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.

BackgroundPenicillium 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.MethodsTotal 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. ResultsWe 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.ConclusionsThe 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.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-015-2043-x) contains supplementary material, whi...

A novel double-stranded RNA virus designated Penicillium digitatum virus 1 (PdV1) was isolated from the citrus fruit rot pathogen P. digitatum (HS-RH1). The full-length cDNA sequence of the dsRNA/PdV1 (5211bp) possesses two partially overlapping open reading frames, which encode a coat protein (CP) and a putative RNA-dependent RNA polymerase (RdRp), respectively. Phylogenetic analysis based on multiple alignments of the amino acid sequences of the RdRp and CP indicated that PdV1 tentatively belongs to the genus Victorivirus in the Totiviridae family. Electron micrographs of negatively stained viral particles purified from the peak fraction of sucrose density gradient centrifugation showed spherical particles ~35nm in diameter. Transfection experiments with purified virions indicated that PdV1 could reduce the vegetative growth and virulence of P. digitatum strain HS-F6. In summary, we report the first isolation and characterization of a mycovirus from P. digitatum that contributes to the hypovirulence phenotypes of the host strain.

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 EC50 values to P. italicum isolates. The results could provide references for de novo antifungal drug design.

To date, partitiviruses, including gammapartitiviruses, have been extensively studied in various fungal hosts but have not been reported in Penicillium digitatum (also called green mold, the pathogenic fungus infecting citrus). In the present work, we isolated and molecularly characterized a double-stranded RNA (dsRNA) partitivirus from citrus green mold, which we have named "Penicillium digitatum gammapartitivirus 1" (PdGV1). The bisegmented genome of PdGV1 contains two dsRNA segments (dsRNA1 and dsRNA2) with a length of 1795 bp and 1622 bp, respectively. Each of the two genomic dsRNAs contains a single open reading frame encoding a putative RNA-dependent RNA polymerase (RdRp) and a coat protein (CP), respectively. Phylogenetic analysis based on RdRp and CP sequences showed that PdGV1 clustered with mycoviruses belonging to the genus Gammapartitivirus, family Partitiviridae, e.g., Penicillium stoloniferum virus S. The 5'- and 3'-untranslated regions (UTRs) of the PdGV1 genomic dsRNAs both contained unique conserved RNA motifs that have never been found in any other partitivirus. This is the first report of a new gammapartitivirus that infects the citrus-pathogenic fungus P. digitatum.