Next Article in Journal
Correction: Wu et al. Rust Fungi on Medicinal Plants in Guizhou Province with Descriptions of Three New Species. J. Fungi 2023, 9, 953
Previous Article in Journal
The Difference in Diversity between Endophytic Microorganisms in White and Grey Zizania latifolia
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mining the Penicillium expansum Genome for Virulence Genes: A Functional-Based Approach to Discover Novel Loci Mediating Blue Mold Decay of Apple Fruit

by
Dianiris Luciano-Rosario
1,*,†,
Hui Peng
2,†,
Verneta L. Gaskins
1,
Jorge M. Fonseca
1,
Nancy P. Keller
3,4 and
Wayne M. Jurick II
1,*
1
Food Quality Laboratory, USDA-ARS, Beltsville, MD 20705, USA
2
Everglades Research and Education Center, Horticultural Sciences Department, University of Florida, Belle Glade, FL 33430, USA
3
Department of Medical Microbiology and Immunology, University of Wisconsin, Madison, WI 53706, USA
4
Department of Plant Pathology, University of Wisconsin, Madison, WI 53706, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
J. Fungi 2023, 9(11), 1066; https://doi.org/10.3390/jof9111066
Submission received: 26 September 2023 / Revised: 26 October 2023 / Accepted: 28 October 2023 / Published: 1 November 2023
(This article belongs to the Section Fungi in Agriculture and Biotechnology)

Abstract

:
Blue mold, a postharvest disease of pome fruits, is caused by the filamentous fungus Penicillium expansum. In addition to the economic losses caused by P. expansum, food safety can be compromised, as this pathogen is mycotoxigenic. In this study, forward and reverse genetic approaches were used to identify genes involved in blue mold infection in apple fruits. For this, we generated a random T-DNA insertional mutant library. A total of 448 transformants were generated and screened for the reduced decay phenotype on apples. Of these mutants, six (T-193, T-275, T-434, T-588, T-625, and T-711) were selected for continued studies and five unique genes were identified of interest. In addition, two deletion mutants (Δt-625 and Δt-588) and a knockdown strain (t-434KD) were generated for three loci. Data show that the ∆t-588 mutant phenocopied the T-DNA insertion mutant and had virulence penalties during apple fruit decay. We hypothesize that this locus encodes a glyoxalase due to bioinformatic predictions, thus contributing to reduced colony diameter when grown in methylglyoxal (MG). This work presents novel members of signaling networks and additional genetic factors that regulate fungal virulence in the blue mold fungus during apple fruit decay.

1. Introduction

Penicillium expansum is a necrotrophic ascomycete and the causal agent of blue mold disease, a postharvest rot of apples and pears [1]. This filamentous fungus infects apples through wounds created during harvesting, handling, or packing, and the symptoms are most prevalent after storage, as this species can grow prolifically at very low temperatures. The prevalence of blue mold causes significant economic losses in the apple industry. For example, in the United States, blue mold is estimated to cause USD 50–250-million-dollar losses annually [1]. P. expansum not only threatens the apple industry, but also poses food safety concerns, as it produces patulin, which is difficult to eliminate from contaminated food products due to not being degraded by pasteurization [2,3]. The presence of patulin in processed food products is regulated by the FDA, the European Union, and the United Nations World Health Organization, which impose slightly different threshold levels, ranging from 10 to 50 ppb. Blue mold disease management relies primarily on fungicide applications, resulting in limited effectiveness due to the development of fungicide-resistant strains [4].
Due to the challenges in blue mold disease management and patulin reduction strategies, previous research has largely focused on understanding patulin biosynthesis and decontamination. While the molecular biology of P. expansum regarding virulence has also been investigated, only a few genes have been identified as significant and potential targets for disease-management potential. Many studies have focused on conserved virulence factors studied across the ascomycota fungi (e.g., creA, veA, laeA), and others have relied on mining-omics datasets [5,6,7,8]. While these approaches are high yielding, a forward genetics approach has not been implemented to date in Penicillium expansum. Historically, forward genetics approaches have produced many new advances toward gene discovery in diverse biological processes, including understanding the genetic mechanisms of disease caused by fungal plant pathogens [9,10]. One approach towards identifying and discovering genes of interest is through a random mutagenesis approach. Previous research in molecular plant pathology has successfully led to the identification of important loci for pathogenesis and/or virulence in diverse species [9,11,12,13,14,15].
This is the first study to generate a P. expansum library of random T-DNA insertional mutants using Agrobacterium-mediated transformation (AMT). We screened 448 of these transformants in apples for virulence defects and identified six single-spore T-DNA mutants (T-193, T-275, T-588/T-711, T-434, and T-625) that showed reduced lesion diameter upon apple inoculation when compared to the wild-type strain. They exhibited little to no growth penalty in vitro and contained a single T-DNA insertion in the genome. After performing TAIL-PCR, we identified the interrupted locus for each transformant. We obtained deletion mutants for two of these loci (∆t-625 and ∆t-588) and a knockdown strain (t-434KD). Strains ∆t-625 and t-434KD did not phenocopy the T-DNA mutants upon apple inoculations in three different cultivars. Interestingly, ∆t-588 showed reduced lesion diameter when compared to the control strain in Honeycrisp apples. None of these mutants showed reduced nor enhanced tolerance to Congo Red or Sorbitol at different concentrations. We hypothesize that the gene interrupted in T-588 likely functions as a glyoxalase, as ∆t-588 shows reduced colony diameter when compared to the control strains upon exposure to methylglyoxal. Hence, a new dimension of the virulence picture has been added to the complex process of decay, in which it is crucial for the fungus to effectively detoxify and/or metabolize toxic byproducts of the primary metabolism. The T-DNA mutant library also provides a functional platform for the fungal genetics community that can be screened for a particular phenotype in vitro or in vivo to understand the genetic basis of any given biological question.

2. Materials and Methods

2.1. Fungal Propagation and Media

Penicillium expansum strain R19, which was preserved in the USDA-ARS Food Quality Lab (FQL), was used in this study. Fungal growth conditions in culture and inoculation of apple fruit (‘Golden Delicious’ and ‘Ginger Gold’) were conducted following the methods used in a previous study by [16]. For the single-gene-deletion mutant generation, auxotrophic strains TDL 9.1 and TWW 12.1 (Pe21 background) were utilized. The strains were cultivated for conidial production in potato dextrose agar (PDA) (Becton, Dickinson and Company, Franklin Lakes, NJ, USA) or glucose minimal media (GMM) (10 g/L glucose, 1× Nitrate Salts, 1× Trace Elements, pH 6.5) at 25 °C for 5 to 7 days. GMM was also prepared using ammonium salts for culturing the knockdown strain TDL 37.1.

2.2. AMT of Penicillium expansum R19

The Agrobacterium tumefaciens-mediated transformation (AMT) of Penicillium expansum R19 was conducted following the protocol described by [17], with minor modifications. The Agrobacterium tumefaciens strain EHA105 was used as a T-DNA acceptor, and the binary plasmid pPK2 carrying the hygromycin B resistance gene (Hph) (Figure S1) was used as a T-DNA donor [18]. The plasmid was transferred into EHA105 cells using the heat-shock method [19]. A single colony of A. tumefacies carrying the plasmid pPK2 was cultured in 5 mL of YEB medium (5 g/L beef extract, 1 g/L yeast extract, 5 g/L peptone, 5 g/L sucrose, 0.5 g/L MgCl2) with kanamycin (100 μg mL−1) and rifampicin (50 μg mL−1) at 28 °C on a rotary shaker at 200 rpm until the OD600 reached 0.6–0.8. Cells from 2 mL of the culture were collected by centrifuge at 4000× g for 5 min and resuspended with 500 μL of liquid-induction medium (IM, 0.6 mg mL−1 MgSO4·7 H2O, 0.3 mg mL−1 NaCl, 0.01 mg mL−1 CaCl2, 0.01 mg mL−1 FeSO4, 0.4 mg mL−1 NH4NO3, 0.5% glycerol, 40 mM MES, 0.5 μg mL−1 H3BO3, 0.5 μg mL−1 ZnSO4·7H2O, 0.5 μg mL−1 CuSO4·5H2O, 0.5 μg mL−1 MnSO4·H2O, 0.5 μg mL−1 Na2MoO4·2H2O, and 0.2 mg mL−1 glucose, 1 mM potassium phosphate buffer, pH 5.8). The supernatant was discarded after the sample was centrifuged at 4000× g for 5 min. The bacterial pellets were diluted to OD600 = 0.15 with the IM supplemented with 200 μM acetosyringone, and then they were grown at 28 °C for 8–12 h on a rotary shaker at 200 rpm until the OD600 reached 0.3. P. expansum spores were harvested from 10-day-old PDA plates and diluted to 106 mL−1 with sterilized water. The fungal-spore suspension and the induced A. tumefaciens culture were gently mixed by pipette at a ratio of 1:1, and set at room temperature for 2 h without shaking. A total of 200 μL of mixture was spread onto the surface of a nitrocellulose filter (0.45 μm) placed on an IM agar plate (liquid IM plus 1.0 mg mL−1 glucose and 1.5% agar) supplemented with 200 μM acetosyringone. The plates were incubated at 25 °C for 60 h. After cocultivation, the membranes were transferred onto a selection medium (PDA) containing 300 ng mL−1 cefotaxime and 100 ng mL−1 hygromycin B. After 4–5 days, visible colonies of transformants were isolated and streaked onto fresh selection plates. Consequent single colonies were selected for further experiments.

2.3. Fruit Wounding and Fungal Inoculation

Apple fruits (‘Golden Delicious’, ‘Ginger Gold’, ‘Fuji’, and ‘Honeycrisp’) were obtained from either the Penn State Fruit Research and Extension Center (Biglersville, PA) or local organic produce vendors. When possible, organic fruits were used, and when that was not possible, fruit without postharvest fungicide treatments or preharvest sprays with postharvest efficacy were procured. Pathogen inoculation was carried out with blemish- and disease-free fruits at a similar size following the published [16] method with minor modifications. Fruits were thoroughly washed with tap water and dried with a paper towels before sanitization with a spray of 75% ethanol. Fruits were dried with tissues and then wounded using a cylindrical wounding tool (3 mm deep × 3 mm diameter). The wounds were inoculated with 10 μL of 104 conidia/mL suspension of P. expansum immediately after wounding. Twenty minutes after inoculation, boxes were loaded with inoculated fruits, covered with lids, and stored at room temperature. The lesion expansion of blue mold decay was measured at 4, 7, and 10 days, or every day for 7 days postinoculation (DPI).

2.4. Standard Molecular Methods

TAIL-PCR

Genomic DNA of the candidate transformants was extracted from mycelia using the Sigma Extract-N-Amp™ Plant PCR Kit (Sigma-Aldrich, St. Louis, MO, USA) according to the manufacturer’s instructions. Thermal asymmetric interlaced PCR (TAIL-PCR) was applied to identify the locus of T-DNA insertion. Degenerated primers included AD1 (5′-wagtgnagwancanaga-3′) and AD2 (5′-ntcgastwtsgwgtt-3′). Based on the DNA sequence of the T-DNA, a total of six specific primers in nested positions close to the left or right border of the T-DNA were designed and synthesized. Specific primers for the left border were LB1 (5′-gtgcctaatgagtgagctaactccc-3′), LB2 (5′-ctcacattaattgcgttgcgctc-3′), and LB3 (5′-gagcaattcggcgttaattcagt-3′). Specific primers for the right border were RB1 (5′-tggcactggccgtcgttttacaac-3′), RB2 (5′-aacgtcgtgactgggaaaaccct-3′), and RB3 (5′-cccttcccaacagttgcgca-3′). For the first-round PCR, components in 10 μL of reaction volume included 0.05 unit of HS ExTaq (TaKaRa Bio, Mountain View, CA, USA), 2 mM dNTP mix, 1 ng of genome DNA, 1 ppM T-DNA-specific primer (LB1 or RB1), and 10 ppM degenerated primer mix (AD1 and AD4). TAIL-PCR was performed according to the protocol used by [20], with some modification. For the secondary TAIL-PCR, the 10 μL of the reaction mix included 1 μL of a 1:100 dilution of the primary PCR products, 1 ppM LB2 or RB2 primer, and the composition of other reagents was the same as that of the primary PCR reaction mix. For the third TAIL-PCR, the 10 μL of reaction mix included 1 μL of a 1:100 dilution of the secondary PCR products and 1 ppm LB3 or RB3 primer, and other components were the same as those in the primary PCR reaction. Products of the third-round PCR were separated on 1% agarose gel with ethidium bromide. Single, prominent bands were purified with the Qiagen® DNA Gel Extraction Kit (Qiagen, Frederick, MD, USA), and consequent DNA fragments were sequenced using T-DNA-specific primer (LB3 or RB3).

2.5. Identification of T-DNA Flanking Insertion

Genomic DNA of the strain R19 (wild-type) and the T-DNA mutants were extracted from mycelia using the Qiagen DNeasy Plant Mini Kit (Qiagen, Frederick, MD, USA). The T-DNA-induced mutation was identified through thermal asymmetric interlaced PCR (TAIL-PCR) according to a published protocol [20]. A total of thirteen primers, including three T-DNA left-border ones (Ppk2-LB1, 2, and 3), three right-border ones (Ppk2-LB1, 2, and 3), and seven arbitrary-degenerate ones (AD1-8), were used in the TAIL-PCR. Three rounds of PCR were run on a Bio-Rad PCR machine (CFX96 Touch™ Real-Time PCR Detection System) with EX Taq DNA Polymerase (TaKaRa Bio, Mountain View, CA, USA) using the product of the previous PCR as a template for the next. A common arbitrary primer and nested T-DNA-specific primers were employed in the PCR in a consecutive manner. Candidate bands from the second or third PCR were purified with the QIAquick PCR Purification Kit (Qiagen, Frederick, MD, USA), and resulting DNA fragments were subcloned with the Invitrogen TOPO TA Cloning Kit (Invitrogen, Frederick, MD, USA). These were sequenced with an ABI 3730 DNA Sequencer within the Iowa State University DNA Facility (Ames, IA, USA).

2.6. Gene Expression Analysis

Total RNA was isolated from frozen tissues using the RNeasy Plant Mini Kit (Qiagen, Frederick, MD, USA) according to the manufacturer’s instructions. Genomic DNA contamination was removed with TURBO DNase (Invitrogen, Frederick, MD, USA) on the column following the manufacturer’s instruction. Purified RNA was quantified with a NanoDrop 1000 spectrophotometer (Thermo Scientific, Waltham, MA, USA). The first strand of cDNA was synthesized with quantified RNA using iScriptTM Transcription Supermix (Bio-RAD, Hercules, CA, USA) following the manual. The primer pairs, CACGGTCTTGCCACTTGCTCGT/TCGTTTCGAGTACGTCGCGCTG and ACGTCGTCCCCATCTACGAG/GCTCAGCGAGGATCTTCATC were used for the expression analysis of PeHSP and the housekeeping gene actin (Genbank accession number XM_016739901). Quantitative PCR was performed with the CFX96 Touch™ Real-Time PCR Detection System (Bio-RAD, Hercules, CA, USA) using GoTaq® qPCR Master Mix (Promega, Durham, NC, USA) and the following program: 95 °C for 2 min, followed by 45 cycles of 95 °C for 5 s and 60 °C for 15 s. A melt–curve analysis was performed to check the specificities of the PCR products. Relative expression of a given locus was normalized to the expression of actin and shown in fold changes (lowest value = 1) using the cycle threshold (Ct) 2−ΔΔCt method [21]. The significance of the expression difference in various tissues at different timepoints was determined by Student’s t-test (p = 0.05) (Microsoft Excel 2010). All PCR experiments were performed in triplicate.

2.7. Southern Blot Analysis

Genomic DNA of the strain R19 and T-DNA mutants were isolated from mycelia using the Qiagen DNeasy Plant Mini Kit (Qiagen, Frederick, MD, USA) and digested with restriction endonuclease PAC I, separated by gel electrophoresis with 1.3% TAE agarose gel, and blotted onto Amersham Hybond-N+ nylon membranes (GE Healthcare Life Sciences, Laurel, MD, USA). A DNA fragment of 283 bp, including a portion of Hph DNA and a portion of Trp terminator, was used as the DNA hybridization probe. This fragment was amplified from the pPK2 plasmid with the primer pair TTCGGGCGTACACAAATCGCCC/AGGCACTCTTTGCTGCTTGGACA. The probes were labeled and detected using the DIG High Prime DNA Labeling and Detection Starter Kit II (Roche, Branchburg, NJ, USA) following the manual. Chemiluminescent signals on the membrane were recorded by Kodak X-ray film (8 × 10 inches).

2.8. Gene Deletion Construct Design

To obtain constructs for generating deletion mutants, primers were designed as described by [22,23] to obtain three distinct flanks. One of these flanks corresponded to the selectable marker pyrG from Aspergillus fumigatus, and the other two flanks corresponded to 1.5 kb upstream and downstream of the open reading frame (ORF) for the desired gene of interest. We used PCR to obtain independent amplicon fragments with primer, as listed in Table S1. After gel purification, double-joint PCR was performed and used as a template DNA for obtaining the constructs by a final PCR reaction. Regarding the knockdown constructs, the same principle was implemented, except that the primers were designed to insert the Aspergillus nidulans promoter niad/niia region upstream of the gene-of-interest’s ORF.

2.9. Polyethylene Glycol-Mediated Protoplast Transformation

The transformation process was carried out as described by [23,24], with some modifications. To obtain protoplasts, 109 spores of strain TDL 9.1 or TWW 12.1 were incubated in 250 mL liquid GMM supplemented with yeast extract (YE), (1 g/L) uridine (0.56 g/L), and uracil (1.26 g/L) for 13.5 h at 25 °C and 280 rpm. After this, the cultures were centrifuged in bottles at 11,500× g for 2 min. After decanting the supernatant, the pellet was washed with sterile dH2O and transferred to a 20 mL tube, which was centrifuged at 25,000× g for 1 min. The supernatant was discarded, and approximately 200 µL of the pellet was transferred to the protoplast solution containing 90 mg of lysing enzymes from Trichoderma harzianum (Sigma-Aldrich, St. Louis, MO, USA), 60 mg of yatalase (TakaraBio, Kusatsu, Shiga, Japan) and 30 mL of Osmotic Medium (1 M MgSO4, 10 mM sodium phosphate buffer, pH 5.8). The solution was incubated in an Erlenmeyer flask for 4 h at 100 rpm at 28 °C. Then, the protoplasts were transferred to a 50 mL conical tube, and 10 mL of trapping buffer (0.6 M sorbitol, 10 mM Tris-HCl, pH 7.5) were added on top. This solution was centrifuged for 15 min at 1300× g. The remaining interphase was collected in a 15 mL conical tube and gently washed by adding 10 mL of STC buffer (1.2 M sorbitol, 10 mM CaCl2, 10 mM Tris-HCl, pH 7.5). This solution was then centrifuged for 7 min at 7800× g. The supernatant was discarded, and the pellet was resuspended in 1 mL of STC, then transferred to a 1.5 mL tube, where it was centrifuged for 15 s in a tabletop centrifuge at room temperature. The supernatant was discarded, and the pellet was resuspended in 300 µL of STC buffer. After this, 100 µL of the protoplast suspension, 15 µL of the purified construct, and 85 µL of STC buffer were incubated on ice for 50 min. After this period, 1.25 mL of a 60% PEG solution (60% PEG 4000, 50 mM CaCl2, 50 mM Tris-HCl, pH 7.5) was added to the tubes and incubated at room temperature for 20 additional minutes. Finally, 5 mL of STC buffer was added to the mixture, which was used to plate the transformants using an overlay method over sorbitol minimal media (SMM). Once the transformants grew, colonies were transferred to individual GMM plates and screened using PCR and Southern Blot.

2.10. Stress Plate Assays

GMM was supplemented with Congo Red (Sigma-Aldrich, St. Louis, MO, USA) at a final concentration of 10 µg/mL, 25 µg/mL, 100 µg/mL, or 200 µg/mL, sorbitol at a final concentration of 1 M, 1.2 M, or 1.5 M, and/or methylglyoxal (Sigma-Aldrich, St. Louis, MO, USA) at a final concentration of 0.1%. The plates were point-inoculated with total 100 spores, and incubated at 25 °C for 5 days for the Congo Red and Sorbitol or 12 days for the methylglyoxal.

2.11. Data Analysis

Statistical significance was analyzed with the t-test at each time point or 2-way ANOVA method, as appropriate. Multiple comparison of values was done with the Tukey HSD test. (* p < 0.05, ** p < 0.01,*** p < 0.001).

3. Results

3.1. Identification of Penicillium expansum AMT Virulence Mutants

A collection of 448 transformants was generated through Agrobacterium-mediated transformation (AMT) using the T-DNA of pPK2hph as the insertional mutagen. All transformants were assessed for their impact on virulence in bioassays using three apples per transformant. Out of the 448 transformants, 50 putative mutants displayed a strong reduction (>20%) in virulence. To specifically isolate virulence-related genes using the candidate mutant pool, only those mutants that displayed virulence reduction on apples greater than a growth-rate change in culture medium were retested, again using three apples per transformant. In this way, 12 virulence mutants with a reproducible virulence defect were identified (Figure 1). Cultivation of these transformants on a rich medium revealed that T-DNA insertion has a specific impact on the morphology of fungi, such as the growth rate and color (Figure 1A). The colony of R19 has a blue ring with white edges. All transformants except for T163 showed various shades of brown pigmentation in media. Among the 12 mutants, T489 had the lowest growth-rate change in the culture medium but displayed 44.4% virulence reduction when compared to wild-type R19 (Figure 1B).

3.2. Determination of T-DNA Integration Pattern

Using a probe derived from the T-DNA, DNA-blotting with the genomic DNA of six selected transformants (T-193, T-275, T-434, T-588, T-625, T-711) was performed. As shown in Figure 2, no band was observed in the lane carrying the DNA of wild-type P. expansum (R19), while one unique band appeared on each lane of the transformants, clearly showing a single insertion of the T-DNA in the genome of examined transformants. Except for T-588 and T-711 (which were inserted into the same locus), other transformant bands have different sizes, suggesting the positional diversity of the T-DNA insertion in those transformants within the P. expansum R19 genome (Figure 2).

3.3. Identification of Loci Interrupted by T-DNA Insertion

Thermal asymmetric interlaced polymerase chain reaction (TAIL-PCR) was implemented on selected virulence mutants to isolate DNA sequences flanking the T-DNA insertion. Using the combination of three T-DNA-specific primers and four different degenerate primers, LB and RB flanking sequences were obtained from transformants T-193 (PEXP_001700), T-434 (PEXP_051540), T-588 (PEXP_003290), T-625 (PEXP_008080), and T-711 (PEXP_003290), while only the LB flanking sequence was identified from T-275 (PEXP_016530). TAIL-PCR was not able to amply the RB of T-DNA in T275, indicating a possible excision failure of the right border during the integration of T-DNA into the R19 genome. DNA sequencing revealed flanking sequences of a border of T-DNA were exclusive for each transformant, indicating the uniqueness of the T-DNA integration (Table S3). Interestingly, T-588 and T-711 were found to share identical RB and LB flanking sequences, suggesting they are derived from the same event of T-DNA integration. We thus did not pursue T-711 for further study. Nevertheless, these flanking sequences allowed us to localize the T-DNA integration sites and identify putative integration-affected genes by comparing the adjacent genome sequence of the T-DNA inserting locus in transformants and P. expansum R19 (Figure 2B, Table 1).
For mutant T-193, a putative ORF located downstream 1518 bp to the insertion locus was found to be the coding sequence closest to the T-DNA. Blasting against the NCBI protein database revealed that the putative protein exhibited the highest percent identity (99.82%) to a U8 small nucleolar RNA-associated protein 15 (RAP15) in Penicillium viridicatum. The putative protein has the UTP15_C motif conserved in other RAP15 members, and is thus named PeRAP15. For T-275, T-DNA was found to be inserted 33 bp upstream of a putative start codon of an ORF. The putative protein encoded by the ORF showed the highest similarity (87.3%) to the glycyl-tRNA synthetase (GlyRS) in Penicillium digitatum according to the BLAST results against the NCBI protein database. No conserved motif was identified for this protein sequence (GlyRS) or from other species that also contained the putative protein. This ORF was thus named PeGlyRS. For T434, T-DNA was inserted into the coding area of a putative gene. Its amino acid sequence displayed high similarity to N-acetylglucosaminyl transferases (GPI1). Alignment of this sequence to other GPI1s showed the sequence contained the conserved domain. Consequently, this putative gene was named PeGPI1. For T588, T-DNA was found to be in the coding sequence of a putative ORF. Blasting against the NCBI database showed that the putative protein had the highest identity to glyoxalases (GloIs) from other species. Like other GloIs, this protein also contained a conserved VOC motif responsible for the isomerization of the spontaneously formed hemithioacetal adduct between GSH and 2-oxoaldehydes (such as methylglyoxal) into S-2-hydroxyacylglutathione [25]. For T-625, T-DNA was inserted into a putative ORF that codes a protein showing high similarity to the heat-shock protein 40 (HSP40) cochaperone Jid1 in P. digitatum. The protein, together with other HSP40s, shared a common DNA-J domain, a typical motif of a chaperone [26]. This mutant, T-625, was investigated and characterized by [16], and shown to be involved in secretion processes, toxin production, and virulence in apple fruits.

3.4. Expression Profiles for Genes up- and Downstream of the T-DNA Insertion

To assess how the T-DNA insertion affects the transcription of candidate genes, we investigated the expression of candidate genes in both wild-type (R19) and mutant fungi from liquid-medium culture and P. expansum-decayed apple fruit tissue at 10 DPI. All candidate genes, including RAP15, GloI, GlyRS, HSP40, and GPI1, are transcribed in fungi from both in vitro culture and decayed apple tissue, indicating they are active in vivo (Figure 3). Among the five genes, RAP15, GlyRS, HSP40, and GPI1 displayed significantly lower expression in the decayed apple tissue than in the culture medium, while GloI showed five-times-higher expression in decayed tissue than in the culture medium. In the culture at 10 DPI, transcription of all candidate genes was heavily reduced (>62%) in transformants. In decayed tissues at 10 DPI, transcription of four genes (RAP15, GloI, GlyRS, and HSP40) was significantly lower (>79%) in transformants, and only the expression of GPI1 was decreased when comparing apple and media culture conditions. This observation suggests T-DNA insertion in both noncoding and coding areas can disrupt the transcription of these genes. Considering that the T-DNA insertion occurred in the coding region in three transformants, including T434, three genes, including GPI1, were expected to have aberrant-sized transcripts, or may not be transcribed due to the T-DNA insertion. Hence, these directly impacted loci serve as functional null mutants.

3.5. Targeted Single-Gene Deletion or Knockdown Strain Phenotypes

To confirm that the observed T-DNA phenotypes are due to the interruption of the identified genes, we aimed to generate deletion strains for all the identified loci. We successfully obtained deletion mutants for the genes that had T-DNA insertions T-625 and T-588. Since obtaining deletion strains for the genes that contained T-DNA insertions in T-193, T-262, or T-434 was not possible, we sought to obtain knockdown strains. We successfully generated a knockdown strain for T-434, but not for T-193 or T-262. The genotypes for the obtained strains are listed in Table S2. Thus, we hypothesize that the deletion or silencing of these loci results in a lethal phenotype. After obtaining the deletion and knockdown strains, we assessed the radial growth and germination initiation rates of these when compared to the control strain, and found no statistically significant differences between our deletion mutants and the control strain (Figure 4A,B). We observed that, when grown on ammonia as a nitrogen source, t-434KD exhibits a slight but statistically significant increase in radial growth when compared to the control strain in the same media. Upon assessing the germination initiation rate for t-434KD using ammonia as the nitrogen source, we did not observe any significant difference between the control and knockdown strains (Figure 4A,B).

3.6. ∆T-588 Strain Shows Reduced Lesion Decay Diameter In Vivo

After generating the deletion and knockdown mutants, these strains were assessed for variation in lesion development in ‘Fuji’, ‘Golden Delicious’, and ‘Honeycrisp’ apples. The primary objective was to test the hypothesis that these genes were involved in virulence during apple fruit decay, as evidenced by the initial screening performed using the T-DNA mutants. Interestingly, neither ∆t-625 nor t-434KD showed reduced lesion diameter when compared to the control strain on any of the tested apple cultivars (Figure 5). As for ∆t-588, it was observed that there was reduced lesion diameter in ‘Honeycrisp’ apples when compared to the control strain. The same trend holds for ∆t-588 in both ‘Fuji’ and ‘Golden Delicious’, although the differences for these cultivars are not statistically significant. Unexpectedly, we observed that t-434KD showed a slight but significant increase in lesion diameter in ‘Fuji’ and ‘Honeycrisp’ cultivars when compared to the control strain (Figure 5). In addition, we subjected these mutants to cell wall and membrane stressors, Congo Red and sorbitol, in different concentrations. None of the strains showed increased or reduced susceptibility to the assessed stressors (Figure 6).

3.7. ∆t-588 Shows Reduced Colony Diameter When Exposed to Methylglyoxal

While bioinformatically assessing possible gene function for the identified loci, it was noted that the locus interrupted in the ∆t-588 mutant contained a glyoxalase domain. Therefore, it was hypothesized that PEXP_003290 could function in detoxifying methylglyoxal. To test this, ∆t-588 and the control strain were cultured in GMM and GMM supplemented with methylglyoxal at a final concentration of 0.1% [27]. After incubating the plates for 14 days, it was observed that the ∆t-588 colony had decreased radial growth when compared to the control strain (Figure 7).

4. Discussion

To develop novel management strategies for blue mold decay, one needs to understand the molecular underpinnings and biochemical pathways that modulate virulence in P. expansum. In this study, a forward and reverse genetics approach was used to discover and characterize new genes important in P. expansum virulence during apple fruit decay. After screening 448 random T-DNA insertion mutants for reduced lesion diameter in apples, we obtained 6 mutants (T-191, T-434, T-262, T-588, T625, and T-711) and identified 5 loci that are involved in P. expansum virulence. Upon identification of the disrupted genes, we aimed to obtain deletion mutants for each of the loci and were successful in obtaining two deletion mutants (∆t-625 and ∆t-588) and one knockdown strain (t-434KD). In general, we found similar technical trends, patterns, and limitations with the AMT system as in other studies [9,28]. For example, a common technical challenge in AMT is its low efficiency and left-border excision failures, which was encountered in this study [28]. In addition, deletion or knockdown strains were unable to be generated for some of the identified loci, which is a common occurrence that has been observed in other fungal pathogens, like Magnaporthe oryzae, Fusarium oxysporum, and Botrytis cinerea.

4.1. UTP15 (T-DNA Mutant T-193)

The U3 small nucleolar RNA-associated protein 15 (UTP15) is involved in the nucleolar processing of pre-18S rRNA. This protein is typically found together with the other two WD proteins, i.e., CIRH1A and WDR43, and are found in the nuclear matrix fraction of HeLa cells. They are contained in the t-UTP subcomplex of SSU processomes, and have been shown to bind directly with each other and to bind tightly to the FC, independent of rRNA transcription, and are involved in ribosome biogenesis [29]. Interestingly, in a forward genetic screening in zebrafish, a mutation in utp15 could affect vascular patterning in a p53-dependent manner in zebrafish embryos [30]. In another study, UTP15 was shown to be upregulated when A. nidulans was treated with abscisic acid (ABA) [31]. Aside from the discovery of this locus in other fungal systems, our work is key in nature regarding this finding in Penicillium spp.

4.2. Glycyl-tRNA Synthetase (T-DNA Mutant T-275)

Glycyl-tRNA synthetase (GlyRS) can charge tRNAs with glycine, and thus plays an important role in protein synthesis, which makes this enzyme ubiquitous in all cell types [32]. It has been reported that GlyRS interacts with elongation factor 1-delta (EEF1D), which participates in the enzymatic delivery of aminoacyl tRNAs to the ribosome as a subunit of the elongation factor-1 complex [33]. GlyRS has been shown to exert diverse functions, like acting as gene regulators, involvement in intron splicing, and angiogenesis, in addition to their role in translation [34]. Conserved signature indels (CSIs) were identified in GlyRS protein sequences across Xanthomonadales and a limited number of Beta- or Gamma-proteobacteria that mostly act as plant or human pathogens, indicating the potential involvement of GlyRS in pathogenesis [35]. Translative efforts have been directed at the development of new antibiotics, which would inhibit the activity of aminoacyl-tRNA synthetases, including GlyRS from pathogens such as Enterococcus faecalis, Staphylococcus aureus, Streptococcus pneumoniae, Helicobacter pylori, Mycobacterium tuberculosis, or Candida albicans [36,37]. Recently, the characterization and function of GlyRS proteins in plant-specific fungal pathogens have not been reported. Our finding indicates the potential involvement of the GlyRS protein in the virulence process during apple decay in P. expansum. Systematic functional analysis of GlyRS protein genes could shed additional light on the rationale of this class of proteins in the fungal virulence of P. expansum and other fungal pathogens with varying biology and lifestyle.

4.3. Glyoxalase (T-DNA Mutants T588 and T711, and ∆t-588 Mutant)

Glyoxalase I (GloI) is part of the glyoxalase system that catalyzes the detoxification of the endogenous reactive metabolite, methylglyoxal (MG), which is produced as a byproduct of glycolysis [38]. GloI is responsible for the isomerization of the hemithioacetal formed spontaneously from α-oxoaldehyde and GSH to S-2-hydroxyacylglutathione derivatives RCOCH(OH)-SG® RCH(OH)CO-SG in the cytosol of cells [39]. This reaction decreases the steady-state concentrations of the metabolic byproduct methylglyoxal and other reactive acyclic α-oxoaldehydes that mediate glycation reactions [40,41]. Glycation is a nonenzymatic glycosylation that impairs the functioning of biomolecules by bonding sugar molecules to a protein or lipid molecule. This process is implicated in many age-related chronic diseases, such as cardiovascular and Alzheimer’s disease [41,42,43,44]. The primary function of glyoxalase I is the detoxification of α-oxoaldehydes as part of the enzymatic defense against glycation. In addition, glyoxalases have been implicated in bacterial and fungal virulence. Glyoxalase A (gloA) was found to have a role in the virulence activation of Listeria monocytogenes, and similarly to ∆t-588, gloA mutants were sensitive to exogenous MG [45]. In C. albicans, the deletion of glx3 results in an avirulent phenotype along with decreased biofilm filamentation [46]. However, the role of glyoxalase in Penicillium–apple fruit interactions has not previously been ascertained. It is envisioned that the discovery of this new locus, encoding a glyoxalase enzyme, is important for the fungus to tolerate the generation of toxic byproducts from the central metabolism (e.g., methylglyoxal) that are generated during decay. Hence, this finding adds a new, logical perspective to the concept of fungal virulence factors, in that it is critical for the pathogen to detoxify metabolites from the primary metabolism to enable sufficient self-protection and to ensure optimal lesion/decay during the apple fruit colonization.

4.4. GPI1 (T-DNA Mutant T-434 and t-434KD)

Glycosylphosphatidylinositol (GPI) is a glycolipid that acts as an anchor to many cell-surface proteins, and in fungi, can also be found in their cell walls [47,48]. GPI anchors are bound to the C-terminus of proteins through post-translational modifications [48]. GPI-anchored proteins have diverse roles in cell biology and include cell–cell interactions, mammalian or protozoan antigens, and enzymes, among others [47]. Synthesis of GPI and its transfer to proteins comprise a series of 23 steps that take place in the endoplasmic reticulum (ER) [49]. In mammals, the first GPI synthesis step is catalyzed by a protein complex of six core subunits that include PIGQ or its homolog in Saccaromyces cerevisiae GPI1 [49]. In S. cerevisiae, gpi1-deficient strains show a temperature-sensitive phenotype, are unable to synthesize GPI in vitro, and, when cultured in permissive temperatures, have a separation defect after budding [50]. Due to the importance of GPI, a previous study proposed using GPI anchor biosynthesis as a potential target for developing antifungals [51]. This idea was explored through the development of host-induced gene silencing (HIGS) against many plant-pathogen fungal species and was evidenced by a patent [52]. While we could not generate a knockout strain, temperature was not considered as a selection criterion for screening transformants in this study. Hence, it cannot be ruled out that the possibility of obtaining a complete null mutant using temperature and other parameters that center around GPI biology and biochemistry. The assays within this study’s evaluation knockdown strain suggests that the level of GPI1 is sufficient to foster GPI biosynthesis, as no impact was observed in the growth nor virulence defect, in contrast to its T-DNA mutant counterpart. This is the first time a GPI-biosynthesis-related gene has been studied in Penicillium spp., and these findings shed light on additional aspects of protein post-translational modifications and promising antifungal targets as new areas for decay control.

4.5. HSP40 (T-DNA Mutant T-625 and ∆t-625)

Hsp40 proteins are divided into three classes, I, II, and III, which differ in their structural patterns [53,54]. The major role of Hsp40s are to regulate Hsp70’s activity through an ATP-dependent interaction [55,56,57]. Hsp70 and Hsp40 are involved in protein folding, transportation, and localization, and thus act in many cellular activities, including fugal–plant interactions [55,56,57,58]. For instance, Hsp70 in Colletotrichum gloeosporioides displayed higher expression in the infected coffee plant [59]. Hsp70 proteins of M. oryzae were also found to accumulate in the endoplasmic reticulum (ER) of fungi and were responsible for disease development in rice plants [60]. These proteins were observed during the cell-wall degradation of host plants, and thereby assisted virulence when a plant pathogenic fungus, Rhizoctonia solani, was attacking soybean seedlings [61].
Besides acting as an ATP-dependent chaperone, Hsp40 has also been found to independently assist viruses to invade host cells. In mammalian cells, double-stranded RNA-dependent protein kinase (PKR) is a key regulator of the antiviral innate immune response [62]. In P. expansum, Hsp40 exhibited lower expression (about 0.12-fold) in infected apple tissue than in culture at 10 dpi (Figure 3). This unexpected observation might be caused by the delay in examination regarding the transcriptional level of Hsp40. We only examined the expression of Hsp40 with highly decayed tissue at 10 dpi in this study. However, a previous study showed that the possible roles of this gene involve modulating protein secretion through the ER and vesicle-mediated transport [16].

4.6. Loci That Were Recalcitrant to Gene Deletion/Functional Analysis

Deletion mutants for the genes disrupted in T-DNA mutants T-191 and T-262 could not be obtained during this study. It can be hypothesized that these genes may be essential due to both their bioinformatic prediction and their recalcitrant behavior upon transformation with both deletion and knockdown constructs. Previous studies have also documented a lack of success in generating targeted deletion mutants after a random mutagenesis approach. For example, in another study, a random ATM was performed in B. cinerea, and the authors could not generate targeted deletion mutants for one out of three characterized genes [9]. In addition, two of these genes are involved in RNA processing, which is an essential cellular process; so, even though there may be redundancy, we hypothesize that removing any of these components can be lethal to the organism. One of the benefits of ATM is the possibility of studying essential genes due to the random insertion, which can provide a T-DNA mutant knockdown strain, as evidenced in this study.
While it is unfortunate that other loci did not phenocopy the T-DNA mutants and/or resulted in loci that are presumed lethal, this information provides abundant opportunities for further investigation. Mechanisms behind why the mutants did not phenocopy remain elusive and may be due to a variety of factors. These may include the transcriptional, translational, and post-translational impacts of the locus and/or the surrounding loci from the T-DNA insertion [9]. Perhaps applying comparative genomic and/or transcriptomic approaches to the strains to determine pleiotropic or off-target effects of the T-DNA insertion may be warranted. Previous studies in yeast, where the entire KO library was sequenced, showed compensatory mechanisms were activated upon deletion of a single locus [63]. Perhaps this may be a parallel situation in our system that warrants a more detailed omics-based investigation. Our approach to utilizing the genes with presumed lethal phenotypes can be further explored using dsRNA interference in single and combinations to not only assess gene function but also to determine if these loci can be targeted using dsRNAs as a form of control strategy which has been implemented in other cropping systems [64,65].

5. Conclusions

In general, this study provides new insights into the complex biological picture of P. expansum–apple fruit interactions and the genes that are involved in the decay process. As shown by the diversity of loci discovered from this study, the virulence and the decay process represent a complex, well-orchestrated system that has multiple overlapping systems involving secretion, membrane components, RNA processing, protein synthesis, and metabolite detoxification. For the first time, a random T-DNA mutant library was generated, and identified five loci important to virulence in this species. Although some of these loci seem to be essential, identification through this method was possible. Two deletion mutants and one knockdown strain were generated that were evaluated for virulence. Locus PEXP_003290 was found to phenocopy the T-DNA mutant and supports the hypothesis that this locus acts as a glyoxalase, which is important for full virulence during apple fruit decay. The generation of this random T-DNA transformant library also serves as a resource for investigating other aspects of P. expansum biology. Future studies should focus on obtaining a more comprehensive library by increasing the total number of transformants. It can then be utilized to investigate other aspects of P. expansum biology, like differential production of specialized metabolites, tolerance or sensitivity to various fungicides, the metabolism of substrates at cold temperatures, spore germination defects, protein secretion, and natural product synthesis. Any of these avenues that explore the functional genetic basis of a given biological process regulating the decay process are envisioned to set the foundation for future blue mold decay strategies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/jof9111066/s1, Table S1: List of primers for generating deletion and knock down mutants; Table S2: List of Strains used in this study; Figure S1: T-DNA of pPK2 hph; Figure S2: Southern Blot confirmation for deletion and knock-down strains and design schematics. Table S3: List of BLAST results for loci interrupted by T-DNA insertion.

Author Contributions

Conceptualization, D.L.-R., H.P. and W.M.J.II; Formal analysis, D.L.-R., H.P., V.L.G. and W.M.J.II; Funding acquisition, W.M.J.II; Investigation, D.L.-R., H.P., V.L.G. and W.M.J.II; Methodology, D.L.-R., H.P., V.L.G. and W.M.J.II; Project administration, D.L.-R., H.P. and W.M.J.II; Resources, N.P.K. and W.M.J.II; Supervision, J.M.F., N.P.K. and W.M.J.II; Validation, D.L.-R., H.P., V.L.G. and W.M.J.II; Visualization, D.L.-R., H.P., V.L.G. and W.M.J.II; Writing—original draft, D.L.-R. and H.P.; Writing—review and editing, D.L.-R., H.P., V.L.G., J.M.F., N.P.K. and W.M.J.II. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by United States Department of Agriculture (National Program 303 Plant Diseases).

Institutional Review Board Statement

Not Applicable.

Informed Consent Statement

Not Applicable.

Data Availability Statement

Data is contained within the article and supplementary material.

Acknowledgments

This research was supported entirely by base funds provided by the USDA-ARS National Program 303 Plant Diseases to Lead Scientist WMJ II, and in part by an appointment to the Agricultural Research Service (ARS) Research Participation Program administered by the Oak Ridge Institute for Science and Education (ORISE) through an interagency agreement between the U.S. Department of Energy (DOE) and the U.S. Department of Agriculture (USDA). ORISE is managed by ORAU under DOE contract number DE-SC0014664. All opinions expressed in this manuscript are the authors’ and do not necessarily reflect the policies and views of the USDA, DOE, or ORAU/ORISE.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Luciano-Rosario, D.; Keller, N.P.; Jurick, W.M. Penicillium expansum: Biology, Omics, and Management Tools for a Global Postharvest Pathogen Causing Blue Mould of Pome Fruit. Mol. Plant Pathol. 2020, 21, 1391–1404. [Google Scholar] [CrossRef]
  2. Ioi, J.D.; Zhou, T.; Tsao, R.; F Marcone, M. Mitigation of Patulin in Fresh and Processed Foods and Beverages. Toxins 2017, 9, 157. [Google Scholar] [CrossRef] [PubMed]
  3. Yu, L.; Qiao, N.; Zhao, J.; Zhang, H.; Tian, F.; Zhai, Q.; Chen, W. Postharvest Control of Penicillium expansum in Fruits: A Review. Food Biosci. 2020, 36, 100633. [Google Scholar] [CrossRef]
  4. Jurick, W.M.; Choi, M.-W.; Gaskins, V.L.; Peter, K.A.; Cox, K.D. Would You Like Wood or Plastic? Bin Material, Sanitation Treatments, and Bin Inoculum Levels Impact Blue Mold Decay of Stored Apple Fruit. Plant Dis. 2023, 107, 1177–1182. [Google Scholar] [CrossRef] [PubMed]
  5. Kumar, D.; Barad, S.; Chen, Y.; Luo, X.; Tannous, J.; Dubey, A.; Glam Matana, N.; Tian, S.; Li, B.; Keller, N.; et al. LaeA Regulation of Secondary Metabolism Modulates Virulence in Penicillium expansum and Is Mediated by Sucrose. Mol. Plant Pathol. 2017, 18, 1150–1163. [Google Scholar] [CrossRef]
  6. El Hajj Assaf, C.; Snini, S.P.; Tadrist, S.; Bailly, S.; Naylies, C.; Oswald, I.P.; Lorber, S.; Puel, O. Impact of veA on the Development, Aggressiveness, Dissemination and Secondary Metabolism of Penicillium Expansum. Mol. Plant Pathol. 2018, 19, 1971–1983. [Google Scholar] [CrossRef]
  7. Tannous, J.; Kumar, D.; Sela, N.; Sionov, E.; Prusky, D.; Keller, N.P. Fungal Attack and Host Defence Pathways Unveiled in Near-Avirulent Interactions of Penicillium expansum creA Mutants on Apples. Mol. Plant Pathol. 2018, 19, 2635–2650. [Google Scholar] [CrossRef]
  8. Sánchez-Torres, P.; Vilanova, L.; Ballester, A.R.; López-Pérez, M.; Teixidó, N.; Viñas, I.; Usall, J.; González-Candelas, L.; Torres, R. Unravelling the Contribution of the Penicillium expansum PeSte12 Transcription Factor to Virulence during Apple Fruit Infection. Food Microbiol. 2018, 69, 123–135. [Google Scholar] [CrossRef]
  9. Giesbert, S.; Schumacher, J.; Kupas, V.; Espino, J.; Segmüller, N.; Haeuser-Hahn, I.; Schreier, P.H.; Tudzynski, P. Identification of Pathogenesis-Associated Genes by T-DNA–Mediated Insertional Mutagenesis in Botrytis Cinerea: A Type 2A Phosphoprotein Phosphatase and an SPT3 Transcription Factor Have Significant Impact on Virulence. Mol. Plant-Microbe Interact. 2012, 25, 481–495. [Google Scholar] [CrossRef]
  10. Michielse, C.B.; van Wijk, R.; Reijnen, L.; Cornelissen, B.J.; Rep, M. Insight into the Molecular Requirements for Pathogenicity of Fusarium Oxysporum f. Sp. Lycopersici through Large-Scale Insertional Mutagenesis. Genome Biol. 2009, 10, R4. [Google Scholar] [CrossRef]
  11. Yu, M.; Yu, J.; Hu, J.; Huang, L.; Wang, Y.; Yin, X.; Nie, Y.; Meng, X.; Wang, W.; Liu, Y. Identification of Pathogenicity-Related Genes in the Rice Pathogen Ustilaginoidea Virens through Random Insertional Mutagenesis. Fungal Genet. Biol. 2015, 76, 10–19. [Google Scholar] [CrossRef] [PubMed]
  12. Xu, Y.; Ao, K.; Tian, L.; Qiu, Y.; Huang, X.; Liu, X.; Hoy, R.; Zhang, Y.; Rashid, K.Y.; Xia, S.; et al. A Forward Genetic Screen in Sclerotinia sclerotiorum Revealed the Transcriptional Regulation of Its Sclerotial Melanization Pathway. Mol. Plant Microbe Interact. 2022, 35, 244–256. [Google Scholar] [CrossRef] [PubMed]
  13. Deng, S.; Wang, C.; Zhang, X.; Lin, L. Bidirectional Promoter Trapping T-DNA for Insertional Mutagenesis in Verticillium Dahliae. Can. J. Microbiol. 2014, 60, 445–454. [Google Scholar] [CrossRef] [PubMed]
  14. Maruthachalam, K.; Klosterman, S.J.; Kang, S.; Hayes, R.J.; Subbarao, K.V. Identification of Pathogenicity-Related Genes in the Vascular Wilt Fungus Verticillium dahliae by Agrobacterium tumefaciens-Mediated T-DNA Insertional Mutagenesis. Mol. Biotechnol. 2011, 49, 209–221. [Google Scholar] [CrossRef] [PubMed]
  15. Lv, B.; Zheng, L.; Liu, H.; Tang, J.; Hsiang, T.; Huang, J. Use of Random T-DNA Mutagenesis in Identification of Gene UvPRO1, A Regulator of Conidiation, Stress Response, and Virulence in Ustilaginoidea Virens. Front. Microbiol. 2016, 7, 2086. [Google Scholar] [CrossRef]
  16. Jurick, W.M.; Peng, H.; Beard, H.S.; Garrett, W.M.; Lichtner, F.J.; Luciano-Rosario, D.; Macarisin, O.; Liu, Y.; Peter, K.A.; Gaskins, V.L.; et al. Blistering1 Modulates Penicillium expansum Virulence Via Vesicle-Mediated Protein Secretion. Mol. Cell. Proteom. 2020, 19, 344–361. [Google Scholar] [CrossRef]
  17. Utermark, J.; Karlovsky, P. Genetic Transformation of Filamentous Fungi by Agrobacterium Tumefaciens. Protoc. Exch. 2008. [Google Scholar] [CrossRef]
  18. Mullins, E.D.; Chen, X.; Romaine, P.; Raina, R.; Geiser, D.M.; Kang, S. Agrobacterium-Mediated Transformation of Fusarium oxysporum: An Efficient Tool for Insertional Mutagenesis and Gene Transfer. Phytopathology 2001, 91, 173–180. [Google Scholar] [CrossRef]
  19. Höfgen, R.; Willmitzer, L. Storage of Competent Cells for Agrobacterium Transformation. Nucleic Acids Res. 1988, 16, 9877. [Google Scholar] [CrossRef]
  20. Singer, T.; Burke, E. High-Throughput TAIL-PCR as a Tool to Identify DNA Flanking Insertions. Methods Mol. Biol. 2003, 236, 241–272. [Google Scholar] [CrossRef]
  21. Livak, K.J.; Schmittgen, T.D. Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods 2001, 25, 402–408. [Google Scholar] [CrossRef] [PubMed]
  22. Lim, F.Y.; Sanchez, J.F.; Wang, C.C.C.; Keller, N.P. Toward Awakening Cryptic Secondary Metabolite Gene Clusters in Filamentous Fungi. Methods Enzymol. 2012, 517, 303–324. [Google Scholar] [CrossRef] [PubMed]
  23. Nickles, G.; Ludwikoski, I.; Bok, J.W.; Keller, N.P. Comprehensive Guide to Extracting and Expressing Fungal Secondary Metabolites with Aspergillus fumigatus as a Case Study. Curr. Protoc. 2021, 1, e321. [Google Scholar] [CrossRef] [PubMed]
  24. Greco, C.; Pfannenstiel, B.T.; Liu, J.C.; Keller, N.P. Depsipeptide Aspergillicins Revealed by Chromatin Reader Protein Deletion. ACS Chem. Biol. 2019, 14, 1121–1128. [Google Scholar] [CrossRef]
  25. Thornalley, P.J. Glyoxalase I—Structure, Function and a Critical Role in the Enzymatic Defence against Glycation. Biochem. Soc. Trans. 2003, 31, 1343–1348. [Google Scholar] [CrossRef]
  26. Perales-Calvo, J.; Muga, A.; Moro, F. Role of DnaJ G/F-Rich Domain in Conformational Recognition and Binding of Protein Substrates. J. Biol. Chem. 2010, 285, 34231–34239. [Google Scholar] [CrossRef]
  27. Long, X.; He, N.-M.; Tan, L.-X.; Yang, Y.-H.; Zhou, J.-P.; Liu, Z.-Y.; Mo, M.-H.; Liu, T. Methylglyoxal Has Different Impacts on the Fungistatic Roles of Ammonia and Benzaldehyde, and Lactoylglutathione Lyase Is Necessary for the Resistance of Arthrobotrys oligospora to Soil Fungistasis. Front. Cell. Infect. Microbiol. 2021, 11, 640823. [Google Scholar] [CrossRef]
  28. Li, G.; Zhou, Z.; Liu, G.; Zheng, F.; He, C. Characterization of T-DNA Insertion Patterns in the Genome of Rice Blast Fungus Magnaporthe Oryzae. Curr. Genet. 2007, 51, 233–243. [Google Scholar] [CrossRef]
  29. Sato, M.; Araki, N.; Kumeta, M.; Takeyasu, K.; Taguchi, Y.; Asai, T.; Furukawa, K.; Horigome, T. Interaction, Mobility, and Phosphorylation of Human Orthologues of WD Repeat-Containing Components of the Yeast SSU Processome t-UTP Sub-Complex. Biochem. Cell Biol. 2013, 91, 466–475. [Google Scholar] [CrossRef]
  30. Mouillesseaux, K.; Chen, J.-N. Mutation in Utp15 Disrupts Vascular Patterning in a P53-Dependent Manner in Zebrafish Embryos. PLoS ONE 2011, 6, e25013. [Google Scholar] [CrossRef]
  31. Xu, G.; Yang, S.; Meng, L.; Wang, B.-G. The Plant Hormone Abscisic Acid Regulates the Growth and Metabolism of Endophytic Fungus Aspergillus nidulans. Sci. Rep. 2018, 8, 6504. [Google Scholar] [CrossRef] [PubMed]
  32. Nie, A.; Sun, B.; Fu, Z.; Yu, D. Roles of Aminoacyl-tRNA Synthetases in Immune Regulation and Immune Diseases. Cell Death Dis. 2019, 10, 901. [Google Scholar] [CrossRef] [PubMed]
  33. Sang Lee, J.; Gyu Park, S.; Park, H.; Seol, W.; Lee, S.; Kim, S. Interaction Network of Human Aminoacyl-tRNA Synthetases and Subunits of Elongation Factor 1 Complex. Biochem. Biophys. Res. Commun. 2002, 291, 158–164. [Google Scholar] [CrossRef]
  34. Rubio Gomez, M.A.; Ibba, M. Aminoacyl-tRNA Synthetases. RNA 2020, 26, 910–936. [Google Scholar] [CrossRef] [PubMed]
  35. Naushad, H.S.; Gupta, R.S. Phylogenomics and Molecular Signatures for Species from the Plant Pathogen-Containing Order Xanthomonadales. PLoS ONE 2013, 8, e55216. [Google Scholar] [CrossRef]
  36. Szymański, M.; Deniziak, M.; Barciszewski, J. The New Aspects of Aminoacyl-tRNA Synthetases. Acta Biochim. Pol. 2000, 47, 821–834. [Google Scholar] [CrossRef]
  37. Pang, L.; Weeks, S.D.; Van Aerschot, A. Aminoacyl-tRNA Synthetases as Valuable Targets for Antimicrobial Drug Discovery. Int. J. Mol. Sci. 2021, 22, 1750. [Google Scholar] [CrossRef]
  38. Morgenstern, J.; Campos Campos, M.; Nawroth, P.; Fleming, T. The Glyoxalase System—New Insights into an Ancient Metabolism. Antioxidants 2020, 9, 939. [Google Scholar] [CrossRef]
  39. Thornalley, P.J. The Glyoxalase System in Health and Disease. Mol. Asp. Med. 1993, 14, 287–371. [Google Scholar] [CrossRef]
  40. Rabbani, N.; Thornalley, P.J. Methylglyoxal, Glyoxalase 1 and the Dicarbonyl Proteome. Amino Acids 2012, 42, 1133–1142. [Google Scholar] [CrossRef]
  41. Rabbani, N.; Xue, M.; Weickert, M.O.; Thornalley, P.J. Multiple Roles of Glyoxalase 1-Mediated Suppression of Methylglyoxal Glycation in Cancer Biology-Involvement in Tumour Suppression, Tumour Growth, Multidrug Resistance and Target for Chemotherapy. Semin. Cancer Biol. 2018, 49, 83–93. [Google Scholar] [CrossRef] [PubMed]
  42. Thornalley, P.J.; Battah, S.; Ahmed, N.; Karachalias, N.; Agalou, S.; Babaei-Jadidi, R.; Dawnay, A. Quantitative Screening of Advanced Glycation Endproducts in Cellular and Extracellular Proteins by Tandem Mass Spectrometry. Biochem. J. 2003, 375, 581–592. [Google Scholar] [CrossRef] [PubMed]
  43. Münch, G.; Deuther-Conrad, W.; Gasic-Milenkovic, J. Glycoxidative Stress Creates a Vicious Cycle of Neurodegeneration in Alzheimer’s Disease—A Target for Neuroprotective Treatment Strategies? J. Neural Transm. Suppl. 2002, 62, 303–307. [Google Scholar] [CrossRef]
  44. Münch, G.; Mayer, S.; Michaelis, J.; Hipkiss, A.R.; Riederer, P.; Müller, R.; Neumann, A.; Schinzel, R.; Cunningham, A.M. Influence of Advanced Glycation End-Products and AGE-Inhibitors on Nucleation-Dependent Polymerization of Beta-Amyloid Peptide. Biochim. Biophys. Acta 1997, 1360, 17–29. [Google Scholar] [CrossRef]
  45. Anaya-Sanchez, A.; Feng, Y.; Berude, J.C.; Portnoy, D.A. Detoxification of Methylglyoxal by the Glyoxalase System Is Required for Glutathione Availability and Virulence Activation in Listeria Monocytogenes. PLoS Pathog. 2021, 17, e1009819. [Google Scholar] [CrossRef]
  46. Cabello, L.; Gómez-Herreros, E.; Fernández-Pereira, J.; Maicas, S.; Martínez-Esparza, M.C.; de Groot, P.W.J.; Valentín, E. Deletion of GLX3 in Candida albicans Affects Temperature Tolerance, Biofilm Formation and Virulence. FEMS Yeast Res. 2019, 19, foy124. [Google Scholar] [CrossRef]
  47. Paulick, M.G.; Bertozzi, C.R. The Glycosylphosphatidylinositol Anchor: A Complex Membrane-Anchoring Structure for Proteins. Biochemistry 2008, 47, 6991–7000. [Google Scholar] [CrossRef]
  48. Kinoshita, T. Glycosylphosphatidylinositol (GPI) Anchors: Biochemistry and Cell Biology: Introduction to a Thematic Review Series. J. Lipid Res. 2016, 57, 4–5. [Google Scholar] [CrossRef]
  49. Liu, S.-S.; Liu, Y.-S.; Guo, X.-Y.; Murakami, Y.; Yang, G.; Gao, X.-D.; Kinoshita, T.; Fujita, M. A Knockout Cell Library of GPI Biosynthetic Genes for Functional Studies of GPI-Anchored Proteins. Commun. Biol. 2021, 4, 1–11. [Google Scholar] [CrossRef]
  50. Leidich, S.D.; Orlean, P. Gpi1, a Saccharomyces Cerevisiae Protein That Participates in the First Step in Glycosylphosphatidylinositol Anchor Synthesis. J. Biol. Chem. 1996, 271, 27829–27837. [Google Scholar] [CrossRef]
  51. McLellan, C.A.; Whitesell, L.; King, O.D.; Lancaster, A.K.; Mazitschek, R.; Lindquist, S. Inhibiting GPI Anchor Biosynthesis in Fungi Stresses the Endoplasmic Reticulum and Enhances Immunogenicity. ACS Chem. Biol. 2012, 7, 1520–1528. [Google Scholar] [CrossRef] [PubMed]
  52. Stahl, D.; Temme, N.; Deising, H.; Garcia Oliveira, E.; Schafer, W. Development of Fungal Resistant Crops by HIGS (Host-Induced Gene Silencing) Mediated Inhibition of Gpi-Anchored Cell Wall Protein Synthesis. U.S. Patent Application No. 16/317,335, 7 July 2017. [Google Scholar]
  53. Cheetham, M.E.; Caplan, A.J. Structure, Function and Evolution of DnaJ: Conservation and Adaptation of Chaperone Function. Cell Stress Chaperones 1998, 3, 28–36. [Google Scholar] [CrossRef] [PubMed]
  54. Hennessy, F.; Cheetham, M.E.; Dirr, H.W.; Blatch, G.L. Analysis of the Levels of Conservation of the J Domain among the Various Types of DnaJ-like Proteins. Cell Stress Chaperones 2000, 5, 347–358. [Google Scholar] [CrossRef] [PubMed]
  55. Shen, Y.; Meunier, L.; Hendershot, L.M. Identification and Characterization of a Novel Endoplasmic Reticulum (ER) DnaJ Homologue, Which Stimulates ATPase Activity of BiP in Vitro and Is Induced by ER Stress. J. Biol. Chem. 2002, 277, 15947–15956. [Google Scholar] [CrossRef]
  56. Craig, E.A.; Huang, P.; Aron, R.; Andrew, A. The Diverse Roles of J-Proteins, the Obligate Hsp70 Co-Chaperone. Rev. Physiol. Biochem. Pharmacol. 2006, 156, 1–21. [Google Scholar] [CrossRef]
  57. Steel, R.; Doherty, J.P.; Buzzard, K.; Clemons, N.; Hawkins, C.J.; Anderson, R.L. Hsp72 Inhibits Apoptosis Upstream of the Mitochondria and Not through Interactions with Apaf-1. J. Biol. Chem. 2004, 279, 51490–51499. [Google Scholar] [CrossRef]
  58. Kampinga, H.H.; Craig, E.A. The HSP70 Chaperone Machinery: J Proteins as Drivers of Functional Specificity. Nat. Rev. Mol. Cell Biol. 2010, 11, 579–592. [Google Scholar] [CrossRef]
  59. Chen, Z.J.; Ribeiro, A.; Silva, M.C.; Santos, P.; Guerra-Guimarães, L.; Gouveia, M.; Fernandez, D.; Rodrigues, C.J. Heat Shock-Induced Susceptibility of Green Coffee Leaves and Berries to Colletotrichum gloeosporioides and Its Association to PR and Hsp70 Gene Expression. Physiol. Mol. Plant Pathol. 2003, 63, 181–190. [Google Scholar] [CrossRef]
  60. Yi, M.; Chi, M.-H.; Khang, C.H.; Park, S.-Y.; Kang, S.; Valent, B.; Lee, Y.-H. The ER Chaperone LHS1 Is Involved in Asexual Development and Rice Infection by the Blast Fungus Magnaporthe oryzae. Plant Cell 2009, 21, 681–695. [Google Scholar] [CrossRef]
  61. Lakshman, D.K.; Roberts, D.P.; Garrett, W.M.; Natarajan, S.S.; Darwish, O.; Alkharouf, N.; Pain, A.; Khan, F.; Jambhulkar, P.P.; Mitra, A. Proteomic Investigation of Rhizoctonia solani AG 4 Identifies Secretome and Mycelial Proteins with Roles in Plant Cell Wall Degradation and Virulence. J. Agric. Food Chem. 2016, 64, 3101–3110. [Google Scholar] [CrossRef]
  62. Connor, J.H.; Lyles, D.S. Inhibition of Host and Viral Translation during Vesicular Stomatitis Virus Infection. eIF2 Is Responsible for the Inhibition of Viral but Not Host Translation. J. Biol. Chem. 2005, 280, 13512–13519. [Google Scholar] [CrossRef] [PubMed]
  63. Puddu, F.; Herzog, M.; Selivanova, A.; Wang, S.; Zhu, J.; Klein-Lavi, S.; Gordon, M.; Meirman, R.; Millan-Zambrano, G.; Ayestaran, I.; et al. Genome Architecture and Stability in the Saccharomyces Cerevisiae Knockout Collection. Nature 2019, 573, 416–420. [Google Scholar] [CrossRef] [PubMed]
  64. Wytinck, N.; Ziegler, D.J.; Walker, P.L.; Sullivan, D.S.; Biggar, K.T.; Khan, D.; Sakariyahu, S.K.; Wilkins, O.; Whyard, S.; Belmonte, M.F. Host Induced Gene Silencing of the Sclerotinia sclerotiorum ABHYDROLASE-3 Gene Reduces Disease Severity in Brassica Napus. PLoS ONE 2022, 17, e0261102. [Google Scholar] [CrossRef] [PubMed]
  65. Degnan, R.M.; McTaggart, A.R.; Shuey, L.S.; Pame, L.J.S.; Smith, G.R.; Gardiner, D.M.; Nock, V.; Soffe, R.; Sale, S.; Garrill, A.; et al. Exogenous Double-Stranded RNA Inhibits the Infection Physiology of Rust Fungi to Reduce Symptoms in Planta. Mol. Plant Pathol. 2023, 24, 191–207. [Google Scholar] [CrossRef]
Figure 1. Morphology and virulence of T-DNA transformants. (A) Growth of T-DNA transformants on apples and PDA medium. Golden Delicious apples were inoculated with R19 or transformant spores and then stored in the dark at room temperature for 10 days. PDA plates were added R19 or transformant spores and then stored in an incubator at 25 °C for 6.5 days. (B) Virulence and growth reduction to wild-type strain R19 in transformants. The data represent the mean value (±SD) of three replicates, and comparisons were assessed using unpaired t-test with the Welch correction.
Figure 1. Morphology and virulence of T-DNA transformants. (A) Growth of T-DNA transformants on apples and PDA medium. Golden Delicious apples were inoculated with R19 or transformant spores and then stored in the dark at room temperature for 10 days. PDA plates were added R19 or transformant spores and then stored in an incubator at 25 °C for 6.5 days. (B) Virulence and growth reduction to wild-type strain R19 in transformants. The data represent the mean value (±SD) of three replicates, and comparisons were assessed using unpaired t-test with the Welch correction.
Jof 09 01066 g001
Figure 2. DNA gel-blotting analysis of T-DNA. (A) WT (R19) genomic DNA was digested with restriction endonucleases PAC I. Molecular markers are indicated on the left. (B) Diagram of P. expansum genes interrupted by T-DNA insertion. A total of five transformants were used for localizing T-DNA insertion through thermal asymmetric interlaced polymerase chain reaction (TAIL-PCR). Five genes were found potentially irrupted by T-DNA insertion in coding sequences (CDS) or untranslated regions (UTR), including GloI (glyoxalases I), GlyRS (glycyl-tRNA synthetase), HSP40 (heat-shock protein 40), GPI1 (N-acetylglucosaminyl transferases), and RAP15 (U8 small nucleolar RNA-associated protein 15). The gene structure figure was generated with Gene Structure Display Server (v2.0; http://gsds.cbi.pku.edu.cn/index.php accessed on 2 February 2019).
Figure 2. DNA gel-blotting analysis of T-DNA. (A) WT (R19) genomic DNA was digested with restriction endonucleases PAC I. Molecular markers are indicated on the left. (B) Diagram of P. expansum genes interrupted by T-DNA insertion. A total of five transformants were used for localizing T-DNA insertion through thermal asymmetric interlaced polymerase chain reaction (TAIL-PCR). Five genes were found potentially irrupted by T-DNA insertion in coding sequences (CDS) or untranslated regions (UTR), including GloI (glyoxalases I), GlyRS (glycyl-tRNA synthetase), HSP40 (heat-shock protein 40), GPI1 (N-acetylglucosaminyl transferases), and RAP15 (U8 small nucleolar RNA-associated protein 15). The gene structure figure was generated with Gene Structure Display Server (v2.0; http://gsds.cbi.pku.edu.cn/index.php accessed on 2 February 2019).
Jof 09 01066 g002
Figure 3. Expression patterns of genes potentially responsible for the virulence mutation of transformants in PDA culture and decayed apple tissues. Total RNA was extracted from PDA culture and decayed apple tissues after 10 days post inoculation (DPI). The data represent the mean value (±SD) of three replicates.
Figure 3. Expression patterns of genes potentially responsible for the virulence mutation of transformants in PDA culture and decayed apple tissues. Total RNA was extracted from PDA culture and decayed apple tissues after 10 days post inoculation (DPI). The data represent the mean value (±SD) of three replicates.
Jof 09 01066 g003
Figure 4. Targeted mutant morphology and physiology. (A) Mutant appearance after cultivation in GMM plates for 7 days using initial inoculum of 106 spores at 25 °C. Knockdown mutant was assessed in GMM and GMM-amended media (using NH4 as a nitrogen source). (B) Radial growth measurements 7 days post inoculation in GMM, and germination initiation assay on liquid minimal media (LMM) 21 h post inoculation at 25 °C. * p < 0.05, ** p < 0.01, *** p < 0.001.
Figure 4. Targeted mutant morphology and physiology. (A) Mutant appearance after cultivation in GMM plates for 7 days using initial inoculum of 106 spores at 25 °C. Knockdown mutant was assessed in GMM and GMM-amended media (using NH4 as a nitrogen source). (B) Radial growth measurements 7 days post inoculation in GMM, and germination initiation assay on liquid minimal media (LMM) 21 h post inoculation at 25 °C. * p < 0.05, ** p < 0.01, *** p < 0.001.
Jof 09 01066 g004
Figure 5. Apple inoculation assays for targeted mutant strains. Apple lesions were measured for 7 days in Fuji, Golden Delicious, and Honeycrisp cultivars. ** p < 0.01, *** p < 0.001.
Figure 5. Apple inoculation assays for targeted mutant strains. Apple lesions were measured for 7 days in Fuji, Golden Delicious, and Honeycrisp cultivars. ** p < 0.01, *** p < 0.001.
Jof 09 01066 g005
Figure 6. Congo Red and sorbitol susceptibility plate assay. The mutants were assessed at different concentrations of GMM amended with (A) Congo Red or (B) sorbitol. GMM using ammonium salts was used to evaluate t-434KD. Images were obtained 5 days post inoculation.
Figure 6. Congo Red and sorbitol susceptibility plate assay. The mutants were assessed at different concentrations of GMM amended with (A) Congo Red or (B) sorbitol. GMM using ammonium salts was used to evaluate t-434KD. Images were obtained 5 days post inoculation.
Jof 09 01066 g006
Figure 7. Methylglyoxal (MG) plate assay. Deletion strain Δt-588 was assessed for in vitro susceptibility of 0.1% MG in GMM-amended plates using an initial inoculum of 100 spores. (A) Plates were imaged and the (B) colony diameters quantified 12 days post inoculation.
Figure 7. Methylglyoxal (MG) plate assay. Deletion strain Δt-588 was assessed for in vitro susceptibility of 0.1% MG in GMM-amended plates using an initial inoculum of 100 spores. (A) Plates were imaged and the (B) colony diameters quantified 12 days post inoculation.
Jof 09 01066 g007
Table 1. T-DNA mutant list and annotation description.
Table 1. T-DNA mutant list and annotation description.
T-DNA Mutant NameLocus IDAnnotation Description
T-193PeRAP15 (PEXP_001700)U3 small nucleolar RNA-associated protein 15, C-terminal
T-275PeGlyRs (PEXP_016530)hypothetical protein
T-588/T-711PeGIO (PEXP_003290)hypothetical protein/VOC containing domain/serine rich protein
T-434PeGPI1 (PEXP_051540)phosphatidylinositol N-acetylglucosaminyltransferase subunit Q/GPI1
T-625Blistering1 (PEXP_008080)DNA-J-domain-containing protein
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Luciano-Rosario, D.; Peng, H.; Gaskins, V.L.; Fonseca, J.M.; Keller, N.P.; Jurick, W.M., II. Mining the Penicillium expansum Genome for Virulence Genes: A Functional-Based Approach to Discover Novel Loci Mediating Blue Mold Decay of Apple Fruit. J. Fungi 2023, 9, 1066. https://doi.org/10.3390/jof9111066

AMA Style

Luciano-Rosario D, Peng H, Gaskins VL, Fonseca JM, Keller NP, Jurick WM II. Mining the Penicillium expansum Genome for Virulence Genes: A Functional-Based Approach to Discover Novel Loci Mediating Blue Mold Decay of Apple Fruit. Journal of Fungi. 2023; 9(11):1066. https://doi.org/10.3390/jof9111066

Chicago/Turabian Style

Luciano-Rosario, Dianiris, Hui Peng, Verneta L. Gaskins, Jorge M. Fonseca, Nancy P. Keller, and Wayne M. Jurick, II. 2023. "Mining the Penicillium expansum Genome for Virulence Genes: A Functional-Based Approach to Discover Novel Loci Mediating Blue Mold Decay of Apple Fruit" Journal of Fungi 9, no. 11: 1066. https://doi.org/10.3390/jof9111066

APA Style

Luciano-Rosario, D., Peng, H., Gaskins, V. L., Fonseca, J. M., Keller, N. P., & Jurick, W. M., II. (2023). Mining the Penicillium expansum Genome for Virulence Genes: A Functional-Based Approach to Discover Novel Loci Mediating Blue Mold Decay of Apple Fruit. Journal of Fungi, 9(11), 1066. https://doi.org/10.3390/jof9111066

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop