Development and performance assessment of a luminex xMAP® direct hybridization assay for the detection and identification of indoor air fungal contamination

Considered as a public health problem, indoor fungal contamination is generally monitored using classical protocols based on culturing. However, this culture dependency could influence the representativeness of the fungal population detected in an analyzed sample as this includes the dead and uncultivable fraction. Moreover, culture-based protocols are often time-consuming. In this context, molecular tools are a powerful alternative, especially those allowing multiplexing. In this study a Luminex xMAP® assay was developed for the simultaneous detection of 10 fungal species which are most frequently in indoor air and that may cause health problems. This xMAP® assay was found to be sensitive, i.e. its limit of detection is ranging between 0.05 and 0.01 ng of gDNA. The assay was subsequently tested with environmental air samples which were also analyzed with a classical protocol. All the species identified with the classical method were also detected with the xMAP® assay, however in a shorter time frame. These results demonstrate that the Luminex xMAP® fungal assay developed in this study could contribute to the improvement of public health and specifically to the indoor fungal contamination treatment.


Introduction
Currently, indoor airborne fungal contamination is suggested to be associated with public health problems [1]. Even if indoor air fungal contaminants could be allergenic or could have an implication in respiratory diseases, such as asthma, wheezing or rhinitis [2][3][4][5][6][7] the scientific evidence for the causal link between these molds and adverse health effects is still poorly documented. The use of classical methods involving culturing and microscopic visualization in routine monitoring analysis could be pointed out as one of the reasons explaining this lack of evidence in scientific literature. Even if culture-dependent tools are useful and well documented, these techniques are a1111111111 a1111111111 a1111111111 a1111111111 a1111111111 Moreover, the use of Luminex xMAP 1 technology has the advantage that the sample size requirements are reduced and that no skilled mycologist is needed to perform the microscopicbased identification analysis.

Fungal strains and DNA isolation
All the fungal species and strains used in this study are listed in Table 1. All of them were purchased from the BCCM/IHEM collection located at the Scientific Institute of Public Health in Brussels (WIV-ISP, Belgium).
The culture and extraction protocol were both previously described in Libert et al. (2015) [10]. Briefly, all the strains were incubated with constant agitation in a S10 Sabouraud liquid medium (Biorad, Temse, Belgium) at 25˚C during 3 to 10 days depending on the species' culture conditions. Then, after a centrifugation of 1 min at 12 000g to eliminate all of the Sabouraud liquid, 0.25 ml of acid washed glass beads (Sigma Aldrich, Diegem, Belgium) were added to the wet sample (300 mg) transferred into cryotubes and put at -80˚C during 40 minutes. After that, the samples were freeze-dried overnight with a freeze-dryer Epsilon 1-6D (Martin Christ, Osterode am Harz, Germany) and subsequently bead-beaten (1 minute, maximal speed) with a Mini bead beater (Biospec Products, OK, USA).
Finally, an adapted phenol chloroform (24:1) protocol (Ashktorab and Cohen 1992) was applied to extract DNA, which was then purified with the Qiagen CTAB genomic Tip-20 kit (Qiagen Benelux-B.V., KJ Venlo, the Netherlands) and eluted with 100μl Gibco 1 DNase, RNase, protease free water (Life Technologies, Gent, Belgium). The purity and the amount of extracted DNA were evaluated with a Nanodrop 1 2000 (Thermo Scientific, Wilmington, USA).DNA integrity was verified on a 2% agarose gel.
The PCR amplifications were carried out with following thermal cycler programme: 30 s at 98˚C (initial denaturation); 35 cycles of 30 s at 98˚C (denaturation), of 1 min at 55˚C (annealing) and 1 min at 72˚C (extension); a final extension at 72˚C during 10 min and a final hold at 16˚C.

Probe selection
All probes were selected from literature (qPCR Taqman probes transferred into Luminex probes, Table 1) and adapted if needed in order to have a similar range of length. Probe quality estimation (hairpins, ΔG values. . .) was performed with Visual OMP version 7.8.42.0 (DNA software, Washington, USA). An in silico analysis of each probe was done with the Thermoblast tool from Visual OMP version 7.8.42.0 (DNA software, Washington, USA).
All probes were manufactured by Eurogentec (Liège, Belgium) and tagged with a 5'-end amino modifier C12, followed by a RP-Cartridge-Gold purification.

Probe coupling to xMAP ® beads
According to the Luminex xMAP 1 technology, each specific probe is coupled to a specific set of beads. The coupling protocol used in this study is based on the Luminex recommendations for carbodiimide coupling of amine-modified oligonucleotides to MagPlex magnetic carboxylated microspheres (beads) [25]. The final concentration of the coupled beads was 12500 beads/ μl. These working stocks of coupled beads were stored in the dark at 4˚C until their use.
The hybridization mix contained per reaction 33 μl of a specific coupled bead set (76 beads/ μl), 12 μl of Tris-EDTA buffer, pH 8 (Sigma-Aldrich, Diegem, Belgium) and 5 μl of fresh PCR amplicons. The hybridization reaction was performed on an Aeris thermal cycler (Esco, Barnsley, The Netherlands) according to the following protocol i.e., a first step at 96˚C during 1 min 30 s and a second one at 58˚C during 30 min. Before a third incubation step of 5 min at 58˚C, 25 μl of reporter mix composed of 4 μg/ml of SAPE (Streptavidin-R-Phycoerythrin) (Life Technologies, Gent, Belgium) and 1X TMAC buffer (Sigma-Aldrich, Diegem, Belgium), were added to each sample.
Finally, all analyses were performed on a MAGPIX device (Luminex Corporation, Austin, USA) equipped with the xPONENT for MAGPIX v4.2 software (Luminex Corporation, Austin, USA). The protocol applied in all runs was performed at 58˚C with a minimum of 50 beads counted for each bead region. A wash of each sample was also carried out on the machine during the analysis.
Because each bead set has a unique spectral address (distinct red color code) and each PCR amplicon hybridized to the probe bound to the beads is labeled with SAPE, the fluorescence intensity (red and green) gives information on the amount of beads per region (bead set) and on the amount of beads bound to a PCR amplicon. This last information is given by the median fluorescence intensity value (MFI) and is defined for each region (bead set).

Data analysis and interpretation
The data analysis and the interpretation of the results were based on Wuyts et al. (2015) [26]. At the end of each run, the bead counts were checked to verify whether the bead count was homogenous for all of the coupled bead sets. If this was not the case, the run was repeated. Then, the median fluorescence intensity (MFI) of each target was used to calculate a signal-tonoise ratio (SN) with the following formula i.e., where SN target a is the signal-to-noise ratio observed for the set of coupled beads selected for the detection of the target a; MFI sample target a corresponds to the MFI value (collected by the xPONENT software) observed for the target in a specific sample; and MFI NTC target a is the MFI value obtained for the NTC target a According to Wuyts et al. (2015) [26], a result was considered as positive if SN target a ! 3.00. For each SN ratio close to the limit (i.e. close to 3), a t-test (95% confidence) was performed with the SN ratio obtained for the negative control. If the difference is significant ( Ã ), the results were considered as positive. If no difference was observed with the negative control, the data were considered as negative.

Specificity test
The specificity of each coupled beads-probes was tested with 3 different assays: a simplex analysis where each probe was tested only on its targeted species (i.e., one bead set, one species) and a multiplex analysis where each species was subjected to each coupled bead set at the same time. Finally, a mix of gDNA extracted from different species was tested.
Simplex analysis (DNA of 1 species, 1 set of beads). The probe specificity was firstly tested during a simplex analysis. This test consisted of the one by one analysis of all the targeted species with their specific coupled beads-probe set. So, 10 specific hybridization mixes were made, i.e. one mix per set of beads and one mix per species, and analyzed in duplicate during 3 independent runs. All analyses were performed with 5 μl of PCR amplicon, obtained as described above. The repetitions were done with DNA template extracted from independent cultures. For each mix and run, one non-template-control PCR reaction (NTC PCR ) i.e., gDNA replaced by water, was introduced in order to evaluate the background linked to the analysis.
Multiplex analysis (DNA of 1 species, multiple sets of beads). The second step of the specificity evaluation consisted of the multiplex analysis. In this test, the hybridization mix contained every of the coupled bead sets, and this mix was tested on each species in duplicate in 3 independent runs. Every analysis was performed with 5 μl of PCR amplicon, obtained for each species with 5 ng of gDNA from pure culture, as elaborated above. The repetitions were done with DNA template extracted from independent cultures. A NTC PCR was added to all runs.
DNA mix analysis (DNA of multiple species, multiple sets of beads). Subsequently, a Luminex analysis was performed on the product of a PCR reaction containing a mix of gDNA extracted from different species. Nine different PCR mixes were made ( Table 2). Mix 1 contained 5 ng of gDNA from all targeted species as mentioned in the Table 2 and was considered as a positive control.
For the mixes 2 to 9, gDNA of all species targeted in the test was added in the PCR mix, except for one, whereby the missing species changed for each mix ( Table 2). All analyses of mixes were done in duplicate in 3 independent runs. For every analysis, a NTC PCR was added.

Sensitivity evaluation: Limit of detection
To evaluate the sensitivity of the Luminex assay developed in this study, a serial dilution of gDNA of each targeted species was made to determine the limit of detection (LOD). Because no guidelines exist on the development and the performance assessment of molecular methods for fungal detection, the LOD estimation performed in this study is based on the workflows elaborated for the validation of molecular methods for the detection of GMO and food-pathogens. In these fields, the LOD is defined as the lowest concentration of an analyte which is detected with a probability of 95% [27,28].
In order to estimate this LOD, 9 dilutions from 10 to 0.001 ng of gDNA of each species were amplified independently in a PCR reaction according to the above described PCR protocol. All of the PCR reactions were analyzed in duplicate in 3 independent runs.

Proof of concept with real-life environmental samples
The environmental testing was performed on real-life air samples previously collected from contaminated houses [10]. The protocols for the sampling, the DNA extraction and microscopic determination were previously described in Libert et al. (2015) [10]. All the samples were collected in duplicate. Briefly, air samples were collected in contaminated habitats using the Coriolis μ air sampler (Bertin Techologies, Montigny-Le-Bretonneux, France) which collects contaminants into a liquid (15 ml) compatible with the classical culture-based approach as well as with molecular-based tools. A first group of samples was centrifuged during 15 min A.alternata 3320 at 5000 g in order to extract the DNA using an adapted phenol-chloroform DNA extraction protocol (including bead-beating and freeze-drying steps, followed by a CTAB Tip20 purification protocol [10]). In parallel, the second group of samples was analyzed according to the culture-based protocol [10,12] which included an incubation step of 5 days at 25˚C for mesophilic and 2 days at 45˚C for thermophilic fungi and a microscopic visualization for species determination.
The Luminex assay was applied on samples that were all collected in the same house in 4 different rooms (i.e. bathroom, bedroom, kitchen and living-room). The analysis was performed in duplicate with 3 independent runs according to the protocols described above. The PCR amplifications were done with 5 ng of gDNA extracted from each sample. In each run a NTC PCR and 8 positive controls (i.e., gDNA extracted from pure cultures of A. alternata, A. versicolor, A. fumigatus, C. cladosporioides, C. herbarum, P. chrysogenum, S. chartarum and U. botrytis) were added. It should be noted that neither A. creber nor A. sydowii positive controls were added, because their detection was done with the same probe than that for the detection of A. versicolor.
For each SN ratio close to the limit (i.e. close to 3), a t-test (confidence 95%) was performed with the higher negative results considered as a negative control. If the difference is significant ( Ã ), the results were considered as positive. If the no difference was observed with the negative control, the data were considered as negative.
To verify that no inhibitors were present in the environmental samples extracts and to confirm that if the species would be present, it could be detected with the Luminex fungal assay, a final amount of 5 ng of gDNA from pure culture of each species not-detected during the Luminex assay and the classical monitoring were spiked into the DNA extracted from the 4 environmental samples. The spiked DNA extract comes from the same culture than that was used to prepare the DNA used as positive controls. The Luminex analysis was subsequently repeated.

Probe selection
This Luminex assay aims at the detection of the 10 fungal species most frequently found in indoor air in Belgium and in Europe and that may cause health problems i.e., A. alternata, A. creber, A. fumigatus, A. sydowii, A. versicolor, C. cladosporioides, C. herbarum, P. chrysogenum, S. chartarum and U. botrytis [2,22,23]. The Luminex xMAP 1 technology used in this study is based on the direct DNA hybridization to a specific probe coupled with a unique set of beads. Specific probes to be coupled to the beads were selected in the internal transcribed spacer (ITS) regions of the ribosomal DNA. These ribosomal regions have the advantages to be conserved and to show, at the genus level, few polymorphisms [29,30,31], allowing the specific detection of particular species. Another advantage is the fact that the small subunit (SSU), 5.8 S and the large subunit (LSU) from the rRNA genes flanking the ITS 1 and 2 regions, are sufficiently conserved among species to design some universal primers such as the primers ITS1 and ITS2 and ITS3 and ITS4 [24] to create the PCR amplicons to be hybridized to the probes. The probes were designed based on qPCR Taqman probes available in literature. Six probes i.e., AaltP2.2, UloP1, StachP2, AfumP1, CherbP1 and VersP1, detect amplicons in ITS 1, while Pchris1 and CcladP1 detect amplicons the ITS 2 region ( Table 1). All of them are specific to their target except the probe VersP1 from the EPA [15,17] which can also hybridize to the ITS 1 region of A. creber and A. sydowii. The probes were evaluated in silico and adapted if needed in order to have a similar range of length (Table 1). Indeed, the stabilization of the formation of the hybridization complex between probe and PCR amplicon is assured by the addition of TMAC which reinforces AT base-pairs [19]. Consequently, the hybridization efficiency is more influenced by the length of the probe than by the nucleotide composition [19]. According to the in silico testing, a consensus length for an optimal detection of each target was found to be between 29 and 31 nucleotides (Table 1).

Specificity test
The specificity of each probe was tested in three steps i.e., a simplex analysis, a multiplex analysis and the multiplex analysis of a mix of DNA. With the simplex analysis, it was verified whether the protocol and probe can detect the amplicon from the targeted species only (i.e., one species, one probe, one bead set for detection of one species); secondly, with the multiplex analysis, it was investigated if no aspecific annealing occurred when a specific amplicon was put into a mix of probes, the specific probe included (i.e., one species, multiple sets of beads mixed for detection of one species); thirdly, with the multiplex analysis of a mix of DNA, it was verified whether no incorrect detection happened when a mix of amplicons from different species was analyzed with a mix of different beads (i.e., multiple species mixed, multiple sets of beads mixed for multiple specific detection).
Simplex and multiplex analysis. During these tests, the PCR amplicon of each species was detected by its own specific coupled probe-bead set in every of the 6 repetitions and this both in the simplex analysis (Table 3) as well as in the multiplex analysis ( Table 4). The simplex analysis yielded average SN ratios ranging between 3.50 ±0.25 (for A. sydowii) and 22.69 ±1.54 (for U. botrytis) ( Table 3). The average SN ratios obtained in the multiplex analysis ranged between 3.52 ±0.05 and 27.15 ±0.18 for C. cladosporioides and U. botrytis, respectively (Table 4). During these analyses, no false positives as well as no false negatives were observed (Tables 3 and 4), following the criteria defined for obtaining a positive result (i.e., MFI ratio !3.00).
DNA mix analysis. Subsequently, an analysis was performed on 9 different DNA mixes (Table 2), with a design allowing to verify that the detection of each species is still possible in the presence of other species and with a mix of set of beads. All the positive controls were correctly detected for each repetition (6/6) with a SN ratio ranging between 4.31 ±0.92 for A. fumigatus and 17.82 ±0.63 for A. versicolor. The first mix contained all of the targeted species. The Luminex analysis resulted in an SN ratio !3.00 (i.e. all positive) for all of the expected species and probes for each repetition (6/6 positive for all), with a lowest SN ratio observed at 4.69 ±0.68 (for A. fumigatus) and the highest SN ratio observed at 16.98 ±2.85 (for A. versicolor) ( Table 5). In the other mixes (mixes 2 to 9) each time one species was omitted (Tables 2 and 5). The gDNA of each target species introduced in the mixes was each time detected (6/6) with SN ratios ranging from 3.37 ±0.25 for the detection of S. chartarum in the mix 7 to 26.10 ±1.37 for P. chrysogenum in the mix 4 ( Table 5), indicating that no false negatives were obtained. No positive Luminex signal was obtained for any coupled bead sets for which no corresponding specific gDNA was added to the mix (SN ratios ranged between 0.92 ±0.20 and 2.62 ±0.25) ( Table 5). This means that for the 6 repetitions, no false positives were observed during this test.

Sensitivity: Limit of detection
The sensitivity of this Luminex assay was determined with 9 points of serial dilutions (from 10 to 0.001 ng of gDNA for each of the targeted species). The LOD was 0.05 ng for A. alternata, A. creber, A. sydowii, A. fumigatus, C. herbarum, P. chrysogenum and S. chartarum and 0.01 ng for A. versicolor, C. cladosporioides and U. botrytis ( Table 6).

Proof of concept using environmental samples
Following the performance assessment, DNA extracted from environmental indoor air samples collected from different rooms inside a contaminated house was analyzed using the fungal Luminex assay. This gave the opportunity to test the molecular tool with real-life samples containing a mix of a priori unknown fungal species at unknown concentration. Therefore, in this study, the proof of concept using environmental samples allows not only testing the performance (detection of DNA at very low or variable concentrations, detection in a mix of species, etc.) of the developed tool but it allows also demonstrating that the developed tool can be used with uncharacterized strains (i.e., strains not coming from a culture collection) of targeted species which could be present in an environmental air sample, and this using DNA extracted from the air samples without prior cultivation step.
The 4 samples used in this study were previously analyzed with classical methods [10] and only 3 species were determined i.e., A. versicolor, Cladosporium sp. and P. chrysogenum. P. chrysogenum was observed in each sample, while A. versicolor was retrieved in 2 of them (i.e., samples from bathroom and living room) and C. cladosporioides in sample 4 (i.e., sample collected in the bathroom) only. Infertile mycelia were also observed in the living room and the kitchen (Table 7).
These 3 species were also detected in the same rooms by the fungal Luminex assay performed on the extracted DNA from these 4 real-life samples (  respectively in the bathroom and the living room and ranged between 3.77 ± 0.09 in the kitchen and 14.83 ±0.22 for P. chrysogenum in the bathroom (Table 8). No other species were detected (Table 8). Because the probe VersP1 is not specific to A. versicolor (Tables 1, 2, 3 and 4), it is important to note that the taxa detected in the sample could also be A. creber and A. sydowii. However, according to the microscopic visualization and based on the fact that A. versicolor is more frequently observed in indoor air from Belgium than A. creber and A. sydowii, the detected species was considered as A. versicolor (Table 8).
In order to verify that no inhibition occurred during the analysis of the environmental samples, and to verify that if a species would have been present, it could be detected in the 4 environmental samples, a spike of all species not detected by the Luminex assay was performed into every environmental sample, according to the results obtained for the first part of the proof of concept analysis. Therefore, based on the results obtained for the analysis of the 4 environmental samples, A. versicolor was spiked into the samples 1 and 2, but not into samples 3 and 4. C. cladosporioides was spiked into the samples 1 to 3 and not in sample 4. Finally, as it was detected in each sample, no spike of P. chrysogenum was performed. Every gDNA spiked into the 4 samples was properly detected (Table 8). Indeed, in the 4 samples each species, spiked or not (as already present), was detected (SN ratio !3.00).

Discussion
For years, fungal indoor air contamination is considered as a public health problem, even if today no substantiated scientific evidence on the causal link exists. This lack of evidence is principally due to a scarcity of data on the full composition of indoor airborne fungal community. Indeed, most of the protocols used today for fungal contamination monitoring are culture-dependent. This implies some limitations, such as that they are not able to reflect the complete indoor fungal community representativeness, as the dead and uncultivable fraction is not detected with these methods. Moreover, they are time-consuming. To improve the collection of data in terms of time and completeness, molecular methods, such as qPCR, have been developed for the detection and identification of the indoor fungal community. However, even if the efficiency of most of these qPCR methods has been well established, their multiplex capacities are still too limited. With its ability to simultaneously analyze up to 50 different targets, the use of the Luminex xMAP 1 technology on a MagPix instrument could significantly upgrade the indoor fungal contamination monitoring as was previously demonstrated for the diagnosis of some relevant fungal and other pathogens from clinical samples [20,21,26,[32][33][34]. This study presents the first Luminex xMAP 1 assay developed for the monitoring of 10 fungal species most frequently found in indoor air in Europe and that may cause health problems, i.e., A. alternata, A. creber, A. fumigatus, A. sydowii, A. versicolor, C. cladosporioides, C. herbarum, P. chrysogenum, S. chartarum and U. botrytis [2,22,23]. While A. creber, A. versicolor, A. sydowii and P. chrysogenum are typical indoor species, A. alternata, A. fumigatus and C. herbarum are known to be outdoor species. Despite their outdoor sources, these 3 species are commonly detected in indoor air samples. Especially observed in the forms of spores, they arrive in indoor environment through the draft and ventilation system (i.e., windows, ventilation or air-conditioning system). Some of these species such as A. alternata, A. versicolor, P. chrysogenum, U. botrytis or S. chartarum are known to have implications in the worsening of respiratory diseases or allergies [5,[35][36][37][38][39]. The detection of these multiple species is therefore important in the context of health issues. The Luminex xMAP 1 technology gives the opportunity to perform in a single run a multiplex analysis able to detect multiple species with high specificity. This specificity is itself defined by the specificity of the probes bound to the each set of beads. These beads are the key elements of this technology.
To have a species-specific detection, the probes were designed to target the ITS1 or ITS2 regions from the ribosomal DNA, which are recognized as the most suitable region for the detection of fungi due to their low intra-species variablility [40,41]. There exist universal primer pairs (i.e. ITS1/2 and ITS 3/4 [24]) to amplify these regions in all fungi. Therefore, targeting of ITS1 and ITS2 also allowed to easily design a multiplex detection based on species-specific probes. Indeed, the targeting of the ITS regions and the use of one set of universal primers for each region allow to reduce the number of step. Consequently, the time and amount of sample needed for the generation of a species-specific amplicon were also decreased. Actually, the PCR amplifications could be done in a duplex PCR reaction with universal primers for fungal gDNA amplification (primers ITS1, 2, 3 and 4). Due to their universal character, the use of these primers also gives the possibility to extend the panel of species to be detected in the future, by adding specific probe-coupled bead sets. This will even more improve the powerful fungal monitoring based on this xMAP 1 assay. Indeed, the specific detection of the generated amplicons is based on the hybridization of these amplicons to specific probes coupled to Luminex bead-sets. Every probe used in this study originates from a previously developed TaqMan 1 qPCR method (Table 1) [17,20]. Therefore, the specificity of each probe used in this study was previously validated [17,20]. According to the Luminex' recommandations, the optimal probe length used for an xMAP 1 analysis is defined between 15 and 20 nucleotides [19]. However, all probes used in this study, obtained from literature, are larger than 20 nucleotides. For most of them, according to the in silico analysis, their specificity was negatively impacted by a size reduction. Therefore, in this study the length defined for the best detection of each species targeted in our mulitplex assays was found to be around 30 nucleotides. To optimize the TMAC-based hybridization step, which is probe-size dependent, all the probes were adapted to have a length between 28 and 31 nuclotides.
The current assay showed a specific detection of each targeted species, as observed during the simplex and multiplex assays. The species-specific detection was further demonstrated with the DNA mixes test. No false negative nor false positives were observed. The remark has to be made that the developed fungal Luminex assay targets 10 different fungal species. Seven of them are detected uniquely with specific probes. The 3 remaining closely related species (i.e., A. creber, A. sydowii and A. versicolor) are detected by one and the same probe, VersP1. This implicates that although the 3 species can be detected using the Luminex assay, no discrimination can be made between them. The poor discrimination between these 3 Versicolores (A. creber, A. sydowii and A. versicolor) species was already observed in other studies [10,15]. To improve the specific identification of these 3 species, an additional marker could be added such as the gene coding for β-tubulin or mycotoxin genes [20,42]. However, the use of a new marker requires the optimization of the PCR workflow (multiplex optimization or addition of an amplification step) and probably of the hybridization temperature in the xMAP 1 workflow, which depends on the probe length. Indeed, with the use of TMAC, the specificity of the annealing depends only on the hybridization temperature, and not on the composition of the probe. Therefore, to avoid any problem of specificity, adding markers (and thus probes) to a workflow will require the optimization of this parameter. When using multiple markers to target closely related species such as the 3 Versicolores in this study, a decision-tree-based workflow could be applied after the Luminex assay to perform the discrimination based on the combination of the signals obtained for each marker. With such a decision-tree-based workflow, the xMAP 1 assay developed in this study could already be used as a first screening without modifications. If a positive result is observed for the VersP1 probe, a second analysis could be performed with some new markers in order to identify the species. This kind of decision tree already exists in other fields e.g., for the detection of genetically modified plants [43]. Alternatively, another molecuar assay could be used, e.g. based on high resolution melting qPCR for the discrimination of closely related species [44]. Being able to discriminate between these 3 species will contribute to our understanding of the impact of indoor fungal contamination and health problems, as currently, the presence of these 3 Aspergillus species in indoor air and their difference in impact on health could not yet be evaluated.
In addition to its specificity, the xMAP 1 assay developed in this study is also sensitive with a LOD ranging between 0.05 and 0.01 ng of gDNA. According to the limited information available on the genome size of each targeted species, these LOD correspond to a range from 2.51 theoretical genomic copy numbers for C. cladosporioides and U. botrytis, to 19.84 theoretical genomic copy numbers for A. fumigatus (Table 9). It should be noted that intraspecies ITS copy number variations observed in some fungal species such as A. alternata, A. fumigatus, A. versicolor or C. cladosporioides could affect the LOD [10,13,45,46], thereby impeding the theoretical genomic copy number variation estimation. Therefore, it has been decided to determine the LOD in mass, instead of in genomic copy numbers.
The proof of concept using environmental samples showed that the xMAP 1 technology can be used for the fungal monitoring in indoor environment. Indeed, the 3 species detected with the classical methods (i.e., A. versicolor, C. cladosporioides and P. chrysogenum) were also detected with the xMAP 1 technology and this without prior cultivation of the sample. Also, the species not detected with the classical protocol were not detected with the xMAP 1 assay. This test indicates that our Luminex assay can be used on mixed and unbalanced concentrations of fungal species without giving the problem of false detection. This observation is supported by the spike test results which showed that the species present in the gDNA extracted from the environmental samples (not spiked) could still be detected even when a high amount of spiked DNA of other species was present. The spike test results also demonstrated that no inhibitors were present in the DNA mix, as the spiked DNA could be detected by the Luminex assay performed on the environmental samples. If the species would have been present in the environmental sample, it would have been detected by the Luminex assay. So, if the species was not detected, it was not present in the environmental samples, or it was present at a concentration below the LOD. According to the observations based on culturing, besides P. chrysogenum, also some infertile mycelia were observed with classical analysis in samples from the kitchen and the bedroom while the xMAP 1 assay only detected P. chrysogenum (Tables 7 and  8). According to these results and those from the spike test, these undetermined taxa do not belong to the targeted species, except for P. chrysogenum. As no DNA sequencing was performed on the infertile mycelia observed on plate, we cannot exclude that these infertile mycelia belong to an untargeted species of our xMAP 1 assay. However, if needed, probes for additional taxa, once determined, can be easily added to the assay, as elaborated above. Additionally, a 'general ITS' probe, detecting all fungal species, could be added to the assay. If none of the specific probes give a Luminex signal, but there is fungal DNA present in the sample (as detected by the general probe), other methodologies, including mass parallel sequencing, could be applied to further characterize the sample, if needed. Nevertheless, this environmental test delivered the proof of concept for the use of the xMAP 1 technology, which is culture-independent and less time-consuming for the analysis of real-life samples i.e., 3 days for the xMAP 1 technology (sampling, DNA extraction included) compared to 5 to 21 days for the classical analysis (sampling and the culturing included).
Our fungal assay has been developed as a detection tool. Because the Luminex xMAP 1 method contains at least one step of PCR amplification, the tool developed in this study should not be considered as a quantitative one. Indeed, due to the exponential amplification of the DNA, the quantification cannot be done accurately, i.e. the detected fluorescence does not translate directly to the number of DNA molecules in the sample. It can merely be considered as a semi-quantitative tool. Additionally, there exists a species-and even strain-dependent copy number variation of the ITS marker, which is still poorly documented as previously discussed [10]. This also complicated the quantification based on this marker. To develop quantitative tools, more studies would have to be made to find a marker that is very well conserved all along the fungal kingdom, and with a constant copy number, preferably a single copy gene, and which could be used for the quantification. Hereto, there is a need for more whole genome sequence data of fungal species (and this for multiple strains per species) to screen the fungal genome in order to find a conserved single copy marker or a unique single copy marker per species, as was previously found for A. fumigatus [45]. Another solution could also be found in the use of new technologies such as digital PCR, which allows absolute quantification, as was done for Candida albicans in blood [50]. However, as digital PCR requires an a priori selection of the fungal species, an absolute quantification of each contaminant cannot be made without a prior screening of the diversity. This screening step, performed before the quantification step, could be done using the Luminex xMAP 1 fungal assay or a metagenomics NGS analysis.
To conclude, this study reported on a fast, specific and sensitive Luminex xMAP 1 assay targeting 10 important fungal contaminants frequently observed in indoor air and that could have health impacts. The use of the xMAP 1 technology allows a culture-independent analysis with a reduced turn-around-time compared to the classical protocols. The use of this molecular multiplex tool to investigate the indoor contamination could improve the monitoring of fungal diversity. The improvement of data on the fungal population in buildings will contribute to the knowledge concerning their impact on public health, especially on respiratory diseases.

Ethical approval
This article does not contain any studies with human participants or animals performed by any of the authors.

Author Contributions
Conceptualization: XL AP FB NR SDK.