High-throughput screening of filamentous fungi using nanoliter-range droplet-based microfluidics

Filamentous fungi are an extremely important source of industrial enzymes because of their capacity to secrete large quantities of proteins. Currently, functional screening of fungi is associated with low throughput and high costs, which severely limits the discovery of novel enzymatic activities and better production strains. Here, we describe a nanoliter-range droplet-based microfluidic system specially adapted for the high-throughput sceening (HTS) of large filamentous fungi libraries for secreted enzyme activities. The platform allowed (i) compartmentalization of single spores in ~10 nl droplets, (ii) germination and mycelium growth and (iii) high-throughput sorting of fungi based on enzymatic activity. A 104 clone UV-mutated library of Aspergillus niger was screened based on α-amylase activity in just 90 minutes. Active clones were enriched 196-fold after a single round of microfluidic HTS. The platform is a powerful tool for the development of new production strains with low cost, space and time footprint and should bring enormous benefit for improving the viability of biotechnological processes.

been used with filamentous fungi. Indeed, growing and screening filamentous fungi in droplets is challenging due to the apical growth of the hyphae after spore germination to form an expanding branched mycelia network 20 . As we show here, filamentous fungi cannot be screened in picoliter droplets in microfluidic systems: the hyphal tips, where most protein secretion occurs 21 , rapidly exit the droplets, leading to uncontrolled coalescence. This limits the incubation period to only a few hours, which is too short to allow secretion of enough enzyme for screening.
This work proposes a novel approach for the HTS of filamentous fungi using nanoliter-range droplet-based microfluidics tools. Single spores can be encapsulated in ~10 nl droplets (volume ~3 orders of magnitude larger than those typically used for sorting bacteria and yeast) which can be incubated and sorted based on fluorescence. This system allowed 24 h incubation of filamentous fungi (Aspergillus niger) and screening of ~7,000 fungi. h −1 based on the activity of a secreted enzyme (α-amylase). A whole-genome UV-mutated A. niger library was enriched 196-fold for fungi secreting α-amylase after one round of microfluidic HTS. Sorted clones were further analyzed using a robotic microplate based system at ~400 clone.h −1 : 98.1% exhibited α-amylase activity. The HTS platform can easily be adapted to screen for other enzymatic activities.

Results
Nanoliter-range droplet-based microfluidic modules. Using nanoliter size droplets is essential for screening of filamentous fungi. Typically, in flasks, the time from which Aspergillus spores germinate to the stage where they start to display detectable enzymatic activity is ~24 h and optimum levels of activity are usually only reached after several days of incubation, depending on the secreted enzyme. The requirement for long incubation time combined with the rapid growth of the fungal hyphae makes screening filamentous fungi in picoliter droplets a challenging task: even in 250 pl droplets the hyphal tips exit the droplets in ~15 h, causing uncontrolled coalescence: (Fig. 1a, Supplementary Movie S1). In contrast, encapsulating single spores in 18 nl droplets allows growth of the branched mycelial network for up to 24 h confined in the droplet (Fig. 1b) with the hyphal tips exiting the droplets only after incubation for 32 h (Fig. 1c). Furthermore, after 24 h incubation, secreted α-amylase was easily detectable in the droplets using a fluorogenic substrate (Fig. 1d).
Increasing droplet volume from the pl to the nl range affects droplet formation, stability and dielectrophoretic sorting. It was straightforward to produce highly monodisperse (polydispersity ≤ 5% by volume) 10-20 nl droplets at 80-90 droplets.s −1 using a flow-focusing 22 device with a constriction orifice 225 μm wide and 250 μm deep ( Supplementary Fig. S1 and Movie S2), an oil phase of Novec HFE-7500 fluorinated oil containing 2.5% (w/w) KryJeffD 900 (a triblock copolymer fluorosurfactant) 23 and an aqueous phase of synthetic growth medium. Droplets were found to be highly stable and could be collected, incubated and reloaded. Analysis of the size distribution of 10 nl droplets after 24 h incubation at 30 °C in the form of a concentrated, creamed emulsion in the glass capillary and re-injection into the sorting chip ( Supplementary Fig. S1) indicated that ≤ 1.0% of droplets were coalesced ( Supplementary Fig. S2).
We developed a fluorescence-activated droplet sorting (FADS) device specifically adapted to sort nanoliter droplets. In dielectrophoretic droplet sorters 12,15,24 , droplets flow in carrier oil towards a Y-shaped junction. With no electric field, all drops flow into the waste channel which offers lower hydrodynamic resistance than the second, collect channel. To direct droplets into the collect channel, on chip electrodes are energized, creating an electrical field gradient, which generates a dielectrophoretic force (DEP) acting on the droplets. In order to be sorted, DEP forces must displace the droplet by a critical distance orthogonal to the flow, d o , in the time, t p , it takes the droplet to traverse the electrical field gradient. The DEP force on a spherical droplet is given by equation (1)  in which, |F DEP | is the magnitude of the DEP force, R is the radius of the droplet being subjected to the DEP force, Re [f CM ] is the Claussius-Mosscoti factor (which itself depends on ε  ⁎ p and ε  ⁎ m : the complex permittivity of the droplet and the carrier fluid, respectively) and ∇ E 2 denotes the electrical field gradient. Assuming that d 0 = 2R, and that, in order to prevent false positives, no more than one droplet can be in the electric field in the sorting chamber at any one time, the theoretical maximum sorting frequency, f sort , can be calculated (see Supplementary  Information): (2) sort m CM 2 where η designates the viscosity of the carrier oil. Equation (2) shows that f sort is directly proportional to R and to ∇ E 2 , and that increasing R apparently favors DEP sorting as lower ∇ E 2 is required. However, larger droplets are more easily split by electric fields 26 . The electric field, E max , above which droplets split, when Maxwell stress surmounts the resistance to deformation due to interfacial tension, is given by the expression: Therefore, the maximum sorting frequency is a function of both the magnitude, E, and inhomogeneity, ∇ E 2 , of the electric field in the sorter and depends on the volume of the drops and the size and geometry of the sorter.
We investigated two different electrode configurations: the first design had both charged and ground electrodes on the same side of the sorter (Fig. 2a), as used for fluorescence-activated sorting of pl droplets 12,15 and the second design had the charged and ground electrodes on the opposite sides of the sorter (Fig. 2b). Finite element analysis (FEA) simulations indicated that in the same-side electrode configuration, although the DEP force was high between the positive and ground electrodes close to the channel wall adjacent to the electrodes, it drops off rapidly across the width of the channel ( Fig. 2a and Supplementary Fig. S3). In contrast, using the cross-side electrode configuration the DEP force extends across the entire width of the sorting channel ( Fig. 2b and Supplementary Fig. S3). By combining a model of the electric field distribution with equations (2) and (3) it was possible to calculate the maximum droplet sorting frequency limited by electrosplitting, f e , as a function of droplet volume ( Fig. 2c; Supplementary Information). With both electrode configurations f e is inversely proportional to the droplet volume. However, if the droplet volume is > 158 pl (~67 μm diameter), the cross-side electrode configuration allows the highest sorting rate, whereas with droplets of < 158 pl the same-side electrode configuration allows the highest sorting rate.
The superiority of the cross-side electrode configuration for sorting 10 nl droplets was confirmed experimentally. Sorting of 10 nl droplets was achieved when AC electric field pulses of 1.4 kV pp , 30 kHz, 30 ms, were applied with the cross-side electrode configuration (Supplementary Movie S3). In contrast, the same electric field pulses could not sort droplets with the same-side electrode configuration (Supplementary Movie S4). The efficiency of the sorting device with the cross-side electrode configuration was further investigated by sorting 10 nl droplets of a binary emulsion comprising low red fluorescent droplets containing 1 μM sulforhodamine B (negative; 65%) and highly red fluorescent droplets containing 30 μM sulforhodamine B and blue ink (positive; 35%) ( Fig. 2d). The false positive and false negative error rates were determined either by color imaging of the emulsion before and after sorting or by high speed video analysis of droplets trajectories during the sorting process at different throughputs. When droplets were sorted at 4 s −1 , the emulsion collected in the positive channel was composed of 100% positive (blue) droplets while the emulsion collected in the negative channel contained 99.9% of negative droplets ( Fig. 2e; Supplementary Fig. S3). Based on the video analysis of the sorting process (Supplementary Movie S5), the false positive and false negative error rates were found to less than 0.07% and equal to 0.07%, respectively (Fig. 2f). These low error rates validate the ability of the cross-side electrode configuration device to efficiently sort 10 nl droplets based on fluorescence. The few false negative errors were always due to instability in the reloading regime of the droplets, resulting in temporary modification of the spacing at the sorting junction. The device could operate at up to 21 droplets.s −1 with acceptable false positive and false negative error rates ( Fig. 2f; Supplementary Movie S6), and was also able to sort 20 nl droplets with similar efficiency ( Supplementary Fig. S3). This is close to the maximum theoretical sorting frequency, f e , for 10 nl droplets calculated from the model (Fig. 2c), which was 3 droplet.s −1 for same-side electrode configuration and 46 droplet.s −1 for the cross-side electrode configuration.
Droplet-based microfluidics screening of filamentous fungi. The screening platform for Aspergillus niger (Fig. 3a) was composed of two distinct microfluidic devices ( Supplementary Fig. S1). A drop maker allowed the encapsulation of single spores with a fluorogenic α-amylase substrate (Supplementary Movie S2). Droplets of 18 nl were produced at ~80 Hz by hydrodynamic flow focusing of a spore suspension with a fluorinated oil phase containing 2.5% (w/w) KryJeffD 900 fluorosurfactant. Fluorinated oils have been shown to facilitate respiratory gas delivery to both prokaryotic and eukaryotic cells in culture 27 . The number of spores per droplet follows a Poisson distribution 28 and was controlled by adjusting the initial density of the spore suspension to give an average number of spores per drop λ, of 0.1-0.3. The droplets were collected in a glass capillary and incubated for 24 h at a 30 °C to allow germination, hyphal growth and α-amylase secretion within the droplets. After incubation, the droplets were re-injected into a sorting device with cross-side configuration electrodes (Supplementary Fig. S1 and Movie S7) which was used to sort droplets based on α-amylase activity (Supplementary Movie S8).
The fluorogenic substrate was starch labeled with multiple BODIPY ® FL fluorophores which are auto-quenched until the starch is hydrolyzed, resulting in an increase in green fluorescence (Fig. 3b). Sorted fungi were recovered, sporulated and further characterized. The microfluidic platform was applied to screen for Aspergillus niger strains which overproduce α-amylase from large whole-genome mutated libraries. The Aspergillus niger strain O58, a Soufflet α-amylase producing reference strain was exposed either to a chemical agent (2% N-Methyl-N′ -Nitro-N-Nitrosoguanidine [MNNG]) or to UV light (254 nm, 1.30 mW.cm −2 ) to introduce mutations in the entire genome. The diversity of chemical and UV libraries was 9.10 6 and 5.10 4 , respectively. The α-amylase activity of the O58 strain, the chemical library (2 and 3 h exposure) and the UV library (30, 60 and 90 s exposure) were analysed with the microfluidic platform ( Fig. 4a and Supplementary Fig. S4). The phenotypic signature of the O58 strain presents an active population [6.2% of drops (theoretically 9.5% with λ = 0.1); 2.5 ± 17.6% RFU] well separated from empty drops or drops showing no α-amylase activity [93.8% of drops; 1 ± 7% RFU] (Fig. 4a). Exposure to MNNG for 2 h resulted in an only slightly modified phenotypic signature [23% of active drops (theoretically 19.7% with λ = 0.22); 2.42 ± 15.7% RFU] (Fig. 4a). In contrast, the phenotypic signature was heavily impacted by 60 s UV-exposure, and the library predominantly comprised variants exhibiting no or low α-amylase activity and only a small number of variants with phenotypes giving equal or higher green fluorescence (α-amylase activity) than the O58 strain [2.6% of active drops (theoretically 19.7% with λ = 0.22); 1.93 ± 50.3% RFU] (Fig. 4a).
Based on these results, the UV mutagenesis method (60 s exposure) was chosen for the subsequent screening steps. The high mutation rates with UV mutagenesis result in high mortality (60 s UV exposure results in a survival rate of 2.6% ( Supplementary Fig. S4)) and a large number of inactive clones. However, a large fraction of inactive clones is also statistically associated, in a random whole-genome mutagenesis strategy, to a higher probability of targeting the desired pathway or gene.
A library was constructed using UV-mutagenesis (60 s exposure) of the O58 mother strain. The library was encapsulated in 18 nl droplets with the fluorogenic α-amylase substrate. The number of viable spores per droplets (21% occupancy; λ = 0.24) was controlled by adjusting the spore density based on the survival curve ( Supplementary Fig.  S4). After 24 h incubation at 30 °C, droplets were sorted at 10 Hz based on fluorescence, selecting droplets containing α-amylase producing fungi ( Fig. 4b; green region) (Movie S8). Around 5.10 4 droplets (~10 4 fungi) were screened in about 90 min, 1.45% were sorted (~750 droplets) and 616 fungi were recovered on an agar plate. The enrichment of the screening steps was evaluated by analyzing sorted fungi at the single clone level (~7,000 fungi) (Fig. 4c). The activity distributions show an important (196-fold) enrichment for active fungi: more than 97% of sorted fungi showed α-amylase activity, the enrichment being limited by co-encapsulation events during the emulsification step 12 .
Microtiter-plate and shake-flask fermentation assays. The sorted fungi population was screened using a microtitre plate-based robotic platform on an industrial relevant solid medium based on rapeseed meal. 1296 colonies were picked into 24-well plates, grown for 5 days and secreted α-amylase activity was measured in 96-well plates using the BODIPY ® FL based substrate ( Fig. 5a and Supplementary Fig. S5). 144 replicates of the O58 fungi were screened as positive controls. In total, 1,272 sorted fungi were analysed (87.2% of the sorted population): 98.9% were founded to be active while only 1.9% did not show detectable α-amylase activity. Most of the sorted fungi showed an activity level similar to O58 (47.3%) while 37.3% showed a higher activity (> 1σ O58 ). We measured the skewness of each distribution as the Pearson coefficient of skewness. We obtained a value of − 0.5 for the O58 distribution, meaning that the distribution is not significantly skewed compared to a normal distribution (p = 0.0124; α = 0.01), and a value of 1.05 for the sorted fungi distribution, indicating that the distribution has a significant positive skew (p < 0.00001; α = 0.01) (Supplementary Information).
We analysed in replicate 10 strains from the plate screen with α-amylase activity ≥ 3 by larger scale shake-flask fermentation. This format is closer to industrially relevant conditions and normally precedes further evaluation and scale-up in fermenters. The α-amylase activity was measured after 6 days of fermentation in PGS medium (100 ml) and compared to the O58 strain activity (Fig. 5b). All the strains tested except one displayed activity similar to or higher than O58, the best (M8) being 2.3-fold higher than the parental strain. The activities measured did not correlate exactly with those from microtitre plates, however, probably due the different production conditions (solid-state fermentation on rapeseed meal versus shake-flask fermentation in liquid synthetic PGS medium).

Discussion
We have developped and applied an efficient new approach for the high-throughput screening of filamentous fungi based on secreted enzyme activities. Using droplet-based microfluidic tools adapted to manipulate nanoliter-volume droplets, the technique allows fast enrichment of large population of fungi and can be integrated to industrial screening processes to speed up the discovery of new production strains of biotechnological interest.
Compared to automated colony pickers and liquid handling systems, the microfluidic platform offers higher throughputs, lower reagent consumption and reduced spatial footprint, all of which contribute to lower screening costs (Supplementary Table S1): screening 10 4 A. niger variants directly using the robotic microtitre plate-based screening platform would have taken more than 16 days and cost $8,770, whereas screening 10 4 variants using the microfluidic system took less than 24 h and cost $14 (and the subsequent screening of 1,272 sorted variants using the microtitre plate-based system took only 8 days and cost $1,100), costs being estimated based on consumables only. In addition to inherent throughput and scale advantages related to the microfluidics format, this approach also greatly simplifies the manipulation of filamentous fungi: off-chip fungi handling is limited to the manipulation of spores and operations on mycelia are performed by manipulating droplets.
Combining microfluidic and robotic microtitre-plate HTS is an attractive strategy. Here, single fungal spores were compartmentalized in droplets together with a fluorogenic α-amylase substrate, and the droplets and the germinated fungi they contain were sorted after 24 h incubation, triggered on fluorescence. However, with such an endpoint measurement, it is difficult to distinguish reliably between variants exhibiting different levels of α-amylase activity as fungal growth, enzyme secretion and the fluorogenic α-amylase reaction are coupled, and in many droplets all (or most) of the substrate was consumed. However, sorting all variants with a fluorescent signal equal to or greater than the mother strain using the microfluidic system efficiently removed variants showing reduced or no α-amylase activity: before sorting, 83% of clones showed little or no α-amylase activity, but a single round of microfluidic sorting resulted in a 196-fold enrichment of active clones. The throughput of microfluidic HTS enables the screening of highly mutated libraries, which are mainly composed of inactive clones, but where there is also a higher probability of finding beneficial mutations, as the higher mutation rate increases the chance of hitting the target gene(s).
Furthermore, the microfluidic assay could be made more quantitative by uncoupling fungal germination, growth and enzyme secretion from the enzymatic assay. This could be achieved by adding the fluorogenic substrate to droplets after incubation to allow fungal growth and enzyme secretion using droplet electrocoalescence 29 or picoinjection 30 . It may also be possible to increase the throughput above ~7,000 fungi.h −1 . The maximum experimental sorting rate (21 Hz) is close to the calculated theoretical maximum sorting rate limited by DEP and electrocoalescence (46 Hz). However, the sorting frequency was not limited by DEP forces or electrosplitting, but by purely fluidic constraints: higher reinjection frequencies resulted in droplets breaking when oil was added to space the reinjected droplets. Higher sorting rates could potentially be achieved by reducing the oil spacer flow rate, a strategy which, combined with other sorter design improvements, allowed sorting of pl volume droplets at 30 kHz 31 . Throughput could be further increased, at least to a certain extent, by increasing λ (at the expense of increasing the number of false positives due to co-compartmentalization of more than one variant in the same droplet 12 ).
The microfluidic system is not limited to screening for α-amylase activity: filamentous fungi can potentially be screened for the production of other enzymes or metabolites. The only restriction is that there must be a fluorogenic assay for the desired enzymatic activity or metabolite, and the metabolite and fluorescent product must not exchange between droplets over the time of the experiment. The use of fluorinated oils for the continuous phase mitigates against exchange, since non-fluorinated molecules are highly insoluble in fluorinated oils 32 . However, exchange can also be mediated by micellar transport 33 . Nevertheless, many microtitre-plate-based assays can be transposed directly to droplets, as was the case with the α-amylase assay used here, and direct, or coupled enzymatic assays can be used. Many other fluorogenic assays can also be simply adapted for microfluidic systems by chemically modifying the fluorogenic substrate to increase the hydrophilicity of the fluorescent product 34 , which reduces partitioning into micelles in the continuous phase and exchange between droplets 33 . For example, Figure 5. Microtiter-plate and shake-flask fermentation assays. (a) Relative α-amylase activity of the O58 strain (144 replicates) and the strains selected using the microfluidic screening (1,272 clones), analyzed using a robotic microtitre plate screening system. The standard deviation σ of the O58 strain is indicated. (b) Histogram displaying the α-amylase secretion level of 10 of the most active fungi from the plate screen and the wild type O58 strain after 6 days of cultivation at 30 °C in 100 ml of PGS medium. Mutants (gray) and O58 (red) were ranked from the lowest to the highest α-amylase activity. All fermentations were performed in duplicate and error bars correspond to ± 1 standard deviation.
Scientific RepoRts | 6:27223 | DOI: 10.1038/srep27223 droplet-based microfluidic assays have been developed to screen for bacteria and yeast producing cellulases, based on either directly screening for exogluconase (cellobiohydrolase) activity using a fluorogenic assay 35 or using a coupled fluorogenic enzyme assays to measure endogluconase activity 16 , and have been used for the ultrahigh-thoughput bioprospection of natural cellulolytic bacteria 35 . In the future, novel assays based on, for example, changes in osmotic pressure 36,37 or mass spectroscopy [38][39][40] might also be developed.
The ability to screen large filamentous fungi populations in a few hours at low cost should bring enormous benefit for identification of new biotechnologically interesting filamentous fungi. Besides this promising biotechnological application, this screening platform would also be applicable to more fundamental studies in combination with whole genome sequencing to help identify mutations and thereby regulatory networks responsible for improved phenotypes. The understanding of those genetic factors might be used to optimize filamentous fungi via inverse metabolic engineering 41 .

Methods
Fabrication of microfluidic devices. Poly-(dimethylsiloxane) (PDMS) microfluidic devices were fabricated as previously described 42  Optical setup, data acquisition and control system for droplet analysis. The optical setup used to monitor microfluidics experiments was previously described, as well as the data acquisition and control system 34 .
In addition, to allow the sorting of a particular droplet, the data acquisition card provided a signal to a model 623B high-voltage amplifier (Trek Inc.) connected to the electrodes of the microfluidic device.

Fluorescence-activated droplet sorting validation.
A binary emulsion comprising 10 nl droplets was produced by mixing two emulsions produced successively using the dropmaker device. The first droplet population contained 1 μM of sulforhodamine B in YGC medium (Yeast extract 5 g.l −1 , Glucose 20 g.l −1 , Chloramphenicol 0.1 g.l −1 ) and the second contained 30 μM of sulforhodamine B and blue ink colorant in YGC medium. The binary emulsion was reloaded into the sorting device and sorted as a function of red fluorescence to select only those droplets containing 30 μM of sulforhodamine B. Both sorted and unsorted droplets were collected. The performance of the sorter was evaluated before and after sorting either by imaging the emulsion on a glass slide using a color camera (Nikon DS-Fi2) or by video analysis of the sorting process using a high-speed camera (Mikroton Eosens MC1310).
Filamentous fungi strain and culture conditions. The Aspergillus niger O58 strain was obtained from the Ets J. Soufflet strain collection. The O58 strain was plated on a PDA (AES Laboratoire) plate and grown at 30 °C for 10 to 15 days for optimal sporulation. Spores were suspended in 20 ml of PGS medium (Glucose 10 g.l −1 , Pancreatic peptone 6 g.l −1 , MgSO 4 .7H 2 O 0.5 g.l −1 , KH 2 PO 4 0.5 g.l −1 , FeSO 4 .7H 2 O 0.5 mg.l −1 ). The spore suspension was stirred for 10 min, filtered through sterilized gauze to remove mycelium residues and stirred again for 10 min. The spore density was estimated using a Thoma cell counting chamber. UV mutation. After sporulation, spores were suspended in 0.01% Triton X-100 (20 ml). The spore suspension was stirred for 10 min, filtered through sterilized gauze and stirred again for 10 min with glass beads to Scientific RepoRts | 6:27223 | DOI: 10.1038/srep27223 homogenize. The spore concentration was adjusted to 3.10 6 spores per ml. 3 ml of the spore suspension was spread on a glass petri dish and exposed to 6 W UV light (Dymax PC-2000/38003) at 254 nm for 60 s at a distance of 5 cm from the lamp (1.30 mW.cm −2 ) to introduce mutations randomly throughout the whole fungi genome. The procedure was repeated three times to UV treat a total of 9 ml of spore suspension. The UV-treated spore suspension was centrifuged 10 min at 2500 g and the supernatant was discarded. The pellet was suspended in PGS medium.

MNNG mutation.
As for UV mutation, a spore suspension was prepared. The spore concentration was adjusted to 1.10 7 spores per ml. 5 ml of 0.2 M phosphate buffer (pH 6.2) containing 6% (w/w) of N-Methyl-N′ -Nitro-N-Nitrosoguanidine (MNNG) was added to 10 ml of spore suspension. The suspension was stirred at 30 °C for the desired mutagenesis time. Mutated spores were then centrifuged 5 min at 4500 g and the pellet was suspended in 5 ml of Triton X-100 (1 g.l −1 ). The procedure was repeated three times to wash away the mutagenesis agent. The pellet was finally suspended in PGS medium.
Characterization of α-amylase activity of O58 strain and libraries. The spore suspension was diluted to the appropriate spore to droplet ratio to obtain typically 0.1-0.2 spore per droplet in PGS medium containing the BODIPY ® FL-labeled DQ TM starch substrate (50 μg.ml −1 ; Enzchek Ultra Amylase Assay Kit; Life Technologies) and Sulforhodamine 101 (5 μM). Spores were encapsulated in 18 nl droplets using the dropmaker device. The emulsion was collected in a glass capillary, incubated off-chip for 24 h at 30 °C and then reloaded in a FADS device for fluorescence analysis.
Microfluidic screening of libraries for α-amylase production. The library spore suspension was encapsulated in 18 nl droplets as described in the previous paragraph. The emulsion was collected in a glass capillary, incubated off-chip for 24 h at 30 °C and then reloaded in a FADS device for fluorescence analysis and selection of droplets containing active fungus at 10 droplets.s −1 by applying AC field pulses (30 kHz; 1600 V pp ; 35 ms). Sorted fungi were spread on PDA plates, grown for 48 h at 30 °C, counted, and then grown for an additional 6 days to allow sporulation. For conservation, spores were suspended in 20 ml of water, counted on a Thoma cell counting chamber and frozen at − 80 °C in 15% glycerol at a density of at least 10 6 spores per ml.