Diversity and Ice Nucleation Activity of Microorganisms Collected With a Small Unmanned Aircraft System (sUAS) in France and the United States

Many microbes relevant to crops, domestic animals, and humans are transported over long distances through the atmosphere. Some of these atmospheric microbes catalyze the freezing of water at higher temperatures and facilitate the onset of precipitation. We collected microbes from the lower atmosphere in France and the United States with a small unmanned aircraft system (sUAS). 55 sampling missions were conducted at two locations in France in 2014 (an airfield in Pujaut, and the top of Puy de Dôme), and three locations in the U.S. in 2015 (a farm in Blacksburg, Virginia, and a farm and a lake in Baton Rouge, Louisiana). The sUAS was a fixed-wing electric drone equipped with a remote-operated sampling device that was opened once the aircraft reached the desired sampling altitude (40–50 meters above ground level). Samples were collected on agar media (TSA, R4A, R2A, and CA) with and without the fungicide cycloheximide. Over 4,000 bacterial-like colonies were recovered across the 55 sUAS sampling missions. A positive relationship between sampling time and temperature and concentrations of culturable bacteria was observed for sUAS flights conducted in France, but not for sUAS flights conducted in Louisiana. A droplet freezing assay was used to screen nearly 2,000 colonies for ice nucleation activity, and 15 colonies were ice nucleation active at temperatures warmer than −8°C. Sequences from portions of 16S rDNA were used to identify 503 colonies from 54 flights to the level of genus. Assemblages of bacteria from sUAS flights in France (TSA) and sUAS flights in Louisiana (R4A) showed more similarity within locations than between locations. Bacteria collected with sUAS on TSA in France and Virginia were significantly different across all levels of classification tested (P < 0.001 for class, order, family, and genus). Principal Coordinates Analysis showed a strong association between the genera Curtobacterium, Pantoea, and Pseudomonas from sUAS flights in Virginia, and Agrococcus, Lysinibacillus, and Paenibacillus from sUAS flights in France. Future work aims to understand the potential origin of the atmospheric microbial assemblages collected with sUAS, and their association with mesoscale atmospheric processes.

Many microbes relevant to crops, domestic animals, and humans are transported over long distances through the atmosphere. Some of these atmospheric microbes catalyze the freezing of water at higher temperatures and facilitate the onset of precipitation. We collected microbes from the lower atmosphere in France and the United States with a small unmanned aircraft system (sUAS). 55 sampling missions were conducted at two locations in France in 2014 (an airfield in Pujaut, and the top of Puy de Dôme), and three locations in the U.S. in 2015 (a farm in Blacksburg, Virginia, and a farm and a lake in Baton Rouge, Louisiana). The sUAS was a fixed-wing electric drone equipped with a remote-operated sampling device that was opened once the aircraft reached the desired sampling altitude (40-50 meters above ground level). Samples were collected on agar media (TSA, R4A, R2A, and CA) with and without the fungicide cycloheximide. Over 4,000 bacterial-like colonies were recovered across the 55 sUAS sampling missions. A positive relationship between sampling time and temperature and concentrations of culturable bacteria was observed for sUAS flights conducted in France, but not for sUAS flights conducted in Louisiana. A droplet freezing assay was used to screen nearly 2,000 colonies for ice nucleation activity, and 15 colonies were ice nucleation active at temperatures warmer than −8 • C. Sequences from portions of 16S rDNA were used to identify 503 colonies from 54 flights to the level of genus. Assemblages of bacteria from sUAS flights in France (TSA) and sUAS flights in Louisiana (R4A) showed more similarity within locations than between locations. Bacteria collected with sUAS on TSA in France and Virginia were significantly different across all levels of classification tested (P < 0.001 for class, order, family, and genus). Principal Coordinates Analysis showed a strong association between the genera Curtobacterium, Pantoea, and Pseudomonas from sUAS flights in Virginia, and Agrococcus, Lysinibacillus, and Paenibacillus from sUAS flights in France. Future work aims to understand the potential origin of the atmospheric microbial assemblages collected with sUAS, and their association with mesoscale atmospheric processes.

INTRODUCTION
Many microorganisms that impact the health of crops, domestic animals, and humans are transported over long distances through the atmosphere (Schmale and Ross, 2015). The dispersal of these microbes in the atmosphere is comprised of a complex series of events (Isard and Gage, 2001). Such events include microbial life history stages (e.g., aerosolization, escape from terrestrial or aquatic environments, survival while airborne, deposition onto a suitable environment, etc.), that interact with biotic and abiotic factors (e.g., rainfall, air currents, UV radiation, etc.) that control dispersal (Isard et al., 2005). These events and factors change in time and space across large-scale biological and meteorological gradients, forming a framework to understand the aerobiological processes governing the spread of microorganisms in the atmosphere (Aylor, 1986;Schmale and Ross, 2015).
Some atmospheric microbes catalyze the freezing of water at higher temperatures and may facilitate the onset of precipitation (Burrows et al., 2009;Murray et al., 2012;Morris et al., 2014). Such ice-nucleating microorganisms likely play an important role in the water cycle (Sands et al., 1982;Morris et al., 2014). Bacteria in the atmosphere may be associated with cloud formation (Ahern et al., 2007;Joly et al., 2013;Hu et al., 2017). A number of studies have described the ice nucleation activity (INA) of microorganisms collected from different components of the water cycle (Burrows et al., 2009;Murray et al., 2012;Morris et al., 2013). It has been hypothesized that the INA phenotype emerged as an adaptation to harsh conditions in the atmosphere; freeze tolerance and/or microbes facilitate their own deposition from the atmosphere (with low nutrient concentration) to land (rich in nutrients) through the production of ice-nucleating proteins or macromolecules .
An average of 2 × 10 4 fungal spores and bacteria m −3 was reported over many continental regions (Spracklen and Heald, 2014). Residence times of about 3-17 days could be enough to allow the transport of microbes over thousands of kilometers (Mayol et al., 2017). Islands, such as those in the West and Central Pacific, may act as stepping stones for the long-range transport of microorganisms across oceans (Mayol et al., 2017). The modeling of Lagrangian coherent structures (LCS) may help describe the long distance transport of microorganisms within and among habitats (Tallapragada et al., 2011).
New technologies are needed to study the movement of microorganisms in the atmosphere. Previous studies (Schmale et al., 2008(Schmale et al., , 2012Lin et al., 2013Lin et al., , 2014 have used unmanned aircraft systems (UASs) to study the transport of microorganisms tens to hundreds of meters above the ground. These UASs are equipped with unique devices for collecting microbes in the atmosphere during flight. Autonomous systems enable teams of UASs to perform complex atmospheric sampling tasks, and coordinate flight missions with one another colleagues (Techy et al., 2010). Data collected with UASs can be used to validate and improve disease forecasting models along highways in the sky, connecting transport scales across farms, states, and continents (Schmale and Ross, 2015). Though terrestrial environments are often considered a major contributor to atmospheric microbial aerosols, little is known about aquatic sources of microbes. Droplets containing microorganisms can aerosolize from the soil and water surface, liberating them into the atmosphere (Pietsch et al., 2015;Joung et al., 2017). Aerosolized droplets may be produced in several ways, including bubble-bursting and fragmentation from breaking waves (Wu, 1981). Fluorescent dyes have been released as surrogate pathogenic microorganisms, and tracked with teams of unmanned systems in the water and in the air (Powers et al., 2018a).
The overall goal of our work was to measure the diversity and ice nucleation activity of microbes in the lower atmosphere. The specific objectives of our work were to: (1) collect microbes from the lower atmosphere at several sites with a small unmanned aircraft system (sUAS), (2) quantify culturable microorganisms from sUAS flights using different types of agar media to maximize the diversity of microorganisms detected, (3) determine the frequency of INA microbes collected with sUAS at −8 • C, (4) identify bacteria from sUAS missions using sequences from portions of 16S rDNA, and (5) explore associations of microbial diversity to the different flights (55), locations (Virginia, Louisiana, and France), and media (R4A, TSA, R2A, and CA).

MATERIALS AND METHODS
GPS locations of field studies are reported below in decimal degrees (DD). Field studies were conducted at two locations in France in 2014; an airfield in Pujaut (43.994737,4.750301), and the top of Puy de Dôme (45.773260,2.964211). Field studies were also conducted at three locations in the U.S. in 2015; a farm in Blacksburg, Virginia (37.196049 Table 1). These locations were selected based on their proximity to potential terrestrial and aquatic sources of ice nuclei.

Small Unmanned Aircraft System
A small unmanned aircraft system (sUAS) was used to collect microorganisms in the lower atmosphere. The sUAS was an electric fixed-wing drone (the Clouds Fly AXN) with a remotecontrolled microbe-sampling device (Figure 1, top). The sUAS was equipped with the Hitec Sensor Station (HTS-SS) including pitot airspeed (HTS-AS) and GPS (HTS-GPS) sensors. The microbe-sampling device consisted of a single 60 mm verticallymounted petri plate containing agar media (Figure 1, top). This is a modification of the design used by Schmale et al. (2008), Schmale et al. (2012), and Lin et al. (2013Lin et al. ( , 2014. The sampling device was closed during takeoff and landing and was opened by remote control from the ground once the sUAS was at the target sampling altitude and airspeed. The device remained open for the duration of the sampling interval which was usually 10 minutes, as detailed in Table 1. Immediately following sample collection, the exposed plate containing the agar media was removed from the sampling device and stored in a small plastic container for transport to the laboratory. CFUs per flight were converted to CFU/m 3 of air sampled as described previously (Aylor et al., 2011;Lin et al., 2014).
Frontiers in Microbiology | www.frontiersin.org FIGURE 1 | A small unmanned aircraft system (sUAS) was used to collect samples 50 m above ground level. The microbe-sampling device consisted of a single 60 mm vertically-mounted petri plate containing agar media. The sampling device was closed during takeoff and landing, and was opened by remote control from the ground once the sUAS was at the target sampling altitude and airspeed (top and middle panels). The device remained open for the duration of the 10-min sampling interval, and was closed at the end of the sampling interval prior to landing. Immediately following sample collection, the exposed plate containing the agar media was removed from the sampling device and stored in a small plastic container for transport to the laboratory. Plates were incubated for 3-4 days at room temperature, and colony forming units (CFUs) were counted and colonies (bottom panel) were picked and transferred to new media. Left plate was from sUAS mission LA5, and right plate was from sUAS mission LA6.

Collection and Culturing of Microorganisms
Agar collection plates from the sUAS sampling missions contained TSA, R2A, R4A, and CA media with or without 50 mg/L cycloheximide. These media were selected based on their use in previous published reports Hanlon et al., 2017), and in an effort to maximize the potential diversity of microorganisms cultured from environmental samples. Plates were incubated for 3-4 days at room temperature, and colony forming units (CFUs) were counted and the plates were photographed (Figure 1, bottom). Colonies were picked with sterile toothpicks and streaked onto the same media they were collected on to obtain pure cultures in a laminar flow hood under sterile working conditions. Pure cultures were stored at −80 • C in 20% glycerol. Individual microbes from France flights were inoculated onto TSA stabs in 1.5 mL Eppendorf tubes for transport back to the laboratories at Virginia Tech, USA. The French collection was cultured on fresh TSA plates prior to INA screening. Flight plates from the LSU and Blacksburg locations were incubated at room temperature, approximately 22 • C, for 4-7 days and then stored at 4 • C. Microbes were moved as individual isolates to fresh streak plates, or chosen directly from flight plates for INA screening. Samples were then picked as individual microbes, from streak plates or flight plates, and inoculated directly into the 96-well plate for droplet freezing assays.

Weather Data
Historical hourly weather data for each of the regions where the sUAS sampling was conducted were obtained from https://www.wunderground.com/ ( Table 2). For Louisiana, weather data were obtained from the Baton Rouge Metro-Ryan airport (https://www.wunderground.com/history/airport/ KBTR/). For Virginia, weather data were obtained from the Virginia Tech airport (https://www.wunderground.com/history/ airport/KBCB/). For Pujaut, France weather data were obtained from the Caumont airport (https://www.wunderground.com/ history/airport/LFMV/). For Puy de Dôme, France weather data were obtained from the Clermont Ferrand airport (https://www. wunderground.com/history/airport/LFLC/).
Associations of concentrations of culturable bacteria (CFUs/m 3 ) with each of the weather variables for sUAS flights in France and Louisiana (where the same media was used across each of the flights at each location) were examined using regression analysis.

Ice Nucleation Assays
Ice nucleation assays were performed in 96-well plates. Pure cultures of microbes were transferred from culture plates with a sterile toothpick to 140 µL of DI water. Plates were covered with pierceable foil (LMT-ALUMA-II, Phenix Research Products Inc, Candler, NC) to prevent cross-contamination across wells and incubated at 4 • C for 1 h. The foil was pierced with a multichannel pipet to load the cooling bath, and a second pierceable foil was placed on top after 12 µL sample aliquots were loaded. Droplet freezing assays were performed on a PARAFILM R M (Sigma P6543, 20 in. × 50 ft) float assembly in a Lauda Alpha RA 12 (LCKD 4908) cooling bath (LAUDA-Brinkmann, LP, Delran, NJ, 08075) with ethylene glycol coolant fluid (Air gas RAD64000246). Droplets (12 µL) of the microbial suspension were loaded in duplicate, and the temperature of the cooling bath was decreased stepwise every 2 minutes from −2 • to −12 • C, or until frost spontaneously formed on the float. Microbes were considered to be INA if at least one drop froze at or before −8 • C. These putative INA strains identified from the original test were grown individually on their respective media type and stored in 20% glycerol at −80 • C. These isolates were plated from frozen glycerol stocks to test for INA a second time ( Table 3).

Identification of Microorganisms Using Portions of 16S rDNA
All cultures (regardless of INA phenotype) from cryogenic storage were streaked on R2A plates to ensure isolates were pure prior to sequencing. When a mix of different colonies was isolated from a cryogenic sample, additional streaks were performed to ensure the purity of each individual sample. Up to 15 isolates per flight were randomly selected (using assignments based on a random number generator) from each flight collection for sequence-assisted identification. Colony PCR was used to amplify a portion of 16S rDNA. Different colony PCR strategies were performed on the cultures including: (1) transferring the colony directly to a PCR tube to serve as a template, (2) transferring a tooth-pick microbial biomass to 100 µL of distilled water, and then using 5 µL of this dilution as the template for PCR, (3) heating 10 µL of the cited biomass solution and taking 2 µL as the template, and (4) using a DNA extraction kit (Puregene Yeast/Bact. Kit B; Qiagen #1042607) with microbial glycerol stock as the template for PCR.
A portion of the 16S rDNA gene was amplified using the primers 518 F (5 ′ CCA GCA GCC GCG GTA ATA CG 3 ′ ) and 1491R (5 ′ ATC GGY TAC CTT GTT ACG ACT TC 3 ′ ) (Turner et al., 1999) and http://www.earthmicrobiome.org/empstandardprotocols/16s/). GoTaq R Green Master Mix (Promega M712) was used with 2-5 µL of template in a 20 µL PCR reaction. The parameters used in the thermal cycler were: 1 × cycle denaturation 95 • C for 10 min, 5 × cycle denaturation 95 • C for 1 min, annealing 53 • C for 30 s, extension 72 • C for 1min 40 s, 25 × cycle 94 • C for 30 s, 53 • C for 20 s, 72 • C for 1 min, 1 × cycle 72 • C for 10 min, hold 4 • C). Resulting PCR products were loaded on an 1% TBE agarose gel with Ethidium bromide and visualized on a UV transilluminator. Enzymes rSAP and ExoI were used to clean-up the PCR products, and the cleaned products were sequenced by Eton Biosciences (AB1 3,730 × 1 DNA Sequencer, Eton Biosciences, 104 T.W. Alexander Drive, Bldg 4A, RTP, NC 27709, USA). The 4 Peaks program (Version 1.8. The Netherlands) was used to trim the resulting sequences with a standard quality (QC) of 30 or higher. The trimmed sequences then were queried against the NCBI nr (non-redundant) database to identify the microorganisms to the level of genus. Sequences were also tested to find chimeras using the DECIPHER program, Version 2.0 (Wright et al., 2012). Accession numbers were obtained from GenBank following submission using the BankIt submission tool at NCBI (ncbi.nlm.nih.gov).

Analysis of Microbial Diversity Among sUAS Collections
Identities from sequences from portions of 16S rDNA were used to associate each isolate with class, order, family, and genus. Media effects (CA, R2A, TSA) were considered for flights conducted in Virginia, since this was the only location for which three media types were used. Location effects were ascertained by comparing isolates from TSA media for France and Virginia. Null hypothesis tests of identical assemblage composition among treatments (i.e., media types and location) were run using PERMANOVA (Anderson, 2001) for the taxonomic levels of class, order, family, and genus using 10,000 permutations. Patterns of genera composition were graphically summarized using PCoA (Principal Coordinates Analysis, i.e., Metric multidimensional scaling) (Kruskal, 1964a) (Kruskal, 1964b). To account for presence/absence nature of data, Jaccard's presence/absence index (Jarccard, 1908) was used to create required dissimilarity matrices for analyses. Vector fitting (Oksanen, 2015) was used to depict the strength of association between the ordination projection and the genera that underlay the projection. Vector fitting tests for the null hypothesis of no association were implemented using 1,000 permutations.

RESULTS
Fifty-five sampling missions were conducted at two locations in France in 2014, and three locations in the U.S. in 2015. Over 4,000 colonies were recovered across the 55 sUAS sampling missions ( Table 1).

Concentrations of Microbes
CFUs per flight were converted to CFU/m 3 of air sampled as described previously (Aylor et al., 2011;Lin et al., 2014). For flights conducted in Virginia, concentrations of microbes ranged from 2 to 53 (Flight CF23) CFUs/m 3 (Table 1, Figure 2). For flights conducted in Louisiana, concentrations of microbes ranged from 3 to 79 (Flight LA19) CFUs/m 3 (Table 1, Figure 2). For flights conducted in France, CFUs ranged from 0 to 43 (Flight FRA09) CFUs/m 3 (Table 1, Figure 2). Concentrations of microbes were generally lower for media containing the fungicide cylohexamide (CHX) (Figure 2). For sUAS flights in France, there was a significant positive relationship between time of the sUAS flight (r = 0.61, P < 0.01) (Figure 3) and temperature during the sUAS flight (r = 0.71, P < 0.01) and CFU/m 3 . There was a significant negative relationship between relative humidity during the sUAS flight and CFU/m 3 (r = 0.52, P = 0.02). There was also a significant negative relationship between relative humidity and temperature (r = 0.67, P < 0.01). There was no significant relationship between windspeed and CFU/m 3 (r = 0.27, P = 0.26).
For sUAS flights in Louisiana, there was no significant relationship between any of the variables tested (time of the sUAS flight, temperature, relative humidity, and windspeed) and concentrations of viable bacteria (CFU/m 3 ) (P > 0.05 for all variables tested).

INA Strains
The ice nucleation activity (INA) of the isolates collected from the sUAS flights was determined using a droplet freezing assay. A total of 1,963 isolates was tested for INA, and 15 isolates (∼1% of the total screened) were INA at −8 • C following the initial screen. However, only one of these isolates (LA6P5) showed INA following a second screen from cryogenic cultures. This isolate (LA6P5) was identified as Xanthomonas (GenBank accession number MG924929) ( Table 4). None of the control droplets containing DI water froze at −12 • C (the lower temperature limit) in any of the assays.

Diversity of Microbes in sUAS Collections
Graphical summaries of taxa in location/media combinations are given for class (Figure 4) and for genus (Figure 5). Analyses on the effect of media for sUAS flights in Virginia showed evidence of assemblage variation at the level of genus (P = 0.05). No evidence of assemblage variation with media was observed at the levels of class (P = 0.34), order (P = 0.34), and family (P = 0.12). Three genera were significantly associated with projection scatter. These were: (1) Curtobacterium, which was associated with CA and TSA (which showed tight clustering), (2) Microbacterium which was associated with three of the CA samples, and (3) Massilia which was associated with three of the R2A samples.
Atmospheric assemblages of bacteria collected with sUAS on TSA in France and Virginia were significantly different across all levels of classification tested (P < 0.001 for class, order, family, and genus). Principal Coordinates Analysis showed a strong association between the genera Curtobacterium, Pantoea, and Pseudomonas from sUAS flights in Virginia, and Agrococcus, Lysinibacillus, and Paenibacillus from sUAS flights in France (Figure 6).
Mean Jaccard genus-based dissimilarities of sUAS flights in France (TSA) (Figure 7) and sUAS flights in Louisiana (R4A) (Figure 8) showed more similarity within locations than between locations; 0.74 for France|France and 0.78 for Louisiana|Louisiana, and 0.82 for France|Louisiana (Koleff et al., 2003). This is also supported with PERMANOVA testing of the H 0 that there are no true community differences between France (TSA media) and Louisiana (R4A media) (P < 0.0001).

Flight
Strain ID Genus GenBank accesion number Strain LA6P5 showed INA following a second screen from cryogenic cultures.

DISCUSSION
Little is known about the transport, diversity, and ice nucleation activity of microorganisms in the lower atmosphere. We collected microbes from the lower atmosphere in France and the United States with a small unmanned aircraft system (sUAS). At least one culturable microorganism was recovered from all of the 55 sUAS sampling missions, with concentrations of culturable   air (Proctor, 1935). Mayol et al. (Mayol et al., 2017) reported airborne concentrations of about 8,000 per m 3 over the North Atlantic Ocean.
Over 4,000 CFUs were recovered across the 55 sUAS sampling missions on four different types of agar media (TSA, R4A, R2A, and CA) with and without the fungicide cycloheximide. Schmale and colleagues have reported collections of Fusarium from the atmosphere using large UAS (Tallapragada et al., 2011;Schmale et al., 2012;Lin et al., 2013Lin et al., , 2014, but most of these collections used a selective medium (a modified Nash-Snyder medium containing the fungicide PCNB and a cocktail of antibiotics) which biased the recovery of fungi in this genus. To our knowledge, our work represents the first published report of microbial exploration in the lower atmosphere using a handlaunched drone. Since our sUAS can be launched by hand, this work demonstrates that such sUAS can be deployed in remote locations lacking an established runway for takeoff and landing (a requirement for large fixed wing drones), thereby opening new perspectives for exploration of tropospheric aerobiology and its dynamics.
For sUAS flights in France, we observed a positive relationship between sampling time and temperature and concentrations of culturable bacteria. Thus, greater concentrations of bacteria were collected during UAS flights in France that were conducted later in the day and during higher temperatures. During daylight hours, the Earth's surface is heated by the sun, and there is an upward transfer of heat (Oke, 1987). This transfer of heat generates atmospheric turbulence and vertical mixing, stretching the planetary boundary layer kilometer distances (Oke, 1987). We speculate that our observed associations between temperature and sampling time with atmospheric concentrations of bacteria are driven in part by vertical mixing of bacteria from the surface of the Earth. Maldonado-Ramirez et al. (2005) collected viable spores of the fungus Fusarium graminearum over 158 UAS sampling flights in New York, USA over four years (Maldonado-Ramirez et al., 2005). Though these authors did not observe a significant difference between spore concentrations among UAS flights conducted during different times of the day, these authors did observe significant differences among UAS flights conducted during clear, cloudy, and rainy conditions. For sUAS flights in Louisiana, we did not observe a significant relationship between temperature and sampling time with concentrations of culturable bacteria (P > 0.05 for all variables). Temperatures were much greater during sUAS flights in Louisiana than in France, suggesting potential transport and/or viability limitations during the sampling conditions encountered in Louisiana. Environmental factors such as ultraviolet radiation, relative humidity, and temperature are known to have a significant impact on the viability of microorganisms in the atmosphere and thus may impact their ability to survive transport over long distances (Isard and Gage, 2001;Schmale and Ross, 2015). The identification and tracking of overlapping backtrajectories in high probability regions (Schmale and Ross, 2015) and Lagrangian coherent structures (LCS) (Tallapragada et al., 2011) could help generate hypotheses concerning potential sources and transport distances of the microbes collected during our sUAS missions, but ecological and environmental factors impacting their viability during transit must also be considered.
Atmospheric assemblages of bacteria from sUAS flights in France (TSA) and sUAS flights in Louisiana (R4A) showed more similarity within locations than between locations. Moreover, bacteria collected with sUAS on TSA in France and Virginia were significantly different across all levels of classification tested (P < 0.001 for class, order, family, and genus). Cuthbertson et al. (2017) suggested that day of sampling (temporal) is more important than location (spatial) with regards to sample diversity, most likely due to changes in meteorological conditions such as wind direction which appeared to produce distinct assemblages at the level of phylum. Bowers et al. (2009) sampled microbes from the atmosphere at a high-elevation site in Colorado, and reported that microbial abundances, diversity, and composition were fairly stable over a two-week sampling period (Bowers et al., 2009). Another study by Bowers et al. (2011), showed that the composition of the microbes in the sample was linked to the source of the sample (soil, plant leaf, freshwater, etc.) (Bowers et al., 2011). These authors also argued that land-use type was significantly related to the bacterial assemblage composition in the atmosphere, where the most abundant phyla were Proteobacteria, Actinobacteria and Firmicutes. Interestingly, these three phyla have also been reported from environmental samples collected in Japan (together with Cyanobacteria and Acidobacteria) (Hu et al., 2017), France (together with Bacteroidetes) (Amato et al., 2017), Svalbard (close to the North Pole) (Cuthbertson et al., 2017), and from our sUAS collections reported in this manuscript. Amato et al. (2017) reported diverse microbial assemblages in cloud water samples collected at the Puy de Dôme in France, potentially from a diverse range of sources. Teasing out contributions of local and more distant sources to the sUAS assemblages would require a coordinated sampling effort where multiple sUAS would sample at multiple sites simultaneously. Such an effort is needed to dissect local (the Pujaut, France site was surrounded by farmland, and Louisiana and Virginia sites were on farmland) and more regional sources which could be attributed to larger land-use trends.
Principal Coordinates Analysis showed a strong association between the genera Curtobacterium, Pantoea, and Pseudomonas from sUAS flights in Virginia, and Agrococcus, Lysinibacillus, and Paenibacillus from sUAS flights in France. We speculate that these associations are likely due to differences in potential sources and/or weather patterns in these regions (Monteil et al., 2014), but this was not tested in the current study. Such speculation is reasonable given potential linkages between microbial diversity in collections and the attribution of land-use types reported by others (Bowers et al., 2011;Monteil et al., 2014). The weather data we report surrounding our sUAS collections ( Table 2) suggest that there may be potential associations between microbial diversity, land-use, and weather, at least for France (where we observed significant relationships between temperature and time of sampling).
A droplet freezing assay was used to screen nearly 2,000 colonies for ice nucleation activity, and 15 colonies were INA at temperatures warmer than −8 • C following the initial screen. Only one of these isolates (Xanthomonas) was INA following a second screen. Failor et al. (2017) suggested that discrepancies in INA between initial screens of bacteria and repeated screens (following storage and re-plating) could be attributed to growth conditions not conducive to INA. Hanlon et al. (2017) performed INA assays on microbes collected from bulk rain samples from a series of simulated rain events, and found about 1.5% (121/8,331) of the colonies tested following an initial screen (Hanlon et al., 2017). However, only 0.4% (34/8,331) of microbial ice nucleators were INA following a second screen. Powers et al. (2018b) reported about 1% (7/720) of aquatic microbes collected from a freshwater lake in Virginia, USA (Powers et al., 2018b). Likewise, Lindeman and colleagues (Lindemann et al., 1982) showed that on average about 1.7% of the bacteria recovered in air samples above crops were INA and were not greater than 4% of the total viable bacteria detected in the air. On the other hand, about 16% (7/44) of the strains isolated from cloud water on the Puy de Dôme in France were INA (Joly et al., 2013). Failor et al. (2017) isolated ice+ bacteria from bulk rain samples (1.8% of the strains tested, or 593/33,134) from 23 different precipitation events on seven media types, and observed mostly ice+ strains in the class Gammaproteobacteria . Therefore, the relative abundance of INA bacteria in the total populations that we detected in air in this study were the same order of magnitude as in previous studies of the air, lake water and rainfall, but are an order of magnitude lower than in cloud water (Lindemann et al., 1982;Joly et al., 2013;Failor et al., 2017;Powers et al., 2018b).
Although the quantities of INA bacteria detected here seem very small, when translated into the processes that loft bacteria into the atmosphere and up to cloud heights, these amounts are likely to be very significant relative to cloud processes. Over crops, Lindemann and colleagues (Lindemann et al., 1982) observed net upward bacterial flux of up to 500 bacteria m −2 s −1 and Carotenuto and colleagues (Carotenuto et al., 2017) observed up to twice that rate for bacteria and fungi combined. Upward flux is sustained by sensible heat flux from sunlight that warms the ground. This usually occurs markedly between about 10:00 to 16:00 on sunny summer days. During this 6-h period, an upward flux of 500 bacteria m −2 s −1 would contribute a cumulative total of about 10 10 bacteria into the air over a 1 ha field (and therefore 10 8 INA bacteria if 1% of these were INA). This suggests that seemingly low concentrations of bacteria in the air (including those observed by sUAS) could nevertheless be of considerable importance as they accumulate over time, and that colossal sampling efforts are needed to capture and quantify them.
Future work aims to understand the potential origin of the atmospheric microbial assemblages collected with sUAS, and their association with mesoscale atmospheric processes. Such work could make use of known sources of ice-nucleating microorganisms in a series of release and recapture experiments (e.g., Prussin et al., 2014), and/or coordinated sampling with sUAS and ground-based collection schemes in periods surrounding major meteorological events such as hurricanes.

AUTHOR CONTRIBUTIONS
CJ-S cultured microbes, collected and analyzed 16S sequence data, and led the writing of the manuscript. RH conducted ice nucleation assays, managed the strain collection, provided technical advice related to the sequencing, analyzed data, and assisted in writing the manuscript. KA conducted statistical analyses using the 16S data and assisted in writing the manuscript. CP constructed the sUAS sampling devices, assisted in field collections with sUAS in Virginia and Louisiana, and assisted in writing the manuscript. CM helped coordinate field collections in France, analyzed data, and assisted in writing the manuscript. DS managed the project, designed experiments, was the pilot in command (PIC) for all sUAS missions, cultured microbes, analyzed data, and assisted in writing the manuscript.