Front. Plant Sci., 15 May 2017
Sec. Plant Breeding

Comparative Transcriptomics Reveals Differential Gene Expression Related to Colletotrichum gloeosporioides Resistance in the Octoploid Strawberry

Feng Wang1†, Feng Zhang2,3†, Mengmeng Chen1, Zhiheng Liu1, Zhihong Zhang2, Junfan Fu1 and Yue Ma1,2*
  • 1College of Plant Protection, Shenyang Agricultural University, Shenyang, China
  • 2College of Horticulture, Shenyang Agricultural University, Shenyang, China
  • 3Institute of Crop Science, Chinese Academy of Agricultural Sciences, Beijing, China

The strawberry is an important fruit worldwide; however, the development of the strawberry industry is limited by fungal disease. Anthracnose is caused by the pathogen Colletotrichum gloeosporioides and leads to large-scale losses in strawberry quality and production. However, the transcriptional response of strawberry to infection with C. gloeosporioides is poorly understood. In the present study, the strawberry leaf transcriptome of the ‘Yanli’ and ‘Benihoppe’ cultivars were deep sequenced via an RNA-seq analysis to study C. gloeosporioides resistance in strawberry. Among the sequences, differentially expressed genes were annotated with Gene Ontology terms and subjected to pathway enrichment analysis. Significant categories included defense, plant–pathogen interactions and flavonoid biosynthesis were identified. The comprehensive transcriptome data set provides molecular insight into C. gloeosporioides resistance genes in resistant and susceptible strawberry cultivars. Our findings can enhance breeding efforts in strawberry.


Cultivated strawberry (Fragaria × ananassa Duch) is an important fruit worldwide and is popular with consumers because of its flavor, fragrance and nutritional value. Anthracnose is one of the most destructive diseases in strawberry cultivation areas. In strawberry, anthracnose is mainly caused by three fungal pathogens: Colletotrichum gloeosporioides, Colletotrichum fragariae, and Colletotrichum acutatum (Peres et al., 2005). C. gloeosporioides typically causes petiole and stolon lesions and crown rot in strawberry and can also produce symptoms in the fruit (McInnes et al., 1992). With the expansion of strawberry planting areas in China, anthracnose is causing increasingly significant harm to China’s strawberry industry. The strawberry greenhouse planting area is gradually expanding in Northern China, and the spread of strawberry anthracnose disease is significantly impacting the development of the strawberry industry in this region.

In strawberry, significant efforts have been devoted to revealing the molecular components associated with anthracnose resistance. Subtractive hybridization (Casado-Díaz et al., 2006), microarray (Guidarelli et al., 2011), and proteomics (Fang et al., 2012) have been used to identify the genes associated with anthracnose infection. Li et al. (2013) identified a set of NB-LRR genes that display ecotype-specific responses to C. gloeosporioides inoculation via the genome-wide isolation and characterization of NB-LRR genes in woodland strawberry. Zhang et al. (2016) found that NB-LRR transcripts accumulate in response to SA pretreatment and that SA directly inhibits the germination of C. gloeosporioides. Amil-Ruiz et al. (2016) revealed that specific aspects of SA- and JA-dependent signaling pathways are activated in strawberry upon interaction with C. acutatum. However, the components of the strawberry defense network remain largely unknown or poorly understood, and the exact mechanisms have not been identified.

In this study, we performed disease assessment of the four main strawberry cultivars (‘Yanli,’ ‘Toyonoka,’ ‘Sachinoka’ and ‘Benihoppe’) in northern China following infection with C. gloeosporioides. We compared the dynamic changes in gene expression during the process of infection with C. gloeosporioides. We generated a comparative transcriptome between ‘Yanli’ and ‘Benihoppe,’ with three different libraries for each cultivar. These data deepen our understanding of the complex genetic and molecular mechanisms of strawberry defense.

Materials and Methods

Plant Materials

Plant tissue-cultured material from strawberry cultivars ‘Yanli,’ ‘Toyonoka,’ ‘Sachinoka’ and ‘Benihoppe’ was grown in a greenhouse at the Fruit Molecular Biology Laboratory of Shenyang Agricultural University. The experiment was organized into two groups, with 20 strains per variety in each group. The plants were selected for uniform size and color and an absence of visual defects. The plants were managed under identical conditions. Collected leaves were frozen immediately in liquid nitrogen and stored at -80°C until further use.

Colletotrichum gloeosporioides Culture and Infection

The pathogenic fungi used for the experiments were grown in shaking Potato Dextrose (PD) cultures for 7–10 days at a culture temperature of 25°C and then were filtered through four layers of gauze into a diluent. The accession number of the internal transcribed spacer (ITS) sequence of the C. gloeosporioides strain is KT224457. After incubation, the suspension population was counted using a bacterial counting chamber and adjusted to 106 spores ml-1. Then, the conidial suspension was transferred to a watering can, and the plant leaves and stems were sprayed to cover their surfaces with a uniform layer of the suspension. The inoculated leaves were subsequently incubated at 25°C for disease development. The results were recorded after 0, 24, 72, 96, and 120 h.

Disease Assessment

External symptoms of disease were indicated by any sign of disease spots on the plants. The disease incidence and disease index (DI) on each leaf were observed and recorded after incubation. When the spot zone beyond the wounded site on a leaf was greater than 1 mm wide, the site was scored as infected. The disease incidence was expressed as the percentage (%) of infected leaves among the total leaves. The disease index of each leaf was rated according to the percentage of observed spots on a scale of 0–5: (0) no disease spots; (1) <5% disease spot area; (2) 5–10% disease spot area; (3) 10–15% disease spot area; (4) 15–20% disease spot area; and (5) >20% disease spot area.

DI=Σ(number of diseased plants in this index×disease index)total number of leaves×maximum disease index×100               (1)

RNA Isolation and Quality Assessment

Total RNA was isolated from the plant leaves using the modified CTAB method as described by Chang et al. (2007). The RNA samples were treated with DNase (TaKaRa, Japan) for 4 h. The integrity of the RNA samples was examined using an Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, United States).

cDNA Library Preparation and Illumina Sequencing

For the RNA-seq experiments, the strawberry cultivars ‘Yanli’ (resistant) and ‘Benihoppe’ (susceptible) were used. The samples included 0-h-untreated ‘Yanli’ (R1), 72-h-untreated ‘Yanli’ (R2) 72-h-infected (with C. gloeosporioides) ‘Yanli’ (R3), 0-h-untreated ‘Benihoppe’ (S1), 72-h-untreated ‘Benihoppe’ (S2) and 72-h-infected ‘Benihoppe’ (S3). Leaves were randomly collected for each treatment. A total of six pools were submitted to the OEbiotech Company (Shanghai, China) for cDNA library preparation and sequencing. The paired-end library preparation and sequencing were performed following standard Illumina methods using a DNA sample kit (#FC-102-1002, Illumina). The cDNA library was sequenced on the HiSeqTM 2500 Illumina sequencing platform.

Mapping Reads to the Reference Genome

Clean reads were preprocessed by removing the adaptor sequences and discarding the empty and low-quality sequences. The reads were then mapped to the octoploid strawberry genome FAN_r1.11 (Hirakawa et al., 2013) using Tophat with the default parameters.

Functional Annotation

We annotated genes based on a set of sequential BLAST searches (Altschul et al., 1997) to identify the most descriptive annotation for each sequence. The genes were compared with sequences in the National Center for Biotechnology Information (NCBI) non-redundant (Nr) protein databases2 and the Swiss-Prot protein database3. The genes were also queried against the Plant Resistance Gene database4.

Digital Gene Expression Analysis

Gene expression levels in the RNA-seq analysis were measured as reads per kilobase of exon model per million mapped reads (RPKM). DESeq software5 was used to identify differentially expressed genes (DEGs) via pairwise comparisons, and the results of all statistical tests were adjusted for multiple testing with the Benjamini—Hochberg false discovery rate (FDR < 0.01). Sequences were deemed to be significantly differentially expressed if the adjusted P-value obtained by this method was <0.05 and at least a two-fold change (>1 or <-1 in log 2 ratio value) in the RPKM was observed between the two libraries.

GO and KEGG Pathway Enrichment Analysis

The functions of the DEGs were characterized using Gene Ontology (GO) terms with the agriGO database (Du et al., 2010). The second method used for DEG characterization was based on a Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis. Pathways with Q-values < 0.05 were considered significantly enriched in DEGs. The Q-value was used to determine the P-value threshold (Benjamini and Hochberg, 1995).

qRT-PCR Assay

First-strand cDNA was synthesized from one microgram of total RNA using the PrimeScript RT reagent kit with gDNA Eraser (TaKaRa, Japan) according to the manufacturer’s instructions.

To determine the expression of DEGs, quantitative reverse transcription PCR (qRT-PCR) was performed on the iQ5 Real-Time PCR Detection System (Bio-Rad, United States) using SYBR Premix Ex Taq (TaKaRa, Japan). The cycling conditions were as follows: 40 cycles of denaturation at 95°C for 30 s, annealing at 58°C for 45 s, and extension at 72°C for 30 s. A dissociation curve was generated at the end of each PCR cycle to verify that a single product had been amplified using software provided with the iQ5 Real-Time PCR Detection System. Each gene was analyzed in triplicate, and then the average threshold cycle (CT) was calculated per sample. All of the data were normalized to the level of the Fv26S internal transcript control. Relative fold changes in gene expression were calculated using the comparative Ct method (2-ΔΔCT) (Liu C. et al., 2012; Liu T. et al., 2012). Primers used in the PCR analysis are listed in Supplementary Table S3.

Statistical Analysis

The data were statistically analyzed using DPS 7.05 software. After conducting an analysis of variance (ANOVA), Duncan’s test was used to identify significant differences between mean values.


Differential Responses of Strawberry Cultivars to C. gloeosporioides

In northern China, ‘Yanli,’ ‘Toyonoka,’ ‘Sachinoka’ and ‘Benihoppe’ are the main strawberry cultivars. The different cultivars responded differently to the pathogen C. gloeosporioides. We observed the leaves of ‘Yanli,’ ‘Toyonoka,’ ‘Sachinoka’ and ‘Benihoppe’ after 72 h of infection with C. gloeosporioides. Disease symptoms appeared in the strawberry cultivars (Figure 1A). The disease incidence and index values were determined at this time point, and statistical analysis was then performed (Figures 1B,C). The disease incidence and index values were the lowest for ‘Yanli’ and the highest for ‘Benihoppe’ among the four cultivars.


FIGURE 1. Disease symptoms on the leaves of strawberry after C. gloeosporioides infection. (A) Disease symptoms on ‘Yanli,’ ‘Toyonoka,’ ‘Sachinoka’ and ‘Benihoppe’ after 72 h of C. gloeosporioides infection. (B) Disease incidence on the leaves of the strawberry cultivars after 72 h of C. gloeosporioides infection. (C) Disease index values of the strawberry cultivars after 72 h of C. gloeosporioides infection. ∗∗p < 0.05.

Transcriptomes of Resistant and Susceptible Strawberry Cultivars Infected with C. gloeosporioides

To identify the genes that are responsive to C. gloeosporioides infection in strawberry and compare gene expression patterns between resistant and susceptible cultivars, six libraries were constructed from Colletotrichum gloeosporioides-infected leaf tissues and from non-vaccinated leaf tissues as a control. A total of 16.12 G and 16.05 G reads were sequenced from the ‘Yanli’ (R1, R2, and R3) and ‘Benihoppe’ (S1, S2, and S3) cultivars, respectively, and 206433226 clean reads were obtained from the six libraries (Table 1).


TABLE 1. Number of reads based on Illumina sequencing data in libraries.

The reads were compared with the octoploid strawberry genome FAN_r1.1.6 Gene expression levels in the RNA-seq analysis were measured as RPKM (Supplementary Table S1). The analyses revealed pairwise alterations of gene expression. R3 did not cluster with R1 and R2, and S3 did not cluster with S1 and S2 (Figure 2). The heat map, drawn from all gene count data from both genome references based on Euclidean distances, depicts the relationships of all of the transcriptomes.


FIGURE 2. Global changes in gene expression in the resistant strawberry cultivar ‘Yanli’ and the susceptible strawberry cultivar ‘Benihoppe’ during C. gloeosporioides infection. Overall clustering of six samples using all normalized count data calculated by Cuffdiff. R1: 0-h-untreated ‘Yanli,’ R2: 72-h-untreated ‘Yanli,’ R3: 72-h-infected ‘Yanli’; S1: 0-h-untreated ‘Benihoppe,’ S2: 72-h-untreated ‘Benihoppe,’ S3: 72-h-infected ‘Benihoppe.’

The genes were annotated based on their similarity to known or putative sequences in web databases. Among the 230838 genes, 163033 (70.63%) had at least one significant match to a gene model according to the BLAST search (Table 2). The expressed genes were categorized according to GO and KEGG pathway analyses. All of the genes were queried against the Plant Resistance Gene database7 using BLASTp, and a total of 5227 putative R genes were identified. A list of the identified R genes is provided in Supplementary Table S2.


TABLE 2. Summary of strawberry gene annotations.

Transcriptomic Analysis of C. gloeosporioides-responsive Genes in Strawberry

To identify the C. gloeosporioides-responsive genes in resistant and susceptible strawberry cultivars, R3 was compared with R1 and R2, with a union set of 2707 DEGs as C. gloeosporioides-responsive genes in ‘Yanli’ (Figure 3A), and S3 was compared with S1 and S2, with a union set of 1354 DEGs as C. gloeosporioides-responsive genes in ‘Benihoppe’ (Figure 3B). Then, the C. gloeosporioides-responsive genes were compared between ‘Yanli’ and ‘Benihoppe’ to reveal genes that were highly conserved or divergent in response to infection (Figure 3C). In total, we identified 2481 unigenes in ‘Yanli’ (Supplementary Table S4), 1128 unigenes in ‘Benihoppe’ (Supplementary Table S5), and 226 unigenes that intersected.


FIGURE 3. Global changes in C. gloeosporioides-responsive genes in resistant and susceptible strawberry cultivars. (A) The number of common and specific DEGs in ‘Yanli.’ (B) The number of common and specific DEGs in ‘Benihoppe.’ (C) The number of common and specific C. gloeosporioides-responsive genes in ‘Yanli’ and ‘Benihoppe.’

The expression patterns of these unigenes were analyzed in ‘Yanli’ and ‘Benihoppe’ during infection with C. gloeosporioides according to their FPKM values. The unigenes responsive to C. gloeosporioides infection in ‘Yanli’ but not responsive in ‘Benihoppe’ are shown in Figure 4A; those responsive in ‘Benihoppe’ but not ‘Yanli’ are shown in Figure 4C. To functionally classify these unigenes, we observed significant enrichment of these genes in multiple GO categories. The top GO terms with unigenes in ‘Yanli’ included cell wall organization process, flavonoid biosynthetic process, cutin biosynthetic process, response to biotic stimulus process and plant-type secondary cell wall biogenesis. The top KEGG pathways with unigenes in ‘Yanli’ included Flavonoid biosynthesis, Phenylpropanoid biosynthesis, Gap junction, Plant–pathogen interaction and Pathogenic Escherichia coli infection (Figure 4B). In ‘Benihoppe,’ 1128 unigenes were enriched in GO terms including response to heat, defense response to other organism, response to stress, and response to cold. The top KEGG pathways with unigenes in ‘Benihoppe’ included Protein processing in endoplasmic reticulum, MAPK signaling pathway, and Estrogen signaling pathway (Figure 4D).


FIGURE 4. Specific C. gloeosporioides-responsive genes in ‘Yanli’ or Benihoppe. (A) Expression patterns of C. gloeosporioides-responsive genes specific to ‘Yanli’ among the six pools. (B) KEGG enrichment analyses of C. gloeosporioides-responsive genes specific to ‘Yanli.’ (C) Expression patterns of C. gloeosporioides-responsive genes specific to ‘Benihoppe’ among the six pools. (D) KEGG enrichment analyses of C. gloeosporioides-responsive genes specific to Benihoppe.’

To confirm the RNA-seq results, some genes were validated using qRT-PCR. In the plant–pathogen interaction pathways, nine strawberry genes (CERK1, EDS1, RPM1, CNGCs, CDPK, Rboh, CaMCML, WRKY33 and PR1) were up-regulated. The gene expression patterns (nine DEGs in the plant–pathogen interaction pathways) determined by qRT-PCR showed trends similar to those found in the RNA-seq data, which suggested that our transcriptome analysis was reliable (Figure 5).


FIGURE 5. Relative expression values according to the qRT-PCR analysis of the up-regulated strawberry genes in the plant–pathogen interaction pathways between mock and infected leaves at 72 h of ‘Yanli’ (A) and ‘Benihoppe’ (B). The qPCR results were generally in accordance with the gene expression profiles in the transcriptome data. p < 0.1.


This study was designed to identify strawberry genes related to defense against anthracnose. C. gloeosporioides is the primary infectious species in China. C. gloeosporioides typically causes petiole and stolon lesions and crown rot in strawberry and can also result in symptoms in the fruit (McInnes et al., 1992). ‘Yanli,’ ‘Toyonoka,’ ‘Sachinoka’ and ‘Benihoppe’ are the main strawberry cultivars in northern China. We conducted a disease assessment and found that disease symptoms appeared within 72 h, with ‘Yanli’ being a resistance accession and ‘Benihoppe’ being a susceptible accession. However, when a reference genome is available, the sequencing reads are aligned primarily by mapping onto the sequenced reference genome. The genome of F. × ananassa was firstly sequenced in 2013 (Hirakawa et al., 2013). We infected the ‘Yanli’ and ‘Benihoppe’ strawberry cultivars with C. gloeosporioides and compared their transcriptomes.

A total of 32 G of raw sequence data were generated via Illumina sequencing of strawberry leaves with or without C. gloeosporioides infection. In response to pathogen attacks, plants have evolved complex signaling and defense pathways. The transcriptomic analysis indicated that infection C. gloeosporioides elicited a wide range of responses. Many of the identified genes encode proteins with demonstrated roles in resistance and defense functions.

In plants, innate immunity is triggered via the response of PRRs to PAMPs (pathogen-associated molecular patterns), thereby providing the first line of inducible pathogen-associated molecular PTI. Chitin is one of the most common PAMPs in fungi and is highly conserved in this group (Boller and Felix, 2009; Chen and Ronald, 2011; Gust et al., 2012). CERK1 belongs to the PRR protein family. The stimulation of PRRs leads to PTI. These genes may be involved in PTI in strawberry against C. gloeosporioides.

The second class of receptors involves gene-for-gene-type interactions, which are used by the plant’s extracellular receptors to recognize pathogen virulence molecular effectors, thereby resulting in ETI. ETI is often initiated by a subset of R genes. The evolution of new plant disease-resistance genes is subject to strong selection pressure, and these genes play key roles in recognizing proteins expressed by specific pathogen-derived AVR genes and provide stable resistance for many generations (Ellis et al., 2000). In this study, 5227 R genes were identified in strawberry (Supplementary Table S2).

Plants can induce transcription factors to DNA sequences as a method of adjusting to biological and abiotic stresses (Atkinson and Urwin, 2012). Pathogen-response genes, such as WRKY, play crucial roles in plant responses. The WRKY family has been identified in the plant kingdom, and different WRKY families are found in higher plants (Ulker and Somssich, 2004). In strawberry, Casado-Díaz et al. (2006) found that the strawberry Fawrky-1 gene may be associated with the concerted gene activation cascade that occurs within the mechanisms of plant defense. In this study, certain DEGs encoding WRKY exhibited different expression patterns, and we found that WRKY33 transcripts were accumulated in response to C. gloeosporioides (Figure 4).

Plants can produce pathogenesis-related (PR) proteins in pathological-related environments. Various studies have suggested that PR genes play important roles in disease resistance. The expression level of PR1 was found to increase with pathogen infection (Palm and Medzhitov, 2009). PR1 expression increases during infection by Meloidogyne incognita in Arabidopsis thaliana (Hamamouch et al., 2011). The over-expression of CABPR1 in Nicotiana tabacum was found to be consistent with a significant increase in resistance following inoculation with pathogenic fungi (Sarowar et al., 2005). We found that PR1 expression was up-regulated in response to infection, thus supporting the SA pathway response.

In ‘Yanli,’ C. gloeosporioides-responsive unigenes were enriched in Plant–pathogen interaction pathway (Figure 4B). In this pathway, genes including CNGCs, CDPL, and Rboh were up-regulated. Calcium-mediated signaling was also activated because the ROS pathway is responsible for the hyper-sensitive response. All of the strawberry CNGC genes were up-regulated after infection, indicating that Ca2+ participates in signal transduction during early stage strawberry defense responses. Calcium is an essential second messenger in the signal transduction pathways that regulate plant responses to fungi (Zipfel, 2009). Ca2+ can activate RBOH protein, which is involved in the production of reactive oxygen species in vitro (Sagi and Fluhr, 2001). Our results indicated that RBOH was up-regulated (Figure 4), which is consistent with previous observations in wheat, cotton, cucumber, and Brassica rapa after infection by the fusarium wilt fungal pathogen (Dowd et al., 2004; Chen et al., 2015). This result indicates that higher ROS levels occurred in the strawberry leaves after infection.

Among ‘Yanli,’ ‘Toyonoka,’ ‘Sachinoka’ and ‘Benihoppe,’ ‘Yanli’ was resistant and ‘Benihoppe’ was susceptible (Figure 1). Intriguingly, the expression of most genes differed significantly between ‘Yanli’ and ‘Benihoppe.’ There were 4013 DEGs between R1 and S1, 3914 between R2 and S2, and 3895 between R3 and S3. The C. gloeosporioides-responsive genes also differed significantly (Figure 3). The expression of most genes differed significantly between ‘Yanli’ and ‘Benihoppe,’ indicating that ‘Yanli’ and ‘Benihoppe’ responded differently to the pathogen C. gloeosporioides. Furthermore, the C. gloeosporioides-responsive genes differed significantly and were enriched in different pathways, supporting the hypothesis that these cultivars may have different networks of defense mechanisms mediated by a diverse array of signaling molecules. More work is needed to understand the complex network of defense signaling pathways in strawberry. Our RNA-seq data provide molecular information relating to the defense response in strawberry.


In conclusion, RNA-seq was utilized to reveal the transcriptomic response of the strawberry cultivars ‘Yanli’ and ‘Benihoppe’ inoculated with C. gloeosporioides, the causal agent of anthracnose. The DEGs were confirmed by qRT-PCR, which demonstrated the accuracy and reliability of the RNA-seq data. GO and KEGG enrichment analyses established that plant–pathogen interaction pathways play key roles in strawberry resistance to C. gloeosporioides. In strawberry, resistance to a variety of pathogens has been reported to be mostly polygenic and quantitatively inherited, making it difficult to associate molecular markers with disease-resistance genes. Our results further our understanding of the strawberry immune system and may facilitate future disease control through biotechnological and breeding strategies.

Author Contributions

ZL, JF, and YM planned and organized the study. ZZ provided the plant material. FW and MC performed the Colletotrichum gloeosporioides culture, infection and RNA isolation and quality assessment. FZ performed qRT-PCR measurements. FW, FZ, and YM analyzed the data. FW, FZ, MC, and YM drafted the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


The authors acknowledge the financial support from the China National Natural Science Foundation (No. 31101525), the China Postdoctoral Science Foundation Grant (No. 2014M561251), and the Cultivation Plan for Youth Agricultural Science and Technology Innovative Talents of Liaoning Province (No. 2014045), Millions Talents Project Fundation of Liaoning Province (2017069).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2017.00779/full#supplementary-material

TABLE S1 | Gene IDs and expression data.

TABLE S2 | R genes.

TABLE S3 | Primers used for qRT-PCR.

TABLE S4 | Colletotrichum gloeosporioides-responsive genes specific to ‘Yanli’.

TABLE S5 | Colletotrichum gloeosporioides-responsive genes specific to ‘Benihoppe’.


  1. ^ http://strawberry-garden.kazusa.or.jp//
  2. ^ http://www.ncbi.nlm.nih.gov
  3. ^ http://www.expasy.ch/sprot
  4. ^ http://prgdb.crg.eu/
  5. ^ http://bioconductor.org/packages/release/bioc/html/DESeq.html
  6. ^ http://strawberry-garden.kazusa.or.jp/
  7. ^ http://prgdb.crg.eu/


Altschul, S. F., Madden, T. L., Schäffer, A. A., Zhang, J., Zhang, Z., Miller, W., et al. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389–3402. doi: 10.1093/nar/25.17.3389

PubMed Abstract | CrossRef Full Text | Google Scholar

Amil-Ruiz, F., Garrido-Gala, J., Gadea, J., Blanco-Portales, R., Muñoz-Mérida, A., Trelles, O., et al. (2016). Partial activation of SA- and JA-defensive pathways in strawberry upon Colletotrichum acutatum interaction. Front. Plant Sci. 7:1036. doi: 10.3389/fpls.2016.01036

PubMed Abstract | CrossRef Full Text | Google Scholar

Atkinson, N. J., and Urwin, P. E. (2012). The interaction of plant biotic and abiotic stresses: from genes to the field. J. Exp. Bot. 63, 3523–3543. doi: 10.1093/jxb/ers100

PubMed Abstract | CrossRef Full Text | Google Scholar

Benjamini, Y., and Hochberg, Y. (1995). Controlling the false discovery rate: a practical and powerful approach to multiple testing. J. R. Stat. Soc. Series B Methodol. 57, 289–300.

Google Scholar

Boller, T., and Felix, G. (2009). A renaissance of elicitors: perception of microbe-associated molecular patterns and danger signals by pattern-recognition receptors. Annu. Rev. Plant Biol. 60, 379–406. doi: 10.1146/annurev.arplant.57.032905.105346

PubMed Abstract | CrossRef Full Text | Google Scholar

Casado-Díaz, A., Encinas-Villarejo, S., de los Santos, B., Schiliró, E., Yubero-Serrano, E. M., Amil-Ruíz, F., et al. (2006). Analysis of strawberry genes differentially expressed in response to Colletotrichum infection. Physiol. Plant 128, 633–650. doi: 10.1111/j.1399-3054.2006.00798.x

CrossRef Full Text | Google Scholar

Chang, L., Zhang, Z., Yang, H., Li, H., and Dai, H. (2007). Detection of strawberry RNA and DNA viruses by RT-PCR using total nucleic acid as a template. J. Phytopathol. 155, 431–436. doi: 10.1111/j.1439-0434.2007.01254.x

CrossRef Full Text | Google Scholar

Chen, J., Pang, W., Chen, B., Zhang, C., and Piao, Z. (2015). Transcriptome analysis of Brassica rapa near-isogenic lines carrying clubroot-resistant and –susceptible alleles in response to Plasmodiophora brassicae during early infection. Front. Plant Sci. 6:1183. doi: 10.3389/fpls.2015.01183

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., and Ronald, P. C. (2011). Innate immunity in rice. Trends Plant Sci. 16, 451–459. doi: 10.1016/j.tplants.2011.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Dowd, C., Wilson, I. W., and McFadden, H. (2004). Gene expression profile changes in cotton root and hypocotyl tissues in response to infection with Fusarium oxysporum f. sp. vasinfectum. Mol. Plant Microbe Interact. 17, 654–667. doi: 10.1094/MPMI.2004.17.6.654

PubMed Abstract | CrossRef Full Text | Google Scholar

Du, Z., Zhou, X., Ling, Y., Zhang, Z., and Su, Z. (2010). agriGO: a GO analysis toolkit for the agricultural community. Nucleic Acids Res. 38, W64–W70. doi: 10.1093/nar/gkq310

PubMed Abstract | CrossRef Full Text | Google Scholar

Ellis, J., Dodds, P., and Pryor, T. (2000). The generation of plant disease resistance gene specificities. Trends Plant Sci. 5, 373–379. doi: 10.1016/S1360-1385(00)01694-0

CrossRef Full Text | Google Scholar

Fang, X., Chen, W., Xin, Y., Zhang, H., Yan, C., Yu, H., et al. (2012). Proteomic analysis of strawberry leaves infected with Colletotrichum fragariae. J. Proteomics 75, 4074–4090. doi: 10.1016/j.jprot.2012.05.022

PubMed Abstract | CrossRef Full Text | Google Scholar

Guidarelli, M., Carbone, F., Mourgues, F., Perrotta, G., Rosati, C., Bertolini, P., et al. (2011). Colletotrichum acutatum interactions with unripe and ripe strawberry fruits and differential responses at histological and transcriptional levels. Plant Pathol. 60, 685–697. doi: 10.1111/j.1365-3059.2010.02423.x

CrossRef Full Text | Google Scholar

Gust, A. A., Willmann, R., Desaki, Y., Grabherr, H. M., and Nürnberger, T. (2012). Plant Lysm proteins: modules mediating symbiosis and immunity. Trends Plant Sci. 17, 495–502. doi: 10.1016/j.tplants.2012.04.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Hamamouch, N., Li, C., Seo, P. J., Park, C. M., and Davis, E. L. (2011). Expression of Arabidopsis pathogenesis-related genes during nematode infection. Mol. Plant Pathol. 12, 355–364. doi: 10.1111/j.1364-3703.2010.00675.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirakawa, H., Shirasawa, K., Kosugi, S., Tashiro, K., Nakayama, S., Yamada, M., et al. (2013). Dissection of the octoploid strawberry genome by deep sequencing of the genomes of Fragaria species. DNA Res. 21, 169–181. doi: 10.1093/dnares/dst049

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J., Zhang, Q. Y., Gao, Z. H., Wang, F., Duan, K., Ye, Z. W., et al. (2013). Genome-wide identification and comparative expression analysis of NBS–LRR-encoding genes upon Colletotrichum gloeosporioides infection in two ecotypes of Fragaria vesca. Gene 527, 215–227. doi: 10.1016/j.gene.2013.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C., Li, C., Liang, D., Wei, Z., Zhou, S., Wang, R., et al. (2012). Differential expression of ion transporters and aquaporins in leaves may contribute to different salt tolerance in Malus species. Plant Physiol. Biochem. 58, 159–165. doi: 10.1016/j.plaphy.2012.06.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, T., Liu, Z., Song, C., Hu, Y., Han, Z., She, J., et al. (2012). Chitin-induced dimerization activates a plant immune receptor. Science 336, 1160–1164. doi: 10.1126/science.1218867

PubMed Abstract | CrossRef Full Text | Google Scholar

McInnes, T. B., Black, L. L., and Gatti, J. M. (1992). Disease-free plants for management of strawberry anthracnose crown rot. Plant Dis. 76, 260–264. doi: 10.1094/PD-76-0260

CrossRef Full Text | Google Scholar

Palm, N. W., and Medzhitov, R. (2009). Pattern recognition receptors and control of adaptive immunity. Immunol. Rev. 227, 221–233. doi: 10.1111/j.1600-065X.2008.00731.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Peres, N. A., Timmer, L. W., Adaskaveg, J. E., and Correll, J. C. (2005). Lifestyles of Colletotrichum acutatum. Plant Dis. 89, 784–796. doi: 10.1094/PD-89-0784

CrossRef Full Text | Google Scholar

Sagi, M., and Fluhr, R. (2001). Superoxide production by plant homologues of the gp91(phox) NADPH oxidase. Modulation of activity by calcium and by tobacco mosaic virus infection. Plant Physiol. 126, 1281–1290. doi: 10.1104/pp.126.3.1281

PubMed Abstract | CrossRef Full Text | Google Scholar

Sarowar, S., Kim, Y. J., Kim, E. N., Kim, K. D., Hwang, B. K., Islam, R., et al. (2005). Overexpression of a pepper basic pathogenesis-related protein 1 gene in tobacco plants enhances resistance to heavy metal and pathogen stresses. Plant Cell Rep. 24, 216–224. doi: 10.1007/s00299-005-0928-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Ulker, B., and Somssich, I. E. (2004). WRKY transcription factors: from DNA binding towards biological function. Curr. Opin. Plant Biol. 7, 491–498. doi: 10.1016/j.pbi.2004.07.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q. Y., Zhang, L. Q., Song, L. L., Duan, K., Li, N., Wang, Y. X., et al. (2016). The different interactions of Colletotrichum gloeosporioides with two strawberry varieties and the involvement of salicylic acid. Hortic. Res. 3, 16007. doi: 10.1038/hortres.2016.7

PubMed Abstract | CrossRef Full Text | Google Scholar

Zipfel, C. (2009). Early molecular events in PAMP-triggered immunity. Curr. Opin. Plant Biol. 12, 414–420. doi: 10.1016/j.pbi.2009.06.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: strawberry, Colletotrichum gloeosporioides, transcriptome, resistance gene, quantification of gene expression

Citation: Wang F, Zhang F, Chen M, Liu Z, Zhang Z, Fu J and Ma Y (2017) Comparative Transcriptomics Reveals Differential Gene Expression Related to Colletotrichum gloeosporioides Resistance in the Octoploid Strawberry. Front. Plant Sci. 8:779. doi: 10.3389/fpls.2017.00779

Received: 21 September 2016; Accepted: 25 April 2017;
Published: 15 May 2017.

Edited by:

Vasileios Fotopoulos, Cyprus University of Technology, Cyprus

Reviewed by:

Seonghee Lee, University of Florida, United States
Sotiris Tjamos, Agricultural University of Athens, Greece

Copyright © 2017 Wang, Zhang, Chen, Liu, Zhang, Fu and Ma. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yue Ma, mayuechong@126.com

These authors have contributed equally to this work and share first authorship.