{'fraction_non_alphanumeric': 0.030241111565181854, 'fraction_numerical': 0.01882003742498871, 'mean_word_length': 5.332606919095614, 'pattern_counts': {'":': 0, '<': 0, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 2, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 2, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Lake microbial communities are resilient after a whole-ecosystem disturbance

Disturbances act as powerful structuring forces on ecosystems. To ask whether environmental microbial communities have capacity to recover after a large disturbance event, we conducted a whole-ecosystem manipulation, during which we imposed an intense disturbance on freshwater microbial communities by artificially mixing a temperate lake during peak summer thermal stratification. We employed environmental sensors and water chemistry analyses to evaluate the physical and chemical responses of the lake, and bar-coded 16S ribosomal RNA gene pyrosequencing and automated ribosomal intergenic spacer analysis (ARISA) to assess the bacterial community responses. The artificial mixing increased mean lake temperature from 14 to 20 °C for seven weeks after mixing ended, and exposed the microorganisms to very different environmental conditions, including increased hypolimnion oxygen and increased epilimnion carbon dioxide concentrations. Though overall ecosystem conditions remained altered (with hypolimnion temperatures elevated from 6 to 20 °C), bacterial communities returned to their pre-manipulation state as some environmental conditions, such as oxygen concentration, recovered. Recovery to pre-disturbance community composition and diversity was observed within 7 (epilimnion) and 11 (hypolimnion) days after mixing. Our results suggest that some microbial communities have capacity to recover after a major disturbance.


Introduction
Disturbances, defined as events that trigger a discrete change in the physical or chemical environment that may affect the local community (Glasby and Underwood, 1996), have a key role in the structuring of ecosystems (Sousa, 1984;White and Pickett, 1985). Characteristics and ecological consequences of disturbances vary immensely (Fraterrigo and Rusak, 2008), hindering advances in the theoretical and predictive models (Sousa, 1984;White and Jentsch, 2001). As there is an inverse relationship between disturbance frequency and magnitude (large events seldom occur, small events occur frequently (Sousa, 1984;Romme et al., 1998;Turner and Dale, 1998;White and Jentsch, 2001), extreme disturbances are inherently difficult to study in nature. Nonetheless, observations after large disturbances are crucial for understanding the susceptibility of ecosystems to extreme events such as those expected under altered climate scenarios (for example, Katz and Brown (1992) and works cited).
Microbial communities are at the heart of all ecosystem functions, and thus their responses to disturbances may influence ecosystem recovery (Allison and Martiny, 2008). Some have advised that microbial variables be included in predictive models of ecosystem change (for example, Arhonditsis and Brett, 2004;McGuire and Treseder, 2010;Sarmento et al., 2010). But, before this can occur, it is necessary to understand the robustness of microbial assemblages to disturbances (Allison and Martiny, 2008). In ecology, community robustness is comprised of resistance, defined as the ability to withstand change in the face of a disturbance, and resilience, the pace of recovery (if any) following a disturbance (Pimm, 1984).
Water column mixing is a disturbance to microbial communities because it disrupts the physicalchemical gradients created by thermal stratification known to define niches for microorganisms (for example, Heaney and Talling, 1980;Vincent et al., 1984;Ovreas et al., 1997;Cytryn et al., 2000;Fenchel and Finlay, 2008). From an ecological perspective, bacterioplankton responses to mixing may provide unique insight into understanding disturbances. As it disrupts gradients in the water column and occurs seasonally, mixing affects spatial and temporal drivers of microbial communities. Further, both microbial and environmental responses to mixing are tractable. High-throughput sequencing and fingerprinting methods can be used to assess changes in microbial community diversity, while environmental sensors can be used to quantify spatial and temporal changes in the environment. The alignment of these measurements allows a comprehensive perspective of microbial responses to an ecosystem-level disturbance. Therefore, lake mixing and microbial communities together create an interesting model system for understanding disturbance ecology.
The work described here builds on a series of studies focused on the response of freshwater microbial communities to mixing (Shade et al., 2007(Shade et al., , 2010a(Shade et al., , b, 2011Jones et al., 2008). In an earlier study, we observed that patterns of bacterial community succession in a temperate eutrophic lake were linked to spring and fall mixing events (Shade et al., 2007), suggesting an importance of lake mixing for community dynamics. In a second study of a small, darkly stained sub-tropical lake, we observed a surprising recovery among bacterial communities following mixing events caused by typhoons, with repeatable trajectories of community change across typhoon events and across years Shade et al., 2010b). The results of these two surveys prompted us to ask whether bacteria in temperate lakes that mix seasonally would be similarly resilient to an episodic mixing event. Therefore, we conducted an enclosure experiment in a temperate humic lake to evaluate the effect of mixing and specific mixingassociated environmental changes on communities in the epilimnion (near surface stratum) and hypolimnion (bottom stratum) . Both epilimnion and hypolimnion communities returned to their control composition within 10 days, demonstrating community robustness despite differences in initial compositions and in disturbance characteristics. Finally, we conducted a reciprocal transplant experiment to ask how epilimnion and hypolimnion bacterial communities responded to environmental conditions in the opposite strata (Shade et al., 2010a). Consequently, we identified possible 'generalist' community members, persistent in both epilimnion and hypolimnion, that could serve as pioneer species after mixing. The results from the transplant experiment suggested that lake bacterial communities harbor members that initiate post-disturbance succession and eventual recovery.
Following from the results of this series of studies, we asked whether freshwater bacterial communities would be robust (resistant and resilient) to an artificial mixing event that occurred during summer stratification, which is a time when communities are compositionally distinct from communities that occur during seasonal spring and fall mixing events. Mixing at the peak of summer stratification is a disturbance that would not happen naturally, and therefore represents an extreme event for summer microbial communities. Thus, this disturbance experiment provided a unique opportunity to measure both community and ecosystem responses to a large, unlikely disturbance in a relatively wellstudied system. Summer mixing also allowed us to observe changes in the bacterial community separate from other seasonal changes that are concurrent with spring and autumn mixing (for example, ice-off). Further, a whole-ecosystem experiment avoided limitations of mesocosm and enclosure experiments (Schindler, 1998). One common criticism of mesocosms is that they fail to reproduce conditions in nature and provide limited insight, and so it was possible that the community responses we previously observed in mesocosm experiments might not occur in the whole-lake. For these reasons, we wanted to pursue a whole-ecosystem experiment to understand microbial community responses to an unlikely, extreme mixing event.
Here we report the bacterial community responses to the disturbance of an entire small temperate lake, North Sparkling Bog (Boulder Junction, WI, USA, Supplementary Table 1). We assessed bacterioplankton communities using 16S ribosomal RNA gene tag pyrosequencing to gain an unprecedented depth of sampling of the lake microbial diversity.

Materials and Methods
Experiment, site, environmental sensors and sampling The study lake, North Sparkling Bog (46100'17.47' N,89142'18.37' W), is darkly stained with high dissolved organic carbon (average 9.5 mg l À 1 ). Though this is a relatively shallow lake (4.5 m maximum depth), the darkly stained color and surface area (0.46 ha) to depth ratio (1022 m) causes strong stratification during most of the ice-free season, which creates strong vertical gradients in temperature and dissolved oxygen, and overall high stability. North Sparkling Bog is a Global Lakes Ecological Observatory Network site. It has a buoy instrumented with environmental sensors (see below); sensor data are available online (www.gleon.org) and lake characteristics are given in Supplemental  Table 1. Whole-lake destratification was performed by mixing epilimnion and hypolimnion waters with the aid of buoyant forces . Briefly, two rubber membranes with inflatable tubing on the outer rim (8.25 m diameter) were oscillated in the water column. When the tubing was filled with compressed air, the membranes rose to the lake surface. After air was released from the tubing, the membrane sank to the bottom. This repeated process of inflation and deflation caused turbulence and internal waves, and eventually resulted in complete mixing.
The lake was stratified in oxygen and temperature before the artificial mixing that began on 02 July 2008. An autoprofiling sonde (Yellow Springs Instruments, Yellow Springs, OH, USA; 6600 V2-4) measured hourly profiles of dissolved oxygen and temperature, and these observations were used to assess progression of lake mixing (data not shown). The buoy was also outfitted with a thermistor string (Apprise Technologies, Inc., Duluth, MN, USA) and a sensor that measured aqueous and atmospheric carbon dioxide concentrations every 2 h (Browne, 2004). The mixing manipulation continued until the surface and bottom-water temperatures differed by less than the accuracy of the thermistors (0.2 1C), on 10 July 2008 (the last day of mixing, referred to throughout as Mix), signifying that complete mixing had been achieved.
Samples were collected weekly or twice weekly during June, July and August of 2007, 2008 and 2009, respectively, for standard water chemistry and microbial community analyses. Manual measurements of oxygen and temperature were collected at every sampling time (Yellow Springs Instrument 550 A) and used to identify the thermocline. Integrated epilimnion and hypolimnion samples were collected as previously described . In 2008, integrated thermal-layer samples were collected more frequently (usually daily) during and after the whole-lake manipulation. In 2008 during the experimental period, we also collected depth-discrete epilimnion (0 m) and hypolimnion (4 m) water samples using a peristaltic pump that was rinsed between each sample. Depth-discrete samples were collected in this manner even when the lake was completely mixed (and therefore no epilimnion or hypolimnion layers existed). Thus, on the day that complete mixing was achieved (Mix), we considered the 0 m and 4 m water samples to represent epilimnion and hypolimnion water, respectively, even though the temperature throughout the water column was homogeneous. Collected water was preserved for water chemistry (total and dissolved nitrogen, total and dissolved phosphorus, dissolved organic and inorganic carbon, water color) using the North Temperate Lakes Long Term Ecological Research standard protocols (http://lter.limnology.wisc.edu/). Water for bacterial community analysis was filtered onto a 0.2 mm polyethersulfone membrane (Pall, New York, NY, USA), frozen at À 20 1C at the field station, and then transferred to À 80 1C at the main laboratory in Madison.
During 2008, we also measured concentrations of several terminal electron acceptors and donors that are hallmarks of distinct biogeochemical transformations mediated by microorganisms under anoxic conditions (for example, iron reduction, sulfate reduction and methanogenesis), throughout the ice-free season. Data presented here were collected on 01 July (day À 9), 10 July (Mix), 13 July (day 3), 15 July (day 5), 17 July (day 7), 19 July (day 9), 21 July (day 11) and 30 July (day 20). All parameters were measured at 4 m. We also measured five-point water column profiles of carbon dioxide and methane. For the carbon dioxide and methane profiles, one sample was always collected at 0 m, one sample was always collected at 4 m to match the microbial samples, and three additional samples were collected at intermediate depths that best characterized the inflection points in dissolved oxygen over the oxycline. These middle three sample locations varied with the thermal and oxygen structure of the water column, and therefore were not the same depths across all sampling times. When the water column was mixed (no oxygen inflection point was observed), methane and carbon dioxide samples were collected equidistant from the surface and bottom at 0, 1, 2, 3 and 4 m.

Water chemistry
Water chemistry analyses were performed by the North Temperate Lakes Long Term Ecological Research, according to their standard protocols (http://lter.limnology.wisc.edu/), except for iron, sulfide, carbon dioxide and methane. Samples for dissolved Fe were acidified with HCl to pHo2 and analyzed colorimetrically using ferrozine (Stookey, 1970). Samples for dissolved sulfide were stabilized by addition of 1% (wt vol À 1 ) zinc acetate (to convert dissolved sulfide to ZnS), and analyzed colorimetrically via the methylene blue method (Cline, 1969). Carbon dioxide and methane were measured using gas chromatography (Shimadzu GC-14A, Kyoto, Japan), equipped with a flame ionization detector and a methanizer, which allowed for analysis of both CH 4 and CO 2 (which is converted to CH 4 by a methanizer) in a single injection.
Rarefied Faith's phlyogenetic diversity (PD, (Faith, 1992)) and rarefied richness were calculated in QIIME. Pielou's Evenness, J, was calculated based on Shannon diversity, H (J ¼ H/log (richness)). The Grubb's outlier test (outliers package (Komsta, 2011)) in the R environment for statistical computing (R Development Core Team, 2010) was used to determine if the day of complete mixing values of richness, Pielou's evenness or PD were significantly different from the rest of the observations in the series. Grubb's test calculates a score for the outlier by subtracting the proposed outlier value from the set mean and then dividing by the s.d. (G statistic), and also by comparing the ratio of variances of the set with and without the proposed outlier included (U statistic). A paired, two-sided t-test was used to determine if epilimnion and hypolimnion trajectories of richness, Pielou's evenness or PD were significantly different from each other.
Non-metric multidimensional scaling analysis was performed with UniFrac weighted distances (Lozupone et al., 2007) vegan package in R (Oksanen et al., 2011). The results are not sensitive to choice of metric: in this case, weighted UniFrac performed similarly to unweighted UniFrac and also Bray-Curtis similarity, suggesting that for this ecosystem patterns in community membership and community structure are highly correlated (see Supplemental Online Materials). The R environment was used for calculation of evenness and creation of bar charts and heatmaps by phyla. The vegan package in R was used for correspondence analysis, Mantel tests and Procrustes analysis.
For the evaluation of the species abundance distributions (Magurran, 2007), we fitted log-normal models to the ranked abundance of each taxon in the pyrosequencing data, assuming an underlying Poisson distribution (vegan package in R). Null, preemptive, log-normal, Zipf and Zipf-Mandelbrot models were each fit to the data (Wilson, 1991). Null and pre-emptive models did not produce good fits, as evaluated by Akaike information criterion, Bayesian information criterion and deviance. However, the remaining models had comparable fits. The log-normal model was selected because it made the fewest assumptions about the underlying ecology of the data. Species abundance distributions for other models reproduced the same shift in the community structure at the day of mixing.
To examine bacterial community change before, during and after the mixing event, we used automated ribosomal intergenic spacer analysis (ARISA) to fingerprint bacterial communities from integrated epilimnion and integrated hypolimnion water samples (Fisher and Triplett (1999) collected during June, July and August, with modifications provided in Shade et al. (2011)). This resulted in ARISA profiles for the following samples: 17 in 2007, 60 samples in 2008 and 17 in 2009. Note that this integrated sampling strategy was different than what was used for the depth-discrete samples that were pyrosequenced. The demarcation of epilimnion and hypolimnion samples was determined by thermal profiles measured at the time of sampling, and the thermocline was included in the epilimnion sample. Thus, integrated epilimnion and integrated hypolimnion together included the whole-lake column.
ARISA profiles from integrated thermal-layer samples were used to determine to what extent physical displacement of OTUs via water motion could explain observed changes in community composition, during the July 2008 manipulation. We compared day À 10 samples (30 June 2008) with Mix (10 July 2008) to assess the proportion of Mix community members that could be attributed to OTUs present in the epilimnion or hypolimnion before the water column was mixed. Any community members that were present in the Mix sample but were not observed in the pre-mixing samples were designated as 'Unknown origin.' We used a Venn analysis, pie charts and hierarchical clustering with a heatmap of OTU occurrences to visualize differences in pre-mix and Mix communities (R Development Core Team, 2010).
Next, we performed simulations of community richness and evenness to determine the probability of observing the Mix community given the pre-mix composition. For the simulation, we used re-sampling without replacement of integrated epilimnion and integrated hypolimnion communities in 'mixing scenarios' of the following proportions (epilimnion/hypolimnion contribution): 10/90, 20/80, 30/70, 40/60, 50/50, 60/40, 70/30, 80/20 and 90/10 (Supplemental Methods). For 1000 resamplings at each scenario, we report a P-value to describe how often we encountered the observed Mix richness and evenness in the simulated re-sampling. Using the ARISA data from the integrated water samples, we could evaluate whether the community composition at Mix could be explained by different degrees of physical displacement of community members. Richness was compared rounded with the nearest whole number, and evenness was compared with three significant digits. More details for this analysis are provided in Supplemental Materials.
We also used ARISA fingerprints to compare interannual dynamics in North Sparkling Bog the year before (2007, June through October) and the year after the manipulation (2009, June and July). These samples were fingerprinted from integrated epilimnion and integrated hypolimnion water. If the lake bacterial communities recovered after the manipulation, we would expect the dynamics after the manipulation to be generally consistent across years because of the seasonality previously observed in North Temperate bog bacterial communities Shade et al., 2008). To assess community variability, we calculated permuted analysis of betadispersion of pairwise Bray-Curtis similarities (vegan package (Oksanen et al., 2011) and (Anderson et al., 2006)). This analysis calculates dispersion from a mean compositional centroid within a categorical group. In this case, groups were months within years (for example, July 2007). As epilimnion and hypolimnion communities vary differently through time (Shade et al., 2007), epilimnion and hypolimnion samples were analyzed separately.
ARISA fingerprints were further used to calculate a rate of compositional change (a.k.a. similaritydecay) in the community. The slope of the log-linear regression between Bray-Curtis dissimilarity and days between observations is a rate of compositional change (Nekola and White, 1999;Soininen et al., 2007); this is similar to the strategy used to assess rate of change in the previous mesocosm experiment . We compared rates across the 2007-2009 June-August time series to the 2008 post-mixing time series (July-August), and also to the rate of change in temperature over the same time period. The temperature rate was assessed in the same way, by regressing the change in temperature against the time between observations.

Results
Physical and chemical changes after mixing Artificial mixing continued for 8 days until temperature was the same throughout the water column. As compared with the previous ice-free season (2007), mixing disrupted the normal physical and chemical structure of the lake, creating a pulse disturbance and exposing bacterial communities to markedly different environmental conditions. These different conditions included increased bottomwater temperature from 6 to 201C and increased dissolved oxygen concentration from below detection (o0.01 mg l À 1 ) to near 3 mg l À 1 (Figures 1a and b). The vertical oxygen gradient was homogenized gradually, but then returned within 2 days to its pre-manipulation anoxic conditions (12 July 2008, day 2) (Figures 1c, d).
We consider 'recovery' as the point in which the biology or chemistry matched the pre-mixing time point, day À 9. During the recovery period of the lake, stratification was relatively weak and more susceptible to weather-induced mixing, which is not typical for this typically strongly stratified lake ( Figure 1a). We observed a spontaneous secondary Lake microbial communities are resilient A Shade et al mixing event 4 days after treatment had ceased (day 4). This secondary event occurred because overall the manipulation warmed the lake and created nearly isothermal conditions. The lack of thermal stability in the water column made the lake unusually susceptible to weather events that encourage mixing, such as temperature changes or wind.
To the best of our knowledge (and based on inspection of thermal profiles collected from the time that instrumented buoys were first deployed on similar lakes, available at http://www.gleon.org), a spontaneous, full water column mixing event has never been observed in mid-summer in North Sparkling Bog.
To understand any possible functional microbial responses that may be attributed to lake mixing, we measured the concentrations of various hypolimnion respiration products over the mixing experiment. Measurements of hypolimnion iron (Fe 2 þ and Fe Total ), sulfur (H 2 S) and nitrogen (NH 4 ) respiration products decreased to a minimum on the day complete mixing was achieved (Mix), and then rapidly increased, generally following the overall response of dissolved oxygen (Figure 2a). SO 4 had a 'lag' response that correlated with the post-mixing increase in H 2 S. Additionally, Fe 3 þ and NO 3 did not vary substantially over the course of the experiment, and were both present in comparably low concentrations.
We also considered the dynamics of gases CO 2 and CH 4 , two microbial products, over the course of the mixing experiment. CO 2 concentrations in the surface waters during the mixing experiment were markedly higher than CO 2 concentrations in the month prior and were more comparable with concentrations observed in the autumn during fall mixing (Figure 2b). Methane and CO 2 , two products of microbial respiration that often accumulate in the hypolimnia of stratified lakes, were homogenous throughout the water column at low concentrations at Mix, but re-accumulated at the lowest depth (4 m) by day 20 post-mixing (Figures 2c-e).
Comparing patterns in ARISA fingerprinting to patterns in pyrosequencing We first compared bacterial community composition measurements made with automated ribosomal intergenic spacer analysis (ARISA) to those generated via tag pyrosequencing from the same depthdiscrete 0 m and 4 m July 2008 water samples, to ensure that they revealed the same community patterns. Using Bray-Curtis similarity for the ARISA fingerprints and weighted UniFrac distance for the pyrosequencing data, we observed similar patterns and results across both data sets (Mantel Pearson's R ¼ 0.79, P ¼ 0.001 on 1000 permutations). We found a similar agreement across data sets by Procrustes correlation between correspondence analyses sample scores (R ¼ 0.75, P ¼ 0.003 on 1000 permutations). This analysis implied that community patterns generally were consistent across these two c-e were measured using standard analyses of depth-discrete, manually-collected samples.

Changes in bacterial diversity after mixing
To characterize community robustness, we analyzed important ecological aspects of community diversity: richness (number of distinct taxa), evenness (equitability of taxa), structure (representation of dominant and rare taxa), PD, the number and relationship of phylogenetic lineages) and composition (relative abundance of specific taxa defined as phylotypes at the 97% sequence similarity level), using the tag pyrosequence data. Changes in community composition, loss of species richness and gain in species richness relative to the pre-mix conditions (possible bloom or growth of populations) were used as indicators of community responses to the disturbance. July 2008 (the month of the mixing experiment) samples from water collected at depth 0 m (epilimnion) and depth 4 m (hypolimnion) were used to assess immediate community responses to mixing. We predicted that the epilimnion (0 m) and hypolimnion (4 m) communities would be indistinguishable under mixed conditions. However, richness and evenness were different, whereas PD was comparable across depths ( Figure 3). This suggested that despite thermal homogeneity of the lake, each stratum maintained a separate community structure, possibly driven by species sorting by environmental drivers other than temperature. Such a situation could arise if the time scale of vertical water movement (likely driven by convective mixing rather than wind, (Read et al., 2012)) was slower than growth or loss rates.
For each aspect of diversity evaluated and for each depth (0 and 4 m) of the lake, the bacterial communities recovered within days to pre-perturbation states (Figure 3, Figure 4). There was an opposite response in diversity in each of the two stratum ( Figure 3). The epilimnion communities became richer and included taxa that were more phylogenetically divergent after the mixing (Mix). All measures of epilimnion diversity peaked on day 3 Log-normal models fit to ranked taxa abundance to assess changes in the rarity and dominance of members in the (a) epilimnion (0 m) and (b) hypolimnion (4 m). Red arrows below the x axis show the shift in community structure after mixing, and black arrows show recovery. Each line color is a community structure at a different sampling time. In the epilimnion, communities experienced a second, greater shift in structure on day 3, shown by the two red arrows.
Lake microbial communities are resilient A Shade et al post-mixing, and then returned to pre-mixing levels ( Figure 3, open symbols). Pielou's evenness on day 3 post-mixing was detected as significantly different from the remaining July 2008 time series by the Grubb's outlier test (G ¼ 2.04, U ¼ 0.002, Po0.001).
Together, these data suggest that the maximal change of the epilimnion community to the mixing experiment occurred not on the day of complete mixing, but after a 3-day lag.
In comparison, the hypolimnion communities were less rich (Grubb's G ¼ 1.97, U ¼ 0.072, P ¼ 0.006), contained fewer phylotypes (G ¼ 1.93, U ¼ 0.102, P ¼ 0.012), and were less even (G ¼ 3.145, U ¼ 0.017, Po0.001) on Mix, but rapidly recovered to pre-mixing levels by day 3 (Figure 3, closed symbols). These analyses suggest that the hypolimnion community changed maximally around Mix. As an overall comparison of the diversity time series, epilimnion and hypolimnion communities were different in evenness and PD, but they were not different in richness, suggesting that though the number of OTUs observed at either depth were comparable, each depth contained a different representation of taxa (evenness T ¼ 2.93 Po0.05; PD T ¼ À 3.64, Po0.05; richness P ¼ 0.27) Species abundance distributions were used to investigate changes in the bacterial communities' structures (Magurran, 2007;McGill et al., 2007).
There was a strong shift in the structures of both epilimnion and hypolimnion communities (Figure 4). The epilimnion became more even (taxa shared similar relative abundances) during and immediately after the manipulation and the largest shift was between Mix and day 3 (Figure 4a). The hypolimnion became less even (loss of richness and dominance by a few taxa, Figure 4b) during the manipulation, and largely recovered by day 3. The recovery of the shape of the species abundance distributions further supports community recovery.

Changes in bacterial community composition after mixing
We next assessed changes in epilimnion and hypolimnion composition (defined as the relative abundances of OTUs) before and after lake mixing, using tag pyrosequence data ( Figure 5). The epilimnion community experienced a large change between day À 9 pre-mixing and Mix, a second smaller change between Mix and day 3 post mixing, and then another large change between days 3 and 7 post-mixing (Figure 5a). After day 7, the post-mixing epilimnion community was most similar in composition to the pre-mixing state, suggesting epilimnion community recovery (weighted UniFrac distance between epilimnion day À 9 and day 7 ¼ 0.13, Supplementary Figure 1a). The hypolimnion community in North Sparkling Bog changed directionally until it returned to a composition similar to pre-mixing by day 11 postmixing (Figure 5b; weighted UniFrac distance between hypolimnion day À 9 and day 11 ¼ 0.17; Supplementary Figure 1b). As described above, the epilimnion and hypolimnion communities became less distinct during mixing (weighted UniFrac distance ¼ 0.15), but then diverged quickly as thermal stratification was re-established (Figures 5c, d). Epilimnion and hypolimnion communities remained 15% different by weighted UniFrac distance, though the dissimilarity between communities at the two depths reached a minimum at Mix (Figure 5d). Further, similarity between the epilimnion day 20 and day À 9 communities was greater than the similarity between the same days in the hypolimnion, suggesting that the epilimnion community recovered more completely than the hypolimnion.
There were distinct responses to mixing among phylogenetically coherent groups defined at the phylum or Proteobacteria class level (Figure 6). In the epilimnion, Gamma-Proteobacteria bloomed briefly on day 3 post-mixing, possibly responding to nitrogen or other nutrients up-welled from the hypolimnion after a lag for growth ( Supplementary  Figure 2), suggesting that there were community changes that could be attributed to biology. In the hypolimnion, populations of Acidobacteria and Delta-Proteobacteria were below detection at Mix, Phylum-level (and Proteobacterial class level) changes in epilimnion (0 m) and hypolimnion (4 m) composition through the lake mixing experiment. Relative abundances of phylotyped reads are shown, identified to their closest match in the Ribosomal Database Project. Asterisks indicate that no members associated with the phylum were detected. In the x axis, orange circles identify the pre-mix time point (day À 9), green circles identify the day complete mixing was achieved (Mix) and yellow circles identify the recovery time point (day 20). See Supplementary Figure 3 for a different perspective of phylum-level dynamics, normalized by total occurrences within each phylum to visualize maximum changes. Lake microbial communities are resilient A Shade et al but re-appeared by day 3. Several members of these latter phyla are strict anaerobes, and may have either died or were protected within a micro-niche when the hypolimnion was oxygenated. Bacteriodetes also decreased in the hypolimnion and gradually recovered within 20 days. Other phyla, such as hypolimnion Chlorobi and epilimnion Firmicutes, had dynamics that did not obviously correspond to the mixing experiment.
To compare similar patterns of change across phyla, we standardized each phylum's dynamics to its maximum observed abundance (for example, a dominant phylum that changed twofold would be scaled the same as a rare phylum that also changed twofold, Supplementary Figure 3). In the epilimnion, members of less-abundant phyla generally responded to mixing more dramatically than dominant phyla. In the hypolimnion, some rare phyla (for example, Gamma-Proteobacteria and Firmicutes) bloomed after the mix, while some abundant phyla (for example, Bacteroidetes and Acidobacteria) dropped below detection on Mix and recovered by day 3 post-mixing.
Understanding community changes at mixing: displacement and ecology To determine the separate contributions of physical displacement of cells and biological processes (that is, active migration or growth of OTUs), we used ARISA fingerprints of integrated epilimnion and integrated hypolimnion water samples (Figure 7). ARISA OTUs that were not detected before the manipulation (day À 10, note that this is 1 day earlier that the day À 9 time point used for pyrosequencing) were unique to the 'Epi Mix' and 'Hypo Mix' communities ( Figure 7a). These unique Mix OTUs were of unknown origin, and represented 11% of the epilimnion and 12% of the hypolimnion communities ( Figure 7b). Further, the epilimnion mix and hypolimnion mix communities were more similar to each other than they were to either the epilimnion or the hypolimnion community before mixing (Figure 7c).
Finally, we performed re-sampling of the day À 10 communities in various scenarios of mixing to assess how likely it was to observe the Mix communities given the day À 10 composition. We re-sampled communities one thousand times, including epilimnion and hypolimnion proportions of 10/90, 20/80, 30/70, 40/60 and 50/50, and assessed richness (assessed by the number of unique ARISA fragment lengths) and Pielou's evenness for each re-sampling. For these tests (a total of nine mixing scenarios of different proportions of epilimnion and hypolimnion contribution), we never observed the actual ARISA-based richness and evenness after 1000 resamplings at each mixing scenario (Supplementary Figure 4). Together, these results suggest that there were changes in composition that could not be attributed only to physical displacement of cells, but that may be instead explained by ecological processes.

Long-and short-term effects on bacterial community dynamics
To evaluate whether the mixing experiment caused prolonged changes in bacterial community dynamics, we assessed seasonal variability using community ARISA fingerprints from 2007 through July 2009. Because bacterial communities in north temperate bog lakes were previously shown to change seasonally Shade et al., 2008) we used permuted analysis of betadispersion of Bray-Curtis similarities to address whether there were overall monthly differences in variability between the year of mixing and the other 2 years. We detected no differences in overall variability across months and years (epilimnion P ¼ 0.35 on 999 permutations, hypolimnion P ¼ 0.10, Supplementary Figure 5). This suggests that the mixing experiment in July 2008 did not have prolonged or permanent effects on bacterial community dynamics.
We then asked whether the mixing manipulation influenced the bacterial community rate of change, by comparing rates observed in the epilimnia and hypolimnia from 2007, 2008 and 2009. We calculated the rate of change (similarity-decay) by regressing the Bray-Curtis dissimilarity (calculated using ARISA profiles) against the time interval between samples, for all possible pairs of samples within each time series, and recording the slope of the regression (Supplementary Figure 6). The hypolimnion communities sampled in 2008 post-mixing changed faster than any of the other communities. The bacterial community rate of change was correlated with the rate of change in temperature (r ¼ 0.94, P ¼ 0.005, n ¼ 6). Taken together, these findings suggest that the elevated temperature in the hypolimnion contributed to a faster rate of community change than is normally observed in the cold bottom waters of this bog lake.

Discussion
We employed a whole-ecosystem manipulation to ask whether microbial communities could recover after a large disturbance, and to ask how responses of microbial communities compared with changes in environmental conditions. Bacterial communities were not resistant to summer mixing, but were highly resilient. The communities in the epilimnion and hypolimnion both changed as a result of the disturbance, but recovered by days 7 and 11 postdisturbance, respectively. Though microbial communities in the epilimnion and hypolimnion had contrasting responses to the disturbance (for example, epilimnion increased in diversity and the hypolimnion decreased), they had similar patterns  Figure 7 Comparison of composition across integrated epilimnion and integrated hypolimnion communities fingerprinted with ARISA, collected 2 days before the mixing manipulation began (day À 10) and on the day that complete water column mixing was achieved (Mix). (a) ARISA OTUs detected in and shared across the epilimnion and hypolimnion on day À 10 and Mix. Numbers in parentheses show the proportion of detected OTUs belonging to the group. (b) For the epilimnion and hypolimnion Mix communities, the proportion of OTUs that originated in the epilimnion on day À 10, OTUs that originated in the hypolimnion on day À 10, OTUs that were associated with both layers on day À 10, and OTUs that were uniquely observed on the day of Mix (labeled as 'unknown origin'). (c) Heatmap and cluster analysis of community composition of epilimnion and hypolimnion day À 10 and epilimnion and hypolimnion Mix. OTUs are in columns (ARISA OTU IDs are provided along the x axis) and rows are separate communities. OTUs are clustered based on similar abundance and occurrence patterns across the four samples. of recovery. Although the manipulation seemed to increase the rate of community change in the hypolimnion during summer 2008, the communities did not exhibit prolonged changes in dynamics after the mixing experiment, as assessed by comparing overall variability in the year of mixing to ARISA fingerprints from 2007 through July 2009. The richness and evenness of the integrated epilimnion and integrated hypolimnion at Mix were not the same as any simulated mixing scenarios, suggesting that vertical transport of cells in the water column was not enough to explain differences between the two layers (between day À 10 and Mix). We also observed changes in individual phylotype and OTU relative abundances that could not be explained by cell transport in the water column. For example, at mixing, Alpha-Proteobacteria tags increased in relative abundance in the hypolimnion, while other tags previously abundant in the hypolimnion fell below detection (for example, Acidobacteria, Delta-Proteobacteria). Several members of these latter phyla are strict anaerobes, and were either protected in micro-niches or possibly died when the hypolimnion was oxygenated. Interestingly, the Beta-Proteobacteria decreased in the epilimnion as a proportion of total tags, but changed minimally in the hypolimnion. Thus, we infer that biological processes (for example, growth, death, revival from a dormant state) contributed to the observed temporal dynamics. For example, dormant microorganisms persisting in a sediment seed bank may have a role in community resilience after disturbance (Jones and Lennon, 2010), but further study is needed to determine the role of inactive members, as well as extremely rare members, to overall community recovery.

Hypolimnion response and recovery
Hypolimnion communities changed more drastically than epilimnion in composition, richness and phylogenetic diversity. The recovery of hypolimnion composition coincided with the depletion of dissolved oxygen and the re-accumulation of anaerobic respiration products such as methane, carbon dioxide and other reduced molecules (for example, ferrous iron, ammonium and sulfide). As would be expected, accumulation of sulfide coincided with a decrease in sulfate, several days after Mix. Taken together, these observations are consistent with the transitions from aerobic to anaerobic respiration that we would predict from first-principles of biogeochemistry (for example, electron-donor and -acceptor pairs). Thus, the microorganisms and their respiration products helped to create the environmental conditions that eventually harbored a recovered hypolimnion community. This demonstrates a fundamental feedback loop for the hypolimnion community: microorganisms help to create the strong oxidation-reduction gradients, and this gradient then dictates the observed community robustness to mixing by acting as an environmental filter that drives species sorting.
In our previous mesocosm experiment, the hypolimnion community also recovered as oxygen concentration decreased, despite that temperature immediately recovered to pre-disturbance conditions . Notably, the hypolimnion temperature was approximately 6 1C during the mesocosm experiment (the same as the ambient hypolimnion temperature) and the community recovered after 10 days, while in the whole-lake manipulation the temperature was 20 1C (elevated as compared with the pre-disturbance 6 1C) and the community recovered after 11 days. This indicates that oxygen concentration is more important than temperature in determining hypolimnion community composition. This does not exclude the possibility that other, unmeasured environmental changes also were important in driving the microbial community, but oxygen proved to be the most important among those that we measured. However, elevated temperatures appeared to accelerate the hypolimnion community changes after mixing. Recent estimates of bacteriovore grazing rates in a similar small bog lake suggested bacterial growth rates were on the order of 0.33 day even at 9 o C (Hahn et al., 2012). This is consistent with the community composition in the recovering hypolimnion turning over during a period of days.
Summer lake mixing as a novel, large disturbance We have provided several lines of evidence that our mixing experiment represented both a novel and extreme disturbance to lake bacteria. The first is the drastic changes in temperature and oxygen, key drivers of epilimnion and hypolimnion aquatic bacterial communities, respectively . The temperature in the hypolimnion increased from 6 to 20 1C. The historical ice-free hypolimnion temperature has likely never exceeded 10 1C. The oxygen in the hypolimnion was increased from 0 to near 3 mg/l. Stable bog lakes such as North Sparkling Bog remain anaerobic in the bottom waters with the exception of seasonal mixing (fall and spring). Both of these environmental changes were major deviations from seasonal norms. Second, the known seasonality of the microbial communities in North Sparkling Bog  provides additional evidence for the novelty of the mixing event. Because North Sparkling Bog bacterial communities typically follow a temperate lake seasonal trajectory, we have observed that community composition in July is very different from the composition observed during spring and fall mixing. Therefore, the mixing disturbance was novel to summer 'July' communities. Finally, the duration of the mixing event presents an extreme disturbance for this ecosystem. After 8 days, complete mixing was achieved in the lake, whereas seasonal spring and autumn mixing occur after weeks to months of gradual cooling of the epilimnion and periodic entrainment of these cooled waters into the hypolimnion, until temperature homogeneity is reached through to the sediment of the lake. The experimental mixing occurred over a much shorter period than seasonal mixing.

Conclusions
Microorganisms as sentinels of change Our results have implications for understanding microbiological responses to environmental perturbations at both fine and coarse temporal scales. In this study, the microbial communities changed after the disturbance, exhibiting sensitivity to a pulse perturbation that occurred over a relatively short time scale (1 week). Further, the clear connection between biogeochemical conditions and community composition suggests that there are direct linkages between changes in environmental forcing, biogeochemical conditions, and microbial communities. The communities recovered within days, and the analysis of variability from 2007 to 2009 demonstrated that the monthly baseline in temporal community variation was relatively stable and may represent a deterministic point of return. Work by our group Shade et al., 2007;Shade et al., 2008) and others (Crump and Hobbie, 2005;Fuhrman et al., 2006;Nelson, 2008;Andersson et al., 2010) has shown repeatable seasonal patterns in microbial communities inhabiting rivers, oceans, seas and lakes, providing cumulative evidence of a consistent microbial community 'baseline' that spans many aquatic systems. In ecology, increases in ecosystem variance have been suggested as indicators of ecological regime shifts (for example, Carpenter and Brock, 2006). Therefore, gradual changes in aquatic microbial communities that exceed normal variability may correspond to press disturbances, such as those expected with gradual global change (for example, increases in mean temperature). Meanwhile, more dramatic microbial shifts may signal pulse disturbances, such as extremes in precipitation or temperature associated with expected increases in climate variability (Folland et al., 2002). The current work and these previous studies together indicate that consistent time series of aquatic microbial communities may serve as robust sentinels over multiple time scales, from local pulse disturbances to global presses.