Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

RNA-Seq derived identification of differential transcription in the chrysanthemum leaf following inoculation with Alternaria tenuissima

BMC Genomics volume 15, Article number: 9 (2014) Cite this article

Abstract

Background

A major production constraint on the important ornamental species chrysanthemum is black spot which is caused by the necrotrophic fungus Alternaria tenuissima. The molecular basis of host resistance to A. tenuissima has not been studied as yet in any detail. Here, high throughput sequencing was taken to characterize the transcriptomic response of the chrysanthemum leaf to A. tenuissima inoculation.

Results

The transcriptomic data was acquired using RNA-Seq technology, based on the Illumina HiSeq™ 2000 platform. Four different libraries derived from two sets of leaves harvested from either inoculated or mock-inoculated plants were characterized. Over seven million clean reads were generated from each library, each corresponding to a coverage of >350,000 nt. About 70% of the reads could be mapped to a set of chrysanthemum unigenes. Read frequency was used as a measure of transcript abundance and therefore as an identifier of differential transcription in the four libraries. The differentially transcribed genes identified were involved in photosynthesis, pathogen recognition, reactive oxygen species generation, cell wall modification and phytohormone signalling; in addition, a number of varied transcription factors were identified. A selection of 23 of the genes was transcription-profiled using quantitative RT-PCR to validate the RNA-Seq output.

Conclusions

A substantial body of chrysanthemum transcriptomic sequence was generated, which led to a number of insights into the molecular basis of the host response to A. tenuissima infection. Although most of the differentially transcribed genes were up-regulated by the presence of the pathogen, those involved in photosynthesis were down-regulated.

Background

Chrysanthemum (Chrysanthemum morifolium Ramat.) is the second most commercially valuable ornamental species after rose [1, 2]. A serious production constraint is represented by black spot disease (causative pathogen the necrotrophic fungus Alternaria tenuissima (Fr.) Wiltsh) [3]. The disease is most damaging during humid, warm conditions, which makes it a year-round problem for greenhouse-based production [3, 4]. Severe infections damage the commercial value of the plant, as they cause leaf necrosis, and reduce the quantity and quality of the flowers [5]. Little is known regarding the chrysanthemum/black spot host-pathogen interaction. However, in the host-pathogen system involving the model species Arabidopsis thaliana and the related pathogen A. brassicicola, the pathogenesis-related protein PR4 is significantly up-regulated [6]. In the tomato/A. alternata system, the ethylene (ET), jasmonate (JA) and salicylic acid (SA) signalling pathways are all activated as part of the host response [7], while the response of mint to A. alternata infection features many proteins related to stress and defence [8]. Finally, Egusa et al. (2009) [9] have shown that the transcription of the genes PGIP (polygalacturonase inhibiting protein) and PPO (polyphenol oxidase) is induced in the leaf of the Japanese pear when challenged by A. alternata. Plant defence responses are first activated in the organs located at the site of infection but are then extended to the uninfected systemic (distal) organs, activating a systemic acquired resistance (SAR) which is effective against a broad spectrum of pathogens in the whole plant [10–12]. However, limited studies have examined the induction of SAR in chrysanthemum/black spot host-pathogen interaction, so far. In the present study, the systemic responses are expected by surveying gene expression profiles in the noninfected (systemic) leaves.

RNA-Seq technology has been developed to enable the simultaneous sequencing of very large numbers of short reads, and in so doing has revolutionized the qualitative and quantitative analysis of the transcriptome [13–16]. When applied to cotton infected with a wilt pathogen, of the >32,000 genes identified by mapping the reads to a genomic sequence assembly, over 3,000 were found to be up- or down-regulated as part of the defence response [16]. Similarly, an analysis of the banana-Fusarium oxysporum interaction successfully demonstrated the up-regulation of genes involved in hormone synthesis, pathogenesis-related genes, transcription factors and signalling/regulatory genes [17]. Finally, the lettuce-Botrytis cinerea interaction has been shown to feature the induction of genes involved in the phenylpropanoid pathway and in terpenoid synthesis, as well as a global down-regulation of genes responsible for photosynthesis [18].

In present study, we aimed to 1) elucidate the localized responses to the infections on the inoculated site by comparing libraries generated from mock-inoculated and inoculated leaves; 2) to describe the systemic response by comparing libraries generated from neighbouring leaves from mock-inoculated and inoculated leaves. The present study reports the outcome of a RNA-Seq based analysis of the chrysanthemum-A. tenuissima interaction. The experiment has yielded information regarding the identity of the genes which are either up- or down-regulated as part of the defence response. The majority of the differentially transcribed (DT) genes were involved in either pathogen recognition, reactive oxygen species detoxification, cell wall modification or phytohormone signalling, but also a range of transcription factors, belonging to various families were identified. Validation of the RNA-Seq data was provided by subjecting a set of 23 of the DT genes to quantitative RT-PCR (qPCR).

Results

Analysis of RNA-Seq libraries

The major characteristics of the four libraries (Figure 1) are summarized in Table 1 and Additional file 1: Figure S1. The number of raw reads per library ranged from ~7.2 to ~7.6 million, and the total number of base pairs sequenced from 352,864,582 to 370,177,360 (Table 1) (Accession No. for library A SRS464569; Accession No. for library B SRS480632; Accession No. for library C SRS480633; Accession No. for library D SRS480635). After removal of reads including adaptor sequence, reads in which >10% of the bases were uncertain, a total of, respectively 7,524,234, 7,248,778, 7,201,318 and 7,554,640 clean reads were obtained, corresponding to 368,687,466, 355,190,122, 352,864,582 and 370,177,360 base pairs (Table 1). The proportion of clean reads was >99.30% in each library (Additional file 1: Figure S1).

Figure 1
figure 1

The source of RNA libraries prepared from mock- and pathogen-inoculated chrysanthemum plants. A and B: mock inoculation, C and D: inoculation with an A. tenuissima spore suspension. A, C: The first and fourth true leaves without any treatment. B: the second and third mock-inoculated true leaves. D: the second and third pathogen-inoculated true leaves. Leaves were harvested at 0 h, 6 h, 24 h, 48 h and 72 h after inoculation.

Full size image
Table 1 Summary of read mapping
Full size table

Read mapping

A reference gene database [Raw sequence data were deposited in the NCBI Sequence Read Archive database (http://trace.ncbi.nlm.nih.gov/Traces/sra_sub/sub.cgi?) under accession number SRP029991] which included all known Chrysanthemum morifolium unigene sequences was used to map the RNA-Seq reads. Based on the chosen criteria, 69.75% of the clean reads recognized sequences in this database (Table 1). On a per library basis, the proportions of the clean reads uniquely mapped to the database were, respectively, 46.81%, 45.82%, 46.14% and 44.56%. In addition, the proportion of the clean reads from library D uniquely mapped to the publicly available A.tenuissima database was 0.02%. The number of genes identified increased with the number of reads, but above 6,000,000 reads no further genes were detected, implying full saturation of the transcriptome (Figure 2).

Figure 2
figure 2

Sequencing saturation analysis of each library. Sequencing saturation in the four libraries of A, B, C and D. The number of different genes detected rose as the read number was increased.

Full size image

GO classification of DT genes

Of the 217 genes classified as DT genes in the contrast between library A (leaf 1 and 4 of plants mock-inoculated on leaf 2 and 3) and B (mock-treated leaf 2 and 3), 106 could be assigned a GO classification; the equivalent number for the A vs C (leaf 1 and 4 of plants inoculated by the pathogen on leaf 2 and 3) contrast was 418 out of 659, for the B vs D (pathogen-infected leaf 2 and 3) contrast 1,057 out of 1,705, and for the C vs D contrast 294 out of 494 (Additional file 2: Table S1, Additional file 3: Table S2, Additional file 4: Table S3 and Additional file 5: Table S4, Additional file 6: Table S5, Additional file 7: Table S6, Additional file 8: Table S7 and Additional file 9: Table S8). Furthermore, of the 659 genes classified as DT genes in the contrast between library A and C (A vs C), 469 (71.2%) behaved similarly between A vs C and B vs D; the equivalent number for the B vs D contrast was 469 out of 1705 (27.5%) (Additional file 7: Table S6 and Additional file 8: Table S7). For the A vs B contrast, seven genes were categorized as “cellular component”, six as “molecular function” and 14 as “biological process”; the respective distributions in A vs C, B vs D and C vs D were ten, ten and 20, nine, 11 and 20, and nine, nine and 18 (Figure 3). The frequency of DT genes was highest in the contrast B vs D. The commonest molecular functions of the DT genes in this contrast were binding and catalytic activity; in terms of cellular component, most were associated with cells, cell parts, macromolecular complexes, membranes, membrane parts, organelles and organelle parts; finally, in terms of biological process, the majority were associated with cellular processes, metabolic processes and the response to stimulus (Figure 3). Signalling responses to or mediated by JA, SA and ET were well represented, particularly those active in JA-mediated signalling. Some transcription factors and cell wall modification genes were also differentially transcribed, as were genes involved in secondary metabolism (phenylpropanoid pathway and terpenoid synthesis). Among the biological processes well represented among the down-regulated genes were those involved in photosynthesis. In the contrast between B and D, a small number of up-regulated genes belonged to the categories ‘cell killing’, ‘positive regulation of biological process’, ‘membrane-enclosed lumen’ and ‘receptor activity’. A higher number of genes identified in this contrast were associated with ‘response to stimulus’ and ‘signalling’ than in the contrasts A vs B, A vs C and C vs D (Figure 3).

Figure 3
figure 3

Gene Ontology (GO) classifications of DT genes. DT genes were annotated in three categories: biological process, cellular component and molecular function. Y-axis (right) represents the number of DT genes in each category; Y-axis (left) represents the percentage of a specific category of DT genes within that main category. Panels A, B, C and D (left) represents DT genes in the contrast between library A (leaf 1 and 4 of plants mock-inoculated on leaf 2 and 3) and B (mock-treated leaf 2 and 3) (A-VS-B) (right) , library A and C (leaf 1 and 4 of plants inoculated by the pathogen on leaf 2 and 3) (A-VS-C) (right), library B and D (pathogen-infected leaf 2 and 3) (B-VS-D) (right), library C and D (C-VS-D) (right), respectively.

Full size image

Changes in transcription level

The distribution of unigene coverage in each sample was analysed as a way of evaluating the quality of the RNA-Seq dataset (Figure 4). The term “gene coverage” reflects the proportion of the full gene sequence represented by RNA-Seq reads. For most of the unigenes, gene coverage was >50%. The transcription level of each unigene (Additional file 10: Table S9) was derived from the number of relevant reads recovered following Mortazavi et al. (2008) [19]. DT genes (Additional file 6: Table S5, Additional file 7: Table S6, Additional file 8: Table S7, Additional file 9: Table S8) were identified using an algorithm developed by Audic et al. (1997) [20]. Between A and B, 20 genes were up- and 197 down-regulated, between A and C, the totals were 562 and 97, respectively, between B and D, 1,181 and 524, respectively and C and D, 245 and 249, respectively (Figure 5).

Figure 4
figure 4

Distribution of gene coverage analysis of each library. Distribution of gene coverage in libraries A, B, C and D.

Full size image
Figure 5
figure 5

The number of DT unigenes identified in each library contrast.

Full size image

Transcription factors, cell wall modification genes and genes involved in JA and SA signalling were all regulated by A. tenuissima inoculation

The transcription data indicated that infection by A. tenuissima regulated a number of transcription factors and genes associated with pathogenesis and JA and SA signalling. Three members of the GRAS-type transcription factor family were more abundantly transcribed in the inoculated than in the mock-inoculated leaves. When qPCR was performed on a selection of the 23 of the DT genes to validate the conclusions drawn from the RNA-Seq analysis, all the genes behaved as predicted. Most of the genes examined (the exception was the MYB transcription factor CL9570) were induced after A. tenuissima inoculation (Table 2). The qPCR analysis suggested that genes involved in cell wall modification, JA and SA signalling and transcription factors comprised a network of interactions, providing the host with a capacity to fine-tune its disease response.

Table 2 Genes differentially transcribed in ‘Zaoyihong’ leaves in response to mock-inoculation (B) and A. tenuissima inoculation (D)
Full size table

Discussion

Global patterns of transcription in response to infection by A. tenuissima

The chrysanthemum genome is polyploid and large, so has not as yet benefited from comprehensive and integrated genomic and transcriptomic sequence analysis. The molecular basis of its defence response against pathogen infection is currently poorly understood, but the advent of high-throughput sequencing technology now allows an unprecedented opportunity to explore it. About 70% of the reads in each of the four RNA-Seq libraries were mappable back to known transcripts (Table 1), a proportion which is somewhat lower than achieved in the lettuce/B. cinerea system [18], probably reflecting the more comprehensive status of the lettuce transcriptome. The 30% of reads which were not mappable are presumably associated with as yet unidentified transcripts [16]. The validation through qPCR of the transcripts identified as regulated by A. tenuissima infection showed that the RNA-Seq method is well suited for the analysis of transcription induced as part of the defence response in chrysanthemum (Table 2). In all, 659 DT genes were identified in the A vs C contrast and 1,705 in B vs D (Figure 5, Tables S7, S8). A comparison with the outcomes of the lettuce/B. cinerea and A. thaliana/B. cinerea interactions [18, 21, 22] showed that only two genes were up-regulated in all three systems, namely Lsa004290.1, At1g74360 and Unigene7965_All (encoding an LRR protein kinase) and Lsa016859.1, At4g17500 and Unigene17395_All (ERF1). During the early phase of both the lettuce/B. cinerea and the A. thaliana/B. cinerea interactions, genes in the ET pathway are heavily involved in the defence response, a finding which was replicated in the chrysanthemum/A. tenuissima interaction. The RNA-Seq data further show that a large number of genes are involved in the host response between six and 72 h post inoculation, including several genes involved in the JA and SA pathways. Some of these are discussed in more detail below.

Pathogen recognition-related genes modulated by A. tenuissima infection

Plant pattern recognition receptors (PRRs) perceive microbe-associated molecular patterns (MAMPs), a set of molecular signatures encompassing whole classes of microbes. This recognition initiates a basal level of immunity (termed MAMP-triggered immunity) [23]. Receptor-like kinases (RLKs), which form the largest plant receptor family, are PRRs localized at the plasma membrane [24]. The involvement of RLKs in certain host-pathogen interactions has been well documented experimentally [18]. The RNA-Seq-based transcriptomic analysis of the lettuce/B. cinerea interaction has revealed that several types of RLK are differentially transcribed [18], and the same phenomenon was recorded in the chrysanthemum/A. tenuissima interaction. The damage caused by microbes can induce the host to synthesize MAMP-like products, termed damage-associated molecular patterns (DAMPs) [23]. Oligogalacturonides (OGs) released from the plant cell wall activate the DAMP-associated response. The induction of an OG receptor wall-associated kinase (WAK) has been suggested as being necessary for the survival A. thaliana challenged by pathogen infection [25]. Five WAK-like kinases (encoded by Unigene52017_All, Unigene36001_All, Unigene55763_All, Unigene49198_All and Unigene12436_All) were among the DT genes detected in the B vs D contrast (Additional file 11: Table S10), and a sixth (Unigene49198_All) was identified in the A vs C contrast (Additional file 12: Table S11). Another study showed that WAKs were involved in the immune responses against B. cinerea. Furthermore, transgenic plants overexpressing WAK1 conferred resistance to B. cinerea in A. thaliana[26].

Genes encoding a second class of PRR, the leucine-rich repeat RLKs (LRR-RLK s), are known to be involved in both basal and cultivar-specific host immunity [27]. In A. thaliana, BRI1 encodes an LRR-RLK which forms a heterodimeric complex with a second LRR-RLK called BAK1 acting as a negative regulator of microbial infection-induced cell death [27]. Bak1 mutants develop spreading necrosis after the triggering of apoptosis by infection with B. cinerea[27]. The transcription of two LRR-RLK (Unigene20925_All and Unigene27322_All), four BRI-like (Unigene55939_All, Unigene15368_All, Unigene29292_All, Unigene18133_All), six BAK1 genes (Unigene37501_All, Unigene36228_All, Unigene16709_All, Unigene16958_All, Unigene14705_All, and Unigene27008_All), and two somatic embryogenesis receptor kinase (SERK) genes (Unigene22508_All, Unigene14416_All) was modulated by A. tenuissima infection in the B vs D contrast (Additional file 13: Table S12), as were an additional two BRI-like (Unigene15368_All and Unigene4877_All), a somatic embryogenesis receptor kinase (SERK) genes (Unigene22508_All), and four BAK1 genes (Unigene36228_All, Unigene4877_All, Unigene14705_All and Unigene27008_All) in the A vs C contrast (Additional file 12: Table S11). All the above genes were up-regulated except for Unigene18133_All. In the lettuce/B. cinerea interaction, one BRI-like (Lsa034184.1) gene was up-regulated, but no BAK1 genes [18].

The products of genes encoding cysteine-rich receptor-like kinases (CRKs) are RLKs which contain an extracellular cysteine-rich repeat domain. These genes are reportedly activated by oxidative stress, pathogen attack and exposure to SA [18, 28–30]. Although less widely researched, they have been characterized to be involved in the pathogen defense and programmed cell death in A. thaliana[28, 30]. Here, three CRKs w ere identified as DT genes (Unigene55939_All, Unigene15489_All and Unigene14705_All) in the B vs D comparison (Additional file 14: Table S13) and a further three (Unigene4299_All, Unigene14705_All and Unigene15489_All) in the A vs C comparison (Additional file 15: Table S14). Thus, PRR-RLKs are clearly involved in the chrysanthemum host response to A. tenuissima infection. The early activation of these genes may reflect their transcription as an attempt by the host to recognize MAMPs/DAMPs.

ROS (reactive oxygen species) detoxification genes modulated by A. tenuissima infection

MAMPs are known to trigger the production of ROS in response to pathogen infection, largely derived from NADPH oxidase activity (commonly referred to as “respiratory burst oxidase homologues” (rboh)) [17, 31]. Two chrysanthemum rboh homologues were recognized: Unigene300_All (homologue of rbohD) and Unigene45792_All (rbohF); both were differentially transcribed in the B vs D (Additional file 16: Table S15) and A vs C contrasts (Additional file 17: Table S16). The lettuce rbohD homologue (Lsa002796.1) was induced 48 h after infection with B. cinerea, but its rbohF homologue (Lsa018309.1) was not up-regulated [18]. In A. thaliana, both rbohD and rbohF are required for ROS detoxification [32]. Two further pathogen-inducible α-dioxygenases (Unigene32071_All and Unigene12359_All) (Additional file 16: Table S15) were differentially transcribed in the B vs D contrast, and one (Unigene32071_All) in A vs C (Additional file 17: Table S16); this class of gene was also up-regulated in the lettuce/B. cinerea system; its product is involved in protecting the cell against oxidative stress [18].

Genes associated with photosynthesis were mostly down-regulated by A. tenuissima infection

DT genes involved in photosynthesis were uniformly down-regulated in the B vs D contrast, with the sole exception of Unigene2020_All (Additional file 18: Table S17). In the contrast C vs D, four photosynthesis-related DT genes were detected, and all were down-regulated (Additional file 19: Table S18). The response mirrors the outcomes in the lettuce/B. cinerea and lettuce/Verticillium dahliae systems [18, 33], as well as in other documented plant-pathogen interactions [34–39]. Unigene6198 _All, predicted to be involved in the determination of the plant’s circadian clock, was down-regulated in both the B vs D (Additional file 18: Table S17) and C vs D contrasts (Additional file 19: Table S18). Similar examples of the suppression of such genes by pathogen infection have been described in both lettuce [18] and A. thaliana[40]. In the case of A. thaliana, B. cinerea infection appears to dampen the oscillating transcription of certain core clock components, leading to the suggestion that the pathogen attempts to dampen the host’s defence response, because a set of genes associated with plant immune responses was revealed to be regulated by the plant’s circadian clock [41].

Genes associated with cell wall protection affected by A. tenuissima infection

Both PGIP and PPO were up-regulated in response to A. tenuissima infection (Table 2). PGIP and PPO are known to respond to various cues [42–48], including the presence of A. solani[42], A. triticina[49], A. macrospora[50], Sclerospora graminicola[51] and Colletotrichum lindemuthianum[52]. PGIP transcript abundance increased over the period 6–24 h after inoculation with A. tenuissima (Table 2). Many fungi secrete endo-polygalacturonase, an enzyme which degrades the polysaccharides present in the plant cell wall. Host genes involved in the early defence response include those which help to reinforce the cell wall, and thereby inhibit pathogen entry. PGIP’s role in mediating resistance to A. alternata infection has been shown in both rough lemon and Japanese pear [9, 53], while the growth of B. cinerea is restricted in transgenic tomato heterologously expressing pear PGIP[54]. PPO transcription in chrysanthemum was enhanced after inoculation (Table 2), similar to what has been observed in the leaf of Japanese pear inoculated with A. alternata[9]. In tomato, the constitutive expression of PPO increases host resistance to Pseudomonas syringae, while its down-regulation enhances susceptibility [55, 56].

JA and SA signalling pathway-related genes involved in the response to A. tenuissima infection

CmJAZ (Unigene11800 _All), CmMYC2 (Unigene3689 _All), CmVSP (CL10952), CmNPR1 (Unigene23699_All) and CmTGA (Unigene52251_All) were all induced by A. tenuissima infection (Table 2). NPR1 is a major component of SA signalling, functioning as a co-activator of the TGA transcription factors known to regulate the transcription of various SA-responsive genes [57]. JAZ (jasmonate ZIM domain) proteins repress JA signalling by binding to transcriptional regulators such as MYC2. The degradation of JAZ proteins relieves the JAZ-mediated repression of the JA signalling pathway and thereby activates a large number of JA responsive genes [57, 58]. JA and SA both play an important part in host defence against herbivorous insects and microbial pathogens. In A. thaliana, VSP2 and PR-1 are, respectively, JA- and SA-responsive proteins [58]. Overall, a number of genes within the JA and SA signalling pathways are involved in the defence response of chrysanthemum against A. tenuissima infection.

Transcription factors responding to A. tenuissima infection

Transcription factors are central to the control of the timing and placement of defence response gene expression [59]. Their mode of action is to first recognize and then bind to regulatory elements located in the promoter region of their target genes, thereby activating or de-activating their transcription. Here, five classes of transcription factor (MYB, AP2/ERF, WRKY, NAC and GRAS) were identified among the DT genes responding to A. tenuissima infection. In A. thaliana, AtMYB30, AtMYB44 and AtMYB96 are all involved in the triggering of apoptosis and therefore resistance against biotrophic bacterial pathogens such as Pseudomonas syringae[60–64]. AtMYB108 is required for resistance against B. cinerea and A. brassicicola[65]. The transcription of AtMYB58 has been associated with secondary cell wall formation [66]. Here, one CmMYB copy (Unigene6575_All and a homologue of AtMYB58), was prominently transcribed from 48 h after inoculation with A. tenuissima (Table 2). Therefore, it was speculated that CmMYB (Unigene6575_All) may be involved in defence response to the necrotrophic fungus A. tenuissima through the regulation of secondary cell wall biosynthesis. Modifications to the plant cell wall were already recognized as the potential mechanism of resistance, which was previously provided by the reports on the response of plants to fungal challenges [67]. A second CmMYB (Unigene16113_All, homologous to AtMYB74) was also up-regulated by A. tenuissima infection (Table 2). In A. thaliana, AtMYB74 reportedly responds to salinity stress and the exogenous supply of abscisic acid, ETH and JA [68], but not as yet to pathogen infection. Therefore, it was speculated that this gene may be involved in defense response to the necrotrophic fungus A. tenuissima and may be also associated with ABA and JA signaling in Chrysanthemum. Other transcription factors in the AP2/ERF, WRKY, NAC and GRAS families were abundantly transcribed in the leaf following A. tenuissima inoculation. Many reports have also indicated that four families of transcription factors: MYB proteins, ethylene-responsive-element-binding factors (ERF), WRKY proteins and NAC proteins link to plant stress responses, such as pathogens [59, 60, 69–72]. Little evidence has been provided to date regarding the participation of GRAS transcription factors in the defence response.

In conclusion, it was clear that infection with A. tenuissima induced a wide range of genes in the chrysanthemum leaf. The response involved a complex set of interactions between pathogenesis-related genes, genes in the JA and SA signalling pathway and transcription factors. A more detailed understanding of the identity of these genes will help to unravel the molecular basis of the defence response of chrysanthemum to A. tenuissima infection, and eventually lead to the recognition of candidates for the targeted genetic improvement of chrysanthemum.

Conclusions

In this study, we characterized the leaf transcriptome of chrysanthemum and provided the comparative DT genes involeved in the interaction between chrysanthemum and A. tenuissima. These findings provide a substantial contribution to existing sequence resources of chrysanthemum, and a strong basis for further characterization of gene expression profiles in the interaction of chrysanthemum and A. tenuissima. The majority of the DT genes were those involved in pathogen recognition, reactive oxygen species detoxification, cell wall modification, phytohormone signalling, and transcription factors belonging to various families were also identified, which will improve our understanding of the molecular mechanisms underlining direct response and induced systemic resistance of chrysanthemum to A. tenuissim.

Methods

Plant materials, A. tenuissima inoculum preparation and inoculation

The chrysanthemum variety ‘Zaoyihong’ was obtained from the Chrysanthemum Germplasm Resource Preserving Centre, Nanjing Agricultural University, China. Uniform cuttings were propagated in sand, and rooted seedlings transplanted into a 2:1 mixture of garden soil and vermiculite without fertilizer supplementation. The plants were grown under a 16 h photoperiod with a day/night temperature of, respectively, 25°C and 18°C. The relative humidity was maintained at 68-75% [3]. A. tenuissima conidia were isolated from diseased chrysanthemum plants, and cultured on potato dextrose agar at 25°C. An aqueous suspension of 106 spores per ml was prepared with a few drops of Triton X-100 added as a wetting agent [3]. The surface of the second and third true leaves of 20 day old root cuttings was punctured with a needle (approximately 0.30 mm diameter), and a 10 μl droplet of spore suspension was placed on the puncture site. Mock treatments comprised 10 μl droplets of sterile distilled water. After inoculation, the plants were held at 100% relative humidity and 25°C in the dark for 24 h, and then illuminated with 120 μmol m-2 s-1 cool white fluorescent light with a 12 h photoperiod. The leaves of three seedlings were sampled for each treatment at 0 h, 6 h, 24 h, 48 h and 72 h after inoculation. The samples collected at defined time points of each treatment were pooled for RNA-seq.

RNA extraction

Four separate libraries (A-D, see Figure 1) were prepared. Extracts of the second and the third true leaves of mock-treated and pathogen-infected plants gave rise to, respectively, libraries B and D, while extracts of the first and the fourth noninfected (systemic) true leaves of plants gave rise to libraries A and C. Total RNA was isolated using a Total RNA Isolation System (Takara, Japan), following the manufacturer’s recommendations. The quality of the total RNA (RNA Integrity Number > 6.5 and 28S:18S > 1.0) was verified using a 2100 Bioanalyzer RNA Nano chip device (Agilent, Santa Clara, CA, USA) and its concentration ascertained using an ND-1000 spectrophotometer (NanoDrop, Wilmington, DE). The standards applied were 1.8 ≤ OD260/280 ≤ 2.2 and OD260/230 ≥ 1.8. At least 10 μg RNA was pooled in an equimolar fashion from each of the three sampled plants [73].

cDNA library construction and Illumina sequencing

Each total RNA extract was first treated with RNase-free DNase I (TaKaRa, Dalian, China) to remove contaminating DNA, and the mRNA content was concentrated by capturing on magnetic oligo (dT) beads. The mRNA was fragmented to a size of ~200 bp using a fragmentation buffer, and the resulting fragments used to synthesize the first cDNA strand by priming with random hexamers. The second strand was generated using a SuperScript Double-Stranded cDNA Synthesis kit (Invitrogen, Camarillo, CA), purified via magnetic beads, the ends repaired and a single adenine base added to the 3′ ends. Sequencing adaptors were then ligated to the fragments, and agarose gel electrophoresis used to select the range of fragments suitable for PCR amplification. Sequencing using an Illumina HiSeq™ 2000 platform was performed at the Beijing Genomics Institute (Shenzhen, China; http://www.genomics.cn/index.php), following the manufacturer’s protocols.

Treatment of sequence data

Raw reads were saved as .fastq files, and filtered to remove adaptor sequences, reads in which the proportion of non-called bases was > 10% and reads in which low quality (≤ 5) bases represented > 50% of the read. The remaining reads were mapped onto the set of chrysanthemum unigene sequences using SOAPaligner/SOAP2 [74]. A maximum of two mismatches was permitted for the purpose of alignment. The frequency of occurrence of individual reads was normalized to RPKM (reads per kb per million reads) [19]. Differential transcription between pathogen-inoculated and mock samples was based on the log2 ratio of the two RPKM values. All raw RNA-Seq data have been deposited at the sequence read archive (SRA) of NCBI (Additional file 20: Table S19).

Identification of DT genes

Following Audic and Claverie (1997) [20], a stringent algorithm was developed to identify DT genes. The FDR (false discovery rate) provides a criterion to determine the P-value threshold in multiple tests and analyses by manipulating the FDR value. Here, differential transcription was declared provided that the P-value was < 0.05, the FDR ≤ 0.001 and the absolute value of log2 induction ratios of treated samples compared with mock-treatment or control (CK) was ≥ 1.0. Standard gene ontology (GO) was used to describe DT gene functionality, and a hypergeometric test was used to map the DT genes to GO terms based on the BGI WEGO (Web Gene Ontology Annotation Plot, http://wego.genomics.org.cn/cgi-bin/wego/index.pl).

Quantitative real-time PCR (qPCR) validation

The transcription of 23 selected candidate genes was determined using quantitative real time PCR (qPCR). The samples collected at different time points were pooled, which were used for qPCR analysis, contaminating DNA removed by RNase-free DNase I treatment and the first cDNA strand synthesized using a Super RT kit (BioTeke, Beijing, China). A set of gene-specific primer pairs (sequences given in Table 3) was designed using Primer3 software [18, 75]. qPCRs were based on SYBR_Green I (TOYOBO’, Japan) implemented in a Rotor-Gene 3000 device (Corbett, Australia). The chrysanthemum EF1α gene was used as a reference. Each 25 μl qPCR reaction contained 10 μl SYBR Green PCR master mix, 0.2 μM of each primer and 10 ng cDNA, and the amplification regime consisted of an initial denaturation of 95°C/60 s, followed by 40 cycles of 95°C/15 s, 55°C/15 s, 72°C/20 s. Transcript abundances are given as the mean ± SE of three replicates. Relative transcription levels were calculated using the 2-△△CT method [76].

Table 3 Primers of quantitative reverse transcription-polymerase chain reaction for validation of the digital gene expression data
Full size table

References

  1. Shibata M: Importance of genetic transformation in ornamental plant breeding. Plant Biotechnol. 2008, 25 (1): 3-8. 10.5511/plantbiotechnology.25.3.

    Article  Google Scholar 

  2. Teixeira da Silva JA, Shinoyama H, Aida R, Matsushita Y, Raj SK, Chen F: Chrysanthemum Biotechnology: Quo vadis?. Crit Rev Plant Sci. 2013, 32 (1): 21-52. 10.1080/07352689.2012.696461.

    Article  CAS  Google Scholar 

  3. Xu G, Liu Y, Chen S, Chen F: Potential structural and biochemical mechanisms of compositae wild species resistance to Alternaria tenuissima. Russ J Plant Physl+. 2011, 58 (3): 491-497. 10.1134/S1021443711030216.

    Article  CAS  Google Scholar 

  4. Deng Y, Chen S, Chang Q, Wang H, Chen F: The chrysanthemum × Artemisia vulgaris intergeneric hybrid has better rooting ability and higher resistance to alternaria leaf spot than its chrysanthemum parent. Sci Hortic-Amsterdam. 2012, 134: 185-190.

    Article  Google Scholar 

  5. Mirkova E, Konstantinova P: First report of Alternaria leaf spot on gerbera (Gerbera jamesonii H. Bolux ex JD Hook) in Bulgaria. J Phytopathol. 2003, 15 (6): 323-328.

    Article  Google Scholar 

  6. Mukherjee AK, Carp MJ, Zuchman R, Ziv T, Horwitz BA, Gepstein S: Proteornics of the response of Arabidopsis thaliana to infection with Alternaria brassicicola. J Proteomics. 2010, 73 (4): 709-720. 10.1016/j.jprot.2009.10.005.

    Article  CAS  PubMed  Google Scholar 

  7. Jia C, Zhang L, Liu L, Wang J, Li C, Wang Q: Multiple phytohormone signalling pathways modulate susceptibility of tomato plants to Alternaria alternata f. sp. lycopersici. J Exp Bot. 2013, 64 (2): 637-650. 10.1093/jxb/ers360.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  8. Sinha R, Chattopadhyay S: Changes in the leaf Proteome profile of Mentha arvensis in response to Alternaria alternata infection. J Proteomics. 2011, 74 (3): 327-336. 10.1016/j.jprot.2010.11.009.

    Article  CAS  PubMed  Google Scholar 

  9. Egusa M, Ochi H, Tsuge T, Otani H, Kodama M: Identification of putative defense-related genes in Japanese pear against Alternaria alternata infection using suppression subtractive hybridization and expression analysis. J Gen Plant Pathol. 2009, 75 (2): 119-124. 10.1007/s10327-009-0150-8.

    Article  CAS  Google Scholar 

  10. Kulye M, Liu H, Zhang Y, Zeng H, Yang X, Qiu D: Hrip1, a novel protein elicitor from necrotrophic fungus, Alternaria tenuissima, elicits cell death, expression of defence-related genes and systemic acquired resistance in tobacco. Plant Cell Environ. 2012, 35 (12): 2104-2120. 10.1111/j.1365-3040.2012.02539.x.

    Article  CAS  PubMed  Google Scholar 

  11. Yano A, Suzuki K, Uchimiya H, Shinshi H: Induction of hypersensitive cell death by a fungal protein in cultures of tobacco cells. Mol Plant Microbe In. 1998, 11 (2): 115-123. 10.1094/MPMI.1998.11.2.115.

    Article  CAS  Google Scholar 

  12. Garcia-Brugger A, Lamotte O, Vandelle E, Bourque S, Lecourieux D, Poinssot B, Wendehenne D, Pugin A: Early signaling events induced by elicitors of plant defenses. Mol Plant Microbe In. 2006, 19 (7): 711-724. 10.1094/MPMI-19-0711.

    Article  CAS  Google Scholar 

  13. Iandolino A, Nobuta K, Goes da Silva F, Cook DR, Meyers BC: Comparative expression profiling in grape (Vitis vinifera) berries derived from frequency analysis of ESTs and MPSS signatures. BMC Plant Biol. 2008, 8 (1): 1-16. 10.1186/1471-2229-8-1.

    Article  Google Scholar 

  14. Vizoso P, Meisel LA, Tittarelli A, Latorre M, Saba J, Caroca R, Maldonado J, Cambiazo V, Campos-Vargas R, Gonzalez M: Comparative EST transcript profiling of peach fruits under different post-harvest conditions reveals candidate genes associated with peach fruit quality. BMC Genomics. 2009, 10: 423-440. 10.1186/1471-2164-10-423.

    Article  PubMed Central  PubMed  Google Scholar 

  15. Wang QQ, Liu F, Chen XS, Ma XJ, Zeng HQ, Yang ZM: Transcriptome profiling of early developing cotton fiber by deep-sequencing reveals significantly differential expression of genes in a fuzzless/lintless mutant. Genomics. 2010, 96 (6): 369-376. 10.1016/j.ygeno.2010.08.009.

    Article  CAS  PubMed  Google Scholar 

  16. Xu L, Zhu L, Tu L, Liu L, Yuan D, Jin L, Long L, Zhang X: Lignin metabolism has a central role in the resistance of cotton to the wilt fungus Verticillium dahliae as revealed by RNA-Seq-dependent transcriptional analysis and histochemistry. J Exp Bot. 2011, 62 (15): 5607-5621. 10.1093/jxb/err245.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  17. Li CY, Deng GM, Yang J, Viljoen A, Jin Y, Kuang RB, Zuo CW, Lv ZC, Yang QS, Sheng O, et al: Transcriptome profiling of resistant and susceptible Cavendish banana roots following inoculation with Fusarium oxysporum f. sp. cubense tropical race 4. BMC Genomics. 2012, 13: 374-384. 10.1186/1471-2164-13-374.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  18. de Cremer K, Mathys J, Vos C, Froenicke L, Michelmore RW, Cammue B, de Coninck B: RNAseq-based transcriptome analysis of Lactuca sativa infected by the fungal necrotroph Botrytis cinerea. Plant Cell Environ. 2013, 36: 1992-2007.

    CAS  PubMed  Google Scholar 

  19. Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5 (7): 621-628. 10.1038/nmeth.1226.

    Article  CAS  PubMed  Google Scholar 

  20. Audic S, Claverie JM: The significance of digital gene expression profiles. Genome Res. 1997, 7 (10): 986-995.

    CAS  PubMed  Google Scholar 

  21. Ferrari S, Galletti R, Denoux C, de Lorenzo G, Ausubel FM, Dewdney J: Resistance to Botrytis cinerea induced in Arabidopsis by elicitors is independent of salicylic acid, ethylene, or jasmonate signaling but requires PHYTOALEXIN DEFICIENT3. Plant Physiol. 2007, 144 (1): 367-379. 10.1104/pp.107.095596.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  22. Mulema JM, Denby KJ: Spatial and temporal transcriptomic analysis of the Arabidopsis thaliana–Botrytis cinerea interaction. Mol Biol Rep. 2012, 39 (4): 4039-4049. 10.1007/s11033-011-1185-4.

    Article  CAS  PubMed  Google Scholar 

  23. Boller T, Felix G: A Renaissance of Elicitors: Perception of Microbe-Associated Molecular Patterns and Danger Signals by Pattern-Recognition Receptors. Annu Rev Plant Biol. 2009, 60: 379-406. 10.1146/annurev.arplant.57.032905.105346.

    Article  CAS  PubMed  Google Scholar 

  24. Toer M, Lotze MT, Holton N: Receptor-mediated signalling in plants: molecular patterns and programmes. J Exp Bot. 2009, 60 (13): 3645-3654. 10.1093/jxb/erp233.

    Article  CAS  Google Scholar 

  25. He ZH, He DZ, Kohorn BD: Requirement for the induced expression of a cell wall associated receptor kinase for survival during the pathogen response. Plant J. 1998, 14 (1): 55-63. 10.1046/j.1365-313X.1998.00092.x.

    Article  CAS  PubMed  Google Scholar 

  26. Brutus A, Sicilia F, Macone A, Cervone F, de Lorenzo G: A domain swap approach reveals a role of the plant wall-associated kinase 1 (WAK1) as a receptor of oligogalacturonides. Proc Natl Acad Sci U S A. 2010, 107 (20): 9452-9457. 10.1073/pnas.1000675107.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  27. Kemmerling B, Schwedt A, Rodriguez P, Mazzotta S, Frank M, Qamar SA, Mengiste T, Betsuyaku S, Parker JE, Müssig C: The BRI1-associated kinase 1, BAK1, has a brassinolide-independent role in plant cell-death control. Curr Biol. 2007, 17 (13): 1116-1122. 10.1016/j.cub.2007.05.046.

    Article  CAS  PubMed  Google Scholar 

  28. Czernic P, Visser B, Sun WN, Savoure A, Deslandes L, Marco Y, van Montagu M, Verbruggen N: Characterization of an Arabidopsis thaliana receptor-like protein kinase gene activated by oxidative stress and pathogen attack. Plant J. 1999, 18 (3): 321-327. 10.1046/j.1365-313X.1999.00447.x.

    Article  CAS  PubMed  Google Scholar 

  29. Chen KG, Du LQ, Chen ZX: Sensitization of defense responses and activation of programmed cell death by a pathogen-induced receptor-like protein kinase in Arabidopsis. Plant Mol Biol. 2003, 53 (1): 61-74.

    Article  CAS  PubMed  Google Scholar 

  30. Chen KG, Fan BF, Du LQ, Chen ZX: Activation of hypersensitive cell death by pathogen-induced receptor-like protein kinases from Arabidopsis. Plant Mol Biol. 2004, 56 (2): 271-283. 10.1007/s11103-004-3381-2.

    Article  CAS  PubMed  Google Scholar 

  31. Torres MA, Jones JD, Dangl JL: Reactive oxygen species signaling in response to pathogens. Plant Physiol. 2006, 141 (2): 373-378. 10.1104/pp.106.079467.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  32. Torres MA, Dangl JL, Jones JD: Arabidopsis gp91phox homologues AtrbohD and AtrbohF are required for accumulation of reactive oxygen intermediates in the plant defense response. Proc Natl Acad Sci USA. 2002, 99 (1): 517-522. 10.1073/pnas.012452499.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  33. Klosterman SJ, Anchieta A, Garcia-Pedrajas MD, Maruthachalam K, Hayes RJ, Subbarao KV: SSH reveals a linkage between a senescence-associated protease and Verticillium wilt symptom development in lettuce (Lactuca sativa). Physiol Mol Plant P. 2011, 76 (1): 48-58. 10.1016/j.pmpp.2011.05.003.

    Article  CAS  Google Scholar 

  34. Milli A, Cecconi D, Bortesi L, Persi A, Rinalducci S, Zamboni A, Zoccatelli G, Lovato A, Zolla L, Polverari A: Proteomic analysis of the compatible interaction between Vitis vinifera and Plasmopara viticola. J Proteomics. 2012, 75 (4): 1284-1302. 10.1016/j.jprot.2011.11.006.

    Article  CAS  PubMed  Google Scholar 

  35. Gyetvai G, Sonderkaer M, Goebel U, Basekow R, Ballvora A, Imhoff M, Kersten B, Nielsen KL, Gebhardt C: The transcriptome of compatible and incompatible interactions of potato (Solanum tuberosum) with Phytophthora infestans revealed by DeepSAGE analysis. PLoS One. 2012, 7 (2): e31526-10.1371/journal.pone.0031526.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  36. Bonfig KB, Gabler A, Simon UK, Luschin-Ebengreuth N, Hatz M, Berger S, Muhammad N, Zeier J, Sinha AK, Roitsch T: Post-translational derepression of invertase activity in source leaves via down-regulation of invertase inhibitor expression is part of the plant defense response. Mol Plant. 2010, 3 (6): 1037-1048. 10.1093/mp/ssq053.

    Article  CAS  PubMed  Google Scholar 

  37. Swarbrick PJ, Schulze-Lefert P, Scholes JD: Metabolic consequences of susceptibility and resistance (race-specific and broad-spectrum) in barley leaves challenged with powdery mildew. Plant Cell Environ. 2006, 29 (6): 1061-1076. 10.1111/j.1365-3040.2005.01472.x.

    Article  CAS  PubMed  Google Scholar 

  38. Scharte J, Schon H, Weis E: Photosynthesis and carbohydrate metabolism in tobacco leaves during an incompatible interaction with Phytophthora nicotianae. Plant Cell Environ. 2005, 28 (11): 1421-1435. 10.1111/j.1365-3040.2005.01380.x.

    Article  CAS  Google Scholar 

  39. Chou H-M, Bundock N, Rolfe SA, Scholes JD: Infection of Arabidopsis thaliana leaves with Albugo candida (white blister rust) causes a reprogramming of host metabolism. Mol Plant Pathol. 2000, 1 (2): 99-113. 10.1046/j.1364-3703.2000.00013.x.

    Article  CAS  PubMed  Google Scholar 

  40. Windram O, Madhou P, McHattie S, Hill C, Hickman R, Cooke E, Jenkins DJ, Penfold CA, Baxter L, Breeze E, et al: Arabidopsis defense against Botrytis cinerea: chronology and regulation deciphered by high-resolution temporal transcriptomic analysis. Plant Cell. 2012, 24 (9): 3530-3557. 10.1105/tpc.112.102046.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  41. Wang W, Barnaby JY, Tada Y, Li H, Tör M, Caldelari D, du Lee FXD, Dong X: Timing of plant immune responses by a central circadian regulator. Nature. 2011, 470 (7332): 110-114. 10.1038/nature09766.

    Article  CAS  PubMed  Google Scholar 

  42. Thipyapong P, Steffens JC: Tomato polyphenol oxidase (differential response of the polyphenol oxidase F promoter to injuries and wound signals). Plant Physiol. 1997, 115 (2): 409-418.

    CAS  PubMed Central  PubMed  Google Scholar 

  43. Ferrari S, Sella L, Janni M, de Lorenzo G, Favaron F, D’Ovidio R: Transgenic expression of polygalacturonase-inhibiting proteins in Arabidopsis and wheat increases resistance to the flower pathogen Fusarium graminearum. Plant Biol. 2012, 14: 31-38.

    Article  CAS  PubMed  Google Scholar 

  44. Shetty SM, Chandrashekar A, Venkatesh YP: Promoter analyses and transcriptional profiling of eggplant polyphenol oxidase 1 gene (SmePPO1) reveal differential response to exogenous methyl jasmonate and salicylic acid. J Plant Physiol. 2012, 169: 718-730. 10.1016/j.jplph.2012.01.007.

    Article  CAS  PubMed  Google Scholar 

  45. Bhonwong A, Stout MJ, Attajarusit J, Tantasawat P: Defensive role of tomato polyphenol oxidases against cotton bollworm (Helicoverpa armigera) and beet armyworm (Spodoptera exigua). J Chem Ecol. 2009, 35 (1): 28-38. 10.1007/s10886-008-9571-7.

    Article  CAS  PubMed  Google Scholar 

  46. Mahanil S, Attajarusit J, Stout MJ, Thipyapong P: Overexpression of tomato polyphenol oxidase increases resistance to common cutworm. Plant Sci. 2008, 174 (4): 456-466. 10.1016/j.plantsci.2008.01.006.

    Article  CAS  Google Scholar 

  47. Wang J, Constabel CP: Polyphenol oxidase overexpression in transgenic Populus enhances resistance to herbivory by forest tent caterpillar (Malacosoma disstria). Planta. 2004, 220 (1): 87-96. 10.1007/s00425-004-1327-1.

    Article  CAS  PubMed  Google Scholar 

  48. Thipyapong P, Hunt MD, Steffens JC: Systemic wound induction of potato (Solanum tuberosum) polyphenol oxidase. Phytochemistry. 1995, 40 (3): 673-676. 10.1016/0031-9422(95)00359-F.

    Article  CAS  Google Scholar 

  49. Tyagi M, Kayastha AM, Sinha B: The role of peroxidase and polyphenol oxidase isozymes in wheat resistance to Alternaria triticina. Biol Plantarum. 2000, 43 (4): 559-562. 10.1023/A:1002821715719.

    Article  CAS  Google Scholar 

  50. Bashan Y: Phenols in cotton seedlings resistant and susceptible to Alternaria macrospora. J Phytopathol. 1986, 116: 1-10. 10.1111/j.1439-0434.1986.tb00888.x.

    Article  CAS  Google Scholar 

  51. Niranjan RAJS, Sarosh B, Shetty H: Induction and accumulation of polyphenol oxidase activities as implicated in development of resistance against pearl millet downy mildew disease. Funct Plant Biol. 2006, 33 (6): 563-571. 10.1071/FP06003.

    Article  Google Scholar 

  52. Campos A, Ferreira A, Hampe M, Antunes I, Brancao N, Silveira E, Osorio V, Augustin E: Peroxidases and polyphenol oxidase activity in bean anthracnose resistance. Pesqui Agropecu Bras. 2004, 39: 637-643. 10.1590/S0100-204X2004000700004.

    Article  Google Scholar 

  53. Gotoh Y, Nalumpang S, Isshiki A, Utsumi T, Gomi K, Yamamoto H, Akimitsu K: A cDNA encoding polygalacturonase-inhibiting protein induced in citrus leaves by polygalacturonase of Alternaria citri. J Gen Plant Pathol. 2002, 68 (1): 57-61. 10.1007/PL00013054.

    Article  CAS  Google Scholar 

  54. Powell AL, van Kan J, ten Have A, Visser J, Greve LC, Bennett AB, Labavitch JM: Transgenic expression of pear PGIP in tomato limits fungal colonization. Mol Plant Microbe In. 2000, 13 (9): 942-950. 10.1094/MPMI.2000.13.9.942.

    Article  CAS  Google Scholar 

  55. Li L, Steffens JC: Overexpression of polyphenol oxidase in transgenic tomato plants results in enhanced bacterial disease resistance. Planta. 2002, 215 (2): 239-247. 10.1007/s00425-002-0750-4.

    Article  CAS  PubMed  Google Scholar 

  56. Thipyapong P, Hunt MD, Steffens JC: Antisense downregulation of polyphenol oxidase results in enhanced disease susceptibility. Planta. 2004, 220 (1): 105-117. 10.1007/s00425-004-1330-6.

    Article  CAS  PubMed  Google Scholar 

  57. Ballaré CL: Jasmonate-induced defenses: a tale of intelligence, collaborators and rascals. Trends Plant Sci. 2011, 16 (5): 249-257. 10.1016/j.tplants.2010.12.001.

    Article  PubMed  Google Scholar 

  58. Leon-Reyes A, van der Does D, de lange ES, Delker C, Wasternack C, van Wees CSM, Ritsema T, Pieterse CMJ: Salicylate-mediated suppression of jasmonate-responsive gene expression in Arabidopsis is targeted downstream of the jasmonate biosynthesis pathway. Planta. 2010, 232 (6): 1423-1432. 10.1007/s00425-010-1265-z.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  59. Singh KB, Foley RC, Oñate-Sánchez L: Transcription factors in plant defense and stress responses. Curr Opin Plant Biol. 2002, 5 (5): 430-436. 10.1016/S1369-5266(02)00289-3.

    Article  CAS  PubMed  Google Scholar 

  60. Dubos C, Stracke R, Grotewold E, Weisshaar B, Martin C, Lepiniec L: MYB transcription factors in Arabidopsis. Trends Plant Sci. 2010, 15 (10): 573-581. 10.1016/j.tplants.2010.06.005.

    Article  CAS  PubMed  Google Scholar 

  61. Vailleau F, Daniel X, Tronchet M, Montillet J-L, Triantaphylidès C, Roby D: A R2R3-MYB gene, AtMYB30, acts as a positive regulator of the hypersensitive cell death program in plants in response to pathogen attack. Proc Natl Acad Sci U S A. 2002, 99 (15): 10179-10184. 10.1073/pnas.152047199.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  62. Shim JS, Jung C, Lee S, Min K, Lee YW, Choi Y, Lee JS, Song JT, Kim JK, Choi YD: AtMYB44 regulates WRKY70 expression and modulates antagonistic interaction between salicylic acid and jasmonic acid signaling. Plant J. 2013, 73 (3): 483-495. 10.1111/tpj.12051.

    Article  CAS  PubMed  Google Scholar 

  63. Zou B, Jia Z, Tian S, Wang X, Gou Z, Lü B, Dong H: AtMYB44 positively modulates disease resistance to Pseudomonas syringae through the salicylic acid signalling pathway in Arabidopsis. Funct Plant Biol. 2012, 40 (3): 304-313.

    Article  Google Scholar 

  64. Seo PJ, Park CM: MYB96-mediated abscisic acid signals induce pathogen resistance response by promoting salicylic acid biosynthesis in Arabidopsis. New Phytol. 2010, 186 (2): 471-483. 10.1111/j.1469-8137.2010.03183.x.

    Article  CAS  PubMed  Google Scholar 

  65. Mengiste T, Chen X, Salmeron J, Dietrich R: The BOTRYTIS SUSCEPTIBLE1 gene encodes an R2R3MYB transcription factor protein that is required for biotic and abiotic stress responses in Arabidopsis. Plant Cell. 2003, 15 (11): 2551-2565. 10.1105/tpc.014167.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  66. Zhou J, Lee C, Zhong R, Ye ZH: MYB58 and MYB63 are transcriptional activators of the lignin biosynthetic pathway during secondary cell wall formation in Arabidopsis. Plant Cell. 2009, 21 (1): 248-266. 10.1105/tpc.108.063321.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  67. Ramírez V, Agorio A, Coego A, García-Andrade J, Hernández MJ, Balaguer B, Ouwerkerk PB, Zarra I, Vera P: MYB46 modulates disease susceptibility to Botrytis cinerea in Arabidopsis. Plant Physiol. 2011, 155 (4): 1920-1935. 10.1104/pp.110.171843.

    Article  PubMed Central  PubMed  Google Scholar 

  68. Yanhui C, Xiaoyuan Y, Kun H, Meihua L, Jigang L, Zhaofeng G, Zhiqiang L, Yunfei Z, Xiaoxiao W, Xiaoming Q: The MYB transcription factor superfamily of Arabidopsis: expression analysis and phylogenetic comparison with the rice MYB family. Plant Mol Biol. 2006, 60 (1): 107-124. 10.1007/s11103-005-2910-y.

    Article  PubMed  Google Scholar 

  69. Du H, Yang SS, Feng BR, Tang YX, Huang YB, Liang Z: Genome-wide analysis of the MYB transcription factor superfamily in soybean. BMC Plant Biol. 2012, 12 (1): 106-127. 10.1186/1471-2229-12-106.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  70. Zheng Z, Qamar SA, Chen Z, Mengiste T: Arabidopsis WRKY33 transcription factor is required for resistance to necrotrophic fungal pathogens. Plant J. 2006, 48 (4): 592-605. 10.1111/j.1365-313X.2006.02901.x.

    Article  CAS  PubMed  Google Scholar 

  71. Lin R, Zhao W, Meng X, Wang M, Peng Y: Rice gene OsNAC19 encodes a novel NAC-domain transcription factor and responds to infection by Magnaporthe grisea. Plant Sci. 2007, 172 (1): 120-130. 10.1016/j.plantsci.2006.07.019.

    Article  CAS  Google Scholar 

  72. Xia N, Zhang G, Liu XY, Deng L, Cai GL, Zhang Y, Wang XJ, Zhao J, Huang LL, Kang ZS: Characterization of a novel wheat NAC transcription factor gene involved in defense response against stripe rust pathogen infection and abiotic stresses. Mol Biol Rep. 2010, 37 (8): 3703-3712. 10.1007/s11033-010-0023-4.

    Article  CAS  PubMed  Google Scholar 

  73. Wang H, Jiang J, Chen S, Qi X, Peng H, Li P, Song A, Guan Z, Fang W, Liao Y, et al: Next-generation sequencing of the chrysanthemum nankingense (Asteraceae) transcriptome permits large-scale unigene assembly and SSR marker discovery. PLoS One. 2013, 8 (4): e62293-10.1371/journal.pone.0062293.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

  74. Li R, Yu C, Li Y, Lam TW, Yiu SM, Kristiansen K, Wang J: SOAP2: an improved ultrafast tool for short read alignment. Bioinformatics. 2009, 25 (15): 1966-1967. 10.1093/bioinformatics/btp336.

    Article  CAS  PubMed  Google Scholar 

  75. Rozen S, Skaletsky H: Primer3 on the WWW for general users and for biologist programmers. Methods Mol Biol. 2000, 132: 365-386.

    CAS  PubMed  Google Scholar 

  76. Song A, Lou W, Jiang J, Chen S, Sun Z, Guan Z, Fang W, Teng N, Chen F: An isoform of eukaryotic initiation factor 4E from Chrysanthemum morifolium Interacts with chrysanthemum virus B coat protein. PLoS One. 2013, 8 (3): e57229-10.1371/journal.pone.0057229.

    Article  CAS  PubMed Central  PubMed  Google Scholar 

Download references

Acknowledgment

This work was supported by funding from the Program for New Century Excellent Talents in University of Chinese Ministry of Education (NCET-10-0492, NCET-12-0890, NCET-11-0669), the National Natural Science Foundation of China (31272196), the Natural Science Fund of Jiangsu Province (BK201164, BK2012773), 948 Project of Ministry of Agriculture (2011-G17, 2013-S13), Fund for Independent Innovation of Agricultural Sciences in Jiangsu Province [CX(12)2020], the Program for Hi-Tech Research, Jiangsu, China, (BE2012350,BE2011325), and the Youth Science and Technology Innovation Fund from Nanjing Agricultural University (KJ2011009).

Author information

Authors and Affiliations

  1. College of Horticulture, Nanjing Agricultural University, Nanjing, 210095, China

    Huiyun Li, Sumei Chen, Aiping Song, Haibin Wang, Weimin Fang, Zhiyong Guan, Jiafu Jiang & Fadi Chen

  2. Jiangsu Province Engineering Lab for Modern Facility Agriculture Technology and Equipment, Nanjing, 210095, China

    Huiyun Li & Fadi Chen

Authors
  1. Huiyun Li
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Sumei Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Aiping Song
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Haibin Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Weimin Fang
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Zhiyong Guan
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Jiafu Jiang
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Fadi Chen
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding authors

Correspondence to Jiafu Jiang or Fadi Chen.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

LHY and SAP contributed to RACE PCR, real-time PCR, bioinformatics analysis and writing of the manuscript. CFD, JJF and CSM conceived of the study, and participated in its design and contributed to revisions of the manuscript. FWM and GZY participated in experiment materials preparation. WHB helped with the RNA extraction. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2013_5635_MOESM1_ESM.tiff

Additional file 1: Figure S1: Composition of the raw reads in the four RNA libraries. “Clean” reads are those remaining after removal of adaptor sequences, reads in which the proportion of non-called bases was >10% and reads in which low quality (≤ 5) bases represented >50% of the reads. The numbers in parentheses indicate the percentage of each type of read present. (TIFF 79 KB)

12864_2013_5635_MOESM2_ESM.xls

Additional file 2: Table S1: GO classification of the genes differentially transcribed in the contrast between libraries A and B. (XLS 33 KB)

12864_2013_5635_MOESM3_ESM.xls

Additional file 3: Table S2: GO classification of the genes differentially transcribed in the contrast between libraries A and C. (XLS 52 KB)

12864_2013_5635_MOESM4_ESM.xls

Additional file 4: Table S3: GO classification of the genes differentially transcribed in the contrast between libraries B and D. (XLS 94 KB)

12864_2013_5635_MOESM5_ESM.xls

Additional file 5: Table S4: GO classification of the genes differentially transcribed in the contrast between libraries C and D. (XLS 45 KB)

12864_2013_5635_MOESM6_ESM.xls

Additional file 6: Table S5: Genes differentially transcribed in the contrast between libraries A and B. The criteria applied for assigning significance were: P- value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(B/A)| ≥ 1. Genes listed in descending order of absolute |log2Ratio(B/A)|. GeneIDs retrieved from the Chrysanthemum Reference Sequence Database. Annotation of unigene sequences performed using BlastX (E <10). The “GeneLength” column gives the length of exon sequence. A- and B expression: frequency of unigene transcripts in libraries A and B, respectively. A- and B-RPKM: reads per kb per million reads for each unigene in libraries A and B, respectively. Log2 Ratio(B/A): the ratio between the RPKM in B and the RPKM in A. Up-Down-Regulation(B/A), P-value and FDR of each gene are also shown. KEGG: annotation according to the KEGG database by BLAST. Blast nr: identification of homologues in GenBank. GO Component, GO Function and Go Process: ontology information of Cellular Components, Molecular Function and Biological Processes of Gene-corresponding GO terms. “-”: no hit. (XLS 106 KB)

12864_2013_5635_MOESM7_ESM.xls

Additional file 7: Table S6: Genes differentially transcribed in the contrast between libraries A and C. The criteria applied for assigning significance were: P- value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(C/A)| ≥ 1. Genes listed in descending order of absolute |log2Ratio(C/A)|. GeneIDs retrieved from the Chrysanthemum Reference Sequence Database. Annotation of unigene sequences performed using BlastX (E <10). The “GeneLength” column gives the length of exon sequence. A- and C expression: frequency of unigene transcripts in libraries A and C, respectively. A- and C–RPKM: reads per kb per million reads for each unigene in libraries A and C, respectively. Log2 Ratio(C/A): the ratio between the RPKM in C and the RPKM in A. Up-Down-Regulation(C/A), P-value and FDR of each gene are also shown. KEGG: annotation according to the KEGG database by BLAST. Blast nr: identification of homologues in GenBank. GO Component, GO Function and Go Process: ontology information of Cellular Components, Molecular Function and Biological Processes of Gene-corresponding GO terms. “-”: no hit. (XLS 286 KB)

12864_2013_5635_MOESM8_ESM.xls

Additional file 8: Table S7: Genes differentially transcribed in the contrast between libraries B and D. The criteria applied for assigning significance were: P- value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(D/B)| ≥ 1. Genes listed in descending order of absolute |log2Ratio(D/B)|. GeneIDs retrieved from the Chrysanthemum Reference Sequence Database. Annotation of unigene sequences performed using BlastX (E <10). The “GeneLength” column gives the length of exon sequence. B- and D expression: frequency of unigene transcripts in libraries B and D, respectively. B- and D-RPKM: reads per kb per million reads for each unigene in libraries B and D, respectively. Log2 Ratio(D/B): the ratio between the RPKM in D and the RPKM in B. Up-Down-Regulation(D/B), P-value and FDR of each gene are also shown. KEGG: annotation according to the KEGG database by BLAST. Blast nr: identification of homologues in GenBank. GO Component, GO Function and Go Process: ontology information of Cellular Components, Molecular Function and Biological Processes of Gene-corresponding GO terms. “-”: no hit. (XLS 688 KB)

12864_2013_5635_MOESM9_ESM.xls

Additional file 9: Table S8: Genes differentially transcribed in the contrast between libraries C and D. The criteria applied for assigning significance were: P- value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(D/C)| ≥ 1. Genes listed in descending order of absolute |log2Ratio(D/C)|. GeneIDs retrieved from the Chrysanthemum Reference Sequence Database. Annotation of unigene sequences performed using BlastX (E <10). The “GeneLength” column gives the length of exon sequence. C– and D expression: frequency of unigene transcripts in libraries C and D, respectively. C- and D-RPKM: reads per kb per million reads for each unigene in libraries C and D, respectively. Log2 Ratio(D/C): the ratio between the RPKM in D and the RPKM in C. Up-Down-Regulation(D/C), P-value and FDR of each gene are also shown. KEGG: annotation according to the KEGG database by BLAST. Blast nr: identification of homologues in GenBank. GO Component, GO Function and Go Process: ontology information of Cellular Components, Molecular Function and Biological Processes of Gene-corresponding GO terms. “-”: no hit. (XLS 222 KB)

12864_2013_5635_MOESM10_ESM.xls

Additional file 10: Table S9: The transcription level of each unigene derived from the number of relevant reads recovered in the four libraries. The “GeneLength” column gives the length of exon sequence. (XLS 8 MB)

12864_2013_5635_MOESM11_ESM.doc

Additional file 11: Table S10: The differential transcription of Wall-associated receptor kinase-like (WAK-like) genes in the contrast B vs D. The criteria applied for assigning significance were: P-value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(D/B)| ≥ 1. RPKM: reads per kb per million reads. (DOC 116 KB)

12864_2013_5635_MOESM12_ESM.doc

Additional file 12: Table S11: The differential transcription of Wall-associated receptor kinase-like (WAK-like), brassinosteroid insensitive 1 (BRI-like), somatic embryogenesis receptor kinase (SERK), and BRI1-associated receptor kinase 1 (BAK1) genes in the contrast A vs C. The criteria applied for assigning significance were: P-value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(C/A)| ≥ 1. RPKM: reads per kb per million reads. (DOC 120 KB)

12864_2013_5635_MOESM13_ESM.doc

Additional file 13: Table S12: The differential transcription of leucine-rich repeat receptor-like kinase (LRR-RLK), brassinosteroid insensitive 1 (BRI-like), BRI1-associated receptor kinase 1 (BAK1), and somatic embryogenesis receptor kinase (SERK) genes in the contrast B vs D. The criteria applied for assigning significance were: P-value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(D/B)| ≥ 1. RPKM: reads per kb per million reads. (DOC 130 KB)

12864_2013_5635_MOESM14_ESM.doc

Additional file 14: Table S13: The differential transcription of cysteine-rich receptor-like protein kinase (CRKs) genes in the contrast B vs D. The criteria applied for assigning significance were: P-value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(D/B)| ≥ 1. RPKM: reads per kb per million reads. (DOC 114 KB)

12864_2013_5635_MOESM15_ESM.doc

Additional file 15: Table S14: The differential transcription of cysteine-rich receptor-like protein kinase (CRKs) genes in the contrast A vs C. The criteria applied for assigning significance were: P-value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(C/A)| ≥ 1. RPKM: reads per kb per million reads. (DOC 115 KB)

12864_2013_5635_MOESM16_ESM.doc

Additional file 16: Table S15: The differential transcription of respiratory burst oxidase, and alpha-dioxygenase genes in the contrast B vs D. The criteria applied for assigning significance were: P-value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(D/B)| ≥ 1. RPKM: reads per kb per million reads. (DOC 116 KB)

12864_2013_5635_MOESM17_ESM.doc

Additional file 17: Table S16: The differential transcription of respiratory burst oxidase, and alpha-dioxygenase genes in the contrast A vs C. The criteria applied for assigning significance were: P-value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(C/A)| ≥ 1. RPKM: reads per kb per million reads. (DOC 114 KB)

12864_2013_5635_MOESM18_ESM.doc

Additional file 18: Table S17: The differential transcription of photosynthesis and circadian rhythm-related genes in the contrast B vs D. The criteria applied for assigning significance were: P-value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(D/B)| ≥ 1. RPKM: reads per kb per million reads. (DOC 158 KB)

12864_2013_5635_MOESM19_ESM.doc

Additional file 19: Table S18: The differential transcription of photosynthesis and circadian rhythm-related genes in the contrast C vs D. The criteria applied for assigning significance were: P-value < 0.05, FDR ≤ 0.001, and estimated absolute |log2Ratio(D/C)| ≥ 1. RPKM: reads per kb per million reads. (DOC 118 KB)

Additional file 20: Table S19: High-throughput sequencing metadata. (XLS 61 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver ( https://creativecommons.org/publicdomain/zero/1.0/ ) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Cite this article

Li, H., Chen, S., Song, A. et al. RNA-Seq derived identification of differential transcription in the chrysanthemum leaf following inoculation with Alternaria tenuissima. BMC Genomics 15, 9 (2014). https://doi.org/10.1186/1471-2164-15-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-15-9

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us