- Research article
- Open Access
Comparative transcriptional profiling of orange fruit in response to the biocontrol yeast Kloeckera apiculata and its active compounds
BMC Genomics volume 17, Article number: 17 (2016)
The yeast Kloeckera apiculata strain 34–9 is an antagonist that shows biological control activity against the postharvest fungal pathogens of citrus. An antifungal compound, 2-phenylethanol (PEA), has been identified from the extract of K. apiculata. To better understand the molecular processes underlying the response of citrus fruit tissue to K. apiculata, the extract and PEA, microarray analyses were performed on navel oranges using an Affymetrix Citrus GeneChip.
As many as 801, 339 and 608 differentially expressed genes (DEGs) were identified after the application of K. apiculata, the extract and PEA, respectively. In general, K. apiculata induced the expression of defence-related genes. In addition to chitinase and β-1,3-glucanase, genes involved in ethylene (ET), jasmonic acid (JA), calcium signalling, MAPK signalling and phenylalanine metabolism were induced. In contrast, monodehydroascorbate reductase, superoxide dismutase (SOD), catalase (CAT), peroxidase (POD) and carotenoid biosynthesis genes were down-regulated. The expression profiles for the extract- and PEA-treated samples were similar to that found for yeast (sharing 57.4 % DEGs), with a significant increase in the transcript levels of defence-related genes.
This study provides a global picture of the gene expression changes in navel oranges after the application of the antagonist yeast K. apiculata, its extract and PEA. The interpretation of the DEGs revealed new insight into the molecular processes that regulate the defence responses in orange tissue.
The biological control of postharvest pathogenic fungi using microbial antagonists has recently emerged as a promising alternative to the use of synthetic fungicides [1–4]. Over 30 yeasts have been isolated and investigated for their biological control efficacy against postharvest fruit diseases. Kloeckera apiculata strain 34–9, a yeast isolated from the epiphytes of citrus roots , has been shown to suppress postharvest fungal pathogens in citrus, e.g., Penicillium digitatum and Penicillium italicum, the causal agent of green and blue mold, respectively [6, 7].
Knowledge regarding the modes of action of biological control agents is essential for developing appropriate commercial formulations and application methods to maximize the potential use of biological control. Several routes have been proposed to explain the action mechanism of biological control agents. The yeast-induced defence response of fruit has been considered a potential means to suppress infection with plant pathogens, and growing evidences have supported this point of view [1–4]. Pichia guilliermondii and Candida famata enhanced the accumulation of phytoalexins, scoparone and scopoletin in citrus wound tissues [8, 9]. Nantawanit et al.  reported that P. guilliermondii induced a defence response in chili fruit against Colletotrichum capsici. The yeast Aureobasidium pullulans can induce the accumulation of chitinase, β-1,3-glucanase and peroxidase in apple fruit . Candida saitoana induced defence responses in apple fruit . Rhodosporidium paludigenum induced resistance and defence-related responses against P. digitatum in citrus fruit . The biocontrol capability of Pichia caribbica was based on the activation of defence-related enzymes in peaches .
Reactive oxygen species (ROS) and the phytohormones signalling pathway have been shown to regulate the yeast response processes . Castoria et al.  indicated that the ability to tolerate high levels of ROS production in fruit tissue is an essential characteristic of effective yeast antagonists. Macarisin et al.  reported that yeasts on the surfaces of fruit produced H2O2 and O2 −; O2 − acted as a global regulator to activate the fruit defence responses. The application of P. membranefaciens to citrus fruit enhanced the levels of H2O2 and O2 − in the host tissue . In contrast, fungal pathogens suppress the host tissue defence responses by acidifying the fruit with organic acids, such as citric and gluconic acid [19, 20]. The acidification of the tissue might suppress host cells’ production of H2O2, and enhance the sensitivity to pathogen-produced pectolytic enzymes [21–23].
Three phytohormones, salicylic acid (SA), jasmonic acid (JA) with its derivatives (collectively called jasmonates) and ethylene (ET), have been shown to play major roles in regulating defence responses in plants [24–26]. Candida oleophila induced disease resistance by increasing the production of phytoalexin, ET biosynthesis and phenylalanine ammonia-lyase (PAL) activity . Preharvest treatment with SA and methyl jasmonate (MeJA) induced defence-related enzymes in sweet cherries . The integration of antagonistic yeast with SA or JA resulted in a remarkably improved biocontrol efficacy [29, 30]. In general, the SA-signalling pathway is believed to mediate the resistance to biotrophic pathogens, whereas the JA/ET-signalling pathway is thought to be necessary for resistance to necrotrophic pathogens . This regulation consists of the control of a complex regulatory network that connects the different pathways to enable each to assist or antagonize the others as required and fine-tune the defence response to the individual pathogen. Other plant hormones, including abscisic acid (ABA), gibberellin and auxin, act as moderators of the plant immune signalling network and have also been implicated in plant defence .
Many classes of compounds derived from Trichoderma strains, including proteins, peptaibols, oligosaccharides, low-molecular-weight compounds and small secondary metabolites, can elicit plant defence responses , such as the expression of pathogenesis-related (PR) proteins, the induction of lignification and ROS-accumulation. Several studies have reported global changes in fruit gene expression in response to adverse stress [33–37] and antagonist yeasts [38–41] by using proteomic and transcriptomic analyses. Many fruit defence response genes and proteins were identified that may increase fruit resistance; however, little is known regarding the molecular basis of functional compounds from antagonist yeast underlying the induction of host responses. An antifungal compound 2-phenylethanol (PEA) was previously identified from the extract of K. apiculata . The present study was undertaken to provide a systematic view of the citrus response to the yeast K. apiculata and its functional compounds by using an Affymetrix Citrus GeneChip.
Global changes in citrus gene expression profiles
To obtain an overall picture of the gene regulation, biocontrol yeast, the ether extract and active compound were used to treat citrus and two independent microarray analyses were performed for each treatment. To reduce experimental variation, two sets of six fruit exocarps were harvested from each treated and untreated (water control) fruit. After the removal of low-quality and internal reference probes, a total of 20,083 transcripts were reliably detected in the microarray analysis. Microarray analysis gene changes in citrus exocarp revealed as many as 801, 339 and 608 differentially expressed genes (DEGs) that showed a significant (P < 0.05) change in expression (≥1.3-fold) after 24 h of incubation with K. apiculata, the extract and PEA (Fig. 1). We further analysed these genes in subsequent experiments.
All DGEs were aligned against the Arabidopsis database by using the Citrus HarvEST software, and detailed descriptions of the sequences are shown in Additional file 1. GO categories were assigned to the 1052 DEGs using the Blast2GO program (http://www.blast2go.org). The DEGs were categorized into 22 groups based on their biological processes, as shown in Fig. 2. Response to stimulus (267), metabolic process (495), cellular process (501), pigmentation (208), biological regulation (225), multicellular organismal process (142) and developmental process (135) were the major categories. Categories based on the cellular component revealed that the responsive genes were mainly related to cell (564), cell part (564), organelle (440) and organelle part (159). With regard to molecular function, the DEGs were classified as catalytic activity (385), binding (470), transcription regulator activity (56), transporter activity (53), molecular transducer activity (40), electron carrier activity (36), enzyme regulator activity (19), antioxidant activity (11), structural molecule activity (9), translation regulator activity (5) and nutrient reservoir activity (2). Of these, the plastid and intracellular organelle were the major sub-cellular organelles involved in the citrus response.
The responsive genes were further assessed using KEGG pathway analysis (http://www.genome.jp/kegg/) (Additional file 1). The most represented pathways are phenylpropanoid biosynthesis (26), limonene and pinene degradation (18), ABC transporters (14), proteasome (3), lysosome (5), oxidative phosphorylation (9), flavonoid biosynthesis (20), the regulation of autophagy (3), calcium signalling pathway (13), apoptosis (18), fatty acid metabolism (7), MAPK signalling pathway (25), phenylalanine, tyrosine and tryptophan biosynthesis (7), citrate cycle (TCA cycle) (5), flavone and flavonol biosynthesis (14), starch and sucrose metabolism (13), arachidonic acid metabolism (4), phenylalanine metabolism (7), ascorbate and aldarate metabolism (6) and carotenoid biosynthesis (5). Most of these pathways were consistent with biological processes that were already identified by GO analysis. Some of these pathways were related to the defence response based on previous knowledge, such as phenylpropanoid biosynthesis and the calcium signalling pathway [26, 31, 32, 36].
Change pattern of gene expression in citrus in response to K. apiculata
Of the 801 DEGs in orange exocarp tissue treatment with K. apiculata, 56 % of the annotated genes were down-regulated and 44 % were up-regulated. Furthermore, the microarray data for the probes of the significant dataset were mapped to Arabidopsis using the MapManBin software (http://ppdb.tc.cornell.edu/dbsearch/searchacc.aspx). The data obtained from this analysis are presented in Fig. 3 and Additional file 2. Major and minor CHO metabolism (10), glycolysis (3), fermentation (1), TCA (3), mitochondrial electron transport/ATP synthesis (3), lipid metabolism (17), amino acid metabolism (12), redox (6), nucleotide metabolism (9), DNA (10) and proteins (82) associated with the cell (11) showed more down-regulated than up-regulated genes, while more up-regulated genes were found in PS (5), N-metabolism (2), hormone metabolism (17) and secondary metabolism (14) associated with the cell wall (9) in response to K. apiculata treatment. A number of new genes that are potentially related to defence responses were identified in this study. Based on microarray and previous data [26, 31, 32, 36], hormone, reactive oxygen species (ROS), lipid, secondary metabolite, cell wall, stress, phenyalanine metabolism related genes were selected for further analysis. Figure 4 and Additional file 3 summarize the changes in these defence-related genes.
The first noticeable pathway is the hormone metabolism pathway. In total, 17 differentially expressed genes are involved in hormone metabolism, including the ethylene (ET)-signalling pathway of eight ethylene response factor (ERF) genes (Cit.18086.1.S1_at, Cit.22763.1.S1_s_at, Cit.2675.1.S1_s_at, Cit.2677.1.S1_at, Cit.17142.1.S1_s_at, Cit.18673.1.S1_at, Cit.20640.1.S1_at, Cit.16845.1.S1_at); jasmonic acid (JA)-signalling pathway of one hydroperoxide lyase (HPL) (Cit.10444.1.S1_at) and five allene oxide synthase (AOS) genes (Cit.905.1.S1_at, Cit.6011.1.S1_at, Cit.23585.1.S1_at, Cit.31140.1.S1_at, Cit.996.1.S1_s_at) and five abscisic acid (ABA)-signalling pathway genes (Cit.13424.1.S1_at, Cit.5225.1.S1_at, Cit.10675.1.S1_at, Cit.13166.1.S1_at, Cit.34429.1.S1_s_at) (Fig. 4 and Additional file 3). Most of these genes were highly expressed in yeast-treated orange tissue, such as ERF (Cit.2677.1.S1_at), which was up-regulated 4.0-fold according to the microarray data. In addition, two polyamine (polyamine oxidase), three auxin-responsive and five gibberellic acid (GA; gibberellin receptor, gibberellin oxidase) genes were down-regulated.
Reactive oxygen species (ROS) accumulation has been well studied for biocontrol yeast-induced defence responses in fruits [15–18]. The second group of metabolic pathways is involved in the redox and antioxidation pathway. In total, five genes involved in antioxidant biosynthesis were down-regulated, including monodehydroascorbate reductase (Cit.3320.1.S1_s_at, Cit.3318.1.S1_at), superoxide dismutase (SOD, Cit.5267.1.S1_at), catalase (CAT, Cit.8351.1.S1_s_at) and peroxidase (POD, Cit.8515.1.S1_s_at) (Fig. 4 and Additional file 3). For example, SOD (Cit.5267.1.S1_at) was down-regulated 5.5-fold according to the microarray data, and the qRT-PCR data were consistent with these results, demonstrating that the level of SOD was 2.2-times lower in response to K. apiculata-treatment than in CK. Moreover, Cytochrome P450 plays an important role in the redox pathway and has been well characterized [43–45]. Over nine different cytochrome P450 genes were detected in our microarray data, such as monooxygenase/p-coumarate 3-hydroxylase, monooxygenase 83B1 and ent-kaurenoate oxidase.
The third group of metabolic pathways consists of signalling pathway and pathogenesis-related (PR) proteins. The signalling pathway genes included 13 genes for calcium and 25 genes for MAPK signalling, most of which were up-regulated in response to K. apiculata application. The other genes encoding for chitinase (Cit.15242.1.S1_at, 1.8-fold) and β-1,3-glucanase (Cit.10558.1.S1_s_at, 1.4-fold) were stimulated by K. apiculata application (Fig. 4 and Additional file 3). In addition, five different disease resistance protein genes (TIR-NBS-LRR class) were also found in our study.
The fourth group of significant K. apiculata-responsive genes included secondary metabolic processes, such as the phenylpropanoid pathway, limonene and pinene degradation, flavone and flavonol biosynthesis and carotenoid biosynthesis (Fig. 3). These genes were induced following K. apiculata application. A significant increase in the expression of the genes encoding for chalcone-flavanone isomerase (Cit.17011.1.S1_s_at), cinnamoyl-CoA reductase (Cit.13313.1.S1_s_at), violaxanthin de-epoxidase (Cit.30844.1.S1_s_at) and shikimate 5-dehydrogenase (Cit.25466.1.S1_at) was observed. These genes are mainly involved in to lignin and flavanol biosynthesis. Moreover, genes involved in the biosynthesis of the other secondary metabolites, such as carotenoid and terpenes, were down-regulated, including p-coumarate 3-hydroxylase (Cit.30567.1.S1_at), 3-chloroallyl aldehyde dehydrogenase (Cit.30574.1.S1_s_at), ent-kaurenoate oxidase (Cit.13587.1.S1_at) and carotenoid isomerase (Cit.29769.1.S1_s_at).
Several families of transcription factors, including the WRKY, R2R3-MYB, bHLH (basic helix-loop-helix) and WD40 genes, showed significant transcriptional changes in response to K. apiculata application, as revealed by the microarray data.
Comparative analysis of gene expression in citrus between different treatment
To further analyze the response of the orange exocarp tissue to the extract and PEA, a total of 339 and 608 genes were identified. The expression profiles in response to the extract and PEA were similar to that found for K. apiculata; 57.4 % of the 803 DEGs in response to K. apiculata were also altered in the extract/PEA treatments (Fig. 1). The common up-regulated and down-regulated genes in citrus between K. apiculata-PEA and special in K. apiculata were listed in Additional file 4. The distribution of the genes among the GO and KEGG functional categories indicated that a large number of defence-related genes were also included in these DEGs. Orange exocarp tissue responded similarly to K. apiculata, the extract and PEA (Figs. 3 and 4).
Verification of microarray data by qRT-PCR analyses
To confirm that the DEGs identified by the microarray gene expression were indeed differentially expressed, 20 genes were selected based on their biological significance for confirmation in a biologically independent experiment using qRT-PCR, including SOD-, ET-, JA- and ABA-related genes, which were detected in the microarray data and bioinformatic analyses. The relative transcript abundance patterns for the stress application were compared using the transcriptome data. The results of the qRT-PCR experiments revealed that most of the genes showed the same expression pattern as the microarray data (Fig. 5), such as Cit.2677.1.S1_at, which was 3.2-times higher in response to K. apiculata-treatment than in CK in qRT-PCR data, and 4.0-times higher for the microarray data.
Pathogenesis-related (PR) proteins activity
Of the PR proteins, chitinase and β-1,3-glucanase are two of the most fully characterized enzymes that are capable of hydrolyzing the polymers of fungal cell walls  Furthermore, the accumulation of chitinase and β-1,3-glucanase is important in retarding fungal growth and decreasing the spoilage of fruits caused by fungal pathogens. The level of chitinase and β-1,3-glucanase was observed to be significantly higher after inoculation of fruit with yeast, the extract and PEA compared with the control (Fig. 6ab), which is consistent with the microarray data. In addition, some PR-proteins could act in cell wall reinforcement by catalyzing lignification, such as PR-9 . As illustrated in Fig. 7, the application of K. apiculata and the extract to discs of orange peel increased the lignin content in fruit peel relative to the water control peel discs 60 h after application.
Hydrogen peroxide (H2O2) level in orange tissue
Reactive oxygen species (ROS) burst has been shown to regulate the yeast response processes . Our study showed that K. apiculata treatment resulted in a high level of intracellular H2O2 when applied to oranges (Fig. 6c); this level decreased dramatically 12 h after the application of yeast to citrus fruit, although the statistic analysis showed that there were significant differences (P < 0.05) between K. apiculata-treatment and control at the point of 24 h, 36 h and 48 h. The extract and PEA did not enhance the level of H2O2.
Polyamine level in orange tissue
Polyamines, mainly diamine putrescine (Put), triamine spermidine (Spd) and tetraamine spermine (Spm), act as an important source of H2O2 production and have been suggested to be involved in the response to pathogen attack or responsible for enhanced disease resistance in higher plants . In yeast-treated citrus, the level of Put, Spd and Spm were observed to be lower than normal control, especially at 24 h (Fig. 6). In the extract and PEA-treated citrus, there was not a rule can be followed for Spm, Put and Spd (Fig. 6d, e, f).
Interactions between postharvest yeast biological control agents and host tissue have recently been widely studied in various fruits, including citrus, apples, pears; however, the molecular mechanisms are poorly understood. Gene expression profiling via the use of microarrays has been recognized as a powerful approach to obtain an overall view of gene expression and the physiological processes involved in the response to a particular stimulus [48–50]. In this study, we used a microarray to identify global changes in gene expression that occur in orange fruit exocarp tissues following the application of the yeast biological control agent K. apiculata, its extract and PEA. A large number of newly discovered and interesting genes encoding transcription and post-transcription factors were included in these DEGs, indicating that these genes may be key regulators that control the defence response by activating or repressing numerous genes. Additionally, a number of putative homologs of genes for host resistance were also found.
A variety of genes involved in the response to biotic and abiotic stresses, signalling, defence, hormones and secondary metabolism were identified in K. apiculata-treated orange exocarp tissue (Fig. 3). These findings imply that complex biochemical and molecular processes are involved in the reaction of fruit host tissue to the introduction of yeast cells, which have the potential to influence the efficacy of the biocontrol agent.
Plant ROS-signalling pathways have been shown to play essential roles in the regulation of host defence response processes [40, 51]. Previous data showed that the production of ROS by yeast antagonists may serve as a signal to trigger an oxidative burst in host tissue, leading to the activation of host defence mechanisms . In contrast, fungal pathogens, such as Penicillium expansum and P. digitatum, suppress host cell defence responses by inhibiting the production of H2O2 in host cells [19, 20]. In the current study, a significant accumulation of H2O2 was observed in the host tissue after the application of K. apiculata to cells at 12 h (Fig. 6c). These finding were consistent with previous reports [17, 41]. The intensity of ROS production peaked shortly after the application of yeast cells to intact fruit, with a concomitant accumulation of H2O2 in the host fruit tissue itself . Overall, these results support the notion that the intense production of ROS in fruit tissue induced by yeast cells plays a major role in the early stages of the application of K. apiculata.
Plant mitochondria have been reported as a source of the oxidative burst . In our microarray data, 39 genes were finally categorized as mitochondrial by GO categories, including respiratory chain complex and electron transport. In addition, MapManBin identified three and six genes as being involved in mitochondrial electron transport/ATP synthesis and redox ((reduction-oxidation) reactions, respectively. The oxidative burst in fruit host tissue is likely responsible for the disordered of energy metabolism in mitochondria. ROS accumulation in yeast-treated tissue is also accompanied by a decrease in the expression levels of genes encoding for ROS-detoxifying enzymes, including monodehydroascorbate reductase, SOD, CAT and POD (Fig. 4). Changes in antioxidant gene expression, which lead to an increase in the ROS levels, and the activation of defence mechanisms have been supported by several reports [15, 18, 40, 51]. ROS contribute to the activation of plant defence by inducing changes in gene expression, including the redox regulation of transcription factors, production of PR proteins, ET synthesis and cell death .
Plant hormones play pivotal roles in the regulation of the defence signalling network . The signalling pathways crosstalk in an antagonistic or synergistic manner, providing the plant with a powerful capacity to finely regulate its immune response . SA, Jas and ET are recognized as major defence hormones . Other hormones, including ABA, auxin and gibberellins, affect the SA-JA-ET backbone of the plants immune signalling network, resulting in positive or negative effects on biotrophic and necrotrophic pathogens [56–58]. In our microarray data, we detected a significant increase in JA-signalling (8 genes) and ET-signalling (6 genes) gene expression in yeast-treated orange exocarp tissue. This finding was further supported by qRT-PCR data (Fig. 5), and the decreased Spd and Put levels (Fig. 6df) also supported these results. ET and polyamines have a common precursor, and they appear to have opposing physiological roles . In addition, our results identified a difference in the expression of 27 genes related to other hormones (polyamine, gibberellins, auxins and ABA). Their signalling pathways may have indirect effects on plant immunity by antagonistically or synergistically interacting with the SA-JA-ET backbone of the plant immune signalling network . The data supported the K. apiculata-induced citrus defence response via the JA/ET-signalling pathway.
Following K. apiculata treatment, 57.4 % of DEGs showed the same pattern of change as was found following the extract/PEA treatments (Fig. 1 and Additional file 4). This result indicated that the extract and PEA have the potential to influence the efficacy of the biocontrol agent of K. apiculata by inducing defence response genes. Furthermore, the ability of K. apiculata to induce the defence response of citrus is partly related to the extract and PEA. The global expression profiles were similar in response to treatment with K. apiculata, the expression of defence-associated genes being greatly enhanced under the extract- and PEA-treatment (Figs. 3 and 4). Meanwhile, we noticed that application of the yeast strain enhanced larger effects than application of the extract or PEA, which may be interpreted the effects of K. apiculata involved in both biotic and abiotic stress. Only abiotic stress of K. apiculata was partly achieved via secretion secondary metabolite, such as PEA.
As a metabolite of L-Phe, PEA may negatively regulate the biosynthetic pathway and indirectly influence the production of a fruit’s L-Phe-derived metabolites . These metabolites have protective and regulatory functions in plants and can be categorized into three broad groups: phytoalexins (flavonoids, isoflavanones), phytoanticipins and signalling molecules (e.g., SA). Flavonoid glycosides serve as potential modulators of cell division, while flavonoids serve as regulators of auxin transport and SA acts as a regulator of both local and systemic pathogen-induced defence gene activation, the oxidative burst and pathogen-induced cell death [56, 61]. All L-Phe-related genes in the K. apiculata and PEA treatments (Cit.16303.1.S1_at, Cit.1280.1.S1_s_at, Cit.29769.1.S1_s_at, Cit.6742.1.S1_s_at, Cit.9171.1.S1_at, Cit.9944.1.S1_at, Cit.9944.1.S1_x_at, Cit.17011.1.S1_s_at, Cit.12979.1.S1_at) shared the same up- and down-regulated pattern as the control, except for Cit.15355.1.S1_at.
Our study provides a global picture of the gene expression changes that result from the application of the yeast biological control agent K. apiculata, its extract and PEA on citrus fruit (Fig. 8). The microarray data revealed a large number of genes that were reported to be involved in the defence response. The functional categorization of the DEGs revealed the involvement of a number of important pathways, including oxidative phosphorylation, phenylpropanoid biosynthesis, mitochondrial electron transport/ATP synthesis, MAPK signalling, calcium signalling and hormone cross communication, in regulating of the defence process. In addition, similar global expression profiles were acquired with the expression of defence-associated genes between K. apiculata, the extract and PEA treatments.
Fruit material and biocontrol agent
Olinda Valencia oranges (Citrus sinensis L. Osbeck) were harvested at commercial maturity (28 April) from adult trees grown in Yichang City, Hubei Province, China. Fruits without physical injuries and infections were selected based on uniformity in size. Prior to use, the fruits were disinfected with 2 % (v/v) NaOCl solution for 2 min, rinsed with tap water and air-dried. A strain of K. apiculata 34–9 was isolated from the epiphytes of citrus roots . The strain was grown in YPD medium (1 % yeast extract, 2 % peptone, 2 % dextrose and 2 % agar).
Citrus RNA extraction and microarray analysis
The extract was obtained from the cell-free culture of K. apiculata as described previously . Citrus fruits soaked for 5 min in 1.0 × 108 cells/mL K. apiculata (KA), 1530 μg/mL PEA, the extract (1000 × dilute) and water control (CK), then air dried and the samples were then placed on plastic cases for 24 h. The fresh exocarp of citrus was separated with a knife after washing the fruit with water, after which it was directly frozen in liquid nitrogen and stored at −80 °C. Each sample consisted of the pooled exocarp of six fruits. Two biological replicates were used for each line.
The total RNA was extracted as described previously . The Affymetrix GeneChip One-cycle Target Labeling Kit (Affymetrix, Santa Clara, CA; http://www.affymetrix.com/) was used for expression analysis. The GeneChip Citrus Genome Array (platform: GPL5731) contains 30,171 probe sets representing up to 33,879 citrus transcripts based on EST sequences obtained from several citrus species and citrus hybrids. The arrays were performed according to the manufacturer’s recommended protocols. Microarray experiments were designed to comply with MIAME guidelines . The differentially expressed genes (DEGs) were selected and functionally annotated as described in Gallego-Giraldo et al. . We used the classical ttest to identify DEGs and defined p-value < 0.05 to be statistically significant. The details of the citrus cDNA microarray data were submitted to NCBI under GEO accession numbers GSE45680.
Quantitative real-time RT-PCR (qRT-PCR)
The total RNA was reverse transcribed into first-strand cDNA using the M-MLV first strand kit (Invitrogen, USA) according to the manufacturer’s instructions. Twenty genes were chosen for confirmation by qRT-PCR with SYBR®GREEN Master Mix (Toyobo, Osaka, Japan). Primers for the chosen genes were designed with the primer express software (Applied Biosystems, USA) and are presented in Table 1. A qRT-PCR assay for gene expression analysis was performed on a Roche 480 Real-time PCR System (Roche Molecular Systems, Belleville, USA) using the beta-actin (ACTB) gene as an endogenous control according to Yan et al. . Briefly, the primers for the target gene and ACTB were diluted in the SYBR Mix, and 10 μL of the reaction mix were added to each well. The reactions were performed with an initial incubation at 50 °C for 2 min and at 95 °C for 1 min followed by 40 cycles of 95 °C for 15 s, 60 °C for 20 s and 72 °C 10 s. The levels of gene expression were analysed with a LightCycler®480. Zero template controls were included for each primer pair. Each PCR reaction was carried out in triplicate for three repeats, and the data were presented as the means ± SD.
Plant endogenous H2O2, lignin and enzyme activity analyses
The disinfected fruits were inoculated on their circumference using an inoculating needle (5.0 mm). For inoculations, a 10 μL aliquot of the yeast suspension at a concentration of 1.0 × 108 cells/mL was dropped onto each prick. After air drying, the fruits were stored in enclosed plastic trays to maintain a high humidity (approximately 95 %). The plastic cases were maintained at 28 °C for the indicated periods. To measure the elicitation effect, the tissue surrounding each wound of fruit was collected at hour 0, 12, 24, 36, 48, 60 after treatment, and immediately immersed in liquid nitrogen and stored at −80 °C until use. A 10-g sample (fresh weight; FW) of the exocarp was ground into a powder in liquid nitrogen.
The concentration of H2O2 was assayed using H2O2 assay kits (Nanjing Jiancheng Bioengineering Institute, China) according to the manufacturer’s instructions. The enzyme activities were determined by a Shimadzu UV-1800 spectrophotometer (Shimadzu, Japan). The activities of chitinase and β-1,3-glucanase were measured as described previously . The lignin content was quantified using the method described by Syros et al. .
Quantification of free polyamines by high-performance liquid chromatography (HPLC)
The free polyamines were quantified using a method previously described in Liu and Moriguchi  and Fu et al. . Samples were prepared and collected as described above. A sample of peel tissues (0.5 g) was homogenized in 5 mL of 5 % cold perchloric acid (PCA) for 30 min on ice. The supernatant was transferred to a new tube after centrifugation at 12000 rpm (4 °C) for 15 min; the resulting pellet was reconstituted with 5 mL of 5 % PCA and maintained on ice for 30 min before centrifugation at the same conditions. The supernatant was mixed, and 500 μL of it was benzoylated. The supernatant was mixed with 10 mL of benzoyl chloride and 1 mL of 2 mol NaOH. The resultant solution was vortexed for 30 s and then incubated for 25 min in a water bath at 37 °C. The benzoylated polyamines were then leached with 2 mL of ethyl ether, vacuum dried in a concentrator (Eppendorf 5301, Germany) and re-dissolved with 100 μL of methanol (HPLC grade). The benzoyl-polyamines (20 μL) were analysed using an Agilent 1200 HPLC systems (Santa Clara, CA, USA) equipped with a C18 reversed phase column (4.6 mm × 150 mm, particle size 5 μm) and a UV-detector according to Shi et al.  with minor modification. The column was eluted at 1 mL/min, with a programmed gradient of solvents (methanol/water), changing from 60 to 95 % in 23 min. Chromatograms were scanned at 230 nm. The polyamines were quantified in triplicate.
All the statistical analyses in this study were conducted using the Statistical Program SPSS 13.0 for windows (SPSS Inc, Chicago, IL). Analysis of variance (ANOVA) was performed and Duncan’s multiple range test was used for means separation. The statistical significance in this experiment is all applied at the level P < 0.05.
Janisiewicz WJ, Korsten L. Biological control of postharvest diseases of fruits. Annu Rev Phytopathol. 2002;40:411–41.
Droby S, Wisniewski M, Macarisin D, Wilson C. Twenty years of postharvest biocontrol research: is it time for a new paradigm? Postharvest Biol Technol. 2009;52:137–45.
Sharma RR, Singh D, Singh R. Biological control of postharvest diseases of fruits and vegetables by microbial antagonists: a review. Biol Control. 2009;50:205–21.
Jamalizadeh M, Etebarian HR, Aminian H, Alizadeh A. A review of mechanisms of action of biological control organisms against post-harvest fruit spoilage. EPPO Bulletin. 2011;41:65–71.
Long CA, Zhang W, Deng BX. Biological control of Penicillium italicum of citrus and Botrytis cinerea of grape by strain 34–9 of Kloeckera apiculata. Eur Food Res Technol. 2005;221:197–201.
González-Candelas L, Alamar S, Sánchez-Torres P, Zacarías L, Marcos JF. A transcriptomic approach highlights induction of secondary metabolism in citrus fruit in response to Penicillium digitatum infection. BMC Plant Biol. 2010;10:194.
Liu P, Luo L, Long CA. Characterization of competition for nutrients in the biocontrol of Penicillium italicum by Kloeckera apiculata. Biol Control. 2013;67:157–62.
Arras G. Mode of action of an isolate of Candida famata in biological control of Penicillium digitatum in orange fruits. Postharvest Biol Technol. 1999;8:191–8.
Rodov V, Ben-Yehoshua S, D’hallewin S, Castia T. Accumulation of phytoalexins scoparone and scolopetin in citrus fruits subjected to various postharvest treatments. Acta Hortic. 1994;381:517–23.
Nantawanit N, Chanchaichaovivat A, Panijpan B, Ruenwongsa P. Induction of defense response against Colletotrichum capsici in chili fruit by the yeast Pichia guilliermondii strain R13. Biol Control. 2010;52:145–52.
Ippolito A, EI Ghaouth A, Wilson CL, Wisniewski M. Control of postharvest decay of apple fruit by Aureobasidium pullulans and induction of defense responses. Postharvest Biol Techol. 2000;19:265–72.
El Ghaouth A, Wilson CL, Wisniewski M. Control of postharvest decay of apple fruit with Candida saitoana and induction of defense responses. Phytopathology. 2003;93:344–8.
Lu LF, Lu HP, Wu CQ, Fang WW, Yu C, Ye CZ. Rhodosporidium paludigenum induces resistance and defense-related responses against Penicillium digitatum in citrus fruit. Postharvest Biol Techol. 2013;85:196–202.
Xu BT, Zhang HY, Chen KP, Xu Q, Yao Y, Gao H. Biocontrol of postharvest Rhizopus decay of peaches with Pichia caribbica. Curr Microbiol. 2013;67:255–61.
Tian SP, Qin GZ, Li BQ. Reactive oxygen species involved in regulating fruit senescence and fungal pathogenicity. Plant Mol Biol. 2013;82:593–602.
Castoria R, Caputo L, De Curtis F, De Cicco V. Resistance of postharvest biocontrol yeasts to oxidative stress: a possible new mechanism of action. Phytopathology. 2003;93:564–72.
Macarisin D, Droby S, Bauchan G, Wisniewski M. Superoxide anion and hydrogen peroxide in the yeast antagonist–fruit interaction: a new role for reactive oxygen species in postharvest biocontrol? Postharvest Biol Technol. 2010;58:194–202.
Luo Y, Zhou YH, Zeng KF. Effect of Pichia membranaefaciens on ROS metabolism and postharvest disease control in citrus fruit. Crop Prot. 2013;53:96–102.
Prusky D, McEvoy JL, Saftner R, Conway WS, Jones R. The relationship between host acidification and virulence of Penicillium spp. on apple and citrus fruit. Phytopathology. 2004;94:44–51.
Barad S, Horowitz SB, Moscovitz O, Lichter A, Sherman A, Prusky D. A Penicillium expansum glucoseoxidase–encoding gene, GOX2, is essential for gluconic acid production and acidification during colonization of deciduous fruit. Mol Plant-Microbe Interact. 2012;25:779–88.
Eshel D, Miyara I, Ailing T, Dinoor A, Prusky D. pH regulates endoglucanase expression and virulence of Alternaria alternata in persimmon fruit. Mol Plant-Microbe Interact. 2002;15:774–9.
Macarisin D, Cohen L, Eick A, Rafael G, Belausov E, Wisniewski M, et al. Penicillium digitatum suppresses production of hydrogen peroxide in host tissue during infection of citrus fruit. Phytopathology. 2007;97:1491–500.
Miyara I, Shafran H, Davidzon M, Sherman A, Prusky D. pH regulation of ammonia secretion by Colletotrichum gloeosporioides and its effect on appressorium formation and pathogenicity. Mol Plant-Microbe Interact. 2010;23:304–16.
Mauch-Mani B, Mauch F. The role of abscisic acid in plant–pathogen interactions. Curr Opin Plant Biol. 2005;8:409–14.
Zhao Y, Wei T, Yin KQ, Chen Z, Gu H, Qu L, et al. Arabidopsis RAP2.2 plays an important role in plant resistance to Botrytis cinerea and ethylene responses. New Phytol. 2012;195:450–60.
Koornneef A, Pieterse CMJ. Cross talk in defense signaling. Plant Physiol. 2008;146:839–44.
Droby S, Vinokur V, Weiss B, Cohen L, Daus A, Golaschmidt EE, et al. Induction of resistance to Penicillium digitatum in grapefruit by the yeast biocontrol agent Candida oleophila. Phytopathology. 2002;92:393–9.
Yao HJ, Tian SP. Effects of pre- and post-harvest application of salicylic acid or methyl jasmonate on inducing disease resistance of sweet cherry fruit in storage. Postharvest Biol Tech. 2005;35:253–62.
Yu T, Chen JS, Chen RL, Huang B, Liu DH, Zheng XD. Biocontrol of blue and gray mold diseases of pear fruit by integration of antagonistic yeast with salicylic acid. Int J Food Microbiol. 2007;116:339–45.
Cao SF, Zheng YH, Wang KT, Tang SS, Rui HJ. Effect of yeast antagonist in combination with methyl jasmonate treatment on postharvest anthracnose rot of loquat fruit. Biol Control. 2009;50:73–7.
Gimenez-Ibanez S, Solano R. Nuclear jasmonate and salicylate signaling and crosstalk in defense against pathogens. Front Plant Sci. 2013;4:72.
Shoresh M, Harman GE, Mastouri F. Induced systemic resistance and plant responses to fungal biocontrol agents. Annu Rev Phytopathol. 2010;48:21–43.
Zhu A, Li W, Ye J, Sun X, Ding Y, Cheng Y, et al. Microarray expression profiling of postharvest Ponkan mandarin (Citrus reticulata) fruit under cold storage reveals regulatory gene candidates and implications on soluble sugars metabolism. J Integr Plant Biol. 2011;53:358–74.
Qin GZ, Wang YY, Cao BH, Wang WH, Tian SP. Unraveling the regulatory network of the MADS box transcription factor RIN in fruit ripening. Plant J. 2012;70:243–55.
Yu K, Xu Q, Da X, Guo F, Ding Y, Deng X. Transcriptome changes during fruit development and ripening of sweet orange (Citrus sinensis). BMC Genomics. 2012;13:10.
Yun Z, Gao HJ, Liu P, Liu SZ, Luo T, Jin S, et al. Comparative proteomic and metabolomic profiling of citrus fruit with enhancement of disease resistance by postharvest heat treatment. BMC Plant Biol. 2013;13:44.
Zheng ZL, Zhao YH. Transcriptome comparison and gene coexpression network analysis provide a systems view of citrus response to ‘Candidatus Liberibacter asiaticus’ infection. BMC Genomics. 2013;14:27.
Chan Z, Qin G, Xu X, Li B, Tian S. Proteome approach to characterize proteins induces by antagonist yeast and salicylic acid in peach fruit. J Proteome Res. 2007;6:1677–88.
Jiang F, Zheng X, Chen J. Microarray analysis of gene expression profile induced by the biocontrol yeast Cryptococcus laurentii in cherry tomato fruit. Gene. 2009;430:12–6.
Hershkovitz V, Ben-Dayan C, Raphael G, Pasmanik-Chor M, Liu J, Belausov E, et al. Global changes in gene expression of grapefruit peel tissue in response to the yeast biocontrol agent Metschnikowia fructicola. Mol Plant Pathol. 2012;13:338–49.
Hershkovitz V, Sela N, Taha-Salaime L, Liu J, Rafael G, Kessler C, et al. De-novo assembly and characterization of the transcriptome of Metschnikowia fructicola reveals differences in gene expression following interaction with Penicillium digitatum and grapefruit peel. BMC Genomics. 2013;14:168.
Liu P, Cheng YJ, Yang M, Liu YJ, Chen K, Long CA, et al. Mechanisms of action for 2-phenylethanol isolated from Kloeckera apiculata in control of Penicillium molds of citrus fruits. BMC Microbiol. 2014;14:242.
Mizutani M, Ohta D. Diversification of P450 genes during land plant evolution. Annu Rev Plant Biol. 2010;61:291–315.
Guttikonda SK, Trupti J, Bisht NC, Chen H, An YQC, Pandey S, et al. Whole genome co-expression analysis of soybean cytochrome P450 genes identifies nodulationspecific P450 monooxygenases. BMC Plant Biol. 2010;10:243.
Kushiro T, Okamoto M, Nakabayashi K, Yamagishi K, Kitamura S, Asami T, et al. The Arabidopsis cytochrome P450 CYP707A encodes ABA 8′-hydroxylases: key enzymes in ABA catabolism. EMBO J. 2004;23:1647–56.
Van Loon LC, Rep M, Pieterse CMJ. Significance of inducible defense-related proteins in infected plants. Annu Rev Phytopathol. 2006;44:135–62.
Fu XZ, Chen CW, Wang Y, Liu JH, Moriguchi T. Ectopic expression of MdSPDS1 in sweet orange (Citrus sinensis Osbeck) reduces canker susceptibility: involvement of H2O2 production and transcriptional alteration. BMC Plant Biol. 2011;11:55.
Feng C, Chen M, Xu CJ, Bai L, Yin XR, Li X, et al. Transcriptomic analysis of Chinese bayberry (Myrica rubra) fruit development and ripening using RNA-Seq. BMC Genomics. 2012;13:19.
Fu XZ, Gong XQ, Zhang YX, Wang Y, Liu JH. Different transcriptional response to Xanthomonas citri subsp. citri between kumquat and sweet orange with contrasting canker tolerance. PLoS One. 2012;7:e41790.
Zamboni A, Zanin L, Tomasi N, Pezzotti M, Pinton R, Varanini Z, et al. Genome-wide microarray analysis of tomato roots showed defined responses to iron deficiency. BMC Genomics. 2012;13:101.
Ballester AR, Lafuente MT, Forment J, Gadea J, De Vos RCH, Bovy AG, et al. Transcriptomic profiling of citrus fruit peel tissues reveals fundamental effects of phenylpropanoids and ethylene on induced resistance. Mol Plant Pathol. 2011;12:879–97.
Tiwari BS, Belenghi B, Levine A. Oxidative stress increased respiration and generation of reative oxygen species, resulting in ATP depletion, opening of mitochondiral permeability transition, and programmed cell death. Plant Physiol. 2002;128:1271–81.
Yang K, Liu Y, Zhang S. Activation of a mitogen-activated protein kinase pathway is involved in disease resistance in tobacco. Proc Natl Acad Sci U S A. 2001;98:741–6.
Pieterse CMJ, Leon-Reyes A, Van der Ent S, Van Wees SCM. Networking by small-molecule hormones in plant immunity. Nat Chem Biol. 2009;5:308–16.
Pieterse CMJ, Van der Does D, Zamioudis C, Leon-Reyes A, Van Wees SCM. Hormonal modulation of plant immunity. Ann Rev Cell Dev Biol. 2012;28:489–521.
Loreti E, Povero G, Novi G, Solfanelli C, Alpi A, Perata P. Gibberellins, jasmonate and abscisic acid modulate the sucrose-induced expression of anthocyanin biosynthetic genes in Arabidopsis. New Phytol. 2008;179:1004–16.
Song S, Qi T, Huang H, Ren Q, Wu D, Chang C, et al. The Jasmonate-ZIM domain proteins interact with the R2R3-MYB transcription factors MYB21 and MYB24 to affect Jasmonate-regulated stamen development in Arabidopsis. Plant Cell. 2011;23:1000–13.
An D, Yang J, Zhang P. Transcriptome profiling of low temperature-treated cassava apical shoots showed dynamic responses of tropical plant to cold stress. BMC Genomics. 2012;13:64.
Serrano M, Martinez-Madrid MC, Romojaro F. Ethylene biosynthesis and polyamine and ABA levels in cut carnations treated with aminotriazole. J Am Soc Hort Sci. 1999;124:81–5.
Tieman DM, Loucas HM, Kim JY, Clark DG, Klee HJ. Tomato phenylacetaldehyde reductases catalyze the last step in the synthesis of the aroma volatile 2-phenylethanol. Phytochemistry. 2007;68:2660–9.
Walter S, Nicholson P, Doohan FM. Action and reaction of host and pathogen during Fusarium head blight disease. New Phytol. 2010;185:54–66.
Liu YZ, Liu Q, Tao NG, Deng XX. Efficient isolation of RNA from fruit peel and pulp of ripening navel orange (Citrus sinensis Osbeck). J Huazhong Agr Univ. 2006;25:300–4.
Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, et al. Minimum information about a microarray experiment (MIAME)–toward standards for microarray data. Nat Genet. 2001;29:365–71.
Gallego-Giraldo L, Jikumaru Y, Kamiya Y, Tang Y, Dixon RA. Selective lignin downregulation leads to constitutive defense response expression in alfalfa (Medicago sativa L.). New Phytol. 2011;190:627–39.
Yan JW, Yuan FR, Long GY, Qin L, Deng ZN. Selection of reference genes for quantitative real-time RT-PCR analysis in citrus. Mol Biol Rep. 2012;39:1831–8.
Syros T, Yupsanis T, Zafiriadis H, Economou A. Activity and isoforms of peroxidases, lignin and anatomy, during adventitious rooting in cuttings of Ebenus cretica L. J Plant Physiol. 2004;161:69–77.
Liu JH, Moriguchi T. Changes in free polyamine titers and expression of polyamine biosynthetic genes during growth of peach in vitro callus. Plant Cell Rep. 2007;26:125–31.
Shi J, Fu XZ, Peng T, Huang XS, Fan QJ, Liu JH. Spermine pretreatment confers dehydration tolerance of citrus in vitro plants via modulation of antioxidative capacity and stomatal response. Tree Physiol. 2010;30:914–22.
The research was financially supported by the National Basic Research Program of China (973 program; no: 2013CB127100), the National Natural Science Foundation of China (grant nos. 31171773, 30972062 and 31401831) and the Modern Agriculture (Citrus) Technology System (CARS-27).
The authors declare that they have no competing interests.
CL and PL designed the research; PL, CK, GL and XY performed research and analysed the data; PL and XY wrote the paper. All authors read and approved the final manuscript.
List of differentially expressed genes and GO and KEGG analyses between K. apiculata, the extract and phenylethanol treatment. The table contained information of the differential expressed genes with expression difference, and genes differential expressed at 0.05 significance level. (XLS 1869 kb)
List and MapMan analysis of differentially expressed genes in citrus in response to K. apiculata, the extract and phenylethanol, respectively. (XLS 264 kb)
List of defence-related differentially expressed genes in citrus under K. apiculata treatment. (XLS 115 kb)
The common or in special up-regulated or down-regulated expressed genes in citrus between K. apiculata, PEA and the extract treatment. (XLS 301 kb)
About this article
Cite this article
Liu, P., Chen, K., Li, G. et al. Comparative transcriptional profiling of orange fruit in response to the biocontrol yeast Kloeckera apiculata and its active compounds. BMC Genomics 17, 17 (2016) doi:10.1186/s12864-015-2333-3
- Biological control
- Kloeckera apiculata
- Plant defence