- Research article
- Open Access
Transcriptional profiling analysis of Penicillium digitatum, the causal agent of citrus green mold, unravels an inhibited ergosterol biosynthesis pathway in response to citral
© The Author(s). 2016
- Received: 15 February 2016
- Accepted: 18 July 2016
- Published: 11 August 2016
Green mold caused by Penicillium digitatum is the most damaging postharvest diseases of citrus fruit. Previously, we have observed that citral dose-dependently inhibited the mycelial growth of P. digitatum, with the minimum inhibitory concentration (MIC) of 1.78 mg/mL, but the underlying molecular mechanism is barely understood.
In this study, the transcriptional profiling of the control and 1/2MIC-citral treated P. digitatum mycelia after 30 min of exposure were analyzed by RNA-Seq. A total of 6355 genes, including 2322 up-regulated and 4033 down-regulated genes, were found to be responsive to citral. These genes were mapped to 155 KEGG pathways, mainly concerning mRNA surveillance, RNA polymerase, RNA transport, aminoacyl-tRNA biosynthesis, ABC transporter, glycolysis/gluconeogenesis, citrate cycle, oxidative phosphorylation, sulfur metabolism, nitrogen metabolism, inositol phosphate metabolism, fatty acid biosynthesis, unsaturated fatty acids biosynthesis, fatty acid metabolism, and steroid biosynthesis. Particularly, citral exposure affected the expression levels of five ergosterol biosynthetic genes (e.g. ERG7, ERG11, ERG6, ERG3 and ERG5), which corresponds well with the GC-MS results, the reduction in ergosterol content, and accumulation of massive lanosterol. In addition, ERG11, the gene responsible for lanosterol 14α-demethylase, was observed to be the key down-regulated gene in response to citral.
Our present finding suggests that citral could exhibit its antifungal activity against P. digitatum by the down-regulation of ergosterol biosynthesis.
- Penicillium digitatum
- Ergosterol biosynthesis
Significant losses can occur after harvest during storage and marketing of citrus fruit primarily due to green mold, caused by Penicillium digitatum, and secondarily by blue mold and sour rot caused by P. italicum and Geotrichum citri-aurantii, respectively . Among them, P. digitatum is the most economically devastating pathogen causing about 90 % of the total loss of postharvest citrus fruit . Currently, the control of postharvest pathogens depends mainly on synthetic chemical application, but extensive use of synthetic chemical has caused the emergence of fungicide-resistant populations, therefore seeking natural and effective microbiocidal agents as alternatives has been immediate areas of research focus [3, 4].
Citral, one of the volatile constituents in plant essential oils, has been demonstrated to have strong antifungal activity against P. digitatum, P. italicum, and G. citri-aurantii [4–6]. Fumigation of oranges with citral (20, 60 or 150 mL/L in absorbent pads) in a closed system, following application of conidia (20 μL of 106 conidia/mL) to puncture wounds, delayed the onset of sour rot at room temperature by 7–10 days and at 5 °C, by 13–30 days . Previously, we have observed that citral inhibited the myclial growths of P. italicum and G. citri-aurantii in a dose dependent manner [7, 8], and the application of wax enriched with citral significantly decreased the incidence of green mold after 6 days of storage at 25 ± 2 °C . Therefore, it might be an alternative as fungicide in controlling postharvest diseases in citrus fruit.
The antifungal mechanism of volatile compounds has been attributed to its capacity to disturb the cellular membrane, interfere with the cellular metabolism, react with active sites of enzymes, or act as H+ carriers [10, 11]. In our previous studies, citral was found to destroy the membrane permeability and integrity of P. italicum and G. citri-aurantii by causing significant losses in total lipids or ergosterol contents [7, 8]. In addition, citral at a minimum inhibitory concentration (MIC, 1.78 mg/mL) evidently altered the mitochondrial morphology and inhibited the citrate cycle (TCA cycle) of P. digitatum . However, information about the inhibitory mechanism of citral on P. digitatum at molecular level is rather limited, and thus, requires further studies.
Recently, numerous reports regarding the global gene expression in response to essential oils or their volatile components in fungal have been conducted. Parveen et al.  found that the cell wall- and membrane-related genes were major targets of α-terpinene against Saccharomyces cerevisiae, showing that genes associated with protein analysis, carbohydrate metabolism, and transcription were repressed, whereas genes involved in ergosterol biosynthesis, phospholipid biosynthesis, cell wall organization and detoxification were activated. Rao et al.  reported that carvacrol can stimulate the expression levels of genes involved in alternate metabolic and energy pathways, stress response, autophagy, drug efflux, but repress genes mediating ribosome biogenesis and RNA metabolism. A microarray analysis revealed that 1 % (v/v) p-anisaldehyde can induce the expression levels of genes related to sulphur assimilation, aromatic aldehydes metabolism, and secondary metabolism in S. cerevisiae, suggesting that the normal metabolism of aromatic aldehydes was influenced . In another report, 90 min of exposure to 128 μg/mL thymol resulted in a reduction in intracellular thiamine concentration in yeast by impairing thiamin metabolism. In contrast, many genes involving in maintaining cell membrane and cell wall integrity, sulfur assimilation, methionine biosynthesis, and production of AdoMET were drastically induced . These findings provide some novel clues for the understanding of the antifungal mechanism of plant essential oils.
RNA-Seq is an increasingly attractive method for whole-genome expression studies in many biological systems. As compared with traditional microarrays, despite their instrumental role in profiling the global expression of genes, this technique extends the possibilities of transcriptome studies to the analysis of previously unidentified genes and splice variants due to the advantages such as no need to provide a referenced genomic sequence, a fully coverage of three types of RNA, an unlimited dynamic range of quantification at reduced technical variability, the declining cost of sequencing, and so on [17, 18]. It has been widely used to unravel the pathogenetic mechanisms or diseases control in citrus fruit [19–21]. Therefore, our present study is to analyze the global gene expression profiles in P. digitatum mycelia with or without citral treatment by RNA-Seq, in an effort to explore the underlying molecular mechanism and to find some key metabolic pathways or genes involved in.
P. digitatum was isolated from infected citrus fruit and preserved on potato dextrose agar (PDA) at 25 ± 2 °C. Two hundred micro liter fungal suspensions (5 × 105 cfu/mL) were added to 40 mL potato dextrose broth (PDB) and incubated in a moist chamber at 28 ± 2 °C for 72 h. The mycelia were vacuum-filtered and weighted at a 6 h interval to make a growth curve.
Based on the result of growth curve, 1 g wet mycelia in logarithmic metaphase (48 h of culture) were added to 20 mL PBS (pH 6.8) and incubated with 1/2MIC (half of minimum inhibitory concentration; 0.89 mg/mL) or MIC of citral for 0, 30, 60, and 120 min. Samples without citral were severed as a control. The resulting mycelia were vacuum-filtered, weighted, and recorded to select optimal concentration and time for the next analysis. All collected P. digitatum mycelia were grinded to powders in liquid nitrogen and stored at −80 °C until further use.
RNA isolation, integrity examination and RNA-Seq library preparation
Total RNAs from control, 1/2 MIC or MIC citral-treated samples after 30 min of exposure were extracted with TRIzol regent (Invitrogen, USA) according to the manufacturer’s instruction and treated with RNase-free DNase I (Takara Biotechnology, China). RNA integrity was determined by a 2100 bioanalyzer (Agilent, USA). Poly (A) mRNA from control and 1/2MIC citral-treated samples, designated as CK30 and T30, respectively, was isolated with oligo-dT beads and then treated with the fragmentation buffer. The cleaved RNA fragments were then transcribed into first-strand cDNA using reverse transcriptase and random hexamer primers, followed by second strand cDNA synthesis using DNA polymerase I and RNaseH. The double stranded cDNA was further subjected to end-repair using T4 DNA polymerase, Klenow fragment, and T4 Polynucleotide kinase, followed by a single A base addition using Klenow 3′ to 5′ exo-polymerase. It was then ligated with adapter or index adapter using T4 quick DNA ligase. Adaptor ligated fragments were selected according to the size and the desired range of cDNA fragments were excised from the gel. PCR was performed to selectively enrich and amplify the fragments. Finally, after validating on an Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-Time PCR System, the cDNA library was sequenced on a flow cell using Illumina HiSeq2000 ™ .
Assembly and functional annotation
Transcriptome de novo assembly of the short reads was carried out using Trinty software . After filtration of the low quality reads, clean reads were randomly clipped into 25-bp k-mers and assembled using de Bruijin graph and Trinty software. Furthermore, the reads producing large fragment without N is named contig and the result sequences of trinity is called unigenes. Then unigenes were clustered using TGI Clustering tools to obtain the finally unigenes . Finally, blastx alignment (e value < 0.00001) between unigenes and protein databases (non-redundant protein (NR), Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG) and the Clusters of Orthologous Groups of proteins (COG)) were performed, and the best aligning results were selected as a priority order of NR, Swiss-Prot, KEGG and COG to determine the sequence orientation. When a unigene was not aligned to the above databases, ESTScan software  was used to predict its coding regions and to determine its sequence orientation. Gene ontology (GO) terms annotations of unigenes were performed with the Blast2GO program according to the NR annotation. The COG and KEGG pathway annotations were analyzed using Blastall software against the cluster of orthologous groups’ database and the kyoto encyclopedia of genes and genomes database.
Differential expression of unigenes
Where, N1 is the total fragments number of the CK30, and N2 is total fragments number of T30; gene A holds x fragments in CK30 and i fragments in T30. The larger difference of the expression level between the CK30 and T30 is calculated with FDR ≤ 0.001 and log2 Ratio ≥ 1.
GO annotation and pathway analysis of differentially expressed genes (DEGs)
Where, in GO analysis, N is the number of all genes with GO annotation; n is the number of DEGs in N; M is the number of all genes that are annotated to the certain GO terms; m is the number of DEGs in M. The calculated p value goes through Bonferroni Correction, taking corrected-P value ≤ 0.05 as a threshold. In KEGG analysis, N is the number of all genes with KEGG annotation, n is the number of DEGs in N, M is the number of all genes annotated to specific pathways, and m is number of DEGs in M.
Real-time Fluorescence Quantitative PCR (RTFQ-PCR) analysis
Primer pair sequences designed for validation of differentially expressed genes in CK30 and T30 treatment P. digitatum using Real-time Fluorescence Quantitative PCR (RTFQ-PCR)
Primer sequence (5′-3′, forward/reverse)
Sterols compositions by gas chromatography-mass spectrometry (GC-MS)
The sterols compositions of P. digitatum cells exposure to citral at a concentration of 1/2MIC for 0, 30, 60, and 120 min were determined by GC-MS method [28, 29] with some modifications. The 2-day-old mycelia from 50 mL PDB were collected and centrifuged at 4000 g for 15 min. Then the samples were dried with a vacuum freeze drier for 4 h. About 0.5 g of dry mycelia were homogenized with liquid nitrogen and were suspended in 4 mL of freshly prepared 25 % (w/v) NaOH and 8 mL of absolute ethanol, vortexed for 30 sec and sonicated for 10 min, and then saponified at 80 °C for 1 h. The mixtures were extracted with 2 mL petroleum ether for three times and washed by saturated NaCl solution twice. The organic upper layer was transferred into a new 10.0 mL plastic microcentrifuge safe-lock tube containing 40 ± 2 mg anhydrous sodium sulphate. The samples were then vacuum concentrated. Each residue was dissolved in 950 μL of methyl tert-butyl ether and 50 μL of silylation reagent mixture N-methyl-N-trimethylsilyltrifluoroacetamide (MSTFA) and N-trimethylsilylimidazole (TSIM) (9:1) was added. The samples were gently shaken and stored for completion of the silylation reaction at room temperature for at least 30 min, before being subjected to GC–MS analysis. The samples without citral treatment were used as a control.
The analytical GC was carried out on a Shimadzu QP2010 plus gas chromatograph (Shimadzu, Kyoto, Japan) equipped with flame ionization detector (FID). A non-polar cross-linked fused-silica capillary column, HP-5MS (30 m × 0.25 mm × 0.25 μm; Agilent, Santa Clara, CA), was used. The oven temperature was held at 50 °C for 1 min, programmed at a rate of 50 °C/min to 260 °C, then increased to 300 °C at a rate of 4 °C/min where it remained for 10 min. The carrier gas was helium (1.3 mL/min). The injector temperature was 250 °C, detector temperature 310 °C and the volume injected was 2 μL. MS analysis was carried out on the same chromatograph equipped with a Shimadzu QP 2010 GC/MS system, ionisation voltage 70 eV, ion source temperature 200 °C, mass range m/z 50–700, scanning interval 0.5 s and scanning speed 1000 amu/s. The sterol TMS ethers were identified by comparison with commercial references, the NIST™ database or data from literature.
Determination of total ergosterol and lanosterol contents
Total ergosterol and lanosterol contents of P. digitatum cells exposure to citral at a concentration of 1/2MIC for 0, 30, 60, and 120 min were determined by high-performance liquid chromatography (HPLC) as reported by our previous study , whereas the detected wavelength for lanosterol was set at 210 nm. The samples without citral treatment were used as a control.
All data were expressed as the mean ± SD by measuring three independent replicates. Analysis of variance using one-way ANOVA followed by Duncan’s test was performed to test the significance of differences between means obtained among the treatments at the 5 % level of significance using SPSS statistical software package release 16.0 (SPSS Inc., Chicago, IL, USA).
The growth of P. digitatum
To determine the effects of citral on the growth of P. digitatum, the wet weight of P. digitatum mycelia exposed to citral was determined. As shown in Fig. 1b, citral markedly inhibited the mycelial growth of P. digitatum during the whole period, whereas that of the control group remained stable. After 30 min of exposure, the wet weight of 1/2MIC and MIC mycelia were 0.83 ± 0.03 g and 0.71 ± 0.06 g, respectively, which were significantly lower than that of the control (0.99 ± 0.06 g, P < 0.05). This inhibition was more obvious with the prolonging of the culture time. Considering that a RNA integrated number (RIN) value higher than 6.5 was the prerequisite for the RNA-Seq analysis, we further determined the RIN values in citral-treated or control samples at 30 min incubation. Results showed that the RIN values in 1/2 MIC or MIC citral -treated and control groups were 6.7, 2.6, and 7.9, respectively (data not show), suggesting that MIC citral -treated samples were not suitable for further analysis. Therefore, only samples treated with 1/2MIC citral (T30) or without citral treatment (CK30) were chosen for RNA-Seq analysis.
Overall transcriptome profiles P. digitatum
Summary data of reads in untreated strains (CK30) and citral-treated strains (T30) of P. digitatum mycelia transcriptomes
CK30 & T30
Total clean reads
Total clean nucleotides (nt)
Total number of contigs
Mean length of contigs (nt)
The number of unigenes
Mean length of unigenes (nt)
The number of all-unigenes
Mean length of unigenes
Unigenes annotation against NR
Unigenes annotation against NT
Unigenes annotation against Swiss-Prot
Unigenes annotation against KEGG
Unigenes annotation against COG
Unigenes annotation against GO
A total of 16,898 all-unigenes were mapped to the referenced canonical pathways in KEGG, 8714 of them were annotated to 168 KEGG pathways. Among of which, 2550, 1087, and 898 unigenes were distributed to metabolic pathway, biosynthesis of secondary metabolites, and microbial metabolism in diverse environments, respectively. Few unigenes correlated with carotenoid biosynthesis (1 unigene), phosphotransferase system (PTS, 1 unigene), and sesquiterpenoid and triterpenoid biosynthesis (2 unigenes) were also found.
Differences in overall gene expressions and pathway between CK30 and T30
The differences in gene expressions between CK30 and T30 were further analyzed (Additional file 3: Table S3). The variations of gene length and total reads distributions in gene expression levels were eliminated by FPKM calculation. The different gene expressions were calculated by false discovery rate (FDR ≤ 0.001) and the absolute value of log2Ration ≥ 1 as a threshold. A total of 6355 differential expression genes were detected between CK30 and T30, including 2322 up-regulated and 4033 down-regulated genes (Table 2). The distribution of down-regulated genes was nearly twice as that of up-regulated genes, which might be caused by the inhibition of P. digitatum mycelial growth.
Enrichment pathway analysis of DEGs in P. digitatum
DEGs genes with pathway annotation (3232)
All genes with pathway annotation (8714)
64 (1.98 %)
134 (1.54 %)
Ribosome biogenesis in eukaryotes
78 (2.41 %)
171 (1.96 %)
Alanine, aspartate and glutamate metabolism
40 (1.24 %)
83 (0.95 %)
Starch and sucrose metabolism
175 (5.41 %)
426 (4.89 %)
91 (2.82 %)
213 (2.44 %)
Fatty acid biosynthesis
28 (0.87 %)
62 (0.71 %)
21 (0.65 %)
45 (0.52 %)
MAPK signaling pathway
88 (2.72 %)
221 (2.54 %)
Amino sugar and nucleotide sugar metabolism
98 (3.03 %)
248 (2.85 %)
Fructose and mannose metabolism
90 (2.78 %)
230 (2.64 %)
Phosphatidylinositol signaling system
22 (0.68 %)
53 (0.61 %)
16 (0.5 %)
38 (0.44 %)
38 (1.18 %)
98 (1.12 %)
80 (2.48 %)
210 (2.41 %)
Inositol phosphate metabolism
30 (0.93 %)
79 (0.91 %)
121 (3.74 %)
326 (3.74 %)
Citrate cycle (TCA cycle)
21 (0.65 %)
56 (0.64 %)
Linoleic acid metabolism
36 (1.11 %)
97 (1.11 %)
64 (1.98 %)
173 (1.99 %)
mRNA surveillance pathway
55 (1.7 %)
150 (1.72 %)
29 (0.9 %)
80 (0.92 %)
Pentose phosphate pathway
17 (0.53 %)
48 (0.55 %)
94 (2.91 %)
267 (3.06 %)
86 (2.66 %)
247 (2.83 %)
24 (0.74 %)
73 (0.84 %)
32 (0.99 %)
96 (1.1 %)
42 (1.3 %)
127 (1.46 %)
20 (0.62 %)
64 (0.73 %)
47 (1.45 %)
144 (1.65 %)
Biosynthesis of unsaturated fatty acids
14 (0.43 %)
53 (0.61 %)
Fatty acid metabolism
43 (1.33 %)
147 (1.69 %)
29 (0.9 %)
109 (1.25 %)
1 (0.03 %)
11 (0.13 %)
33 (1.02 %)
224 (2.57 %)
Genes related to stress response
To survive from drastic environment changes, eukaryotes must rapidly adjust their gene expression programs to develop a series of protective mechanisms [30, 31]. In this study, considerable alterations in gene expression levels related to cell rescue, defense, and detoxification were influenced by citral. Genes associated with the over-represented pathway, such as ribosome biogenesis in eukaryotes, mRNA surveillance pathway, RNA polymerase, RNA transport and aminoacyl-tRNA biosynthesis, were partially repressed, indicating a lesion in the translational activity of cells. This result, to a large extent, accounts for the low RIN value (2.6) in P. digitatum cells treated with MIC of citral for 30 min. In addition, one gene responsible for mitochondrial ribosomal protein (MRPS5) was found to be down-regulated, indicating a blocking in mitochondrial translation and the loss of mitochondrial functions . Correspondingly, a repression in energy-related pathways including TCA cycle, glycolysis/gluconeogenesis, oxidative phosphorylation, sulfur metabolism, and nitrogen metabolism, was also observed. Most genes in these pathways, such as SDH1, frdA, CS, IDH1, HK, PFK, PK, fbcH, ATPeF1G, ATPeV0C, cysQ, cysJ, aspA, nasB, and glnA, were down-regulated. Previously, we have demonstrated that citral treatment could evidently alter the mitochondrial morphology, reduce the ATP content and inhibit the TCA cycle of P. digitatum . These results supported the hypothesis that the down-regulation of energy-demanding processes is a very important approach to help the cells escaping from energy crisis under extreme environmental conditions . Interestingly, significant changes in genes involved in the environmental stress response were also observed, which was consistent with a previous report describing the global gene expressions in S. cerevisiae cells treated with thymol .
Genes involved in multidrug resistance
The typical gene subfamilies involved in multidrug resistance in fungal, such as multi-drug resistance (MRD, ABCB subfamily), multi-drug resistance-associated protein (MRP, ABCC subfamily), and pleiotropic drug resistance (PDR, ABCG subfamily), were affected by citral. A slight increase in the MRD gene PMR1 expression was observed in samples exposed to citral. This gene was an important determinant of resistance to demethylation inhibitor-resistant strains, which is normally suggested to protect the cells against triflumizole . Genes belonged to plasma membrane ATP-binding cassette (ABC) transporters, such as PDR5 and SNQ2, were also up-regulated. These two genes were involved in a robust induction of drug efflux mechanisms, and their expression levels in S. cerevisiae cells have been demonstrated to be strongly induced by thymol . Interestingly, the expression levels of five MRP genes including MRP1, MRP4, MRP5, MRP7, and MRP13, were down-regulated. This type of genes generally transports xenobiotic compounds or toxic metabolites that conjugated to glutathione, glucuronide, or sulfur, and the over-expression of MRP1 gene can enhance the multidrug resistance in eukaryotic cells . Therefore, the decrease in expression levels of MRP genes suggested that the multidrug resistance in P. digitatum cell might be impaired by citral.
Genes associated with cell integrity
Several targets including cell wall, cell membrane, mitochondrion, and genetic material, have been proposed to account for the antifungal activity of essential oils or their volatile compositions [3, 12–16, 36]. The cell wall is an extracellular layer outside the cell membrane which protects the cell against mechanical damage, osmotic strength and determines the cell shape. We found that several pathways associated with cell wall biogenesis, such as amino sugar and nucleotide sugar metabolism, starch and sucrose metabolism, were drastically influenced by citral. CHI1, Fks1, two important cell wall-related genes encoding for chitin synthase and 1,3-β-glucan synthase, respectively, were repressed. This finding underlined that citral can block the formation of cell wall architecture and thereby improve the sensitivity of cells to citral. In contrast, Parveen et al.  reported that theses two genes were dramatically induced in S. cerevisiae cell exposed to α-terpinene, which in turn activated cell wall compensatory mechanism to overcome the terpinene toxicity. This discrepancy might be, to a large extent, caused by different strains, or by different antifungal mechanisms, as suggested by other reports [36, 37].
The lipophilic nature of terpenoid enables it preferentially enter into the lipid membrane, which results in the increased membrane fluidity and eventually leads to the increase of membrane permeability [38, 39]. The expression levels of some genes involved in cell membrane-related pathways, such as inositol phosphate metabolism, fatty acid biosynthesis, biosynthesis of unsaturated fatty acids, fatty acid metabolism, and steroid biosynthesis, were changed upon citral treatment. Among them, INO1, PI4K, CDIPT, accC, FAS1, FAS2, FAB1, ERG7, ERG11, ERG6, ERG3, ERG5 were down-regulated, suggesting that the integrity of cell membrane might be destroyed. Oppositely, expression of FAD6 was increased by 12.4 folds. These results were consistent with our previous study that citral could induce the disruption of membrane integrity by causing a decrease in total lipid contents of pathogen cells .
Ergosterol is one of the principal sterol components in fungal membrane. Generally, a decrease in ergosterol contents results in osmotic disturbances, disruption of cell growth and proliferation [7, 39–41]. Some chemical fungicides widely used in controlling the green mold of citrus fruit, such as imazalil, prochloraz, and triflumizole, exhibited their antifungal mode by blocking the ergosterol biosynthesis, which can give rise to the disruption of cell structure and function, even to the death of cell [42–44]. Notably, the expression levels of five genes in ergosterol biosynthetic pathway encoding for lanosterol synthase (ERG7), lanosterol 14α-demethylase (ERG11), C-24 sterol methyl-transferase (ERG6), C-5 sterol desaturase (ERG3), and C-22 sterol desaturase (ERG5), were down-regulated by 2.7, 2.1, 2.0, 2.1 and 2.0 folds, respectively.
Effect of citral on gene expression levels in ergosterol biosynthesis
Analysis of sterols composition by GC-MS
The sterols composition of P. digitatum by GC-MS
Major ions [m/z]
1/2MIC citral-treated group
Ergosterol and lanosterol contents
The ergosterol and lanosterol contents in 1/2MIC citral-treated P. digitatum cells correspond well with the transcriptional regulation of ERG7, ERG11, ERG6, ERG3, and ERG5 genes. After the initial addition of citral, the gene expressions of ERG11, ERG6, ERG3, and ERG5 were partially inhibited and these led to the decrease in total ergosterol content before 30 min of exposure. Correspondingly, the absence of ergosterol might induce the gene expression of ERG7 at this time, as convinced by the up-regulation of this gene and the accumulation of lanosterol content. The gene expressions of ERG7 in both groups have the same changing tendency after 60 min of exposure, but the expression levels in 1/2MIC citral-treated samples were significantly higher (P < 0.05) than those in control group. This up regulation, together with the reduction of ergosterol content, might stimulate the gene expressions of ERG6, ERG3 and ERG5 to provide more ergosterol to meet the demand of cells, which is thought to be a general response to decreased ergosterol levels . This regulatory manner is similar to a previous finding that described in S. cerevisiae cells exposed to α-terpinene . Nevertheless, the gene expressions of ERG6 and ERG3 in 1/2MIC group were still lower than the control groups, which were supposed to contribute a lot to the inhibition of ergosterol biosynthesis, as illustrated by other documents [41, 55–57]. The gene expression of ERG5 in 1/2MIC group was lower than that in the CK group at 120 min of treatment, which might be caused by the reduction in its precursor (ergosta-5,7,24(28)-trienol). Interestingly, the expression levels of ERG11 continuously declined in 1/2MIC citral-treated groups and were always lower than those of CK group. Considering that the down-regulation of ERG11 can block the ergosterol biosynthesis and lead to the decrease of ergosterol content to inhibit the growth of fungi [43, 47], our present finding implies that ERG11 is a critical target of citral against P. digitatum.
In conclusion, our present research revealed that the global gene expressions involved in many important pathways in P. digitatum hyphae, including ABC transport, steroid biosynthesis, amino sugar and nucleotide sugar metabolism, starch and sucrose metabolism, TCA cycle, oxidative phosphorylation, RNA degradation, and ribosome, can be greatly affected by citral exposure. RTFQ-PCR verified that citral can affect the expression level of five genes (ERG7, ERG11, ERG6, ERG3 and ERG5) in ergosterol biosynthetic pathway of P. digitatum, showing that ERG11 is a critical targeted gene. In addition, the ergosterol and lanosterol contents in citral-treated samples correspond well with the coordinated regulation of these five genes and GC-MS results. These results suggest that citral could exhibit its antifungal activity against the mycelial growth of P. digitatum by disrupting ergosterol biosynthesis.
MIC, minimum inhibitory concentration; PDA, potato dextrose agar; PDB, potato dextrose broth; DEGs, differentially expressed genes; FPKM, fragments per Kb million fragments; FDR, false discovery rate; RTFQ-PCR, real-time fluorescence quantitative PCR; RIN, RNA integrated number; CDS, coding sequence; NR, non-redundant protein; NT, non-redundant nucleotide sequence; KEGG, Kyoto Encyclopedia of Genes and Genomes; COG, clusters of orthologous groups of proteins; GO, gene ontology; NCBI, National Center for Biotechnology Information; GC-MS, gas chromatography-mass spectrometry; HPLC, high-performance liquid chromatography; MSTFA, N-methyl-N-trimethylsilyltrifluoroacetamide; TSIM, N-trimethylsilylimidazole; MRD, multi-drug resistance; MRP, Multi-drug resistance-associated protein; PDR, Pleiotropic drug resistance; 1/2MIC, Half of minimum inhibitory concentration
This study was supported by the National Natural Science Foundation of China (Nos.31572172 and 31271964) and Research Foundation of Education Bureau of Hunan Province (Nos.15A 181 and 12B126).
Availability of data and materials
The datasets generated during and/or analyzed during the current study are not publicly available due to the follow-up study about these datasets is under way but are available from the corresponding author on reasonable request.
NT designed research; QO performed research; NT and QO analyzed data; and QO, NT and GJ wrote the paper. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Smilanick JL, Mansour MF, Gabler FM, Sorenson D. Control of citrus postharvest green mold and sour rot by potassium sorbate combined with heat and fungicides. Postharvest Biol Tec. 2008;47:226–38.View ArticleGoogle Scholar
- Macarisin D, Cohen L, Eick A, Rafael G, Belausov E, Wisniewski M, Droby S. Penicillium digitatum suppresses production of hydrogen peroxide in host tissue during infection of citrus fruit. Phytopathology. 2007;97:1491–500.View ArticlePubMedGoogle Scholar
- Shao X, Cheng S, Wang H, Yu D, Mungai C. The possible mechanism of antifungal action of tea tree oil on Botrytis cinerea. J Appl Microbiol. 2013;114:1642–9.View ArticlePubMedGoogle Scholar
- Wuryatmo E, Able AJ, Ford CM, Scott ES. Effect of volatile citral on the development of blue mould, green mould and sour rot on navel orange. Australas Plant Path. 2014;43:403–11.View ArticleGoogle Scholar
- Droby S, Eick A, Macarisin D, Cohen L, Rafael G, Stange R, McColumb G, Dudaic N, Nasserd A, Wisniewski M, et al. Role of citrus volatiles in host recognition, germination and growth of Penicillium digitatum and Penicillium italicum. Postharvest Biol Tec. 2008;49:386–96.View ArticleGoogle Scholar
- Wuryatmo E, Klieber A, Scott ES. Inhibition of citrus postharvest pathogens by vapor of citral and related compounds in culture. J Agr Food Chem. 2003;51:2637–40.View ArticleGoogle Scholar
- Tao NG, OuYang QL, Jia L. Citral inhibits mycelial growth of Penicillium italicum by a membrane damage mechanism. Food Control. 2014;41:116–21.View ArticleGoogle Scholar
- Zhou HE, Tao NG, Jia L. Antifungal activity of citral, octanal and α–terpineol against Geotrichum citri–aurantii. Food Control. 2014;37:277–83.View ArticleGoogle Scholar
- Fan F, Tao NG, Jia L, He XL. Use of citral incorporated in postharvest wax of citrus fruit as a botanical fungicide against Penicillium digitatum. Postharvest Biol Tec. 2014;90:52–5.View ArticleGoogle Scholar
- Bajpai VK, Sharma A, Baek KH. Antibacterial mode of action of Cudrania tricuspidata fruit essential oil, affecting membrane permeability and surface characteristics of food–borne pathogens. Food Control. 2013;32:582–90.View ArticleGoogle Scholar
- Ultee A, Bennik MHJ, Moezelaar R. The phenolic hydroxyl group of carvacrol is essential for action against the food–borne pathogen Bacillus cereus. Appl Environ Microb. 2002;68:1561–8.View ArticleGoogle Scholar
- Zheng SJ, Jing GX, Wang X, Ouyang QL, Jia L, Tao NG. Citral exerts its antifungal activity against Penicillium digitatum by affecting the mitochondrial morphology and function. Food Chem. 2015;158:76–81.View ArticleGoogle Scholar
- Parveen M, Hasan MK, Takahashi J, Murata Y, Kitagawa E, Kodama O, Iwahashi H. Response of Saccharomyces cerevisiae to a monoterpene: evaluation of antifungal potential by DNA microarray analysis. J Antimicrob Chemother. 2004;54:46–55.View ArticlePubMedGoogle Scholar
- Rao A, Zhang YQ, Muend S, Rao R. Mechanism of antifungal activity of terpenoid phenols resembles calcium stress and inhibition of the TOR pathway. Antimicrob Agents Chemother. 2010;54:5062–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Yu L, Guo N, Yang Y, Wu XP, Meng RZ, Fan JW, Wang XL, Liu JB, Deng XM. Microarray analysis of p–anisaldehyde–induced transcriptome of Saccharomyces cerevisiae. J Ind Microbiol Biot. 2010;37:313–22.View ArticleGoogle Scholar
- Bi X, Guo N, Jin J, Liu J, Feng H, Shi J, Xiang H, Wu X, Dong J, Hu H, et al. The global gene expression profile of the model fungus Saccharomyces cerevisiae induced by thymol. J Appl Microbiol. 2010;108:712–22.View ArticlePubMedGoogle Scholar
- Wang LK, Feng ZX, Wang X, Wang XW, Zhang XG. DEGseq: an R package for identifying differentially expressed genes from RNA–seq data. Bioinformatics. 2010;26:136–8.View ArticlePubMedGoogle Scholar
- Van Verk MC, Hickman R, Pieterse CM, Van Wees SC. RNA–Seq: revelation of the messengers. Trends Plant Sci. 2013;18:175–9.View ArticlePubMedGoogle Scholar
- Hershkovitz V, Sela N, Taha–Salaime L, Liu J, Rafael G, Kessler C, Aly R, Levy M, Wisniewski M, Droby S. 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Aritua V, Achor D, Gmitter FG, Albrigo G, Wang N. Transcriptional and microscopic analyses of citrus stem and root responses to Candidatus Liberibacter asiaticus infection. PLoS One. 2013;8.Google Scholar
- Liu P, Cheng YJ, Yang M, Liu YJ, Chen K, Long CA, Deng XX. Mechanisms of action for 2–phenylethanol isolated from Kloeckera apiculata in control of Penicillium molds of citrus fruits. BMC Microbiol. 2014;14:242.View ArticlePubMedPubMed CentralGoogle Scholar
- Liu B, Jiang GF, Zhang YF, Li JL, Li XJ, Yue JS, Chen F, Liu HQ, Li HJ, Zhu SP. Analysis of transcriptome differences between resistant and susceptible strains of the citrus red mite Panonychus citri (Acari: Tetranychidae). PLoS One. 2011;6:e28516.View ArticlePubMedPubMed CentralGoogle Scholar
- Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng QD, et al. Full–length transcriptome assembly from RNA–Seq data without a reference genome. Nat Biotechnol. 2011;29:644–52.View ArticlePubMedPubMed CentralGoogle Scholar
- Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, et al. TIGR Gene Indices clustering tools (TGICL), a software system for fast clustering of large EST datasets. Bioinformatics. 2003;19:651–2.View ArticlePubMedGoogle Scholar
- Iseli C, Jongeneel CV, Bucher P. ESTScan: a program for detecting, evaluating, and reconstructing potential coding regions in EST sequences. Proc Int Conf Intell Syst Mol Biol. 1999;99:138–48.Google Scholar
- Audic S, Claverie JM. The significance of digital gene expression profiles. Genome Res. 1997;7:986–95.PubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using Real–Time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25:402–8.View ArticlePubMedGoogle Scholar
- Müller C, Staudacher V, Krauss J, Giera M, Bracher F. A convenient cellular assay for the identification of the molecular target of ergosterol biosynthesis inhibitors and quantification of their effects on total ergosterol biosynthesis. Steroids. 2013;78:483–93.View ArticlePubMedGoogle Scholar
- Alcazar-Fuoli L, Mellado E, Garcia-Effron G, Lopez JF, Grimalt JO, Cuenca-Estrella JM, Rodriguez-Tudela JL. Ergosterol biosynthesis pathway in Aspergillus fumigatus. Steroids. 2008;73:339–47.View ArticlePubMedGoogle Scholar
- Legay G, Marouf E, Berger D, Neuhaus JM, Mauch–Mani B, Slaughter A. Identification of genes expressed during the compatible interaction of grapevine with Plasmopara viticola through suppression subtractive hybridization (SSH). Eur J Plant Pathol. 2011;129:281–301.View ArticleGoogle Scholar
- Li JW, Liu J, Zhang H, Xie CH. Identification and transcriptional profiling of differentially expressed genes associated with resistance to Pseudoperonospora cubensis in cucumber. Plant Cell Rep. 2011;30:345–57.View ArticlePubMedGoogle Scholar
- Zhang L, Zhang Y, Zhou YM, An S, Zhou YX, Cheng J. Response of gene expression in Saccharomyces cerevisiae to amphotericin B and nystatin measured by microarrays. J Antimicrob Chemother. 2002;49:905–15.View ArticlePubMedGoogle Scholar
- Gasch AP, Spellman PT, Kao CM, Carmel–Harel O, Eisen MB, Storz G, Botstein D, Brown PO. Genomic expression programs in the response of yeast cells to environmental changes. Mol Biol Cell. 2002;11:4241–57.View ArticleGoogle Scholar
- Nakaune R, Adachi K, Nawata O, Tomiyama M, Akutsu K, Hibi T. A novel ATP–binding cassette transporter involved in multidrug resistance in the phytopathogenic fungus Penicillium digitatum. Appl Environ Microb. 1998;64:3983–8.Google Scholar
- Paumi CM, Chuk M, Snider J, Stagljar I, Michaelis S. ABC transporters in Saccharomyces cerevisiae and their interactors: new technology advances the biology of the ABCC (MRP) subfamily. Microbiol Mol Biol R. 2009;73:577–93.View ArticleGoogle Scholar
- Bakkali F, Averbeck S, Averbeck D, Idaomar M. Biological effects of essential oils–a review. Food Chem Toxicol. 2008;46:446–75.View ArticlePubMedGoogle Scholar
- Kalemba D, Kunicka A. Antibacterial and antifungal properties of essential oils. Curr MedChem. 2003;10:813–29.Google Scholar
- Burt S. Essential oils: their antibacterial properties and potential applications in foods—a review. Int J Food Microbiol. 2004;94:223–53.View ArticlePubMedGoogle Scholar
- Khan A, Ahmad A, Akhtar F, Yousuf S, Xess I, Khan LA, Manzoor N. Ocimum sanctum essential oil and its active principles exert their antifungal activity by disrupting ergosterol biosynthesis and membrane integrity. Res Microbiol. 2010;161:816–23.View ArticlePubMedGoogle Scholar
- Onyewu C, Blankenship JR, Del Poeta M, Heitman J. Ergosterol biosynthesis inhibitors become fungicidal when combined with calcineurin inhibitors against Candida albicans, Candida glabrata, and Candida krusei. Antimicrob Agents Chemother. 2003;47:956–64.View ArticlePubMedPubMed CentralGoogle Scholar
- Sun XP, Wang JY, Feng D, Ma ZH, Li HY. PdCYP51B, a new putative sterol 14a–demethylase gene of Penicillium digitatum involved in resistance to imazalil and other fungicides inhibiting ergosterol synthesis. Appl Microbiol Biot. 2011;91:1107–19.View ArticleGoogle Scholar
- Hamamoto H, Hasegawa K, Nakaune R, Lee YJ, Makizumi Y, Akutsu K, Hibi T. Tandem repeat of a transcriptional enhancer upstream of the sterol 14α-demethylase gene (CYP51) in Penicillium digitatum. Appl Environ Microb. 2000;66:3421–6.View ArticleGoogle Scholar
- Liu J, Yuan YZ, Wu Z, Li N, Chen YL, Qin TT, Geng HH, Xiong L, Liu D. A novel sterol regulatory element–binding protein gene (sreA) identified in Penicillium digitatum is required for prochloraz resistance, full virulence and erg11 (cyp51) regulation. PLoS One. 2015;10:e0117115.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang JL, Yu JH, Liu J, Yuan YZ, Li N, He MQ, Qi T, Hui G, Li X, Liu DL. Novel mutations in CYP51B from Penicillium digitatum involved in prochloraz resistance. J Microbiol. 2014;52:762–70.View ArticlePubMedGoogle Scholar
- Nes WD. Biosynthesis of cholesterol and other sterols. Chem Rev. 2011;111:6423–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Nes WD, Hanners PK, Parish EJ. Control of fungal sterol C-24 transalkylation: importance to developmental regulation. Biochem Biophys Res Commun. 1986;139:410–5.View ArticlePubMedGoogle Scholar
- Liu J, Wang SQ, Qin TT, Li N, Niu YH, Li DD, Yuan YZ, Geng H, Xiong L, Liu D. Whole transcriptome analysis of Penicillium digitatum strains treatmented with prochloraz reveals their drug–resistant mechanisms. BMC Genomics. 2015;16:1.View ArticleGoogle Scholar
- Jensen–Pergakes KL, Kennedy MA, Lees ND, Barbuch R, Koegel C, Bard M. Sequencing, disruption, and characterization of the Candida albicans sterol methyltransferase (ERG6) gene: drug susceptibility studies in erg6 Mutants. Antimicrob Agents Chemother. 1998;42:1160–7.PubMedPubMed CentralGoogle Scholar
- Veen M, Stahl U, Lang C. Combined overexpression of genes of the ergosterol biosynthetic pathway leads to accumulation of sterols in Saccharomyces cerevisiae. FEMS Yeast Res. 2003;4:87–95.View ArticlePubMedGoogle Scholar
- da Silva BN, Nakassugi LP, Oliveira JFP, Kohiyama CY, Mossini SAG, Grespan R, Nerilo SB, Mallmann CA, Filho BAA, Machinski Jr M. Antifungal activity and inhibition of fumonisin production by Rosmarinus officinalis L. essential oil in Fusarium verticillioides (Sacc.) Nirenberg. Food Chem. 2015;166:330–6.View ArticleGoogle Scholar
- Tian J, Wang YZ, Zeng H, Li ZY, Zhang P, Tessema A, Peng X. Efficacy and possible mechanisms of perillaldehyde in control of Aspergillus niger causing grape decay. Int J Food Microbiol. 2015;202:27–34.View ArticlePubMedGoogle Scholar
- Xu X, Gao Y, Sun L. Free and esterified triterpene alcohol composition of bee pollen from different botanical origins. Food Res Int. 2012;48:650–6.View ArticleGoogle Scholar
- Chang WQ, Zhang M, Li Y, Li XB, Gao YH, Xie ZY, Lou HX. Lichen endophyte derived pyridoxatin inactivates Candida growth by interfering with ergosterol biosynthesis. Biochim Biophys Acta. 2015;1850:1762–71.View ArticlePubMedGoogle Scholar
- Njombolwana NS, Erasmus A, Fourie PH. Evaluation of curative and protective control of Penicillium digitatum following imazalil application in wax coating. Postharvest Biol Tec. 2013;77:102–10.View ArticleGoogle Scholar
- Kagan IA, Michel A, Prause A, Scheffler BE, Pace P, Duke SO. Gene transcription profiles of Saccharomyces cerevisiae after treatment with plant protection fungicides that inhibit ergosterol biosynthesis. Pestic Biochem Phys. 2005;82:133–53.View ArticleGoogle Scholar
- Sanglard D, Ischer F, Parkinson T, Falconer D, Bille J. Candida albicans mutations in the ergosterol biosynthetic pathway and resistance to several antifungal agents. Antimicrob Agents Chemother. 2003;47:2404–12.View ArticlePubMedPubMed CentralGoogle Scholar
- Khodavandi A, Alizadeh F, Vanda NA, Karimi G, Chong PP. Possible mechanisms of the antifungal activity of fluconazole in combination with terbinafine against Candida albicans. Pharm Biol. 2014;52:1505–9.View ArticlePubMedGoogle Scholar