Transcriptome analysis of fungicide-responsive gene expression profiles in two Penicillium italicum strains with different response to the sterol demethylation inhibitor (DMI) fungicide prochloraz

Abstract Background Penicillium italicum (blue mold) is one of citrus pathogens causing undesirable citrus fruit decay even at strictly-controlled low temperatures (< 10 °C) during shipping and storage. P. italicum isolates with considerably high resistance to sterol demethylation inhibitor (DMI) fungicides have emerged; however, mechanism(s) underlying such DMI-resistance remains unclear. In contrast to available elucidation on anti-DMI mechanism for P. digitatum (green mold), how P. italicum DMI-resistance develops has not yet been clarified. Results The present study prepared RNA-sequencing (RNA-seq) libraries for two P. italicum strains (highly resistant (Pi-R) versus highly sensitive (Pi-S) to DMI fungicides), with and without prochloraz treatment, to identify prochloraz-responsive genes facilitating DMI-resistance. After 6 h prochloraz-treatment, comparative transcriptome profiling showed more differentially expressed genes (DEGs) in Pi-R than Pi-S. Functional enrichments identified 15 DEGs in the prochloraz-induced Pi-R transcriptome, simultaneously up-regulated in P. italicum resistance. These included ATP-binding cassette (ABC) transporter-encoding genes, major facilitator superfamily (MFS) transporter-encoding genes, ergosterol (ERG) anabolism component genes ERG2, ERG6 and EGR11 (CYP51A), mitogen-activated protein kinase (MAPK) signaling-inducer genes Mkk1 and Hog1, and Ca2+/calmodulin-dependent kinase (CaMK) signaling-inducer genes CaMK1 and CaMK2. Fragments Per Kilobase per Million mapped reads (FPKM) analysis of Pi-R transcrtiptome showed that prochloraz induced mRNA increase of additional 4 unigenes, including the other two ERG11 isoforms CYP51B and CYP51C and the remaining kinase-encoding genes (i.e., Bck1 and Slt2) required for Slt2-MAPK signaling. The expression patterns of all the 19 prochloraz-responsive genes, obtained in our RNA-seq data sets, have been validated by quantitative real-time PCR (qRT-PCR). These lines of evidence in together draw a general portrait of anti-DMI mechanisms for P. italicum species. Intriguingly, some strategies adopted by the present Pi-R were not observed in the previously documented prochloraz-resistant P. digitatum transcrtiptomes. These included simultaneous induction of all major EGR11 isoforms (CYP51A/B/C), over-expression of ERG2 and ERG6 to modulate ergosterol anabolism, and concurrent mobilization of Slt2-MAPK and CaMK signaling processes to overcome fungicide-induced stresses. Conclusions The present findings provided transcriptomic evidence on P. italicum DMI-resistance mechanisms and revealed some diversity in anti-DMI strategies between P. italicum and P. digitatum species, contributing to our knowledge on P. italicum DMI-resistance mechanisms.