Pathogenic fungi including Penicillium digitatum and Penicillium italicum are the main destructive pathogens in the citrus industry, causing great losses during postharvest process. To our knowledge, only one mycovirus from P. digitatum has been reported, and the prevalence of such mycoviruses against citrus postharvest pathogenic fungi and their genotyping were still under investigation. In the present study, we showed that 39 of 152 Penicillium isolates from main citrus-growing areas in China were infected with various mycoviruses belonging to polymycoviruses, Narna-like viruses, and families Totiviridae, Partitivirdae and Chrysoviridae. The next generation sequencing (NGS) towards virus genome library and the following molecular analysis revealed two novel mycoviruses Penicillium digitatum polymycovirus 1 (PdPmV1) and Penicillium digitatum Narna-like virus 1 (PdNLV1), coexisting in P. digitatum strain HS-RH2. The fungicide-resistant P. digitatum strains HS-F6 and HS-E9 coinfected by PdPmV1 and PdNLV1 exhibited obvious reduction in triazole drug prochloraz resistance by mycelial growth analysis on both PDA plates and citrus fruit epidermis with given prochloraz concentration. This report at the first time characterized two novel mycoviruses from P. digitatum and revealed the mycovirus-induced reduction of fungicide resistance.

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.

Penicillium digitatum is the most destructive postharvest pathogen of citrus fruits, causing fruit decay and economic loss. Additionally, control of the disease is further complicated by the emergence of drug-resistant strains due to the extensive use of triazole antifungal drugs. In this work, an orthologus gene encoding a putative sterol regulatory element-binding protein (SREBP) was identified in the genome of P. digitatum and named sreA. The putative SreA protein contains a conserved domain of unknown function (DUF2014) at its carboxyl terminus and a helix-loop-helix (HLH) leucine zipper DNA binding domain at its amino terminus, domains that are functionally associated with SREBP transcription factors. The deletion of sreA (ΔsreA) in a prochloraz-resistant strain (PdHS-F6) by Agrobacterium tumefaciens-mediated transformation led to increased susceptibility to prochloraz and a significantly lower EC₅₀ value compared with the HS-F6 wild-type or complementation strain (COsreA). A virulence assay showed that the ΔsreA strain was defective in virulence towards citrus fruits, while the complementation of sreA could restore the virulence to a large extent. Further analysis by quantitative real-time PCR demonstrated that prochloraz-induced expression of cyp51A and cyp51B in PdHS-F6 was completely abolished in the ΔsreA strain. These results demonstrate that sreA is a critical transcription factor gene required for prochloraz resistance and full virulence in P. digitatum and is involved in the regulation of cyp51 expression.

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.