Skip to main content
  • Research article
  • Open access
  • Published:

Oligonucleotide microarray for the identification of potential mycotoxigenic fungi

Abstract

Background

Mycotoxins are secondary metabolites which are produced by numerous fungi and pose a continuous challenge to the safety and quality of food commodities in South Africa. These toxins have toxicologically relevant effects on humans and animals that eat contaminated foods. In this study, a diagnostic DNA microarray was developed for the identification of the most common food-borne fungi, as well as the genes leading to toxin production.

Results

A total of 40 potentially mycotoxigenic fungi isolated from different food commodities, as well as the genes that are involved in the mycotoxin synthetic pathways, were analyzed. For fungal identification, oligonucleotide probes were designed by exploiting the sequence variations of the elongation factor 1-alpha (EF-1 α) coding regions and the internal transcribed spacer (ITS) regions of the rRNA gene cassette. For the detection of fungi able to produce mycotoxins, oligonucleotide probes directed towards genes leading to toxin production from different fungal strains were identified in data available in the public domain. The probes selected for fungal identification and the probes specific for toxin producing genes were spotted onto microarray slides.

Conclusions

The diagnostic microarray developed can be used to identify single pure strains or cultures of potentially mycotoxigenic fungi as well as genes leading to toxin production in both laboratory samples and maize-derived foods offering an interesting potential for microbiological laboratories.

Background

Mycotoxins are fungal toxins which pose a threat to human, animal and plant health. These toxins can cause acute or chronic toxicity in humans and animals that eat contaminated foods or crops, depending on the quantities produced and consumed [1]. It is estimated that 25% of all food commodities produced on earth are contaminated with mycotoxins due to the fact that fungi develop on these commodities [2]. A study done in South Africa by Rabie et al. [3] showed that mycotoxins such as aflatoxins, beauvericin, deoxynivalenol, moniliformin, trichothecene and zearalenone are contaminants of food commodities. The most known and studied group of mycotoxins in South Africa are fumonisins which have been associated with oesophageal cancer in humans and the cause of leucoencephalomalacia (LEM) in horses, mules and donkeys [4]. It is thus necessary to eliminate or reduce the presence of mycotoxins in the food chain.

An important step in controlling contaminants in the food production chain is by identifying food-borne fungi. The conventional methods used for the detection of fungal contamination are based on phenotypic and physiological characteristics that make use of standard culture and biochemical/serological tests. However, these methods are very time-consuming, laborious and do not detect mycotoxins. Recently, a variety of molecular methods have been used for fungal pathogen identification and for their potential to produce mycotoxins [5]. Molecular methods were used for Aspergillus species differentiation using Southern blot hybridization assays [6] and PCR-based restriction fragment length polymorphisms [7]. Most assays that have been developed included PCR-based methods that exploited the highly conserved ribosomal RNA gene sequences for the design of species-specific primers [8] as well as generic PCR detection assays developed for genes involved in the biosynthesis of some mycotoxins [9, 10]. Although these assays are an improvement compared to conventional methods, the overall throughput is still limited. Only a limited number of diagnostic regions can be identified for a single organism at a time. If all potentially mycotoxigenic fungi must be included, these assays become laborious and expensive. The use of integrated platforms that combine identification and typing methods for several fungi would facilitate the rapid and accurate identification of possible mycotoxigenic fungi in food commodities.

The microarray technique allows the rapid and parallel characterization of a range of organisms and has the intrinsic ability to perform multiplexed and low-volume biological assays. This technique has been increasingly used for diagnostic purposes as it has the ability to detect more than one parameter at a time [11, 12]. Leinberger et al. [13] exploited the polymorphisms of the internal transcribed regions in the ribosomal RNA cassette for the microarray-based detection and identification of Candida and Aspergillus species. In a similar experiment, DeSantis et al. [14] generated a 62358-probe oligonucleotide of small subunit ribosomal RNA (ssu rRNA) for the detection of 18 different orders of microbes from environmental samples and novel variants exhibiting mutations in their ssu rRNA. Microarrays have also been successfully used to study the expression levels of mycotoxin gene clusters. Schmidt-Heydt and Geisen [15] developed a microarray which contained oligonucleotide probes for the biosynthesis pathways of fumonisin, aflatoxin, ochratoxin, patulin and trichothecene. These studies showed that it is possible to use a range of organisms for hybridization to an array containing many different oligonucleotide probes. The microarray technique is thus analogous to performing many PCR reactions and hybridization reactions at the same time and has the advantage of being versatile [16].

The aim of this study was to develop a diagnostic microarray for the identification of single strains of food-borne fungi that are most prevalent in South African food commodities, and to detect the ability of these fungi to produce mycotoxins in laboratory and food samples. A total of 40 food-borne fungi isolated from different foods that belong to the genera Alternaria, Aspergillus, Bipolaris, Claviceps, Curvularia, Diplodia, Drechslera, Eurotium, Fusarium, Penicillium and Pithomyces, were used. For fungal discrimination, the polymorphisms of the internal transcribed spacer (ITS) regions and the elongation factor 1- alpha (EF-1 α) gene were exploited for the design of the oligonucleotide probes. The specificity of a probe was increased in some instances by substituting an oligonucleotide with a high affinity DNA analogue known as locked nucleic acid (LNA). A locked nucleic acid nucleotide analogue consists of a 2'-O,4'-C methylene bridge and locks the LNA structure into a rigid bicyclic formation and displays unprecedented hybridization affinity towards complementary DNA and RNA [17]. It is most disruptive, and thus gives a better signal, in a centre position. For the detection of fungi that can produce mycotoxins, oligonucleotide probes for the genes leading to mycotoxin production were selected from public databases and included in the oligonucleotide array. The combination of ITS, EF-1 α and mycotoxin genes on the same array was evaluated for the potential of the array to identify the forty fungal isolates and the genes involved in pathways leading to toxin production.

Results

Probe design

A 96-probe oligonucleotide microarray was constructed for the simultaneous detection and identification of potentially mycotoxigenic fungi. Probes for the array were designed by exploiting the polymorphisms of the internal transcribed spacer (ITS) regions of the rRNA complex. Amplification of fungal DNA with the universal fungal primers ITS1 and ITS4 and subsequent sequence analysis allowed the differentiation of most of the fungal species studied. Several unique polymorphisms (sequence data can be found in GenBank with accession numbers [GenBank:FJ864706, GenBank:FJ864709, GenBank:FJ864710, GenBank:FJ864708, GenBank:FJ864711, GenBank:FJ864703, GenBank:FJ864704, GenBank:FJ864705, GenBank:FJ864707, and GenBank:FJ864712]) could be identified within the PCR products generated for each fungal species. However, amplification of the Fusarium species showed no significant differences between the sequences of the PCR products generated with the ITS primers. Therefore, the elongation factor 1-alpha (EF-1 α) gene was used for the identification of polymorphisms in Fusarium species and for the design of unique species- or genus-specific probes. A total of 38 probes could be designed or identified for Alternaria, Aspergillus, Penicillium and Stenocarpella species from the ITS regions and 22 probes could be designed or identified for the Fusarium species from the EF-1 α gene (Table 1). This probe set was then extended by searching public databases for additional probes for the ITS regions and the EF-1 α gene. No unique probes could be designed for Drechslera species, Eurotium chevalieri, Fusarium sambucinum, F. semitectum, Penicillium funiculosum, P. rugulosum and Pithomyces chartarum. The Fusarium and Penicillium strains share many sequence similarities with the other species used in this study. This rendered the development of species-specific oligonucleotide probes more difficult. For the strains Pithomyces and Eurotium no unique polymorphisms could be identified that could be used for the design of unique probes.

Table 1 Probe sequences and names of species- and toxin- specific genes for different fungal isolates

The public databases were also used to identify oligonucleotide probes specific for genes leading to toxin production. The probes selected for the biosynthesis genes leading to the regulation of fumonisins, aflatoxin, nivalenol, deoxynivalenol and tricothecenes were 18 - 22 nucleotides in length. A total of 23 toxin-specific probes were identified and spotted onto the array. The oligonucleotide sequences and their specificities are shown in Table 1. Most probes used for the final array construction were oligonucleotide probes identified in public databases as the probe sequences were diverse and minimal cross-hybridization was obtained. Some sequence data is available upon request.

Optimization of labeling and hybridization conditions

To avoid amplification bias and to get a more uniform genetic locus representation, targets were labeled using a random approach that does not involve amplification. All labeled target DNA positively hybridized to the array (Figure 1) showing fluorescent net signal intensities ranging from 2000 to 6000 intensity units demonstrating efficient hybridization of the target DNA. The hybridization conditions were further tested to get the optimal discrimination of target species and genes leading to toxin production without having unspecific signal intensities by determining the optimal PCR annealing temperature for fungal DNA using the probes in Table 1. Aspergillus clavatus and A. versicolor were used for this purpose as they showed cross-hybridization to other species-specific probes in the initial experiment. This was expected as the ITS region of both species are very similar. An increase in hybridization temperature from 42°C to 53°C showed that there is nearly no cross-hybridization between these two species and there was no decrease in net signal intensity (results not shown). Although the ITS sequences are quite similar for both fungal species, high hybridization efficiencies were obtained with net signal intensities of about 2000 signal units for A. clavatus and of about 3500 signal units for A. versicolor (Figure 2A). In general, it was also observed that the optimal probe annealing temperatures for PCR amplifications was about 5°C higher than the optimal probe hybridization temperature (results not shown). The probes and their optimal annealing temperatures are listed in Table 1.

Figure 1
figure 1

Sections of fluorescent images showing DNA hybridized to the array. Sections of fluorescent images after hybridization of target DNA to the diagnostic array. A. (Top) Hybridization profile of Aspergillus versicolor; (Middle) Penicillium corylophilum; (Bottom) P. expansum. B. The arrangement of a few oligonucleotide probes within the indicated fields of a section of the array. Oligonucleotide probe names were used to indicate the field. Each column represents four replicates of the same spot.

Figure 2
figure 2

Relative intensities of hybridized DNA. Relative intensities after hybridization of labeled target DNA to the array. Each experiment was done in triplicate and the medians and their standard deviations were calculated for each spot on the array. Only positive hybridization results are shown. A. Relative intensities of fungal strains hybridizing to probes designed from the internal transcribed (ITS) regions of Alternaria, Aspergillus, Penicillium and Stenocarpella species. B. Relative intensities of Fusarium species hybridizing to probes designed from the elongation factor (EF). C. Relative intensities of Alternaria, Aspergillus, Penicillium and Stenocarpella species hybridizing to their relevant mycotoxin genes. D. Relative intensities of Fusarium species hybridizing to their relevant mycotoxin genes.

Specificity and functionality of the microarray

The specificity of the array was tested by using the forty precharacterized fungal isolates listed in Table 2. The hybridization of fungal isolate to the array gave insight into the affinity of test probes for their correct target and the effect of multiple versus single diagnostic probes/species. The hybridization of each fungal isolates for 16 - 24 hours at 53°C resulted in different hybridization patterns for the different fungal strains (Figure 1) with relative intensities indicating the level of hybridization of each target to the probe (Figure 2). Thirty-two test samples showed high affinity for their probes producing a best match result. It was possible to positively identify the test organisms with at least one probe due to the presence of multiple diagnostic probes with fluorescent net signal intensities ranging from 2992 to 6000 intensity units. SNR values obtained from the relative intensities of hybridized DNA indicated in the graph, gave a clear indication whether a spot was present (SNR>/= 3.0) or absent (SNR<3.0). Weak cross-hybridization was observed for Aspergillus clavatus and A. niger, but these fungal isolates could be positively identified due to the multiple probes on the array. Although the multiple probes per species used for the array construction showed big differences in hybridization efficiencies with some probes showing no hybridization, at least one oligonucleotide showed high hybridization efficiency for most of the fungal isolates tested and could be used for species- or toxin-specific gene identification. Eight species could not be positively identified as they did not reveal specific hybrization patterns (Table 3).

Table 2 Fungal cultures used in this study, their potential mycotoxins and the host of the fungus
Table 3 Fungal species identified (+) or not identified (-) with diagnostic chip developed

The three controls (ITS1, ITS3, ITS4) which were specific for universal fungal sequences served as internal standards to ensure that the parameters (labelling and hybridization) were similar across experiments. A similar intensity of controls across slides indicated that the relative signal intensities of probes are also similar across slides. Further, some probes in this study were modified to contain locked nucleic acids (LNAs) in at least two selected single nucleotide polymorphisms (SNP) sites per fragment. SNP's were found to be most effective, and thus gave better signal, if they were in a centre position. A probe with multiple polymorphisms along the probe length, regardless of position or modification at the polymorphic site, showed less cross-hybridization (results not shown) which is consistent with the data obtained by You et al. [18].

The functionality of the microarray was tested by hybridizing precharacterized fungal isolates to the array. Twenty-five fungal isolates were characterized for the presence of mycotoxin genes by growing them at 25°C for 1 week, extracting genomic DNA and PCR-amplified the DNA of each individual fungal isolate using the toxin-specific oligonucleotide probes that were used for array construction. Different species showed different amplifications of toxin-producing genes (Table 4). These results indicated which fungal isolates have the potential to produce mycotoxins and hybridized to probes specific for genes leading to toxin production on the array. The amplicons obtained were consistent with the signal intensities obtained when samples were hybridized to the array (Figure 2C-D). The microarray chip developed was also tested for its ability to detect genes leading to mycotoxin production without any knowledge about the identity of the fungal isolate. In this study, Fusarium anthophilum was used to test this approach as no species-specific probes were present on the slide. The hybridization of this fungus to the fum5 F and fum5 R probes (Figure 2C-D) indicated that the fungus is able to produce fumonisins confirming that mycotoxin-producing genes can be detected. It should be noted that the presence of a gene in the genome does not mean that a gene is transcribed and expressed.

Table 4 Fungal species screened and scored for for presence (+) or absence (-) of mycotoxin genes with PCR

Validation of the array

The performance and reproducibility of the array was tested starting from independently extracted fungal DNA from eight blind fungal samples that were hybridized to the array. Binary scores obtained from the array were compared to the binary scores from replicate experiments. Repeatability of the binary scores obtained from the hybridizations from replicate experiments of the same fungi were on average 95%. The results obtained were also compared in each case to the identity obtained for the same culture grown by standard laboratory procedures and to the correlation of the PCR product amplified from the same sample with the positively identified oligonucleotide probes. The same procedure was followed for the mycotoxin biosynthesis genes. The identities of the amplicons and the identities of the fungi obtained by standard methods showed that the array was able to identify the fungi and mycotoxin genes correctly; seven of the eight fungal isolates could be identified up to the species level (Table 3). Fusarium sambucinum could not be identified to species level due to the absence of species-specific signals. In all cases the genes leading to mycotoxin production could be identified.

Discussion

The identification and detection of fungi has become increasingly dependent on molecular characterization. Methods such as Southern blot hybridization assays, restriction fragment length polymorphism analysis and PCR-based assays exploiting the internal transcribed spacer (ITS) and elongation factor 1-alpha (EF-1 α) regions are all effective for the detection and identification of food-borne fungi. However, all these methods can identify only a single organism at a time. Suitable detection methods, anticipating mycotoxin risks, are needed to ensure a safe food production chain and eliminate the risk factors. Oligonucleotide microarrays have a high multiplexing capacity and have proved to be an efficient approach to overcome these limitations. This technology offers an identification process based on sequence confirmation through hybridization [16] and has the ability to analyze many samples simultaneously.

In this study, a microarray chip was developed for the identification of fungal isolates, as well as biosynthesis genes of the most important mycotoxins that can threaten the safety and quality of food products especially those derived from maize. For the design of genus- and species-specific probes the ITS regions of the rRNA gene cassette were exploited. These coding regions show a high degree of variation [19] and analysis of the fungal ITS alignments revealed significant differences among the different fungi. However, analysis of the ITS regions of Fusarium species showed that they have similar sequences which could have cross hybridized on the array, making it non-specific. Kane et al. [20] found that in 50mer oligonucleotide arrays, cross-hybridization occurred between fragments of relatively low sequence similarity. The highly repetitive DNA content of plant genomes resulted in cross-hybridization of DNA fragments to printed-probe DNA and the overall spot intensity of many probes was increased. Therefore, the EF regions were used for the design of species-specific probes for Fusarium species. For some probes with similar sequences the chances of cross hybridization were minimized by substituting a single oligonucleotide in the probe sequence using a high affinity DNA analogue known as locked nucleic acid (LNA) at three specific points to increase the specificity and the Tm of a probe. The LNAs were inserted at a single nucleotide polymorphism (SNP) site for improved performance of the probe. Letowski et al. [21] found that probes containing polymorphisms toward the centre of the probe showed a higher discrimination power. If LNAs are to be included then they must be inserted in a triplicate series around the centre of the probe. Further, G-T mismatch sites must be avoided and should preferably be inserted at sites where adenine is the identity of the base [18].

Cross hybridization has also been reported in several microarray-based species detection studies where single regions were used for identification. Anthony et al. [22] found that in oligonucleotide arrays, cross-hybridization occured between Listeria species and it was necessary to include additional probes to the array. In a similar study done by Volokhov et al [23], E. coli and Salmonell a isolates produced indistinguishable hybridization profiles when single probes were used. However, they showed that multiple probes improve the sensitivity of the array when compared with the single diagnostic probes that could be unsuitable for a group of closely related organisms. In this study, the probes spotted onto the array were a mixture of single and multiple probes for each species that were either genus-, species-specific or specific for genes leading to toxin production. When multiple probe sequences were used the discriminatory power of the array increased as a sample hybridized to at least one probe of the multiple probes on the array. In addition, probes for the array construction were designed around a Tm of 56°C so that all probes would hybridize under similar conditions. Considering that microarray experiments are non-equilibrium measurements, it is desirable that microarray probes exhibit uniform thermodynamic properties, which many probe design tools aim to achieve by demanding a narrow distribution of the probe melting temperature Tm [24]. The fungi hybridizing to the diagnostic array may, however, represent a taxon or haplotype that was not included in the array design. In some of the species complexes included in this study several haplotypes of ITS1 and/or TEF1a genes may be found suggesting that probes my fail to detect some of the haplotypes. The cross hybridization that was observed between A. clavatus and A. niger indicates that more strains need to be studied and additional probes still need to be designed to discriminate between these two species. This also applies to the eight fungal species that could not be identified to species level.

The random labeling strategy used in this study was applied to diminish secondary structures [25] and to have an efficient target. Previous studies suggested that amplification products of large samples resulted in poor hybridization and target PCR amplification resulted in amplification bias [26]. Although high levels of amplification are desirable for PCR assays, this feature is less critical for microarrays as only limited probe is available on the array surface [16]. As target genomic DNA was not a limiting resource in this study, a random approach that omits the target amplification step prior to DNA hybridization proved to be efficient for the sensitive detection of fungi. This approach ensured that there is an equal amount of target sequences available for dye coupling and thus their representation on the array was balanced. This makes the microarray an attractive tool for single strain fungal infections compared to morphological identification. Zheng et al [27] identified the three fungal pathogens, Candida, Cryptococcus neoformans and Aspergillus directly from 27 clinical specimens using a microarray. However the ability of the present microarray to reliably detect mixed infections and single copy genes such as TEF1a was not established. It is also likely that in a sample containing multiple fungi, the fast-growing fungi are extracted in greater concentrations than the slow-growing fungi making the identification of all the fungi present in the sample not possible.

The microarray developed was also evaluated for its ability to detect genes leading to toxin production without prior knowledge of the fungus that produced it. Determination of toxin producing genes is often of a greater concern than the identification of the exact fungal species. Although our understanding of the biosynthesis of mycotoxins is incomplete several genes have been identified. Often more than one gene plays a key role in the biosynthetic pathway and it is important to include as many genes as possible on the microarray chip for proper identification of toxin-producing fungi. In this study, two unknown Fusarium species that show no species-specific signals on the array were identified to genus level by their positive hybridization to the mycotoxin-producing genes tri5, tri13 and fum5 F and fum5 R, respectively. However, based only on the hybridization signal it was not possible to predict whether the respective mycotoxin was produced. This could have been achieved if cDNA would have been used as a target in the array hybridization where differentially expression of mycotoxin genes would have indicated mycotoxin production. Schmidt-Heydt and Geisen [15] used RNA to detect the activation of gene clusters under conditions conducive for the biosynthesis of trichothecenes, fumonisin, ochratoxin, aflotoxin and patulin. However, they found that the biosynthesis of secondary metabolites, like mycotoxins, is dependent on environmental conditions like substrate, pH, temperature and water activity [28] and thus mycotoxins are not always expressed.

Conclusions

With the multiplexing capacity as one of the important features of microarrays, the method developed in the present study can be used to detect more than one parameter at a time, namely fungal species and genes involved in pathways leading to toxin production. A total of 32 fungi could be identified and their potential to produce mycotoxins could be determined. This study describes the omission of the target amplification step of target DNA prior to hybridization in a DNA-based microarray experiment. The results indicated that the random labeling technique could provide enough labeled target DNA for the direct detection of a single fungal infection from infected maize kernels using the microarray In the long term, the developed microarray chip could be used to hybridize DNA and cDNA labeled with different Cy dyes for the simultaneous detection of fungal identity and toxin involved genes. The genomic DNA would determine the fungal identity and the cDNA would determine whether genes for mycotoxin biosynthesis are expressed. The cDNA approach can also be useful to determine which gene clusters are expressed under conditions conducive for the biosynthesis of trichothecenes, fumonisin, ochratoxin, aflatoxin and patulin.

Methods

Fungal cultures and DNA extraction

A total of forty food-borne fungi posing a health threat in South Africa were obtained from the Agricultural Research Council culture collection (ARC), Pretoria, South Africa and are listed in Table 1. Up to two isolates of each taxon were used depending on availability. Further, eight blind samples were taken at random from the forty fungi to validate the array. Fungal strains were grown on 1.5% malt extract agar at 25°C for 1-2 weeks. Total genomic fungal DNA was extracted following the DNA extraction protocol described by Raeder and Broda [29] and column-purified using the QIAquick PCR Purification Kit (QIAGEN). Total genomic DNA of inoculated maize kernels was isolated by the same protocol.

Design and identification of unique oligonucleotide probes

For DNA sequencing and subsequent oligonucleotide design, the full ITS regions of all 40 fungi were amplified using universal fungal primers ITS1 (5'-TCCGTAGGTGAACCTGCGG-3') and ITS4 (5'-TCCTCCGCTTATTGATATGC-3') as described by White et al [19]. The PCR amplifications were performed in a 25 μl volume containing 0.4 μM of each universal primer, 1.5 mM MgCl2, 0.2 mM of each dNTP, 0.5 U Taq polymerase and 1 × reaction buffer (Bioline) and 8 ng template DNA. The PCR amplification consisted of 35 cycles of denaturation at 94°C for 30 sec, primer annealing at 50°C for 45 sec, and primer extension at 72°C for 1 min; an initial denaturation step of 94°C for 5 min, and a final extension at 72°C for 5 min. Amplicons were then sent for sequencing to Inqaba Biotec (Pretoria, South Africa). Sequenced fragments were aligned using ClustalX [30] in order to identify polymorphisms. For species that showed no significant polymorphisms of the internal transcribed spacer (ITS) regions of the rRNA complex, a different conserved region, namely the elongation factor 1- alpha (EF-1 α) gene, was amplified using primers EF1 and EF2 [31]. The monoplex PCR amplifications were performed in a 20 μl volume containing 0.4 mM of each forward and reverse primer, 0.25 mM of each dNTP, 1× reaction buffer (Bioline), 0.5 U Taq polymerase (Bioline) and 6 ng template DNA. The PCR amplification consisted of 30 cycles of denaturation at 94°C for 30 sec, primer annealing at 57°C for 45 sec, and primer extension at 72°C for 1 min; an initial denaturation of 94°C for 5 min, and a final extension of 72°C for 7 min. The amplified fragments were then column-purified (QIAquick PCR purification Kit, QIAGEN GmbH) and sent for sequencing to Inqaba Biotec (Pretoria, South Africa). The ITS and EF-1 α sequences were then submitted to GenBank [GenBank:FJ864706, GenBank:FJ864709, GenBank:FJ864710, GenBank:FJ864708, GenBank:FJ864711, GenBank:FJ864703, GenBank:FJ864704, GenBank:FJ864705, GenBankFJ864707, and GenBank:FJ864712] and served as targets for the design of multiple probes for each species which are able to discriminate between the forty fungal isolates. The sequences of the conserved regions were aligned using the ClustlX software [30], manual adjustments were made and areas of interspecies variation were identified. These regions were used for the design of genus- and species-specific probes of various lengths (14-25 bases) and within a narrow range of the melting temperature of 56°C (± 5°C). All oligonucleotide probes were designed using the Primer Designer 4 package (Version 4.2, Scientific and Educational Software, Cary, NC). The probe set was then extended by searching public databases (NCBI and EMBL) for genus-or species-specific oligonucleotide probes. The specificity of each oligonucleotide was assessed by conducting BLAST searches and only unique oligonucleotide probes were chosen to be printed onto the array. Probes with similar sequences were then made more stringent by the insertion of LNA (locked nucleic acid) to increase the specificity of each oligonucleotide and to get a set of probes with similar hybridization efficiencies (Table 2). Multiple probes were identified, if possible, for each species or toxin. The probes identified and designed were synthesised by Inqaba Biotech, Pretoria (Pretoria, South Africa).

In addition, the public databases were used to identify toxin-specific probes for genes leading to toxin production for each of the 40 fungi. To test the optimal annealing temperature for array hybridization, monoplex PCR amplifications were carried out for all the probes identified. The PCR amplifications were performed in a 25 μl volume containing 0.4 μM of each oligonucleotide, 1.5 mM MgCl2, 0.2 mM of each dNTP, 0.5 U Taq polymerase and 1 × reaction buffer (Bioline) and 5 ng template DNA. The PCR amplification consisted of 30 cycles of denaturation at 94°C for 30 sec, oligonucleotide specific annealing temperatures varying from 55°C to 60°C for 45 sec depending on the primer used, and extension at 72°C for 1 min; an initial denaturation step at 94°C for 5 min, and a final extension step at 72°C for 5 min. Aliquots of amplicons were resolved on 1% agarose gels.

Array construction

Arrays were constructed from 86 uniquely designed species- and toxin-specific oligonucleotide probes. Equal volumes (10 μl each) of 100 pmol/ml oligonucleotide and 100% DMSO were transferred into a 384-well plate (Amersham PharmaciaBiotech) and stored at -20°C. Sixteen replicates of each oligonucleotide were printed onto Vapour Phase Coated Glass Slides (Amersham Pharmacia Biotech) using a Molecular Dynamics Gen III spotter at the African Centre for Gene Technologies (ACGT) Microarray Facility, University of Pretoria, Pretoria, South Africa http://fabinet.up.ac.za/microarray. Following printing, the slides were allowed to dry overnight at 45-50% relative humidity. Spotted DNA was then bound to the slides by UV cross-linking at 250 mJ and baked at 80°C for 2 h. The DNA internal transcribed spacer oligonucleotides ITS1, ITS3 and ITS4 served as controls for global normalization and were spotted at concentrations of 50 ng/μl, 100 ng/μl, 150 ng/μl and 200 ng/μl onto the array.

Labeling of target DNA

For target labeling, DNA was extracted from the forty fungi listed in Table 1 using the DNA extraction procedure described before. Extracted DNA was precipitated in 90% ethanol and 0.9 mM NaAc (pH 5.2) to exclude low-molecular-weight fragments. The precipitate was collected by centrifugation at 3,600 g for 30 min. Two micrograms of DNA was labelled with Cy3 or Cy5 by using a Cy™Dye Post-labelling Reactive Dye Pack (GE Healthcare, UK). Each labelling reaction contained DNA diluted in 5 μl 0.2 M Na2CO3 (pH9) and 2.5 μl Cy5 mono NHS ester 4000 pmol dye resuspended in 12 μl DMSO. The reactions were incubated at room temperature for 90 minutes in the dark. After labeling, the dye coupling reaction was column-purified using the QiaQuick PCR purification (QIAGEN GmbH).

Hybridization

Vacuum-dried Cy5-labeled target and 0.3 pmol of the Cy3-labeled control probes were resuspended in 40 μl of hybridisation mixture containing 50% formamide (SIGMA), 25% 2× hybridization buffer (Amersham Pharmacia Biotech), and 25% deionized water. This mixture was denatured at 95°C for five minutes and stored on ice for hybridization. The hybridization solution was pipetted onto a glass slide, covered with a cover slip (24 × 60 mm, No.1, Marienfeld, Germany) and inserted into a custom-made hybridization chamber (N.B. Engineering Works, Pretoria, South Africa). The hybridization was performed overnight at 53°C. After hybridization, the slides were washed twice in 2× SSC and 0.2% SDS at 37°C for 6 minutes, once in 0.2× SSC and 0.2% SDS at room temperature for 5 minutes and twice in 0.075× SSC at room temperature for 5 min. The slides were rinsed in de-ionised water for 2 s and dried by centrifugation at 1000 × g for 5 minutes.

Data acquisition and processing

Oligonucleotide arrays were scanned with a GenePix 4000B scanner (Molecular Dynamics, USA). The mean pixel intensity of each array that resulted from the individual hybridizations was quantified with the Array Vision 6.0 software (Imaging Research Inc., Molecular Dynamics, USA). Individual net signal intensities were obtained by subtracting the local background from the raw spot intensity value. Irregular spots were manually flagged for removal. Further data analysis was performed in the Microsoft Excel software (Microsoft, Richmond, Washington). Anomalous spots not detected through manual inspection were flagged for removal, if the signal intensity of a spot varied more than 10% from the mean of the sixteen replicates on each slide. Signal intensities of the sixteen replicates were then averaged and intensity values were normalized across slides by global regression on the spot intensity data of the internal transcribed spacer oligonucleotides ITS1, ITS3 and ITS4, which were used as a reference for normalization of all spot intensity data (reference design). The net signal intensity of each spot was divided by the median signal intensity of the sixteen replicates and spots with an SNR ((Signal median - Background median) × Standard deviation Background) value below the median were removed from the analysis [32]. Each spot was then either assigned a 1 (present, SNR>/= 3.0) or a 0 (absent, SNR<3.0) according to the median SNR value. The probes with the highest SNR value were considered to be the best target-probe match.

The data discussed above has been deposited at NCBI Gene Expression Omnibus (GEO) [33] and is accessible through GEO series accession number GSE19227.

Reproducibility of the array

The reproducibility of the array was tested using fungal DNA that was independently extracted from eight blind fungal samples obtained from the Forestry and Agricultural Biotechnology Institute, Pretoria. Fungal DNA was labeled and hybridized to the diagnostic chip. For each hybridization experiment, one technical replicate (using independent labeling reactions) was performed, each replication consisting of a reverse labelling experiment. Data analysis was done as described above and binary scores were obtained. Signal intensity values of replicate hybridizations were plotted against each other in Microsoft Excel to verify that the independent fungal samples showed the same scoring pattern. The results were also compared in each case to the identity obtained for the same culture grown by standard laboratory procedures. In addition, the probes positively identified were used for PCR amplification of the eight samples and the results obtained for the array were confirmed with the PCR product amplified from the same sample. The BLAST program was used to obtain the identities of the amplicons. The same procedure was followed for the mycotoxin biosynthesis genes.

Abbreviations

DNA:

deoxyribonucleic acid

EF-1 α:

elongation factor 1-alpha

ITS:

internal transcribed spacer

LNA:

locked nucleic acids

PCR:

polymerase chain reaction

SNP:

single nucleotide polymorphisms

ssu rRNA:

small subunit ribosomal RNA.

References

  1. Barrett JR: Mycotoxins: Of molds and maladies. Environ Health Perspectives. 2000, 108: A20-A27. 10.2307/3454282.

    Article  CAS  Google Scholar 

  2. Mellor S: Problem of Mycotoxins and some solutions. Pig progress. 2003, 5: 12-15.

    Google Scholar 

  3. Rabie CJ, Marais GJ: Toxigenic fungi and mycotoxins in South African foods and feeds. 2000, Report to the Department of Health, Pretoria

    Google Scholar 

  4. Peraica M, Radic B, Lucic A, Pavlovic M: Toxic effects of mycotoxin in humans. Bulletin WHO. 1999, 77: 754-756.

    CAS  Google Scholar 

  5. Niessen L: PCR-based diagnosis and quantification ofmycotoxin producing fungi. Int J Food Microbiol. 2007, 119: 38-46. 10.1016/j.ijfoodmicro.2007.07.023.

    Article  CAS  PubMed  Google Scholar 

  6. Hebart HJ, Loffler J, Meissner C, Serey F, Schmidt D, Bohme A, Martin H, Engel A, Bunje D, Kern WV, Schumacher U, Kanz L, Einsele H: Early detection of Aspergillus infection after allogeneic stem cell transplatation by polymerase chain reaction screening. J Infec Dis. 2000, 181: 1713-1719. 10.1086/315435.

    Article  CAS  Google Scholar 

  7. Mirhendi H, Diba K, Kordbacheh P, Jalalizand N, Makimura K: Identification of pathogenic Aspergillus species by a PCR-restriction enzyme method. J Med Microbiol. 2007, 56: 1568-1570. 10.1099/jmm.0.47319-0.

    Article  CAS  PubMed  Google Scholar 

  8. Mishra PK, Fox RTV, Culham A: Development of a PCR-based assay for rapid and reliable identification of pathogenic Fusaria. FEMS Microbiol Lett. 2003, 218: 329-332. 10.1111/j.1574-6968.2003.tb11537.x.

    Article  CAS  PubMed  Google Scholar 

  9. Waalwijk C, Lee van der T, de Vries I, Hesselink T, Arts J, Kema GHJ: Synteny in toxigenic Fusarium species: the fumonisin gene cluster and the mating type region as examples. Eur J Plant Pathol. 2004, 110: 533-544. 10.1023/B:EJPP.0000032393.72921.5b.

    Article  CAS  Google Scholar 

  10. Paterson RRM, Archer S, Kozakiewicz Z, Lea A, Locke T, O'Grady E: A gene probe for the patulin metabolic pathway with potential for use in patulin and novel disease control. Biocontrol Science and Technol. 2000, 10: 509-512. 10.1080/09583150050115098.

    Article  Google Scholar 

  11. Kristensen R, Gauthier G, Berdal KG, Hamels S, Remacle J, Holst-Jensen A: DNA microarray to detect and identify trichothecene- and moniliformin-producing Fusarium species. J Appl Microbiol. 2007, 102: 1060-1070.

    CAS  PubMed  Google Scholar 

  12. Uttamchandani M, Neo JL, Ong BNZ, Moochhala S: Applications of microarrays in pathogen detection and biodefence. Trends Biotechnol. 2008, 27: 53-61. 10.1016/j.tibtech.2008.09.004.

    Article  PubMed  Google Scholar 

  13. Leinberger DM, Schumacher U, Autenrieth IB, Bachmann TT: Development of a DNA microarray for detection and identification of fungal pathogens involved in invasive mycoses. J Clinical Microbiol. 2005, 43: 4943-4953. 10.1128/JCM.43.10.4943-4953.2005.

    Article  CAS  Google Scholar 

  14. DeSantis TZ, Stone CE, Murray SR, Moberg JP, Andersen GL: Rapid quantification and taxonomic classification of environmental DNA from both prokaryotic and eukaryotic origins using a microarray. FEMS Microbiol Lett. 2005, 245: 271-278. 10.1016/j.femsle.2005.03.016.

    Article  CAS  PubMed  Google Scholar 

  15. Schmidt-Heydt M, Geisen R: A microarray for monitoring the production of mycotoxins in food. Int J Food Microbiol. 2007, 117: 131-140. 10.1016/j.ijfoodmicro.2007.01.014.

    Article  CAS  PubMed  Google Scholar 

  16. Vora GJ, Meador CE, Stenger DA, Andreadis JD: Nucleic acid amplification strategies for DNA-microarray-based pathogen detection. Appl Environ Microbiol. 2004, 70: 3047-3054. 10.1128/AEM.70.5.3047-3054.2004.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  17. Johnson MP, Haupt LM, Griffiths LR: Locked nucleic acids (LNA) singlenucleotide polymorphism (SNP) genotype analysis and validation using real-time PCR. NAR. 2004, 32: e55-10.1093/nar/gnh046.

    Article  PubMed Central  PubMed  Google Scholar 

  18. You Y, Moreira BG, Behlke MA, Owczarzy R: Design of LNA probes that improve mismatch discrimination. Nucleic Acids Res. 2006, 34: e60-10.1093/nar/gkl175.

    Article  PubMed Central  PubMed  Google Scholar 

  19. White TJ, Bruns T, Lee S, Taylor J: Amplification and direct sequencing of fungal ribosomal RNA genes for phylogenetics. PCR Protocols: A Guide to Methods and Applications. Edited by: Innis MA, Gelfand DH, Sninsky JJ, White TJ. 1990, San Diego: Academic Press Inc, 315-322.

    Google Scholar 

  20. Kane MD, Jatkoe TA, Stumpf CR, Lu J, Thomas JD, Madore SJ: Assessment of the sensitivity and specificity of oligonucleotide (50mer) microarrays. Nucleic Acids Res. 2000, 28: 4552-4557. 10.1093/nar/28.22.4552.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  21. Letowski J, Brousseau R, Masson L: Designing better probes: effect of probe size, mismatch position and number on hybridization in DNA oligonucleotide microarrays. J Microbiol Methods. 2004, 57: 269-278. 10.1016/j.mimet.2004.02.002.

    Article  CAS  PubMed  Google Scholar 

  22. Anthony RM, Brown TJ, French GL: Rapid diagnosis of bacteremia by universal amplification of 23S ribosomal DNA followed by hybridization to an oligonucleotide array. J Clinical Microbiol. 2000, 38: 781-788.

    CAS  Google Scholar 

  23. Volokhov D, Rasooly A, Chumakov K, Chizhikov V: Identification of Listeria species by microarray-based assay. J Clinical Microbiol. 2002, 40: 4720-4728. 10.1128/JCM.40.12.4720-4728.2002.

    Article  CAS  Google Scholar 

  24. Graf A, Gasser B, Dragostis M, Sauer M, Leparc GG, Tuechler T, Kreil DP, Mattanovich D: Novel insights into the unfolded protein response using Pichia pastoris specific DNA microarrays. BMC Genomics. 2008, 9: 390-10.1186/1471-2164-9-390.

    Article  PubMed Central  PubMed  Google Scholar 

  25. Lane S, Everman J, Logea F, Call DR: Amplicon structure prevents target hybridization to oligonucleotide microarrays. Biosensors and Bioelec. 2004, 20: 728-725. 10.1016/j.bios.2004.04.014.

    Article  CAS  Google Scholar 

  26. Southern E, Mir K, Shepinov M: Molecular interactions on microarrays. Nature Genet. 1999, 21: 5-9. 10.1038/4429.

    Article  CAS  PubMed  Google Scholar 

  27. Zheng XZ, Kong F, Halliday C, Chen S, Lau A, Playford G, Sorrell TC: Reverse Line Blot Hybridization Assay for Identification of Medically Important Fungi from Culture and Clinical Specimens. J Clinical Microbiol. 2007, 45: 2872-2880. 10.1128/JCM.00687-07.

    Article  Google Scholar 

  28. Hope D, Aldred D, Magan N: Comparison of environmental profiles for growth and deoxynivalenol production by Fusarium culmorum and F. graminearum on wheat grain. Lett Appl Microbiol. 2005, 40: 295-300. 10.1111/j.1472-765X.2005.01674.x.

    Article  CAS  PubMed  Google Scholar 

  29. Raeder U, Broda P: Rapid preparation of DNA from filamentous fungi. Lett Appl Microbiol. 1985, 1: 17-20. 10.1111/j.1472-765X.1985.tb01479.x.

    Article  CAS  Google Scholar 

  30. Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tool. Nucleic Acids Res. 1997, 25: 4876-4882. 10.1093/nar/25.24.4876.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  31. O'Donnell K, Cigelnik E: Two divergent intragenomic rRNA ITS2 types within a monophyletic lineage of the fungus Fusarium are nonorthologous. Molecular Phylogenet Evol. 1997, 7: 103-116. 10.1006/mpev.1996.0376.

    Article  Google Scholar 

  32. Martens M, Weidner S, Link B, de Vos P, Gillis M, Willems A: 2007. A Prototype Taxonomic Microarray Targeting the rpsA Housekeeping Gene Permits Species Identification within the Rhizobial Genus. Syst Appl Microbiol. 2007, 30 (5): 390-400. 10.1016/j.syapm.2007.01.002.

    Article  CAS  PubMed  Google Scholar 

  33. NCBI Gene Expression Omnibus (GEO).http://ncbi.nlm.nih.gov/geo/

  34. Bluhm BH, Flaherty JE, Cousin MA, Woloshuk CP: Multiplex polymerase chain reaction assay for the differential detection of trichothecene- and fumonisin-producing species of Fusarium in cornmeal. J Food Proteomics. 2002, 65: 1955-1961.

    CAS  Google Scholar 

  35. Martin KJ, Rygiewicz PT: Fungal-specific PCR primers developed for analysis of the ITS region of environmental DNA extracts. BMC Microbiology. 2005, 5: 28-10.1186/1471-2180-5-28.

    Article  PubMed Central  PubMed  Google Scholar 

  36. Barros E, Crampton M, Marais G, Lezar S: A DNA-based method to quantify Stenocarpella maydis in maize. Maydica. 2008, 53: 125-129.

    Google Scholar 

  37. Doohan FM, Parry DW, Jenkinson P, Nicholson P: The use of species-specific PCR-based assays to analyse Fusarium ear blight of wheat. Plant Path. 1998, 47: 197-205. 10.1046/j.1365-3059.1998.00218.x.

    Article  CAS  Google Scholar 

  38. Williams KJ, Dennis JI, Smyl C, Wallwork H: The application of species-specific assays based on the polymerase chain reaction to analyse Fusarium crown rot of durum wheat. Australasian Plant Path. 2002, 31: 119-127. 10.1071/AP01079.

    Article  Google Scholar 

  39. Yoder WT, Christianson LM: Species-specific primers resolve members of Fusarium section Fusarium - Taxonomic status of the edible "Quorn" fungus reevaluated. Fungal Genet and Biol. 1998, 23: 68-80. 10.1006/fgbi.1997.1027.

    Article  CAS  Google Scholar 

  40. Möller EM, Chelkowski J, Geiger HH: Species-specific PCR assays for the fungal pathogens Fusarium moniliforme and Fusarium subglutinans and their application to diagnose maize ear rot disease. J Phytopath. 1999, 147: 497-502. 10.1111/j.1439-0434.1999.tb03856.x.

    Article  Google Scholar 

  41. Nicholson P, Simpson DR, Weston G, Rezanoor HN, Lees AK, Parry DW, Joyce D: Detection and quantification of Fusarium culmorum and Fusarium graminearum in cereals using PCR assays. Physiol and Mol Plant Path. 1998, 53: 17-37. 10.1006/pmpp.1998.0170.

    Article  CAS  Google Scholar 

  42. Mulé G, Susca A, Stea G, Moretti A: A species-specific PCR assay based on the calmodulin partial gene for identification of Fusarium verticillioides, F. proliferatum and F. subglutinans. Eur J Plant Path. 2004, 110: 495-502. 10.1023/B:EJPP.0000032389.84048.71.

    Article  Google Scholar 

  43. 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. 10.1016/S0168-1605(02)00250-7.

    Article  CAS  PubMed  Google Scholar 

  44. Dombrink-Kurtzman MA: The sequence of the isoepoxydon dehydrogenase gene of the patulin biosynthetic pathway in Penicillium species. Antonie van Leeuwenhoek. 2007, 91: 179-189. 10.1007/s10482-006-9109-3.

    Article  CAS  PubMed  Google Scholar 

  45. Dombrink-Kurtzman MA: The isoepoxydon dehydrogenase gene of the patulin metabolic pathway differs for Penicillium griseofulvum and Penicillium expansum. Antonie van Leeuwenhoek. 2005, 89: 1-8. 10.1007/s10482-005-9002-5.

    Article  PubMed  Google Scholar 

  46. Lee L, Han Y-K, Kim K-H, Yun S-H, Lee Y-W: Tri13 and Tri7 Determine Deoxynivalenol- and Nivalenol-Producing Chemotypes of Gibberella zeae. Appl and Environ Microbiol. 2002, 68: 2148-2154. 10.1128/AEM.68.5.2148-2154.2002.

    Article  CAS  Google Scholar 

  47. Nicholson P, Simpson DR, Wilson AH, Chandler E, Thomsett M: Detection and differentiation of trichothecene and enniatin-producing Fusarium species on small-grain cereals. Eur J Plant Path. 2004, 110: 503-514. 10.1023/B:EJPP.0000032390.65641.a7.

    Article  CAS  Google Scholar 

  48. Niessen ML, Vogel RF: Group specific PCR-detection of potential trichothecene-producing Fusarium-species in pure cultures and cereal samples. Syst Appl Microbiol. 1998, 21: 618-631.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

This study benefited from the financial support of the Young Researchers Establishment Fund (YREF). Ms Adriaana Jakobs is thanked for assistance with the identification of the fungal strains used in this study.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Sabine Lezar.

Additional information

Authors' contributions

SL: conceived the study, designed the experiment, microarray study, statistical analysis and drafted the manuscript. EB: participated in the study co-ordination and helped to draft the manuscript. Both authors read and approved the final manuscript.

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Lezar, S., Barros, E. Oligonucleotide microarray for the identification of potential mycotoxigenic fungi. BMC Microbiol 10, 87 (2010). https://doi.org/10.1186/1471-2180-10-87

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2180-10-87

Keywords