The completed genome sequence of the pathogenic ascomycete fungus Penicillium digitatum

Transcriptome Analysis Reveals Potential Regulators of DMI Fungicide Resistance in the Citrus Postharvest Pathogen Penicillium digitatum

Green mold, caused by Penicillium digitatum, is the major cause of citrus postharvest decay. Currently, the application of sterol demethylation inhibitor (DMI) fungicide is one of the main control measures to prevent green mold. However, the fungicide-resistance problem in the pathogen P. digitatum is growing. The regulatory mechanism of DMI fungicide resistance in P. digitatum is poorly understood. Here, we first performed transcriptomic analysis of the P. digitatum strain Pdw03 treated with imazalil (IMZ) for 2 and 12 h. A total of 1338 genes were up-regulated and 1635 were down-regulated under IMZ treatment for 2 h compared to control while 1700 were up-regulated and 1661 down-regulated under IMZ treatment for 12 h. The expression of about half of the genes in the ergosterol biosynthesis pathway was affected during IMZ stress. Further analysis identified that 84 of 320 transcription factors (TFs) were differentially expressed at both conditions, making them potential regulators in DMI resistance. To confirm their roles, three differentially expressed TFs were selected to generate disruption mutants using the CRISPR/Cas9 technology. The results showed that two of them had no response to IMZ stress while ∆PdflbC was more sensitive compared with the wild type. However, disruption of PdflbC did not affect the ergosterol content. The defect in IMZ sensitivity of ∆PdflbC was restored by genetic complementation of the mutant with a functional copy of PdflbC. Taken together, our results offer a rich source of information to identify novel regulators in DMI resistance.

Transcriptome Analysis of mfs2-Defective Penicillium digitatum Mutant to Reveal Importance of Pdmfs2 in Developing Fungal Prochloraz Resistance

Demethylation inhibitors (DMIs), including prochloraz, are popular fungicides to control citrus postharvest pathogens such as Penicillium digitatum (green mold). However, many P. digitatum strains have developed prochloraz resistance, which decreases drug efficacy. Specific major facilitator superfamily (MFS) transporter gene mfs2, encoding drug-efflux pump protein MFS2, has been identified in P. digitatum strain F6 (PdF6) to confer fungal strain prochloraz resistance. However, except for the drug-efflux pump function of MFS2, other mechanisms relating to the Pdmfs2 are not fully clear. The present study reported a transcriptome investigation on the mfs2-defective P. digitatum strain. Comparing to the wild-type strain, the mfs2-defective strain showed 717 differentially expressed genes (DEGs) without prochloraz induction, and 1221 DEGs with prochloraz induction. The obtained DEGs included multiple isoforms of MFS transporter-encoding genes, ATP-binding cassette (ABC) transporter-encoding genes, and multidrug and toxic compound extrusion (MATE) family protein-encoding genes. Many of these putative drug-efflux pump protein-encoding genes had significantly lower transcript abundances in the mfs2-defective P. digitatum strain at prochloraz induction, as compared to the wild-type strain, including twenty-two MFS transporter-encoding genes (MFS1 to MFS22), two ABC transporter-encoding genes (ABC1 and ABC2), and three MATE protein-encoding genes (MATE1 to MATE3). The prochloraz induction on special drug-efflux pump protein genes in the wild-type strain was not observed in the mfs2-defective strain, including MFS21, MFS22, ABC2, MATE1, MATE2, and MATE3. On the other hand, the up-regulation of other drug-efflux pump protein genes in the mfs2-defective strain cannot recover the fungal prochloraz resistance, including MFS23, MFS26, MFS27, MFS31, MFS33, and ABC3 to ABC8. The functional enrichment of DEGs based on Kyoto Encyclopedia of Genes and Genomes (KEGG), Clusters of Orthologous Groups (COG), and euKaryotic Orthologous Groups (KOG) database resources suggested some essential contributors to the mfs2-relating prochloraz resistance, including ribosome biosynthesis-related genes, oxidative phosphorylation genes, steroid biosynthesis-related genes, fatty acid and lipid metabolism-related genes, and carbon- and nitrogen-metabolism-related genes. The results indicated that the MFS2 transporter might be involved in the regulation of multiple drug-efflux pump protein gene expressions and multiple metabolism-related gene expressions, thus playing an important role in developing P. digitatum prochloraz resistance.

Inhibitory Properties of Cinnamon Bark Oil against Postharvest Pathogen Penicillium digitatum In Vitro

Penicillium digitatum is a major postharvest pathogen that threatens the global citrus fruit industry and causes great economic losses annually. In the present study, inhibitory properties of cinnamon bark oil (CBO) against P. digitatum in vitro were investigated. Results indicated that 0.03% CBO could efficiently inhibit the spore germination, germ tube elongation, mycelial growth, colonial expansion and conidial accumulation of P. digitatum. The results of fluorescein diacetate (FDA) and MitoTraker Orange (MTO) staining also proved the suppression effects of CBO against P. digitatum. Meanwhile, CBO could inhibit green mold rots induced by P. digitatum in citrus fruit when the working concentration of CBO exceeded 0.06%. In addition, the expressions of 12 genes critical for the growth and virulence of P. digitatum were also significantly regulated under CBO stress. Through a transcriptomic analysis, a total of 1802 common differentially expressed genes (DEGs) were detected in P. digitatum after 4 h and 8 h of CBO treatment. Most of the DEG products were associated with carbohydrate, amino acid and lipid metabolism. They directly or indirectly led to the disturbance of the membrane and the generation of reactive oxygen species (ROS). Our results may deepen the understanding of antifungal properties of CBO against P. digitatum and provide the theoretical foundation to uncover the antifungal mechanism of CBO at the molecular level.

Discovery and Transcriptional Profiling of Penicillium digitatum Genes That Could Promote Fungal Virulence during Citrus Fruit Infection

Green mold caused by Penicillium digitatum (Pers.:Fr.) Sacc is the most prevalent postharvest rot concerning citrus fruits. Using the subtractive suppression hybridization (SSH) technique, different P. digitatum genes have been identified that could be involved in virulence during citrus infection in the early stages, a crucial moment that determines whether the infection progresses or not. To this end, a comparison of two P. digitatum strains with high and low virulence has been carried out. We conducted a study on the gene expression profile of the most relevant genes. The results indicate the importance of transcription and regulation processes as well as enzymes involved in the degradation of the plant cell wall. The most represented expressed sequence tag (EST) was identified as PDIP_11000, associated with the FluG domain, which is putatively involved in the activation of conidiation. It is also worth noting that PDIP_02280 encodes a pectin methyl esterase, a cell wall remodeling protein with a high expression level in the most virulent fungal strains, which is notably induced during citrus infection. Furthermore, within the group with the greatest representation and showing significant induction in the early stages of infection, regulatory proteins (PDIP_68700, PDIP_76160) and a chaperone (PDIP_38040) stand out. To a lesser extent, but not less relevant, it is worth distinguishing different regulatory proteins and transcription factors, such as PDIP_00580, PDIP_49640 and PDIP_78930.

The completed genome sequence of Pestalotiopsis versicolor, a pathogenic ascomycete fungus with implications for bayberry production

The pathogenic fungus Pestalotiopsis versicolor is a major etiological agent of fungal twig blight disease affecting bayberry trees. However, the lack of complete genome sequence information for this crucial pathogenic fungus hinders the molecular and genetic investigation of its pathogenic mechanism. To address this knowledge gap, we have generated the complete genome sequence of P. versicolor strain XJ27, employing a combination of Illumina, PacBio, and Hi-C sequencing technologies. This comprehensive genome sequence, comprising 7 chromosomes with an N50 contig size of 7,275,017 bp, a GC content ratio of 50.16%, and a total size of 50.80 Mb, encompasses 13,971 predicted coding genes. By performing comparative genomic analysis between P. versicolor and the genomes of eleven plant-pathogenic fungi, as well as three closely related fungi within the same group, we have gained initial insights into its evolutionary trajectory, particularly through gene family analysis. These findings shed light on the distinctive characteristics and evolutionary history of P. versicolor. Importantly, the availability of this high-quality genetic resource will serve as a foundational tool for investigating the biology, molecular pathogenesis, and virulence of P. versicolor. Furthermore, it will facilitate the development of more potent antifungal medications by uncovering potential vulnerabilities in its genetic makeup.

Genomic insights into Aspergillus sydowii 29R-4-F02: unraveling adaptive mechanisms in subseafloor coal-bearing sediment environments

Introduction

Aspergillussydowii is an important filamentous fungus that inhabits diverse environments. However, investigations on the biology and genetics of A. sydowii in subseafloor sediments remain limited.

Methods

Here, we performed de novo sequencing and assembly of the A. sydowii 29R-4-F02 genome, an isolate obtained from approximately 2.4 km deep, 20-million-year-old coal-bearing sediments beneath the seafloor by employing the Nanopore sequencing platform.

Results and Discussion

The generated genome was 37.19 Mb with GC content of 50.05%. The final assembly consisted of 11 contigs with N50 of 4.6 Mb, encoding 12,488 putative genes. Notably, the subseafloor strain 29R-4-F02 showed a higher number of carbohydrate-active enzymes (CAZymes) and distinct genes related to vesicular fusion and autophagy compared to the terrestrial strain CBS593.65. Furthermore, 257 positively selected genes, including those involved in DNA repair and CAZymes were identified in subseafloor strain 29R-4-F02. These findings suggest that A. sydowii possesses a unique genetic repertoire enabling its survival in the extreme subseafloor environments over tens of millions of years.

Transcriptomic Profile of Penicillium digitatum Reveals Novel Aspects of the Mode of Action of the Antifungal Protein AfpB

Antifungal proteins (AFPs) from filamentous fungi are promising biomolecules to control fungal pathogens. Understanding their biological role and mode of action is essential for their future application. AfpB from the citrus fruit pathogen Penicillium digitatum is highly active against fungal phytopathogens, including its native fungus. Our previous data showed that AfpB acts through a multitargeted three-stage process: interaction with the outer mannosylated cell wall, energy-dependent cell internalization, and intracellular actions that result in cell death. Here, we extend these findings by characterizing the functional role of AfpB and its interaction with P. digitatum through transcriptomic studies. For this, we compared the transcriptomic response of AfpB-treated P. digitatum wild type, a ΔafpB mutant, and an AfpB-overproducing strain. Transcriptomic data suggest a multifaceted role for AfpB. Data from the ΔafpB mutant suggested that the afpB gene contributes to the overall homeostasis of the cell. Additionally, these data showed that AfpB represses toxin-encoding genes, and they suggest a link to apoptotic processes. Gene expression and knockout mutants confirmed that genes coding for acetolactate synthase (ALS) and acetolactate decarboxylase (ALD), which belong to the acetoin biosynthetic pathway, contribute to the inhibitory activity of AfpB. Moreover, a gene encoding a previously uncharacterized extracellular tandem repeat peptide (TRP) protein showed high induction in the presence of AfpB, whereas its TRP monomer enhanced AfpB activity. Overall, our study offers a rich source of information to further advance in the characterization of the multifaceted mode of action of AFPs. IMPORTANCE Fungal infections threaten human health worldwide and have a negative impact on food security, damaging crop production and causing animal diseases. At present, only a few classes of fungicides are available due to the complexity of targeting fungi without affecting plant, animal, or human hosts. Moreover, the intensive use of fungicides in agriculture has led to the development of resistance. Therefore, there is an urgent need to develop antifungal biomolecules with new modes of action to fight human-, animal-, and plant-pathogenic fungi. Fungal antifungal proteins (AFPs) offer great potential as new biofungicides to control deleterious fungi. However, current knowledge about their killing mechanism is still limited, which hampers their potential applicability. AfpB from P. digitatum is a promising molecule with potent and specific fungicidal activity. This study further characterizes its mode of action, opening avenues for the development of new antifungals.

Whole-genome and comparative genome analysis of Mucor racemosus C isolated from Yongchuan Douchi

Mucor racemosus is the predominant fungal in the zhiqu stage of the fermentation of Yongchuan Douchi (Mucor-type), which plays an important role in the fermentation process of Yongchuan Douchi. However, there is a lack of information on the genetic analysis of M. racemosus. In this study, we isolated and identified M. racemosus C (accession no JAPEHQ000000000) from Yongchuan Douchi and analyzed the physiological indicators, then genomic information of the strain to perform a comprehensive analysis of its fermentation capacity and safety. M. racemosus C had neutral protease activity up to 68.051 U/mL at 30 °C and alkaline protease activity up to 57.367 U/mL at 25 °C. In addition, comparing the genomic data with the COGs database (NCBI), it was predicted that M. racemosus C undergoes extensive amino acid metabolism, making C suitable for the production of fermented foods (e.g., Douchi, Syoyu, and sufu). Finally, we performed virulence genes and resistance genes analysis, hemolysis experiment, aflatoxins assay, antibiotic resistance assay to evaluate the safety of M. racemosus C, and the results showed that M. racemosus C was safe, non-toxin-producing and non-hemolytic.

Botryosphaeriaceae gene machinery: Correlation between diversity and virulence

The Botryosphaeriaceae family comprises numerous fungal pathogens capable of causing economically meaningful diseases in a wide range of crops. Many of its members can live as endophytes and turn into aggressive pathogens following the onset of environmental stress events. Their ability to cause disease may rely on the production of a broad set of effectors, such as cell wall-degrading enzymes, secondary metabolites, and peptidases. Here, we conducted comparative analyses of 41 genomes representing six Botryosphaeriaceae genera to provide insights into the genetic features linked to pathogenicity and virulence. We show that these Botryosphaeriaceae genomes possess a large diversity of carbohydrate-active enzymes (CAZymes; 128 families) and peptidases (45 families). Botryosphaeria, Neofusicoccum, and Lasiodiplodia presented the highest number of genes encoding CAZymes involved in the degradation of the plant cell wall components. The genus Botryosphaeria also exhibited the highest abundance of secreted CAZymes and peptidases. Generally, the secondary metabolites gene cluster profile was consistent in the Botryosphaeriaceae family, except for Diplodia and Neoscytalidium. At the strain level, Neofusicoccum parvum NpBt67 stood out among all the Botryosphaeriaceae genomes, presenting a higher number of secretome constituents. In contrast, the Diplodia strains showed the lowest richness of the pathogenicity- and virulence-related genes, which may correlate with their low virulence reported in previous studies. Overall, these results contribute to a better understanding of the mechanisms underlying pathogenicity and virulence in remarkable Botryosphaeriaceae species. Our results also support that Botryosphaeriaceae species could be used as an interesting biotechnological tool for lignocellulose fractionation and bioeconomy.

Comparative genomic study of the Penicillium genus elucidates a diverse pangenome and 15 lateral gene transfer events

The Penicillia are known to produce a wide range natural products-some with devastating outcome for the agricultural industry and others with unexploited potential in different applications. However, a large-scale overview of the biosynthetic potential of different species has been lacking. In this study, we sequenced 93 Penicillium isolates and, together with eleven published genomes that hold similar assembly characteristics, we established a species phylogeny as well as defining a Penicillium pangenome. A total of 5612 genes were shared between ≥ 98 isolates corresponding to approximately half of the average number of genes a Penicillium genome holds. We further identified 15 lateral gene transfer events that have occurred in this collection of Penicillium isolates, which might have played an important role, such as niche adaption, in the evolution of these fungi. The comprehensive characterization of the genomic diversity in the Penicillium genus supersedes single-reference genomes, which do not necessarily capture the entire genetic variation.

Genome-wide identification of the Penicillium digitatum bZIP gene family and the roles of one key member, PdatfA

Penicillium digitatum is the most common cause of postharvest decay in citrus fruits around the world. Previous studies revealed that the bZIP gene family plays crucial roles in development, stress adaptation, and pathogenicity in fungi. However, little is known about the bZIP genes in P. digitatum. In this study, we systematically identified the bZIP family in 23 Penicillium species and analyzed their evolutionary relationships. We found that gene loss and gene duplication shaped the evolution of the Penicillium bZIP family. P. digitatum experienced 3 bZIP gene loss events, but with no gene duplication. We subsequently characterized the biological functions of one important member, PdatfA in P. digitatum by constructing the deletion mutant. Results showed that ΔPdatfA exhibited a moderate growth defect, reduced pigmentation, and slightly increased resistance to fungicides iprodione and fludioxonil. However, ΔPdatfA displayed similar rot symptoms to that of the wild type. The ΔPdatfA mycelia were not affected in response to oxidative stress while its conidia showed enhanced resistance due to the upregulation of catalases. Our results provide new insights into the evolution and functions of the bZIP gene family in Penicillium.

Hybrid Assembly Improves Genome Quality and Completeness of Trametes villosa CCMB561 and Reveals a Huge Potential for Lignocellulose Breakdown

Trametes villosa is a wood-decaying fungus with great potential to be used in the bioconversion of agro-industrial residues and to obtain high-value-added products, such as biofuels. Nonetheless, the lack of high-quality genomic data hampers studies investigating genetic mechanisms and metabolic pathways in T. villosa, hindering its application in industry. Herein, applying a hybrid assembly pipeline using short reads (Illumina HiSeq) and long reads (Oxford Nanopore MinION), we obtained a high-quality genome for the T. villosa CCMB561 and investigated its genetic potential for lignocellulose breakdown. The new genome possesses 143 contigs, N50 of 1,009,271 bp, a total length of 46,748,415 bp, 14,540 protein-coding genes, 22 secondary metabolite gene clusters, and 426 genes encoding Carbohydrate-Active enzymes. Our CAZome annotation and comparative genomic analyses of nine Trametes spp. genomes revealed T. villosa CCMB561 as the species with the highest number of genes encoding lignin-modifying enzymes and a wide array of genes encoding proteins for the breakdown of cellulose, hemicellulose, and pectin. These results bring to light the potential of this isolate to be applied in the bioconversion of lignocellulose and will support future studies on the expression, regulation, and evolution of genes, proteins, and metabolic pathways regarding the bioconversion of lignocellulosic residues.