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

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. 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. 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.

[1]  Satoru Miyano,et al.  Open source clustering software , 2004 .

[2]  K. Chung,et al.  The glutathione peroxidase-mediated reactive oxygen species resistance, fungicide sensitivity and cell wall construction in the citrus fungal pathogen Alternaria alternata. , 2016, Environmental microbiology.

[3]  C. P. Semighini,et al.  Quantitative Analysis of the Relative Transcript Levels of ABC Transporter Atr Genes in Aspergillus nidulans by Real-Time Reverse Transcription-PCR Assay , 2002, Applied and Environmental Microbiology.

[4]  Matthew D. Young,et al.  Gene ontology analysis for RNA-seq: accounting for selection bias , 2010, Genome Biology.

[5]  M. Cyert,et al.  Genetic analysis of calmodulin and its targets in Saccharomyces cerevisiae. , 2001, Annual review of genetics.

[6]  J. Meis,et al.  A Novel Y319H Substitution in CYP51C Associated with Azole Resistance in Aspergillus flavus , 2015, Antimicrobial Agents and Chemotherapy.

[7]  S. Katiyar,et al.  Aspergillus fumigatus CYP51 sequence: potential basis for fluconazole resistance. , 2001, Medical mycology.

[8]  K. Chung,et al.  How the Pathogenic Fungus Alternaria alternata Copes with Stress via the Response Regulators SSK1 and SHO1 , 2016, PloS one.

[9]  N. Xu,et al.  The type II Ca2+/calmodulin-dependent protein kinases are involved in the regulation of cell wall integrity and oxidative stress response in Candida albicans. , 2014, Biochemical and biophysical research communications.

[10]  Zhonghua Ma,et al.  PdMfs1, a major facilitator superfamily transporter from Penicillium digitatum, is partially involved in the imazalil-resistance and pathogenicity , 2012 .

[11]  A. H. McKay,et al.  Toxicity and Resistance Potential of Selected Fungicides to Galactomyces and Penicillium spp. Causing Postharvest Fruit Decays of Citrus and Other Crops. , 2012, Plant disease.

[12]  D. Kelly,et al.  Azole fungicides - understanding resistance mechanisms in agricultural fungal pathogens. , 2015, Pest management science.

[13]  Jing Gao,et al.  A newly identified amino acid substitution T123I in the 14&agr;‐demethylase (Erg11p) of Candida albicans confers azole resistance , 2017, FEMS yeast research.

[14]  D. Margosan,et al.  Characterization of fungicide-resistant isolates of Penicillium digitatum collected in California , 2007 .

[15]  Y. van de Peer,et al.  RNA-sequencing of Cercospora beticola DMI-sensitive and -resistant isolates after treatment with tetraconazole identifies common and contrasting pathway induction. , 2016, Fungal genetics and biology : FG & B.

[16]  P. Janssen,et al.  Mutation in cytochrome P-450-dependent 14 alpha-demethylase results in decreased affinity for azole antifungals. , 1990, Biochemical Society transactions.

[17]  V. Bus.  ED50 levels of Penicillium digitatum and P. italicum with reduced sensitivity to thiabendazole, benomyl and imazalil , 1992 .

[18]  E. Kothe,et al.  A transporter for abiotic stress and plant metabolite resistance in the ectomycorrhizal fungus Tricholoma vaccinum , 2015, Environmental Science and Pollution Research.

[19]  T. Gabaldón,et al.  Genome, Transcriptome, and Functional Analyses of Penicillium expansum Provide New Insights Into Secondary Metabolism and Pathogenicity. , 2015, Molecular plant-microbe interactions : MPMI.

[20]  K. Chung,et al.  A Major Facilitator Superfamily Transporter-Mediated Resistance to Oxidative Stress and Fungicides Requires Yap1, Skn7, and MAP Kinases in the Citrus Fungal Pathogen Alternaria alternata , 2017, PloS one.

[21]  G. Holmes,et al.  Sensitivity of Penicillium digitatum and P. italicum to Postharvest Citrus Fungicides in California. , 1999, Phytopathology.

[22]  H. Hamamoto,et al.  Tandem Repeat of a Transcriptional Enhancer Upstream of the Sterol 14α-Demethylase Gene (CYP51) inPenicillium digitatum , 2000, Applied and Environmental Microbiology.

[23]  Dominique Sanglard,et al.  Amino Acid Substitutions in the Cytochrome P-450 Lanosterol 14α-Demethylase (CYP51A1) from Azole-Resistant Candida albicans Clinical Isolates Contribute to Resistance to Azole Antifungal Agents , 1998, Antimicrobial Agents and Chemotherapy.

[24]  Zhonghua Ma,et al.  Gene transcription profiling of Fusarium graminearum treated with an azole fungicide tebuconazole , 2009, Applied Microbiology and Biotechnology.

[25]  H. Judelson,et al.  Transgene-Induced Silencing of the Zoosporogenesis-Specific NIFC Gene Cluster of Phytophthora infestans Involves Chromatin Alterations , 2007, Eukaryotic Cell.

[26]  N. Zhang,et al.  Nonapoptotic Death of Saccharomyces cerevisiae Cells That Is Stimulated by Hsp90 and Inhibited by Calcineurin and Cmk2 in Response to Endoplasmic Reticulum Stresses , 2008, Eukaryotic Cell.

[27]  S Miyano,et al.  Open source clustering software. , 2004, Bioinformatics.

[28]  T. Foster,et al.  ERG2 and ERG24 Are Required for Normal Vacuolar Physiology as Well as Candida albicans Pathogenicity in a Murine Model of Disseminated but Not Vaginal Candidiasis , 2015, Eukaryotic Cell.

[29]  T. Michailides,et al.  The Point Mutation G461S in the MfCYP51 Gene is Associated with Tebuconazole Resistance in Monilinia fructicola Populations in Brazil. , 2017, Phytopathology.

[30]  H. Saito,et al.  Regulation of cross-talk in yeast MAPK signaling pathways. , 2010, Current opinion in microbiology.

[31]  F. Abe,et al.  Mechanistic role of ergosterol in membrane rigidity and cycloheximide resistance in Saccharomyces cerevisiae. , 2009, Biochimica et biophysica acta.

[32]  C. Waalwijk,et al.  MgSlt2, a cellular integrity MAP kinase gene of the fungal wheat pathogen Mycosphaerella graminicola, is dispensable for penetration but essential for invasive growth. , 2006, Molecular plant-microbe interactions : MPMI.

[33]  B. Tudzynski,et al.  BcSAK1, a Stress-Activated Mitogen-Activated Protein Kinase, Is Involved in Vegetative Differentiation and Pathogenicity in Botrytis cinerea , 2006, Eukaryotic Cell.

[34]  Mark D. Robinson,et al.  edgeR: a Bioconductor package for differential expression analysis of digital gene expression data , 2009, Bioinform..

[35]  Lior Pachter,et al.  Sequence Analysis , 2020, Definitions.

[36]  G. Kaatz,et al.  Multidrug efflux pumps of Gram-positive bacteria. , 2016, Drug resistance updates : reviews and commentaries in antimicrobial and anticancer chemotherapy.

[37]  K. Hammond-Kosack,et al.  The adaptation of Fusarium culmorum to DMI Fungicides Is Mediated by Major Transcriptome Modifications in Response to Azole Fungicide, Including the Overexpression of a PDR Transporter (FcABC1) , 2018, Front. Microbiol..

[38]  K. Kuroda,et al.  Activation of signaling pathways related to cell wall integrity and multidrug resistance by organic solvent in Saccharomyces cerevisiae , 2013, Current Genetics.

[39]  Shih-Hsiung Wu,et al.  Site-specific His/Asp phosphoproteomic analysis of prokaryotes reveals putative targets for drug resistance , 2017, BMC Microbiology.

[40]  P. Le Pape,et al.  The amino acid substitution N136Y in Candida albicans sterol 14alpha-demethylase is involved in fluconazole resistance. , 2016, Medical mycology.

[41]  Jin-Qun Huang,et al.  Comparative genomic, transcriptomic and secretomic profiling of Penicillium oxalicum HP7-1 and its cellulase and xylanase hyper-producing mutant EU2106, and identification of two novel regulatory genes of cellulase and xylanase gene expression , 2016, Biotechnology for Biofuels.

[42]  Hongye Li,et al.  The citrus postharvest pathogen Penicillium digitatum depends on the PdMpkB kinase for developmental and virulence functions. , 2016, International journal of food microbiology.

[43]  Khyati Kapoor,et al.  Multidrug resistance in yeast Candida. , 2005, International review of cytology.

[44]  K. Chung,et al.  Resistance to oxidative stress via regulating siderophore-mediated iron acquisition by the citrus fungal pathogen Alternaria alternata. , 2014, Microbiology.

[45]  P. Fourie,et al.  Imazalil resistance in Penicillium digitatum and P. italicum causing citrus postharvest green and blue mould: Impact and options , 2015 .

[46]  J. Heesemann,et al.  AfMkk2 is required for cell wall integrity signaling, adhesion, and full virulence of the human pathogen Aspergillus fumigatus. , 2010, International journal of medical microbiology : IJMM.

[47]  J. Tyndall,et al.  Impact of Homologous Resistance Mutations from Pathogenic Yeast on Saccharomyces cerevisiae Lanosterol 14α-Demethylase , 2017, Antimicrobial Agents and Chemotherapy.

[48]  Hongye Li,et al.  Functional analysis of two sterol regulatory element binding proteins in Penicillium digitatum , 2017, PloS one.

[49]  J. Cleary,et al.  Genome-wide Expression Profiling of the Response to Polyene, Pyrimidine, Azole, and Echinocandin Antifungal Agents in Saccharomyces cerevisiae* , 2003, Journal of Biological Chemistry.

[50]  H. Hamamoto,et al.  A novel ABC transporter gene, PMR5, is involved in multidrug resistance in the phytopathogenic fungus Penicillium digitatum , 2002, Molecular Genetics and Genomics.

[51]  A. Serrano,et al.  8-Dehydrosterols induce membrane traffic and autophagy defects through V-ATPase dysfunction in Saccharomyces cerevisae. , 2015, Biochimica et biophysica acta.

[52]  B. McDonald,et al.  QTL mapping of fungicide sensitivity reveals novel genes and pleiotropy with melanization in the pathogen Zymoseptoria tritici. , 2015, Fungal genetics and biology : FG & B.

[53]  M. Tomiyama,et al.  A Novel ATP-Binding Cassette Transporter Involved in Multidrug Resistance in the Phytopathogenic FungusPenicillium digitatum , 1998, Applied and Environmental Microbiology.

[54]  Tao Cai,et al.  Automated genome annotation and pathway identification using the KEGG Orthology (KO) as a controlled vocabulary , 2005, Bioinform..

[55]  N. Hawkins,et al.  Constraints on the evolution of azole resistance in plant pathogenic fungi , 2013 .

[56]  R. Calderone,et al.  Azole Susceptibility and Transcriptome Profiling in Candida albicans Mitochondrial Electron Transport Chain Complex I Mutants , 2012, Antimicrobial Agents and Chemotherapy.

[57]  M. Stephens,et al.  RNA-seq: an assessment of technical reproducibility and comparison with gene expression arrays. , 2008, Genome research.

[58]  Ravi Kumar,et al.  Calcium/calmodulin-dependent kinases are involved in growth, thermotolerance, oxidative stress survival, and fertility in Neurospora crassa , 2014, Archives of Microbiology.

[59]  Yongze Yuan,et al.  A novel major facilitator superfamily transporter in Penicillium digitatum (PdMFS2) is required for prochloraz resistance, conidiation and full virulence , 2016, Biotechnology Letters.

[60]  Hongye Li,et al.  Os2 MAP kinase-mediated osmostress tolerance in Penicillium digitatum is associated with its positive regulation on glycerol synthesis and negative regulation on ergosterol synthesis. , 2014, Microbiological research.

[61]  Karl Kuchler,et al.  The Yeast Protein Kinase C Cell Integrity Pathway Mediates Tolerance to the Antifungal Drug Caspofungin through Activation of Slt2p Mitogen-Activated Protein Kinase Signaling , 2003, Eukaryotic Cell.

[62]  M. Raymond,et al.  Relative Contributions of the Candida albicans ABC Transporters Cdr1p and Cdr2p to Clinical Azole Resistance , 2009, Antimicrobial Agents and Chemotherapy.

[63]  David E. Levin,et al.  Cell Wall Integrity Signaling in Saccharomyces cerevisiae , 2005, Microbiology and Molecular Biology Reviews.

[64]  S. Mahmoudi,et al.  Expression of Major Efflux Pumps in Fluconazole-Resistant Candida albicans. , 2017, Infectious disorders drug targets.

[65]  Yongze Yuan,et al.  Whole transcriptome analysis of Penicillium digitatum strains treatmented with prochloraz reveals their drug-resistant mechanisms , 2015, BMC Genomics.

[66]  C. Zuo,et al.  Proteomic analysis of conidia germination in Fusarium oxysporum f. sp. cubense tropical race 4 reveals new targets in ergosterol biosynthesis pathway for controlling Fusarium wilt of banana , 2015, Applied Microbiology and Biotechnology.

[67]  R. González-Méndez,et al.  Calcium/calmodulin kinase1 and its relation to thermotolerance and HSP90 in Sporothrix schenckii: an RNAi and yeast two-hybrid study , 2011, BMC Microbiology.

[68]  Shaojie Li,et al.  The Hsp90 Co-chaperones Sti1, Aha1, and P23 Regulate Adaptive Responses to Antifungal Azoles , 2016, Front. Microbiol..

[69]  K. Tietjen,et al.  Structural and Functional Elucidation of Yeast Lanosterol 14α-Demethylase in Complex with Agrochemical Antifungals , 2016, PloS one.

[70]  P. D. Rogers,et al.  Jjj1 Is a Negative Regulator of Pdr1-Mediated Fluconazole Resistance in Candida glabrata , 2018, mSphere.

[71]  N. L. van der Weerden,et al.  Activation of stress signalling pathways enhances tolerance of fungi to chemical fungicides and antifungal proteins , 2014, Cellular and Molecular Life Sciences.

[72]  Thomas D. Schmittgen,et al.  Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. , 2001, Methods.

[73]  L. Myers,et al.  Mediator Tail Module Is Required for Tac1-Activated CDR1 Expression and Azole Resistance in Candida albicans , 2017, Antimicrobial Agents and Chemotherapy.

[74]  P. Sánchez-Torres,et al.  PdSlt2 Penicillium digitatum mitogen-activated-protein kinase controls sporulation and virulence during citrus fruit infection. , 2017, Fungal biology.

[75]  K. Irie,et al.  MKK1 and MKK2, which encode Saccharomyces cerevisiae mitogen-activated protein kinase-kinase homologs, function in the pathway mediated by protein kinase C , 1993, Molecular and cellular biology.

[76]  D. Sanglard,et al.  Activity of Isavuconazole and Other Azoles against Candida Clinical Isolates and Yeast Model Systems with Known Azole Resistance Mechanisms , 2015, Antimicrobial Agents and Chemotherapy.

[77]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[78]  Guoxing Jing,et al.  Transcriptional profiling analysis of Penicillium digitatum, the causal agent of citrus green mold, unravels an inhibited ergosterol biosynthesis pathway in response to citral , 2016, BMC Genomics.

[79]  L. Koymans,et al.  Contribution of mutations in the cytochrome P450 14alpha-demethylase (Erg11p, Cyp51p) to azole resistance in Candida albicans. , 1999, Microbiology.

[80]  Junli Liu,et al.  Heterologous expression of the Monilinia fructicola CYP51 (MfCYP51) gene in Pichia pastoris confirms the mode of action of the novel fungicide, SYP-Z048 , 2015, Front. Microbiol..

[81]  Hilde van der Togt,et al.  Publisher's Note , 2003, J. Netw. Comput. Appl..

[82]  V. Bus,et al.  Occurrence of Penicillium digitatum and P. italicum resistant to benomyl, thiabendazole, and imazalil on citrus fruit from different geographic origins , 1991 .

[83]  W. Melchers,et al.  Molecular Epidemiology of Aspergillus fumigatus Isolates Harboring the TR34/L98H Azole Resistance Mechanism , 2012, Journal of Clinical Microbiology.

[84]  L. Kitchen,et al.  Diethylcarbamazine-related antimicrobial activity in Mycobacterium tuberculosis-infected blood. , 1998, The Journal of antimicrobial chemotherapy.

[85]  B. Wong,et al.  Fluconazole Transport into Candida albicans Secretory Vesicles by the Membrane Proteins Cdr1p, Cdr2p, and Mdr1p , 2010, Eukaryotic Cell.

[86]  M. Raymond,et al.  PDR16‐mediated azole resistance in Candida albicans , 2006, Molecular microbiology.

[87]  T. Tsuchiya,et al.  Multidrug efflux transporters in the MATE family. , 2009, Biochimica et biophysica acta.

[88]  PCR-based detection of sterol demethylation inhibitor-resistant strains of Penicillium digitatum. , 2001, Pest management science.

[89]  P. Sánchez-Torres,et al.  Molecular insights into fungicide resistance in sensitive and resistant Penicillium digitatum strains infecting citrus , 2011 .

[90]  D. Hagiwara,et al.  A Novel Zn2-Cys6 Transcription Factor AtrR Plays a Key Role in an Azole Resistance Mechanism of Aspergillus fumigatus by Co-regulating cyp51A and cdr1B Expressions , 2017, PLoS pathogens.

[91]  G. Kema,et al.  MgMfs1, a major facilitator superfamily transporter from the fungal wheat pathogen Mycosphaerella graminicola, is a strong protectant against natural toxic compounds and fungicides. , 2007, Fungal genetics and biology : FG & B.

[92]  S. Fillinger,et al.  The osmosensing signal transduction pathway from Botrytis cinerea regulates cell wall integrity and MAP kinase pathways control melanin biosynthesis with influence of light. , 2011, Fungal genetics and biology : FG & B.

[93]  A. Sanz,et al.  The CWI Pathway: Regulation of the Transcriptional Adaptive Response to Cell Wall Stress in Yeast , 2017, Journal of fungi.

[94]  H. Schoonbeek,et al.  Fungal transporters involved in efflux of natural toxic compounds and fungicides. , 2000, Fungal genetics and biology : FG & B.

[95]  Novel mutations in CYP51B from Penicillium digitatum involved in prochloraz resistance , 2014, Journal of Microbiology.

[96]  J. Adrio,et al.  Antifungals , 2018, Reactions Weekly.

[97]  D. Diekema,et al.  Contributions of Aspergillus fumigatus ATP-Binding Cassette Transporter Proteins to Drug Resistance and Virulence , 2013, Eukaryotic Cell.

[98]  E. Balzi,et al.  Multidrug resistance in Aspergillus nidulans involves novel ATP-binding cassette transporters , 1997, Molecular and General Genetics MGG.

[99]  W. Schäfer,et al.  The stress-activated protein kinase FgOS-2 is a key regulator in the life cycle of the cereal pathogen Fusarium graminearum. , 2012, Molecular plant-microbe interactions : MPMI.

[100]  W. Nierman,et al.  Transcriptome analysis of Aspergillus fumigatus exposed to voriconazole , 2006, Current Genetics.

[101]  D. Margosan,et al.  Imazalil resistance linked to a unique insertion sequence in the PdCYP51 promoter region of Penicillium digitatum , 2007 .

[102]  D. Denning,et al.  The cdr1B efflux transporter is associated with non-cyp51a-mediated itraconazole resistance in Aspergillus fumigatus. , 2013, The Journal of antimicrobial chemotherapy.

[103]  P. Sánchez-Torres,et al.  The transcription factor PdSte12 contributes to Penicillium digitatum virulence during citrus fruit infection , 2017 .