Penicillium digitatum is one of the most destructive postharvest pathogen of citrus fruits, causing fruit decay and economic loss. The emergence of fungicide-resistant strains made the control of P. digitatum more difficult. While the genome of P. digitatum is available, there has been few reports about its resistant mechanism from the transcriptome perspective and there has been no large-scale functional annotation of the genome using expressed genes derived from transcriptomes. Total RNA of P. digitatum strain HS-F6 (prochloraz-resistant strain) and HS-E3 (prochloraz-susceptible strain) before and after prochloraz-treatment were extracted and sequenced on an Illumina Hiseq 2000 platform. The transcriptome data of four samples were compared and analyzed using differential expression analysis, novel transcripts prediction and alternative splicing analysis, SNP analysis and quantitative real-time PCR. We present a large scale analysis about the transcriptome data of P. digitatum. The whole RNA was extracted from a prochloraz-resistant strain (HS-F6) and a prochloraz-susceptible strain (HS-E3) before and after prochloraz-treatment and sequenced by Illumina technology. A total of more than 100 million reads were generated and de novo assembled into 9760 transcripts that contained annotated genes after quality control and sequence assembling. 6625 single nucleotide variations (SNVs) were identified from the sequences aligned against the reference genome. Gene expression profiling analysis was performed upon prochloraz treatment in HS-F6 and HS-E3, and differential expression analysis was used to identify genes related to prochloraz-response and drug-resistance: there are 224 differentially expressed genes in HS-E3 and 1100 differentially expressed genes in HS-F6 after prochloraz-treatment. Moreover, gene expression profile in prochloraz-resistant strain HS-F6 is quite different from that in HS-E3 before prochloraz-treatment, 1520 differential expression genes were identified between the two strains. Gene ontology (GO) term enrichment and KEGG enrichment were then performed to classify the differential expression genes. Among these genes, there are a lot of transporter encoding genes including 14 MFS (Major Facilitator Superfamily) transporters, 8 ABC (ATP-binding cassette transporter) and 3 MATE (multidrug and toxic compound extrusion family) transporters. Meanwhile, the roles of typical MFS, ABC and MATE proteins in prochloraz resistance were investigated using real-time quantitative PCR. The sequencing-based transcriptome data of P. digitatum demonstrate differences between prochloraz-resistant and prochloraz-susceptible strains with prochloraz-treatment. The differences existed in expressed transcripts, splice isoforms and GO categories, which would contribute to our knowledge on the molecular mechanisms involved in drug resistance of P. digitatum.

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. italicumisolates 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 forP. digitatum(green mold), howP. italicumDMI-resistance develops has not yet been clarified.Results: The present study prepared RNA-sequencing (RNA-seq) libraries for twoP. italicumstrains (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 inP. italicumresistance. These included ATP-binding cassette (ABC) transporter-encoding genes, major facilitator superfamily (MFS) transporter-encoding genes, ergosterol (ERG) anabolism component genesERG2,ERG6andEGR11(CYP51A), mitogen-activated protein kinase (MAPK) signaling-inducer genesMkk1andHog1, and Ca²⁺/calmodulin-dependent kinase (CaMK) signaling-inducer genesCaMK1andCaMK2. 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 twoERG11isoformsCYP51BandCYP51Cand the remaining kinase-encoding genes (i.e.,Bck1andSlt2) 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 forP. italicumspecies. Intriguingly, some strategies adopted by the present Pi-R were not observed in the previously documented prochloraz-resistantP. digitatumtranscrtiptomes. These included simultaneous induction of all majorEGR11isoforms (CYP51A/B/C), over-expression ofERG2andERG6to 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 onP. italicumDMI-resistance mechanisms and revealed some diversity in anti-DMI strategies betweenP. italicumandP. digitatumspecies, contributing to our knowledge onP. italicumDMI-resistance mechanisms.

Penicillium sp. are damaging to a range of foods and fruits including citrus. To date, double-stranded (ds)RNA viruses have been reported in most Penicillium species but not in citrus pathogen P. crustosum. Here we report a novel dsRNA virus, designated as Penicillium crustosum chrysovirus 1 (PcCV1) and isolated from P. crustosum strain HS-CQ15. PcCV1 genome comprises four dsRNA segments, referred to as dsRNA1, dsRNA2, dsRNA3, and dsRNA4, which are 3600, 3177, 3078, and 2808 bp in length, respectively. Sequence analysis revealed the presence of four open reading frames (ORFs) in the PcCV1 genome. ORF1 in dsRNA1 encodes a putative RNA-dependent RNA polymerase (RdRp) and ORF2 in dsRNA2 encodes a putative coat protein (CP). The two remaining ORFs, ORF3 in dsRNA3 and ORF4 in dsRNA4, encode proteins of unknown function. Phylogenetic analysis based on RdRp sequences showed that PcCV1 clusters with other members of the genus Chrysovirus, family Chrysoviridae. Transmission electron microscope (TEM) analysis revealed that the PcCV1 visions are approximately 40 nm in diameter. Regarding biological effects of PcCV1, HS-CQ15 harboring the chrysovirus exhibited no obvious difference in colony morphology under fungicide-free conditions but decreased resistance to demethylation inhibitor (DMI)-fungicide prochloraz, as compared to PcCV1-cured strain. Here we provide the first evidence of a virus present in citrus pathogenic fungus P. crustosum and the chrysovirus-induced change in fungicide-resistance of its host fungus.

Calcium (Ca²⁺)/calmodulin-dependent protein kinases (CaMKs) act as a class of crucial elements in Ca²⁺-signal transduction pathways that regulate fungal growth, sporulation, virulence, and environmental stress tolerance. However, little is known about the function of such protein kinase in phytopathogenic Penicillium species. In the present study, a new CaMK gene from the citrus pathogenic fungus P. italicum, designated PiCaMK1, was cloned and functionally characterized by gene knockout and transcriptome analysis. The open reading frame of PiCaMK1 is 1209 bp in full length, which encodes 402 amino acid residues (putative molecular weight ~45.2 KD) with the highest homologous (~96.3%) to the P. expansum CaMK. The knockout mutant ΔPiCaMK1 showed a significant reduction in vegetative growth, conidiation, and virulence (i.e., to induce blue mold decay on citrus fruit). ΔPiCaMK1 was less sensitive to NaCl- or KCl-induced salinity stress and less resistant to mannitol-induced osmotic stress, indicating the functional involvement of PiCaMK1 in such environmental stress tolerance. In contrast, the PiCaMK1-complemented strain ΔPiCaMK1COM can restore all the defective phenotypes. Transcriptome analysis revealed that knockout of PiCaMK1 down-regulated expression of the genes involved in DNA replication and repair, cell cycle, meiosis, pyrimidine and purine metabolisms, and MAPK signaling pathway. Our results suggested the critical role of PiCaMK1 in regulating multiple physical and cellular processes of citrus postharvest pathogen P. italicum, including growth, conidiation, virulence, and environmental stress tolerance.

Sterol 14α-demethylases from Cytochrome P450 family (CYP51s) are essential enzymes in sterol biosynthesis and well-known as the target of antifungal drugs. The 3D structure of CYP51A from Penicillium italicum (PiCYP51A) was constructed through homology modeling based on the crystal structure of human CYP51A (PDB: 3LD6). Molecular dynamics (MD) simulation was operated to relax the initial model and followed by quality assessment using PROCHECK program. On the basis of the docking information on the currently available CYP51s with the patent demethylase inhibitors (DMIs), pharmacophore-based virtual screening combined with docking analysis was performed to pick out twelve new compounds from ZINC database. Six hits revealed in the ligand database suggested potential ability to inhibit PiCYP51A. Compared to patent fungicide triazolone, the top three lead compounds had similar or higher affinity with the target enzyme, and accordingly, exhibited comparable or lower EC₅₀ values to P. italicum isolates. The results could provide references for de novo antifungal drug design.