Next Article in Journal
Effects of Zinc Acetate Hydrate Supplementation on Renal Anemia with Hypozincemia in Hemodialysis Patients
Next Article in Special Issue
Occurrence and Exposure Assessment of Major Mycotoxins in Foodstuffs from Algeria
Previous Article in Journal
Stability and Chemical Conversion of the Purified Reference Material of Gymnodimine-A under Different Temperature and pH Conditions
Previous Article in Special Issue
Development and Validation of LC-Q-TOF-MS Methodology to Determine Mycotoxin Biomarkers in Human Urine
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Ochratoxin A Defective Aspergillus carbonarius Mutants as Potential Biocontrol Agents

by
Belén Llobregat
,
Luis González-Candelas
and
Ana-Rosa Ballester
*
Instituto de Agroquímica y Tecnología de Alimentos (IATA), CSIC, Calle Catedrático Agustín Escardino 7, 46980 Paterna, Spain
*
Author to whom correspondence should be addressed.
Toxins 2022, 14(11), 745; https://doi.org/10.3390/toxins14110745
Submission received: 3 October 2022 / Revised: 27 October 2022 / Accepted: 28 October 2022 / Published: 31 October 2022
(This article belongs to the Special Issue Occurrence and Determination of Mycotoxins)

Abstract

:
Aspergillus carbonarius is one of the main species responsible for wine, coffee and cocoa toxin contamination. The main mycotoxin produced by this fungus, ochratoxin A (OTA), is a secondary metabolite categorized as a possible carcinogen because of its significant nephrotoxicity and immunosuppressive effects. A polyketide synthase gene (otaA) encodes the first enzyme in the OTA biosynthetic pathway. It is known that the filamentous fungi, growth, development and production of secondary metabolites are interconnected processes governed by global regulatory factors whose encoding genes are generally located outside the gene clusters involved in the biosynthesis of each secondary metabolite, such as the veA gene, which forms part of the VELVET complex. Different fungal strains compete for nutrients and space when they infect their hosts, and safer non-mycotoxigenic strains may be able to outcompete mycotoxigenic strains during colonization. To determine the possible utility of biopesticides based on the competitive exclusion of mycotoxigenic strains by non-toxigenic ones, we used A. carbonarius ΔotaA and ΔveA knockout mutants. Our results showed that during both in vitro growth and infection of grapes, non-mycotoxigenic strains could outcompete the wild-type strain. Additionally, the introduction of the non-mycotoxigenic strain led to a drastic decrease in OTA during both in vitro growth and infection of grapes.
Key Contribution: Non-mycotoxigenic strains may outcompete mycotoxigenic Aspergillus carbonarius strains and lower the ochratoxin A level in infected grapes.

1. Introduction

Mycotoxins are highly poisonous and carcinogenic secondary metabolites produced by several fungi. Aspergillus, Fusarium and Penicillium are the major filamentous fungi associated with the production of the mycotoxins that are hazardous for both human and animal health [1]. Aflatoxins, deoxynivalenol, ochratoxin A (OTA), fumonisins, T-2 and HT-2 toxins, patulin and zearalenone are considered the most relevant mycotoxins from a food/feed safety point of view [2]. The best strategy to lower mycotoxin levels in the food chain is to prevent the growth of the fungus and, if it has already grown, to prevent it from producing the toxin. Different physical (temperature and humidity management), chemical (use of chemical antifungal agents, electrolyzed oxidizing water treatment, etc.), and biological (use of antagonistic microorganisms, plant extracts) methods have been proposed [3]. One of the most important strategies for managing fungal diseases is to employ antagonistic microorganisms. The four main mechanisms deployed by biocontrol microorganisms are: (i) competing for resources and available space; (ii) producing antibiotics; (iii) inducing resistance; and (iv) direct parasitism. In addition to bacteria and yeast, it is also feasible to use filamentous fungi as biocontrol agents, and fungi that do not produce toxins and can compete with mycotoxigenic fungi are especially interesting. The mechanism of mycotoxin biocontrol is primarily attributed to competitive exclusion, but mycotoxins are sometimes reduced by a larger amount, which can be explained by the displacement of mycotoxin-producing fungi on the crop. Some examples are the uses of aflatoxigenic strains of Aspergillus flavus. Different products are commercially available: A. flavus AF36 (AF36 Prevail ®) was the first atoxigenic A. flavus strain-based product to be registered worldwide for the biological control of aflatoxins in fields in 2003 [4]; A. flavus strain NRL18543 is the active ingredient of this product. In 2004, the U.S. Environmental Protection Agency also registered Afla-Guard ®, which contains A. flavus NRRL21882 as the active agent [5], and it completely lacks both aflatoxin and cyclopiazonic acid [6]. A. flavus Og0222, the biocontrol ingredient of AflasafeTM, was developed for African countries and is commercially available for applications in maize and groundnuts, among others. Atoxigenic A. flavus strain MUCL5491 (biocontrol product named AF-X1) has also been approved for use in Italy to control the aflatoxin contamination of maize [7]. These atoxigenic A. flavus-based biocontrol products have been developed worldwide and share the characteristic of using native A. flavus strains that have been selected for the region where they will be employed. These particular strains have been carefully chosen because they do not produce aflatoxins and/or cyclopiazonic acid [6,8]. However, whether they produce other mycotoxins or harmful secondary metabolites is unknown. In this regard, we should consider the impact that global climatic changes can have on environmental factors, such as water, temperature and carbon dioxide, which could result in the production of new mycotoxins and/or higher levels of already known mycotoxins and, therefore, food safety and security concerns could increase [9,10].
OTA is a toxic secondary metabolite produced by the fungi of the Aspergillus and Penicillium genera, among others. Exposure to a low concentration of this mycotoxin may lead to neurotoxicity, hepatotoxicity, immunotoxicity and carcinogenicity. Therefore, OTA has been defined as Group 2B (possible carcinogenic to humans) by the International Agency for Research on Cancer (IARC). Comparative genomics has helped us to identify the genes involved in OTA biosynthesis [11,12,13,14]. The role of different genes has been proved: otaA, previously designated as AcOTApks [15], codes for a polyketide synthase [15,16]; otaB-AcOTAnrps codes for a non-ribosomal peptide synthase [17]; otaD-AcOTAhal encodes a halogenase [18]. A P450 cytochrome oxidase (otaC-AcOTAP450) [12], a bZIP transcription factor (otaR1-AcOTAbZIP) [19], and a cyclase otaY-AcOTASnoaL gene [20] also form part of the OTA biosynthetic cluster. In addition to the specific genes of the OTA biosynthetic pathway, different studies have shown the involvement of global transcription factors in the regulation of secondary metabolites biosynthesis in A. carbonarius, including genes veA and laeA [16,21]. Disruption of either of these two genes significantly reduces OTA production during in vitro and in vivo growth.
Recent works have established that climate change can modify the OTA levels produced by A. carbonarius [9,10,22]. Thus, ensuring that the selected strain does not produce mycotoxins regardless of environmental changes in temperature or humidity is crucial. The working hypothesis of this study is that the application of non-ochratoxigenic strains might lower mycotoxin levels in food and feed. In this work, we focused on two genes for the generation of non-OTA-producing A. carbonarius strains: otaA and veA. While the deletion of the otaA gene should affect only OTA production, the deletion of veA could alter the production of other secondary metabolites [23,24] to pave the way to develop safer non-mycotoxigenic strains. We investigated the putative biocontrol activity of these two mutants and determined the lowering of mycotoxin levels as a consequence of the displacement of the wild-type (wt) strain during the co-cultivation of wt and mutant strains in culture medium and during grape berry infection.

2. Results

2.1. Construction and Characterization of Knockout Mutants

In order to study whether a non-mycotoxigenic strain could outcompete a mycotoxigenic strain, we constructed two different knockout mutants in the OTA-producing A. carbonarius ITEM 5010 background [17]: (i) one based on the deletion of the otaA gene (the mutant denoted as Δpks), the first gene in the OTA biosynthesis; (ii) a second one based on the veA gene (the mutant denoted as ΔveA), and a global regulator of secondary metabolism. The followed gene replacement strategy is shown in Figure S1. Gene replacement plasmids pRFHU2-pks and pRFHU2-veA were obtained by the USER-Friendly cloning strategy [25]. The co-cultivation of the Agrobacterium tumefaciens cells carrying the plasmid with the conidia of A. carbonarius led to hygromycin-resistant colonies appearing approximately 3 days after being transferred to selective PDA plates. PCR amplification was used to check the correct insertion of the T-DNA containing the hygromycin-resistant marker that replaced the gene of interest (Figure S1). To select the knockout mutants with no further T-DNA integrations, the copy number of the integrated T-DNA was determined by quantitative real-time PCR (qPCR) (Table S1). Three knockout mutants and two ectopic mutants for each construction were selected for further analyses.
The phenotypic traits of the wt strain, and the ectopic and knockout mutants regarding growth, conidiation and OTA production, were recorded on PDA plates after 5 days post inoculation (dpi) of growth (Figure 1). The secondary metabolism was impaired in the ΔveA mutants. It resulted in a brownish pigment observed at the back of the colony compared to the wt strain and the Δpks mutants (Figure 1A). Non-statistically significant differences for the growth area were detected among the wt as well as the ectopic and Δpks mutants (Figure 1B). However, the growth of the ΔveA mutants was significantly lower than that of the wt and ectopic mutants (Figure 1E). The conidiation of the three Δpks mutants and the three ΔveA mutants was lower than for the wt strain (Figure 1C,F). The ΔveA mutants displayed a marked reduction in conidiation, approximately 50 % compared to the wt strain. The OTA production analysis showed that none of the Δpks and ΔveA knockout mutants was able to produce the mycotoxin under the assayed conditions (Figure 1D,G). These results were confirmed by HPLC-MS (data not shown).

2.2. Different Stresses Did Not Affect the Growth of the Knockout Mutants

In order to investigate the sensitivity of mutants Δpks and ΔveA to several stresses, the potato dextrose broth (PDB) supplemented with several concentrations of different stressor compounds was tested (Figure 2 and Figure 3). The assay included high osmolarity, tested with sorbitol, osmotic stress induced by sodium salt (NaCl) and compounds that affect the cell wall and cell membranes, such as calcofluor white (CFW) and sodium dodecyl sulfate (SDS), respectively. Oxidative stress was tested with different oxygen peroxide (H2O2) concentrations. The wt strain growth was not affected by increasing the H2O2, CFW or sorbitol concentrations, but it diminished at pH 7.5 compared to the more acidic pHs of 3.0, 4.5 and 6.0. Moreover, its growth was severely affected by adding high NaCl and SDS concentrations. Compared to the wt, the Δpks mutants did not show significant differences for any of the tested compounds or pHs (Figure 2).
No significant differences were observed for the ΔveA mutants compared to the wt when H2O2, SDS or sorbitol were added to culture media (Figure 3). However, the addition of CFW at 1125 and 2250 µg/mL reduced the growth of both ΔveA mutants. Although the addition of 1125 mM NaCl modified their growth, the overall curve profile was the same for the three strains. The pH of the medium affected the growth of the ΔveA knockout mutants, which was higher at pH 3.0 and lower at 7.5 compared to the growth of the wt.

2.3. Competitiveness of Mutants Δpks and ΔveA during In Vitro Growth

In order to assess the ability of mutants Δpks and ΔveA to compete with the wt strain, we co-inoculated the wt and the mutants at different ratios (10 wt:1 Δ, 1 wt:1 Δ, 1 wt:10 Δ) on a 96-well microplate. The control wells were inoculated with either the wt or knockout mutants, and they were cultured under the same conditions. The percentage of each strain at the beginning of the experiment was verified by the colony-counting method, and by inoculating the PDA plates supplemented, or not, with antibiotic (Figure S2A,B). Only the knockout strains were able to develop on the antibiotic-supplemented PDA plates. Each well was collected independently after 7 days of incubation at 28 °C for further analyses. The percentage of each strain under competition was calculated by the colony-counting (Figure S2C,D) and qPCR (Figure 4A,B) methods. Based on the qPCR data, when the mutants and the wt strain were inoculated at the same ratio (1 wt: 1 Δ), mutants were capable of displacing the wt strain at 14.5% wt: 85.5% Δpks, and 25.6% wt: 74.4% ΔveA (Figure 4A,B). Furthermore, mutants were able to displace the wt even under the most unfavorable condition (10 wt:1 Δ ratio). The Δpks mutant succeeded in imposing itself on the wt more intensely than the ΔveA mutant.
Regarding OTA production, the wt strain synthesized OTA, but neither the Δpks (Figure 4C) nor the ΔveA (Figure 4D) mutant produced OTA under the assayed conditions. OTA was not even detected in the 1 wt: 10 Δpks mixture culture. Under the most unfavorable growth condition (10 wt: 1 Δ), OTA production lowered from the expected value of 90.9 % to 62.5 % and 68.0 % when the wt was co-inoculated with the Δpks mutant or the ΔveA mutant, respectively.

2.4. Competitiveness of the Δpks and ΔveA Mutants during Grape Berry Infection

The results of the co-inoculation of the wild-type strain and the mutants in grape berries suggested that both mutants were able to lower the amount of OTA produced by the wt strain (Figure 5). Both mutants Δpks and ΔveA were capable of diminishing OTA production from the expected value of 50.0% to 21.1% and 10.1%, respectively, when the wt strain and mutants were equally inoculated (1 wt:1 Δ). Moreover, the presence of the Δpks mutant reduced OTA production from an expected value of 90.9% to 30.7% under the most adverse condition (10 wt:1 Δ).

3. Discussion

The present study investigated the potential application of non-toxigenic knockout mutants of A. carbonarius for reducing the OTA produced by the ochratoxigenic A. carbonarius ITEM 5010 wt strains. Different commercial products are already available for the application of non-toxigenic A. flavus strains [6,7,26], which supports the notion that competitive exclusion occurs between toxigenic and non-toxigenic strains [27]. Afla-Guard® and AF36 are two commercial biological control products based on non-aflatoxigenic A. flavus strains, which have been approved by the U.S. Environmental Protection Agency for the biological control of A. flavus and aflatoxin contamination in peanut, corn and cottonseed with considerable success [28]. The addition of atoxigenic A. flavus strain A2085 (AF-X1TM) in artificially inoculated maize ears very significantly lowers aflatoxin concentrations [7]. Despite the good results obtained in the control of aflatoxins produced by A. flavus, there are no commercial products available to minimize the OTA produced by A. carbonarius. Numerous research works on this fungus have focused on natural strains that do not produce OTA and on deciphering the genetic causes of non-mycotoxin production. A common putative cluster of genes involved in OTA biosynthesis has been described in A. carbonarius [11,12,13], and OTA production is regulated by different biosynthetic genes, such as genes otaA-AcOTApks and otaB-AcOTAnrps, among others [15,16,17]. Other general regulatory genes, such as veA and laeA, impact mycotoxin production [16,21]. Previous research has demonstrated that co-inoculation of an OTA-producing A. carbonarius strain with a non-OTA-producing A. carbonarius strain, which harbors a Y728H mutation in the otaA gene, can lower OTA levels [29]. If, as suggested by the authors, failed OTA production in this strain is due to a single point mutation, then there is a chance that this mutation can be reversed and lead back to an OTA-producer strain. It is well-known that climate change can affect the production of mycotoxins, and an increase in CO2 on a grape-based matrix often stimulates the development of A. carbonarius and OTA production [22]. With the present climate change scenario, it is critical to ensure that mycotoxins are never synthesized despite changing environmental conditions.
In this work, we specifically wished to determine whether the use of non ochratoxigenic strains obtained by deleting either the otaA genes or veA genes could reduce OTA contamination in foods when competing with their parental mycotoxigenic strain. The deletion of the A. carbonarius otaA gene led to reduced conidiation (Figure 1C) and, as expected, no OTA production was detected during either in vitro growth on rich medium (Figure 1D) or grape berry infection (Figure 5). Other authors have also reported similar results about Δpks mutants’ lack of OTA production [15,16]. Nevertheless, when growing on a minimal medium at 25 °C, Gallo et al. [15] did not observe any phenotypic differences between the wt A. carbonarius ITEM 5010 and the ΔAcOTApks strain. When grown on yeast extract sucrose medium at pH 4.0 and 28 °C, A. carbonarius isolate NRRL 368 and the Acpks-knockout mutant showed comparable growth, sporulation and conidial germination patterns [16]. The different strains utilized in these research works, the incubation temperature and the culture medium are all variables that may influence how the fungus behaves. However, none of the mutants produced OTA, which confirms that this gene is essential for the synthesis of this mycotoxin.
The heterodimeric VELVET complex VelB/VeA/LaeA is a fungus-specific protein complex that controls development, secondary metabolism and pathogenicity. In A. carbonarius, VeA gene deletion resulted in decreased conidiation, which indicates that this gene is essential for the strain’s ability to produce conidia (Figure 1) [21]. In contrast, VeA acts as a negative regulator of asexual development in Aspergillus cristatus and Aspergillus nidulans, among others [30,31]. Although it would appear that this gene plays a species-specific role in the regulation of conidia formation, its role in the production of several mycotoxins is clear: although VeA generally plays a positive role in the regulation of secondary metabolite production, some secondary metabolites are negatively regulated by VeA [24,32]. In line with that previously shown by Crespo-Sempere et al. [21], OTA production was affected by the deletion of the veA gene in A. carbonarius (Figure 1). This gene also positively regulates OTA biosynthesis in Aspergillus niger and Aspergillus ochraceus [33,34], and it plays a role in the production of aflatoxins in Aspergillus flavus [35] and fumonisin by Fusarium verticillioides [36]. Light, pH, and osmotic and oxidative stresses are a few examples of environmental stressors that VeA responds to. For instance, in the dark, the OTA production in the AcΔveA mutant decreased to a greater extent compared to the wt A. carbonarius when both were grown under light conditions. Nevertheless, conidia production similarly diminished in both situations compared to the wt [21]. A. carbonarius exhibits different behavior depending on the pH of the environment. At pH 7.0, it produces less OTA than at pH 4.0 [37]. Our results showed that the growth of both the Δpks mutant and the wt was marginally lower under the neutral pH 7.5 conditions than under more acidic conditions (Figure 2). However, the ΔveA mutant developed more quickly under the most acidic situation that we examined but less rapidly under neutral conditions (Figure 3). Previous studies have shown that the disruption of the veA gene reduces tolerance to oxidative stress in A. flavus and A. niger [34,38], but it has no effect on A. fumigatus [39]. Our study demonstrated that loss of genes veA and otaA in A. carbonarius did not change the behavior of the mutants in a liquid medium with hydrogen peroxide employed as an oxidant inducer. This finding suggests that veA performs no function in preventing oxidative stress. We investigated how osmotic stressors, such as NaCl and sorbitol, affect growth in A. carbonarius, and whether genes otaA and veA could impact this. Our results revealed that osmotic stress caused by high NaCl concentrations leads to diminished fungal growth. However, the growth in both mutants was unaffected compared to the wt, which is similar to that described in A. flavus, where the response to osmotic is barely affected by VeA [38].
The parallel differentiation and quantification of A. carbonarius toxigenic and non-toxigenic strains in mixed cultures were completed by qPCR following a similar strategy to that used to discriminate non aflatoxigenic biocontrol strains and aflatoxigenic strains of A. flavus by droplet digital PCR and qPCR [40,41]. The qPCR results (Figure 4) were the equivalent to those obtained by the traditional plate-counting approach using PDA supplemented, or not, with Hygromicin B as a selective marker for mutants (Figure S2). Even under the most unfavorable initial inoculation conditions (10 wt:1 Δ), both mutants Δpks and ΔveA were able to displace the wt strain. A low percentage of the non-toxigenic A. carbonarius was able to significantly lower the OTA level in both culture medium (Figure 4) and infected grapes (Figure 5). The overall results herein described indicate that the knockout mutants reduced OTA production to a greater extent than growth. It is worth noting that the co-inoculation of grapes with a low percentage of the Δpks mutant resulted in a more pronounced OTA reduction than with the ΔveA mutant. This suggests that the ΔveA mutant displays less fitness than the Δpks mutant during grape infection in the presence of the wt strain. This differential fitness could be related to the role of the missing proteins. As previously mentioned, OtaA is a polyketide synthase that contributes specifically to OTA biosynthesis, while protein VeA is involved in global regulation and affects several functions, such as sexual development, secondary metabolism, and even pathogenicity. As VeA is a global regulator of secondary metabolism, it could be advantageous to utilize ΔveA mutants if the goal of using biocontrol agents is to eliminate the possibility of any mycotoxin being present, and not only OTA. Further experiments should be run to test this hypothesis.
The preliminary results of this study using these mutants as biocontrol agents are promising for the control of OTA. Similarly to what happens with A. flavus [7], it could be feasible to employ non ochratoxigenic mutants to outcompete toxigenic A. carbonarius strains to reduce OTA contamination in food and feed.

4. Materials and Methods

4.1. Fungal Strains and Culture Conditions

Aspergillus carbonarius ITEM 5010 was used as the wt strain to construct knockout mutants for genes otaA and veA. All the strains were grown on potato dextrose agar (PDA, Difco-BD Diagnostics, Sparks, MD, USA) or PDB, with or without the corresponding antibiotic. Cultures were incubated at 28 °C for 7–14 days. Conidia were scraped off agar with a sterile spatula, suspended in distilled water and titrated by a hemacytometer. Plasmid vectors were cloned and propagated in Escherichia coli DH5α grown in Luria–Bertani medium (LB; bacto tryptone 10 g, yeast extract 5 g, NaCl 5 g, agar 15 g, per liter) supplemented with 25 µg/mL of kanamycin at 37 °C. Agrobacterium tumefaciens AGL-1 was grown in the LB medium supplemented with 20 µg/mL of rifampicin, 100 µg/mL of kanamycin and 75 µg/mL of carbenicillin at 28 °C.

4.2. Construction and Verification of the A. carbonarius Knockout Mutants

In order to construct the otaA and veA gene replacement plasmids, the 5’ and 3’ flanking regions (≈1.5 kb) were amplified from A. carbonarius ITEM 5010 genomic DNA and cloned into plasmid vector pRF-HU2 [25]. All the primers pairs were designed by the Primer3 software [42]. The amplification of 10 ng of genomic DNA was performed using High-Taq DNA polymerase (Bioron GmbH, Ludwigshafen, Germany) according to the manufacturer’s instructions and using primers pairs: veA-VA-O1/veA-VA-O2 and OTApks_O1/OTApks_O2 for the upstream regions; veA-VA-A3/veA-VA-A4 and OTApks_A3/OTApks_A4 for the downstream regions (Table S2). The cycling conditions consisted of: 94 °C for 3 min, 35 cycles of 94 °C for 15 s, 58 °C for 20 s and 72 °C for 1 min 30 s and 72 °C for 10 min.
Binary vector pRF-HU2 was designed to be used with the USER (Uracil-Specific Excision Reagent) friendly cloning technique (New England Biolabs, Ipswich, MA, USA), as previously described [43]. First to obtain plasmids pRFHU2-VEA and pRFHU2-PKS (Figure S1), the amplified upstream and downstream fragment regions were mixed with digested vector pRF-HU2 (ratio of 30:30:120 ng) and the USER enzyme mix, and then, they were incubated according to the manufacturer’s conditions. An aliquot of the mixture was directly used for the transformation of the E. coli DH5α chemically-competent cells. Kanamycin-resistant transformants were screened by PCR using primers pairs RF-5/RF-2 and RF-1/RF-6 (Table S2). Proper fusion was confirmed by DNA sequencing.
The transformation of A. tumefaciens was completed by introducing both plasmids independently in electrocompetent A. tumefaciens AGL-1 cells with a Gene Pulser apparatus (Bio-Rad, Richmond, CA, USA) [44]. Then, the final A. carbonarius transformation was completed as previously described [44] by mixing equal volumes of a conidial suspension of A. carbonarius (104, 105 and 106 conidia/mL) and IMAS-induced A. tumefaciens cultures, and spreading them onto paper filter layered on agar plates containing IMAS. After co-cultivation at 24 °C for 40 h, membranes were transferred to PDA plates containing hygromycin B (HygB, 100 µg/mL, InvivoGen, San Diego, CA, USA) to select fungal transformants and cefotaxime (200 µg/mL, Sigma-Aldrich, St. Louis, MO, USA) as an inhibitory growth agent of A. tumefaciens cells. Finally, A. carbonarius ITEM 5010 hygromycin-resistant colonies appeared after 3 to 4 days of incubation at 28 °C, which allowed monosporic cultures to be subsequently obtained.
Genomic DNA extraction was completed as formerly described [43]. The PCR analysis of the transformants was used to confirm the disruption of genes veA and otaA: screening the deletion of the veA gene (primers pairs: veA-VI/veA-VJ), the deletion of the otaA gene (primers pairs: OTApks-F3/OTApks-R4), and the insertion of selection marker HygB (primers pairs: HMBF1/HMBR1) (Table S2). Finally, to check the number of T–DNA copies integrated into the genome of transformants, qPCR was carried out as previously described [44]. The primers pairs used to determine the T-DNA copy number of genes veA and otaA were designed in the upstream region of the genes with primers pairs veA-VG/veA-VH, and OTApks-R4/OTApks-F8, respectively. The gene selected as a reference was the non-ribosomal peptide synthetase (nrps) gene (ID: 132610), which was detected with primer pairs AcNRP_F/AcNRP_R (Table S1) [21]. Then, qPCR reactions were performed as formerly reported [44] in a LightCycler 480 Instrument (Roche Diagnostics, Mannheim, Germany) equipped with the LightCycler SW 1.5 software and by calculating the number of T–DNA copies integrated into the genome of each transformant according to Pfaffl [45].

4.3. Characterization of the Knockout Mutants: Mycelial Growth, Conidiation and Growth under Stress Conditions

For the phenotypical characterization, 5 µL of conidial suspension (1 × 105 conidia/mL) was placed in the center of PDA plates and incubated for up to 5 days at 28 °C. Plates were scanned daily, and the area of growth was analyzed with ImageJ 1.53q (Wayne Rasband, National Institutes of Health, Bethesda, MD, USA). To determine conidia production, three plugs (5 mm diameter) were collected from the center, middle and border of the colony, and 500 µL of methanol was added to the three plugs. The mixture was shaken in an Omni bead Ruptor 24 (Omni International Inc., Kennesaw, GE, USA) for 1 min and at speed 4. Conidia were counted with a hemacytometer and methanol extracts were stored at −80 °C to determine OTA production.
Growth profiles were performed on the 96-well PDB plates supplemented with different concentrations of compounds H2O2, SDS, CFW, NaCl and sorbitol, and at distinct pHs (3.0, 4.5, 6.0, and 7.5). Spores were collected after 7 days of growth from the pointed-centrally inoculated PDA plates, and concentration was adjusted using a hemocytometer. To determine growth profiles, the 96-well PDB plates containing 100 µL of medium were inoculated in triplicate at a final concentration of 1 × 105 conidia/mL, and they were incubated at 24 °C for up to 7 days. Absorbance at 600 nm was measured automatically at 2 hourly intervals with a FLUOstar Omega (BMG Labtech). Growth curve analyses were performed in R 4.1.1. The area under the curve (AUC) used to describe the growth of the fungi under different stress conditions was calculated at 7 dpi for each replicate based on the spline fit model included in the ‘grofit’ package [46].

4.4. In Vitro Competition Assays

Under controlled laboratory conditions, the competition assays involving the ochratoxigenic A. carbonarius ITEM 5010 strain and the non-mycotoxigenic mutants were performed on 96-well PDB plates. One knockout mutant for each gene (ΔveA10b and Δpks8c) was selected. Spore suspension concentrations were adjusted to 1 × 105 conidia/mL. Five 10:0, 10:1, 1:1, 1:10, and 0:10 mix ratios (the wt strain vs. Δ) were used to inoculate PDB on 96-wells plates. Cultures were grown for 7 days at 28 °C in the dark. Each competition ratio was performed in 15 wells and was repeated twice on different days. At the end of the competition assays, five wells were used to estimate the growth percentage of each competing strain. Then plates were stored at −20 °C for further analyses.
Each strain’s percentage of growth was calculated by two different approaches: (i) diluting the content of the well with 1 mL of water, counting the conidia concentration in at least five wells and plating 100 µL of 2 × 103 conidia/mL onto PDA plates and the PDA plates containing 100 µg/mL of hygromycin. After 24–48 h at 28 °C, incipient colonies were counted, and the percentages of the wt and knockout mutants were determined for each ratio. Only the knockout mutants were able to grow on the DA plates containing hygromycin; (ii) quantifying fitness by qPCR. At least three wells per ratio were used for DNA extraction. The DNAs of A. carbonarius ITEM5010 and the knockout mutants were used as a proxy to determine each strain’s growth percentage by employing specific primers: (i) AcveA-VI/AcveA-VJ, and OTApks-F5/OTApks-R6 to detect gene veA and gene otaA, respectively, in the wt strain; (ii) HPH3F/HPH4R for the hygromycin-resistant marker in the knockout mutants; (iii) AcNRP_F/AcNRP_R to detect the reference gene [21]. A qPCR analysis was performed as previously described.
To determine whether non-mycotoxigenic mutants could inhibit ochratoxin A production by the wt strain A. carbonarius ITEM 5010, OTA was independently extracted from at least three co-inoculated PDB wells in 2 mL tubes containing three 2.7 mm steel beads with 500 µL of methanol with the help of an Omni Bead Ruptor 24 (Omni International, Inc., Kennesaw, GE, USA). Samples were centrifuged and filtered, and methanol extracts were stored at −80 °C until further analyses.

4.5. Artificial Inoculation of Grape Berries

For each strain, nine grape berries were surface-sterilized with 2% sodium hypochlorite for 5 min, rinsed with tap water and air-dried. Conidia, collected by scraping the surface of 7-day-old colonies grown on PDA plates, were resuspended in water and adjusted to 1 × 105 conidia/mL. Five (wt strain vs. deletant) 10:0, 10:1, 1:1, 1:10 and 0:10 mix ratios were prepared. Aliquots of 5 µL of the conidial suspensions were inoculated on the grape berries that had been previously wounded with a needle at a 5 mm depth. After 10 days at 28 °C in the dark, three biological replicates containing all three grape berries were collected, frozen in liquid nitrogen and ground to a fine powder. Pooled samples were stored at −80 °C until further use. To determine OTA levels, 0.5 g of frozen tissues was mixed with 250 µL of methanol and three steal beads, as previously described, and the methanol extracts were stored at −80 °C to determine OTA production.

4.6. OTA Quantification

OTA was detected by a previously described high-pressure chromatography (HPLC) method [44] with minor modifications. Briefly, methanolic extracts were filtered and injected into HPLC. OTA detection and quantification were completed by an HPLC system (ACQUITY Arc Sys Core 1–30 cm, Waters Co., Milford, CT, USA) equipped with a Waters temperature control module, a Waters 2475 fluorescence detector (excitation wavelength of 330 nm and emission wavelength of 460 nm), and a Kinetex biphenyl column (4.6 mm × 150 mm, 5 µm, Phenomenex Inc., Torrance, CA, USA). Twenty microliters of each extract were injected. The mobile phase was acetonitrile: water: acetic acid (57:41:2, v/v/v) with isocratic elution for 15 min at a flow rate of 1 mL/min. Working standard solutions were prepared by appropriate diluting the known volumes of the stock solution with methanol, and they were used to obtain calibration curves in the chromatographic system. OTA was expressed as a percentage of the wt strain.
Samples were also analyzed by mass spectrometry with a Waters Acquity UPLC I-Class System (Waters Corporation, Milford, MA, USA) coupled to a Bruker Daltonics QToF-MS mass spectrometer (maXis impact series, resolution ≥ 55,000 FWHM Bruker Daltonics, Bremen, Germany) using ESI for the positive (ESI (+)) ionization mode. UPLC separation was performed in an ACE-Excel C18-PFP column (3.0 mm × 100 mm, 1.7 µm size particle) at a flow rate of 0.3 mL/min. Separation was completed using water with formic acid at 0.01% as the weak mobile phase (A) and methanol with formic acid at 0.01% as the strong mobile phase (B). The gradient started with 20% of B at 0 min, which was maintained for 3 min. Next, the proportion was set at 40% and linearly increased to 80% until minute 22.50. Eluent B was set at 20% and kept constant until minute 30. Nitrogen was used as the desolvation gas with a flux of 9 L/min and as the nebulizing gas with a flux of 2.0 bar. The drying temperature was 200 °C and the column temperature was 40 °C. The voltage source was 4.0 kV for ESI (+). The MS experiment was completed by employing HR-QTOF-MS, applying 24 eV for ESI (+) and using broadband collision-induced dissociation (bbCID). The MS data were acquired within an m/z range of 50–1200 Da.
The external calibrant solution was delivered by a KNAUER Smartline Pump 100 with a pressure sensor (KNAUER, Berlin, Germany). The instrument was externally calibrated before each sequence with 10 mM sodium formate solution. The mixture was prepared by adding 0.5 mL of formic acid and 1.0 mL of 1.0 M sodium hydroxide to an isopropanol/Milli Q water solution (1:1, v/v). UPLC-QToF-MS analyses were performed at the Metabolomics Platform of CEBAS-CSIC, Campus Universitario de Espinardo, 30100 Espinardo, Murcia (Spain).

4.7. Statistical Analysis

The ANOVA analysis was used to determine if there were significant differences among means. Tukey’s test was carried out to determine if significant (p < 0.05) differences occurred between individual treatments (Statpoint).

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/toxins14110745/s1. Figure S1: Deletion of genes otaA and veA in Aspergillus carbonarius. (A) Physical maps of plasmids pRFHU2-PKS and pRFHU2-VEA. (B) Diagram of the deletion cassette used to replace the target region (gene of interest) in the wild-type (wt) strain by homologous recombination by generating mutants Δpks and ΔveA. The primers used in the construction and analysis of both plasmids are shown. (C) Amplification band patterns of the different polymerase chain reactions (PCR) of the wt ITEM 5010 (wt), two ectopic mutants (preceded by letter ‘e’) and three knockout mutants (preceded by symbol ‘Δ’). Figure S2: Competitiveness of knockout mutants Δpks (A–C, pink bars) and ΔveA (B–D, blue bars) against the mycotoxigenic wt strain A. carbonarius ITEM 5010 (black bars) at day 0 (A,B) and 7 days post inoculation (C,D). Competitiveness was determined by counting colonies on PDA and the PDA supplemented with 100 µg/mL of hygromycin. Only knockout mutants can grow on the PDA supplemented with antibiotic. Values are the mean of at least three biological replicates and error bars represent the standard error of the mean (SEM). Table S1: Estimation of the number of T-DNA copies integrated into the genome of the mutants. Table S2: List of the primers used in this study.

Author Contributions

Conceptualization, L.G.-C. and A.-R.B.; methodology, B.L. and A.-R.B.; formal analysis, B.L. and A.-R.B.; writing—original draft preparation, B.L. and A.-R.B.; writing—review and editing, L.G.-C. and A.-R.B.; supervision, L.G.-C. and A.-R.B.; funding acquisition, L.G.-C. and A.-R.B. All authors have read and agreed to the published version of the manuscript.

Funding

The work was partially cofunded by Projects RTI2018-093392-A-I00 and PID2021-126005OB-I00 (MCIN /AEI /10.13039/501100011033 / FEDER, UE) and by Project AGROALNEXT/2022/028 (GVA, PRTR, MRR, NextGenerationEU). The Ramón y Cajal postdoctoral contract (RYC-2017-22009, MINECO and FSE) to A.-R.B. is also acknowledged, as are contracts JAEINT_19_01896 (CSIC) and PRE2019-089326 (MICINN/AEI/FSE, UE) to B.L.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors acknowledge the technical assistance of Adrián González in the final working stages. The authors acknowledge the support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bennett, J.W.; Klich, M. Mycotoxins. Clin. Microbiol. Rev. 2003, 16, 497–516. [Google Scholar] [CrossRef] [Green Version]
  2. Eskola, M.; Kos, G.; Elliott, C.T.; Hajšlová, J.; Mayar, S.; Krska, R. Worldwide Contamination of Food-Crops with Mycotoxins: Validity of the Widely Cited ‘FAO Estimate’ of 25%. Crit. Rev. Food Sci. Nutr. 2020, 60, 2773–2789. [Google Scholar] [CrossRef]
  3. Liu, Y.; Yamdeu, J.H.G.; Gong, Y.Y.; Orfila, C. A Review of Postharvest Approaches to Reduce Fungal and Mycotoxin Contamination of Foods. Compr. Rev. Food Sci. Food Saf. 2020, 19, 1521–1560. [Google Scholar] [CrossRef]
  4. EPA. US Environmental Protection Agency Office of Pesticide Programs. In Biopesticide Registration Action Document: Aspergillus flavus AF36; PC Code 006456; EPA: Washington, DC, USA, 2003. [Google Scholar]
  5. EPA. US Environmental Protection Agency Office of Pesticide Programs. In Biopesticide Registration Action Document Aspergillus flavus (NRRL 21882) (PC Code 006500); EPA: Washington, DC, USA, 2004. [Google Scholar]
  6. Dorner, J.; Lamb, M. Development and Commercial Use of Afla-Guard®, an Aflatoxin Biocontrol Agent. Mycotoxin Res. 2006, 22, 33–38. [Google Scholar] [CrossRef]
  7. Mauro, A.; Garcia-Cela, E.; Pietri, A.; Cotty, P.J.; Battilani, P. Biological Control Products for Aflatoxin Prevention in Italy: Commercial Field Evaluation of Atoxigenic Aspergillus flavus Active Ingredients. Toxins 2018, 10, E30. [Google Scholar] [CrossRef] [Green Version]
  8. Ehrlich, K.C.; Cotty, P.J. An Isolate of Aspergillus flavus Used to Reduce Aflatoxin Contamination in Cottonseed Has a Defective Polyketide Synthase Gene. Appl. Microbiol. Biotechnol. 2004, 65, 473–478. [Google Scholar] [CrossRef]
  9. Kolawole, O.; Meneely, J.; Petchkongkaew, A.; Elliott, C. A Review of Mycotoxin Biosynthetic Pathways: Associated Genes and Their Expressions under the Influence of Climatic Factors. Fungal Biol. Rev. 2021, 37, 8–26. [Google Scholar] [CrossRef]
  10. Perrone, G.; Ferrara, M.; Medina, A.; Pascale, M.; Magan, N. Toxigenic Fungi and Mycotoxins in a Climate Change Scenario: Ecology, Genomics, Distribution, Prediction and Prevention of the Risk. Microorganisms 2020, 8, 1496. [Google Scholar] [CrossRef]
  11. Gil-Serna, J.; García-Díaz, M.; González-Jaén, M.T.; Vázquez, C.; Patiño, B. Description of an Orthologous Cluster of Ochratoxin A Biosynthetic Genes in Aspergillus and Penicillium Species. A Comparative Analysis. Int. J. Food Microbiol. 2018, 268, 35–43. [Google Scholar] [CrossRef]
  12. Ferrara, M.; Gallo, A.; Perrone, G.; Magistà, D.; Baker, S.E. Comparative Genomic Analysis of Ochratoxin A Biosynthetic Cluster in Producing Fungi: New Evidence of a Cyclase Gene Involvement. Front. Microbiol. 2020, 11, 581309. [Google Scholar] [CrossRef]
  13. Castellá, G.; Bragulat, M.R.; Puig, L.; Sanseverino, W.; Cabañes, F.J. Genomic Diversity in Ochratoxigenic and Non Ochratoxigenic Strains of Aspergillus carbonarius. Sci. Rep. 2018, 8, 5439. [Google Scholar] [CrossRef] [Green Version]
  14. Wang, Y.; Wang, L.; Wu, F.; Liu, F.; Wang, Q.; Zhang, X.; Selvaraj, J.N.; Zhao, Y.; Xing, F.; Yin, W.-B.; et al. A Consensus Ochratoxin A Biosynthetic Pathway: Insights from the Genome Sequence of Aspergillus ochraceus and a Comparative Genomic Analysis. Appl. Environ. Microbiol. 2018, 84, e01009-18. [Google Scholar] [CrossRef] [Green Version]
  15. Gallo, A.; Knox, B.P.; Bruno, K.S.; Solfrizzo, M.; Baker, S.E.; Perrone, G. Identification and Characterization of the Polyketide Synthase Involved in Ochratoxin A Biosynthesis in Aspergillus carbonarius. Int. J. Food Microbiol. 2014, 179, 10–17. [Google Scholar] [CrossRef]
  16. Maor, U.; Barda, O.; Sadhasivam, S.; Bi, Y.; Levin, E.; Zakin, V.; Prusky, D.B.; Sionov, E. Functional Roles of LaeA, Polyketide Synthase, and Glucose Oxidase in the Regulation of Ochratoxin A Biosynthesis and Virulence in Aspergillus carbonarius. Mol. Plant Pathol. 2021, 22, 117–129. [Google Scholar] [CrossRef]
  17. Gallo, A.; Bruno, K.S.; Solfrizzo, M.; Perrone, G.; Mulè, G.; Visconti, A.; Baker, S.E. New Insight into the Ochratoxin A Biosynthetic Pathway through Deletion of a Nonribosomal Peptide Synthetase Gene in Aspergillus carbonarius. Appl. Environ. Microbiol. 2012, 78, 8208–8218. [Google Scholar] [CrossRef] [Green Version]
  18. Ferrara, M.; Perrone, G.; Gambacorta, L.; Epifani, F.; Solfrizzo, M.; Gallo, A. Identification of a Halogenase Involved in the Biosynthesis of Ochratoxin A in Aspergillus carbonarius. Appl. Environ. Microbiol. 2016, 82, 5631–5641. [Google Scholar] [CrossRef] [Green Version]
  19. Gerin, D.; Garrapa, F.; Ballester, A.-R.; González-Candelas, L.; De Miccolis Angelini, R.M.; Faretra, F.; Pollastro, S. Functional Role of Aspergillus carbonarius AcOTAbZIP Gene, a BZIP Transcription Factor within the OTA Gene Cluster. Toxins 2021, 13, 111. [Google Scholar] [CrossRef]
  20. Ferrara, M.; Gallo, A.; Cervini, C.; Gambacorta, L.; Solfrizzo, M.; Baker, S.E.; Perrone, G. Evidence of the Involvement of a Cyclase Gene in the Biosynthesis of Ochratoxin A in Aspergillus carbonarius. Toxins 2021, 13, 892. [Google Scholar] [CrossRef]
  21. Crespo-Sempere, A.; Marín, S.; Sanchis, V.; Ramos, A.J. VeA and LaeA Transcriptional Factors Regulate Ochratoxin A Biosynthesis in Aspergillus carbonarius. Int. J. Food Microbiol. 2013, 166, 479–486. [Google Scholar] [CrossRef]
  22. Cervini, C.; Verheecke-Vaessen, C.; Ferrara, M.; García-Cela, E.; Magistà, D.; Medina, A.; Gallo, A.; Magan, N.; Perrone, G. Interacting Climate Change Factors (CO2 and Temperature Cycles) Effects on Growth, Secondary Metabolite Gene Expression and Phenotypic Ochratoxin A Production by Aspergillus carbonarius Strains on a Grape-Based Matrix. Fungal Biol. 2021, 125, 115–122. [Google Scholar] [CrossRef]
  23. Assaf, C.E.H.; Snini, S.P.; Tadrist, S.; Bailly, S.; Naylies, C.; Oswald, S.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] [Green Version]
  24. Wiemann, P.; Brown, D.W.; Kleigrewe, K.; Bok, J.W.; Keller, N.P.; Humpf, H.U.; Tudzynski, B. FfVel1 and FfLae1, Components of a Velvet-like Complex in Fusarium fujikuroi, Affect Differentiation, Secondary Metabolism and Virulence. Mol. Microbiol. 2010, 77, 972–994. [Google Scholar] [CrossRef] [Green Version]
  25. Frandsen, R.J.; Andersson, J.A.; Kristensen, M.B.; Giese, H. Efficient Four Fragment Cloning for the Construction of Vectors for Targeted Gene Replacement in Filamentous Fungi. BMC Mol. Biol. 2008, 9, 70. [Google Scholar] [CrossRef] [Green Version]
  26. Dorner, J.W. Biological Control of Aflatoxin Contamination in Corn Using a Nontoxigenic Strain of Aspergillus flavus. J. Food Prot. 2009, 72, 801–804. [Google Scholar] [CrossRef]
  27. Cotty, P.J.; Bayman, P. Competitive Exclusion of a Toxigenic Strain of Aspergillus flavus by an Atoxigenic Strain. Phytopathology 1993, 83, 1283. [Google Scholar] [CrossRef]
  28. Lewis, M.H.; Carbone, I.; Luis, J.M.; Payne, G.A.; Bowen, K.L.; Hagan, A.K.; Kemerait, R.; Heiniger, R.; Ojiambo, P.S. Biocontrol Strains Differentially Shift the Genetic Structure of Indigenous Soil Populations of Aspergillus flavus. Front. Microbiol. 2019, 10, 1738. [Google Scholar] [CrossRef] [Green Version]
  29. Castellá, G.; Bragulat, M.R.; Cigliano, R.A.; Cabañes, F.J. Transcriptome Analysis of Non-Ochratoxigenic Aspergillus carbonarius Strains and Interactions between Some Black Aspergilli Species. Int. J. Food Microbiol. 2020, 317, 108498. [Google Scholar] [CrossRef] [Green Version]
  30. Tan, Y.; Wang, H.; Wang, Y.; Ge, Y.; Ren, X.; Ren, C.; Wang, Y.; Ren, X.; Liu, Y.; Liu, Z. The Role of the VeA Gene in Adjusting Developmental Balance and Environmental Stress Response in Aspergillus cristatus. Fungal Biol. 2018, 122, 952–964. [Google Scholar] [CrossRef]
  31. Bayram, Ö.; Krappmann, S.; Ni, M.; Bok, J.W.; Helmstaedt, K.; Valerius, O.; Braus-Stromeyer, S.; Kwon, N.-J.; Keller, N.P.; Yu, J.-H.; et al. VelB/VeA/LaeA Complex Coordinates Light Signal with Fungal Development and Secondary Metabolism. Science 2008, 320, 1504–1506. [Google Scholar] [CrossRef]
  32. Bok, J.W.; Soukup, A.A.; Chadwick, E.; Chiang, Y.-M.; Wang, C.C.C.; Keller, N.P. VeA and MvlA Repression of the Cryptic Orsellinic Acid Gene Cluster in Aspergillus nidulans Involves Histone 3 Acetylation. Mol. Microbiol. 2013, 89, 963–974. [Google Scholar] [CrossRef]
  33. Wang, G.; Zhang, H.; Wang, Y.; Liu, F.; Li, E.; Ma, J.; Yang, B.; Zhang, C.; Li, L.; Liu, Y. Requirement of LaeA, VeA, and VelB on Asexual Development, Ochratoxin A Biosynthesis, and Fungal Virulence in Aspergillus ochraceus. Front. Microbiol. 2019, 10, 2759. [Google Scholar] [CrossRef] [Green Version]
  34. Zhang, J.; Chen, H.; Sumarah, M.W.; Gao, Q.; Wang, D.; Zhang, Y. VeA Gene Acts as a Positive Regulator of Conidia Production, Ochratoxin a Biosynthesis, and Oxidative Stress Tolerance in Aspergillus niger. J. Agric. Food Chem. 2018, 66, 13199–13208. [Google Scholar] [CrossRef]
  35. Amaike, S.; Keller, N.P. Distinct Roles for VeA and LaeA in Development and Pathogenesis of Aspergillus flavus. Eukaryot. Cell 2009, 8, 1051–1060. [Google Scholar] [CrossRef] [Green Version]
  36. Myung, K.; Zitomer, N.C.; Duvall, M.; Glenn, A.E.; Riley, R.T.; Calvo, A.M. The Conserved Global Regulator VeA Is Necessary for Symptom Production and Mycotoxin Synthesis in Maize Seedlings by Fusarium verticillioides. Plant Pathol. 2012, 61, 152–160. [Google Scholar] [CrossRef] [Green Version]
  37. Barda, O.; Maor, U.; Sadhasivam, S.; Bi, Y.; Zakin, V.; Prusky, D.; Sionov, E. The PH-Responsive Transcription Factor PacC Governs Pathogenicity and Ochratoxin A Biosynthesis in Aspergillus carbonarius. Front. Microbiol. 2020, 11, 210. [Google Scholar] [CrossRef]
  38. Baidya, S.; Duran, R.M.; Lohmar, J.M.; Harris-Coward, P.Y.; Cary, J.W.; Hong, S.-Y.; Roze, L.V.; Linz, J.E.; Calvo, A.M. VeA Is Associated with the Response to Oxidative Stress in the Aflatoxin Producer Aspergillus flavus. Eukaryot. Cell 2014, 13, 1095–1103. [Google Scholar] [CrossRef] [Green Version]
  39. Dhingra, S.; Andes, D.; Calvo, A.M. VeA Regulates Conidiation, Gliotoxin Production and Protease Activity in the Opportunistic Human Pathogen Aspergillus fumigatus. Eukaryot Cell 2012, 11, 1531–1543. [Google Scholar] [CrossRef] [Green Version]
  40. Schamann, A.; Schmidt-Heydt, M.; Geisen, R. Analysis of the Competitiveness between a Non-Aflatoxigenic and an Aflatoxigenic Aspergillus flavus Strain on Maize Kernels by Droplet Digital PCR. Mycotoxin Res. 2022, 38, 27–36. [Google Scholar] [CrossRef]
  41. Garcia-Lopez, M.T.; Luo, Y.; Ortega-Beltran, A.; Jaime, R.; Moral, J.; Michailides, T.J. Quantification of the Aflatoxin Biocontrol Strain Aspergillus flavus AF36 in Soil and in Nuts and Leaves of Pistachio by Real-Time PCR. Plant Dis. 2021, 105, 1657–1665. [Google Scholar] [CrossRef]
  42. Untergasser, A.; Cutcutache, I.; Koressaar, T.; Ye, J.; Faircloth, B.C.; Remm, M.; Rozen, S.G. Primer3—New Capabilities and Interfaces. Nucleic Acids Res. 2012, 40, e115. [Google Scholar] [CrossRef]
  43. Crespo-Sempere, A.; López-Pérez, M.; Martínez-Culebras, P.V.; González-Candelas, L. Development of a Green Fluorescent Tagged Strain of Aspergillus carbonarius to Monitor Fungal Colonization in Grapes. Int. J. Food Microbiol. 2011, 148, 135–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Gerin, D.; González-Candelas, L.; Ballester, A.-R.; Pollastro, S.; De Miccolis Angelini, R.; Faretra, F. Functional Characterization of the Alb1 Orthologue Gene in the Ochratoxigenic Fungus Aspergillus carbonarius (AC49 Strain). Toxins 2018, 10, 120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Pfaffl, M.W. A New Mathematical Model for Relative Quantification in Real-Time RT–PCR. Nucleic Acids Res. 2001, 29, e45. [Google Scholar] [CrossRef] [PubMed]
  46. Kahm, M.; Hasenbrink, G.; Lichtenberg-Fraté, H.; Ludwig, J.; Kschischo, M. Grofit: Fitting Biological Growth Curves with R. J. Stat. Softw. 2010, 33, 1–21. [Google Scholar] [CrossRef]
Figure 1. Phenotypic traits of the wild-type strain of A. carbonarius ITEM 5010 (denoted as ‘wt’, black bars), two ectopic mutants (denoted as ‘E’, gray bars) and three knockout mutants (color-filled bars) of genes otaA (BD) and veA (EG), denoted as Δpks and ΔveA, respectively. (A) The front and back colony views of the different strains point-inoculated on PDA plates in the dark at 7 days post-inoculation. Growth area (B,E), conidiation (C,F) and OTA production (D,G) were tested on the PDA plates centrally point-inoculated with 5 uL of 1 × 105 conidia/mL. After incubation at 28 °C for 5 days, colonies were scanned to determine the growth area with the ImageJ software. Three plugs were collected from the center, middle and inner parts of the colony for conidia determination and OTA extraction purposes. Values were normalized to those of the wt growing under the same conditions. Error bars represent the standard error of the mean of at least three biological replicates. The bars with different letters in the same panel are statistically different as determined by the one-way ANOVA and Tukey’s test (p < 0.05). nd. not detected under the assayed conditions.
Figure 1. Phenotypic traits of the wild-type strain of A. carbonarius ITEM 5010 (denoted as ‘wt’, black bars), two ectopic mutants (denoted as ‘E’, gray bars) and three knockout mutants (color-filled bars) of genes otaA (BD) and veA (EG), denoted as Δpks and ΔveA, respectively. (A) The front and back colony views of the different strains point-inoculated on PDA plates in the dark at 7 days post-inoculation. Growth area (B,E), conidiation (C,F) and OTA production (D,G) were tested on the PDA plates centrally point-inoculated with 5 uL of 1 × 105 conidia/mL. After incubation at 28 °C for 5 days, colonies were scanned to determine the growth area with the ImageJ software. Three plugs were collected from the center, middle and inner parts of the colony for conidia determination and OTA extraction purposes. Values were normalized to those of the wt growing under the same conditions. Error bars represent the standard error of the mean of at least three biological replicates. The bars with different letters in the same panel are statistically different as determined by the one-way ANOVA and Tukey’s test (p < 0.05). nd. not detected under the assayed conditions.
Toxins 14 00745 g001
Figure 2. Growth of the wild-type strain A. carbonarius ITEM 5010 (■) and two Δpks mutants (△ represents Δpks8a, and ▽ denotes Δpks27c) in the presence of different H2O2 (A), CFW (B), NaCl (C), SDS (D) and sorbitol (E) concentrations and at distinct pHs (F). Growth was determined as the area under the curve (AUC) after 7 days of incubation at 24 °C. The two-way ANOVA, followed by Tukey’s test (p < 0.05), was performed to determine the significant growth of the Δpks8a (denoted as *) and Δpks27c (denoted as Φ) knockout mutants compared to the wild type. Values represent the mean ± standard error of the mean of three biological replicates. The experiment was repeated twice.
Figure 2. Growth of the wild-type strain A. carbonarius ITEM 5010 (■) and two Δpks mutants (△ represents Δpks8a, and ▽ denotes Δpks27c) in the presence of different H2O2 (A), CFW (B), NaCl (C), SDS (D) and sorbitol (E) concentrations and at distinct pHs (F). Growth was determined as the area under the curve (AUC) after 7 days of incubation at 24 °C. The two-way ANOVA, followed by Tukey’s test (p < 0.05), was performed to determine the significant growth of the Δpks8a (denoted as *) and Δpks27c (denoted as Φ) knockout mutants compared to the wild type. Values represent the mean ± standard error of the mean of three biological replicates. The experiment was repeated twice.
Toxins 14 00745 g002
Figure 3. Growth of the wt strain A. carbonarius ITEM 5010 (■) and two ΔveA (△ represents ΔveA10b, and ▽ denotes ΔveA12a) in the presence of different H2O2 (A), CFW (B), NaCl (C), SDS (D) and sorbitol (E), concentrations and at distinct pHs (F). Growth was determined as the area under the curve (AUC) after 7 days of incubation at 24 °C. The two-way ANOVA, followed by Tukey’s test (p < 0.05, indicated as *), was performed to establish whether significant differences existed between the knockout mutants and the wild type. Values represent the mean and standard error of the mean of three biological replicates. The experiment was repeated twice.
Figure 3. Growth of the wt strain A. carbonarius ITEM 5010 (■) and two ΔveA (△ represents ΔveA10b, and ▽ denotes ΔveA12a) in the presence of different H2O2 (A), CFW (B), NaCl (C), SDS (D) and sorbitol (E), concentrations and at distinct pHs (F). Growth was determined as the area under the curve (AUC) after 7 days of incubation at 24 °C. The two-way ANOVA, followed by Tukey’s test (p < 0.05, indicated as *), was performed to establish whether significant differences existed between the knockout mutants and the wild type. Values represent the mean and standard error of the mean of three biological replicates. The experiment was repeated twice.
Toxins 14 00745 g003
Figure 4. Competitiveness of the Δpks (A) and Δvea (B) knockout mutants against the mycotoxigenic wild-type (wt) strain A. carbonarius ITEM 5010 on potato dextrose broth (PDB). Competition assays were performed at the 1:0, 10:1, 1:1, 1:10 and 0:1 (wt:Δ) ratios. The expected values for the wt strain are indicated as gray dashed lines. The percentage of each strain was estimated by qPCR. Values are the mean of at least three biological replicates and error bars represent the standard error of the mean (SEM). OTA production during the competition assays of (C) the wt strain vs. the Δpks mutant, and the (D) wt vs. Δvea mutant. Values are normalized based on the OTA production by the wt strain.
Figure 4. Competitiveness of the Δpks (A) and Δvea (B) knockout mutants against the mycotoxigenic wild-type (wt) strain A. carbonarius ITEM 5010 on potato dextrose broth (PDB). Competition assays were performed at the 1:0, 10:1, 1:1, 1:10 and 0:1 (wt:Δ) ratios. The expected values for the wt strain are indicated as gray dashed lines. The percentage of each strain was estimated by qPCR. Values are the mean of at least three biological replicates and error bars represent the standard error of the mean (SEM). OTA production during the competition assays of (C) the wt strain vs. the Δpks mutant, and the (D) wt vs. Δvea mutant. Values are normalized based on the OTA production by the wt strain.
Toxins 14 00745 g004
Figure 5. OTA production (A,B) during the infection of grape berries. Co-inoculation of the mycotoxigenic wil-type (wt) strain A. carbonarius ITEM 5010 and the Δpks mutant (A,C), and the wt and ΔveA (B,D). Competitive assays were performed at the 1:0, 10:1, 1:1, 1:10, and 0:1 ratios (wt:Δ). Grape berries were incubated at 28 °C for up to 10 days. Values are the mean of three biological replicates and error bars represent the standard error of the mean (SEM). Values are normalized based on the production of OTA by the wt strain.
Figure 5. OTA production (A,B) during the infection of grape berries. Co-inoculation of the mycotoxigenic wil-type (wt) strain A. carbonarius ITEM 5010 and the Δpks mutant (A,C), and the wt and ΔveA (B,D). Competitive assays were performed at the 1:0, 10:1, 1:1, 1:10, and 0:1 ratios (wt:Δ). Grape berries were incubated at 28 °C for up to 10 days. Values are the mean of three biological replicates and error bars represent the standard error of the mean (SEM). Values are normalized based on the production of OTA by the wt strain.
Toxins 14 00745 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Llobregat, B.; González-Candelas, L.; Ballester, A.-R. Ochratoxin A Defective Aspergillus carbonarius Mutants as Potential Biocontrol Agents. Toxins 2022, 14, 745. https://doi.org/10.3390/toxins14110745

AMA Style

Llobregat B, González-Candelas L, Ballester A-R. Ochratoxin A Defective Aspergillus carbonarius Mutants as Potential Biocontrol Agents. Toxins. 2022; 14(11):745. https://doi.org/10.3390/toxins14110745

Chicago/Turabian Style

Llobregat, Belén, Luis González-Candelas, and Ana-Rosa Ballester. 2022. "Ochratoxin A Defective Aspergillus carbonarius Mutants as Potential Biocontrol Agents" Toxins 14, no. 11: 745. https://doi.org/10.3390/toxins14110745

APA Style

Llobregat, B., González-Candelas, L., & Ballester, A. -R. (2022). Ochratoxin A Defective Aspergillus carbonarius Mutants as Potential Biocontrol Agents. Toxins, 14(11), 745. https://doi.org/10.3390/toxins14110745

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