Ken A Aho
email: ahoken@isu.edu

C F Weber

Brent C Christner

Boris A Vinatzer

Cindy E Morris

R Joyce

K C Failor

J T Werth

A L H Bayless-Edwards

David G Schmale

Schmale Iii

Spatiotemporal patterns of microbial composition and diversity in precipitation

Keywords: microbial ecology, community ecology, bioprecipitation, precipitation, atmosphere, metabarcoding, metacommunity, biogeography, meteorology

de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

Introduction

Of growing interest to ecologists are the properties of precipitation-borne microbes, particularly

Accepted Article

Version postprint

This article is protected by copyright. All rights reserved species and strains that may trigger precipitation through the formation of ice in clouds (i.e., bioprecipitation, [START_REF] Sands | The association between bacteria and rain and possible resultant meteorological implications[END_REF], Morris et al. 2014a[START_REF] Hummel | Simulating the influence of primary biological aerosol particles on clouds by heterogeneous ice nucleation[END_REF]. At temperatures warmer than -38 °C, phase transition from liquid water to ice requires the presence of ice nucleating particles (INPs) [START_REF] Atkinson | Rate of homogenous nucleation of ice in supercooled water[END_REF][START_REF] Kanji | Overview of ice nucleating particles[END_REF][START_REF] Knopf | The role of organic aerosol in atmospheric ice nucleation: a review[END_REF]. Various particles can act as ice nucleators [START_REF] Szyrmer | Biogenic and anthropogenic sources of ice-forming nuclei: a review[END_REF][START_REF] Després | Primary biological aerosol particles in the atmosphere: a review[END_REF][START_REF] Harrison | Not all feldspars are equal: a survey of ice nucleating Accepted Article Version postprint Comment citer ce document[END_REF]; however, those most efficient in catalyzing ice formation at warmer temperatures (> -15 °C) are microbial or biogenic in origin [START_REF] Murray | Ice nucleation by particles immersed in supercooled cloud droplets[END_REF], Morris et al 2013[START_REF] Fröhlich-Nowoisky | Ice nucleation activity in the widespread soil fungus Mortierella alpina[END_REF][START_REF] Weber | Polytrichum commune spores nucleate ice and associated microorganisms increase the temperature of ice nucleation activity onset[END_REF][START_REF] Fröhlich-Nowoisky | Bioaerosols in the Earth system: Climate, health, and ecosystem interactions[END_REF][START_REF] Hill | Impacts of bioaerosols on atmospheric ice nucleation processes[END_REF]). Thus, knowledge concerning microbes in precipitation is potentially relevant to understanding global precipitation itself.

Microbes are aerosolized from all terrestrial and aquatic ecosystems [START_REF] Wilson | A marine biogenic source of atmospheric ice-nucleating particles[END_REF][START_REF] Joung | Bioaerosol generation by raindrops on soil[END_REF], including the surfaces of multicellular organisms, e.g., plants [START_REF] Lindemann | Plants as sources of airborne bacteria, including ice nucleation-active bacteria[END_REF]Burrows et al. 2009a[START_REF] Després | Primary biological aerosol particles in the atmosphere: a review[END_REF], and fluxes of bioaerosols are often generated or affected by anthropogenic activities [START_REF] Taha | Bioaerosol releases from compost facilities: Evaluating passive and active source terms at a green waste facility for improved risk assessments[END_REF][START_REF] Suski | Agricultural harvesting emissions of ice-nucleating particles[END_REF]. Once aerosolized, some microbes are well suited for survival in atmospheric environments. For instance, dormant forms of sporeforming bacteria are commonly isolated from aerosols [START_REF] Brunet | Short-Scale Transport of Bioaerosols[END_REF]. Still other microbial taxa are metabolically active in fog [START_REF] Fuzzi | Fog droplets: An atmospheric source of secondary biological aerosol particles[END_REF]) and super-cooled cloud droplets [START_REF] Sattler | Bacterial growth in supercooled cloud droplets[END_REF], apparently using organic acids and alcohols as growth substrates [START_REF] Amato | An important oceanic source of micro-organisms for cloud water at the Puy de Dôme (France)[END_REF], Vaïtilingom et al. 2010). Wet deposition (washout) serves as the primary mechanism for removing bioaerosols from higher atmospheric layers [START_REF] Brunet | Short-Scale Transport of Bioaerosols[END_REF][START_REF] Hanlon | Microbial ice nucleators scavenged from the atmosphere during simulated rain events[END_REF]. Because aerosolized microbial INPs may trigger ice formation and be-scavenged-by precipitation, they are predicted to have much shorter atmospheric residence times than other microbes (Burrows et al. 2009b[START_REF] Stopelli | Ice nucleation active particles are efficiently lost from precipitating clouds[END_REF]. Thus, the ice nucleation-active (INA) phenotype allows broad dissemination of aerosolized microbes [START_REF] Sands | The association between bacteria and rain and possible resultant meteorological implications[END_REF][START_REF] Fitt | The role of rain in dispersal of pathogen inoculum[END_REF]) while limiting prolonged exposure to damaging stresses (e.g., desiccation, UV radiation) associated with atmospheric transport [START_REF] Attard | Effects of atmospheric conditions on ice nucleation activity of Pseudomonas[END_REF][START_REF] Hill | Impacts of bioaerosols on atmospheric ice nucleation processes[END_REF]. Despite this information, however, the community ecology of these processes is poorly understood. For instance, it is largely unknown if atmospheric spatiotemporal characteristics, including site locations and season, correspond to general trends in precipitation-borne microbial taxa or taxa assemblages (Amato et al. 2017, pg. 8, but see [START_REF] Cáliz | A long-term survey unveils strong seasonal patterns in the airborne microbiome coupled to general and regional atmospheric circulations[END_REF].

Accepted Article

Version postprint the lower stratosphere [START_REF] Griffin | Terrestrial microorganisms at an altitude of 20,000 m in Earth's atmosphere[END_REF][START_REF] Després | Primary biological aerosol particles in the atmosphere: a review[END_REF][START_REF] Bryan | Abundance and survival of microbial aerosols in the troposphere and stratosphere[END_REF].

In general, the concentration of atmospheric microbes decreases with altitude [START_REF] Bryan | Abundance and survival of microbial aerosols in the troposphere and stratosphere[END_REF]. Further, because of frictional and gravitational constraints, bioaerosols in the upper troposphere will generally be in the size range of bacteria (0.5 -1.5 µm dia.; [START_REF] Mandrioli | Evidence for long range transport biological and anthropogenic aerosol particles in the atmosphere[END_REF], with long residence times and the potential for global-scale dispersal [START_REF] Brunet | Short-Scale Transport of Bioaerosols[END_REF]. Larger biological particles, including hyphal fragments and fungal spores (2-40 µm dia.), will generally be limited to the PBL, or even the surface boundary layer (SBL; 10-100 m AGL) which at its bottom has the properties of a plane-mixing layer and is subject to high turbulence and small-scale intermittent eddies [START_REF] Raupach | Coherent eddies and turbulence in vegetation canopies: the mixing-layer analogy[END_REF]. Thus, larger microbes will generally have smaller dispersal distances (< 1 km) and residence times [START_REF] Brunet | Short-Scale Transport of Bioaerosols[END_REF][START_REF] Wilkinson | Modelling the effect of size on the aerial dispersal of microorganisms[END_REF]. As an exception, certain meteorological conditions may allow large dispersal distances (> 100 km) of microbes in the PBL (Aylor et al. 2006, Schmale and[START_REF] Schmale | Highways in the sky: scales of atmospheric transport of plant pathogens[END_REF]. Moreover, the strongest phenotypic correlate of aerosolized particle dispersal is aerodynamic diameter, for which physical size is an imprecise proxy [START_REF] Reponen | Aerodynamic versus physical size of spores: Measurement and implication for respiratory deposition[END_REF]. Small aerodynamic diameters can contribute to the relatively long dispersal distances of some physically large spores, including the urediospores of rust fungi (Morris et al. 2014b).

Accepted Article

Version postprint

Comment citer ce document : Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394

This article is protected by copyright. All rights reserved

A Metacommunity Framework for Microbes in Precipitation

A metacommunity defines a set of communities linked by the dispersal of potentially interacting species [START_REF] Leibold | The metacommunity concept: a framework for multi-scale community ecology[END_REF]. This perspective acknowledges the potential importance of multiscale processes in which localities, that hold local taxa assemblages, are connected to other local assemblages within a region. A number of meteorological dispersal models for aerosolized microbes have been previously applied, generally relying on Eulerian and Lagrangian vector field approaches from multivariate calculus [START_REF] Yamamoto | Particle-size distributions and seasonal diversity of allergenic and pathogenic fungi in outdoor air[END_REF][START_REF] Wilkinson | Modelling the effect of size on the aerial dispersal of microorganisms[END_REF][START_REF] Schmale | Highways in the sky: scales of atmospheric transport of plant pathogens[END_REF].

While neither acknowledging metacommunity theory nor explicitly considering precipitation-borne organisms, these models emphasize the dynamic, dispersal-driven constraints of aerosolized microbes, in which weather patterns serve as multi-scale conduits among local microbial assemblages.

We note that a current lack of metacommunity studies of bioaerosols can be linked partially to difficulties with the term community (sensu Ricklefs 1976[START_REF] Stroud | Is a community still a community? Reviewing definitions of key terms in community ecology[END_REF] in microbial systems -because interactions among co-occurring microbial taxa are often negligible or unknown [START_REF] Boon | Interactions in the microbiome: communities of organisms and communities of genes[END_REF]. This is particularly true for atmospheric microbes; whose densities are much lower than those in terrestrial and aquatic environments. Despite these issues, however, tools for the description and analysis of ecological communities and metacommunities have been extremely useful in clarifying processes and characteristics of microbial assemblages [START_REF] Prosser | The role of ecological theory in microbial ecology[END_REF][START_REF] Boon | Interactions in the microbiome: communities of organisms and communities of genes[END_REF].

Partitioning β-diversity. A primary use of metacommunity theory is the elucidation of mechanisms underlying the structure and composition of taxa at multiple spatiotemporal scales [START_REF] Leibold | The metacommunity concept: a framework for multi-scale community ecology[END_REF]. For this purpose, metrics that quantify compositional heterogeneity among local assemblages, i.e., β-diversity, are often used (e.g., [START_REF] Winegardner | Diatom diversity patterns over the past c. 150 years across the conterminous United States of America: Identifying mechanisms behind beta diversity[END_REF]. Total β-diversity, denoted β Total , can be attributed to two processes: turnover of assemblage taxa, β turn , and assemblage nestedness, β nest [START_REF] Baselga | Partitioning the turnover and nestedness components of beta diversity[END_REF]. β turn and β nest , in turn, result from antithetic mechanisms among assemblages: taxa replacement and taxa loss, respectively [START_REF] Baselga | Partitioning the turnover and nestedness components of beta diversity[END_REF]).

To illustrate these properties of β-diversity partitioning in precipitation-borne microbial assemblages, we consider three hypothetical scenarios in which twelve microbial taxa (e.g., families, genera, species, OTUs) distinguish the bioaerosol composition of three precipitation events (Fig. 1).

In scenario A, the assemblages of precipitation events are completely nested, with lower richness

Accepted Article

This article is protected by copyright. All rights reserved events containing compositional subsets of richer events. Because all compositional heterogeneity among events is due to nestedness, β turn = 0, and β nest = β Total . In scenario B, all three events contain six microbial taxa, and (like scenario A) events have three species in common. However, at each event, the three remaining taxa are replaced by three other taxa and these replacements are distinct for each event. Because all compositional heterogeneity among events is due to turnover, β nest = 0, and β turn = β Total. In scenario C, both nestedness and turnover occur, and both contribute to β Total .

A graphical metacommunity heuristic. with smaller aerodynamic diameters, including bacterial cells, can disperse to higher altitudes [START_REF] Deleon-Rodriguez | Microbiome of the upper troposphere: species composition and prevalence, effects of tropical storms, and atmospheric implications[END_REF][START_REF] Bryan | Abundance and survival of microbial aerosols in the troposphere and stratosphere[END_REF], where they generally occur at lower concentrations (Burrows et al. 2009a[START_REF] Gandolfi | Unravelling the bacterial diversity in the atmosphere[END_REF], Griffin et al. 2017, pg. 155), and travel further horizontally (long narrow red arrows). Because of higher velocity air currents, the composition of high altitude bioparticles will be spatially macroscale-driven. Thus, changes in local microbial assemblages should occur over large distances (gradual color changes in hatching).

Further, because of repeated washout of microbes as storms move in time and space, coupled with isolation from surface environments -preventing subsidization of new taxa, β nest should be relatively high among high altitude assemblages (Fig. 1, scenario A), and assemblages will be poorly associated with particular precipitation sites. In contrast, lower atmospheric altitudes will contain higher concentrations of bioparticles, including relatively large microorganisms (e.g., fungal spores and pollen), with poorer dispersal, and stronger ties to local (mesoscale) conditions on the terrestrial surface. Within the PBL, relatively large microbial fluxes should occur over small horizontal distances (short wide green arrows) due to the higher depositional velocities of larger microorganisms [START_REF] Raupach | Coherent eddies and turbulence in vegetation canopies: the mixing-layer analogy[END_REF] and atmospheric mixing, causing strong spatial differentiation among local aerosolized assemblages (dramatic color changes in hatching). In this case, β turn should be relatively high among local assemblages, as spatiotemporal changes on the earth's surface drive corresponding changes to microbial assemblages in the adjacent atmosphere (Fig 1,scenario B). Wet deposition scrubs/scavenges microbes from overlaying atmospheric layers [START_REF] Hanlon | Microbial ice nucleators scavenged from the atmosphere during simulated rain events[END_REF], causing

Accepted Article

Bacterial 16S and fungal 18S rRNA gene amplicon libraries were sequenced from precipitation samples collected during all four seasons at LSU, and three seasons at ISU (spring, summer and winter), and VT (summer, winter and autumn) (Table 1). Relatively rare heavyprecipitation events contributed to the smaller number of samples at ISU. The 30 precipitation events sampled for 18S rRNA genes constituted a subset of the 71 events sampled for 16S except for one case (LSU-spring) in which a precipitation event was sampled for 18S but not 16S rRNA genes (Table 1). The subset of 18S analyses was due to the depletion of sample materials from the 16S rRNA gene amplification procedures, which generally preceded 18S rRNA gene sequencing. In 29 precipitation events both 18S and 16S rRNA genes were successfully amplified and sequenced.

For each sampling event, precipitation was collected in ~120 L galvanized collection cans that were lined with clean, sterile 94 x 122 cm polypropylene bags (Fisher Scientific, Pittsburgh, PA;[START_REF] Failor | Ice nucleation active bacteria in precipitation are unexpectedly diverse and nucleate ice by employing different mechanisms[END_REF]. A minimum volume of 3L was amalgamated from collection cans and sample collection was conducted using a total of eleven cans, ten for samples and one as a procedural control

Accepted Article

Version postprint

This article is protected by copyright. All rights reserved (see Appendix S1,Fig S1). Immediately following each precipitation event, collected material was transferred to sterile 9L carboys (Appendix S1,Fig S2) or left in the collection bags and stored at 4 °C preceding concentrating by filtration, which was always completed within 48 h of precipitation collection. Precipitation was filtered (1-3 L) onto 0.2 µm, 47 mm Supor PES membrane filters or Millipore polycarbonate filters. Filters were placed into sterile 15 mL or 50 mL falcon tubes and frozen preceding further processing. All samples from Louisiana and Virginia were shipped overnight on dry ice to ISU for consolidation with ISU samples.

Sampling procedures were designed to ensure that 16S and 18S rRNA gene sequences were not obtained from contaminated sampling equipment, sampling processes, or downstream molecular procedures. Procedural controls were created by lining one galvanized steel can with a sterile polypropylene bag (just as was done for precipitation sample collection) and covering the can with an ethanol-cleaned lid (Appendix S1,Fig S1). One L of autoclaved nanopure water was poured into the control bag and this water was filtered in the same manner as the precipitation samples. All filters containing procedural controls were put through the same DNA extraction and gene amplification procedures to assess spurious amplification of 16S and 18S rRNA genes due to procedural and/or reagent contamination.

DNA extraction from precipitation samples

Microbial wet lab procedures, including choices of primers, followed guidelines of the Earth Microbiome Project [START_REF] Thompson | A communal catalogue reveals Earth's multiscale microbial diversity[END_REF]. In preparation for DNA extraction, filters were cut using ethanol cleaned and flamed scissors. All filter cutting was performed in a laminar flow hood and filters were only handled with sterile tweezers. Each filter was cut into the smallest possible pieces (< 5mm in length) just above a bead beating tube from the FastDNA SPIN Kit for Soil (MP Biomedicals, USA). Once in the bead beating tube, the manufacturer's protocol for the FastDNA SPIN Kit for Soil was followed, starting with the addition of MT and phosphate buffer solutions. The manufacturer's protocol was followed with the exception of the following steps: in step 4, a mini bead beater was used to homogenize extracts for 70 seconds, step 5 was performed for 8 minutes, a 15 mL tube was always used for step 7, and 100 µL of DNase/pyrogen-free water was used to elute DNA in step 16. The resulting extracts were further purified with the aid of a silica spin filter, following steps

Accepted Article

This article is protected by copyright. All rights reserved 14-22 of the protocol for the MoBio Power Soil Kit (MoBio Laboratories, Carlsbad, CA). In these steps, 160 µL of solution C4 was added to the DNA extraction. The manufacturer's protocol was followed for other steps, with the exception of adding 25 µL of solution C6 and incubating for 1 minute at room temperature prior to centrifugation. All DNA extracts were stored at -20 °C until PCR amplification of the 16S and 18S rRNA genes was attempted.

Bacterial 16S rRNA gene fragment amplification

The V4 regions of bacterial 16S rRNA gene fragments were PCR-amplified in triplicate from each DNA extract for sequencing. PCR amplification was carried out using barcoded primers following the methods of [START_REF] Caporaso | Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms[END_REF]. The 25 µL PCR mix contained the following: 1X 5-PRIME HotMaster Mix, 0.2 µM 515F, 0.2 µM 806R and nuclease free water (13 µL). Thermal-cycling was carried out using 1-6 µL of DNA template in an Eppendorf PRO S Master Cycler under the following conditions: initial denaturation at 94 °C for 3 minutes, 32 cycles of denaturation at 94 °C for 1 minute, annealing at 50 °C for 1 minute and extension at 72 °C for 1:45, followed by a final extension at 72 °C for 10 minutes. We resolved electrophoretically separated PCR products on a 1% TAE gel stained with ethidium bromide to confirm successful amplifications in each reaction. Triplicate reactions were pooled before purification and concentration using the Qiagen MinElute PCR Cleanup kit (Qiagen, Valencia, California, USA). A 35% guanidine-HCl wash step was performed to ensure that large primer-dimers were removed from the samples. Amplicons were stored at -20 °C preceding sequencing at Idaho State University's Molecular Research Core Facility (Pocatello, ID, USA) on the Illumina MiSeq Platform using the V2 500 bp kit.

Fungal 18S rRNA gene fragment amplification

Fungal 18S rRNA gene fragments from the V9 region were PCR-amplified in triplicate from each DNA extract. Steps for 16S amplification were also followed for 18S amplification with two exceptions. First, the 25 µL PCR mix contained 1X 5-PRIME HotMaster Mix, 0.2 µM 1391F, 0.2 µM EukBr and nuclease free water (13 µL). Second, thermal cycling used the following conditions: initial denaturation at 94 °C for 3 minutes, 32 cycles of denaturation at 94 °C for 0:45, annealing at 57 °C for 1 minute and extension at 72 °C for 1:30, followed by a final extension at 72 °C for 10

Accepted Article

This article is protected by copyright. All rights reserved minutes.

Sequence data analysis

Sequence data were analyzed using the mothur (ver. 1.42.3) software package [START_REF] Schloss | Introducing mothur: open-source, platformindependent, community-supported software for describing and comparing microbial communities[END_REF].

Contigs were assembled and parsed on the basis of unique barcodes attached to the 806R primer [START_REF] Caporaso | Ultra-high-throughput microbial community analysis on the Illumina HiSeq and MiSeq platforms[END_REF]. Sequences that did not contain exact matches to the primer and barcodes utilized in the PCR amplification were discarded. Sequences were filtered for quality with a 50-base sliding window. Sequences not obtaining an average quality score of 25, containing ambiguous bases, homopolymers (> 7 bases), having lengths > 259 bases, or deemed chimeric using the UCHIME algorithm [START_REF] Edgar | UCHIME improves sensitivity and speed of chimera detection[END_REF] were eliminated.

All 16S rRNA gene sequences were aligned to the SILVA bacterial 16S rRNA gene reference (release 132) [START_REF] Pruesse | SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB[END_REF]) truncated to contain the amplified V4 hypervariable region.

Sequences were taxonomically classified using the classify.seqs command with arguments set to default values, with the exception of the minimum consensus confidence value being set at 80%. All 16S rRNA gene sequences that were classified as eukaryotic, archaea, mitochondria, or chloroplast were removed from analyzed datasets. Further, pseudoreplicate 16S rRNA gene samples with fewer than 40,000 total reads, following removal of non-informative OTUs, were eliminated from downstream analyses, and thus were not averaged with other pseudoreplicates to obtain a single independent 16S observation per storm (see Appendix S1: Table S2). The 16S rRNA gene sequence library size was normalized by randomly sampling a total of 44,621 sequences from each library, which was the number of sequences in the smallest library.

The sequences of 18S rRNA gene amplicons for each sample were aligned in mothur to the SILVA 18S rRNA gene reference (release 132), truncated to contain the V9 hypervariable region, with minimum consensus confidence set at 80%. Non-fungal sequences were eliminated from the 18S dataset. Analogous to methods for bacteria, pseudoreplicate fungal samples with fewer than 20,000 reads were eliminated. Average read totals for events were approximately 40% lower for 18S than for 16S rRNA gene sequence libraries. The 18S sequence library size was normalized by randomly sampling a total of 28,586 sequences from each library; the number of sequences in the smallest library. This article is protected by copyright. All rights reserved Both 16S and 18S rRNA gene sequences were clustered into operational taxonomic units (OTUs) based on a sequence dissimilarity of ≤ 0.03, using the Needleman-Wunsch algorithm (Kunin et al. 2010, Smith and[START_REF] Smith | Sequence depth, not PCR replication, improves ecological inference from next generation DNA sequencing[END_REF]. Following the quality control procedures described above, the 16S rRNA gene dataset contained 34,107 OTUs describing 71 independent precipitation event samples, and the 18S rRNA gene dataset contained 16,597 OTUs describing 30 independent precipitation event samples. Both 16S and 18S rRNA gene sequence libraries were deposited into MG-RAST metagenomics RAST (Rapid Annotation using Subsystem Technology) server (https://metagenomics.anl.gov). Accession numbers for sequences are given in Appendix S1 ( § S1.6.1). Altitudes of 100, 1000, 5000, 6000, 8000 and 9000 mAGL were chosen for trajectory examination to represent vertical sections of ecological importance in our metacommunity framework (Fig. 2). The 9000 m AGL upper altitudinal bound was chosen due to its proximity to the homogeneous freezing layer (~ -38°C), where ice nucleation by heterogeneous mechanisms is no longer necessary for the freezing of water molecules.

Accepted Article

Collection of Meteorological Data

Classification of convective vs. stratiform precipitation was based on the methods of [START_REF] Biggerstaff | An improved scheme for convective/stratiform echo classification using radar reflectivity[END_REF], and [START_REF] Anagnostou | A convective/stratiform precipitation classification algorithm for volume scanning weather radar observations[END_REF], in addition to the inspection of tropospheric stability indices from National Weather Service soundings (NWS 2018)

Accepted Article

This article is protected by copyright. All rights reserved https://www.weather.gov/lmk/indices), and infrared and visible satellite imagery obtained from GOES-East (NOAA GOES-East 2018). Stratiform precipitation is defined specifically as precipitation collected from stratus and nimbostratus-like cloud systems independent of trailing stratiform regions from convective storm formations [START_REF] Biggerstaff | An improved scheme for convective/stratiform echo classification using radar reflectivity[END_REF]Listemaa 2000, Houze 2014).

Precipitation events were also assigned to six morphological classes: Air mass, Mesoscale Convective System (MCS), Multicell, Nimbostratus, and Squall.

Precipitation event trajectories were classified into five storm origin categories based on a conventional scheme for North American air masses (AMS 2018). Seventy-two-hour trajectories, which originated over marine surfaces, were classified as maritime (m), whereas those that originated over terrestrial surfaces were classified as continental (c). Likewise, those events that originated at latitudes above the subtropical regions of the Northern Hemisphere (> 23.5°N) were classified as Polar (P), and those originating in the subtropical and tropical regions of the Northern Hemisphere (< 23.5°N) were classified as tropical (T). Thus, six possible air mass combination classes were considered: Pacific Maritime Polar (mP1), Atlantic Maritime Polar (mP2), Pacific Maritime Tropical (mT1), Atlantic Maritime Tropical (mT2), Continental Polar (cP), and Continental Tropical (cT).

Statistical Analysis

General constraints. Meteorological data, and univariate summaries of microbial data, were largely analyzed using conventional univariate statistical approaches including general linear models (e.g. ANOVA, ANCOVA) or robust analogues, and tabular methods. Microbial compositional summaries and hypothesis tests required multivariate approaches. To address missing cells, a single ISU autumn and VT summer observation were imputed for ANOVA-type analyses considering site and season interactions, using the random forest algorithm [START_REF] Stekhoven | MissForest-non-parametric missing value imputation for mixed-type data[END_REF]. Imputed points were assumed to be missing completely at random (MCAR). Lack of balance in site and season levels prompted the use of type II (marginal) sums of squares in ANOVA-type analyses, or the emulation of type II sums of squares by changing the order of predictors in linear models and reporting summaries from the last specified predictor. Family-wise type I error (FWE) was controlled, in pairwise comparisons following omnibus ANOVA-type analyses, using [START_REF] Tukey | Comparing individual means in the analysis of variance[END_REF] HSD method. In multiple pairwise tests subsequent to non-ANOVAs, [START_REF] Holm | A simple sequentially rejective multiple test procedure[END_REF] This article is protected by copyright. All rights reserved was used to control FWE (see Aho 2013, Ch. 10). For indicator species analyses (see below) large families of tests for individual taxa were required. In this case, p-values were adjusted for false discovery rate (FDR) within taxonomic hierarchies (i.e., phylum, class, order, family, genus) using the method of [START_REF] Benjamini | Controlling the false discovery rate: a practical and powerful approach to multiple testing[END_REF]. Control of FDR is more appropriate than control of FWE in exploratory applications involving a large number of null hypothesis tests (cf. Hulsman et al. 2013).

Throughout the paper we use the significance level α = 0.05.

Diversity. We measured bacterial and fungal OTU α-diversity in precipitation events using the Shannon-Wiener index [START_REF] Shannon | A mathematical theory of communication[END_REF]. We partitioned total β-diversity (β Total ) into nestedness (β nest ) and spatial turnover (β turn ) components using the incidence-based pairwise measures given in Baselega (2010), i.e., the means of the β sor , β nes , and β sim dissimilarity matrices, respectively. In the For our study, OTUs were defined using rRNA gene sequences of the V4 and V9 regions for bacteria and fungi, respectively. However, inherent physical differences and putative differences in evolutionary history exist for these gene regions. These discrepancies limit comparability of bacterial and fungal diversity using rRNA gene sequences. For example, in a simpler application involving only eukaryotes, Tragin et al. (2017) found that V4 and V9 regions of 18S rRNA genes in phytoplankton samples had similar levels of α-diversity, but the V9 region identified 20% more OTUs. To address these issues, we constrained fungal and bacterial comparisons in four ways. First, we formally compared bacteria and fungi diversity using only β-partitioning measures (i.e., trends in overall heterogeneity of assemblages due to OTU turnover and nestedness), as reflected in our

Accepted Article

Comment citer ce document : This article is protected by copyright. All rights reserved hypotheses. Second, to address richness biases in the 16S V4 regions and 18S V9 rRNA gene regions, we chose β-partitioning metrics shown to be independent of the richness of OTU assemblages [START_REF] Baselga | Comparing methods to separate components of beta diversity[END_REF]. Third, we corrected for inherent differences in the heterogeneity of 16S and 18S rRNA genes by dividing fungal and bacterial β nest and β turn measures by their respective measures of β Total . Thus, in bacterial/fungal comparisons, both β nest and β turn were scaled to be in [0, 1], and represented the proportion of variability in total β-diversity due to nesting and turnover, respectively. Fourth, formal comparisons of fungi and bacteria β-partitioning measures were based only on the 29 rain events for which both 16S and 18S rRNA gene sequences were obtained. The complete set of 72 rain events were a mix of paired (spatiotemporally dependent) and non-paired events (Table 1). Bootstrapping significance tests were constrained so that bacterial and fungal data were paired by precipitation event preceding sampling with replacement. Corresponding p-values were based on one-tailed (reflecting hypothesized directions) percentile confidence intervals [START_REF] Aho | Foundational and applied statistics for biologists using R[END_REF] for the true difference in β-partitioned values.

PERMANOVA. Permutational Multivariate Analysis of Variance (PERMANOVA; [START_REF] Anderson | A new method for non-parametric multivariate analysis of variance[END_REF] was used to consider the effect of macroscale (large spatial scale and seasonal), and mesoscale (sitespecific) drivers of bacterial and fungal precipitation-borne assemblages (Fig. 2). An unresolved problem in next-generation sequencing (NGS) of PCR rRNA gene amplicons is the occurrence of microorganisms whose genomes contain multiple copies of the rRNA gene, causing particular OTUs to have inflated read abundances (Herrera 2009). To address this issue, we conducted PERMANOVAs using both raw read abundances of OTUs and binary-transformed (presence/absence) OTU data. PERMANOVA p-values were based on 10000 permutations.

To decrease the propensity for single precipitation events or OTUs to influence multivariate analyses, we applied a Wisconsin double transformation [START_REF] Bray | An ordination of the upland forest communities of southern Wisconsin[END_REF] This article is protected by copyright. All rights reserved dissimilarity effectively represents and distinguishes multivariate assemblages given sparse (few nonzero entries) datasets which may contain many samples with nothing in common [START_REF] Aho | Using geometric and non-geometric internal evaluators to compare eight vegetation classification methods[END_REF].

Ordination. The relation of precipitation events in OTU space was depicted using Nonmetric Multidimensional Scaling ordinations (NMDS; [START_REF] Kruskal | Multidimensional scaling[END_REF]. These analyses were based on non-imputed datasets with all recorded OTUs. The final NMDS projections were based on the best (lowest stress) solutions resulting from twenty randomized starting configurations. To enhance interpretability, NMDS axes were centered and rotated to principal components. Vector (quantitative data) and factor (categorical data) fitting [START_REF] Oksanen | vegan: community ecology package. R package version 2[END_REF]) was used to examine correlations between raw meteorological data and NMDS scatter. Permutational approaches were used to test null hypotheses of no association in these procedures using 1000 iterations.

Indicator Species Analysis. Indicator species analysis (ISA; Dufrene and Legendre 1997) using nonimputed data was used to test for association of taxa designations of bacterial and fungal OTUs to season and site levels. All unclassified taxa (i.e., "unclassified", "uncultured", "unknown", "Incertae Sedis", "unnamed") were eliminated from ISA procedures. Specific taxa were required to occur at two or more precipitation events for inclusion in analyses. P-values for tests of the null hypothesis that associations between species and categorical levels were no better than random were based on permutation procedures using 1000 iterations.

Software. All analyses were performed within the R statistical environment (R Core Team 2017) with heavy reliance on the packages asbio [START_REF] Aho | asbio: A collection of statistical tools for biologists[END_REF]) and car (Fox and Weisburg 2011) for implementation of conventional general linear models; vegan [START_REF] Oksanen | vegan: community ecology package. R package version 2[END_REF] for most compositional analyses of assemblages including PERMANOVA, ordinations, and matrix temperature calculations; betapart (Baselga and Orme 2012) for beta-diversity partitioning analyses; labdsv [START_REF] Roberts | labdsv: ordination and multivariate analysis for ecology[END_REF] for ISA; and missForest [START_REF] Stekhoven | MissForest: non-parametric missing value imputation for mixed-type data[END_REF] for imputation.

Results

Meteorological Differences among Sites and Seasons

Accepted Article

Comment citer ce document : This article is protected by copyright. All rights reserved As expected (Fig. 2), higher altitude atmospheric layers traveled further in 72 h back-trajectories (slope ≈ 0.01 m m -1 hr -1 , F 1,420 = 391.3, p < 2.2 × 10 -16 ). Significant differences in distance traveled occurred by site when holding season and elevation constant (F 2,420 = 11.2, p = 0.0002), with Idaho storms traveling the furthest and LSU storms arising more locally. These differences were increasingly apparent at higher altitudes (interaction effect: F 2,420 = p = 0.0007). Strong seasonal effects were also evident (F 3,420 = 41.58, p < 2.2 × 10 -16 ), with shorter summer trajectories and longer trajectories in the winter. No trajectory distance interactions occurred between site and season.

Fisher's exact tests identified significant differences (p ≤ 0.03) in cloud air mass origins (i.e., cP, cT, mP1, mP2, mT1, mT2) between sites at all altitudes but one (8000m) when controlling for FWE using [START_REF] Holm | A simple sequentially rejective multiple test procedure[END_REF] procedure. These differences were largely due to the maritime/tropical origin of most LSU and VT precipitation events, compared to Idaho events, which were mostly polar/continental. Idaho maritime events were never Atlantic, whereas VT and LSU maritime events were frequently Atlantic in origin (Fig. 3). Significantly (p < 0.025) distinct seasonal characteristic for air mass origin were evident at all altitudes above 5000 after controlling for FWE. This finding was largely due to precipitation events from LSU (the best-sampled site) being mostly maritime, except for summer storms which were primarily continental (Fig. 3). Weather morphology classes (Airmass, MCS, Multicell, Nimbostratus, Squall) differed significantly by season (p = 0.008), with nimbostratus events occurring primarily in the winter and fall, but not site (p = 0.389; Fig. 3). Similar patterns were apparent for stratiform/convective classes with significant differences occurring with season (p = 0.003), due to the absence of stratiform storms during summer months, but a lack of strong differences with site (p = 0.102).

Microbial Diversity Varies with Season, and Particularly Sites

Site, season, and season by site interactions were not significantly associated with richness and Shannon diversity for either the bacterial or fungal datasets (Table 2). LSU had significantly higher bacterial β turn and β Total than ISU, and LSU had significantly higher β turn and β Total than ISU in the fungal dataset (Table 2; see Appendix S1: Tables S86 andS87 for bootstrap confidence intervals).

Thus, for both 16S and 18S rRNA gene data, LSU events had higher levels of compositional heterogeneity, particularly compared to ISU, and this heterogeneity was due to relatively high taxa

Accepted Article

This article is protected by copyright. All rights reserved turnover in the bacterial dataset and relatively high taxa nestedness in the fungal dataset.

Fungal and Bacterial β-Diversities are Distinct--β turn higher for Fungi, β nest higher for Bacteria.

Fungal and bacterial β-diversity partitioning components differed. Spatiotemporal turnover was higher for fungal assemblages in precipitation, causing β turn /β Total outcomes for fungi to be significantly higher than for bacteria (p < 0.0001; percentile bootstrap confidence intervals; R = 1000; Fig. 4). Assemblage nesting was higher for bacteria, causing bacterial β nest /β Total values to be significantly higher than for fungi (p < 0.0001; Fig. 4).

Bacterial Composition Varies Subequally with Site and Season. Fungal Composition Varies with Site

PERMANOVA analyses of raw and binary transformed data were in strong agreement (Appendix S1: Tables S24, S25, S56, S57). Thus, we present only analyses of raw reads here (Table 3). For 16S rRNA gene library reads, OTU compositional differences had identical levels of statistical significance for site (F 2,61 = 6.32, p = 0.0001) and season (F 3,61 = 2.7, p = 0.0001). For 18S rRNA gene reads, however, much stronger compositional differences were apparent among sites (F 2,20 = 3.1, p = 0.0002) compared to seasons (F 3,20 = 1.46, p = 0.06). Interactions between site and season were absent in both bacterial and fungal datasets, indicating that seasonal assemblage fluctuations were somewhat consistent across sites (Table 3). In pairwise comparisons of bacterial OTU composition, all sites and seasons were significantly different after controlling FWE (Table 3). In pairwise comparisons of fungal OTU composition, all sites were significantly different (Table 3). Seasons were not compared in the fungal dataset due to a lack of omnibus significance. Notwithstanding the effects of smaller sample sizes and inherent inconsistencies in 16S and 18S rRNA gene heterogeneity, our PERMANOVA results suggest that compositional variations were smaller for fungal assemblages, compared to bacterial assemblages, particularly with respect to seasons.

Ordinations Demonstrate Spatiotemporal Variation in Assemblages. Patterns of variation in microbial metacommunity composition detected in PERMANOVA are reflected in NMDS ordinations of 16S and 18S rRNA gene data with respect to both site (Fig. 5) and season (Fig. 6). Also evident, with vector and factor fitting, were the longer distances traveled by Idaho events over 72 hour

Accepted Article

This article is protected by copyright. All rights reserved back-trajectories, and the continental and Atlantic sources of LSU and VT storms, compared to the Pacific maritime and polar sources of Idaho storms (Fig. 5a,6a). Correlations of the NMDS projections to meteorological variables are given in Table 4. Significant 16S rRNA gene correlates included season (p = 0.025; R 2 = 0.11), but not site (p = 0.154, R 2 = 0.13). Correlations of the projection to storm travel distances strengthened as elevation increased. Specifically, for elevations of 100, 1000, 5000, 8000, and 9000m, R 2 values were: 0.01, 0.06, 0.15, 0.23, 0.22, 0.19, respectively, and p-values were: 0.64, 0.10, 0.02, 0.001, <0.001, and 0.001, respectively (Table 4). These results demonstrate that strong relationships exist between meteorological measures and bacterial OTU composition, particularly storm characteristics at higher altitudes.

In the 18S rRNA gene ordination (Figs. 5b,6b), fewer significant relationships occurred between abiotic variables and the projection of precipitation events; and the ordination solution was inferentially weaker (stress = 0.159), compared to the bacterial ordination (stress = 0.08). Despite these differences, we present two-dimensional solutions for both datasets to facilitate their comparison. The fungal projection was strongly associated with site (p = 0.002; R 2 = 0.28; Table 4) but not season (p = 0.13; R 2 = 0.17; Table 4), reversing trends observed for bacteria (Table 4). Other significant correlates of the fungal ordination were meteorological characteristics of moderate altitudes, i.e., air mass categories at 1000m (p = 0.028; r = 0.25), and 6000m (p = 0.003; r = 0.35), and precipitation type, i.e., rain or ice (p = 0.01, R 2 = 0.15). Thus, the NMDS ordination of bacterial OTUs was more strongly correlated with variation in large spatiotemporal scale measures, e.g., season and higher altitude back trajectory distance, whereas the NMDS ordination of fungal OTUs was more strongly correlated with smaller scale spatiotemporal measures, including particular site locations, and moderate altitude air mass categories (Table 4).

Spatiotemporal Dynamics of Taxonomic Composition

The average consensus confidence level for alignment of 16S rRNA gene sequences to taxa was greater than 98.3% for all taxonomic levels, and equaled 100% for levels above order (Appendix S1: Table S28). For 18S rRNA gene sequences, consensus confidence was greater than 98.9% for all taxonomic levels, and equaled 100% for levels above order. For 16S rRNA gene data, unclassified taxa comprised 4%, 5%, 11%, 13%, and 31% of reads across samples for phylum, class, order, family

Accepted Article

This article is protected by copyright. All rights reserved and genus level designations, respectively. Taxonomic designations based on 18S rRNA gene sequences included a large proportion of unclassifiable OTUs for levels below order. Specifically, unclassified taxa constituted 0%, 9%, 32%, 60%, and 71% of reads for phylum, class, order, family and genus level designations, respectively. As a result, we do not present ISA results for fungal taxonomic levels below order. Complete results for fungal families and genera taxa, however, are given in Appendix S1. A large number of bacterial and fungal taxa were statistically significant indicators of sites and/or seasons (Table 5). Thus, only the highest average abundance bacterial and fungal taxa are summarized graphically (Figs 7-10), and only these taxa are considered in the text below. Complete ISA results and taxonomic summaries are given in Appendix S1.

Bacteria. Six named phyla and one unclassified group comprised 98.4% of the total number of 16S rRNA gene sequence reads. In particular, Proteobacteria, Cyanobacteria, Bacteroidetes, Actinobacteria, Firmicutes, Acidobacteria and "unclassified" made up 58%, 14%, 8.1%, 7.6%, 6.1%, 1.6% and 3.6% of read totals, respectively (Figs. 7,8). Proteobacterial subdivisions were dominant in lower taxonomic hierarchies with Gammaproteobacteria, Alphaproteobacteria and Deltaproteobateria making up 26.1%, 25.9% and 4.9% of class read totals, respectively, and, Pseudomonadales (Gammaproteobacteria), Betaproteobacteriales (Betaproteobacteria), Rhizobiales (Alphaproteobacteria), and Sphingomonadales (Alphaproteobacteria), and Burkholderiales (Betaproteobacteria) making up 9.3%, 8.9%, 8.7%, and 8.5%, respectively, of order read totals.

A large percentage of bacterial taxa were significant study site indicators (18% of phyla, 20%, of classes, 18% of orders, 22% of families, and 21% of genera), although much lower percentages of significant taxa remained after controlling for FDR (8% of phyla, 0%, of classes, 3% of orders, 6% of families, and 5% of genera) (Table 5). With respect to bacterial indicators of seasons, 18% of phyla, 21%, of classes, 17% of orders, 17% of families, and 14% of genera were significant indicators, and 0% of phyla, 1%, of classes, 2% of orders, 3% of families, and 2% of genera were statistically significant after controlling for FDR (Table 5).

ISU had the highest number of significant bacterial site indicators across all taxonomic levels (Appendix S1: Table S45). Ignoring FDR, ISU significant indicators included Actinobacteria (phylum), Actinobacteria (class), Cytophagales and Micrococcales (order), Burkholderiaceae,

Accepted Article

This article is protected by copyright. All rights reserved Sphingobacteriaceae, and Microbacteriaceae (family) and Hymenobacter (genus) (Fig. 7). FDRadjusted significant ISU indicators included Massilia (genus). Ignoring FDR, significant LSU indicators include Archangiaceae and Chroococcidiopsaceae (family), and Methylobacterium (genus).

FDR-adjusted significant LSU indicators included Nostocales (order) and Cloacibacterium (genus).

Ignoring FDR, significant VT indicators included Acidobacteriia (class) and Acidobacteriales (family). Significant FDR-adjusted VT indicators included X1174.901.12 (genus).

The season with the largest number of significant bacterial indicators was winter, followed by summer, spring and autumn (Appendix S1: Table S39). Indicators of autumn included the Pseudomonadales (order), and its subtaxa Pseudomonadaceae (family), and Pseudomonas (genus) (Fig. 8). Indicators of spring included Alphaproteobacteria (class), and its subtaxa: Sphingomonadales (order) Sphingomonadaceae (family) and Sphingomonas (genus). The latter three remained statistically significant after controlling FDR. Indicators of summer included three taxa that remained significant after controlling FDR: Deltaproteobacteria (class), and its subtaxa Myxococcales (order), Cystobacteraceae (family), and Melittangium (genus), along with the genus Cloacibacterium.

Indicators of winter included two taxa that remained significant after controlling FDR: Hymenobacteraceae (Family) and Hymenobacter (genus), along with Clostridia, Thermoleophilia, and Verrucomicrobiae (class), Cytophagales (order), and Pedobacter (genus).

Fungi. A relatively large number of fungal taxa (20% of phyla, 15%, of classes, and 29% of orders) were significant site indicators, and 6% of orders remained statistically significant site indicators after controlling for FDR (Table 5). As with bacterial results, ISU had the highest number of indicative taxa, with 15 significant ISA orders, before controlling for FDR, compared to LSU (five orders) and VT (seven orders) (Appendix S1: Table S76). Fungal taxa were poor predictors of seasons. No phyla, only 12% of classes, and 4% of orders were statistically significant ISA season indicators (Table 5, Fig 10), and only 2% of genera were significant seasonal indicators after controlling for FDR (Table 5). The most distinct season levels were summer and spring, which had eight and seven significant order indicators, respectively, before controlling for FDR (Appendix S1: Table S71).

Fungal indicators for ISU included Cystobasidiomycetes, Arthoniomycetes, and Tremellomycetes (class) and Capnodiales and Filobasidiales (order). The latter was significant after

Accepted Article

This article is protected by copyright. All rights reserved controlling FDR (Fig. 9). LSU indicators included Basidiomycota (phylum). VT indicators included Auriculariales (order). Ignoring FDR, spring was indicated by Arthoniomycetes and Taphrinomycetes (class), and summer was indicated by Ustilaginomycetes (class) and the order Ustilaginales (Fig. 10).

Discussion

Bioaerosols constitute a vital component of the biosphere. However, fundamental ecological components of these systems, including microbial composition of precipitation, and the spatiotemporal dynamics and drivers of these assemblages, remain poorly understood [START_REF] Amato | Active microorganisms thrive among extremely diverse communities in cloud water[END_REF]. The limited number of existing studies of precipitation-borne microbes have: 1) generally not considered taxa other than bacteria, or distinguished INPs, 2) concentrated on nearby or identical sites, and/or few points in time, and 3) have not effectively utilized existing analytical and theoretical frameworks for community and metacommunity ecology [START_REF] Bigg | The results of large-scale measurements of natural ice nuclei[END_REF][START_REF] Bigg | Ice nucleus concentrations in remote areas[END_REF][START_REF] Constantinidou | Atmospheric dispersal of ice nucleation-active bacteria: The role of rain[END_REF][START_REF] Amato | An important oceanic source of micro-organisms for cloud water at the Puy de Dôme (France)[END_REF][START_REF] Christner | Geographic, seasonal, and precipitation chemistry influence on the abundance and activity of biological ice nucleators in rain and snow[END_REF][START_REF] Kourtev | Atmospheric cloud water contains a diverse bacterial community[END_REF], Huffman et al.2013[START_REF] Joly | Ice nucleation activity of bacteria isolated from cloud water[END_REF], Tobo et al. 2013, Šantl-Temkiv et al. 2015[START_REF] Demott | Sea spray aerosol as a unique source of ice nucleating particles[END_REF]).

To address this deficiency, we conducted comparative analyses using what is, to date, the most extensive dataset of large-scale spatiotemporal trends in microbial assemblages distributed in precipitation. These data were considered using a metacommunity heuristic (Fig 2) that predicts bioaerosol dispersal patterns are largely a function of altitude, which in turn serves as an environmental filter based on organism size. The metacommunity framework prompted the evaluation of two hypotheses.

Hypothesis 1: Precipitation-borne bacterial assemblages are correlated with macroscale factors

Bacteria assemblages in precipitation were more strongly associated with macroscale spatiotemporal processes than fungi (H 1 ). This trend was demonstrated through both the relatively strong associations of bacterial assemblages to non-local high-altitude meteorological variables (P 1a ), and the higher levels of nestedness in bacterial assemblages, compared to fungi (P 1b ).

Prediction 1a. Assemblages of bacterial OTUs were strongly associated with spatially distant storm characteristics, compared to fungal OTUs. For instance, strong correlations occurred between, high-

Accepted Article

This article is protected by copyright. All rights reserved elevation air trajectory distances and the bacterial OTU NMDS ordination, although this effect was absent in the fungal projection (Table 4).

Strong bacterial macroscale associations were also evident from the perspective of large-scale temporal (seasonal) drivers. Season explained more variation in the bacterial OTU ordination than site in several multivariate analyses, and this trend was reversed for fungal OTUs (Figs. 4,5; Tables 3,4). This suggests that macroscale seasonal forces are more important than short-interval temporal subsidies in determining the composition of precipitation-borne bacterial assemblages. This finding was reiterated in taxonomic analyses, as a large number of significant bacterial indicators occurred for season compared to fungi (Table 5, Figs. 7 and9). The distinctiveness of ISU bacterial OTU composition (Table 3,Figs. 5,8) is likely due to the large spatial distances of ISU to LSU and VT, and ISU's unique meteorology (e.g., Pacific maritime and polar storm sources), compared to LSU and VT, which are more similar meteorologically (Fig. 3). Seasonal trends in the composition of aerosolized bacteria have been previously reported (e.g., [START_REF] Brodie | Urban aerosols harbor diverse and dynamic bacterial populations[END_REF][START_REF] Bowers | Seasonal variability in bacterial and fungal diversity of the near-surface atmosphere[END_REF][START_REF] Striluk | The effect of season and terrestrial biome Accepted Article Version postprint Comment citer ce document[END_REF], and a small set of studies have observed temporal variations in in the bacterial composition of clouds or precipitation samples taken from the same or nearby locations [START_REF] Kourtev | Atmospheric cloud water contains a diverse bacterial community[END_REF][START_REF] Vaïtilingom | Long-term features of cloud microbiology at the puy de Dôme (France)[END_REF], Huffman et al. 2013, Šantl-Temkiv et al. 2015).

Although all four seasons had significantly different bacterial assemblages in pairwise comparisons (Table 3), the uniqueness of winter is evident through both the small area of its confidence ellipse in the bacterial ordination (Fig. 6), and in the relatively high number of significant ISA winter indicator taxa (Fig. 7). Along with the distinct terrestrial surface conditions of these seasons, a potential driver for this phenomenon is the relatively unique (non-ecotonal) conditions of winter storms. The morphology of winter storms (which were limited to nimbostratus and MCS events) were highly distinct compared to summer, autumnal and spring storms (Fig. 3).

Prediction 1b.

Compared to fungi, the β-diversity of bacterial assemblages at particular sites had a larger β nest component (due to species losses in local assemblages compared to the overall metacommunity species pool) and a smaller β turn component (species replacement in time and space) (Fig. 4). The average bacterial OTU richness of individual precipitation events was only 1,911.4/34,107 = 5.6% of the total richness of the entire 16S rRNA gene dataset, whereas the

Accepted Article

This article is protected by copyright. All rights reserved corresponding fungal OTU richness was much higher: 1377.8/16,597 = 8.3%. This discrepancy may occur because fungal OTUs lost from washout are more readily replaced by other OTUs from local subsidization events.

Richness "hotspots" may contribute to high levels of β nest [START_REF] Knopf | The role of organic aerosol in atmospheric ice nucleation: a review[END_REF]. In our study, the richness sample variance was more than five times greater for bacteria than fungi (1,134,132/205,963 = 5.5) for paired events. The two most bacterially diverse precipitation events occurred at ISU (richness = 3865 OTUs) and LSU (richness = 4253 OTUs). These events were linked to very different meteorological processes (e.g., 72 hr. back trajectory storm origins at 9000m were mP1 and mT1 for ISU and LSU, respectively), suggesting that a variety of aerosolization conditions may result in high levels of bacterial diversity in precipitation.

Nestedness of bacterial assemblages in precipitation may be reinforced by the predicted shorter residence time of INA taxa in the troposphere (Burrows et al. 2009b[START_REF] Stopelli | Ice nucleation active particles are efficiently lost from precipitating clouds[END_REF].

The gammaproteobacterial genus Pseudomonas contains the vast majority of documented INA microbial strains [START_REF] Joly | Ice nucleation activity of bacteria isolated from cloud water[END_REF], and was assigned in our sequence analyses to 203 OTUs (see Methods). When considered alone, these 203 OTUs resulted in much higher levels of precipitation event nestedness (β nest = 0.16; β nest /β Total = 0.32) than when considering all bacterial OTUs (β nest = 0.12; β nest /β Total = 0.16). These results suggest that bacterial assemblages in precipitation events with few Pseudomonas OTUs are compositional subsets of events with many Pseudomonas OTUs-potentially from washout of particular INA strains within storm systems.

The nestedness-prone patterns of Pseudomonas OTUs are graphically depicted in Fig. 11.

Rows (precipitation events) and columns (OTUs) in Fig. 11 A), it would be possible to arrange the matrix so that all black cells were in the region above the designated red curve [START_REF] Atmar | The measure of order and disorder in the distribution of species in fragmented habitat[END_REF]Patterson 1993, Rodríguez-Gironés, and[START_REF] Rodríguez-Gironés | A new algorithm to calculate the nestedness temperature of presence-absence matrices[END_REF]. The matrix temperature index, denoted T, can be used to measure nestedness of multivariate observations, e.g., the microbial composition of precipitation events [START_REF] Atmar | The measure of order and disorder in the distribution of species in fragmented habitat[END_REF]. T decreases with nestedness with T = 0 indicating perfect nestedness. Matrix temperature was calculated for precipitation events using all bacterial OTUs (Fig 11a) and using only the 203 OTUs designated as

Accepted Article

This article is protected by copyright. All rights reserved Pseudomonas (Fig. 11b). We found that the matrix temperature of events when only considering Pseudomonas (T = 4.81) was significantly lower than matrix temperature of events when considering all OTUs (T = 7), based on bootstrap tests (p = 0.008, 1000 iterations, Appendix S1: pg. 97).

To further assess the contributions of individual taxa to bacterial assemblage nestedness, we derived β nest and β turn dissimilarity matrices [START_REF] Baselga | Partitioning the turnover and nestedness components of beta diversity[END_REF] based on the 1698 bacterial genera identified, and generated PCoA ordinations from those matrices (NMDS projections failed to converge for the β turn matrix). The gene library reads for Pseudomonas were more strongly correlated to the PCoA projection of the β nest dissimilarity matrix than for any other genera (R 2 = 0.73, p = 0.001; Supplementary Materials; Table S89). For the remaining bacterial genera, R 2 ≤ 0.38 (Table S89).

Further, Pseudomonas had one of the weakest correlations to the PCoA projection of β turn dissimilarities (R 2 = 0.02, p = 0.553). These results support the premise that patterns of relatively high levels of nestedness and low levels of turnover within precipitation-borne aerosolized bacterial assemblages are driven in large part by selective washout of Pseudomonas OTUs as storm systems develop and age. These conclusions are supported by [START_REF] Stopelli | Ice nucleation active particles are efficiently lost from precipitating clouds[END_REF] who found that "older" clouds (those less recently derived from a storm-generating source) have markedly lower proportions of INA taxa, due to ongoing precipitation events as storms move in space and time.

Hypothesis 2: Precipitation-borne fungal assemblages are correlated with mesoscale processes

In agreement with our heuristic, fungal assemblages were mostly associated with mesoscale spatiotemporal drivers (H 2 ). This outcome was evident in the strong associations of fungal OTU composition to particular spatial locations and moderate altitude characteristics of storms (P 2a ), and the higher turnover (species replacement) contributions to β-diversity in fungi, compared to bacteria (P 2b ). Prediction 2a. Fungal assemblages were strongly associated with factors depicting relatively small changes in space and time. For example, the fungal NMDS projection was more strongly correlated with lower altitude air mass origin classes than the bacterial NMDS ordination (Figs. 5,6; Table 4).

Both PERMANOVA (Table 3) and factor-fitting (Table 4) identified stronger fungal OTU compositional associations for sites compared to continental-scale season effects. Further, ISA

Accepted Article

Comment citer ce document : Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394

This article is protected by copyright. All rights reserved identified few fungal seasonal indicators, and large number of significant site indicators (Table 5).

This outcome is likely due to a number of physical factors including the fact that fungal spores and free-living aerosolized fungal organisms are generally much larger (2-56 µm diam.) than prokaryotic cells (0.2-10 µm diam.) (Amato et al. 2017[START_REF] Griffin | Global-scale atmospheric dispersion of microorganisms[END_REF]). Long-distance atmospheric dispersal of fungal spores is generally considered to be less than 10 km [START_REF] Golan | Long-distance dispersal of fungi[END_REF]. Thus, the fungal composition of precipitation is more likely to be influenced by proximate aerosolization events, including agricultural soil disturbance, and spore dispersal as a part of sexual or asexual reproduction (Wilkinson et al. 2011).

A factor unconsidered by our metacommunity framework is the propensity of precipitation to scrub larger aerosols (> 1.5µm) more efficiently than smaller aerosols [START_REF] Gregory | The microbiology of the atmosphere[END_REF][START_REF] Jung | Approximated solution on the properties of the scavenging gap during precipitation using harmonic mean method[END_REF], although sub-micrometer bioaerosols can be effectively scavenged by precipitation [START_REF] Pratt | In situ detection of biological particles in cloud ice-crystals[END_REF]). Size-dependent scrubbing undoubtedly contributes to the composition of metacommunities in precipitation. Indeed, this factor emphasizes the desirability of INA phenotypes for taxa that may otherwise be washout-resistant and prone to prolonged separation from the earth's surface. Prediction 2b. β turn was the largest contributor to β Total for both bacteria and fungi (Table 2, Fig. 4). Despite this general pattern, however, nestedness and turnover components of β Total differed for bacterial and fungal assemblages. Compared to bacteria, fungal β-diversity at particular sites had higher contributions from turnover (taxa replacement) over time, and lower contributions from nestedness (Fig. 4). This outcome reiterates the importance of mesoscale drivers for fungal assemblages in precipitation, with taxa replacement corresponding to locality-specific temporal variations and environmental filtering (cf. [START_REF] Douda | Legacy of post-glacial colonization affects β-diversity: Insights into local community assembly processes[END_REF]. Macroscale effects in airborne bacteria species sorting and opposing mesoscale mechanisms for fungal assemblages likely contribute to pronounced spatiotemporal variations in bacterial and fungal metacommunities conjoined in precipitation (e.g., Myers et al. 2012).

Taxonomic Composition of Sequence Libraries

Bacteria. Our findings agree with previous studies of aerosolized microbes that have reported dominance of Actinobacteria, Cyanobacteria, Firmicutes, Bacteroidetes, and Proteobacteria,

Accepted Article

Comment citer ce document : This article is protected by copyright. All rights reserved particularly Alphaproteobateria and Gammaproteobacteria (Lindeman et al. 1986, Bowers et al. 2013, Striluk et al. 216, many others). These findings are also in general agreement with the limited number of studies that have examined the microbes in clouds and precipitation [START_REF] Vaïtilingom | Long-term features of cloud microbiology at the puy de Dôme (France)[END_REF], Huffman et al. 2013, Šantl-Temkiv et al. 2015, Amato et al. 2017). Reflecting our work, cloud studies from puy de Dôme, France report high proportions of culturable bacteria in cloud water samples related to the genera Pseudomonas and Sphingomonas, along with notable levels of Methylobacterium, Massilia, Bacillus, Streptomycetes, Curtobacterium, and Arthrobacter [START_REF] Vaïtilingom | Long-term features of cloud microbiology at the puy de Dôme (France)[END_REF]).

Among numerous bacterial indicators for ISU was the soil-associated actinobacterial genus Arthrobacter (Fig. 7). High abundance of this taxon may be due to aerosolization from agricultural lands surrounding ISU. LSU was strongly indicated by the genus Cloacibacterium, a wastewaterassociated, gram-negative, non-motile, facultative anaerobe from the family Flavobacteriaceae. VT was strongly indicated by the Edaphobacter, a gram-negative obligatively aerobic genus in Acidobacteria.

Strongly significant spring indicator genera included Massilia and Sphingomonas (Fig. 8).

Massilia is a major group of aerobic flagellated rods that include rhizosphere and root colonizing bacteria, and has been previously detected in aerosols [START_REF] Blatny | Airborne Legionella bacteria from pulp waste treatment plant: aerosol particles characterized as aggregates and their potential hazard[END_REF][START_REF] Fahlgren | Diversity of airborne bacteria in samples collected using different devices for aerosol collection[END_REF], including cloud-top samples [START_REF] Vaïtilingom | Long-term features of cloud microbiology at the puy de Dôme (France)[END_REF]) and the stratosphere [START_REF] Bryan | Microbial distributions and survival in the troposphere and stratosphere[END_REF]. Sphingomonas was previously found to be dominant in precipitation sampling events in winter (December and November) and spring (April, March, and June) at Puy de Dôme, France [START_REF] Vaïtilingom | Long-term features of cloud microbiology at the puy de Dôme (France)[END_REF].

INP concentrations in precipitation were previously reported to be highest during late autumn months and lowest in the summer at our LSU site [START_REF] Christner | Geographic, seasonal, and precipitation chemistry influence on the abundance and activity of biological ice nucleators in rain and snow[END_REF]Fig. 1B). Our current work, based on a preponderance of LSU samples, supports this finding. The genus Pseudomonas was a significant ISA indicator for autumn in our study, and was at lowest abundances during the summer (Fig. 8). This outcome, however, was driven in part by two unusual events (both LSU autumnal storms), identifiable as isolated points along the far-right side of the first axis of bacterial NMDS ordination (Fig 5a,6a). Read totals ( ± se) for Pseudomonas at these events (43,867 ± 370) were 𝑥 This article is protected by copyright. All rights reserved for Pseudomonas remained higher for august than other seasons with removal of the two events;

Pseudomonas p-values for august, however, increased from 0.01 to 0.07. Notably, Pseudomonas was the most common bacterial genus cultured during repeated sampling of clouds from puy de Dôme, and P. syringae was found to be the most commonly isolated species from this genus [START_REF] Vaïtilingom | Long-term features of cloud microbiology at the puy de Dôme (France)[END_REF]. P. syringae is the most effective naturallyoccurring INP known [START_REF] Maki | Ice nucleation Accepted Article Version postprint Comment citer ce document[END_REF], and frequently occurs in precipitation events globally [START_REF] Sands | The association between bacteria and rain and possible resultant meteorological implications[END_REF][START_REF] Amato | An important oceanic source of micro-organisms for cloud water at the Puy de Dôme (France)[END_REF][START_REF] Monteil | Soil water flow is a source of the plant pathogen Pseudomonas syringae in subalpine headwaters[END_REF], Du et al. 2017).

Fungi. The dominance of Ascomycota, and particularly Basidiomycota, in our samples (Figs. 9,10) agrees with previous work on aerosolized fungi, including investigations of cloud and precipitation microbiology [START_REF] Amato | An important oceanic source of micro-organisms for cloud water at the Puy de Dôme (France)[END_REF][START_REF] Vaïtilingom | Long-term features of cloud microbiology at the puy de Dôme (France)[END_REF], Huffman et al. 2013). Although nonsignificant, a weak trend of higher Linnaean fungal diversity during august and summer months (Supplemental Material, Table S53) was evident and is supported by a number of aerosolized microbial studies (e.g., [START_REF] Burch | Effects of meteorological conditions on spore plumes[END_REF].

As noted repeatedly, trends in fungal indicators were site (not season) dependent. Among abundant, strongly significant taxa were ISU fungal indicators including the order Capnodiales, an epiphytic Dothideomycetes, and the class Tremellomycetes, a saprophytic Basidiomycetes. VT was strongly indicated by Microbotryales, a Basidiomycetes that includes plant pathogens (Fig 9).

The high fraction of unclassified fungal taxa in our data is potentially tied to our use of the 18S rRNA gene V9 region over ITS approaches [START_REF] Lepère | Diversity, spatial distribution and activity of fungi in freshwater ecosystems[END_REF]. While this approach may hamper taxonomic analyses at levels below order, it did not affect OTU-based inferences with respect to our research hypotheses, which were largely concerned with β-partitioning.

Conclusions

Our work clarifies the mechanisms and constraints of bioaerosols in the Earth-atmosphere system by compiling and analyzing the most extensive spatiotemporal dataset of microbial assemblages in precipitation to date. We considered these data using a novel metacommunity framework that predicts dispersal patterns as a function of altitude, based on reported general limitations to distributions of airborne microorganisms. In accordance with our heuristic, we found that bacterial

Accepted Article

This article is protected by copyright. All rights reserved composition of precipitation was more strongly associated with macroscale-driven storm characteristics compared to fungi assemblages, which were associated with mescoscale drivers, including particular spatial locations and altitude characteristics of storms. These patterns were demonstrated using a number of approaches including PERMANOVA, ordination, and β-partitioning.

In the later analyses, bacteria had higher contributions to total β-diversity from nestedness due to nonrandom species losses. Hypothetically, this may be due to the inefficiency of taxa subsidizations to high altitudes following washout events. Notably, bacterial nestedness appeared to be reinforced by general trends within the genus Pseudomonas, which contains the highest number of documented INA species and strains. Conversely, total β-diversity of fungi had higher contributions from turnover over time --due to species replacement, hypothetically driven by local environmental filtering. This article is protected by copyright. All rights reserved Tables Table 1. Sample numbers for 72 precipitation events: 71 events sampled for 16S rRNA gene libraries and 30 events sampled for 18S rRNA gene libraries. The first number in each cell defines the number of 16S samples and the second number gives the number of 18S samples.

Accepted Article

This article is protected by copyright. All rights reserved No. significant non-FDR-adjusted taxa 0 3 9

No. significant FDR-adjusted taxa 0 0 1

Prop. significant non-FDR-adjusted taxa 0.00 0.12 0.14 Prop. significant FDR-adjusted taxa 0.00 0.00 0.02

Accepted Article

Version

Figures

Accepted Article

  Figure 2 presents a metacommunity framework for precipitation-borne microbes based on the information given earlier. Local microbial assemblages are distinguished with hatching and colored arrows show subsidization fluxes. Atmospheric bioparticles

  Storm structure classifications were based on cloud top data retrieved from The National Weather Service (NWS) archive of Geostationary Operational Environmental Satellite-East (GOES-East) infrared satellite imagery, and Next Generation Radar (NEXRAD). Level II and III radar reflectivities were provided by the National Climatic Data Center's (NCDC) website (https://www.ncdc.noaa.gov/nexradinv/). Troposphere temperature profiles were retrieved from The NWS radiosonde data archive of stations located in Lake Charles, Louisiana (KLC) Slidell Muni, Louisiana (KLI), Pocatello, Idaho (KSFX), and Roanoke Virginia (KFCX). Seventy-two-hour backward trajectories of air masses present over Baton Rouge, LA, Blacksburg, VA, and Pocatello, ID at times of collection were calculated using the NOAA Air Resources Laboratory Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model.

  sequential Bonferroni procedure Accepted Article Comment citer ce document : Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394

  context of our work, β nest represents variation in bacterial and fungal local assemblages due to subsetting of richer assemblages (Fig 1; scenario A), whereas β turn represents local assemblage heterogeneity due to taxa replacement in precipitation events over spatiotemporal gradients (Fig 1; scenario B). By definition, β Total , β nest , and β turn are in the interval [0, 1], with larger proportions indicating greater assemblage heterogeneity, and β Total = β nest + β turn . Bootstrapping, based on random sampling of rain events with replacement, was used to derive standard errors for β Total , β nest , and β turn and to perform null hypothesis tests of equal partitioning components. Controls of microbial assemblage nestedness were examined using the matrix temperature approach of Rodríguez-Gironés, and Santamaría (2006).

  Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394

  to both 16S and 18S OTU by rain event data matrices. Because of computational constraints, random forest imputations of single Idaho and VT missing cells (and concomitant PERMANOVAs) were based on only the 1500 most abundant 16S and 18S OTUs. Bray-Curtis dissimilarity (Bray and Curtis 1957) was used as the underlying resemblance matrices in PERMANOVA and ordinations (see below). Bray-Curtis Accepted Article Comment citer ce document : Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394

  Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394

  are arranged to maximize event/OTU cooccurrences (black cells) in the upper left-hand part of the matrices (see Fig S34 in the Supplemental Material for a drastically simplified representation). Under perfect nestedness (see Fig. 1, scenario

  Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394

  much higher than at other events(2,202 ± 361). Indicator values (seeDufrene and Legendre 1998) Accepted ArticleComment citer ce document : Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394

  postprint Comment citer ce document : Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394This article is protected by copyright. All rights reserved

Figure captions Figure 1 .

 captions1 Figure captions

Figure 2 .

 2 Figure 2. A metacommunity framework for the origin, mixing, transport and deposition of aerosolized and precipitation-borne microorganisms. Local microbial assemblages are distinguished with hatching, and colored arrows show subsidization fluxes of microbial assemblages. PBL = Planetary Boundary Layer, SBL = Surface Boundary Layer.

Figure 3 .

 3 Figure 3. Radar charts showing, for 72h back trajectories, proportions of occurrences of precipitationevent airmass origins at two altitudes, and weather morphologies, by site and season. For origins: mP1 = Pacific maritime polar, mT1 = Pacific maritime tropical, mT2 = Atlantic maritime tropical, cP = continental polar and cT = continental tropical. Note that an additional origin class occurs at 1000m.

Figure 4 .

 4 Figure 4. Comparisons of bacterial and fungal beta-partitioning ratios for paired precipitation samples (n = 29) for which both 16S and 18S rRNA gene data were available. Bars are means and error bars are bootstrap standard errors (R = 500). Triple asterisks indicate significant differences between 16S and 18S ratios based on percentile bootstrap confidence intervals (R = 1000) at α = 0.001.Accepted Article

Figure 6 .

 6 Figure 6. Identical NMDS projections to Fig. 5, but with seasons emphasized. See Fig. 5 caption for explanation of plot components.

Figure 7 .

 7 Figure 7. Most abundant taxa in 16S rRNA gene libraries across sites. Specifically, the bars show the six most abundant (high read) phyla, and the 15 most abundant classes, orders, families and genera.Also shown is an "other" category, comprising all other named taxa, and "unclassified" and uncultured" IDs. Asterisks (*) denote a significant ISA taxon and are accompanied by bold font entries in the legend. Double asterisks (**) denote a significant ISA taxon after controlling for FDR and are accompanied by bold text and underlined entries in the legend.Accepted Article

Figure 9 .

 9 Figure 9. Most abundant taxa in 18S rRNA gene libraries across sites. Specifically, the bars show the four most abundant phyla, and the 15 most abundant classes and orders. Also shown is an "other" category, comprising all other named taxa. See Fig. 7 caption for additional details.

Figure 10 .

 10 Figure 10. Most abundant taxa in 18S rRNA gene libraries across seasons. Specifically, the bars show the four most abundant phyla, and the 15 most abundant classes and orders. Also shown is an "other" category, comprising all other named taxa. See Fig. 7 caption for additional details.

Figure 11 .

 11 Figure 11. Matrix temperature (nestedness) plots (Atmar and Patterson 1993) summarizing: (a) all 34107 16S OTUs, and (b) 203 Pseudomonas OTUs at 71 precipitation events. Blackened cells in matrices indicate the occurrence of an OTU at an event. Matrix rows (precipitation events) and columns (OTUs) are arranged to maximize nestedness by maximizing the number of occurrences above the red curve. Events with fewer Pseudomonas OTUs were generally compositional subsets of the richer events, resulting in relatively high levels of event nestedness compared to OTUs from all genera.

Table 2 .

 2 OTU richness and diversity summaries for 16S (n = 71) and 18S (n = 30) precipitation event data by site and season. Entries are means ± SEM. Rows within site and season sections provide comparisons of levels in site and season. Entries within row with the same letter within the same section are not significantly different after controlling for FWE among site or season levels using

	Tukey's HSD or Bonferroni adjustment of confidence intervals (for β-partitioning tests). Rows within sections without letters indicate that site or season effects were not significant in the corresponding omnibus test. Formal β-partitioning comparisons of 16S and 18S were made using paired precipitation samples (n = 29) in which both 16S and 18S rRNA amplifications could be made and are shown in Fig. 4. Site Season ISU LSU VT Autumn Spring Summer Winter 16S Rich. 2185 ± 531 1849 ± 135 2001 ± 176 1844 ± 209 1760 ± 338 1755 ± 148 2141 ± 192 SW-div. 5.3 ± 0.3 4.5 ± 0.2 4.9 ± 0.3 4.5 ± 0.3 4.3 ± 0.4 4.6 ± 0.1 5 ± 0.2 β turn 0.46 ± 0.07 B 0.58 ± 0.07 A 0.56 ± 0.1 AB 0.6 ± 0.02 0.62 ± 0.03 0.58 ± 0.03 0.6 ± 0.02 β nest 0.21 ± 0.01 0.14 ± 0.02 0.11 ± 0.02 0.15 ± 0.03 0.15 ± 0.02 0.1 ± 0.04 0.11 ± 0.03 β Total 0.67 ± 0.02 B 0.72 ± 0.03 A 0.68 ± 0.02 AB 0.74 ± 0.02 0.77 ± 0.02 0.67 ± 0.03 0.71 ± 0.02 18S Rich. 1215 ± 25.5 1442 ± 110.8 1243.5 ± 165.8 1387 ± 262 1374 ± 159 1424 ± 151 1319 ± 110 SW-div. 4 ± 0.2 3.8 ± 0.2 3.9 ± 0.3 3.8 ± 0.4 3.9 ± 0.3 3.9 ± 0.2 3.8 ± 0.2 β turn 0.64 ± 0.12 0.7 ± 0.0 0.67 ± 0.12 0.71 ± 0.02 0.74 ± 0.02 0.66 ± 0.02 0.74 ± 0.02 β nest 0.01 ± 0.02 B 0.06 ± 0.01 A 0.05 ± 0.02 AB 0.08 ± 0.06 0.03 ± 0.02 0.07 ± 0.06 0.04 ± 0.06 β Total 0.65 ± 0.07 B 0.76 ± 0.02 A 0.73 ± 0.07 AB 0.79 ± 0.09 0.78 ± 0.02 0.73 ± 0.09 0.78 ± 0.09 Accepted Article

Comment citer ce document : Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394

This article is protected by copyright. All rights reserved

Table 3 .

 3 PERMANOVA results for 16S and 18S rRNA gene raw reads. Both omnibus (left side of table) and pairwise comparisons (right side of table) are shown. For omnibus tests, significant results are bolded. For pairwise tests, levels sharing a lettered exponent are not significantly different after controlling FWE using Holm's sequential Bonferroni procedure. Two imputed data points were required to address empty cells in PERMANOVAs (see Methods).

	Accepted Article	16S Site Season Site × Season Error 18S Site Season Site × Season Error	Omnibus PERMANOVA F R 2 p-val 2 6.32 Df 0.14 0.0001 3 2.65 0.09 0.0001 6 1.18 0.08 0.1049 61 2 3.10 0.17 0.0002 3 1.46 0.12 0.0572 6 0.91 0.15 0.6932 20	Pairwise comparisons ISU A , LSU B , VT C Autumn A , Spring B , Summer C , Winter D ISU A , LSU B , VT C Autumn, Spring, Summer, Winter

Comment citer ce document : Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394

This article is protected by copyright. All rights reserved

Table 4 .

 4 Results of vector (quantitative predictors) and factor fits (categorical predictors) of NMDS ordinations of 16S (Figs 5a, 6a) and 18S (Figs 5b, 6b) datasets. Statistically significant results are bolded.

	Accepted Article	Vectors dist. traveled 100m dist. traveled 1000m dist. traveled 5000m dist. traveled 6000m dist. traveled 8000m dist. traveled 9000m Cloud top temp. Factors Site Season Air mass 100m Air mass 1000m Air mass 5000m Air mass 6000m Air mass 8000m Air mass 9000m Weather morphology Convective or Stratiform Precip. type (rain or ice)	R 2 0.012 0.062 0.154 0.234 0.218 0.193 0.018 0.130 0.109 0.201 0.039 0.025 0.071 0.061 0.072 0.142 0.031 0.021	16S p-value 0.636 0.103 0.016 0.001 < 0.001 0.001 0.546 0.154 0.025 0.030 0.590 0.849 0.277 0.301 0.181 0.075 0.105 0.175	R 2 0.095 0.111 0.035 0.084 0.109 0.115 0.008 0.278 0.170 0.191 0.251 0.183 0.345 0.217 0.245 0.166 0.053 0.152	18S	p-value 0.257 0.224 0.627 0.318 0.219 0.204 0.915 0.002 0.125 0.063 0.028 0.217 0.003 0.114 0.059 0.511 0.206 0.007

Comment citer ce document : Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394

This article is protected by copyright. All rights reserved

Table 5 .

 5 Summary of indicator species analyses (ISA). Shown for bacterial and fungal reads, for both site and season, are the number of taxa considered in ISA (those occurring at ≥ 2 events) and the numbers and proportions of significant (α = 0.05) indicator taxa. Significant results with and without adjustment for FDR are given.

	Phylum	Class	Order	Family	Genus

Acknowledgments

This research was made possible by grants from the Division of Environmental Biology of the National Science Foundation (NSF DEB 1241069, NSF DEB 1241161/1643288). We thank J.S.

Lockhart, Q.L. Washburn, S. Zorio, C. Baxter, H. Lavendar, J. Farrar, A. Gilliland, G. Hunt, and J.

Tucher for technical and intellectual support, and ISU MRCF for microbial sequencing and seed grants.

DATA AVAILABILITY

Data are available from the Northwest Knowledge Network: https://doi.org/10.7923/6t5b-y124

Accepted Article

Version postprint

Comment citer ce document : Aho, K. (Auteur de correspondance), Weber, Christner, B., Vinatzer, B., Morris, C. E., Joyce, Failor, K., Werth, BaylessEdwards, Schmale, D. (2020). Spatiotemporal patterns of microbial composition and diversity in precipitation. Ecological Monographs, 90 (1), e01394. , DOI : 10.1002/ecm.1394