Next Article in Journal
Natural Occurrence of Ochratoxin A in Blood and Milk Samples from Jennies and Their Foals after Delivery
Next Article in Special Issue
Efficient Aflatoxin B1 Sequestration by Yeast Cell Wall Extract and Hydrated Sodium Calcium Aluminosilicate Evaluated Using a Multimodal In-Vitro and Ex-Vivo Methodology
Previous Article in Journal
Aflatoxin B1 and Sterigmatocystin Binding Potential of Lactobacilli
Previous Article in Special Issue
Aflatoxin Reduction in Maize by Industrial-Scale Cleaning Solutions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development of PCR, LAMP and qPCR Assays for the Detection of Aflatoxigenic Strains of Aspergillus flavus and A. parasiticus in Hazelnut

by
Sara Franco Ortega
1,
Ilenia Siciliano
1,
Simona Prencipe
2,
Maria Lodovica Gullino
1,2 and
Davide Spadaro
1,2,*
1
Centre of Competence for the Innovation in the Agro-Environmental Sector—AGROINNOVA, University of Turin, via Paolo Braccini 2, I-10095 Grugliasco, Italy
2
Department of Agricultural, Forestry and Food Sciences (DiSAFA), University of Torino, via Paolo Braccini 2, I-10095 Grugliasco, Italy
*
Author to whom correspondence should be addressed.
Submission received: 30 September 2020 / Revised: 23 November 2020 / Accepted: 26 November 2020 / Published: 30 November 2020
(This article belongs to the Collection Aflatoxins)

Abstract

:
Aspergillus flavus and A. parasiticus are two species able to produce aflatoxins in foodstuffs, and in particular in hazelnuts, at harvest and during postharvest phase. As not all the strains of these species are aflatoxin producers, it is necessary to develop techniques that can detect aflatoxigenic from not aflatoxigenic strains. Two assays, a LAMP (loop-mediated isothermal amplification) and a real time PCR with TaqMan® probe were designed and validated in terms of specificity, sensitivity, reproducibility, and repeatability. The capability of the strains to produce aflatoxins was measured in vitro and both assays showed to be specific for the aflatoxigenic strains of A. flavus and A. parasiticus. The limit of detection of the LAMP assay was 100–999 picograms of DNA, while the qPCR detected 160 femtograms of DNA in hazelnuts. Both techniques were validated using artificially inoculated hazelnuts and naturally infected hazelnuts. The qPCR was able to detect as few as eight cells of aflatoxigenic Aspergillus in naturally infected hazelnut. The combination of the LAMP assay, which can be performed in less than an hour, as screening method, with the high sensitivity of the qPCR, as confirmation assay, is able to detect aflatoxigenic strains already in field, helping to preserve the food safety of hazelnuts.
Key Contribution: A LAMP and a qPCR assays were designed, optimized and validated in order to detect aflatoxigenic Aspergillus sp. with high specificity and sensitivity.

1. Introduction

Aflatoxins (AFs) are toxic polyketide secondary metabolites with low molecular weight and hepatocarcinogenic activity. These compounds are poorly soluble in water and highly soluble in moderately polar solvents, such as chloroform, methanol, or dimethylsulfoxide. Regarding their chemical structure, AFs can be divided into two groups identified based on the fluorescence emitted under ultraviolet light: type B (blue fluorescence) with a difuro-coumaro-cyclopentenone structure; type G (green fluorescence) with a six-member lactone ring replacing the cyclopentenone. Eighteen different types of AFs are currently identified [1], with the most important and common aflatoxins in nature being aflatoxin B1 (AFB1), B2 (AFB2), G1 (AFG1), and G2 (AFG2).
Aspergillus section Flavi includes 33 species and most of them are aflatoxin producers, such as Aspergillus flavus Link, A. parasiticus Speare, and A. toxicarius Murak [2,3,4]. However, non-aflatoxigenic species have also been identified, such as A. oryzae, A. sojae, and A. tamarii. Some of the non-aflatoxigenic species, such as A. oryzae and A. sojae, are widely used in the oriental food fermentation process, being considered as the domesticated species of A. flavus and A. parasiticus [5,6]. A. flavus is the main producer of AFB1, AFB2, and cyclopiazonic acid (CPA) [7,8] while the closely related species A. parasiticus and A. toxicarius are important producers of AFB1, AFB2, AFG1 and AFG2. Though, not all A. flavus strains are aflatoxigenic, depending on the functionality or not of the entire aflatoxin biosynthesis gene cluster [9].
Hazelnut (Corylus avelana L.) is an important product for human nutrition that has beneficial effects on human health by reducing cholesterol and preventing of chronic diseases, such as cancers or cardiovascular complications [10], due to the high concentration of monounsaturated fatty acids, especially oleic acid, and polyunsaturated fatty acids [11]. Italy is the second biggest world producer country of hazelnuts after Turkey. Production of hazelnuts in Italy was 132,699 tons in 2018 [12], representing 13% of the worldwide production [13].
Hazelnuts can be contaminated by different aflatoxigenic molds species at harvest, during storage and processing leading to significant economic losses and possible aflatoxins contamination. The harvest time and method significantly impact the presence of Aspergillus spp. on hazelnuts and therefore in the AFs content of the product. In Turkey, hazelnuts can be directly collected from the tree, but it is a time-consuming procedure and the most common harvest technique is the collection of the mature hazelnuts after dropping to the ground, which usually increases the mold contamination on the surface [14]. Traditionally, after harvesting, hazelnuts are sun-dried to decrease the water activity and to avoid mold growth. However, drying during rainy periods may also enhance the risk of aflatoxin contamination [14]. Hazelnuts storage before drying is another relevant factor, as the common practice of using nylon bags, which increase the humidity, moisten the hazelnuts and boost the susceptibility to fungal growth [14].
During postharvest, Aspergillus spp. and Penicillium spp., which are present at low levels in preharvest, become a threat in nuts, including hazelnuts, due to their adaptation to the storage conditions [14,15]. During postharvest storage, dried hazelnuts require low relative humidity and low temperature (5 °C and 65.5% RH), but when these conditions are not met, the risk of fungal contamination of the kernels and aflatoxin content can increase again [13,15,16,17,18].
Effective controls and accurate and efficient diagnosis should be developed to monitor and control Aspergillus spp. contaminations. The prediction of risk levels associated with mycotoxins in food and the compliance with legal limits and health regulations could be improved by the development of a Decision Support Systems (DSS) [19]. In this perspective, assessing the presence of mycotoxins and quantifying them through the different stages of production is essential. The EU has recently established the maximum level of AFs in nuts, including almonds, pistachios, and apricot kernels and hazelnuts [20]. The maximum level of AFB1 and total AFs in hazelnuts subjected to physical treatments is 8 µg/kg and 15 µg/kg, respectively, whilst it is 5 µg/kg for AFB1 and 10 µg/kg for total AFs in hazelnuts for human consumption. In the USA, the FDA has set 20 ng/g as the maximum level of AFs for nuts, whilst Japan has established the same levels of the EU for imported nuts used for direct consumption or as food ingredients [21,22].
The detection of Aspergillus spp. on food has been traditionally based on plate count and culturing methods studying macro- and micro-morphological characteristics. Traditional detection techniques are time-consuming and require well-trained staff. Molecular techniques overcome the risk of misidentification. Many attempts to develop the most specific and sensitive techniques, able to detect and quantify Aspergillus spp., have been performed. PCR-based assays have been developed to distinguish between aflatoxin producer and non-producer Aspergillus spp. based on the detection of several genes of the aflatoxin biosynthesis gene cluster, by using RNA-based assays, end-point PCR or PCR variants [23,24,25,26,27,28,29,30]. The assays developed included multiplex PCR [31], quantitative real time (qPCR) for the Aspergillus genus, Aspergillus section Flavi spp. [6,32,33,34], specific for A. flavus [35,36,37], and assays able to distinguish between A. flavus and A. parasiticus [38]. The loop-mediated isothermal amplification (LAMP) assay [39,40] has been confirmed as a robust method to detect plant pathogens as viruses, bacteria and fungi [41,42,43]. The high specificity of these assays is based on the use of six primers that recognize eight genetic regions in order to amplify a target DNA. More recently, LAMP assays have been designed to differentiate between ochratoxigenic Aspergillus spp. or aflatoxigenic and non-aflatoxigenic, or fumonisin producer strains of A. flavus [44,45,46,47,48,49]. In addition, a novel LAMP assay designed by Niessen et al. [50] was developed on the nor1 gene for the detection of aflatoxin producers within Aspergillus section Flavi, showing high sensitivity and specificity. This methodology is relatively rapid and may compete with the traditional identification because of its specificity and sensitivity. The detection of the LAMP products can be based on the use of indicators such as calcein or neutral red which causes a color change of the reaction from yellow (negative) to pink (positive) [50]. However, colorimetric changes, such as the Niessen et al. [50] assay, can be difficult to interpret, especially for color-blind people or when the light conditions are not optimal (as in the field or inside storage areas).
As the LAMP assay is less sensitive to inhibitors, crude DNA extraction methods can be associated, which can speed up the DNA extraction from hours to minutes and in combination with portable devices, such as the Genie II and III (Optigene, Horsham, UK), which measure the fluorescence of the LAMP products, this assay can be quickly transferred to the food sector.
On the other hand, quantitative PCR is a robust diagnostic tool due to its higher sensitivity which is critical to detect samples with low target concentrations. Developing and validating these two methods to detect aflatoxigenic molds is important to avoid contamination of AFs in the food processing sector. In this study, we developed three different DNA-based methods, an end-point PCR, a LAMP assay, and a qPCR assay, to detect the presence of two of the most common aflatoxigenic species belonging to Aspergillus section Flavi that are important contaminants in hazelnuts, A. flavus and A. parasiticus. The combination of these analyses in different scenarios may help in obtaining a successful reduction in aflatoxin contamination in the tree nuts industry. Previously, LAMP assays were developed to detect aflatoxigenic molds on nuts, without chemical confirmation on the strains. In this work, we developed LAMP and qPCR assays for the detection of aflatoxigenic strains of A. flavus and A. parasiticus on hazelnuts. The results were verified by HPLC-FLD method for aflatoxin production both in vitro and in vivo, to evaluate the matrix effect [51].

2. Results

2.1. Aflatoxin Analysis In Vitro

Ninety strains of Aspergillus section Flavi isolated from hazelnuts (50 strains) (Table 1) and chestnuts (40 strains) (Table 2) and 7 CBS reference strains (Table S1) were analyzed for aflatoxin production in vitro. The identification of hazelnuts isolates was performed by sequencing the calmodulin gene which was later used for inferring a maximum likelihood tree (Figure S1). This gene identified 42 of the 50 isolates from hazelnuts. As eight isolates presented a low sequence quality, they were not included in the phylogenetic tree. However, visual observation identified these eight isolates as part of Aspergillus section Flavi. Most strains (39 out of 50) were identified as A. flavus, and three strains (FS8, BS1 and BS8) as A. parasiticus/A. sojae. Based on micro- and macro-morphological observations, the strains FS8 and BS1 were identified as A. parasiticus, while the strain BS8 was A. sojae.
AFB1 was the most frequently produced mycotoxin, followed by AFB2. In chestnuts strains, the AFB1 was produced by 13 out of 30 (43.3%) of A. flavus and by the A. parasiticus and A. toxicarius strains, whilst the AFB1 was produced by 28 out of 50 hazelnut strains (56.0%, corresponding to four strains that were not identified by the calmodulin gene, 22 A. flavus and two A. parasiticus) and most of these 28 strains were also able to produce the AFB2 (44.0%, corresponding to three strains that were not identified with the calmodulin gene, 17 A. flavus and two A. parasiticus). AFG1 and AFG2 were produced by A. parasiticus and A. toxicarius strains. One strain of A. sojae from hazelnut (BS8) did not produce any aflatoxin [4].
Regarding the leaves, bracts and kernel isolates, out of the 26 strains from hazelnut leaves, 13 strains produced AFB1 (A. flavus and A. parasiticus), and 10 of them also produced AFB2, whilst 11 and 10 strains from bracts produced AFB1 and AFB2, respectively. However, only three and two out of eight strains from kernels produced AFB1 and AFB2, respectively. AFG1 and AFG2 were only detected in one strain from leaves (FS8) and two from bracts (BS1 and BS6).
The highest amount of mycotoxin was produced by the strains identified as A. toxicarius and A. parasiticus isolated from chestnuts and by the strains FS8 and BS1 (A. parasiticus) and strains BS6 and FV5 (A. flavus) isolated from hazelnuts.
The method used for the quantification was validated for linearity, limit of detection (LOD), limit of quantification (LOQ), and matrix effect (ME). The analytical method was tested in a range of concentrations from 0.1 up to 20 ppb. After extraction, all the samples were diluted to remain in the linearity range. ME was calculated for the four AFs in yeast extract-sucrose (YES) culture broth and the result was 59% for AFB1, 57% for AFB2, 74% for AFG1, and 57% for AFG2. In order to minimize the susceptibility related to the matrix effect, the quantification of the analytes was obtained through the calibration curves prepared in the matrix. Calculated LODs were 0.1, 0.2, 0.4, and 0.2 ppb and calculated LOQs were 0.5, 0.7, 1.2, and 0.8 ppb for AFB1, AFB2, AFG1, and AFG2, respectively.

2.2. Identification of the Region for LAMP and qPCR Primer Design. PCR Validation

The presence or absence (deletion) of a PCR product of approximately 308 bp obtained with the primers AfF and AfR, listed in Table 3, was checked by means of agarose gel. The primers aligned with the aflatoxin biosynthetic cluster of different species (Figure S2) corresponding with the end of the last exon (AfF) and the beginning of the last intron (AfR) of the aflatoxin efflux pump (aflT) gene, as shown in the Figure S3. The product was present in 48 strains and reference strains. The product was not present in A. oryzae var. oryzae CBS 46691, and A. flavus var. columnaris (CBS 485.65) (Figure S4). The amplification was unclear for the strain BS1. Comparing with the aflatoxin production, BS6 (A flavus), BS9 (A. flavus), AFLX1 (A. flavus), AFLX2 (A. flavus), AFSP2 (A. flavus), AFSP3 (A. flavus) did not present the band, but were able to produce at least one of the four aflatoxins studied. On the contrary, FP8 (A. flavus), BV10 (A. flavus), BS5 (A. flavus), NP1 (A. flavus), NP6 (A. flavus), AFCAL7 (A. flavus), AFLX7 (A. flavus) and AFLX10 (A. flavus) present the band but no aflatoxin was detected. The sensitivity of the technique was validated using a standard curve with the DNA of the A. flavus FS3 (Figure S5) starting from 110 ng/µL and reaching an LOD of 111 fg/µL.

2.3. LAMP Assay Validation

The six LAMP primers, designed on a region covering 186 bp (product size calculated using the A. parasiticus sequence of the aflatoxigenic biosynthetic cluster, GenBank accession number: AY371490) of the aflatoxin biosynthetic gene cluster (Figure S2), are listed in Table 3. The primers aligned with the end of the last exon and the beginning of the last intron of the aflT gene as indicated in the Figure S3. Positive results of the LAMP assays were obtained only for the strains with the presence of a 308 bp PCR product checked by means of agarose gel. BLASTn analysis of F3, loopF, loopB, and FIP primers showed similarity at 100% with A. flavus NRRL 3357 (GenBank accession number CP0044620 and XM 002379911) that causes aflatoxin contamination on food and feed [48] and with the strains A. flavus AF36, A. flavus AF13, A. flavus BN008, A. flavus AF70, and other A. parasiticus and A. flavus strains (GenBank accession numbers: AY510455, AY510451, AY510452, AY510453, AY371490, AF515601, respectively), as well, they showed 100% similarity with Aspergillus pseudotamarii (GenBank accession number XM_032052844.1) and A. sojae (CP035528.1). BIP and B3 primers showed 100% similarity with A. flavus NRRL 3357, A. flavus AF36, A. flavus AF13, A. flavus AF70, and A. parasiticus and A. flavus (GenBank accession numbers: CP0044620.1, AY510455, AY510451, AY510452, AY510453, CP051029.1, AF515601, respectively) and approximately at 90% with A. flavus BN008 (AY510452.1). AF36 strain was reported as non-aflatoxigenic whilst the strain BN008 was reported to produce both B and G AFs by Ehrlich et al. [49], which is a common feature of A. parasiticus. According to the online software, the primers showed appropriate structural features, which avoid potential secondary structures. In order to verify the specificity of the LAMP primers, the PCR product obtained with the primers AfF and AfR using the DNA of FP8, BV10 and BS5 was sequenced, and the LAMP primers were found in the sequence, with all the primers matching at 100%.The specificity of the manually design assay was validated according to a panel of inclusion and exclusion of 90 Aspergillus spp. isolated from hazelnuts and chestnuts and seven CBS references (Table 1 and Table 2 and Table S1). Other common pathogens of hazelnuts and chestnuts were included in the panel (Table S2). Considering the 97 Aspergillus spp., and the in vitro mycotoxin production, the LAMP assay produced 42 real positives including the species A. subolivaceus CBS 501.65, A. thomii CBS 120.51, A. arachidicola CBS 117610, A. minisclerotigenes CBS 117635 and A. parvisclerotigenus CBS 121.62, 2 false positives, 5 false negatives and 48 real negatives. Therefore, no LAMP amplification was found in most of the non-aflatoxin producer strains and nor was it found on other fungi commonly found in hazelnuts or chestnuts. However, the LAMP was not able to detect the strains BS3 (A. flavus), AFLX1 (A. flavus), AFLX2 (A. flavus), AFSP2 (A. flavus) and AFSP3 (A. flavus) which produced AFs in vitro, while by using the LAMP assay a positive amplification was obtained using the DNA of two A. flavus AFCAL7 and AFLX7 which were not able to produce aflatoxins in vitro.
The time to positive (Tp) ranged from 15 min 30 s to 37 min 17 s whilst the melting temperature ranged from 83.83 to 92.3.67 °C (Table 1 and Table 2). The sensitivity was tested against 10-fold serial dilutions of the DNA of the strains FS7 (A. flavus) and AFLX6 (A. parasiticus), showing a limit of detection of 100-999 pg/µL of DNA for both species (Table 4, Figure S6) with a Tp of the aflatoxigenic strains obtained in the three independent assays in a repetitive and reproducible way.

2.4. qPCR Validation

The primers and TaqMan® probe sequences, which targeted 74 bp (product size calculated according to the A. parasiticus sequence of the aflatoxigenic biosynthetic cluster, GenBank accession number: AY371490) are listed in Table 3 and Figure S2. The primers aligned with 100% similarity against the sequenced product of the PCR performed with the AfF and AfR primers and the DNA of FP8, BV10 and FV5 and produced an expected amplicon of 74 bp in the three case. Regarding the localization in the aflT gene, the primers and probe were designed at the end of the last exon of the aflT gene (Figure S3). BLASTn analysis of the primers forward, reverse and the probe, showed 100% similarity with A. flavus NRRL3357, A. flavus AF36, A. flavus AF13, A. flavus BN008, A. flavus AF70 and other A. parasiticus and A. flavus (GenBank accession numbers: XM002379911.1, AY510455, AY510451, AY510452, AY510453, AY371490, AF515601, respectively).
SYBR™ Green assay was initially used to check the specificity of the primers and the melting temperature of the qPCR product. The Ct of the positive results ranged from 11.58 to 41.14 whilst the melting temperature of the positive results was between 78.74 and 80.97 °C. Three replications of each reaction were carried out, with some of the technical replicates showing positive amplification and, many non-aflatoxigenic strains were amplified, such as some non-aflatoxigenic, i.e., A. flavus strain AFLX9 and AFLX10. Therefore, in order to increase the specificity of the qPCR, a TaqMan® probe was designed. In total, the SYBR green assay produced 33 real positives (aflatoxigenic strains), 24 real negatives (non-aflatoxigenic strains), 27 false positives and 13 false negatives. The sensitivity of the SYBR Green tested through a standard curve with 10-fold serial dilutions of the DNA from the A. flavus FS7, showed a limit of detection of 169 fg/µL of DNA (with only two out of three replicates detected at this dilution). The standard curve equation was y = −1.1313x + 25.624; with R2 = 0.967 using the ln of the DNA concentration, as shown in the Figure S7. The specificity assay for the qPCR with TaqMan® probe showed 42 real positives (aflatoxigenic strains) including 3 out of 7 reference strains amplified (A. subolivaceus CBS 501.65, A. thomii CBS 120.51 and A. minisclerotigenes CBS 117635), whilst A arachidicola CBS 117610, A. parvisclerotigenus CBS 121.62, A. oryzae var. oryzae CBS 466.31 and A. flavus var. columnaris CBS 485.65 reference strains were not detected with this assay. In total, 8 false positives, 5 false negatives and 42 real negatives were detected. The Ct with the TaqMan® probe ranged from 15.6 to 40.37. No amplification was obtained with the DNA of other fungi commonly found in hazelnuts (Table S2). However, the TaqMan® qPCR was not able to detect the aflatoxigenic strains AFLX1, AFLX2 and AFSP2, whilst it gave positive amplification using the DNA of the A. flavus FP8, BV10, BS5, NP1, NP2, NP3, AFCAL7 and AFLX7 which were unable to produce aflatoxins in vitro.
Using the same DNA of the panel of inclusion-exclusion adjusted between 1 and 50 ng/µL, previously used for the LAMP assay, the Ct ranged from 19.80 to 30.10 (Table 1 and Table 2 and Table S1). The sensitivity of this assay was tested with the serial dilutions of the DNA of A. flavus FS7 and A. parasiticus AFLX6, starting with 16.9 ng/µL and 1.83 ng/µL, respectively. The qPCR with TaqMan® probe showed an LOD of 16.9 fg/µL (only two replicated out of three) and 183 fg/µL (two replicates out of three).
The standard curve for A. flavus FS7 was y = −1.4471x + 24.844, R² = 0.9937, slope of 3.5394 and efficiency of 91.66%, whilst the standard curve using the ln of the DNA concentration of AFLX6 was y = −1.49x + 24.124, R2 of 0.987 with a slope of 3.42, and an efficiency of 96.06% (Figure S7).

2.5. In Vivo Validation

The sensitivity of both techniques (LAMP and qPCR with TaqMan® probe) was evaluated using artificially inoculated hazelnuts in order to test the reliability of these techniques on hazelnuts as food matrix. Three different strains, previously characterized by Prencipe et al. [53] corresponding with one non-aflatoxigenic A. flavus (AFSP1), one aflatoxigenic A. flavus (AFSP4) and one aflatoxigenic A. parasiticus (AFLX6) were used to inoculate hazelnuts at three different concentrations. No amplification was shown in any dilution with any type of DNA extraction with the non-aflatoxigenic A. flavus isolate AFSP1.
Regarding the LAMP results, the type of DNA extraction (crude alkaline extraction or commercial kit) influenced the number of positive replicates at each concentration. With the crude DNA extraction, 3/9, 1/9 and 3/9 hazelnuts inoculated with the aflatoxigenic strain AFSP4 were detected at the concentrations of 3.03 × 105, 3.03 × 106, and 3.03 × 107 spores/mL, respectively. The Tp ranged from 32 min 09 s to 58 min 05 s and the melting temperature from 88.76 to 89.21 °C. The approximate number of cells in the crude extract is reported in Table 5 and was calculated according to the standard curve of the TaqMan® qPCR. The number of cells calculated with the TaqMan® qPCR in each mL of crude extraction corresponded to 1.7 × 105 cells, 1.1 × 104 cells and 6.3 × 105 cells, respectively. Instead, hazelnuts artificially inoculated with the strain AFLX6 were only detected in one and two out of nine replicates with an inoculum of 1.79 × 105 and 1.79 × 106 spores/mL, respectively. The number of cells calculated with the TaqMan® qPCR was 5.6 × 105 and 8.4 × 104, respectively. Hazelnuts inoculated with 1.79 × 107 spores/mL were not detected with the LAMP assay. The Tp ranged from 38 min 06 s to 52 min 46 s and the melting temperature from 88.17 to 89.23 °C.
In all cases, the cytochrome oxidase I (COX) LAMP assay was used as internal control and produced positive results in all the replicates with Tp ranging from 4 min 29 s to 15 min 40 s and melting temperature from 84.92 to 85.36 °C.
The number of positive LAMP results increased significantly when the DNA of the hazelnuts was extracted using a commercial kit. This commercial DNA extraction kit was also used to check the sensitivity of the qPCR with TaqMan® assay in vivo. Regarding the hazelnuts inoculated with the AFSP4 strain, nine out of nine replicates were detected with the LAMP assay. In this case, the Tp ranged from 28 min 16 s to 50 min 16 s and the melting temperature ranged from 87.89 to 88.78 °C. The TaqMan® qPCR produced repetitive results for all the concentrations tested with Ct ranging from 25.65 to 33.98 for the hazelnuts inoculated with 3.03 × 105 spores/mL; 30.63 to 32.49 for the hazelnuts inoculated with 3.03 × 106 spores/mL and 24.49 to 28.66 for the hazelnuts inoculated with 3.03 × 107 spores/mL. The Ct was afterwards used to calculate the approximate number of cells in each µL of the DNA extraction by means of a standard curve. The lowest number of cells detected was 53 spores/mL (samples inoculated with 3.03 × 106 spores/mL) and the maximum 6.2 × 104 spores/mL for samples inoculated with 3.03 × 107 spores/mL.
Regarding the hazelnuts inoculated with the strain AFLX6 and extracted with a commercial kit, they were detected in most replicates, 6/9 at the concentration of 1.79 × 105 spores/mL, 8/9 at the concentration of 1.79 × 106 spores/mL and 9/9 at the concentration of 1.79 × 107 spores/mL, with Tp ranging from 30 min 24 s to 45 min 42 s and melting temperature ranging from 88.03 and 88.78 °C. The TaqMan® qPCR produced repetitive results in all the concentrations of this strain tested. The Ct ranged from 26.90 to 39.89 for the hazelnuts inoculated with 1.79 × 105 spores/mL, 29.91 to 37.10 for the hazelnuts inoculated with 1.79 × 106 spores/mL, 25.63 to 31.11 for the hazelnuts inoculated with 1.79 × 107 spores/mL. The lowest number of cells detected was 60 spores/mL and the maximum 1.9 × 105 spores/mL for samples inoculated with 1.79 × 107 spores/mL.

2.6. Naturally Contaminated Hazelnuts

The presence of Aspergillus spp. in the hazelnuts was checked by plating serial dilution into Aspergillus flavus and parasiticus agar (AFPA) medium. The presence of orange-to yellow colonies in the dilution 1 × 105 confirmed the natural contamination of the lot.
Nineteen out of 50 hazelnuts extracted with a commercial kit showed positive amplification using the qPCR with the TaqMan® probe. The Ct ranged from 33.79 to 39.99 which corresponded with 0.0003 to 0.024 ng of DNA and between 8 spores/µL and 612 spores/µL (Table 6). Though, no amplification was obtained with the LAMP assay with any type of DNA extraction, due to a lower sensibility of this technique compared with the qPCR with TaqMan® probe.

2.7. Comparison between PCR, TaqMan® qPCR, LAMP Assay and Aflatoxin Production

A positive correlation was found between the presence of aflT, detected by end-point PCR, with the production of AFs. Out of 97 strains tested in this study, from hazelnuts, chestnuts and CBS collection, only six strains did not present the band, but at least they produced one of the AFs (usually at least AFB1). These six isolates were all A. flavus (BS6, BS9, AFLX1, AFLX2, AFSP2 and AFSP3). These results suggested that in most of the strains, this region of the aflT was essential for the production of AFs, and it was exploited in order to design the LAMP assay and the qPCR with TaqMan®. However, there were also eight strains that presented the band but did not produce any AFs in vitro all identified as A. flavus. The diagnostic specificity and sensitivity were calculated for this assay, considering real positives—aflatoxin producers of at least one of the four aflatoxins investigated— as well as real negatives—non-aflatoxin producers in vitro. The PCR showed +LR of 4.85 and -LR of 0.16 meaning that approximate 1 in 1.2 positive results, are aflatoxigenic strains, whilst, 1 in 1.2 negative results are non-aflatoxigenic strains. The diagnostic sensitivity was 87.2% and the diagnostic specificity was 82.0% (Table 7).
A comparison between the number of real positives, real negatives, false positives and false negatives was carried out also for the LAMP and the qPCR (Table 7). With DNA from pure cultures of 97 strains (Table 1 and Table 2 and Table S1), the LAMP assay identified 42 real positives, 48 real negatives, 5 false negatives (the LAMP did not produce amplification, but the strain produced at least one aflatoxin in vitro) and 2 false positives (the LAMP was positive, but the strains did not produce any aflatoxin in vitro). The LAMP assay showed a diagnostic sensitivity of 89.4% and diagnostic specificity of 96.0%, with +LR of 22 and –LR 0.11 meaning that approximately 1 in 1.0 strains with a positive LAMP result are real aflatoxigenic strains, whilst approximately 1 in 1.1, with negative LAMP result are non-aflatoxigenic strains.
The qPCR assay with SYBR Green produced 33 real positives, 24 real negatives, 27 false positives and 13 false negatives that related to the +LR of 1.36 (1 in 1.8 positive results are aflatoxigenic strains) and –LR of 0.6 (1 in 1.5 negative results are non-aflatoxigenic strains). The diagnostic sensitivity of this assay was 71.7% and the specificity was 47.1%. All these values increased with the TaqMan® qPCR (42 real positives, 42 real negatives, 8 false positives and 5 false negatives). The TaqMan® qPCR assay showed diagnostics sensitivity of 89.4% and diagnostic specificity of 84.0% with +LR of 8.94 and –LR 0.12 meaning that approximately 1 in 1.1 strains with a positive TaqMan® qPCR positive are real aflatoxigenic strains and 1 in 1.2 are non-aflatoxigenic strains when the TaqMan® qPCR gives a negative result (Table 7).
The TaqMan® qPCR assay was not able to detect A. parvisclerotigenus and A. arachidicola, but the LAMP assay showed a quick amplification using the DNA of these two strains.
Using the data from the artificially inoculated hazelnuts, and considering each biological replicate as an individual sample (81 samples in total), the LAMP using the crude extraction method has a diagnostic sensitivity of 18.5% and diagnostic specificity of 100%, with +LR of 999999 and –LR 0.81, meaning that approximately 1 in 1.0 strains with a positive LAMP result are real aflatoxigenic strains, whilst approximately 1 in 2.6 with negative LAMP result are real non-aflatoxigenic strains. Regarding the results of the LAMP with the DNA of the hazelnuts extracted using a commercial kit, the diagnostic sensitivity and specificity were 92.6% and 100%, respectively, with +LR of 999999 and –LR 0.07 meaning that approximately 1 in 1.0 strains with a positive LAMP result are real aflatoxigenic strains, whilst approximately 1 in 1.1 with negative LAMP results are true non-aflatoxigenic strains. Instead, the TaqMan® qPCR assay showed diagnostics sensitivity and specificity of 1, with +LR of 999999 and –LR 0.00 meaning that approximately 1 in 1.0 strains with a positive TaqMan® qPCR are real aflatoxigenic and 1 out of 1 with a negative TaqMan® qPCR results are non-aflatoxigenic strains.
In connection with the deletion of the aflatoxigenic cluster, three strains (A. flavus FP8, BV10 and BS5) showed a positive amplification with the end-point PCR, and the qPCR assay, but no positive result was found with the LAMP.
The A. flavus strains AFLX1, AFLX2 and AFSP2 reported as aflatoxigenic, were not detected by any technique (LAMP or qPCR), but A. flavus AFSP3 was only detected by the TaqMan® qPCR assay and not using the LAMP assay. None of these strains amplified with the end-point PCR primers.

3. Discussion

Presence of aflatoxigenic molds in food and feed is a critical threat to consumer safety and to the industry. The most effective management actions focus on the reduction in aflatoxin contamination in food, while maintaining the food sanitary standards [54].
Riba et al. [55], reported that 51.8% of the Aspergillus section Flavi strains found in roasted hazelnuts were aflatoxigenic, confirmed in our data which showed that 48% of the hazelnut strains were able to produce AFB1 in vitro. In addition, the AFB2 was also produced in most of the strains that produced the AFB1.
Among Aspergillus section Flavi, A. flavus has been reported more predominant than A. parasiticus in hazelnuts [14]. Traditionally, the Aspergillus species identification in food matrices is based on plating methods on selective media that are time-consuming and labor-intensive [56]. A. oryzae, a non-aflatoxigenic species, is the most phylogenetically related species to A. flavus, as A. sojae and A. nomius to A. parasiticus and A. toxicarius. For example, A. parasiticus and A. sojae can be identified accordingly to the color of colonies, the conidial diameter or using a bleomycin media [57]. However, these differences are not always clear due to the changes in the morphological features caused by the environment, and an accurate identification needs a lot of experience to differentiate between the species. In addition, it has been reported that A. flavus and A. parasiticus after long culturing may change their morphological features to appear like A. oryzae and A. sojae and A. nomius, respectively [58]. Currently, sequencing of barcoding genes such as ITS (internal transcribed spacer), beta-tubulin gene, RNA polymerase II gene (RPB2) or calmodulin gene [4,53,59,60] can resolve almost all the species of Aspergillus section Flavi, however, they cannot distinguish between A. flavus and the phylogenetically related and non-aflatoxigenic species A. oryzae, or between A. parasiticus and A. sojae [4]. Although, in the food-producing industry, it is more critical to determinate the aflatoxigenic potential of the strain than identify the species and this requires a combination of multiple procedures such as categorizing the species in function of the micro and macro morphological features, barcoding and the quantification of the secondary metabolites [29,53,59].
Considering the time necessary for culturing and the challenges associated with the morphological identification, other techniques have been used to determine the presence of aflatoxigenic fungi in foodstuffs. PCR-based assays are widely useful in laboratories due to their specificity, sensitivity and accurate detection properties. They are essential to identify non-predominant species without the pre-enrichment step of growing in a selective medium. However, due to the high degree of DNA relatedness between A. flavus and A. oryzae and A. sojae, A. nomius and A. parasiticus [5,61], the design of PCR-based methods for the identification of different Aspergillus spp. can be a challenge. Although barcoding genes (calmodulin or beta tubulin) are widely used for phylogenetic purposes, they do not help in distinguishing between aflatoxigenic and non-aflatoxigenic strains. Other genes involved in the aflatoxin biosynthetic pathway can be exploited, such as the omt-A, or the aflR, that have been used by some authors to design end-point PCR or TaqMan® qPCR to detect and quantify aflatoxin-producing Aspergillus spp. [23,33,36,62]. LAMP assays have been designed in other genes such as alpha amylase 1 (amy1) to detect A. flavus or ATP citrate lyase subunit 1 (acl1) to detect A. nomius and A. parasiticus [44] in Brazil nuts, peanuts and green coffee beans, similarly to the qPCR assays. The aflatoxigenic biosynthetic gene cluster has been already exploited to design LAMP assay methods on genes such as aflP, aflD or nor1 to distinguish between aflatoxigenic and non-aflatoxigenic fungi in different food matrices [46,48,50].
Within the aflatoxin biosynthetic gene cluster, some studies reported the positive relationship between the lack of some regions with the inability to produce aflatoxins. The region norBcypA showed to present some deletions in the non-aflatoxigenic isolates [9], and the total lack of the aflatoxin cluster has been shown in some strains. Chang et al. [9] reported eight patterns of the cluster according to the lack of some regions. Four types have been described according to the deletions identified within norBcypA regions depending on the amino acids affected and the intergenic region. While four patterns have been characterized by presenting a deletion larger than 45 where the entire aflatoxin gene cluster was lost, even affecting downstream sugar utilization gene cluster (nadA-hexA-glcA-sugR).
Designing primers to obtain specific and sensitive diagnostic tests can be a hard task, and many challenges were faced in this study in order to obtain reliable tests which exploited the differences in the aflatoxigenic biosynthetic cluster between aflatoxigenic and non-aflatoxigenic strains. In this case, we profited from a region of the aflT gene [63], whose deletion was positively correlated with the absence of aflatoxin production in chestnuts and hazelnuts strains. This region allowed us to obtain three types of diagnostic assays that could be applied in different scenarios, as a quick test in the laboratory (end-point PCR), as a quick and user-friendly on-site test (LAMP) and for the quantification of the aflatoxigenic strains and confirmation of samples with low amount of inoculum (qPCR with TaqMan® probe). The latter two methodologies were validated according to European Mediterranean Plant Protection Organization (EPPO) PM 7/98 standards in terms of specificity, sensitivity, repeatability and reproducibility with DNA obtained from pure culture and from hazelnuts. For the LAMP assay, in order to obtain specific primers, instead of using software, primers were manually designed following the guidelines in the Eiken Genome Site (http://loopamp.eiken.co.jp/e/lamp/primer.html) in terms of distance between primer regions, GC content and secondary structure. Out of 97 strains used in this study to confirm the specificity of all the techniques, considering also the CBS reference stains, only six did not present the band and did produce AFs. However, three out of these strains gave a positive result with the TaqMan® qPCR but not with the LAMP. Instead, nine strains presented the band and did not produce any aflatoxin.
Quantification of the mycotoxins in the different matrixes usually requires analytical techniques to quantify these secondary metabolites. Mycotoxins are usually produced at low concentrations, in consequence, liquid chromatography associated with mass spectrometry or fluorescence detection has become a unique method for the quantification of mycotoxins due to the high sensitivity and reliability [64]. Differently from previously published assays [50], we confirmed the aflatoxigenic potential of all the strains used in this study with a robust and sensitive analytical method which showed an LOD of 0.1, 0.2, 0.4, and 0.2 ppb and an LOQ of 0.5, 0.7, 1.2, and 0.8 ppb for AFB1, AFB2, AFG1, and AFG2, respectively. The matrix effect of YES broth was also taken into consideration to calculate the LOQ.
Nonetheless, it is possible that AFs production was below the LOQ of the analytical method (although the method was sensitive) in the nine strains that present the PCR band but no AFs production. However, the most likely explanation, that confirms what is reported in literature, is that the in vitro assay could underrate the aflatoxigenic potential of strains. This idea can explain that some of the strains, which did not produce AFs in vitro, resulted in a positive LAMP or qPCR results as it is likely that they produce aflatoxins in vivo [53,65]. As the objective of this study was the design of techniques that can detect or quantify aflatoxigenic fungi, but not to quantify the level of AFs, HPLC analyses were used as a standard method to validate the molecular techniques. Despite, the end-point PCR assays showed a high specificity of 82.0% with an LOD of 110 fg/µL of DNA. In terms of sensitivity, the LAMP assay designed here, showed a limit of detection of 100–999 pg/µL of DNA, like the limit of detection previously described in other assays [66,67]. However, Niessen et al. [50] also developed a LAMP assay that focused on the detection of aflatoxigenic strains of Aspergillus section Flavi, with primers designed on the nor1 gene, with a lower detection limit of 9.03 pg of genomic DNA per reaction and based on the detection of the positive results by a colorimetric change or by a fluorescence dye. Although the LOD of our LAMP was higher than that of other assays, we obtained a more sensitive diagnostic test with a TaqMan® qPCR that detected 169 fg/µL of genomic DNA of A. flavus and 0.0003 ng/µL of DNA of hazelnuts naturally contaminated by aflatoxigenic strains, corresponding with eight spores/µL. The diagnostic sensitivity was 89.4% both for LAMP and qPCR.
Both techniques, TaqMan® qPCR and LAMP, showed high specificity using DNA from pure culture samples, with no amplification of any other fungi commonly found in hazelnuts. In both cases, the diagnostic specificity and sensitivity were higher than in other studies (96% for the LAMP and 84% for the TaqMan® qPCR), overcoming other assays such as the LAMP designed by Douksouna et al. [48], to detect aflatoxigenic molds in rice, with a specificity of 71.5% using pure DNA. In addition, both tests were validated in terms of repeatability and reproducibility following the EPPO PM7/98 validation standards. Melting temperatures of the SYBR (ranging between 79 and 80 °C) and the LAMP (ranging from 88 to 89 °C) assays were also included as key parameters for the interpretation of the results.
Although both LAMP and qPCR were validated to detect the most common Aspergillus species that are hazelnut contaminants (A. flavus and A. parasiticus), we have seen positive amplification with the DNA of other species such as A. toxicarius, A. subolivaceus, A. thomii, A. arachidicola, A. minisclerotigens and A. parvisclerotigenus. Currently, Aspergillus section Flavi, includes 15 species, 5 of them recently described [4], with four well supported lineages. The two main lineages are centered in A. flavus (including A. minisclerotigenes, A. subolivaceus, A. var. columnaris, A. thomii and A. oryzae) and centered in A. parasiticus (including A. toxicarius, A. arachidicola and A. sojae) [4]. Although in vivo validation of these assays were only performed with A. flavus and A. parasiticus spores, the LAMP and TaqMan® qPCR could be potentially used to detect the DNA of other species of the Aspergillus section Flavi.
It is important to highlight an important caveat in the specificity, as some of the LAMP primers showed 100% similarity in silico with non-aflatoxigenic strains such as. A. sojae or A. pseudotamarii. It was confirmed that the LAMP or TaqMan® qPCR did not amplify A. sojae (isolate BS8), however, we only confirmed that both qPCR and LAMP primers did not amplify the DNA of A. tamarii, which shares 90% similarity with A. pseudotamarii and belongs to the same lineage [4]. In addition, LAMP assays are less sensitive to inhibitors present in the matrix than other molecular techniques allowing the use of crude extraction based on the alkaline disruption of the membrane of the cells to extract the DNA [68]. Previous LAMP assays did not test the sensitivity of the technique with real samples, artificially or naturally infected samples [48] However, it is essential to evaluate the potential application of the assays in real samples. No inhibitions were shown using hazelnuts samples with any of the two assays. During the sensitivity assay using artificially inoculated hazelnuts, repetitive results were obtained using the commercial DNA extraction, within the replicates of the aflatoxigenic A. flavus and the A. parasiticus and no amplification was obtained for the non-aflatoxigenic A. flavus confirming the specificity of both techniques in vivo.
However, not all the LAMP tests gave positive amplification with the crude DNA extracts. This type of crude extraction dilutes the samples over 50,000 times, and it has been a suitable type of extraction for another type of matrix such as plant material, seeds or fruit [52,66,69]. However, in this case it did not identify samples with a low concentration of inoculum.
In some situations, the EPPO (European Mediterranean Plant Protection Organization) protocols suggest the combination of methods rather than individual tests [70]. The low admitted levels of aflatoxins in food in the European Union [71] require the measurement of mycotoxin concentration in the different foodstuff matrixes by HPLC as it might be possible to find lots with concentrations that exceed legal limits [72]. The risk to human health of foodstuffs contaminated with aflatoxins makes it essential to detect aflatoxigenic fungi, and therefore, to evaluate the uncertainly of each diagnostic test in terms of false positives and false negatives. In some situations, such as routine diagnosis, the speed of the analysis could be the most important parameter, so that a single technique, such as the LAMP, could be used as a screening method that will speed up further decisions. Then, despite the negative results of the LAMP on hazelnuts using a crude extraction method and a sensitivity of 100–999 pg/µL of DNA, the assay on portable devices such as the Genie II instrument could be used as an initial screening method per batch, avoiding misreading of color changes. We suggest that negative results should then be confirmed with the TaqMan® qPCR assay (with higher sensitivity of 160 fg/µL of DNA or 8 spores/µL using DNA from hazelnuts) or the LAMP using a more purified extraction method, such as a commercial DNA extraction, which is more sensitive. In addition, the use of the end-point PCR primers and the LAMP assay could be used on the DNA of pure cultures of the Aspergillus section Flavi spp. to confirm the aflatoxigenic potential in some hours or in less than 40 min, respectively. However, HPLC methods, as the one used here, remain a standard method for quantifying the level of aflatoxins.

4. Conclusions

In summary, the direct analysis of hazelnuts, instead of DNA extraction from pure cultures obtained 7 to 14 days after plating, will significantly reduce the time between sampling and result obtainment. This will help greatly in the decision-making behind keeping or discarding a lot of nuts in the food industry, due to possible aflatoxin contamination. The methods developed here, the LAMP assay or the TaqMan® qPCR assay, identified with high confidence the aflatoxigenic A. flavus, and A. parasiticus in both pure DNA and hazelnuts. After the validation following an international EPPO standard (PM7/98), both techniques may be used as complementary to each other to identify the presence of aflatoxigenic strains in food matrices. Positive samples could be confirmed by HPLC methods in order to precisely quantify the AFs levels in the food matrix.

5. Materials and Methods

5.1. Isolate Collection and DNA Extraction

Ninety Aspergillus spp. listed in Table 1 and Table 2 were used for the validation of both LAMP and qPCR assays. Fifty strains (Table 1) were isolated from hazelnut leaves (24), bracts (16) and nuts (8) and forty (Table 2) from different steps of the chestnut production chain [57]. Moreover, seven reference strains of Aspergillus spp. purchased from the Westerdijk Fungal Biodiversity Institute were included to validate the results of the LAMP and the qPCR assays (Table S1). These 7 species were selected as they were part of the two main lineages of Aspergillus section Flavi, the one centred in A. flavus and the other centered in A. parasiticus [4]. Single-spore cultures of isolates were grown in YES agar as described by Frisvad [73] for 7 days at room temperature according to the protocol described by Visagie et al. [74]. The mycelium of each isolate was collected from the agar plate and stored at −20 °C. Total genomic DNA was obtained using the E.Z.N.A. Fungal DNA mini kit (OMEGA Bio-Tek, Norcross, GA, USA), according to the manufacturer’s instructions and stored at −20 °C. DNA concentrations were measured using NanoDrop 2.0 (ThermoFisher, Wilmington, DE, USA) and adjusted to 1–50 ng/µL. The same genomic DNA was used to check the specificity of both LAMP and qPCR assays. Molecular identifications of Aspergillus section Flavi obtained from hazelnut trees were performed using the primers cmd5 (5′-CCGAGTACAAGGARGCCTTC-3′) and cmd6 (5′-CCGATRGAGGTCATRACGTGG-3′) [75] to amplify the calmodulin gene. The PCR reactions were performed in 20 µL with approximately 50 ng of genomic DNA, 2.5 mM deoxynucleotide triphosphates, 2.5 mM MgCl2, 0.5 µM of each primer, QIAGEN reaction buffer diluted 1:10 and 1U of Taq DNA polymerase (QIAGEN, Chatsworth, CA, USA) in a T-100 thermal cycler (Bio-Rad Japan, Tokyo, Japan). The protocol for the calmodulin was set as described by Prencipe et al. [53]. The PCR products were checked by gel electrophoresis at 1% agarose. PCR products were sequenced using the external service of GATC services (Costanza, Germany) in both directions. Contigs of forward and reverse sequences were carried out using the DNA Baser program (Heracle BioSoft SRL, Pitești, Romania).
Alignment of the sequences was performed in CLUSTALW (multiple sequence alignment) using MEGA 6 program [76] with default parameters. The calmodulin sequences of some the isolates listed in the Table S3, previously identified by Frisvad et al. [4] were included in the alignment, as well as the sequence of some chestnut isolates obtained from GenBank (isolates AFSP6, AFCAL9 and AFSP8). The sequences were then edited manually and trimmed obtaining a data set of 461 bp.
The best-fit model (K80 +G) was determined using the tool Model Selection in Topali v2.5 program [77] for maximum likelihood that was carried out using the same software with 100 bootstrap replicates to evaluate the stability. Cladogram tree topologies were adjusted using FigTree v1.4.3 and represented in the Figure S1.

5.2. In vitro Aflatoxin Production

The aflatoxigenic potential of the strains listed in Table 1 and Table 2 and the 7 reference strains (Table S1) was evaluated in an in vitro assay with HPLC coupled fluorescence detection. The protocol used to induce aflatoxin production was based on the in vitro assay described by Prencipe et al. [53]. The strains were grown on YES agar [71] for 1 week at 30 °C in the dark. Three agar plugs (1 cm diameter) were used to inoculate 50 mL of YES broth [74], subsequently incubated for one week at 30 °C in the dark.
The culture broth was filtered through a Whatman filter paper, grade 1 (VWR International, Milan, Italy), to check the aflatoxin production. The broth (10 mL) was extracted with ethyl acetate (Sigma Aldrich, Steinheim, Germany) performing the liquid–liquid extraction twice. The organic extracts were mixed and evaporated to dryness in a rotatory evaporator. The residue was dissolved in 0.5 mL of methanol-water (50:50 v/v) for HPLC analysis. All the samples were diluted to remain in the linearity range. All the strains were grown in triplicate under the same conditions, and AFs were extracted twice.
Analyses were performed using an Agilent 1260 Series (Agilent Technologies, Santa Clara, CA, USA) equipped with a fluorescent detector and an electrochemical cell (KOBRA® CELL, R-biopharm, Darmstadt, Germany) for the post-column derivatization. Chromatographic separation was performed on a Pursuit XRs Ultra C18 (100 mm × 2.0 mm, 2.8 µm, Varian) column in isocratic mode (80:20) using H2O:CH3OH:CH3CN (50:40:10 v/v/v) acidified with 0.035% HNO3 and 119 mg/L KBr in H2O at flow rate of 1 mL/min for 15 min. Quantification of AFs was performed by external calibration method using standard curves of the four aflatoxins (AFG1, AFG2, AFB1, AFB2, Sigma Aldrich, crystallized form) at the concentration range of 0.1 to 25 ppb. The quantity of AFs reported is the average of six replicates (three biological replicates and two technical replicates per biological one).

5.3. HPLC-FLD Method Validation

The developed analytical method was evaluated for recovery, limit of detection (LOD), limit of quantification (LOQ) and matrix effect (ME) for the four aflatoxins in accordance with Matuszewski et al. [78] and Kos et al. [79]. Analytical standards of AFB1, AFB2, AFG1, and AFG2 (Sigma-Aldrich) were used to prepare two sets of samples: (i) standard solution dissolved in mobile phase in a range of concentrations from 0.1 up to 25 ppb; (ii) spiked matrix: the extraction was performed from YES broth media, afterwards, the four aflatoxin standards were added in the range from 0.1 up to 25 ppb. Signal-to-noise method was used to determine limits of detection (LOD) and quantification (LOQ) for the mobile phase and the matrix. S/N ratio 3:1 was used for the determination of LOD, while 10:1 for LOQ.

5.4. Confirmation of Deletion in the Aflatoxin Gene Cluster. PCR LOD

Two manually designed primers (AfF and AfR) were used to check the presence or absence of the region in the aflatoxin biosynthetic gene cluster where the LAMP and qPCR primers were later designed. The primers were designed flanking the region in which the LAMP primers were designed amplifying a region of approximately 308 bp. Possible secondary structures of the primers were also tested as explained above.
A PCR reaction using AfF and AfR of each strain used in this study (Table 1 and Table 2), was performed in 20 µL with approximately 50 ng of genomic DNA, 2.5 mM deoxynucleotide triphosphates, 2.5 mM MgCl2, 0.5 µM of each primer, Qiagen reaction buffer diluted 1:10 and 1U of Taq DNA polymerase (Qiagen, Düsseldorf, Germany) in a T-100 thermal cycler (Bio-Rad Japan, Tokyo, Japan). The protocol included an initial denaturation step at 95 °C for 3 min and 34 cycles, including 95 °C for 30 s, 59.5 °C for 30 s and 72 °C for 1 min, with a final extension at 72 °C for 7 min. The PCR products were checked by gel electrophoresis at 2% agarose (VWR International, Milan, Italy). Limit of detection of the technique was calculated using 10-fold serial dilutions of the A. flavus FV3 (initial concentration 110 ng/µL).

5.5. LAMP and TaqMan® Assay Primers Design

The primers for LAMP and qPCR were designed on the biosynthetic aflatoxin gene cluster, specifically aligning to the aflT gene (transmembrane protein; at the end of the last exon and the beginning of the last intron), in order to detect and quantify the DNA of the aflatoxigenic Aspergillus section Flavi strains, with special interest in A. flavus and A. parasiticus. Due to the difficulty to find specific primers with any software, the six LAMP primers were manually designed: two external primers (F3 and B3), two loop primers (F-loop and B-loop) and two internal primers (BIP and FIP) following the instructions described on-line [80,81]; two primers and a TaqMan® probe were designed for quantitative PCR in the same region (Figures S2 and S3). A multilocus alignment of the complete biosynthetic aflatoxin gene cluster was performed using the sequences of A. flavus (GenBank accession number AY510453.1, AY510455.1, AY510452.1), A. oryzae (GenBank accession number AB196490.1), A. parasiticus (GenBank accession number AY371490.1) and A. nomius (GenBank accession number AY540454.1) as shown in the Figure S2. The alignment was carried out using CLUSTALW (multiple sequence alignment) using MEGA 6 program [74] with default parameters. Possible secondary structures, hairpins, and self-annealing were checked for six LAMP primers, the TaqMan® probe and the two qPCR primers using OligoCalc (http://biotools.nubic.northwestern.edu/OligoCalc.html). BLASTn analysis of all the primers and the TaqMan® probe (National Center for Biotechnology Information, https://blast.ncbi.nlm.nih.gov/) was also carried out to check the specificity. In order to verify the LAMP and qPCR products, the products of the PCR amplicon obtained with DNA of A. flavus FP8, BV10 and BS3 and with the AfF and AfR primers were sequenced using services of Macrogen Europe (Amsterdam, The Netherlands). The primers of the LAMP, qPCR and the probe were then aligned against the three sequences. In order to understand the localization of the primers, they were aligned against the sequence (exons and introns) of the aflT from A. flavus (GenBank: AF515601.1). All primers were synthesized and HPLC purified by Eurofins (UK). Probe for the qPCR assay was synthesized by Metabion International AG (Planegg, Germany).

5.6. LAMP and qPCR Reactions

The LAMP reaction was carried out with 1 µL of DNA for the strains listed in Table 1 and Table 2, Tables S1 and S2 with concentration adjusted to around 1–50 µg/µL or 1 µL of crudely extracted DNA of artificially inoculated and real sample hazelnuts, as explained below, or 1 µL of the DNA extracted using E.Z.N.A. Plant DNA kit (Omega Bio-tek, Inc., Norcross, GA, USA) from inoculated and real sample hazelnuts.
The LAMP reaction was performed in 25 µL including 0.3 µmol/L of the external primers F3 and B3, 1.2 µmol/L of each loop primer (F-loop and B-loop) and 2.3 µmol/L of each internal primer (FIP and BIP) and 1× of the Isothermal Mastermix ISO-004® which contains a ds-DNA binding fluorescent dye (FAM detection channel) (Optigene Ltd., Horsham, UK). The reaction was performed as described by Franco Ortega et al. [43] for 60 min at 65 °C, recording the melting temperature whilst cooling from 95 to 70 °C at 0.05 °C/s using a Genie II® instrument (Optigene Ltd.) and a StepOne thermal cycler (Applied Biosystems, Loughborough, UK). Positive controls, which included two strains able to produce aflatoxins (the AFB1 and AFB2 producer strain FS7 and previously identified as A. flavus and the A. parasiticus strain AFLX6 fungal DNA) and two negative controls, one non-aflatoxigenic strain (A. flavus AFSP1) and one with water, were included in each run. The average time to positive (Tp) of three technical replicates was used to calculate the Tp of each sample. Cytochrome oxidase I (COX) LAMP reaction was also performed using the protocol described by [52] as an internal control to confirm the presence of hazelnut DNA after the crude DNA extraction.
The qPCR assay was initially performed using the two designed qPCR primers (AflF and AflR) with SYBR™ Green to determine the melting temperature of the amplicons, and afterwards with the same two primers (AflF and AflR) and the TaqMan® Afl-probe designed for this study.
In the SYBR Green experiment, the reaction was carried out with 0.3 µmol/L of each primer, 1× of the Power SYBR™ Green Master Mix (Thermo Fisher) and 1 µL of the DNA of each isolate listed in Table 1 and Table 2, Tables S1 and S2. The reactions were performed using a StepOne instrument (Applied Biosystems) with the following conditions: initial denaturation at 95 °C for 4.5 min, 40 cycles of 15 s at 95 °C and 15 s at 60 °C. The melting temperature using SYBR™ Green was measured while heating from 60 to 95 °C at 0.03 °C/s. Each reaction was performed in triplicate and positive (A. flavus FS7 and A. parasiticus strain AFLX6 fungal DNA) and negative controls (non-aflatoxigenic AFSP1 and water) were included in each run.
The same DNA was used to perform the qPCR with a TaqMan® probe, using a final concentration of 0.3 µmol/L of each primer (AflF and AflR), 0.1 µmol/L of the TaqMan® probe, 1× of TaqMan® Universal PCR MasterMix (Applied Biosystems) and 1 µL of the DNA of each isolate listed in Table 1 and Table 2, Tables S1 and S2 or 1 µL of the hazelnut DNA (inoculated and real samples) extracted by using the E.Z.N.A. Plant DNA kit (Omega Bio-tek, Inc.). This assay was not performed using the DNA from the crude extracted hazelnuts. The reactions were performed using a StepOne thermal cycler with the setup explained above. Positive controls (A. flavus FS7 and A parasiticus AFLX6 fungal DNA) and negative controls (the non -aflatoxigenic strain AFSP1 and water) were included in each run.
The average of three technical replicates of the qPCR assay with TaqMan® probe was used to calculate the number of cells present in the DNA extraction from artificially and naturally infected hazelnuts. The Ct values generated by the TaqMan® Real Time qPCR assay were compared against the standard curve to calculate the amount (ng) of DNA in 1 µl of each positive result. The standard curve, to quantify the aflatoxigenic Aspergillus present in hazelnuts, was built using the DNA of the aflatoxigenic strains: FS7 (concentrations ranging from 16.9 ng/µL to 16.9 fg/µL) and A. parasiticus isolate AFLX6 (concentrations ranging from 1.81 ng/µL to 18.3 fg/µL). The standard curve generated with FS7 was used to quantify hazelnuts inoculated with A. flavus, as the standard curve generated with A. parasiticus AFLX6 was used to quantify hazelnuts inoculated with the latter strain. To quantify naturally infected samples, the standard curve generated with the A. flavus FS7 was used, as A. flavus is the most common aflatoxigenic species found in hazelnuts [14]. In order to calculate the number of cells in the inoculated hazelnuts, the elution volume of the commercial DNA (30 µL) of the inoculated hazelnuts was also taken into consideration, as described in other in vivo tests [67].
In all cases, the approximate number of cells detected was calculated by dividing the DNA quantity by the weight of the genome of Aspergillus flavus strain NRRL3357 (0.0000397 ng).

5.7. LAMP and qPCR Validation

Validation of the LAMP assay and the TaqMan® qPCR was performed according to the EPPO standard PM7/98 in terms of specificity, sensitivity, repeatability, and reproducibility. The specificity testing included all the Aspergillus spp. in Table 1, 2 and Table S1 and other common fungi found in hazelnuts and in chestnuts (Table S2). The sensitivity testing of the LAMP and qPCR assays was carried out using 10-fold serial dilution of DNA of two aflatoxigenic species (A. flavus strains FS7 and A. parasiticus strain AFLX6) and in the case of the LAMP, the amplification was also confirmed through an agarose gel. The sensitivity was also related to the in vitro aflatoxin production and later, confirmed on artificially inoculated hazelnuts with the isolates: AFSP1 (non-aflatoxigenic A. flavus), AFSP4 (aflatoxigenic A. flavus) and AFLX6 (aflatoxigenic A. parasiticus), at different concentrations. All the samples were tested in triplicate, and an average of the obtained results was reported (repeatability). Reproducibility was checked by repeating the experiment in different machines, in different days and by two members of the staff.
Analytical specificity and sensitivity of all the assays (PCR, LAMP and qPCR with SYBR Green and TaqMan® probe), using DNA of 97 Aspergillus spp. were evaluated using the calculator tool Diagnostic Test Calculator [82]. Diagnostic sensitivity (DSe) or true-positive rate and analytical specificity (DSp) or true-negative rate, were calculated using the following formulas: Dse = ΣTP/(ΣTP + FN); Dpe = ΣTN/(ΣTN + FP); where TP (True Positive) represents the aflatoxigenic strains, that showed a positive amplification with LAMP or TaqMan® qPCR, TN (True Negative) denotes the non-aflatoxigenic strains that gave a negative result using LAMP or TaqMan® qPCR, FN (False Negatives) indicates samples which gave a negative result using the LAMP or TaqMan® qPCR assay but able to produce at least one of the 4 AFs, and FP (False Positives) are the samples showing positive amplification of LAMP or TaqMan® qPCR but did not produce any AFs. The likelihood positive ratio (LR+) and the likelihood negative ratio (LR−) were calculated according to the ratio between DSe and the FP-rate and the FP-rate and Dse independently for the LAMP and for the qPCR. In addition, these parameters were also calculated for the artificially inoculated hazelnuts, considering as separate assays the LAMP results obtained with crude and the commercial DNA extraction. In this case, each individual technical replicate was used for the calculation.

5.8. Validation Using Artificially Inoculated Hazelnut

A batch of Italian hazelnuts kindly provided by Ferrero S.p.A. (Alba, Italy) was used to perform the sensitivity test of both LAMP and qPCR techniques in order to check the lowest detectable level of the fungi.
The hazelnuts were previously washed in water with 1% Tween 80 (Sigma Aldrich, Gillingham, UK) for 5 min, treated with sodium hypochlorite (2%), rinsed with distilled water and then air-dried for 30 min.
At the same time, spore-suspensions of three different Aspergillus spp. used for the specificity test were prepared. The non-aflatoxigenic A. flavus strain AFSP1, the aflatoxigenic A. flavus strain AFSP4 and the aflatoxigenic A. parasiticus strain AFLX6, were previously grown in YES agar at 30 °C for 10 days, as previously described. Spores-suspension with water and Tween 1% was obtained from the YES-agar plates. The spore concentrations were adjusted by means of a hematocytometer, obtaining 3.09 × 107 for AFSP1, 3.03 × 107 for AFSP4 and 1.79 × 107 for AFLX6. Hazelnuts were artificially inoculated under sterile conditions with 10-fold serial dilution of the three strains: the non-aflatoxigenic A. flavus AFSP1 was inoculated with concentrations 3.09 × 107, 3.09 × 106, 3.09 × 105 spores/mL, A. flavus AFSP4 was inoculated with concentrations 3.03 × 107, 3.03 × 106, 3.03 × 105 spores/mL and A. parasiticus AFLX6 with concentrations 1.79 × 107, 1.79 × 106, 1.79 × 105 spores/mL in independent assays. The inoculation was performed by immersion of 20 hazelnuts in 50 mL of spore suspension for 30 min in agitation (100 rpm) at room temperature. Hazelnuts were then air-dried overnight.
The DNA was extracted from six hazelnuts (for each inoculated concentration) with two different approaches. Three hazelnuts were used to perform a crude DNA extraction based on an alkaline extraction modified from the protocol described by Chomczynski and Rymaszewski [83] and subsequently used to perform a LAMP assay. The crude extraction buffer was composed by 1 mL of pH 13 PEG buffer (50 g/L PEG average Mm of 4600; 20 mmol/L KOH; pH 13.5) as described in Franco Ortega et al. [43] and were placed in a 5 mL tubes and homogenized by vigorous manual shaking for one minute. The solutions were then diluted 1:50 using 10 µL of crude extract in 490 µL of water. The other three hazelnuts of each concentration were ground with liquid nitrogen, extracted using E.Z.N.A. Plant DNA kit (Omega Bio-Tek) according to manufacturer’s instructions and eluted in 30 µL of water.
LAMP assay was performed using 1 µL of both DNA extraction methods in three different assays (technical replicates), whilst the TaqMan® qPCR assay was performed in triplicate using 1 µL of the DNA obtained with a commercial extraction kit. The results of both types of extraction in the LAMP testing were then compared.

5.9. Validation Using Real Samples

Fifty hazelnuts were used to check the reliability of both assays, LAMP and qPCR, using naturally contaminated samples. The hazelnut shells were removed and the hazelnuts were extracted with both types of DNA extraction (crude DNA extraction and using the E.Z.N.A. Plant DNA kit eluted in a final volume of 30 µL) as described before. LAMP assay was performed using 1 µL of both DNA extracts in three independent replicates, whilst the qPCR assay was performed in triplicate using only 1 µL of the extract obtained with the commercial DNA extraction kit. The approximate number of cells on the crude DNA extract was calculated considering the average Ct obtained with 1 µL of the hazelnuts extracted with a commercial kit and using the standard curve of the strain FS7. The elution volume (30 µL) was also considered to calculate the approximate number of cells per hazelnut.
Fifteen hazelnuts from the same batch were used to check the presence of Aspergillus spp. using a traditional approach. Hazelnuts were agitated in 100 mL sterile water with 1% Tween solution for 2 h. The suspension was serially diluted (1:10, 1:102, 1:103, 1:104 and 1:105) and plated onto AFPA medium [84] in triplicate. Sampling of the Aspergillus spp. and the incidence was calculated after two days of incubation at 30 °C in the dark. The identification was performed visually according to the presence of orange-yellow pigmentation on the reverse side of the colonies [53].

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2072-6651/12/12/757/s1, Figure S1: Phylogenetic tree based on the calmodulin gene sequences constructed with Maximum Likelihood method using K80 model plus gamma distribution model. Sequences from different species belonging to the Aspergillus section Flavi, identified by Frisvad et al. [4] were included as references, Figure S2: LAMP, qPCR and end-point PCR primers alignment with the aflT gene within the aflatoxin biosynthetic cluster, Figure S3: LAMP, qPCR and end-point PCR primers alignment with the aflT gene of A. flavus (GenBank AF515601.1) showing the exons and introns; Figure S4: Agarose gel of the 308 bp amplicon obtained with primers AfF and AfR of the 97 Aspergillus spp. used in this study, Figure S5: Agarose gel showing the limit of detection of the 308 bp PCR with the primers AfF and AfR. 1: GelPilot Wide Range Ladder (QIAGEN); 2: GelPilot 100 bp Ladder (QIAGEN), 3: A. flavus FV3 (110 ng/µL); 4: A. flavus FV3 (11 ng/µL); 5: A. flavus FV3 (1.1 ng/µL); 6: A. flavus FV3 (110 pg/µL); 7: A. flavus FV3 (11 pg/µL); 8: A. flavus FV3 (1.1 pg/µL); 9: A. flavus FV3 (110 fg/µL); 10: A. flavus FV3 (11 fg/µL); 11: Negative Control with water; 12: GelPilot 100 bp Ladder (QIAGEN); 13: GelPilot Wide Range Ladder (QIAGEN), Figure S6: LAMP amplifications exported from the Genie II® instrument (A), and melting peaks (B) of 10-fold dilutions of the DNA of A. flavus FS7 strain. The amplification was confirmed by running an agarose gel of one of the three replicates of each dilution (C), Figure S7: qPCR standard curves with SYBR Green (A), and TaqMan ® probe standard curves using the DNA of A. flavus (FS7) (B) and A. parasiticus (AFLX6) (C), Table S1: Reference strains used in this study to check the specificity of the LAMP assay (Tp, time to positive) and the qPCR assays (Ct). The PCR results reported the presence or absence of a 308 bp-amplicon. Aflatoxin production is reported in ppm. Neg = negative, Pos = positive; n.d.: not detected; Melting: melting temperature, Table S2: Other species used in the specificity panel of the LAMP assay (Tp, time to positive) and the qPCR assay (Ct). The PCR results reported the presence or absence of a 308 bp-amplicon. Standard deviation was included for each Tp and melting temperature. Tp: time to positive. Anneal: melting temperature, Table S3: Reference strains used for the identification of the hazelnut isolates. The strains were previously identified by Frisvad et al. [4]. GenBank accession number report the sequence for the calmodulin used to infer a phylogenetic tree that identifies the hazelnut isolates.

Author Contributions

Conceptualization, D.S. and S.F.O.; methodology, S.F.O. and I.S.; validation, S.F.O., I.S., and S.P.; data curation, S.F.O., I.S. and S.P.; writing—Original draft preparation, S.F.O.; writing—Review and editing, S.F.O., I.S., S.P.; M.L.G. and D.S.; supervision, D.S. and M.L.G.; funding acquisition, D.S. and M.L.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Piedmont Region with the project “Food Digital Monitoring—FDM” (Framework Program MIUR—Piedmont Region, Action 3 “Smart Factory”, Call “Smart Factory Technological Platform”) and by CRT Foundation with the project “INNOCHEST—Innovative technologies to guarantee the quality and safety of Piedmontese chestnut”.

Acknowledgments

The authors gratefully acknowledge Ferrero S.p.A for providing the hazelnuts used for this study.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Iram, W.; Anjum, T.; Iqbal, M.; Ghaffar, A.; Abbas, M.; Khan, A.M. Structural analysis and biological toxicity of aflatoxins B1 and B2 degradation products following detoxification by Ocimum basilicum and Cassia fistula aqueous extracts. Front. Microbiol. 2016, 7, 1105. [Google Scholar] [CrossRef] [PubMed]
  2. Cotty, P.J.; Bayman, P.; Egel, D.S.; Elias, K.S. Agriculture, Aflatoxins and Aspergillus. In The Genus Aspergillus; Peberdy, J.F., Ed.; Springer: Boston, MA, USA, 1994; p. 127. [Google Scholar]
  3. Cheraghali, A.A.; Yazdanpanah, H. Interventions to control aflatoxin contamination in pistachio nuts: Iran experience. J. Food Saf. 2010, 30, 382–397. [Google Scholar] [CrossRef]
  4. Frisvard, J.C.; Hubka, V.; Ezekiel, C.N.; Hong, S.B.; Nováková, A.; Chen, A.J.; Arzanlou, M.; Larsen, T.O.; Sklenář, F.; Mahakarnchanakul, W.; et al. Taxonomy of Aspergillus section Flavi and their production of aflatoxins, ochratoxins and other mycotoxins. Stud. Mycol. 2019, 93, 1–63. [Google Scholar] [CrossRef]
  5. Kurtzman, C.P.; Smiley, M.J.; Robnett, C.J.; Wicklow, D.T. DNA relatedness among wild and domesticated species in Aspergillus flavus group. Mycologia 1986, 78, 955–959. [Google Scholar] [CrossRef]
  6. Passone, M.A.; Rosso, L.C.; Ciancio, A.; Etcheverry, M.G. Detection and quantification of Aspergillus section Flavi spp. in stored peanuts by real-time PCR of nor-1 gene, and effects of storage conditions on aflatoxin production. Int. J. Food Microbiol. 2010, 138, 276–281. [Google Scholar] [CrossRef] [PubMed]
  7. Dorner, J.W.; Cole, R.J.; Lomax, L.G.; Gosser, H.S.; Diener, U.L. Cyclopiazonic acid production by Aspergillus flavus and its effects on broiler chickens. Appl. Environ. Microbiol. 1983, 46, 698–703. [Google Scholar] [CrossRef] [Green Version]
  8. Payne, G.A.; Brown, M.P. Genetics and physiology of aflatoxin biosynthesis. Annu. Rev. Phytopathol. 1998, 36, 329–362. [Google Scholar] [CrossRef]
  9. Chang, P.K.; Horn, B.W.; Dorner, J.W. Sequence breakpoints in the aflatoxin biosynthesis gene cluster and flanking regions in non-aflatoxigenic Aspergillus flavus isolates. Fungal Genet. Biol. 2005, 42, 914–923. [Google Scholar] [CrossRef]
  10. Durak, I.; Koksal, I.; Kacmaz, M.; Buyukkocak, S.; Cimen, B.M.Y.; Ozturk, H.S. Hazelnut supplementation enhances plasma antioxidant potential and lowers plasma cholesterol levels. Clin. Chim. Acta. 1999, 284, 113–115. [Google Scholar] [CrossRef]
  11. Alasalvar, C.; Amaral, J.S.; Shahidi, F. Functional lipid characteristics of Turkish Tombul hazelnut (Corylus avellana L.). J. Agric. Food Chem. 2006, 54, 10177–10183. [Google Scholar] [CrossRef]
  12. FAOSTAT. Available online: www.fao.org/faostat/en/#data/QC (accessed on 20 September 2020).
  13. Valente, S.; Meloni, G.R.; Prencipe, S.; Spigolon, N.; Somenzi, M.; Fontana, M.; Gullino, M.L.; Spadaro, D. Effect of drying temperatures and Exposure Times on Aspergillus flavus growth and aflatoxin production on artificially inoculated hazelnuts. J. Food. Prot. 2020, 83, 1241–1247. [Google Scholar] [CrossRef] [PubMed]
  14. Ozay, G.; Seyhan, F.; Pembeci, C.; Saklar, S.; Yilmaz, A. Factors influencing fungal and aflatoxin levels in Turkish hazelnuts (Corylus avellana L.) during growth, harvest, drying and storage: A 3-year study. Food Addit. Contam. Part. A Chem. Anal. Control. Expo. Risk Assess. 2008, 25, 209–218. [Google Scholar] [CrossRef] [PubMed]
  15. Nakai, V.K.; de Oliveira Rocha, L.; Gonçalez, E.; Fonseca, H.; Ortega, E.M.M.; Corrêa, B. Distribution of fungi and aflatoxins in a stored peanut variety. Food Chem. 2008, 106, 285–290. [Google Scholar] [CrossRef]
  16. Joffe, A.Z. Relationships between Aspergillus flavus, A. niger and some other fungi in the mycoflora of groundnut kernels. Plant. Soil. 1969, 31, 57–64. [Google Scholar] [CrossRef]
  17. Bayman, P.; Baker, J.L.; Mahoney, N.E. Aspergillus on tree nuts: Incidence and associations. Mycopathologia 2003, 155, 161–169. [Google Scholar] [CrossRef]
  18. Singh, P.K.; Shukla, A.N. Survey of mycoflora counts, aflatoxin production and induced biochemical changes in walnut kernels. J. Stored Prod. Res. 2008, 44, 169–172. [Google Scholar] [CrossRef]
  19. Battilani, P. Recent advances in modeling the risk of mycotoxin contamination in crops. Curr. Opin. Food Sci. 2016, 11, 10–15. [Google Scholar] [CrossRef]
  20. Commission of the European Communities. Commission Regulation (EC) No 165/2010 of 26 February 2010 amending regulation (EC) No 1881/2006 setting maximum levels for certain contaminants in foodstuffs as regards aflatoxins. Off. J. Eur. Union 2010, 50, 8–12. [Google Scholar]
  21. FDA. Available online: https://www.fda.gov/regulatory-information/search-fda-guidance-documents/guidance-industry-action-levels-poisonous-or-deleterious-substances-human-food-and-animal-feed (accessed on 6 November 2020).
  22. Fukushihoken. Available online: https://www.fukushihoken.metro.tokyo.lg.jp/shokuhin/eng/kabi/kabidoqa.html#:~:text=In%20Japan%2C%20a%20limit%20of,all%20foods%20since%20October%202011 (accessed on 6 November 2020).
  23. Shapira, R.; Paster, N.; Eyal, O.; Menasherov, M.; Mett, A.; Salomon, R. Detection of aflatoxigenic molds in grains by PCR. Appl. Environ. Microbiol. 1996, 62, 3270–3273. [Google Scholar] [CrossRef] [Green Version]
  24. Chen, R.S.; Tsay, J.G.; Huang, Y.F.; Chiou, R.Y.Y. Polymerase chain reaction-mediated characterization of molds belonging to the Aspergillus flavus group and detection of Aspergillus parasiticus in peanut kernels by a multiplex polymerase chain reaction. J. Food Prot. 2002, 65, 840–844. [Google Scholar] [CrossRef]
  25. Somashekar, D.; Rati, E.; Chandrashekar, A. PCR-restriction fragment length analysis of aflR gene for differentiation and detection of Aspergillus flavus and Aspergillus parasiticus in maize. Int. J. Food Microbiol. 2004, 93, 101–107. [Google Scholar] [CrossRef] [PubMed]
  26. Scherm, B.; Palomba, M.; Serra, D.; Marcello, A.; Migheli, Q. Detection of transcripts of the aflatoxin genes aflD, aflO, and aflP by reverse transcription-polymerase chain reaction allows differentiation of aflatoxin-producing and non-producing isolates of Aspergillus flavus and Aspergillus parasiticus. Int. J. Food Microbiol. 2005, 98, 201–210. [Google Scholar] [CrossRef] [PubMed]
  27. Criseo, G.; Racco, C.; Romeo, O. High genetic variability in non-aflatoxigenic A. flavus strains by using Quadruplex PCR-based assay. Int. J. Food Microbiol. 2008, 125, 341–343. [Google Scholar] [CrossRef] [PubMed]
  28. Degola, F.; Berni, E.; Spotti, E.; Ferrero, I.; Restivo, F.M. Facing the problem of “false positives”: Re-assessment and improvement of a multiplex RT-PCR procedure for the diagnosis of A. flavus mycotoxin producers. Int. J. Food Microbiol. 2009, 129, 300–305. [Google Scholar] [CrossRef] [PubMed]
  29. Rodrigues, P.; Venâncio, A.; Kozakiewicz, Z.; Lima, N. A polyphasic approach to the identification of aflatoxigenic and non-aflatoxigenic strains of Aspergillus Section Flavi isolated from Portuguese almonds. Int. J. Food Microbiol. 2009, 129, 187–193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Ahmad, M.M.; Ahmad, M.; Ali, A.; Hamid, R.; Javed, S.; Abdin, M.Z. Detection of Aspergillus flavus and Aspergillus parasiticus from aflatoxin-contaminated peanuts and their differentiation using PCR-RFLP. Ann. Microbiol. 2014, 64, 1597–1605. [Google Scholar] [CrossRef]
  31. Criseo, G.; Bagnara, A.; Bisignano, G. Differentiation of aflatoxin-producing and non-producing strains of Aspergillus flavus group. Lett. Appl. Microbiol. 2001, 33, 291–295. [Google Scholar] [CrossRef]
  32. Goebes, M.D.; Hildemann, L.M.; Kujundzic, E.; Hernandez, M. Real-time PCR for detection of the Aspergillus genus. J. Environ. Monit. 2007, 9, 599–609. [Google Scholar] [CrossRef]
  33. Bintvihok, A.; Treebonmuang, S.; Srisakwattana, K.; Nuanchun, W.; Patthanachai, K.; Usawang, S.A. Rapid and sensitive detection of aflatoxin-producing fungus using an optimized polymerase chain reaction (PCR). Toxicol. Res. 2016, 32, 81–87. [Google Scholar] [CrossRef]
  34. Mitema, A.; Okoth, S.; Rafudeen, S.M. The Development of a qPCR Assay to measure Aspergillus flavus biomass in maize and the use of a biocontrol strategy to limit aflatoxin production. Toxins 2019, 11, 179. [Google Scholar] [CrossRef] [Green Version]
  35. Mayer, Z.; Bagnara, A.; Färber, P.; Geisen, R. Quantification of the copy number of nor-1, a gene of the aflatoxin biosynthetic pathway by real-time PCR, and its correlation to the CFU of Aspergillus flavus in foods. Int. J. Food Microbiol. 2003, 82, 143–151. [Google Scholar] [CrossRef]
  36. Shweta, S.; Madhavan, S.; Paranidharan, V.; Velazhahan, R. Detection of Aspergillus flavus in maize kernels by conventional and real-time PCR assays. Int. Food Res. J. 2013, 20, 3329–3335. [Google Scholar]
  37. Mahmoud, M.A. Detection of Aspergillus flavus in stored peanuts using Real-Time PCR and the expression of aflatoxin genes in toxigenic and atoxigenic A. flavus isolates. Foodborne Pathog. Dis. 2015, 12, 289–296. [Google Scholar] [CrossRef] [PubMed]
  38. Luo, Y.; Gao, W.; Doster, M.; Michailides, T.J. Quantification of conidial density of Aspergillus flavus and A. parasiticus in soil from almond orchards using real-time PCR. J. Appl. Microbiol. 2009, 106, 1649–1660. [Google Scholar] [CrossRef]
  39. Notomi, T.; Okayama, H.; Masubuchi, H.; Yonekawa, T.; Watanabe, K.; Amino, N.; Hase, T. Loop-mediated isothermal amplification of DNA. Nucleic Acids Res. 2000, 28, E63. [Google Scholar] [CrossRef] [Green Version]
  40. Nagamine, K.; Hase, T.; Notomi, T. Accelerated reaction by loop-mediated isothermal amplification using loop primers. Mol. Cell. Probes 2002, 16, 223–229. [Google Scholar] [CrossRef]
  41. Tomlinson, J.A.; Ostoja-Starzewska, S.; Adams, I.P.; Miano, D.W.; Abidrabo, P.; Kinyua, Z.; Alicai, T.; Dickinson, M.J.; Peters, D.; Boonham, N.; et al. Loop-mediated isothermal amplification for rapid detection of the causal agents of cassava brown streak disease. J. Virol. Methods. 2013, 191, 148–154. [Google Scholar] [CrossRef]
  42. Hodgetts, J.; Hall, J.; Karamura, G.; Grant, M.; Studholme, D.J.; Boonham, N.; Karamura, E.; Smith, J.J. Rapid, specific, simple, in-field detection of Xanthomonas campestris pathovar musacearum by loop-mediated isothermal amplification. J. Appl. Microbiol. 2015, 119, 1651–1658. [Google Scholar] [CrossRef]
  43. Franco Ortega, S.; Tomlinson, J.; Hodgetts, J.; Spadaro, D.; Gullino, M.L.; Boonham, N. Development of loop-mediated isothermal amplification assays for the detection of seedborne fungal pathogens, Fusarium fujikuroi and Magnaporthe oryzae, in rice seeds. Plant. Dis. 2018, 102, 1549–1558. [Google Scholar] [CrossRef] [Green Version]
  44. Luo, J.; Vogel, R.F.; Niessen, L. Development and application of a loop-mediated isothermal amplification assay for rapid identification of aflatoxigenic molds and their detection in food samples. Int. J. Food Microbiol. 2012, 159, 214–224. [Google Scholar] [CrossRef]
  45. Luo, J.; Taniwaki, M.H.; Iamanaka, B.T.; Vogel, R.F.; Niessen, L. Application of loop-mediated isothermal amplification assays for direct identification of pure cultures of Aspergillus flavus, A. nomius, and A. caelatus and for their rapid detection in shelled Brazil nuts. Int. J. Food Microbiol. 2014, 172, 5–12. [Google Scholar] [CrossRef] [PubMed]
  46. Liu, P.; Li, B.; Yin, R.; Weng, Q.; Chen, Q. Development and evaluation of ITS and aflP based LAMP assays for rapid detection of Aspergillus flavus in food samples. Can. J. Microbiol. 2014, 60, 579–584. [Google Scholar] [CrossRef] [PubMed]
  47. Al-Sheikh, H.M. LAMP-PCR detection of ochratoxigenic Aspergillus species collected from peanut kernel. Genet. Mol. Res. 2015, 14, 634–644. [Google Scholar] [CrossRef] [PubMed]
  48. Douksouna, Y.; Kwallah, A.O.; Nyerere, A.; Runo, S.; Ambang, Z. Application and Evaluation of the Loop-Mediated Isothermal Amplification assay for the Detection of Aflatoxigenic Fungi Contaminants of Rice Grains in Kenya. Plant Pathol Microbiol. 2020, 11, 491. [Google Scholar] [CrossRef]
  49. Ferrara, M.; Logrieco, A.F.; Moretti, A.; Susca, A. A loop-mediated isothermal amplification (LAMP) assay for rapid detection of fumonisin producing Aspergillus species. Food Microbiol. 2020, 90, 103469. [Google Scholar] [CrossRef]
  50. Niessen, L.; Bechtner, J.; Fodil, S.; Taniwaki, M.H.; Vogel, R.F. LAMP-based group specific detection of aflatoxin producers within Aspergillus section Flavi in food raw materials, spices, and dried fruit using neutral red for visible-light signal detection. Int. J. Food Microbiol. 2018, 266, 241–250. [Google Scholar] [CrossRef]
  51. Nixon, G.J.; Svenstrup, H.F.; Donald, C.E.; Carder, C.; Stephenson, J.M.; Morris-Jones, S.; Huggett, J.F.; Foy, C.A. A novel approach for evaluating the performance of real time quantitative loop-mediated isothermal amplification based methods. Biomol. Detect. Quantif. 2014, 2, 4–10. [Google Scholar] [CrossRef] [Green Version]
  52. Tomlinson, J.; Dickinson, M.J.; Boonham, N. Rapid Detection of Phytophthora ramorum and P. kernoviae by Two-Minute DNA extraction followed by isothermal amplification and amplicon detection by generic lateral flow device. Phytopathology 2010, 100, 143–1493. [Google Scholar] [CrossRef] [Green Version]
  53. Prencipe, S.; Siciliano, I.; Contessa, C.; Botta, R.; Garibaldi, A.; Gullino, M.L.; Spadaro, D. Characterization of Aspergillus section Flavi isolated from fresh chestnuts and along the chestnut flour process. Food Microbiol. 2018, 69, 159–169. [Google Scholar] [CrossRef]
  54. Magan, N.; Olsen, M. Mycotoxins in Food. Detection and Control; Woodhead Publishing Ltd.: Cambridge, UK; CRC Press LLC.: Boca Raton, FL, USA, 2004. [Google Scholar] [CrossRef]
  55. Riba, A.; Matmoura, A.; Mokrane, S.; Mathieu, F.; Sabaou, N. Investigations on aflatoxigenic fungi and aflatoxins contamination in some nuts sampled in Algeria. Afr. J. Microbiol. Res. 2013, 7, 4974–4980. [Google Scholar] [CrossRef] [Green Version]
  56. Pitt, J.I.; Hocking, A.D.; Glenn, D.R. An improved medium for the detection of Aspergillus flavus and A. parasiticus. J. Appl. Bacteriol. 1983, 54, 109–114. [Google Scholar] [CrossRef] [PubMed]
  57. Klich, M.A.; Mullaney, E.J. Use of a bleomycin-containing medium to distinguish Aspergillus parasiticus from A. sojae. Mycologia 1989, 81, 159–160. [Google Scholar] [CrossRef]
  58. Torres, J.; Guarro, J.; Suarez, G.; Suñe, N.; Calvo, M.A.; Ramírez, C. Morphological changes in strains of Aspergillus flavus link ex fries and Aspergillus parasiticus speare related with aflatoxin production. Mycopathologia. 1980, 72, 171–180. [Google Scholar] [CrossRef]
  59. Samson, R.A.; Visagie, C.M.; Houbraken, J.; Hong, S.B.; Hubka, V.; Klaassen, C.H.; Perrone, G.; Seifert, K.A.; Susca, A.; Tanney, J.B.; et al. Phylogeny, identification and nomenclature of the genus Aspergillus. Stud. Mycol. 2014, 78, 141–173. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Makhlouf, J.; Carvajal-Campos, A.; Querin, A.; Tadrist, S.; Puel, O.; Lorber, S.; Oswald, I.P.; Hamze, M.; Bailly, J.D.; Bailly, S. Morphologic, molecular and metabolic characterization of Aspergillus section Flavi in spices marketed in Lebanon. Sci. Rep. 2019, 9, 5263. [Google Scholar] [CrossRef] [PubMed]
  61. Chang, P.K.; Ehrlich, K.C.; Hua, S.S.H. Cladal relatedness among Aspergillus oryzae isolates and Aspergillus flavus S and L morphotypes. Int. J. Food. Microbiol. 2006, 108, 172–177. [Google Scholar] [CrossRef]
  62. Rodríguez, A.; Rodríguez, M.; Luque, M.I.; Martín, A.; Córdoba, J.J. Real-time PCR assays for detection and quantification of aflatoxin-producing molds in foods. Food Microbiol. 2012, 31, 89–99. [Google Scholar] [CrossRef]
  63. Yu, J.; Chang, P.K.; Ehrlich, K.C.; Cary, J.W.; Bhatnagar, D.; Cleveland, T.E.; Payne, G.A.; Linz, J.E.; Woloshuk, C.P.; Bennett, J.W. Clustered Pathway Genes in Aflatoxin Biosynthesis. Appl. Environ. Microbiol. 2004, 70, 1253–1262. [Google Scholar] [CrossRef] [Green Version]
  64. Krska, R.; Baumgartner, S.; Joseph, R. The state-of-the-art in the analysis of type-A and -B trichothecene mycotoxins in cereals. Fresenius J. Anal. Chem. 2001, 371, 285–299. [Google Scholar] [CrossRef]
  65. Probst, C.; Cotty, P.J. Relationships between in vivo and in vitro aflatoxin production: Reliable prediction of fungal ability to contaminate maize with aflatoxins. Fungal Biol. 2012, 116, 503–510. [Google Scholar] [CrossRef]
  66. Franco Ortega, S.; Bustos Lopez, M.; Nari, L.; Boonham, N.; Gullino, M.L.; Spadaro, D. Rapid Detection of Monilinia fructicola and Monilinia laxa on Peach and Nectarine using Loop-Mediated Isothermal Amplification. Plant. Dis. 2019, 103, 9. [Google Scholar] [CrossRef]
  67. Franco Ortega, S.; Prencipe, S.; Gullino, M.L.; Spadaro, D. New molecular tool for a quick and easy detection of apple scab in the field. Agronomy 2020, 10, 581. [Google Scholar] [CrossRef] [Green Version]
  68. Wilson, I.G. Inhibition and facilitation of nucleic acid amplification. Appl. Environ. Microbiol. 1997, 63, 3741–3751. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Franco Ortega, S.; Tomlinson, J.; Gilardi, G.; Spadaro, D.; Gullino, M.L.; Garibaldi, A.; Boonham, N. Rapid detection of Fusarium oxysporum f. sp. lactucae on soil, lettuce seeds and plants using loop-mediated isothermal amplification. Plant Pathol. 2018, 67, 1462–1473. [Google Scholar] [CrossRef] [Green Version]
  70. OEPP/EPPO Bulletin. PM 7/76 Use of EPPO Diagnostic Standards. EPPO Bull. 2018, 48, 373–377. [Google Scholar] [CrossRef] [Green Version]
  71. ISO-ISO 16050:2003. Available online: www.iso.org/standard/29628.html (accessed on 20 September 2020).
  72. Lund, F.; Frisvad, J.C. Penicillium verrucosum in wheat and barley indicates presence of ochratoxin A. J. Appl. Microbiol. 2003, 95, 1117–1123. [Google Scholar] [CrossRef] [Green Version]
  73. Frisvad, J.C. Physiological criteria and mycotoxin production as aids in identification of common asymmetric Penicillia. Appl. Environ. Microb. 1981, 41, 568–579. [Google Scholar] [CrossRef] [Green Version]
  74. Visagie, C.M.; Houbraken, J.; Frisvad, J.C.; Hong, S.B.; Klaassen, C.H.W.; Perrone, G.; Seifert, K.A.; Varga, J.; Yaguchi, T.; Samson, R.A. Identification and nomenclature of the genus Penicillium. Stud. Mycol. 2014, 78, 343–371. [Google Scholar] [CrossRef] [Green Version]
  75. Hong, S.B.; Go, S.J.; Frisvard, J.C.; Samson, R.A. Polyphasic taxonomy of Aspergillus fumigatus and related species. Mycologia 2005, 97, 1316–1329. [Google Scholar] [CrossRef]
  76. Tamura, K.; Stecher, G.; Peterson, D.; Filipski, A.; Kumar, S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol. Biol. Evol. 2013, 30, 2725–2729. [Google Scholar] [CrossRef] [Green Version]
  77. Milne, I.; Wright, F.; Rowe, G.; Marshall, D.F.; Husmeier, D.; McGuire, G. TOPALi: Software for automatic identification of recombinant sequences within DNA multiple alignments. Bioinformatics 2004, 20, 1806–1807. [Google Scholar] [CrossRef] [PubMed]
  78. Matuszewski, B.K.; Constanzer, M.L.; Chavez-Eng, C.M. Strategies for the assessment of matrix effect in quantitative bioanalytical methods based on HPLC-MS/MS. Anal. Chem. 2003, 75, 3019–3030. [Google Scholar] [CrossRef] [PubMed]
  79. Kos, J.; Janić Hajnal, E.; Jajić, I.; Krstović, S.; Mastilović, J.; Šarić, B.; Jovanov, P. Comparison of ELISA, HPLC-FLD and HPLC-MS/MS methods for determination of aflatoxin M1 in natural contaminated milk samples. Acta Chim. Slov. 2016, 63, 747–756. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Eiken Site. Available online: http://loopamp.eiken.co.jp/e/lamp/primer.html (accessed on 20 November 2020).
  81. Guide to LAMP Primer Designing of Primer Explorer v4. Available online: https://primerexplorer.jp/e/v4_manual/pdf/PrimerExplorerV4_Manual_1.pdf (accessed on 20 November 2020).
  82. Diagnostic Test Calculator. Available online: http://araw.mede.uic.edu/cgi-bin/testcalc.pl (accessed on 20 November 2020).
  83. Chomczynski, P.; Rymaszewski, M. Alkaline polyethylene glycol-based method for direct PCR from bacteria, eukaryotic tissue samples, and whole blood. Biotechniques 2006, 40, 454–458. [Google Scholar] [CrossRef] [PubMed]
  84. Corry, J.E.L.; Curtis, G.D.W.; Baird, R.M. Handbook of Culture Media for Food Microbiology; Elsevier Science: Amsterdam/Dordrecht, The Netherlands, 2003; Volume 37, p. 662. [Google Scholar]
Table 1. Strains from hazelnuts used in this study to check the specificity of the loop-mediated isothermal amplification (LAMP) assay (Tp, time to positive) and the TaqMan® qPCR assay (Ct). The PCR results reported the presence or absence of a 308 bp amplicon. Aflatoxin production is reported in ppm. Standard deviation was included for each Tp, Ct, melting temperature and for the four aflatoxins. ±: PCR result was not clear. Neg = negative, Pos = positive; n.d.: not detected; Melting = melting temperature.
Table 1. Strains from hazelnuts used in this study to check the specificity of the loop-mediated isothermal amplification (LAMP) assay (Tp, time to positive) and the TaqMan® qPCR assay (Ct). The PCR results reported the presence or absence of a 308 bp amplicon. Aflatoxin production is reported in ppm. Standard deviation was included for each Tp, Ct, melting temperature and for the four aflatoxins. ±: PCR result was not clear. Neg = negative, Pos = positive; n.d.: not detected; Melting = melting temperature.
NameSpeciesPCR LAMP AssayqPCR TaqMan® Probe (Ct)qPCR with SYBR GreenAflatoxin Production (ppm)
Tp (min:s)Melting (°C)CtMelting (°C)AFB1AFB2AFG1AFG2
Strains from Leaves
FS2A. flavusPos23:21 ± 03:1089.56 ± 0.0925.68 ± 0.1212.98 ± 0.079.79 ± 0.0038.8 ± 0.72.3 ± 0.2n.d.n.d.
FS3A. flavusNegNeg84.13 ± 0.265Neg15.66 ± 0.079.94 ± 0.21n.d.n.d.n.d.n.d.
FS4A. flavusNegNeg84.03 ± 0.093Neg34.04 ± 5.164.74 ± 10.53n.d.n.d.n.d.n.d.
FS5A. flavusNegNeg84.18 ± 0.086Neg22.60 ± 1.079.64 ± 0.00n.d.n.d.n.d.n.d.
FS6A. flavusNegNeg84.23 ± 0.088Neg23.93 ± 0.679.34 ± 0.11n.d.n.d.n.d.n.d.
FS7A. flavusPos21:33 ± 01:5289.61 ± 0.1720.38 ± 0.4623.10 ± 0.879.19 ± 0.218.5 ± 1.51.6 ± 0.5n.d.n.d.
FS8A. parasiticusPos20:47 ± 02:0789.51 ± 0.0120.11 ± 0.36Neg65.78 ± 0.4247,491 ± 50717,262 ± 11375,443 ± 169828,789 ± 669
FS9A. flavusPos19:19 ± 01:1089.21 ± 0.1523.49 ± 0.9337.71 ± 0.965.03 ± 0.7443.2 ± 8.42.2 ± 0.6n.d.n.d.
FS10Aspergillus sp. NegNeg84.13 ± 0.301Neg18.18 ± 0.779.79 ± 0.11n.d.n.d.n.d.n.d.
FV1A. flavusPos22:37 ± 03:2789.51 ± 0.1526.57 ± 0.25NegNeg238 ± 89.4n.d.n.d.n.d.
FV2A. flavusPos20:14 ± 02:1389.60 ± 0.0826.79 ± 0.4037.10 ± 0.165.19 ± 0.2127.6 ± 3.4n.d.n.d.n.d.
FV3A. flavus.Pos19:43 ± 01:2089.61 ± 0.0925.13 ± 0.7329.30 ± 0.279.19 ± 0.00314 ± 31.31.8n.d.n.d.
FV4A. flavusNegNeg84.33 ± 0.090Neg34.44 ± 6.766.82 ± 8.53n.d.n.d.n.d.n.d.
FV5A. flavusPos22:57 ± 02:2189.50 ± 0.0124.09 ± 0.4017.20 ± 0.179.78 ± 0.1146,996 ± 41779871 ± 288n.d.n.d.
FV7A. flavusNegNeg84.22 ± 0.092Neg33.54 ± 3.479.34 ± 0.63n.d.n.d.n.d.n.d.
FV8A. flavusPos20:11 ± 04:4189.66 ± 0.0124.40 ± 0.2730.18 ± 1.779.49 ± 0.00133 ± 8.8n.d.n.d.n.d.
FV9A. flavusNegNeg84.29 ± 0.009Neg29.58 ± 0.479.34 ± 0.11n.d.n.d.n.d.n.d.
FV10A. flavusNegNeg84.29 ± 0.010Neg12.68 ± 0.479.78 ± 0.00n.d.n.d.n.d.n.d.
FP1A. flavusNegNeg84.08 ± 0.081Neg37.52 ± 0.866.66 ± 9.70n.d.n.d.n.d.n.d.
FP2A. flavusPos21:40 ± 02:2889.61 ± 0.0922.21 ± 0.3130.56 ± 0.279.48 ± 0.11106 ± 7.313.7 ± 1.5n.d.n.d.
FP3A. flavusPos21:31 ± 01:2289.55 ± 0.0922.40 ± 0.2813.70 ± 0.079.77 ± 0.11495 ± 10112.9 ± 4.4n.d.n.d.
FP4Aspergillus sp. Pos22:35 ± 02:2589.66 ± 0.0121.13 ± 0.9114.08 ± 0.279.77 ± 0.11369 ± 37.111.8 ± 0.6n.d.n.d.
FP5Aspergillus sp. NegNeg84.18 ± 0.176Neg37.96 ± 1.766.34 ± 9.60n.d.n.d.n.d.n.d.
FP6Aspergillus sp. NegNeg84.33 ± 0.085Neg15.75 ± 0.079.47 ± 0.00n.d.n.d.n.d.n.d.
FP8A. flavusPosNeg84.24 ± 0.09836.79 ± 0.4034.3 ± 1.179.62 ± 0.00n.d.n.d.n.d.n.d.
FP9A. flavusPos20:52 ± 01:3289.51 ± 0.1622.39 ± 0.1413.59 ± 0.079.93 ± 0.0012.9 ± 1.12.9 ± 1n.d.n.d.
Strains from Bracts
BV2Aspergillus sp. Pos21:38 ± 01:4989.43 ± 0.1124.21 ± 0.6926.00 ± 0.279.33 ± 0.32493 ± 10812.5 ± 4.1n.d.n.d.
BV4A. flavusPos22:10 ± 02:0489.43 ± 0.3125.47 ± 0.5813.85 ± 0.479.62 ± 0.11255 ± 51.52.5 ± 1.2n.d.n.d.
BV5A. flavusPos22:58 ± 02:3489.66 ± 0.1622.12 ± 0.51Neg66.36 ± 0.0047.0 ± 7.21.2 ± 0.2n.d.n.d.
BV6Aspergillus sp. Pos21:39 ± 02:2289.43 ± 0.1021.25 ± 0.3821.91 ± 0.079.62 ± 0.0034.1 ± 0.84.0 ± 0.8n.d.n.d.
BV8A. flavusPos21:28 ± 03:4389.36 ± 0.2522.84 ± 0.2833.73 ± 4.966.51 ± 9.1728.3 ± 4.84.5 ± 0.5n.d.n.d.
BV9A. flavusPos22:59 ± 02:4689.41 ± 0.2921.97 ± 0.5024.52 ± 0.379.63 ± 0.11386 ± 61.20.8 ± 0.3n.d.n.d.
BV10A. flavusPosNeg84.13 ± 0.01632.45 ± 0.6330.12 ± 0.679.63 ± 0.21n.d.n.d.n.d.n.d.
BS1A. parasiticus±26:52 ± 04:5289.43 ± 0.1123.68 ± 0.1238.00 ± 0.065.76 ± 0.2110,051 ± 78415,086 ± 377492,163 ± 466238,362 ± 1360
BS3A. flavusPosNeg84.13 ± 0.00140.04 ± 0.3312.72 ± 0.280.08 ± 0.009.5n.d.n.d.n.d.
BS4A. flavusPos37:13 ± 00:0088.9930.78 ± 0.1636.58 ± 0.365.76 ± 0.1125.9 ± 22.8 ± 0.9n.d.n.d.
BS5A. flavusPosNeg84.29 ± 0.01431.72 ± 0.2729.46 ± 0.3580.63 ± 0.08n.d.n.d.n.d.n.d.
BS6A. flavusNeg19:35 ± 02:1489.51 ± 0.2122.01 ± 0.45Neg64.72 ± 0.5310,011 ± 3222264 ± 173n.d.n.d.
BS7Aspergillus sp. NegNeg84.20 ± 0.11NegNeg64.70 ± 0.42n.d.n.d.n.d.n.d.
BS8A. sojaeNegNeg84.13 ± 0.22Neg34.15 ± 0.379.62 ± 0.11n.d.n.d.n.d.n.d.
BS9A. flavusNeg22:30 ± 02:5489.41 ± 0.2120.92 ± 0.25Neg65.30 ± 0.7423.5 ± 4.71.9 ± 0.4n.d.n.d.
BS10A. flavusNegNeg84.06 ± 0.09NegNeg65.31 ± 0.74n.d.n.d.n.d.n.d.
Strains from Kernel
NP1A. flavusPosNegNeg27.80 ± 11.9433.15 ± 0.4078.28 ± 0.00n.d.n.d.n.d.n.d.
NP2Aspergillus sp. NegNegNeg26.08 ± 0.9422.41 ± 0.1578.51 ± 0.11n.d.n.d.n.d.n.d.
NP3A. flavusNegNegNegNegNegNegn.d.n.d.n.d.n.d.
NP4A. flavusPos29:17 ± 00:2287.34 ± 0.0928.55 ± 0.4623.58 ± 0.1480.74 ± 3.267.3 ± 1.32.5 ± 1.2n.d.n.d.
NP5A. flavusPos22:53 ± 02:0888.97 ± 0.0021.07 ± 1.8516.65 ± 0.1378.66 ± 0.11137 ± 1.6n.d.n.d.n.d.
NP6A. flavusPosNegNeg34.03 ± 0.1729.22 ± 0.4778.58 ± 0.21n.d.n.d.n.d.n.d.
NP7A. flavusNegNegNegNegNegNegn.d.n.d.n.d.n.d.
NP8A. flavusPos22:40 ± 02:4389.22 ± 0.1522.79 ± 0.9517.11 ± 0.1378.80 ± 0.1146.9 ± 4.645.7 ± 0.1n.d.n.d.
Table 2. Strains from chestnuts used in this study to check the specificity of the LAMP assay (Tp, time to positive) and the TaqMan® qPCR assay (Ct). The PCR results reported the presence or absence of a 308 bp amplicon. Aflatoxin production is reported in ppm. Standard deviation was included for each Tp, Ct, melting temperature and for the four aflatoxins. Neg = negative, Pos = positive; n.d.: not detected; Melting = melting temperature.
Table 2. Strains from chestnuts used in this study to check the specificity of the LAMP assay (Tp, time to positive) and the TaqMan® qPCR assay (Ct). The PCR results reported the presence or absence of a 308 bp amplicon. Aflatoxin production is reported in ppm. Standard deviation was included for each Tp, Ct, melting temperature and for the four aflatoxins. Neg = negative, Pos = positive; n.d.: not detected; Melting = melting temperature.
StrainSpeciesPCR ResultLAMP AssayqPCR qPCR with SYBR GreenAflatoxin Production (ppm)
Tp (min:s)Melting (°C)TaqMan®CtMelting (°C)AFB1AFB2AFG1AFG2
Probe (Ct)
XNA. flavusNegNegNegNegNegNegn.d.n.d.n.d.n.d.
AF2A. flavusNegNegNegNegNeg66.68 ± 1.26n.d.n.d.n.d.n.d.
AF3A. flavusNegNeg84.14 ± 0.150Neg30.20 ± 0.079.34n.d.n.d.n.d.n.d.
AF4A. flavusNegNeg84.24 ± 0.089NegNeg67.57 ± 1.05n.d.n.d.n.d.n.d.
AF11A. flavusNegNegNegNegNeg65.93 ± 2.21n.d.n.d.n.d.n.d.
AF12A. flavusNegNeg84.29 ± 0.012NegNeg67.12 ± 0.42n.d.n.d.n.d.n.d.
AF17A. flavusNegNeg84.29 ± 0.006Neg30.58 ± 1.279.04n.d.n.d.n.d.n.d.
AF18A. flavusNegNeg84.30 ± 0.000Neg3765.79 ± 10.22n.d.n.d.n.d.n.d.
AF20A. flavusNegNeg84.20 ± 0.174NegNeg65.49n.d.n.d.n.d.n.d.
AF24A. flavusNegNegNegNeg27.10 ± 0.579.34n.d.n.d.n.d.n.d.
AFCAL2A. flavusPos23:12 ± 00:3089.66 ± 0.0122.92 ± 1.5515.33 ± 0.680.08 ± 0.00152 ± 264.3 ± 0.3n.d.n.d.
AFCAL3A. flavusPos25:09 ± 04:1289.67 ± 0.0126.34 ± 0.5218.95 ± 0.880.23 ± 0.2114.2 ± 2.3n.d.n.d.n.d.
AFCAL4A. flavusPos22:16 ± 03:1289.82 ± 0.2620.44 ± 0.9813.58 ± 0.380.08 ± 0.1116.1 ± 1.4n.d.n.d.n.d.
AFCAL6A. flavusPos21:03 ± 03:1389.51 ± 0.0118.28 ± 0.6711.88 ± 0.380.23 ± 0.112.5 ± 0.2n.d.n.d.n.d.
AFCAL7A. flavusPos22:30 ± 00:3389.67 ± 0.0126.00 ± 1.6718.35 ± 0.080.08 ± 0.00n.d.n.d.n.d.n.d.
AFCAL9A. flavusNegNeg89.36 ± 2.93Neg32.08 ± 3.479.78 ± 0.00n.d.n.d.n.d.n.d.
AFLX1A. flavusNegNeg84.15 ± 0.151Neg30.60 ± 0.179.78150 ± 25.46.0 ± 0.1n.d.n.d.
AFLX2A. flavusNegNeg84.15 ± 0.150NegNeg65.3343.1 ± 12.94.4 ± 0.9n.d.n.d.
AFLX5A. flavusNegNeg84.30 ± 0.006NegNeg64.43 ± 0.32n.d.n.d.n.d.n.d.
AFLX7A. flavusPos22:14 ± 03:2089.56 ± 0.0921.89 ± 0.74NegNegn.d.n.d.n.d.n.d.
AFLX8A. flavusPos21:14 ± 04:3889.61 ± 0.0933.22 ± 0.8331.03 ± 0.379.63 ± 0.00817 ± 128201 ± 39.8n.d.n.d.
AFLX9A. flavusNegNeg84.14 ± 0.158Neg32.61 ± 8.066.07n.d.n.d.n.d.n.d.
AFLX10A. flavusPosNeg84.12 ± 0.10Neg26.9261.31 ± 3.89n.d.n.d.n.d.n.d.
AFSP1A. flavusNegNeg84.29 ± 0.006Neg32.88 ± 2.979.34n.d.n.d.n.d.n.d.
AFSP2A. flavusNegNeg84.35 ± 0.377Neg38.72 ± 1.566.97 ± 0.9566.97 ± 0.955.6 ± 0.7n.d.n.d.
AFSP3A. flavusNegNeg84.20 ± 0.17439.56 ± 0.6528.96 ± 0.079.345.9 ± 0.1n.d.n.d.n.d.
AFSP4A. flavusPos21:38 ± 02:4089.61 ± 0.0825.59 ± 0.4216.34 ± 0.379.79 ± 0.32177.5 ± 8.94.1 ± 1.1n.d.n.d.
AFSP5A. flavusNegNeg84.29 ± 0.142NegNeg63.09 ± 0.84n.d.n.d.n.d.n.d.
AFSP6A. flavusPos20:58 ± 03:1989.62 ± 0.0823.08 ± 0.4013.25 ± 1.079.19 ± 0.2154.7 ± 4.21.9 ± 0.8n.d.n.d.
AFSP7A. flavusPos23:13 ± 01:3689.72 ± 0.2229.83 ± 0.7622.47 ± 2.779.64 ± 0.1111.2 ± 0.5n.d.n.d.n.d.
AFCAL1A. oryzae var effususNegNeg84.29 ± 0.01Neg37.980.97n.d.n.d.n.d.n.d.
AFCAL5A. oryzae var effususNegNeg84.24 ± 0.08NegNeg65.47 ± 0.21n.d.n.d.n.d.n.d.
AFCAL10A. oryzae var effususNegNeg84.29 ± 0.006NegNeg66.07 ± 0.74n.d.n.d.n.d.n.d.
AFLX3A. oryzae var effususNegNeg84.29 ± 0.006Neg31.8779.64n.d.n.d.n.d.n.d.
AFLX4A. oryzae var effususNegNeg84.24 ± 0.080NegNeg64.88 ± 0.21n.d.n.d.n.d.n.d.
MCAL1A. tamariiNegNeg84.14 ± 0.006Neg28.67 ± 0.379.63n.d.n.d.n.d.n.d.
MCAL2A.tamariiNegNeg84.29 ± 0.004Neg34.99 ± 4.579.48 ± 9.38n.d.n.d.n.d.n.d.
AFCAL8A. parasiticusPos24:10 ± 04:2089.67 ± 0.1521.72 ± 1.3316.52 ± 0.379.93 ± 0.21352 ± 59.468.2 ± 27.5944 ± 169128 ± 25.9
AFLX6A. parasiticusPos20:58 ± 03:3189.57 ± 0.0822.66 ± 0.8718.58 ± 0.379.44 ± 0.28387 ± 36.742.1 ± 8.8662 ± 8.8108 ± 10.1
AFSP8A. toxicariusPos27:49 ± 00:3789.59 ± 0.1128.56 ± 0.4119.24 ± 0.179.64 ± 0.21910 ± 78.2235 ± 24.21142 ± 78.0392 ± 59.2
Table 3. PCR, LAMP and qPCR primers designed and used in this study. Hyphens denote the division between F1 and F2 and B1 and B2 of the FIP and BIP (forward and backward inner) primers.
Table 3. PCR, LAMP and qPCR primers designed and used in this study. Hyphens denote the division between F1 and F2 and B1 and B2 of the FIP and BIP (forward and backward inner) primers.
NameSequenceTechnique
FIP5′-CGAGAATGCTAAGACAGCTGCAAATTAG-CTGGTGGCATATAATGCCGC-3′LAMP
BIP5′-CGCCGTCGGGATCGAATGGAG-CAAAGTTAGCCGGGTACTATCG-3′
loopF5′-GCCACCTGAAAGACGTCCAGCA-3′
loopB5′-AAGCGTTAAGCAGAACAGGTAACG-3′
F35′-CGGAGTACATGGACCAGGTG-3′
B35′-AGAAGACAGAAGAAAAGGTAAGA-3′
COX-FIP5′-ATGGATTTGRCCTAAAGTTTCAGGGCAGGATTTCACTATTGGGT-3′LAMP [52]
COX-BIP5′-TGCATTTCTTAGGGCTTTCGGATCCRGCGTAAGCATCTG-3′
COX-loopF5′-ATGTCCGACCAAAGATTTTACC-3′
COX-loopB5′-GTATGCCACGTCGCATTCC-3′
COX-F35′-TATGGGAGCCGTTTTTGC-3′
COX-B35′-AACTGCTAAGRGCATTCC-3′
AflF5′-CTGGACGTCTTTCAGGT G-3′qPCR
AflR5′-CTTCTCCATTCGATCCCG-3′
Afl-probe5′6-Fam-AM- CCTA ATT TGC AGC TGT CTT AGC ATT CTC- ZNA-4-Eclipse-3′
AfF5′- CCAGGTCTCGATCCAGAGGC -3′PCR
AfR5′- CCGGAACGTGAACTTATGTGCAC -3′
Table 4. Sensitivity testing of the LAMP assay using the DNA of two strains (A. flavus strain FS7 and A. parasiticus strain AFLX6) at different concentrations. Standard deviation was included for each Tp and melting temperature. Tp = time to positive; Melting = melting temperature; n°P = number of positive technical replicates out of 3.
Table 4. Sensitivity testing of the LAMP assay using the DNA of two strains (A. flavus strain FS7 and A. parasiticus strain AFLX6) at different concentrations. Standard deviation was included for each Tp and melting temperature. Tp = time to positive; Melting = melting temperature; n°P = number of positive technical replicates out of 3.
A. flavus FS7A. parasiticus AFLX6
DNA ConcentrationTp (min:s)Melting (°C)n°PTp (min:s)Melting (°C)n°P
10–50 ng/µL22:49 ± 00:5289.51 ± 0.00320:57 ± 03:3189.57 ± 0.083
1–10 ng/µL25:04 ± 00:2189.61 ± 0.09322:50 ± 02:0789.67 ± 0.153
100–999 pg/µL31:37 ± 07:3989.52 ± 0.26328:22 ± 03:0889.61 ± 0.083
10–99 pg/µLNegativeNegative0NegativeNegative0
Table 5. Results of the LAMP assay using crude DNA extracts and DNA extracted with a commercial kit of artificially inoculated hazelnuts at different concentrations with the conidia of three Aspergillus strains (AFSP1, non-aflatoxigenic A. flavus, AFSP4, aflatoxigenic A. flavus and AFLX6, aflatoxigenic A. parasiticus). The crude DNA extraction reliability was evaluated using the COX LAMP assay as internal control. The TaqMan® qPCR was used to calculate the number of cells present in each hazelnut. For each conidial concentration, 3 hazelnuts were inoculated. Tp: LAMP time to positive; Melting: LAMP melting temperature; Neg: negative. n°P = number of positive technical replicates out of 3.
Table 5. Results of the LAMP assay using crude DNA extracts and DNA extracted with a commercial kit of artificially inoculated hazelnuts at different concentrations with the conidia of three Aspergillus strains (AFSP1, non-aflatoxigenic A. flavus, AFSP4, aflatoxigenic A. flavus and AFLX6, aflatoxigenic A. parasiticus). The crude DNA extraction reliability was evaluated using the COX LAMP assay as internal control. The TaqMan® qPCR was used to calculate the number of cells present in each hazelnut. For each conidial concentration, 3 hazelnuts were inoculated. Tp: LAMP time to positive; Melting: LAMP melting temperature; Neg: negative. n°P = number of positive technical replicates out of 3.
Conidial ConcentrationCrude DNA ExtractionCOX AssayCommercial Kit DNA Extraction
Tp (min:s)Melting (°C)n°PCells/Hazelnut (qPCR)Tp (min:s)Melting (°C)Tp (min:s)Melting (°C)n°PTaqMan (Ct)n°PCells/µL (qPCR) (Range of Cells/uL)
AFSP1 (non-aflatoxigenic A. flavus)
3.09 × 105NegNeg0Neg05:43 ± 00:0185.21 ± 0.00NegNeg0Neg0Neg
3.09 × 106NegNeg0Neg05:45 ± 00:4985.22 ± 0.00NegNeg0Neg0Neg
3.09 × 107NegNeg0Neg05:09 ± 00:4285.36 ± 0.00NegNeg0Neg0Neg
AFSP4 (aflatoxigenic A. flavus)
3.03 × 10545:07 ± 12:5888.91 ± 0.1531.7 × 10509:16 ± 05:3385.27 ± 0.0938:34 ± 06:3988.51 ± 0.28930.3 ± 3.295.8 × 103 (2.6 × 104–53)
3.03 × 10639:31 ± 00:0089.21 ± 0.0011.1 × 10405:39 ± 00:4485.21 ± 0.0042:53 ± 04:4788.65 ± 0.09931.6 ± 1.79355 (640–16)
3.03 × 10743:38 ± 05:2589.01 ± 0.1836.3 × 10508:18 ± 05:3685.27 ± 0.0932:59 ± 02:2688.51 ± 0.30926.7 ± 1.792.1 × 104 (6.2 × 104–2.8 × 103)
AFLX6 (aflatoxigenic A. parasiticus)
1.79 × 10544:35 ± 00:0088.77 ± 0.0015.6 × 10505:08 ± 00:2685.14 ± 0.1034:36 ± 04:0788.46 ± 0.12632.3 ± 5.091.9 × 104 (7.9 × 104–8)
1.79 × 10645:26 ± 07:2088.70 ± 0.7528.4 × 10406:37 ± 00:4785.12 ± 0.1739:56 ± 03:3188.46 ± 0.26833.2 ± 3.292.8 × 103 (9.5 × 103–60)
1.79 × 107NegNeg02.4 × 10605:48 ± 00:3985.22 ± 0.1533:36 ± 02:1888.48 ± 0.18927.6 ± 2.098.2 × 104 (1.9 × 10–4.0 × 103)
Table 6. Naturally contaminated hazelnuts used to evaluate the reliability of the TaqMan® qPCR assay. The results listed the Ct, the number of cells in each hazelnut calculated with the TaqMan® qPCR and the number of positive replicates out of three. Standard deviation was included for each Ct and for the approximate number of cells. n°P = number of positive technical replicates out of 3 with the qPCR.
Table 6. Naturally contaminated hazelnuts used to evaluate the reliability of the TaqMan® qPCR assay. The results listed the Ct, the number of cells in each hazelnut calculated with the TaqMan® qPCR and the number of positive replicates out of three. Standard deviation was included for each Ct and for the approximate number of cells. n°P = number of positive technical replicates out of 3 with the qPCR.
Hazelnut NumberCtNumber of Cellsn°P
138.83 ± 0.8320 ± 113
738.39 ± 0.6525 ± 112
839.35 ± 0.2712 ± 22
1037.05 ± 0.2062 ± 83
2239.15 ± 0.1914 ± 22
2339.7481
3237.35 ± 1.5566 ± 602
3337.10 ± 0.2360 ± 103
3434.35 ± 0.61439 ± 1753
3534.40 ± 0.48204 ± 683
3635.44 ± 0.97227 ± 1693
3835.52 ± 1.17219 ± 1393
3936.25 ± 2.09182 ± 1633
4036.61 ± 0.6690 ± 383
4135.52 ± 0.28183 ± 363
4236.04 ± 1.24165 ± 1533
4335.91 ± 0.42142 ± 403
4539.14 ± 0.7916 ± 83
4738.71191
Table 7. Summary of the comparison of the diagnostic sensitivity and specificity of the end-point PCR, LAMP and qPCR with SYBR Green and TaqMan ® probe using DNA of 97 Aspergillus spp.
Table 7. Summary of the comparison of the diagnostic sensitivity and specificity of the end-point PCR, LAMP and qPCR with SYBR Green and TaqMan ® probe using DNA of 97 Aspergillus spp.
PCRLAMPqPCR with SYBR GreenqPCR with TaqMan® Probe
+LR4.85221.365.59
Approximate 1 in 1.2 positive results are aflatoxigenic strainsApproximate 1 in 1.0 positive results are aflatoxigenic strainsApproximate 1 in 1.8 positive results are aflatoxigenic strainsApproximate 1 in 1.2 positive results are aflatoxigenic strains
−LR0.160.110.60.13
Approximate 1 in 1.2 negative results are non-aflatoxigenic strainsApproximate 1 in 1.1 negative results are non-aflatoxigenic strainsApproximate 1 in 1.5 negative results are non-aflatoxigenic strainsApproximate 1 in 1.1 negative results are non-aflatoxigenic strains
Diagnostic sensitivity87.2%89.4%71.7%89.4%
Diagnostic specificity82.0%96%47.1%84%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ortega, S.F.; Siciliano, I.; Prencipe, S.; Gullino, M.L.; Spadaro, D. Development of PCR, LAMP and qPCR Assays for the Detection of Aflatoxigenic Strains of Aspergillus flavus and A. parasiticus in Hazelnut. Toxins 2020, 12, 757. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12120757

AMA Style

Ortega SF, Siciliano I, Prencipe S, Gullino ML, Spadaro D. Development of PCR, LAMP and qPCR Assays for the Detection of Aflatoxigenic Strains of Aspergillus flavus and A. parasiticus in Hazelnut. Toxins. 2020; 12(12):757. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12120757

Chicago/Turabian Style

Ortega, Sara Franco, Ilenia Siciliano, Simona Prencipe, Maria Lodovica Gullino, and Davide Spadaro. 2020. "Development of PCR, LAMP and qPCR Assays for the Detection of Aflatoxigenic Strains of Aspergillus flavus and A. parasiticus in Hazelnut" Toxins 12, no. 12: 757. https://0-doi-org.brum.beds.ac.uk/10.3390/toxins12120757

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