Dong-Qi Xue
Xiu-Ling Chen1
Hong Zhang
Jing-Fu Li- 1College of Horticulture, Northeast Agricultural University, Harbin, China
- 2College of Life Science, Northeast Agricultural University, Harbin, China
Cf-12 is an effective gene for resisting tomato leaf mold disease caused by Cladosporium fulvum (C. fulvum). Unlike many other Cf genes such as Cf-2, Cf-4, Cf-5, and Cf-9, no physiological races of C. fulvum that are virulent to Cf-12 carrying plant lines have been identified. In order to better understand the molecular mechanism of Cf-12 gene resistance response, RNA-Seq was used to analyze the transcriptome changes at three different stages of C. fulvum infection (0, 4, and 8 days post infection [dpi]). A total of 9100 differentially expressed genes (DEGs) between 4 and 0 dpi, 8643 DEGs between 8 and 0 dpi and 2547 DEGs between 8 and 4 dpi were identified. In addition, we found that 736 DEGs shared among the above three groups, suggesting the presence of a common core of DEGs in response to C. fulvum infection. These DEGs were significantly enriched in defense-signaling pathways such as the calcium dependent protein kinases pathway and the jasmonic acid signaling pathway. Additionally, we found that many transcription factor genes were among the DEGs, indicating that transcription factors play an important role in C. fulvum defense response. Our study provides new insight on the molecular mechanism of Cf resistance to C. fulvum, especially the unique features of Cf-12 in responding to C. fulvum infection.
Introduction
Cladosporium fulvum (C. fulvum, syn. Passalora fulva) is a biotrophic pathogen of tomato (Solanum lycopersicum), which causes leaf mold disease (Cooke, 1883). This fungus infects primarily the foliage, and occasionally the petioles and stems (Butler and Jones, 1949; Jones et al., 1997). The infection often results in wilting leaves and defoliation, which reduce fruit yield and quality, and sometimes death of the entire plant (Thomma et al., 2005). C. fulvum has many physiological races, and new physiological races continue to be evolved (Westerink et al., 2004). These physiological races differ in race-specific elicitor proteins encoded by effector genes, and one of which is recognized as the AVR gene. The effector proteins are secreted into the apoplastic space during infection (Nekrasov et al., 2006), and induce either a compatible or incompatible interaction between the fungus and infected plant. An incompatible interaction (chlorosis) occurs when the plant is able to resist the pathogen and prevent infection, while a compatible interaction (necrosis) occurs when the pathogen is able to grow and ramify, causing necrosis to the infected cells (Hammond-Kosack and Jones, 1996). One of the most efficient containments is breeding C. fulvum-resistant tomato cultivars by introducing Cf resistance genes identified from wild Solanum species into cultivated tomato. More than 20 Cf genes have been identified since the discovery of the Cf-1 gene in the 1930s (Lanford, 1937; Kanwar et al., 1980a), and these have been introduced into cultivated tomato (Kerr and Bailey, 1964; Kanwar et al., 1980a,b; Stevens and Rick, 1986; Dickinson et al., 1993; Jones et al., 1993; Joosten and de Wit, 1999; Haanstra et al., 2000; Zhao et al., 2016).
The recognition of plant to pathogen and its subsequent response is a complex and dynamic process (Joosten and de Wit, 1999; Rivas and Thomas, 2005). At least three layers of pathogen recognition/response mechanisms are present in plants. The first one is basal resistance, also known as innate immunity, which can be triggered by microbe-associated molecular patterns such as cell wall components found in microbes. When a pathogen suppresses the basal defense, plants may respond with a hypersensitive response (HR) characterized by deliberate cell death at the site of infection. The third layer of defense is called RNA silencing, where plants recognize and digest the DNA or RNA produced by viruses, making these unusable. C. fulvum penetrates the abaxial side of the leaf and secretes toxic avirulence proteins (Avrs) and extracellular proteins. These specific elicitors are recognized by Cf genes (Lauge et al., 1998) and trigger a hypersensitive response. As a result, these infected parenchyma and epidermal cells collapse and eventually form typical necrotic spots (Steinkamp et al., 1979; Feindt et al., 1981). Studies have unraveled the distinct recognition mechanism between Cf genes and Avrs such as the interaction between Cf-2 and Avr2 (Kruger et al., 2002; Luderer et al., 2002) and between Cf-9 and Avr9 (Koomangersmann et al., 1996). Studies have also revealed that gene expression patterns between Cf-4/Avr4- and Cf-9/Avr9-dependent defense responses have similar gene expression patterns (Romeis et al., 2001; Gabriëls et al., 2006; Nekrasov et al., 2006; Hong et al., 2007; van den Burg et al., 2008).
The Cf-12 gene was identified through the large screening of C. fulvum resistant genes (Kanwar et al., 1980b). It is located at position 31 cM on chromosome 8 of the Lycopersicon esculentum. Cf-12 is efficiently resistant to C. fulvum, but little is known about the molecular mechanism of its defense response. In addition, none of the physiological races of C. fulvum are virulent to Cf-12 carrying plant lines, while many physiological races have been found to be virulent to other Cf genes, including Cf-2, Cf-4, Cf-5, and Cf-9.
RNA-Seq has been widely and successfully applied in biological analysis, particularly in plants (Hong et al., 2007; Varshney et al., 2009; Haas and Zody, 2010) such as wheat (Yang et al., 2015), rice (Bai et al., 2015), maize (Li et al., 2010), cabbage (Wang et al., 2016a), and cucumber (Zhang et al., 2014). In this study, the comprehensive transcriptome analysis of Cf-12-tomato at different infection stages was performed to identify differentially expressed genes (DEGs). The identified DEGs were further verified by qRT-PCR and analyzed using gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). Our results can help identify the key genes and pathways associated with Cf-12-mediated resistance response, and better understand the molecular mechanism of Cf resistance to fungal infection.
Materials and Methods
Plants, C. fulvum Strains, and Pathogen Infection
The resistant tomato line of Cf-12-tomato (CGN7495) and the susceptible line Moneymaker were obtained from the Chinese Academy of Agricultural Sciences (Beijing, China). They were grown in a greenhouse at the Horticultural station of Northeast Agricultural University (Harbin, China). The growing condition was 16-h light and 8-h darkness at 25°C with an ambient humidity of 95%. C. fulvum physiological race 1.2.3 was acquired from tomato-growing regions in Harbin using the single sporangiophore transfer method, as previously described (Hubbeling, 1971; Wong and Wilcox, 2000); and was further axenically propagated on Moneymaker. At the four-five leaf stage, the abaxial surface of the Cf-12-tomato seedlings was inoculated with a suspension of 1 × 107 sporangia per ml (Wang et al., 2007). Leaf samples were harvested at 0–15 days of post-infection (dpi) for microscopic analysis.
Microscopic Observation of C. fulvum in Cf-12 Tomato
In order to observe the interaction process of C. fulvum on Cf-12 tomato, the lactophenol trypan blue staining method was carried out according Franco's approach (Franco et al., 2008). The fungal tissue and dead host cells would be densely stained, while living host cells would impart a translucent and slightly brown color. The leaf samples were harvested at 0–15 dpi, immediately stained, clarified overnight in chloral hydrate solution (2.5 mg/ml) (Keogh et al., 1980), and examined using an EVOS® microscope (ThermoFisher, USA) and an OLYMPUS SZX10 dissecting microscope (Olympus, Japan).
RNA Extraction, Library Preparation, and Sequencing
Total RNA from Cf-12-tomato leaves was extracted and analyzed, as previously described (Fang et al., 2015). The integrity of the isolated RNA was calculated, as previously described (Schroeder et al., 2006); and samples with an RNA integrity number >9.0 were used for libraries construction. The libraries were generated using the NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (NEB, USA). The clustering of the index-coded samples was performed on a cBot Cluster Generation System using the TruSeq PE Cluster Kit v3-cBot-HS (Illumia). After cluster generation, the libraries were sequenced on an Illumina Hiseq 4000 platform conducted by the Novogene Bioinformatics Institute (Beijing, China); and 150-bp paired-end reads were generated.
Quality Control, Mapping, and De novo Assembly
Raw sequence data were processed using Perl scripts (http://www.perl.org), developed by Novogene Bioinformatics Institute (Beijing, China), in order to remove reads that contained adapter fragments and ploy-N stretches (the number of ploy-Ns is >10%). Phred quality scores (Q20: ratio of an error rate 1%, Q30: ratio of an error rate 0.1%) and GC-content were calculated, and only the data with a quality score (Qphred) ≥ 30 (Q30) were used for further analyses.
The reference genome and gene model annotation files were downloaded from the Ensembl Genomes Databases (ftp://ftp.ensemblgenomes.org/pub/release-23/plants/fasta/solanum_lycopersicum/dna/). An index of the reference genome was built using Bowtie v2.2.3 (Langmead and Salzberg, 2012), and paired-end clean reads were aligned to the reference genome using TopHat2 v2.0.12 (Kim et al., 2013). Cufflinks v2.1.1 (Trapnell et al., 2010, 2012) was used to construct and identify both known and novel transcripts from TopHat2 alignment results.
Quantification and Differential Expression Analysis of Transcripts
HTSeq v0.6.1 (EMBL, Heidelberg, Germany) was used to count the read numbers mapped to each gene. Gene expression levels were calculated based on the length of the gene, sequencing depth and read count mapped to this gene using the Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced (FPKM; Trapnell et al., 2010) method. Genes with FPKMs in intervals of 1–3, 3–15, 15–60, and beyond 60 were considered to be expressed at low level, medium level, high level and very high level, respectively.
After calculating the gene expression levels, a differential expression analysis of two conditions or groups was performed using the DESeq R package (1.18.0) (Anders and Huber, 2010). DESeq provides statistical strategies to determine differential gene expression using the negative binomial distribution model (Kij ~NB[μij, ]). The resulting P-values and fold-changes were adjusted using the Benjamini-Hochberg's approach to control the false discovery rate (FDR ≤ DRe5; Benjamini and Hochberg, 1995). Genes with an adjusted P-value (padj) < 0.05 were considered to be significant DEGs.
Validation of DEGs by Quantitative Real-Time PCR
Eighteen DEGs involved in plant disease resistance pathways were validated using quantitative real-time PCR (qRT-PCR). The primer pairs of the selected genes were designed using Primer Premier 6.0 (Premier Biosoft, Canada; Table S1). The tomato actin gene (U60478.1) was used as a reference control. The qRT-PCR was performed using AceQ® qPCR SYBR® Green Master Mix (Vazyme, USA) on an iQ™ 5 Multicolor Real-time PCR Detection System (Bio-Rad, USA). The reaction parameters were as follows: 95°C for 7 min, and 40 cycles of 95°C for 10 s, 58°C for 30 s, and 72°C for 20 s. Each sample was repeated three times, and relative expression levels were evaluated using the 2−ΔΔCt method (Livak and Schmittgen, 2001).
GO and KEGG Enrichment Analysis of DEGs
The GO enrichment analysis of DEGs was performed using the GO-seq based on Wallenius non-central hyper-geometric distribution (Young et al., 2010), in which gene length bias in DEGs were adjusted by Bonferroni correction. GO terms with a corrected P < 0.05 were considered significantly enriched by DEGs. KEGG pathway enrichment (http://www.genome.jp/kegg/) was used to identify significantly enriched signal transduction pathways or metabolic pathways in DEGs. Significantly enriched DEG pathways were identified using the KOBAS 2.0 software (KOBAS, Surrey, UK; Xie et al., 2011), and adjusted by hyper-geometric test and Benjamini-Hochberg FDR correction (FDR ≤ F.05).
Results
Microscopy Observation of C. fulvum Invasion into Cf-12 Tomato Leaves
Light microscopy was used to observe the interaction process between C. fulvum and Cf-12 tomato, or the Moneymaker leaves. A representative image of the C. fulvum mycelium and spores are shown in Figure 1A. Our results revealed that conidiospores germinated at 2 or 3 dpi (Figure 1B), the hypha started penetrating into the stomata at 4 dpi in both the Moneymaker and Cf-2 tomato (Figures 1C,D). No difference was observed between the two cultivars at this stage. Then, the hypha penetrated into the intercellular space and continued from the substomatal cavity into the intercellular space between the spongy mesophyll cells (Thomma et al., 2005), emerging through the stomata at 8 dpi and gradually plugging the stomata at 10 dpi on the Moneymaker cultivar (Figures 1E,F). However, cells surrounding the stomata appeared to die (necrosis) at 8 dpi on Cf-12 tomato (Figure 1G); and a large number of necrotic spots began to appear at 10 dpi (Figure 1H). The area of necrosis was gradually enlarged at 12–15 dpi (Figures 1I,J). In the advanced stages of disease development, the hyphae formed a thick, gray-white mold layer on the abaxial surface of the leaves at 15 dpi on the Moneymaker cultivar (Figure 1K) and formed yellow necrotic spots on the front side of the leaf on Cf-12 tomato (Figure 1L). This observation is in agreement with the results of a previous study (Thomma et al., 2005). Based on this observation, we selected Cf-12 tomato samples at 0, 4, and 8 dpi for mRNA-Seq and qRT-PCR analysis.
Figure 1. Lactophenol trypan blue stained-tomato leaf samples after inoculation with C. fulvum. (A) Morphology of Cladosporium fulvum mycelium and spores; (B) conidiospores germinates (2 or 3 dpi); (C,D) the hypha penetrates into the stomata of the Moneymaker (C) and Cf-2 tomato (D) (4 dpi); (E) the hyphae emerges through the stomata of the Moneymaker cultivar (8 dpi); (F) the hyphae increases and plugs the stomata on the Moneymaker cultivar (10 dpi); (G) cells surrounding the stoma had necrosis lesions on Cf-12 tomato (8 dpi); (H–J) a large number of necrotic spots appeared and the necrosis area increased on Cf-12 tomato (10–15 dpi); (K) thick gray-white mold layer on the abaxial surface of the infected leaf on the Moneymaker cultivar (15 dpi); (L) yellow necrotic spots on Cf-12 tomato. s, spore; gt, germ tube; hy, hypha; p, penetration of the stomata by the hypha; ec, plant epidermal cell walls; v, host vascular tissue; st, stomata; nl, necrotic lesion.
RNA Sequencing and Transcripts Identification
In order to investigate differences in transcriptome between C. fulvum-infected Cf-12 tomatoes (Cf12_B, 4 dpi; Cf12_C, 8 dpi) and non-infected controls (Cf12_A, 0 dpi), RNA from these three samples were sequenced. An average of 59,291,565, 56,131,442 and 70,822,873 raw reads from Cf12_A (0 dpi), Cf12_B (4 dpi), and Cf12_C (8 dpi) were generated, respectively (Table 1, Figure S1). After removing adaptors, low quality reads, duplications and ambiguous reads, an average of 60 million clean reads were obtained for each sample; and 94% of these clean reads were at the Q30 (ratio of error rate ≤ 0.1%) level (Table 1). Furthermore, at least 89% of these clean reads were mapped; of which, ~88% were uniquely mapped and 0.7% were multiple mapped to tomato chromosomes, respectively (Table S2). Pearson correlation coefficient (R2) analysis of the FPKM distribution between different biological replicates revealed a high level of reproducibility of RNA expression patterns (R2 = 0.98–0.99, P < 0.001; Table 2, Figure S2). A large majority of reads were mapped to the genome exon regions (Figure S3), and reads density was positively correlated to chromosome length (Figure S4). Cufflinks analysis revealed a total of 15,395 transcripts from the three samples (Cf12_A, Cf12_B, and Cf12_C), including 14,494 known transcripts (Table S3) and 901 new transcripts (Table S4).
Table 1. Quantitative analysis of raw RNA-seq data.
Table 2. FPKM analysis of gene expression levels.
The Identification of Differentially Expressed Genes between the Infected and Non-infected Cf-12 Tomato
DEGs were identified using the DESeq software with a padj <0.05. For Cf12_B vs. Cf12_A, 9100 DEGs were detected, including 4080 upregulated and 5020 downregulated DEGs (Table S5A, Figure 2); for Cf12_C vs. Cf12_A, 8643 DEGs were detected (3999 upregulated and 4644 downregulated DEGs; Table S5B, Figure 2); and for Cf12_C vs. Cf12_B, 2547 DEGs were identified (1729 upregulated and 818 downregulated DEGs; Table S5C, Figure 2). In addition, 736 DEGs were shared among the three groups (Cf12_B vs. Cf12_A, Cf12_C vs. Cf12_A, and Cf12_C vs. Cf12_B). Approximately 1500 DEGs were shared between Cf12_C vs. Cf12_A and Cf12_C vs. Cf12_B, or Cf12_B vs. Cf12_A and Cf12_C vs. Cf12_B. However, 6426 DEGs were shared between Cf12_B vs. Cf12_A and Cf12_C vs. Cf12_A. This further suggests that a common group of genes were activated or deactivated upon C. fulvum infection (Figure 3). In order to observe the overall changes of the quantity of gene expression, a hierarchical clustering of DEGs using FPKM analysis was performed. Results revealed that after the inoculation of C. fulvum on Cf-12 tomato, genes with low expression quantity in Cf12_A increased in expression (in both Cf12_B and Cf12_C). On the contrary, many genes with high expression quantity in Cf12_A decreased in expression in Cf12_B or Cf12_C (Figure 4). Changes in gene expression quantity were also observed between Cf12_B and Cf12_C (Figure 4), suggesting that the response of Cf-12 tomato to C. fulvum infection changes as time progresses.
Figure 2. Volcano plot showing differentially expressed genes between different libraries. Padj < 0.05 was used as the threshold to judge the significance of the difference in gene expression. Red plots represent upregulated genes; green plots represent downregulated genes; blue plots represent genes with no significant difference.
Figure 3. Venn diagram of the relationship between DEG groups. The numbers indicate the DEG number in each DEG group shown in Table S5.
Figure 4. Hierarchical clustering of DEGs. The blue bands indicate low gene expression quantity, and the red bands represent high gene expression quantity.
The Identification of Gene Expression Patterns of DEGs
In order to identify these similar expression patterns, the relative expression level of DEGs were analyzed by K-means clustering algorithm (Hartigan and Wong, 2013). Clustering analysis revealed that six expression patterns (subclusters) of DEGs were identified (Figure 5). The most prominent group was subcluster_5, in which 5343 genes were upregulated after C. fulvum infection. A similar pattern was observed in subcluster_4, where genes revealed a higher expression level in Cf12_B and Cf12_C; however, the number of genes was much lesser than that in subcluster_5. Subcluster_1, subcluster_2 and subcluster_6 revealed a similar expression pattern, in which most of the genes were downregulated in Cf12_B and Cf12_C; but the downregulated level was relatively small. Subcluster_3 contained 41 genes that were downregulated in Cf12_B, but were upregulated in Cf12_C. These dynamic gene expression patterns further suggest that Cf-12 tomato was resistant to C. fulvum via a highly complex process.
Figure 5. The clustering of DEGs expression patterns. The six expression patterns of DEGs obtained by K-means clustering algorithm is shown, which are represented as upregulated (subcluster_4 and subcluster_5), transient (subcluster_3), and downregulated (subcluster_1, subcluster_2, and subcluster_6). Expression ratios are expressed as Log2.
Validation of RNA-Seq Data by qRT-PCR
In order to validate the RNA-Seq data, qRT-PCR was performed through 18 DEGs using three biological replicates. These 18 genes were selected to reflect some of the functional categories and pathways described below (Section GO and KEGG Enrichment Analysis of DEGs), such as the plant-pathogen interaction pathway and plant hormone signal transduction pathway (Table S1). These qRT-PCR results were compared with the RNA-Seq data. As shown in Figure 6, the trends of these gene expression patterns were consistent and had a strong positive correlation coefficient (R2 = 0.9619), indicating that the RNA-Seq data was reliable.
Figure 6. Correlation of expression levels between RNA-Seq and qRT-PCR.
GO and KEGG Enrichment Analysis of DEGs
In order to further characterize the functions of DEGs, GO enrichment analysis was performed using GOseq. The top 10 enrichment terms of the biological process, cellular component and molecular function were selected, respectively, as the main nodes of the directed acyclic graph. In the biological process category, significant terms were enriched in the phosphorus metabolic process, protein phosphorylation biological process, cellular protein modification process, salicylic acid metabolic process, response to jasmonic acid, and salicylic acid biosynthetic process. These processes were generally associated with disease resistance (Figure S5A). Components associated with membranes such as thylakoid, plastid thylakoid membranes, and photosynthetic membranes (Figure S5B) were significantly enriched in the cellular component category. In the molecular function category, transferase activity, phosphotransferase activity and protein serine/threonine kinase activity were enriched in catalytic activity; and purine nucleoside binding, anion binding, nucleoside binding, adenyl nucleotide binding and ATP binding were enriched in binding (Figure S5C). These binding terms and catalytic activity terms played an important role in signal recognition and signal transduction.
One interesting finding is that many transcription factor genes were significantly upregulated or downregulated in response to C. fulvum. A total of 406 transcription factor-related DEGs were identified and annotated into 60 families (Pérez-Rodríguez et al., 2010; Jin et al., 2014; Wang et al., 2016b; Table S6). These 60 transcription factor families were grouped into two clusters (Figure S6). The majority of DEGs in cluster I revealed an upregulation trend at the early stage (Cf12_B) and a downregulated trend at the later stages (Cf12_C). In cluster II, DEGs exhibited a higher expression at both stages or only at the later stage. For example, one of the MYB family transcription factors (Solyc05g053150.1) was upregulated to more than 5-fold in both Cf12_B vs. Cf12_A and Cf12_C vs. Cf12_A (Table S6). This protein, like most of the 406 transcription factors, remains uncharacterized; and is predicted to be a MYB23-like transcription factor. The downregulation of transcription factors was also identified. For example, the OFP family Solyc09g065350.1 (Wang et al., 2005), another uncharacterized transcription factor, was downregulated to ~5-fold in both infection stages. This result suggests that transcription factors play a critical role in response, especially in the early response, to C. fulvum infection.
The KEGG pathway enrichment analysis was used to identify the biological pathways of incompatible interaction. It was found that the plant-pathogen interaction pathway was significantly enriched (FDR ≤ 0.05). As shown in Figure 7, the number of genes and the rich factor are significantly higher than the other pathways. Many other disease-resistance pathways, including stilbenoid and gingerol biosynthesis, phenylpropanoid biosynthesis, ubiquinone and quinone biosynthesis, flavonoid biosynthesis, pentose phosphate pathway, and nitrogen metabolism were also enriched. In addition, it was found that the transduction plant hormone signals, the biosynthesis of unsaturated fatty acids, fatty acid metabolism, and carbon metabolism participate in Cf-12 tomato response to C. fulvum infection.
Figure 7. Scatter plot of the KEGG pathway enrichment of DEGs. Rich factor is the ratio of the DEG number to the background number in a certain pathway. The size of the dots represents the number of genes, and the color of the dots represents the range of the q-value.
Discussion
In this study, we investigated the transcriptome profiles of Cf-12 tomato in response to C. fulvum infection using RNA-Seq. More than 20,000 transcripts were identified from these three samples (Cf12_A, 0 dpi; Cf12_B, 4 dpi; Cf12_C, 8 dpi), in which 9100 DEGs between Cf12_B and Cf12_A, and 8643 DEGs between Cf12_C and Cf12_A were identified. GO and KEGG analyses revealed that many DEGs and their associated pathways are involved in disease resistance against fungal pathogens such as calcium dependent protein kinases (CDPK), wound salicylic acid inducible protein kinase (SIPK), respiratory burst oxidase homolog protein D/B (Rboh), heat shock protein 90 (Hsp90) and suppressor of G2 allele of SKP1 (SGT1; Piedras et al., 1998; Rivas and Thomas, 2005; Nekrasov et al., 2006; Hong et al., 2007).
In order to cope with the infection of C. fulvum, tomato plants have established a series of defense mechanisms through a complex signal transduction network. The first layer of defense is to recognize C. fulvum by pattern recognition receptors, and initiate the resistance response (Miya et al., 2007). In our study, chitin elicitor receptor kinase 1 (CERK1, Solyc07g049190.2, and Solyc07g049180.2), a pattern recognition protein, was abundantly expressed after infection. It would be interesting to conduct further studies to see whether a higher expression of CERK1 is involved in the activation of chitin signaling.
After the recognition of infection, Cf-12 tomato quickly established complex signal defense pathways such as CDPK. In our study, CDPK (Solyc03g113390.2, Solyc10g074570.1, Solyc02g083850.2, and Solyc10g076900.1) and MEKK1 (Solyc01g104530.2, and Solyc07g053170.2) were expressed higher at the early stage of infection, and subsequently stimulated the respiratory burst oxidase homolog (Rboh, Solyc01g099620.2, and Solyc03g117980.2) at a later stage. This is consistent with previous studies suggesting that these genes play critical roles in Cf-12 tomato response to C. fulvum infection.
The jasmonic acid, brassinosteroid, and ethylene pathways play important roles in the resistance to biotrophic pathogens such as downy mildew and powdery (Ellis and Turner, 2001; Walters et al., 2002). In the present study, we found that the jasmonate-zim-domain gene (JAZ, Solyc12g009220.1), which encodes a major protein in the jasmonic acid signaling pathway, was upregulated after C. fulvum infection. Brassinosteroid insensitive 1-associated receptor kinase 1 (BAK1, Solyc01g104970.2), brassinosteroid-signaling kinase (BSK, Solyc10g085000.1), ethylene receptor (ETR, Solyc06g053710.2) and ethylene response factor 1/2 (ERF1/2, Solyc09g066360.1) were also upregulated; suggesting that similar to previous studies, jasmonic acid, brassinosteroid, and ethylene may play a role in the resistance of Cf-12 tomato to C. fulvum.
Physiological races of C. fulvum that are virulent to Cf-2, Cf-4, Cf-5, and Cf-9 have been reported (Piedras et al., 1998; Romeis et al., 2001; Nekrasov et al., 2006; Hong et al., 2007; Varshney et al., 2009). However, no physiological race virulent to Cf-12 has been identified. In the present study, we revealed that genes encoding the NPR1-like protein (Solyc07g040690.2 and Solyc02g069310.2) and the transcription factor TGA (Solyc11g068370.1 and Solyc06g074320.2) were upregulated, and both of which are involved in the salicylic acid signaling pathway. This could lead to a higher expression of pathogenesis-related protein 1 (PR1) during the disease resistance process. A previous report suggested that the salicylic acid pathway is not required for Cf-2- or Cf-9-dependent resistance to C. fulvum (Brading et al., 2000), suggesting that the salicylic acid pathway might be a unique pathway for Cf-12-dependent resistance.
Responses to infection depend on different gene expression levels, which require various transcription factors. It has been reported that many transcription factors such as WRKY, NAC, MYB, and bZIP families (Singh et al., 2002; Olsen et al., 2005) actively respond to C. fulvum infections. Consistent with previous reports, we found that 35 transcription factors in the MYB family were differentially expressed; and one of which (Solyc05g053150.1) had the highest level of upregulation among all transcription factor DEGs. Many TF families such as TUB, TCP, SET, SBP, PHD, and Orphans have not been reported in the regulation of tomato leaf mold disease resistance, but some of these TFs have been reported in grapevine (Vitis amurensis) against downy mildew (Li et al., 2015) and in ramie against root-lesion nematode infection (Zhu et al., 2014). In addition, majority of the DEGs reported in this study remain uncharacterized. Further studies of these transcription factors related DEGs could greatly help understand the molecular mechanism of the defense response to C. fulvum infection.
In summary, based on our microscopic data and GO/KEGG analysis of RNA-seq data, we believe that the Cf-12 tomato response to C. fulvum infection follows several steps. When the mycelium of C. fulvum grows into the interspace of the stomata and mesophyll cells, the effector proteins secreted by C. fulvum are rapidly recognized by Cf-12 tomato. This triggers downstream defense signaling transductions associated with the Ca2+ channel, as well as several pathways, including pathways involving jasmonic acid, brassinosteroid and ethylene. Then, these defense-related transcription factors (TFs) such as MYB proteins (Moore et al., 2011; van Verk et al., 2011; Puranik et al., 2012) are triggered; which actively regulate downstream resistance pathways. Finally, hypersensitive necrosis occurs locally (Hammond-Kosack and Jones, 1996), and the C. fulvum hyphae is restricted only in infected areas. This study facilitates our understanding of the molecular mechanism of Cf-12 tomato against C. fulvum infection.
Database Link and Accessions
The raw sequencing data of the nine samples have been submitted to the NCBI Sequence Read Archive (SRA, http://www.ncbi.nlm.nih.gov/sra). The accession numbers are: SRR4041970, SRR4041973, SRR4041974, SRR4041975, SRR4042017, SRR4042029, SRR4042030, SRR4042031, and SRR4042032.
Author Contributions
DX, JL, XX, and JJ conceived and designed the experiments; HZ and XFC performed the RNA isolation and qRT-PCR experiments; DX and XLC performed the data analysis; and DX wrote the manuscript. All authors read and approved the final manuscript.
Conflict of Interest Statement
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.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (Grant No. 31572137). And we thank Medjaden Bioscience Limited for scientific editing of this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.02012/full#supplementary-material
Figure S1. Classification of raw reads.
Figure S2. RNA-Seq correlation analysis between samples.
Figure S3. The percentage of reads mapped to genome regions.
Figure S4. Reads density on chromosomes of the tomato reference genome.
Figure S5. Significant enrichment of GO terms. (A) GO significant enrichment in biological process; (B) GO significant enrichment in cellular component; (C) GO significant enrichment in molecular function.
Figure S6. Response of transcription factors to C. fulvum infection. Each line color corresponds to the same gene from Cf12_B/Cf12_A and Cf12_C/Cf12_A. The expression level of the same gene is represented by different colors.
Table S1. Primers used for qPCR.
Table S2. Summary of RNA-Seq data and mapping.
Table S3. Structural optimization of known genes.
Table S4. Results of the new transcription structure annotation.
Table S5. DEGs between two samples. (A) DEGs between Cf12_B and Cf12_A; (B) DEGs between Cf12_C and Cf12_A; (C) DEGs between Cf12_C and Cf12_B.
Table S6. Response of transcription factors to C. fulvum infection.
References
Anders, S., and Huber, W. (2010). Differential expression analysis for sequence count data. Genome Biol. 11:R106. doi: 10.1186/gb-2010-11-10-r106
Bai, B., Wu, J., Sheng, W. T., Zhou, B., Zhou, L. J., Zhuang, W., et al. (2015). Comparative analysis of anther transcriptome profiles of two different rice male sterile lines genotypes under cold stress. Int. J. Mol. Sci. 16, 11398–11416. doi: 10.3390/ijms160511398
Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate - a practical and powerful approach to multiple testing. J. R. Stat. Soc. 57, 289–300.
Brading, P. A., Hammond-Kosack, K. E., Parr, A., and Jones, J. D. (2000). Salicylic acid is not required for Cf-2 - and Cf-9 -dependent resistance of tomato to Cladosporium fulvum. Plant J. Cell Mol. Biol. 23, 305–318. doi: 10.1046/j.1365-313x.2000.00778.x
Butler, E. J., and Jones, S. G. (1949). Tomato Leaf Mould, Cladosporium fulvum Cooke. London: Macmillan.
Dickinson, M. J., Jones, D. A., and Jones, J. D. (1993). Close linkage between the Cf-2/Cf-5 and Mi resistance loci in tomato. Mol. Plant Microbe Interact. 6, 341–347. doi: 10.1094/MPMI-6-341
Ellis, C., and Turner, J. G. (2001). The arabidopsis mutant cev1 has constitutively active jasmonate and ethylene signal pathways and enhanced resistance to pathogens. Plant Cell 13, 1025–1033. doi: 10.1105/tpc.13.5.1025
Fang, S. M., Hu, B. L., Zhou, Q. Z., Yu, Q. Y., and Zhang, Z. (2015). Comparative analysis of the silk gland transcriptomes between the domestic and wild silkworms. BMC Genomics 16:60. doi: 10.1186/s12864-015-1287-9
Feindt, F., Mendgen, K., and Heitefuss, R. (1981). Feinstruktur unterschiedlicher Zellwandreaktionen im Blattparenchym anfälliger und resistenter Rüben (Beta vulgaris L.) nach Infektion durch Cercospora beticola Sacc. J. Phytopathol. 101, 248–264. doi: 10.1111/j.1439-0434.1981.tb03346.x
Franco, F., Dario, M., Dario, C., and Marcello, I. (2008). Chemical-induced resistance against powdery mildew in barley: the effects of chitosan and benzothiadiazole. BioChontrol 53, 387–401. doi: 10.1007/s10526-007-9091-3
Gabriëls, S. H., Takken, F. L., Vossen, J. H., de Jong, C. F., Liu, Q., Turk, S. C., et al. (2006). CDNA-AFLP combined with functional analysis reveals novel genes involved in the hypersensitive response. Mol. Plant Microbe Interact. 19, 567–576. doi: 10.1094/MPMI-19-0567
Haanstra, J. P., Thomas, C. M., Jones, J. D., and Lindhout, P. (2000). Dispersion of the Cf-4 disease resistance gene in Lycopersicon germplasm. Heredity 85(Pt 3), 266–270. doi: 10.1046/j.1365-2540.2000.00754.x
Haas, B. J., and Zody, M. C. (2010). Advancing RNA-Seq analysis. Nat. Biotechnol. 28, 421–423. doi: 10.1038/nbt0510-421
Hammond-Kosack, K. E., and Jones, J. D. (1996). Resistance gene-dependent plant defense responses. Plant Cell 8, 1773–1791. doi: 10.1105/tpc.8.10.1773
Hartigan, J. A., and Wong, M. A. (2013). A K-means clustering algorithm. Appl. Stat. 28, 100–108. doi: 10.2307/2346830
Hong, W., Xu, Y. P., Zheng, Z., Cao, J. S., and Cai, X. Z. (2007). Comparative transcript profiling by cDNA-AFLP reveals similar patterns of Avr4/Cf-4- and Avr9/Cf-9-dependent defence gene expression. Mol. Plant Pathol. 8, 515–527. doi: 10.1111/j.1364-3703.2007.00412.x
Hubbeling, N. (1971). Determination trouble with new races of Cladosporium fulvum Cooke. Meded. Fac. Landbouwwet. Gent. 36, 300–305.
Jin, J., Zhang, H., Kong, L., Gao, G., and Luo, J. (2014). PlantTFDB 3.0: a portal for the functional and evolutionary study of plant transcription factors. Nucleic Acids Res. 42, 1182–1187. doi: 10.1093/nar/gkt1016
Jones, D. A., Dickinson, M. J., Balint-Kurti, P. J., Dixon, M. S., and Jones, J. D. G. (1993). Two complex resistance loci revealed in tomato by classical and RFLP mapping of the Cf-2, Cf-4, Cf-5, and Cf-9 genes for resistance to Cladosporium fulvum. Mol. Plant Microbe Interact. 6, 348–357. doi: 10.1094/MPMI-6-348
Jones, J. B., Jones, J. P., Stall, R. E., and Zitter, T. A. (1997). Compendium of Tomato Diseases. St. Paul. MN: APS Press.
Joosten, M., and de Wit, P. (1999). The tomato-Cladosporium fulvum interaction: a versatile experimental system to study plant-pathogen interactions. Annu. Rev. Phytopathol. 37, 335–367. doi: 10.1146/annurev.phyto.37.1.335
Kanwar, J. S., Kerr, E. A., and Harney, P. M. (1980a). Linkage of Cf-1 to Cf-11 genes for resistance to leaf mold Cladosporium fulvum. Tomato Genet. Coop 30, 20–21.
Kanwar, J. S., Kerr, E. A., and Harney, P. M. (1980b). Linkage of Cf-12 to Cf-24 genes for resistance to leaf mold Cladosporium fulvum Cke Rep. Tomato Genet. Coop 30, 22–23.
Keogh, R. C., Deverall, B. J., and McLeod, S. (1980). Comparison of histological and physiological responses to Phakopsora pachyrhizi in resistant and susceptible soybean. Trans. Br. Mycol. Soc. 74, 329–333. doi: 10.1016/S0007-1536(80)80163-X
Kerr, E. A., and Bailey, D. L. (1964). Resistance to Cladosporium fulvum Cke. obtained from wild species of tomato. Can. J. Bot. 42, 1541–1554. doi: 10.1139/b64-153
Kim, D., Pertea, G., Trapnell, C., Pimentel, H., Kelley, R., and Salzberg, S. L. (2013). TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 14:R36. doi: 10.1186/gb-2013-14-4-r36
Koomangersmann, M., Honee, G. G., and Pjgm, D. (1996). A high affinity binding site for the AVR9 peptide elicitor of Cladosporium fulvum is present on plasma membranes of resitant and susceptible tomato genotypes. Plant Cell 8, 929–938. doi: 10.1105/tpc.8.5.929
Kruger, J., Thomas, C. M., Golstein, C., Dixon, M. S., Smoker, M., Tang, S., et al. (2002). A tomato cysteine protease required for Cf-2-dependent disease resistance and suppression of autonecrosis. Science 296, 744–747. doi: 10.1126/science.1069288
Lanford, A. N. (1937). The parasitism of i Cooke and the genetics of resistance to it. Can. J. Res. C 15, 108–128. doi: 10.1139/cjr37c-008
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi: 10.1038/nmeth.1923
Lauge, R., Joosten, M. H., Haanstra, J. P., Goodwin, P. H., Lindhout, P., and De Wit, P. J. (1998). Successful search for a resistance gene in tomato targeted against a virulence factor of a fungal pathogen. Proc. Natl. Acad. Sci. U.S.A. 95, 9014–9018. doi: 10.1073/pnas.95.15.9014
Li, P., Ponnala, L., Gandotra, N., Wang, L., Si, Y., Tausta, S. L., et al. (2010). The developmental dynamics of the maize leaf transcriptome. Nat. Genet. 42, 1060–1067. doi: 10.1038/ng.703
Li, X., Wu, J., Yin, L., Zhang, Y., Qu, J., and Lu, J. (2015). Comparative transcriptome analysis reveals defense-related genes and pathways against downy mildew in Vitis amurensis grapevine. Plant Physiol. Biochem. 95, 1–14. doi: 10.1016/j.plaphy.2015.06.016
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)). Methods 25, 402–408. doi: 10.1006/meth.2001.1262
Luderer, R., Takken, F. L. W., de Wit, P. J. G. M., and Joosten, M. H. A. J. (2002). Cladosporium fulvum overcomes Cf-2-mediated resistance by producing truncated AVR2 elicitor proteins: cloning and characterization of avirulence gene AVR2. Mol. Microbiology 45, 875–884. doi: 10.1046/j.1365-2958.2002.03060.x
Miya, A., Albert, P., Shinya, T., Desaki, Y., Ichimura, K., Shirasu, K., et al. (2007). CERK1, a LysM receptor kinase, is essential for chitin elicitor signaling in Arabidopsis. Proc. Natl. Acad. Sci. U.S.A. 104, 19613–19618. doi: 10.1073/pnas.0705147104
Moore, J. W., Loake, G. J., and Spoel, S. H. (2011). Transcription dynamics in plant immunity. Plant Cell 23, 2809–2820. doi: 10.1105/tpc.111.087346
Nekrasov, V., Ludwig, A. A., and Jones, J. D. (2006). CITRX thioredoxin is a putative adaptor protein connecting Cf-9 and the ACIK1 protein kinase during the Cf-9/Avr9- induced defence response. FEBS Lett. 580, 4236–4241. doi: 10.1016/j.febslet.2006.06.077
Olsen, A. N., Ernst, H. A., Leggio, L. L., and Skriver, K. (2005). NAC transcription factors: structurally distinct, functionally diverse. Trends Plant Sci. 10, 79–87. doi: 10.1016/j.tplants.2004.12.010
Pérez-Rodríguez, P., Riaño-Pachón, D. M., Corrêa, L. G. G., Rensing, S. A., Kersten, B., and Mueller-Roeber, B. (2010). PlnTFDB: updated content and new features of the plant transcription factor database. Nucleic Acids Res. 38, D822–D827. doi: 10.1093/nar/gkp805
Piedras, P., Hammond-Kosack, K. E., Harrison, K., and Jones, J. D. G. (1998). Rapid, Cf-9- and Avr9-Dependent production of active oxygen species in tobacco suspension cultures. Mol. Plant Microbe Interact. 11, 1155–1166. doi: 10.1094/MPMI.1998.11.12.1155
Puranik, S., Sahu, P. P., Srivastava, P. S., and Prasad, M. (2012). NAC proteins: regulation and role in stress tolerance. Trends Plant Sci. 17, 369–381. doi: 10.1016/j.tplants.2012.02.004
Rivas, S., and Thomas, C. M. (2005). Molecular interactions between tomato and the leaf mold pathogen Cladosporium fulvum. Annu. Rev. Phytopathol. 43, 395–436. doi: 10.1146/annurev.phyto.43.040204.140224
Romeis, T., Ludwig, A. A., Martin, R., and Jones, J. D. (2001). Calcium-dependent protein kinases play an essential role in a plant defence response. EMBO J. 20, 5556–5567. doi: 10.1093/emboj/20.20.5556
Schroeder, A., Mueller, O., Stocker, S., Salowsky, R., Leiber, M., Gassmann, M., et al. (2006). The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol. Biol. 7:3. doi: 10.1186/1471-2199-7-3
Singh, K., Foley, R. C., and Oñate-Sánchez, L. (2002). Transcription factors in plant defense and stress responses. Curr. Opin. Plant Biol. 5, 430–436. doi: 10.1016/S1369-5266(02)00289-3
Steinkamp, M. P., Martin, S. S., Hoefert, L. L., and Ruppel, E. G. (1979). Ultrastructure of lesions produced by Cercospora beticola in leaves of Beta vulgaris. Physiol. Plant Pathol. 15, 13–26. doi: 10.1016/0048-4059(79)90035-3
Thomma, B. P., Van Esse, H. P., Crous, P. W., and De Wit, P. J. (2005). Cladosporium fulvum (syn. Passalora fulva), a highly specialized plant pathogen as a model for functional studies on plant pathogenic Mycosphaerellaceae. Mol. Plant Pathol. 6, 379–393. doi: 10.1111/j.1364-3703.2005.00292.x
Trapnell, C., Roberts, A., Goff, L., Pertea, G., Kim, D., Kelley, D. R., et al. (2012). Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578. doi: 10.1038/nprot.2012.016
Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., van Baren, M. J., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28, 511–515. doi: 10.1038/nbt.1621
van den Burg, H. A., Tsitsigiannis, D. I., Rowland, O., Lo, J., Rallapalli, G., Maclean, D., et al. (2008). The F-box protein ACRE189/ACIF1 regulates cell death and defense responses activated during pathogen recognition in tobacco and tomato. Plant Cell 20, 697–719. doi: 10.1105/tpc.107.056978
van Verk, M. C., Bol, J. F., and Linthorst, H. J. (2011). Prospecting for genes involved in transcriptional regulation of plant defenses, a bioinformatics approach. BMC Plant Biol. 11:88. doi: 10.1186/1471-2229-11-88
Varshney, R. K., Nayak, S. N., May, G. D., and Jackson, S. A. (2009). Next-generation sequencing technologies and their implications for crop genetics and breeding. Trends Biotechnol. 27, 522–530. doi: 10.1016/j.tibtech.2009.05.006
Walters, D., Cowley, T., and Mitchell, A. (2002). Methyl jasmonate alters polyamine metabolism and induces systemic protection against powdery mildew infection in barley seedlings. J. Exp. Bot. 53, 747–756. doi: 10.1093/jexbot/53.369.747
Wang, A., Meng, F., Xu, X., Wang, Y., and Li, J. (2007). Development of molecular markers linked to Cladosporium fulvum resistant gene Cf-6 in tomato by RAPD and SSR methods. Hortscience 42, 11–15.
Wang, S., Wang, C., Zhang, X. X., Chen, X., Liu, J. J., Jia, X. F., et al. (2016a). Transcriptome de novo assembly and analysis of differentially expressed genes related to cytoplasmic male sterility in cabbage. Plant Physiol. Biochem. 105, 224–232. doi: 10.1016/j.plaphy.2016.04.027
Wang, Y., Hoeven, R. S. V. D., Nielsen, R., Mueller, L. A., and Tanksley, S. D. (2005). Characteristics of the tomato nuclear genome as determined by sequencing undermethylated EcoR I digested fragments. Theor. Appl. Genet. 112, 72–84. doi: 10.1007/s00122-005-0107-z
Wang, Y., Lu, W., and Deng, D. (2016b). Bioinformatic landscapes for plant transcription factor system research. Planta 243, 297–304. doi: 10.1007/s00425-015-2453-7
Westerink, N., Joosten, M. H., and De Wit, P. J. (2004). Fungal (A)Virulence Factors at the Crossroads of Disease Susceptibility and Resistance. Binghamton, NY: Food Products Press
Wong, F. P., and Wilcox, W. F. (2000). Distribution of baseline sensitivities to azoxystrobin among isolates of Plasmopara viticola. Plant Dis. 84, 275–281. doi: 10.1094/PDIS.2000.84.3.275
Xie, C., Mao, X., Huang, J., Ding, Y., Wu, J., Dong, S., et al. (2011). KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 39, W316–W322. doi: 10.1093/nar/gkr483
Yang, Z., Peng, Z., Wei, S., Liao, M., Yu, Y., and Jang, Z. (2015). Pistillody mutant reveals key insights into stamen and pistil development in wheat (Triticum aestivum L.). BMC Genomics 16:211. doi: 10.1186/s12864-015-1453-0
Young, M. D., Wakefield, M. J., Smyth, G. K., and Oshlack, A. (2010). Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 11:R14. doi: 10.1186/gb-2010-11-2-r14
Zhang, N., Zhang, H. J., Zhao, B., Sun, Q. Q., Cao, Y. Y., Li, R., et al. (2014). The RNA-seq approach to discriminate gene expression profiles in response to melatonin on cucumber lateral root formation. J. Pineal Res. 56, 39–50. doi: 10.1111/jpi.12095
Zhao, T., Jiang, J., Liu, G., He, S., Zhang, H., Chen, X., et al. (2016). Mapping and candidate gene screening of tomato Cladosporium fulvum-resistant gene Cf-19, based on high-throughput sequencing technology. BMC Plant Biol. 16:51. doi: 10.1186/s12870-016-0737-0
Keywords: Cladosporium fulvum, resistance response, RNA-Seq, differentially expressed genes, Cf-12 tomato
Citation: Xue D-Q, Chen X-L, Zhang H, Chai X-F, Jiang J-B, Xu X-Y and Li J-F (2017) Transcriptome Analysis of the Cf-12-Mediated Resistance Response to Cladosporium fulvum in Tomato. Front. Plant Sci. 7:2012. doi: 10.3389/fpls.2016.02012
Received: 29 August 2016; Accepted: 19 December 2016;
Published: 05 January 2017.
Edited by:
Adi Avni, Tel Aviv University, IsraelReviewed by:
Mahmut Tör, University of Worcester, UKOswaldo Valdes-Lopez, National Autonomous University of Mexico, Mexico
Copyright © 2017 Xue, Chen, Zhang, Chai, Jiang, Xu and Li. 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) or licensor 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.
*Correspondence: Xiang-Yang Xu, xxy709@126.com
Jing-Fu Li, lijf_2005@126.com