Next Article in Journal
Pharmacological Targeting of the ER-Resident Chaperones GRP94 or Cyclophilin B Induces Secretion of IL-22 Binding Protein Isoform-1 (IL-22BPi1)
Next Article in Special Issue
Salicylic Acid Signals Plant Defence against Cadmium Toxicity
Previous Article in Journal
Voluntary Activity Modulates Sugar-Induced Elastic Fiber Remodeling in the Alveolar Region of the Mouse Lung
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Involvement of Salicylic Acid in Anthracnose Infection in Tea Plants Revealed by Transcriptome Profiling

1
Tea Research Institute, Zhejiang University, # 866 Yuhangtang Road, Hangzhou 310058, China
2
Ningbo Huangjinyun Tea Science and Technology Co. Ltd., Yuyao 315412, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2019, 20(10), 2439; https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20102439
Submission received: 9 April 2019 / Revised: 15 May 2019 / Accepted: 16 May 2019 / Published: 17 May 2019
(This article belongs to the Special Issue Salicylic Acid Signalling in Plants)

Abstract

:
Anthracnose is a major leaf disease in tea plant induced by Colletotrichum, which has led to substantial losses in yield and quality of tea. The molecular mechanism with regards to responses or resistance to anthracnose in tea remains unclear. A de novo transcriptome assembly dataset was generated from healthy and anthracnose-infected leaves on tea cultivars “Longjing-43” (LJ43) and “Zhenong-139” (ZN139), with 381.52 million pair-end reads, encompassing 47.78 billion bases. The unigenes were annotated versus Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), National Center for Biotechnology Information (NCBI) non-redundant protein sequences (Nr), evolutionary genealogy of genes: Non-supervised Orthologous Groups (eggNOG) and Swiss-prot. The number of differential expression genes (DEGs) detected between healthy and infected leaves was 1621 in LJ43 and 3089 in ZN139. The GO and KEGG enrichment analysis revealed that the DEGs were highly enriched in catalytic activity, oxidation-reduction, cell-wall reinforcement, plant hormone signal transduction and plant-pathogen interaction. Further studies by quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) and high-performance liquid chromatography (HPLC) showed that expression of genes involved in endogenous salicylic acid biosynthesis and also accumulation of foliar salicylic acid are involved in the response of tea plant to anthracnose infection. This study firstly provided novel insight in salicylic acid acting as a key compound in the responses of tea plant to anthracnose disease. The transcriptome dataset in this study will facilitate to profile gene expression and metabolic networks associated with tea plant immunity against anthracnose.

Graphical Abstract

1. Introduction

Anthracnose is a key leaf disease in the tea plant (Camellia sinensis (L.) O. Kuntze) induced by Colletotrichums [1,2], which has caused a great loss in tea yield and quality owing to its strong infectivity and widespread distribution [3]. Although there have been hundreds of Colletotrichum isolates detected in tea plant [4,5], the dominant species or major pathogen on the tea plant are considered to be Colletotrichum gloeosporioidesis, C. camelliae and C. fructicola [3,4,5], which grow and spread quickly in humid and hot summer and autumn seasons, leading to gray sunken or shrunk necrotic lesions on tea leaves and twigs [6,7]. The loss of tea yield induced by anthracnose was estimated to range from 5% to 20% or even more, depending on the cultivar planted [8,9]. Although plant activators [10], antagonistic bacteria [11] and agronomic techniques such as leaf trimming and rational fertilization have been used to control the anthracnose disease, spraying fungicides is still a major control method in tea fields, resulting in a serious fungicides residue problem. Breeding tea cultivars resistant to anthracnose is considered to be the most effective measure to control anthracnose because there is a great difference in resistance to this disease between various tea cultivars [6,8,9].
Revealing molecular defense responses of tea plant to anthracnose infection and developing a molecular assisted selection (MAS) method are badly needed in the tea plant breeding field. However, little is known about the molecular mechanisms regulating the defense response in tea plants [6] although attempts were made to probe this. Transcriptional analysis and histochemistry revealed that the hypersensitive response (HR) and H2O2 play critical roles in tea plant defense response to C. fructicola [12], and chemical changes of caffeine is considered to be associated with tea–fungi interaction [13]. Nonpathogenic species of Colletotrichum was more vulnerable to catechins and caffeine, and differentiation in secondary metabolites might be an important factor leading to the difference in pathogenicity between cultivars [14]. Different communities of Colletotrichum with little variability within internal transcribed spacer (ITS) and glyceraldehyde 3-phosphate dehydrogenase (GAPDH) regions in the DNA sequence have different virulence [15].
Generally, it was reported that salicylic acid (SA), a secondary metabolite belonging to phenolic compounds, acts as a plant signal or hormone with a regulatory role in a variety of physiological processes under abiotic and biotic stresses, especially as a module in activating disease resistance [16,17,18]. The endogenous SA concentration is related to the conditions of stresses. For instance, SA accumulation, which is responsible for the phenylpropanoid synthesis pathway, was observed in wheat leaves after heavy metal cadmium treatment [19]. In tobacco (Nicotiana tabacum L. cv. Xanthinc), a transient sharp increase in SA induced by ultraviolet UV light, ozone and mosaic virus parallelly led to pathogenesis-related proteins accumulation and enhanced disease resistance [20]. The Arabidopsis thaliana mutants of sid and pad4, which are involved in synthesis of SA and camalexin (a kind of phytoalexin), were defective in SA synthesis, resulting in them being more susceptible to pathogens [21,22]. The defense compounds including pathogenesis-related proteins (PR-1, PR-2, PR-5), camalexin and H2O2 could be induced by endogenous SA [22,23,24]. Transgenic plants with salicylate hydroxylase, which converts SA to other compounds, accumulated almost no SA and lost the resistance to virus [25]. When blocked by an inhibitor of phenylpropanoid (precursor of SA) synthesis, the exogenous SA recovered the ability to plant resistance [23]. All these show that SA indeed plays the key role and an important signaling effect in plant resistance, especially in systemic acquired resistance [25,26]. Due to the fact that there has been no study revealing the relationship of SA to the response of tea plant to diseases, it will be interesting to reveal the responses of SA biosynthesis genes and molecular function of SA in the defense response of tea plant to anthracnose infection.
RNA-seq can be focused on assessing the degree of RNA processing and the types of RNA modification, both of which may play an important role in the disease infection process [27,28]. Information of RNA-seq can be used to profile gene expression levels and to reveal the genes involved in plant defense responses [29,30]. RNA-seq analysis needs no reference genome [31] even though tea plant genome data has been available [32]. Developing a transcriptome assembly dataset generating from anthracnose infected tea leaves will be helpful to reveal the molecular defense responses of tea plant to the anthracnose and also to mine molecular markers for MAS used in resistance breeding.
In the present study, transcriptome sequencing by Trinity was performed on the platform of Illumina HiSeq to construct a de novo transcriptome assembly database generated from anthracnose infected tea leaves of two susceptible tea cultivars “Longjing 43” (LJ43) and “Zhenong 139” (ZN139), in which the unigenes were generated and annotated, and important functional genes and metabolic pathways were also revealed. The data sets will provide references for further study on gene expression profiles, biochemical processes and regulation networks associated with tea plant immunity against anthracnose.

2. Results

2.1. Symptoms of Infected Tea Plants

We inspected the typical anthracnose symptoms in Colletotrichum susceptible tea cultivars “Longjing-43” (LJ43) and “Zhenong-139” (ZN139). The infected leaves were observed to be gray sunken or shrunk necrotic lesions (Figure 1a). In the late stage, the pathogens made the leaves partially withered, fragile, or easily broken, compared to the healthy leaf (Figure 1a left). The infected leaves of both cultivars “LJ43” and “ZN139” showed the same symptoms and the major pathogen isolated from the infected leaves was Colletotrichum gloeosporioides. A detailed view of anthracnose-infected tea plant occurrence in tea fields as well as sample preparation and experiment scheme are shown in Figure 1.

2.2. Transcriptome Profiling

2.2.1. Sequencing Quality, Assembly Characterization and Functional Annotation

Transcriptome sequencing was performed using a de novo assembly (Figure 2) to reveal molecular resistance or defense information in the host tea plant. Based on Illumina sequencing, we generated 381.52 million pair-end reads encompassed 47.78 billion bases, in total, from the two cultivars with two biological replicants. The percentage of N (ambiguous bases), Q20 (reads with mean error rate <1%), Q30 (reads with mean error rate < 0.1%), GC (guanine-cytosine content) were 0.00%, 95%, 89%, 48.5%, respectively. After quality filtering, about 98% clean reads from the raw reads in each sample were reserved for a de novo assembly (Table 1).
After assembling, 335,186 contigs and 109,316 unigenes were generated from cultivar LJ43, while 352,038 contigs and 115,953 unigenes were generated from cultivar ZN139. Details of total length, mean length, N50 (length-weighted median, i.e., 50% contig length), N50 sequence No. (the number of sequences longer than N50), N90 (length-weighted 90% contig length) and N90 sequence No. (the number of sequences longer than N90) were listed in Table 2.
All unigene sequences were aligned against the five databases, i.e., GO (Gene Ontology), KEGG (Kyoto Encyclopedia of Genes and Genomes), Nr (NCBI non-redundant protein sequences), eggNOG (evolutionary genealogy of genes: Non-supervised Orthologous Groups) and Swiss-prot (Table 2). In cultivar LJ43, 45,230 unigenes (41.4%) were annotated by Nr, 24,892 unigenes (22.8%) by GO, 5799 unigenes (5.3%) by KEGG, 43,025 unigenes (39.4%) by eggNOG and 33,895 unigenes (31.0%) by Swiss-prot, while in cultivar ZN139, 46,383 unigenes (40.0%) were annotated by Nr, 25,168 unigenes (21.7%) by GO, 5873 unigenes (5.1%) by KEGG, 44,067 unigenes (38.0%) by eggNOG and 34,860 unigenes (30.1%) by Swiss-prot. A total of 4605 unigenes (4.21%) and 4677 unigenes (4.03%) in two cultivars were annotated in all five reference databases.
Homology searching against the Nr database gave the explanations containing species distribution of the top hits, E-value distribution and similarity distribution (Figure S1). Importantly, about 20.3% unigenes shared the highest homologies to Coffea canephora genes. In the eggNOG analysis, annotated unigenes were classified into 26 items with similar descriptions in both cultivars LJ43 (Figure S2a) and ZN139 (Figure S2b). Based on GO classification, all unigenes were classified to three main categories, i.e., molecular function, cellular component and biological process (Table S1). Catalytic activity (GO: 0003824), binding (GO: 0005488), transporter activity (GO: 0005215) and structural molecule activity (GO: 0005198) were the vital subcategories under the category of molecular function. As for cellular component, cell (GO: 0005623), cell part (GO: 0044464), membrane (GO: 0016020), organelle (GO: 0043226) and membrane part (GO: 0044425) were the primary subcategories. With the respect to biological process category, the top five subcategories were metabolic process (GO: 0008152), cellular process (GO: 0009987), single-organism process (GO: 0044699), biological regulation (GO: 0065007) and localization (GO: 0051179). Moreover, all the annotation results of unigenes for both LJ43 and ZN139 were deposited at the public database Figshare (http://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.7706237.v1, 12 February 2019).

2.2.2. Differential Expression and Enrichment Analysis

Differential expression genes (DEGs) in both cultivars were calculated through anthracnose-infected leaves versus healthy leaves. As shown in Table 3, 1621 DEGs were found in infected leaves of cultivar LJ43, with 1082 up-regulated and 539 down-regulated, while 3089 DEGs were found in infected leaves of cultivar ZN139, with 1527 up-regulated and 1562 down-regulated, compared to the healthy leaves; 755 and 1487 unigenes of the transcripts with GO IDs were mapped by Blast2GO, respectively. Among the DEGs with p < 0.05, 7 were annotated by GO, with 4 for biological process and 3 for cellular component in cultivar LJ43 (Figure 3a), while 26 were annotated by GO, with 8 for the biological process, 7 for the cellular component and 11 for the molecular function in cultivar ZN139 (Figure 3b). It was found that catalytic activity, membrane, oxidation-reduction process, cell periphery and carbohydrate metabolic process were strongly responsible to the infection of Colletotrichum. Other cellular component categories including the extracellular region, cell wall and external encapsulating structure, as well as the secondary metabolic process were next enriched.
Of course, there are also questions to be explained in the present study. As shown in Figure 3, there was great difference in the GO enrichment between cultivars LJ43 and ZN139, but the exact details leading to the difference remains unknown. It was reported that caffeine and catechins such as (−)-epigallocatechin gallate and (+)-catechin may be involved in the resistance of tea plants to anthracnose [13]. LJ43 has higher level of polyphenols than ZN139 although there are no differences in amino acids and caffeine [33]. Difference in tea polyphenols might be one of the many factors leading to difference in GO enrichment between the two cultivars. Differences in gene abundance might also be partially responsible for the difference in GO enrichment.
As we focused on the top 25 KEGG pathways (Figure 4), it was found that the DEGs were enriched in the pathways involving in plant hormone signal transduction (ko04075), plant-pathogen interaction (ko04626), starch and sucrose metabolism (ko00500) and mineral absorption (ko04978) in both cultivars. Besides, flavonoid biosynthesis (ko00941), phenylpropanoid biosynthesis (ko00940), pentose and glucuronate interconversions (ko00040), cutin, suberine and wax (ko00073), biosynthesis metabolism of xenobiotics by cytochrome P450 (ko00980), and drug metabolism-cytochrome P450 (ko00982) also responded to the anthracnose infection (Figure 4).
The differential expression and enrichment analysis results were deposited at the public database Figshare (http://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.7706237.v1, 12 February 2019).

2.2.3. Visualization of Two Vital Plant Metabolic Pathways

Pant hormone signal transduction (ko04075, Figure 5) and plant-pathogen interaction (ko04626, Figure S3) were vital plant metabolic pathways involved in anthracnose infected process. The nomenclature of gene was based on the entry name from KEGG in this paper. The DEGs mapped to these two pathways showed a same trend in both tested tea cultivars. When the two KEGG pathway results were integrated in one figure (Figure 5), it shows that the genes involving in salicylic acid signal transduction pathway, including NPR1 (nonexpressor of pathogenesis-related gene 1), TGA (TGACG motif-binding factor) and pathogenesis-related protein 1 (PR1) were upregulated significantly. PR1 is directly linked to disease resistance. Some plant hormones with signal transduction components, such as auxin (AUX1, AUX/IAA, GH3, SAUR), cytokinine (AHP, A-ARR), gibberellin (GID1, DELLA), jasmonic acid (JAR1) also play an important role in cell enlargement and division, plant growth, induced germination and stress response. As mentioned in pathogen-associated molecular patterns (PAMP) -triggered immunity, RBOH (respiratory burst oxidase), CALM/CLM (calmodulin), WRKY 25/33 and PR1 were found to be upregulated in the anthracnose infected tea leaves. Upstream CNGC (cyclic nucleotide gated channel) involved in calcium ion transferring was downregulated. The changes of gene expression level might be related to HR, cell wall reinforcement, stomatal closure and defense-related gene induction. The expression level of R genes (RPM1, PBS1, EDS1) was upregulated, while the downstream HSP90 (heat shock protein 90kDa beta) was downregulated in effector-triggered immunity.
It is reported that system acquired resistance (SAR) is dependent on salicylic acid signaling and systemic expression of PR genes. Immune system resistance (ISR) depends on ethylene and jasmonic acid (JA) but is not associated with the expression of the PR genes. Both SAR and ISR do result in broad spectrum resistance. Although an antagonistic interaction between SA and JA pathways was revealed, synergistic interactions were also observed [34]. It is considered that there may be a link between JAR1 and SA-related genes. We searched using key words of brassinosteroids, spermine and heat shock transcription factor in the datasets and found that only heat shock transcription factor (HSP90, being down-regulated by infection) contributed to protect from pathogen infection (Figure 5).

2.3. Quantitative Reverse Transcriptase Polymerase Chain Reaction (qRT-PCR) Validation of Salicylic Acid Signaling Related Genes

To further verify the sequencing results, expression levels of genes involving in salicylic acid signal transduction pathway in healthy leaves and anthracnose infected leaves of cultivars LJ43 and ZN139 were tested by qRT-PCR (quantitative reverse transcriptase polymerase chain reaction). Genes ALD1 (AGD2-like defense response gene 1, KEGG Orthology K10206), TGA2 and TGA3 (K14431), PR1 (K13449) in the anthracnose infected leaves showed significant higher expression level than healthy leaves (Figure 6). The expression of NPR1 (K14508) differentiated between cultivar ZN139 and cultivar LJ43, with upregulation in anthracnose infected leaves of LJ43 but no significant difference between healthy leaves (ZN139_H) and infected leaves (ZN139_D) of cultivar ZN139. NPR1 is constitutively expressed in many plants and is activated by modification after infection rather than at the transcriptional level. However, our transcriptome dataset shows that the gene encoding predicted regulatory protein NPR3 isoform X1 (LJ69756_g5) was significantly upregulated in the infected leaves, but it was not found in SA pathway. The reason why NPR1 and NPR3 were upregulated during the anthracnose infection remains to be investigated.
Correlation analysis confirmed that relative expression values generated by the two methods were consistent, with a linear correlation coefficient R-square of 0.878 (Figure 6c), suggesting that endogenous salicylic acid is involved in the response to anthracnose in the tea plant.

2.4. Salicylic Acid Content

Concentrations of free salicylic acid (SA) and bound SA (the inactive storage form 2-O-β-d-glucosylsalicylic acid) [35] were determined by high-performance liquid chromatography (HPLC) (Figure 7) and it showed that levels of free SA, bound SA and total SA in the anthracnose infected leaves were significantly higher than healthy leaves, suggesting that salicylic acid is an important signaling agent responding to anthracnose infection in tea plant.

3. Discussion

Profiling plant immune systems and the pathogen molecules as well as their interaction will provide extraordinary insights into molecular recognition and cell biology of plants [36]. Next-generation molecular approaches such as transcriptomics, metabolomics, whole-genome sequencing and proteomics can be used for this purpose [37]. In the present study, a systematic and accurate transcriptome of healthy leaves and anthracnose-infected leaves of two susceptible tea cultivars was profiled, whose dataset will be of significance for further revealing molecular mechanism of tea plant resistance to anthracnose and mining molecular markers for MAS in tea plant breeding.
Host gene expression fluctuates after pathogen infection and the timing of sampling is important for transcriptome study. However, there have been no studies focused on the expression fluctuation of host genes after pathogen infection. Field investigation showed that the latent period after inoculation of anthracnose pathogens was about 2 weeks [38]. Laboratory tests showed that obvious symptom cold be observed and the resistance capacities of various tea cultivars could be identified after two weeks of wound-inoculation of anthracnose pathogens [8]. In Hangzhou, where we sampled the leaves, the rainy season occurs in June, with humid weather during which anthracnose infection usually takes place. The symptom is observed in early July and so we sampled in early July. Also, tea is a self-sterile plant and all tea cultivars we are planting are cross-pollinated hybrids with a genetically complex background. It would be ideal to obtain a resistant mutant from a disease susceptible tea cultivar. Unfortunately, we have not got one. If a gene showed a same response trend to the infection, it is more definite to conclude that the response has the generality among the susceptible cultivars and is more likely to be related to anthracnose infection when two susceptible cultivars were chosen to be tested.
The present transcriptomic profiling revealed that the oxidation-reduction, catalytic activity, cell wall, membrane and carbohydrate metabolic process in GO analysis (Figure 3), as well as the pathways of Ca2+ signaling and SA signaling in KEGG analysis, were strongly responsible to anthracnose infection. SA signaling pathways with remarkable downstream defense gene expression are validated by qRT-PCR in this paper.
Catalytic activity and oxidation-reduction commonly exist in many biological processes or reactions. Most plant cells generate reactive oxygen species (ROS) [39], including hydrogen peroxide, singlet oxygen and hydroxyl radical [40] during their interaction with pathogens, for example Phytophthora infestans [41] and Colletotrichum lindemuthianum [42]. The generation and accumulation of hydrogen peroxide with a membrane-bound NADPH (nicotinamide adenine dinucleotide phosphate) oxidase is a symptom of the oxidative burst [43], which coincides with programmed cell death (PCD) during HR. The ROS intermediates redox signaling and oxidation-reduction in host plant resistance [44]. Previous study also showed that hydrogen peroxide and HR played a crucial role in tea plant defense [12].
The modification and reinforcement of the cell wall (RCW) as well as deposition of papillae containing callose are considered to be a structural barrier, which can block the penetration of fungal pathogens [42,45]. The detection of hydrogen peroxide at subcellular levels showed that the cross-linking between papillae and HR cells help the pathogen arrest by reinforcing the cell wall apposition [46]. A genetic network of cell-wall damage, regulated by lignin biosynthesis, was characterized in Arabidopsis thaliana (Col) [47]. Plants may have evolved a system of dynamic cell wall remodeling to prevent diseases [48]. Many researches elucidated that the changing of cell wall is a way to hinder pathogens entering the plant cells. Hence, we speculated that the cell wall and callose, with complex chemical structure and physiological function, have direct antifungal effects on facilitating tea plant defense.
KEGG analysis revealed that many DEGs were associated with several pathways, particularly in plant hormone signal transduction (Figure 5) and plant-pathogen interaction (Figure S3) in both cultivars. As shown in Figure 4 and Figure S3, a defense-related gene, HR and PCD could be induced through pathogen-associated molecular patterns triggered immunity (PTI) and effector-triggered immunity (ETI) [36]. Ca2+ is a fateful factor in PTI, resulting in three regulating pathways. One is the ROS-related pathway, in which upstream RBOH with NADPH oxidases [49] is enhanced highly in anthracnose-infected leaves. The generation of ROS acts as a signal to induce HR and RCW. Ca2+ has been reported to activate RBOHs, which is directly phosphorylated responsible to PAMPs perception, at the region of cytosolic N-terminal [50,51]. Ectopic and heterologous expression in the potato (Solanum tuberosum cv Rishiri) suggests that calcium-dependent protein kinases 5 regulates the ROS production via St RBOHB (RBOH) [52]. Another is the calmodulin-related pathway. The expression of calmodulin may result in RCW and stomatal closure. A study on barley (mlo mutations) against mildew showed that plant defense was regulated by Ca2+ dependent interaction with calmodulin and mildew resistance locus O protein, a modulator of plant defense and cell death, in early signaling cascades [53]. Calmodulin detection by enzyme-linked immunosorbent assay (ELISA) showed that calmodulin promoted cell wall regeneration and cell division [54]. The third pathway is related to WRKY transcription factor 25/33, whose high expression induces accumulation of downstream defense proteins encoded by gene PR1. The present results demonstrate that Ca2+ is a key component during plant-pathogen interactions, like that described in a previous study [55].
As for ETI in Figure S3, specific disease resistance genes (R) recognize the effectors that enable pathogens overcoming PTI and ETI [36]. When tea plant cell comes into contact with effectors delivered by fungus, the expression level of nucleotide binding leucine rich repeat domains (NB-LRR) proteins [56] encoded by R genes (RPM1, PBS1, EDS1) increased markedly, resulting in activation of ETI with a threshold of HR, PCD and defense amplification. In contrast, heat shock protein 90 (HSP90) was down-regulated, which was perhaps related to RPM1 function [57]. A study on heat-shock treatment reported that inhibition of HSP90 caused enhancement of disease resistance by released heat shock transcription factor 1 in tomato cultivar Natsunokoma [58], which is consistent with our findings.
It is known that SA, JA and ethylene are three key phytohormones responding to biotic stresses in plant. It is suggested that effective defense against biotrophic or hemibiotrophic pathogens is largely due to PCD in the host, and to associated activation of defense responses regulated by the SA-dependent pathway, while necrotrophic pathogens, benefiting from host cell death, are not limited by SA-dependent defenses, but by JA and ethylene-signaling defense responses [59,60]. Figure 5 shows that the remarkable up-regulation of NPR1, TGA and PR1 genes induced by SA signaling finally facilitates the PR proteins associating with disease resistance to Colletotrichum. Based on SA signaling, the effect of Colletotrichum on tea leaves is likely to be hemibiotrophic. The qRT-PCR testing confirmed that SA signaling was involved in tea resistance to anthracnose, in which notably high expression levels of NPR1, TGA2, TGA3 and PR1 in anthracnose infected leaves were detected in both cultivars (Figure 6). The accumulation of the PR1 protein with antifungal activity enhances plant immunity in the vacuole or extracellular space such as cell wall [26,61]. Among these genes, we demonstrated here that PR1 can be used as a marker for SA-dependent gene expression, being consistent with that in Arabidopsis [62]. NPR1 is a critical regulator and receptor, directly binding to SA [63,64] and a test on yeast indicated that the TGA family was a link between NPR1 and PR1 expression [65]. However, another study showed that there were more genes in SA signaling, not limited to those we confirmed before, such as a high affinity SA-binding protein (SABP2) [66]. WRKY factors act as repressors or activators in system acquired resistance (SAR) network, which are induced by SA [67,68]. The present HPLC analysis revealed that the levels of free SA (active form), bound SA (inactive form) and total SA in tea leaves were significantly increased by Colletotrichum infection (Figure 7). The results are in accordance with those in Arabidopsis thaliana (npr1 and NahG) showing that antimicrobial proteins encoded by the PR1 superfamily were increased with the increase in endogenous SA levels [69], and also are in accordance with those described in two previous studies on tobacco (Nicotiana tabacum L. cv. Xanthinc) [20]. The increase in various forms of SA suggests that the free SA is directly involved in tea plant resistance to anthracnose infection and the bound SA might be a stored form of SA, with a potential for defense resistance. In addition, the expression of ALD1, which is involved in the lysine metabolic pathway, was extremely enhanced in the anthracnose infected leaves. The previous study on Arabidopsis indicated that ALD1 involving in mediation of plant immunity was more similar to SA-triggered immunity [62]. ALD1-triggered pipecolic acid (PA) and N-hydroxypipecolic acid, whose molecular structures are similar to SA, are largely accumulated in the disease infected plant, where they act as a regulator of SAR activated by SA [70]. The cross-links between SA and PA in tea plant defense need to be further investigated. A hypothetical model was clarified in Figure 8 for SA accumulation and related gene expression involving in the SA signaling pathway in tea immunity.
In conclusion, PR1 and endogenous SA, acting as a key compound, clearly play a pivotal role in defense activation of tea immunity to anthracnose. Fully understanding the pathogen molecule pathways and the interaction between Colletotrichum and the tea plant will be helpful to develop agricultural measures to control tea plant anthracnose disease and also to mine molecular markers used in MAS. Despite many recent insights into the tea–Colletotrichum interactions by several assays, there is still a large gap in knowledge of mechanisms for further study, such as focuses on mapping to the reference genome, transgenic tea verification and field testing on exogenous SA, SA agonist or SA inhibitor. Furthermore, many other related topics mentioned above, such as cross talks among signaling pathways, other plant hormones, metabolic processes, phytoalexins and etc., need to be verified in future work.

4. Materials and Methods

4.1. Plant Material

Tea cultivars “Longjing 43” (LJ43) and “Zhenong 139” (ZN139) which are susceptible to anthracnose were used as plant materials in this study. The tea plants were grown in the Experimental Tea Farm of Zhejiang University (Hangzhou, China, 30.3 N, 120.2 E). Leaf samples (10 leaves each sample) were collected in early July 2017 and the leaf slices were cut as described in Figure 1. The responses of a plant to pathogen infections were induced systemically. The healthy leaves were sampled from the individuals which had not completely infected with any part so as to avoid the effects of systematic responses. The leaves were washed with distilled water before cutting the leaf slices to eliminate contamination of RNA from the pathogens and the leaf slices were frozen in liquid nitrogen immediately and then stored at −80 °C until RNA extraction. The healthy leaf samples from LJ43 and ZN139 were labeled as LJ43_H and ZN139_H, and the anthracnose infected leaf samples from LJ43 and ZN139 were labeled as LJ43_D and ZN139_D, respectively (Figure 1). There were two biological replicates for each treatment in the study of de novo transcriptome assembly, and three biological replicates for the studies of qRT-PCR validation and HPLC analysis of salicylic acid.

4.2. Total RNA Extraction, cDNA Library Construction and Deep Sequencing

Total RNA was extracted using RNAprep Pure Kit for Polysaccharides and Polyphenolics rich plant (Cat. no. DP441, Tiangen Biotech Co., Ltd., Beijing, China) according to the instructions of the kit. Concentration and quality of the RNA were assessed using Agilent 2100 Bioanalyzer (Agilent Technologies Inc., Santa Clara, USA), and those with 1.8 ≤ OD260 nm/OD280 nm ≤ 2.1 and RIN ≥ 6.5 were used in subsequent tests. RNA integrity was evaluated by agarose gel electrophoresis. Three micrograms of RNA were used as input material for cDNA preparations.
The cDNA library was constructed according to the protocol of TruSeqRNA Sample Prep Kit (Illumina Inco., San Diego, CA, USA). The quality and quantity of libraries were checked by fluorospectrophotometry (Quant-iTPicoGreen dsDNA Assay Kit, Invitrogen, P7589; Quantifluor-ST fluorometer E6090, Promega, CA, USA) and Agilent 2100 Bioanalyzer with Agilent Bioanalyzer High Sensitivity DNA chip Kit (5067-4626, Agilent Technologies Inc., Santa Clara, USA). Each normalized cDNA library (2 nM) was gradually diluted to 4–5 pM and the sequencing was performed on Illumina HiSeq X Ten Platform in 2 × 150 bp pair-end sequencing mode (Figure 2a) and the original data in FASTQ format (Raw Data) was then generated. The quality assessment of raw data in FASTQ format was done with FastQC [71] (http://www.bioinformatics.babraham.ac.uk/projects/fastqc, 18 January 2018, Q-score ≥ Q20, min read length = 50 bp). The connectors were removed by Cutadapt (Version 1.15) and clean reads were prepared for the further analysis (Figure 2a).

4.3. De Novo Assembly and Functional Annotation

The present tests were carried out in 2017, when the reference genome of tea plant had not been published. Accordingly, de novo assembly was performed to reconstruct the transcriptomes. We used Trinity (Version r20140717, default k-mer = 25 bp) [72,73], consisting of three software modules: Inchworm, Chrysalis and Butterfly, to efficiently process massive clean reads based on robust reconstruction and analysis of De Bruijn Graph (DBG) [72]. Briefly, clean reads with overlap joints were firstly combined to form longer fragments, i.e., contigs. Related pair-end reads and contigs were clustered to form unigenes. Finally, non-redundant unigenes for functional annotation were obtained [73].
The unigene sequences were compared using the local BLASTX [74] program (E-value ≤1e-5) against 5 public databases: Nr [75], GO [76], KEGG [77], Swiss-Prot [78] and eggNOG [79]. GO terms were identified on the platform BLAST2GO [80]. KEGG Pathways and KEGG Orthologs were annotated by KAAS [81] (KEGG Automatic Annotation Server, http://www.genome.jp/tools/kaas/, 18 January 2018).

4.4. Differential Expression Genes (DEGs) Analysis

Firstly, the software RSEM [82] was used to compare read count values on each unigene as original expression, and then standardized as FPKM (expected number of fragments per kilobase of transcript sequence per millions of base pairs sequenced). Then DESeq (Version 1.30.0) was used to analyze DEGs [83]. The screen criteria were expression difference multiple |log2 Fold Change| > 1 (p < 0.05). The bi-directional clustering analysis of DEGs in different samples was performed using R language. Then we mapped all the DEGs to GO databases to reveal the gene enrichment in each term. Furthermore, the DEGs involved in the KEGG pathways including signaling and metabolic pathways were assigned. All the downstream analysis was described in Figure 2b.

4.5. qRT-PCR Validation

To verify the expression levels of DEGs obtained from RNA-seq, qRT-PCR were carried out following the previous protocol [84]. The primary DEGs involved in salicylic acid signal transduction pathway were chosen for validation and specific primers were designed using the PrimerQuest Tool (https://sg.idtdna.com/PrimerQuest/Home/Index, 10 April 2018–23 April 2018). qRT-PCR was run in the machine StepOne Plus (Applied Biosystems, Foster City, CA, USA). The qRT-PCR test was generated using triplicate samples. Each relative expression level was calculated using the 2−ΔΔCT method [85] with β-actin gene as internal control.

4.6. Salicylic Acid Extraction and High-Performance Liquid Chromatography (HPLC) Quantification

The extraction of total SA including free SA and bound SA was performed following the published method [35,86]. The cut leaf slices stored at −80 °C were pulverized in a mortar with liquid nitrogen condition using a pestle. About 0.1 g homogenized powders was mixed with 1 mL pre-cooling 90% methanol, stood overnight, and centrifuged at 8000× g for 10 min. The supernatants were collected while the sediment was re-suspended in 0.5 mL 90% methanol and further extracted for 2 h. The centrifugation and collection were repeated as before. The combined supernatants were evaporated using a SpeedVac machine at 40 °C to 0.3 mL aqueous solution. 20 μL trichloroacetic acid (1 mg/mL) was added to the aqueous solution and mixed by vortex. Partition with 1 mL ethyl acetate: cyclohexane (1:1, v/v) was carried out twice. The organic phase with free SA was dried by nitrogen blowing and dissolved in 0.5 mL mobile phase for HPLC test. The aqueous phase containing bound SA was acid-hydrolyze by adding 0.5 mL HCl (8 M) to release free SA. The same steps for obtaining free SA were implemented again. At this moment, the bound SA was determined in the form of free SA. The contents of free SA and bound SA in the samples were then determined by an Agilent 1100 HPLC system with fluorescence detection. Chromatographic separation was performed using a Kromasil C18 reverse phase column (5 μm, 250 mm × 4.6 mm). The mobile phase was a mixture of 1% acetic acid aqueous solution with 1% phosphoric acid: methanol (40:60, v/v). The flow rate was set at 0.8 mL/min for 50 min in total. The column temperature was set at 35 °C and injection volume was 10 μL. Excitation and emission wavelengths were 294 nm and 426 nm, respectively.
The data of qRT-PCR and HPLC were expressed as mean (mean values of three independent experiments) ± SD (standard deviation). An analysis of variance (ANOVA) test, mainly using least significant difference (LSD), was used for measuring statistical significance by means of SPSS and EXCEL software.

Supplementary Materials

Supplementary materials can be found at https://0-www-mdpi-com.brum.beds.ac.uk/1422-0067/20/10/2439/s1. Raw data files generated by RNA-seq were deposited in the NCBI Sequence Read Archive SRP162639. The files of unigenes were deposited at GenBank Transcriptome Shotgun Assembly TSA GHEL00000000. The annotation and differential expression analysis results of the unigenes were deposited at Figshare (http://0-doi-org.brum.beds.ac.uk/10.6084/m9.figshare.7706237.v1, 12 February 2019). Table S1. All subcategories of GO annotations. Table S2. qRT-PCR primers of selected unigenes. Figure S1. Characteristics of homology searching in Camellia sinensis against the Nr database. Figure S2. eggNOG functional classification of consensus sequence. Figure S3. Visualization of plant-pathogen interaction (ko04626) KEGG pathway for two tea cultivars.

Author Contributions

Y.-L.S., Y.R.L. and X.-Q.Z conceived the project; Y.-L.S., Y.-Y.S., Z.-Y.C. and X.-Y.G. prepared samples and generated the experiments; R.Y., Q.-S.L., X.-M.L., J.-L.L., J.-H.Y. and D.L. obtained and analyzed the RNA-seq data; Y.-L.S. performed bioinformatics analysis and prepared the manuscript; Z.-Y.C. and R.Y. optimized figures and tables; K.-R.W. and L.-J.Z., field management and disease identification; Y.-R.L. and X.-Q.Z. experimental design, writing and polishing the manuscript. All authors read the final manuscript and approved submission.

Funding

This research was financially supported by grants from the Science Technology Department of Zhejiang Province (Project No. 2016C02053-5) and the Science Technology Bureau of Ningbo City (Project No. 2017C10001).

Acknowledgments

The authors acknowledge the Personal Biotechnology Co., Ltd. (Shanghai, China) for offering the sequencing service.

Conflicts of Interest

All of the authors declare no conflict of interest.

Abbreviations

ALD1AGD2-like defense response gene 1
DBGDe Bruijn Graph
DEGsDifferential expression genes
eggNOGEvolutionary genealogy of genes: Non-supervised Orthologous Groups
ETIEffector-triggered immunity
FDRFalse discovery rate
FPKMExpected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced
GOGene Ontology
HPLCHigh-performance liquid chromatography
HRHypersensitive response
ISRImmune system resistance
ITSInternal transcribed spacer
KAASKEGG Automatic Annotation Server
KEGGKyoto Encyclopedia of Genes and Genome
MASMolecular assisted selection
NPR1Nonexpressor of pathogenesis-related gene 1
NrNCBI non-redundant protein sequences
pp-value
PCDProgrammed cell death
PR1Pathogenesis-related gene 1
PTIPathogen-associated molecular patterns triggered immunity
qRT-PCRQuantitative reverse transcriptase polymerase chain reaction
RCWReinforcement of the cell wall
ROSReactive oxygen species
SASalicylic acid
SARSystem acquired resistance
TGATrans-activating TGA Factors

References

  1. Crous, P.W.; Gams, W.; Stalpers, J.A.; Robert, V.; Stegehuis, G. MycoBank: An online initiative to launch mycology into the 21st century. Stud. Mycol. 2004, 19–22. [Google Scholar]
  2. O’Connell, R.J.; Thon, M.R.; Hacquard, S.; Amyotte, S.G.; Kleemann, J.; Torres, M.F.; Damm, U.; Buiate, E.A.; Epstein, L.; Alkan, N.; et al. Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses. Nat. Genet. 2012, 44, 1060. [Google Scholar] [CrossRef]
  3. Cannon, P.F.; Damm, U.; Johnston, P.R.; Weir, B.S. Colletotrichum—Current status and future directions. Stud. Mycol. 2012, 73, 181–213. [Google Scholar] [CrossRef] [PubMed]
  4. Wang, Y.C.; Hao, X.Y.; Wang, L.; Xiao, B.; Wang, X.C.; Yang, Y.J. Diverse Colletotrichum species cause anthracnose of tea plants (Camellia sinensis (L.) O. Kuntze) in China. Sci. Rep. 2016, 6, 35287. [Google Scholar] [CrossRef]
  5. Liu, F.; Weir, B.S.; Damm, U.; Crous, P.W.; Wang, Y.; Liu, B.; Wang, M.; Zhang, M.; Cai, L. Unravelling Colletotrichum species associated with Camellia: Employing ApMat and GS loci to resolve species in the C. gloeosporioides complex. Persoonia 2015, 35, 63–86. [Google Scholar] [CrossRef]
  6. Wang, L.; Wang, Y.C.; Cao, H.L.; Hao, X.Y.; Zeng, J.M.; Yang, Y.J.; Wang, X.C. Transcriptome analysis of an anthracnose-resistant tea plant cultivar reveals genes associated with resistance to Colletotrichum camelliae. PLoS ONE 2016, 11, e0148535. [Google Scholar] [CrossRef]
  7. Guo, M.; Pan, Y.M.; Dai, Y.L.; Gao, Z.M. First report of brown blight disease caused by Colletotrichum gloeosporioides on Camellia sinensis in Anhui Province, China. Plant Dis. 2014, 98, 284. [Google Scholar] [CrossRef] [PubMed]
  8. Yoshida, K.; Takeda, Y. Evaluation of anthracnose resistance among tea genetic resources by wound-inoculation assay. JARQ-Jpn. Agr. Res. Quart. 2006, 40, 379–386. [Google Scholar] [CrossRef]
  9. Yamada, K.; Sonoda, R. A fluorescence microscopic study of the infection process of Discula theae-sinensis in Tea. JARQ-Jpn. Agr. Res. Quart. 2014, 48, 399–402. [Google Scholar] [CrossRef]
  10. Yoshida, K.; Ogino, A.; Yamada, K.; Sonoda, R. Induction of disease resistance in tea (Camellia sinensis L.) by plant activators. JARQ-Jpn. Agr. Res. Quart. 2010, 44, 391–398. [Google Scholar] [CrossRef]
  11. Kim, G.H.; Lim, M.T.; Hur, J.S.; Yum, K.J.; Koh, Y.J. Biological control of tea anthracnose using an antagonistic bacterium of Bacillus subtilis isolated from tea leaves. Plant Pathol. J. 2009, 25, 99–102. [Google Scholar] [CrossRef]
  12. Wang, Y.C.; Hao, X.Y.; Lu, Q.H.; Wang, L.; Qian, W.J.; Li, N.N.; Ding, C.Q.; Wang, X.C.; Yang, Y.J. Transcriptional analysis and histochemistry reveal that hypersensitive cell death and H2O2 have crucial roles in the resistance of tea plant (Camellia sinensis (L.) O. Kuntze) to anthracnose. Horticul. Res. 2018, 5, 18. [Google Scholar] [CrossRef] [PubMed]
  13. Wang, Y.C.; Qian, W.J.; Li, N.N.; Hao, X.Y.; Wang, L.; Xiao, B.; Wang, X.C.; Yang, Y.J. Metabolic changes of caffeine in tea plant (Camellia sinensis (L.) O. Kuntze) as defense response to Colletotrichum fructicola. J. Agric. Food Chem. 2016, 64, 6685–6693. [Google Scholar] [CrossRef]
  14. Lu, Q.; Wang, Y.; Li, N.; Ni, D.; Yang, Y.; Wang, X. Differences in the characteristics and pathogenicity of Colletotrichum camelliae and C. fructicola isolated from the tea plant [Camellia sinensis (L.) O. Kuntze]. Front. Microbiol. 2018, 9, 3060. [Google Scholar] [CrossRef] [PubMed]
  15. Chen, Y.J.; Qiao, W.J.; Zeng, L.; Shen, D.H.; Liu, Z.; Wang, X.S.; Tong, H.R. Characterization, pathogenicity, and phylogenetic analyses of Colletotrichum species associated with brown blight disease on Camellia sinensis in China. Plant Dis. 2017, 101, 1022–1028. [Google Scholar] [CrossRef] [PubMed]
  16. Vlot, A.C.; Dempsey, D.A.; Klessig, D.F. Salicylic acid, a multifaceted hormone to combat disease. Annu. Rev. Phytopathol. 2009, 47, 177–206. [Google Scholar] [CrossRef]
  17. An, C.F.; Mou, Z.L. Salicylic acid and its function in plant immunity. J. Integr. Plant Biol. 2011, 53, 412–428. [Google Scholar] [CrossRef]
  18. Pal, M.; Kovacs, V.; Szalai, G.; Soos, V.; Ma, X.; Liu, H.; Mei, H.; Janda, T. Salicylic acid and abiotic stress responses in rice. J. Agron. Crop. Sci. 2014, 200, 1–11. [Google Scholar] [CrossRef]
  19. Kovacs, V.; Gondor, O.K.; Szalai, G.; Darko, E.; Majlath, I.; Janda, T.; Pal, M. Synthesis and role of salicylic acid in wheat varieties with different levels of cadmium tolerance. J. Hazard. Mater. 2014, 280, 12–19. [Google Scholar] [CrossRef] [Green Version]
  20. Yalpani, N.; Enyedi, A.J.; Leon, J.; Raskin, I. Ultraviolet-Light And Ozone Stimulate accumulation of salicylic-acid, pathogenesis-related proteins and virus-resistance in tobacco. Planta 1994, 193, 372–376. [Google Scholar] [CrossRef]
  21. Zhou, N.; Tootle, T.L.; Tsui, F.; Klessig, D.F.; Glazebrook, J. PAD4 functions upstream from salicylic acid to control defense responses in Arabidopsis. Plant Cell 1998, 10, 1021–1030. [Google Scholar] [CrossRef] [PubMed]
  22. Nawrath, C.; Metraux, J.P. Salicylic acid induction-deficient mutants of Arabidopsis express PR-2 and PR-5 and accumulate high levels of camalexin after pathogen inoculation. Plant Cell 1999, 11, 1393–1404. [Google Scholar] [CrossRef] [PubMed]
  23. Shirasu, K.; Nakajima, H.; Rajasekhar, V.K.; Dixon, R.A.; Lamb, C. Salicylic acid potentiates an agonist-dependent gain control that amplifies pathogen signals in the activation of defense mechanisms. Plant Cell 1997, 9, 261–270. [Google Scholar] [CrossRef]
  24. Loake, G.; Grant, M. Salicylic acid in plant defence-the players and protagonists. Curr. Opin. Plant Biol. 2007, 10, 466–472. [Google Scholar] [CrossRef]
  25. Gaffney, T.; Friedrich, L.; Vernooij, B.; Negrotto, D.; Nye, G.; Uknes, S.; Ward, E.; Kessmann, H.; Ryals, J. Requirement of salicylic-acid for the induction of systemic acquired-resistance. Science 1993, 261, 754–756. [Google Scholar] [CrossRef]
  26. Hammerschmidt, R. Systemic acquired resistance. Adv. Bot. Res. 2009, 51, 173–222. [Google Scholar] [CrossRef]
  27. Westermann, A.J.; Gorski, S.A.; Vogel, J. Dual RNA-seq of pathogen and host. Nat. Rev. Microbiol. 2012, 10, 618–630. [Google Scholar] [CrossRef] [Green Version]
  28. Marguerat, S.; Bahler, J. RNA-seq: From technology to biology. Cell. Mol. Life Sci. 2010, 67, 569–579. [Google Scholar] [CrossRef]
  29. Torres, T.T.; Metta, M.; Ottenwalder, B.; Schlotterer, C. Gene expression profiling by massively parallel sequencing. Genome Res. 2008, 18, 172–177. [Google Scholar] [CrossRef]
  30. Zhang, J.A.; Liang, S.; Duan, J.L.; Wang, J.; Chen, S.L.; Cheng, Z.S.; Zhang, Q.; Liang, X.Q.; Li, Y.R. De novo assembly and characterisation of the transcriptome during seed development, and generation of genic-SSR markers in peanut (Arachis hypogaea L.). BMC Genom. 2012, 13, 90. [Google Scholar] [CrossRef]
  31. Grabherr, M.G.; Haas, B.J.; Yassour, M.; Levin, J.Z.; Thompson, D.A.; Amit, I.; Adiconis, X.; Fan, L.; Raychowdhury, R.; Zeng, Q.D.; et al. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotechnol. 2011, 29, 644. [Google Scholar] [CrossRef] [PubMed]
  32. Wei, C.L.; Yang, H.; Wang, S.B.; Zhao, J.; Liu, C.; Gao, L.P.; Xia, E.H.; Lu, Y.; Tai, Y.L.; She, G.B.; et al. Draft genome sequence of Camellia sinensis var. sinensis provides insights into the evolution of the tea genome and tea quality. Proc. Natl. Acad. Sci. USA 2018, 115, E4151–E4158. [Google Scholar] [CrossRef] [PubMed]
  33. Yang, Y.J.; Liang, Y.R. Clonal Tea Cultivars Records in China; Yang, Y.J., Liang, Y.R., Eds.; Shanghai Scientific and Technical Publishers: Shanghai, China, 2014; pp. 27–98. (In Chinese) [Google Scholar]
  34. Mur, L.A.J.; Kenton, P.; Atzorn, R.; Miersch, O.; Wasternack, C. The outcomes of concentration-specific interactions between salicylate and jasmonate signaling include synergy, antagonism, and oxidative stress leading to cell death. Plant Physiol. 2006, 140, 249–262. [Google Scholar] [CrossRef]
  35. Verberne, M.C.; Brouwer, N.; Delbianco, F.; Linthorst, H.J.M.; Bol, J.F.; Verpoorte, R. Method for the extraction of the volatile compound salicylic acid from tobacco leaf material. Phytochem. Anal. 2002, 13, 45–50. [Google Scholar] [CrossRef]
  36. Jones, J.D.G.; Dangl, J.L. The plant immune system. Nature 2006, 444, 323–329. [Google Scholar] [CrossRef] [Green Version]
  37. Rastogi, G.; Coaker, G.L.; Leveau, J.H.J. New insights into the structure and function of phyllosphere microbiota through high-throughput molecular approaches. FEMS Microbiol. Lett. 2013, 348, 1–10. [Google Scholar] [CrossRef] [Green Version]
  38. Liu, W.; Yuan, D.; Yin, P.; Wang, Z.H.; Guo, G.Y. Research progress on anthracnose of tea plant. Chin. J. Trop. Agric. 2016, 36, 20–26. (In Chinese) [Google Scholar] [CrossRef]
  39. Bolwell, G.P. Role of active oxygen species and NO in plant defence responses. Curr. Opin. Plant Biol. 1999, 2, 287–294. [Google Scholar] [CrossRef]
  40. Banerjee, A.K.; Mandal, A.; Chanda, D.; Chakraborti, S. Oxidant, antioxidant and physical exercise. Mol. Cell. Biochem. 2003, 253, 307–312. [Google Scholar] [CrossRef]
  41. Doke, N. Generation of superoxide anion by potato-tuber protoplasts during the hypersensitive response to hyphal wall components of phytophthora-infestans and specific-inhibition of the reaction by suppressors of hypersensitivity. Physiol. Plant Pathol. 1983, 23, 359–367. [Google Scholar] [CrossRef]
  42. Bolwell, P.P.; Page, A.; Pislewska, M.; Wojtaszek, P. Pathogenic infection and the oxidative defences in plant apoplast. Protoplasma 2001, 217, 20–32. [Google Scholar] [CrossRef]
  43. Lamb, C.; Dixon, R.A. The oxidative burst in plant disease resistance. Annu. Rev. Plant Physiol. Plant Mol. Biol. 1997, 48, 251–275. [Google Scholar] [CrossRef]
  44. Grant, J.J.; Loake, G.J. Role of reactive oxygen intermediates and cognate redox signaling in disease resistance. Plant Physiol. 2000, 124, 21–29. [Google Scholar] [CrossRef]
  45. Ellinger, D.; Naumann, M.; Falter, C.; Zwikowics, C.; Jamrow, T.; Manisseri, C.; Somerville, S.C.; Voigt, C.A. Elevated early callose deposition results in complete penetration resistance to powdery mildew in Arabidopsis. Plant Physiol. 2013, 161, 1433–1444. [Google Scholar] [CrossRef]
  46. ThordalChristensen, H.; Zhang, Z.G.; Wei, Y.D.; Collinge, D.B. Subcellular localization of H2O2 in plants. H2O2 accumulation in papillae and hypersensitive response during the barley-powdery mildew interaction. Plant J. 1997, 11, 1187–1194. [Google Scholar] [CrossRef]
  47. Denness, L.; McKenna, J.F.; Segonzac, C.; Wormit, A.; Madhou, P.; Bennett, M.; Mansfield, J.; Zipfel, C.; Hamann, T. Cell wall damage-induced lignin biosynthesis is regulated by a reactive oxygen species- and jasmonic acid-dependent process in Arabidopsis. Plant Physiol. 2011, 156, 1364–1374. [Google Scholar] [CrossRef]
  48. Bellincampi, D.; Cervone, F.; Lionetti, V. Plant cell wall dynamics and wall-related susceptibility in plant-pathogen interactions. Front. Plant Sci. 2014, 5, 228. [Google Scholar] [CrossRef]
  49. Torres, M.A.; Dangl, J.L.; Jones, J.D.G. Arabidopsis gp91(phox) homologues AtrbohD and AtrbohF are required for accumulation of reactive oxygen intermediates in the plant defense response. Proc. Natl. Acad. Sci. USA 2002, 99, 517–522. [Google Scholar] [CrossRef]
  50. Marino, D.; Dunand, C.; Puppo, A.; Pauly, N. A burst of plant NADPH oxidases. Trends Plant Sci. 2012, 17, 9–15. [Google Scholar] [CrossRef]
  51. Kadota, Y.; Sklenar, J.; Derbyshire, P.; Stransfeld, L.; Asai, S.; Ntoukakis, V.; Jones, J.D.G.; Shirasu, K.; Menke, F.; Jones, A.; et al. Direct regulation of the NADPH oxidase RBOHD by the PRR-associated kinase BIK1 during plant immunity. Mol. Cell 2014, 54, 43–55. [Google Scholar] [CrossRef]
  52. Kobayashi, M.; Ohura, I.; Kawakita, K.; Yokota, N.; Fujiwara, M.; Shimamoto, K.; Doke, N.; Yoshioka, H. Calcium-dependent protein kinases regulate the production of reactive oxygen species by potato NADPH oxidase. Plant Cell 2007, 19, 1065–1080. [Google Scholar] [CrossRef]
  53. Kim, M.C.; Panstruga, R.; Elliott, C.; Muller, J.; Devoto, A.; Yoon, H.W.; Park, H.C.; Cho, M.J.; Schulze-Lefert, P. Calmodulin interacts with MLO protein to regulate defence against mildew in barley. Nature 2002, 416, 447–450. [Google Scholar] [CrossRef]
  54. Sun, D.Y.; Bian, Y.Q.; Zhao, B.H.; Zhao, L.Y.; Yu, X.M.; Shengjun, D. The effects of extracellular calmodulin on cell-wall regeneration of protoplasts and cell-division. Plant Cell Physiol. 1995, 36, 133–138. [Google Scholar]
  55. Lecourieux, D.; Raneva, R.; Pugin, A. Calcium in plant defence-signalling pathways. New Phytol. 2006, 171, 249–269. [Google Scholar] [CrossRef] [Green Version]
  56. Moffett, P. Mechanisms of recognition in dominant r gene mediated resistance. Adv. Virus Res. 2009, 75, 1–299. [Google Scholar] [CrossRef]
  57. Hubert, D.A.; Tornero, P.; Belkhadir, Y.; Krishna, P.; Takahashi, A.; Shirasu, K.; Dangl, J.L. Cytosolic HSP90 associates with and modulates the Arabidopsis RPM1 disease resistance protein. EMBO J. 2003, 22, 5679–5689. [Google Scholar] [CrossRef]
  58. Arofatullah, N.A.; Hasegawa, M.; Tanabata, S.; Ogiwara, I.; Sato, T. Heat shock-induced resistance against pseudomonas syringae pv. tomato (Okabe) Young et al. via Heat shock transcription factors in tomato. Agronomy 2019, 9, 2. [Google Scholar] [CrossRef]
  59. Glazebrook, J. Contrasting mechanisms of defense against biotrophic and necrotrophic pathogens. Annu. Rev. Phytopathol. 2005, 43, 205–227. [Google Scholar] [CrossRef]
  60. Tsuda, K.; Mine, A.; Bethke, G.; Igarashi, D.; Botanga, C.J.; Tsuda, Y.; Glazebrook, J.; Sato, M.; Katagiri, F. Dual regulation of gene expression mediated by extended MAPK activation and salicylic acid contributes to robust innate immunity in Arabidopsis thaliana. PLoS Genet. 2013, 9, e1004015. [Google Scholar] [CrossRef]
  61. Abad, L.R.; DUrzo, M.P.; Liu, D.; Narasimhan, M.L.; Reuveni, M.; Zhu, J.K.; Niu, X.M.; Singh, N.K.; Hasegawa, P.M.; Bressan, R.A. Antifungal activity of tobacco osmotin has specificity and involves plasma membrane permeabilization. Plant Sci. 1996, 118, 11–23. [Google Scholar] [CrossRef]
  62. Song, J.T.; Lu, H.; McDowell, J.M.; Greenberg, J.T. A key role for ALD1 in activation of local and systemic defenses in Arabidopsis. Plant J. 2004, 40, 200–212. [Google Scholar] [CrossRef]
  63. Mou, Z.; Fan, W.H.; Dong, X.N. Inducers of plant systemic acquired resistance regulate NPR1 function through redox changes. Cell 2003, 113, 935–944. [Google Scholar] [CrossRef]
  64. Wu, Y.; Zhang, D.; Chu, J.Y.; Boyle, P.; Wang, Y.; Brindle, I.D.; De Luca, V.; Despres, C. The Arabidopsis NPR1 protein is a receptor for the plant defense hormone salicylic acid. Cell Rep. 2012, 1, 639–647. [Google Scholar] [CrossRef]
  65. Zhou, J.M.; Trifa, Y.; Silva, H.; Pontier, D.; Lam, E.; Shah, J.; Klessig, D.F. NPR1 differentially interacts with members of the TGA/OBF family of transcription factors that bind an element of the PR-1 gene required for induction by salicylic acid. Mol. Plant-Microbe Interact. 2000, 13, 191–202. [Google Scholar] [CrossRef]
  66. Klessig, D.F.; Durner, J.; Noad, R.; Navarre, D.A.; Wendehenne, D.; Kumar, D.; Zhou, J.M.; Shah, J.; Zhang, S.Q.; Kachroo, P.; et al. Nitric oxide and salicylic acid signaling in plant defense. Proc. Natl. Acad. Sci. USA 2000, 97, 8849–8855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Wang, D.; Amornsiripanitch, N.; Dong, X.N. A genomic approach to identify regulatory nodes in the transcriptional network of systemic acquired resistance in plants. PLoS Path. 2006, 2, 1042–1050. [Google Scholar] [CrossRef]
  68. Rushton, P.J.; Somssich, I.E.; Ringler, P.; Shen, Q.X.J. WRKY transcription factors. Trends Plant Sci. 2010, 15, 247–258. [Google Scholar] [CrossRef] [PubMed]
  69. Thomma, B.P.H.J.; Eggermont, K.; Penninckx, I.A.M.A.; Mauch-Mani, B.; Vogelsang, R.; Cammue, B.P.A.; Broekaert, W.F. Separate jasmonate-dependent and salicylate-dependent defense-response pathways in Arabidopsis are essential for resistance to distinct microbial pathogens. Proc. Natl. Acad. Sci. USA 1998, 95, 15107–15111. [Google Scholar] [CrossRef]
  70. Hartmann, M.; Zeier, T.; Bernsdorff, F.; Reichel-Deland, V.; Kim, D.; Hohmann, M.; Scholten, N.; Schuck, S.; Brautigam, A.; Holzel, T.; et al. Flavin monooxygenase-generated n-hydroxypipecolic acid is a critical element of plant systemic immunity. Cell 2018, 173, 456–469. [Google Scholar] [CrossRef] [PubMed]
  71. Brown, J.; Pirrung, M.; McCue, L.A. FQC Dashboard: Integrates FastQC results into a web-based, interactive, and extensible FASTQ quality control tool. Bioinformatics 2017, 33, 3137–3139. [Google Scholar] [CrossRef]
  72. Haas, B.J.; Papanicolaou, A.; Yassour, M.; Grabherr, M.; Blood, P.D.; Bowden, J.; Couger, M.B.; Eccles, D.; Li, B.; Lieber, M.; et al. De novo transcript sequence reconstruction from RNA-seq using the Trinity platform for reference generation and analysis. Nat. Protoc. 2013, 8, 1494–1512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Li, Q.S.; Li, X.M.; Qiao, R.Y.; Shen, E.H.; Lin, X.M.; Lu, J.L.; Ye, J.H.; Liang, Y.R.; Zheng, X.Q. De novo transcriptome assembly of fluorine accumulator tea plant Camellia sinensis with fluoride treatments. Sci. Data 2018, 5, 180194. [Google Scholar] [CrossRef] [PubMed]
  74. Camacho, C.; Coulouris, G.; Avagyan, V.; Ma, N.; Papadopoulos, J.; Bealer, K.; Madden, T.L. BLAST plus: Architecture and applications. BMC Bioinform. 2009, 10, 421. [Google Scholar] [CrossRef] [PubMed]
  75. Pruitt, K.D.; Tatusova, T.; Maglott, D.R. NCBI reference sequences (RefSeq): A curated non-redundant sequence database of genomes, transcripts and proteins. Nucleic Acids Res. 2007, 35, D61–D65. [Google Scholar] [CrossRef]
  76. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef] [PubMed]
  77. Kanehisa, M.; Goto, S.; Kawashima, S.; Okuno, Y.; Hattori, M. The KEGG resource for deciphering the genome. Nucleic Acids Res. 2004, 32, D277–D280. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Boeckmann, B.; Bairoch, A.; Apweiler, R.; Blatter, M.C.; Estreicher, A.; Gasteiger, E.; Martin, M.J.; Michoud, K.; O’Donovan, C.; Phan, I.; et al. The SWISS-PROT protein knowledgebase and its supplement TrEMBL in 2003. Nucleic Acids Res. 2003, 31, 365–370. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Powell, S.; Forslund, K.; Szklarczyk, D.; Trachana, K.; Roth, A.; Huerta-Cepas, J.; Gabaldon, T.; Rattei, T.; Creevey, C.; Kuhn, M.; et al. eggNOG v4.0: Nested orthology inference across 3686 organisms. Nucleic Acids Res. 2014, 42, D231–D239. [Google Scholar] [CrossRef] [PubMed]
  80. Gotz, S.; Garcia-Gomez, J.M.; Terol, J.; Williams, T.D.; Nagaraj, S.H.; Nueda, M.J.; Robles, M.; Talon, M.; Dopazo, J.; Conesa, A. High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36, 3420–3435. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  81. Moriya, Y.; Itoh, M.; Okuda, S.; Yoshizawa, A.C.; Kanehisa, M. KAAS: An automatic genome annotation and pathway reconstruction server. Nucleic Acids Res. 2007, 35, W182–W185. [Google Scholar] [CrossRef]
  82. Li, B.; Dewey, C.N. RSEM: Accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinform. 2011, 12, 323. [Google Scholar] [CrossRef] [PubMed]
  83. Anders, S.; Huber, W. Differential expression analysis for sequence count data. Genome Biol. 2010, 11, R106. [Google Scholar] [CrossRef] [PubMed]
  84. Li, Q.S.; Lin, X.M.; Qiao, R.Y.; Zheng, X.Q.; Lu, J.L.; Ye, J.H.; Liang, Y.R. Effect of fluoride treatment on gene expression in tea plant (Camellia sinensis). Sci. Rep. 2017, 7, 9847. [Google Scholar] [CrossRef] [PubMed]
  85. Livak, K.J.; Schmittgen, T.D. Analysis of relative gene expression data using real-time quantitative PCR and the 2(T)(-Delta Delta C) method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  86. Zawoznik, M.S.; Groppa, M.D.; Tomaro, M.L.; Benavides, M.P. Endogenous salicylic acid potentiates cadmium-induced oxidative stress in Arabidopsis thaliana. Plant Sci. 2007, 173, 190–197. [Google Scholar] [CrossRef]
Figure 1. Sample preparation and experiment scheme. (a) Anthracnose of tea plant occurrence in tea field. Yellow arrow represents healthy leaf (left) and partially enlarged drawing of anthracnose infected leaf (right). (b) A schematic diagram of a healthy leaf (left) and an infected leaf (right). The leaf slices were cut for RNA isolation as described in red boxes. Pairwise comparisons for LJ43 and ZN139 were conducted by arrows. Horizontal arrows mean comparisons between healthy leaves (LJ43_H and ZN139_H) and anthracnose disease infected leaves (LJ43_D & ZN139_D) within cultivar, and the vertical arrow indicates comparisons between cultivars.
Figure 1. Sample preparation and experiment scheme. (a) Anthracnose of tea plant occurrence in tea field. Yellow arrow represents healthy leaf (left) and partially enlarged drawing of anthracnose infected leaf (right). (b) A schematic diagram of a healthy leaf (left) and an infected leaf (right). The leaf slices were cut for RNA isolation as described in red boxes. Pairwise comparisons for LJ43 and ZN139 were conducted by arrows. Horizontal arrows mean comparisons between healthy leaves (LJ43_H and ZN139_H) and anthracnose disease infected leaves (LJ43_D & ZN139_D) within cultivar, and the vertical arrow indicates comparisons between cultivars.
Ijms 20 02439 g001
Figure 2. Overview of RNA-seq analysis workflows. (a) Acquisition of high-quality sequence (clean reads) from 8 normalized RNA samples. (b) The protocol of de novo assembly by Trinity, annotation and downstream analysis for differential expression genes (DEGs).
Figure 2. Overview of RNA-seq analysis workflows. (a) Acquisition of high-quality sequence (clean reads) from 8 normalized RNA samples. (b) The protocol of de novo assembly by Trinity, annotation and downstream analysis for differential expression genes (DEGs).
Ijms 20 02439 g002
Figure 3. Gene Ontology (GO) enrichment test of cultivars LJ43 (a) and ZN139 (b). GO terms are divided into three catalogs, i.e., molecular function (green bars), biological processes (orange bars), cellular component (blue bars), which is shown in Y-axis. X-axis shows the percentage of unigenes in all DEGs aligned to the GO database with Corrected p-value (p) < 0.05.
Figure 3. Gene Ontology (GO) enrichment test of cultivars LJ43 (a) and ZN139 (b). GO terms are divided into three catalogs, i.e., molecular function (green bars), biological processes (orange bars), cellular component (blue bars), which is shown in Y-axis. X-axis shows the percentage of unigenes in all DEGs aligned to the GO database with Corrected p-value (p) < 0.05.
Ijms 20 02439 g003
Figure 4. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment bubble chart of LJ43 (a) and ZN139 (b). The Top 25 of pathway enrichment was revealed with low false discovery rate (FDR) (in red colour, significantly) and high FDR (in green colour, insignificantly). X-axis shows the percentage of DEGs in total number of genes involved in corresponding KEGG pathway. Y-axis expresses KEGG pathway. As well, the size of bubble means the number of DEGs in this item.
Figure 4. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment bubble chart of LJ43 (a) and ZN139 (b). The Top 25 of pathway enrichment was revealed with low false discovery rate (FDR) (in red colour, significantly) and high FDR (in green colour, insignificantly). X-axis shows the percentage of DEGs in total number of genes involved in corresponding KEGG pathway. Y-axis expresses KEGG pathway. As well, the size of bubble means the number of DEGs in this item.
Ijms 20 02439 g004
Figure 5. Visualization of plant hormone signal transduction (ko04075) KEGG pathway for two tea cultivars. The DEGs comparisons of anthracnose infected leaves versus healthy leaves was depicted by boxes in different colors, with red colour representing upregulation, green colour representing downregulation and purple colour representing mixture of upregulation and downregulation. Other boxes mean no DEGs mapping to KEGG ontology terms.
Figure 5. Visualization of plant hormone signal transduction (ko04075) KEGG pathway for two tea cultivars. The DEGs comparisons of anthracnose infected leaves versus healthy leaves was depicted by boxes in different colors, with red colour representing upregulation, green colour representing downregulation and purple colour representing mixture of upregulation and downregulation. Other boxes mean no DEGs mapping to KEGG ontology terms.
Ijms 20 02439 g005
Figure 6. Quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) validation of DEGs and correlation analysis of results between RNA-seq and qRT-PCR. (a) FPKM (expected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced) values present RNA-seq expression of 5 DEGs. The sequence of ALD1 could be found in unigene ID LJ62765_g1 and ZN66762_g1, as NPR1 in LJ68588_g3 and ZN72573_g8, TGA2 in LJ58444_g1 and ZN66591_g2, TGA3 in LJ61836_g1 and ZN62926_g1, PR1 in LJ62639_g1 and ZN62871_g1. All RNA-seq experiments were performed in two biological replicates (n = 2), statistical significance p of which was generated from DEseq software. (b) The relative expression values were firstly normalized using actin as an internal reference, and then made relative to LJ43_H in which the expression value of LJ43_H was fixed as 1. All qRT-PCR experiments were performed in triplicate (n = 3). (c) Correlation analysis of results between RNA-seq and qRT-PCR. FPKM values from RNA-seq were compared to relative expression levels detected by qRT-PCR. Data are represented as mean ± standard deviation (SD). Compared with related healthy leaves, * means p < 0 05, ** means p < 0 01 and ***means p < 0 001.
Figure 6. Quantitative reverse transcriptase polymerase chain reaction (qRT-PCR) validation of DEGs and correlation analysis of results between RNA-seq and qRT-PCR. (a) FPKM (expected number of Fragments Per Kilobase of transcript sequence per Millions base pairs sequenced) values present RNA-seq expression of 5 DEGs. The sequence of ALD1 could be found in unigene ID LJ62765_g1 and ZN66762_g1, as NPR1 in LJ68588_g3 and ZN72573_g8, TGA2 in LJ58444_g1 and ZN66591_g2, TGA3 in LJ61836_g1 and ZN62926_g1, PR1 in LJ62639_g1 and ZN62871_g1. All RNA-seq experiments were performed in two biological replicates (n = 2), statistical significance p of which was generated from DEseq software. (b) The relative expression values were firstly normalized using actin as an internal reference, and then made relative to LJ43_H in which the expression value of LJ43_H was fixed as 1. All qRT-PCR experiments were performed in triplicate (n = 3). (c) Correlation analysis of results between RNA-seq and qRT-PCR. FPKM values from RNA-seq were compared to relative expression levels detected by qRT-PCR. Data are represented as mean ± standard deviation (SD). Compared with related healthy leaves, * means p < 0 05, ** means p < 0 01 and ***means p < 0 001.
Ijms 20 02439 g006
Figure 7. Salicylic acid (SA) content determined by high-performance liquid chromatography (HPLC) in LJ43. All HPLC experiments were performed in triplicate (n = 3). Data are represented as means ± SD. Compared with related healthy leaves, * means p < 0.05, ** means p < 0.01 and *** means p < 0.001.
Figure 7. Salicylic acid (SA) content determined by high-performance liquid chromatography (HPLC) in LJ43. All HPLC experiments were performed in triplicate (n = 3). Data are represented as means ± SD. Compared with related healthy leaves, * means p < 0.05, ** means p < 0.01 and *** means p < 0.001.
Ijms 20 02439 g007
Figure 8. A hypothetical model for SA accumulation and related gene expression in tea plant immunity. Solid lines indicate actual regulations of signaling effect. Dotted lines indicate hypothetical regulations of signaling effect.
Figure 8. A hypothetical model for SA accumulation and related gene expression in tea plant immunity. Solid lines indicate actual regulations of signaling effect. Dotted lines indicate hypothetical regulations of signaling effect.
Ijms 20 02439 g008
Table 1. Summary statistics of raw reads and clean reads.
Table 1. Summary statistics of raw reads and clean reads.
Sample NameLJ43_HLJ43_DZN139_HZN139_DAll Samples
Raw Reads Number (million)38.86 ± 1.1942.94 ± 0.0340.06 ± 3.5237.40 ± 2.27318.52
Total Bases (billion bp)5.83 ± 0.186.44 ± 0.006.01 ± 0.535.61 ± 0.3447.78
N (%)0.00 ± 0.000.00 ± 0.000.00 ± 0.000.00 ± 0.00
Q20 (%)95.29 ± 0.3895.43 ± 0.1395.06 ± 0.2495.32 ± 0.32
Q30 (%)89.38 ± 0.6989.67 ± 0.2788.89 ± 0.4089.47 ± 0.57
GC (%)48.72 ± 0.0149.07 ± 0.1348.79 ± 0.3048.41 ± 0.72
Clean Reads (%)98.04 ± 0.2898.11 ± 0.0597.99 ± 0.1797.99 ± 0.25
Table 2. Information of assembled contigs, unigenes and annotations for two cultivars.
Table 2. Information of assembled contigs, unigenes and annotations for two cultivars.
CultivarLJ43ZN139
ContigTotal Length (bp)102,521,451107,758,129
Sequence Number335,186352,038
Mean Length (bp)306306
N50 (bp)404405
N50 Sequence No.54,93958,050
N90 (bp)150150
N90 Sequence No.251,498264,172
UnigeneTotal Length (bp)65,163,32968,536,138
Sequence Number109,316115,953
Mean Length (bp)596591
N50 (bp)821803
N50 Sequence No.20,02021,541
N90 (bp)261260
N90 Sequence No.81,05186,082
Annotated Unigene Number in DatabaseNr45,230 (41.38%)46,383 (40%)
GO24,892 (22.77%)25,168 (21.71%)
KEGG5799 (5.3%)5873 (5.06%)
eggNOG43,025 (39.36%)44,067 (38%)
Swiss-Prot33,895 (31.01%)34,860 (30.06%)
In all database4605 (4.21%)4677 (4.03%)
Table 3. List of total DEGs in case-control study 1.
Table 3. List of total DEGs in case-control study 1.
ControlCaseUp-Regulated DEGsDown-Regulated DEGsTotal DEGs
LJ43_HLJ43_D10825391621
ZN139_HZN139_D152715623089
1 Healthy leaf were identified as control group, and anthracnose infected leaves as case group.

Share and Cite

MDPI and ACS Style

Shi, Y.-L.; Sheng, Y.-Y.; Cai, Z.-Y.; Yang, R.; Li, Q.-S.; Li, X.-M.; Li, D.; Guo, X.-Y.; Lu, J.-L.; Ye, J.-H.; et al. Involvement of Salicylic Acid in Anthracnose Infection in Tea Plants Revealed by Transcriptome Profiling. Int. J. Mol. Sci. 2019, 20, 2439. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20102439

AMA Style

Shi Y-L, Sheng Y-Y, Cai Z-Y, Yang R, Li Q-S, Li X-M, Li D, Guo X-Y, Lu J-L, Ye J-H, et al. Involvement of Salicylic Acid in Anthracnose Infection in Tea Plants Revealed by Transcriptome Profiling. International Journal of Molecular Sciences. 2019; 20(10):2439. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20102439

Chicago/Turabian Style

Shi, Yun-Long, Yue-Yue Sheng, Zhuo-Yu Cai, Rui Yang, Qing-Sheng Li, Xu-Min Li, Da Li, Xiao-Yuan Guo, Jian-Liang Lu, Jian-Hui Ye, and et al. 2019. "Involvement of Salicylic Acid in Anthracnose Infection in Tea Plants Revealed by Transcriptome Profiling" International Journal of Molecular Sciences 20, no. 10: 2439. https://0-doi-org.brum.beds.ac.uk/10.3390/ijms20102439

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

Article Metrics

Back to TopTop