Peak Abundance of Fatty Acids From Intertidal Biofilm in Relation to the Breeding Migration of Shorebirds

Intertidal biofilm is a thin layer of microbes and meiofauna enmeshed in an extracellular polymeric matrix within and on top of mudflat sediment. This medium provides a dynamic resource for a variety of consumers in estuarine habitats, and is rich in essential fatty acids that birds require for long-distance migration. We measured seasonal changes in biofilm fatty acid content from spring to summer on the Fraser River Estuary, one of the richest and most important ecosystems for migrant and wintering waterbirds in Canada. Fatty acid content in biofilm showed a strong seasonal pattern with a peak in the spring that is associated with the northward migration of Western Sandpipers (Calidris mauri) to their breeding grounds. This peak is linked to the abundance and physiological state, and hence nutritional condition, of epipelic diatoms, which, in turn, depend on a combination of mudflat topography, salinity, temperature, and nutrients that can fluctuate widely with the freshwater inputs from the Fraser River. Specifically, areas with higher elevation (i.e., ∼1 m) had longer periods of exposure to ambient conditions (i.e., light, warm temperatures, gas exchange into/out of biofilms) that facilitated more biofilm growth and higher fatty acid content. Moreover, springtime changes in water chemistry (i.e., salinity/osmotic stresses and nutrients) and temperature facilitated the production of higher overall total lipid/fatty acid contents in the mudflat biofilms compared to summer. Effective conservation of migrating shorebirds depends on the protection of underlying processes at important stopover sites that promote biofilm communities to escalate their production of lipids, including essential fatty acids, during key times of the year.


INTRODUCTION
North American avifauna have suffered a net loss approaching 3 billion birds, or ∼29% relative to 1970 abundance estimates, including steep declines in migratory shorebird species (Rosenberg et al., 2019). The authors attributed this collapse to the "loss of ecosystem integrity, function and services." For migratory shorebirds, a key aspect of ecosystem integrity is the functional quality of the stopover locations that are so vital to the success of their annual migrations. Often these stopover sites are hotspots of biological diversity and productivity and, more often than not, occur in estuarine ecosystems (Butler et al., 2001). Given ongoing declines in bird numbers, it is of paramount importance, from a conservation perspective, to identify and delineate the most critical aspects underpinning the integrity of estuarine ecosystems for migrating shorebirds. Here, we seek to achieve this goal by linking the temporal and spatial dynamics of essential fatty acids, as available in the microphytobenthos, to the needs of migrating shorebirds.
Much of the primary productivity in estuarine environments comes from microphytobenthos (Underwood and Kromkamp, 1999), including biofilm, a thin layer of microbes and meiofauna embedded in an extracellular polymeric matrix that forms on intertidal mudflats. The individual microbial cells grow and reproduce within sediment pores (e.g., clay, silt, and sand) and on the top of the mudflat, and are exposed to ambient conditions during low-tide periods. Mudflat biofilms in the intertidal zone are grazed by marine invertebrates (Herman et al., 2000), fish (Carpentier et al., 2014), and, of key interest here, migratory shorebirds (Elner et al., 2005;Kuwae et al., 2008Kuwae et al., , 2012Mathot et al., 2010;Jardine et al., 2015).
Algal taxa, particularly diatoms, are major sources of essential fatty acids for consumers in aquatic food webs (Kelly and Scheibling, 2012;Quinn et al., 2017). Although many consumers require, to various degrees, specific fatty acids, our focus is on the Western Sandpiper (Calidris mauri). Western Sandpipers extensively graze estuarine intertidal mudflat biofilm (Elner et al., 2005;Kuwae et al., 2008), which serves as a primary source of fatty acids to meet their energetic and physiological demands, particularly during their northward breeding migration each spring (Schnurr et al., 2019). Fats are a preferred migratory fuel because they represent a compact and energy-dense resource compared to carbohydrates and proteins (Guglielmo, 2010). During migration, 90% of the energy required by birds can be derived from lipids, primarily stored as saturated fatty acids (SFA) and monounsaturated fatty acids (MUFA; Guglielmo, 2010), directly or indirectly derived from algae. Further, algae are the main producers (globally) of critically important long-chain polyunsaturated fatty acids (LC-PUFA), including eicosapentaenoic acid (EPA; 20:5n-3) and docosahexaenoic acid (DHA; 22:6n-3) (Hixson et al., 2015). These PUFA are essential to many vertebrates, including birds (Viegas et al., 2017), which benefit from ingesting at least some PUFA (pre-formed) in their diets (Arts et al., 2001).
PUFA are also known to be physiologically important for long distance migration in birds. Consumption of PUFA can increase unsaturation levels of phospholipids in muscle cell membranes (Maillet and Weber, 2007;Weber, 2009), resulting in enhanced membrane fluidity (i.e., decreases in membrane lipid order), which increases permeability, including transmembrane lipid transport (Maillet and Weber, 2006;Weber, 2009). For example, high amounts of EPA and DHA in diets of Semipalmated Sandpipers (Calidris pusilla) enhance oxidative capacity in cellular mitochondria (Maillet and Weber, 2007), positively promoting energy availability required for long distance flights. When muscle cells of Sanderling (Calidris alba) were treated in vitro with n-3 PUFA supplementation, cells increased basal and maximal oxygen consumption (Young, 2019). In addition, PUFA are the precursors of hormones involved in inflammation and cell proliferation, which can facilitate muscle recovery over arduous migration periods (Price, 2010). Therefore, understanding the mechanisms by which estuarine habitats serve to furnish birds with these essential nutrients is essential to effective conservation of migration stopover sites.
Estuarine biofilm communities are dominated by diatoms (Underwood and Paterson, 1993), which can rapidly accumulate fatty acids under specific environmental conditions and growth phases (Schnurr et al., 2019). Under ideal growth conditions diatoms grow and reproduce maximally, causing them to rapidly produce membrane lipids (required for cell division), which consist mostly of phospholipids containing relatively higher proportions of PUFA. When environmental conditions change rapidly (i.e., as a result of an environmental stress), algal cells reduce their rate of cell division and instead store carbon fixed during photosynthesis as triacylglycerol (TAG) -which is comprised of a relatively higher proportion of SFA and MUFA -in lipid droplets located in the cytoplasm (Hu et al., 2008;Yi et al., 2014). For example, marine diatoms in coastal and estuarine ecosystems often experience osmotic stress (Kirst, 1989), and this stress (associated with sudden changes in salinity) induces a lipid/fatty acid accumulation response (Hu and Gao, 2006;Sharma et al., 2012;Mohan and Devi, 2014).
Here, we compared the organic content, and total lipid, chlorophyll-a, and fatty acid contents of intertidal biofilm during spring and summer; two time periods during which shorebirds are migrating northward and southward, respectively, through the Fraser River estuary, British Columbia, Canada. The way in which Western Sandpipers move through the Pacific Flyway differs between the two periods (O'Reilly and Wingfield, 1995;Franks et al., 2014). Northward migration is characterized by huge flocks that rapidly move through to the breeding grounds over narrow time windows (Warnock and Bishop, 1998;Drever et al., 2014). In contrast, southward migration occurs over a protracted period during which birds move through in smaller flocks at a slower pace (Butler et al., 1996). If migrating birds are broadly synchronizing the timing of their northward migration with the peak availability of fatty acids in estuarine biofilm, then seasonal changes in total lipid and fatty acid content should mirror the differences in shorebird use between spring and summer. Given that the realized fatty acid content of estuarine mudflat biofilms is likely a result of complex interactions among a variety of environmental conditions (Schnurr et al., 2019), our aim was to understand how spatial and temporal environmental conditions affect biofilm organic and fatty acid content during these two seasons. Specifically, we investigated how biofilm organic and fatty acid content in biofilms were affected by spatial variables such as mudflat elevation (and how that affects mudflat exposure time), proximity to shore, proximity to the outlet of the Fraser River, and by time-varying factors, including salinity, temperature, nutrient concentrations, and light conditions.

Sampling Dates
Sampling dates were selected to represent the spring (northward) and summer (southward) migration periods in 2017 for Western Sandpipers and Dunlin (Calidris alpina), another common shorebird. Sampling during spring migration was conducted during three campaigns: Campaign 1 occurred April 9-14, Campaign 2 occurred April 21-27, and Campaign 3 occurred May 10-18. Peak northward migration of Western Sandpipers in 2017 occurred around 28 April 2017 (MCD, unpublished data). As the southward migration period occurs over a longer timeframe (i.e., months) than the northward migration (Franks et al., 2014), the two final sampling campaigns were conducted a month apart (August 2-5 and September 5-9) during the summer. These sampling periods occurred during a mix of spring and neap tides (Supplementary Material).

Biofilm Sampling Locations
The Fraser River estuary is located south of the City of Vancouver (49.2827 • N, 123.1207 • W) in southwest British Columbia (Figure 1). It is one of the richest and most important ecosystems for migrant and wintering waterbirds in Canada (Butler and Campbell, 1987). Sampling locations (n = 24) (Figure 1) were determined according to the following criteria: (1) areas that Western Sandpipers are known to utilize for foraging during their stopover periods; and (2) areas across north-south and eastwest transects to understand gradients of proximity to shore and proximity to freshwater river outlets within the estuary. A "spine" of sampling locations was set directly adjacent to provincial government monitoring sites so that we could share their salinity data. Temperature probes were deployed along a north-south transect of sampling locations. East-west transects were perpendicular to the main axis of the "spine" locations to determine whether there was an offshore-onshore spatial gradient in biofilm organic and fatty acid content (Figure 1).

Biofilm Sampling
Biofilm samples (n = 354) were collected during the five campaigns using the same standard procedure. All samples were collected as the mudflat was exposed during ebbing tides, but some were exposed longer (i.e., range of 10-180 min) than others. Upon arrival at each sampling site, a 2.5 × 2.5 m sampling quadrat was set-up on the ground. This quadrat was divided into a sampling grid, whereby the x-and y-axis had 10 equally spaced 25 cm segments. A random number generator was used to determine sampling locations (n = 3) within the quadrat using the segments as guides. Triplicate biofilm samples were taken, whereby a spatula was used to scrape the top ∼2 mm layer of the biofilm at each location, resulting in a ∼1,000 cm 2 sampling area per biofilm sample. These biofilm samples were deposited into 50 ml Falcon TM tubes and immediately stored in a cooler containing dry ice.

Water Sampling for Nutrient Analysis
Water samples were collected for dissolved nutrient analyses by attaching four 100 ml aluminum cans to each stake that marked a sampling location. The aluminum cans were attached ∼12, 25, 35, and 45 cm above the mud surface, so that, upon tidal flooding, each container would fill with water representative of that height and time. Once the tide receded the following day, the cans were removed and contents combined to create a composite sample representative of the flood period on each date at each location. Composite water samples (n = 118) were stored in a 4 • C refrigerator until they were transported to the Pacific Environmental Science Centre, North Vancouver, British Columbia for analysis of macronutrients (soluble reactive phosphorus, nitrate, nitrite, and ammonia). A sub-volume of each sample (n = 118) was sent to the Canada Centre for Inland Waters (CCIW), Burlington, Ontario, for silica analyses.

Temperature and Salinity Data Collection
Temperature and salinity data were continuously collected along the spine (i.e., stations A1-G1; Figure 1). Alpha Mach iBWetLand 22L (model AM015) temperature sensors and data loggers recorded temperature at 15-min intervals, with a resolution of 0.1 • C. Similarly, INW CT2X Conductivity Smart Sensor logged salinity data at 5-min intervals, with a resolution of 0.001 Practical Salinity Units (PSU). Dry zeros were excluded (i.e., when not inundated), and 5th, 50th, and 95th percentiles were calculated for each 24 h and 72-h period around each biofilm sample. Because all these values were strongly inter-correlated (Supplementary Material), we narrowed the list of possible variables by choosing measures related to specific biological hypotheses ( Table 1).

Elevation and Distance Measures
Sample locations were referenced to relevant geographic features using the geospatial tools in Esri ArcMap. Distances of sample locations from the high-water line (HWL) were estimated as the shortest distance from the Canadian Hydrographic Service (CHS) digital HWL to sample location. Distances from the midline of the Fraser River outflow were estimated from the shortest distance between the CHS digital low-water line (that delineated Canoe Passage) to sample location (Figure 1). Elevation measures were obtained as available from LiDAR and Orthophotos taken on the estuary during flights conducted on 23 July 2013. Elevation was determined as an average of the points surrounding each sample location, and corresponds to geodetic elevations in meters, using mean sea level as the zero value.

Organic Content
All frozen (−80 • C) biofilm samples were lyophilized for 2-3 days until completely dry in a freeze-drier (Labconco FreeZone 2.5). Dry biofilm samples were then homogenized using a porcelain mortar and pestle. Ash free dry weight (AFDW) of each sample was determined by taking ∼3 g subsample from the homogenized sample, and ashing them in an oven at 550 • C for 24 h to volatize all organic material, and then re-weighing the samples to determine the difference in mass (i.e., before and after ashing the samples). Organic content was expressed as the percent of the difference from the before-sample.

Chlorophyll-a Content
A subsample of each biofilm sample was used to determine chlorophyll content using a procedure adapted from Arar (1997), and conducted in low lighting conditions to prevent photooxidation of extracted chlorophyll compounds. Depending on the organic fraction of each sample, a specific mass was weighed out to achieve ∼20 mg of organic biomass within the subsample. The chlorophyll within the subsample was then extracted with 10 mL of 9:1 acetone: water in a grinding tube exposed to grinding for 2 min at 500 rpm. Solution and solid biomass contents were transferred from the grinding tube to a centrifugation tube, whereby the sample steeped in a dark refrigerator (4 • C) for 20-23 h. After steeping, samples were centrifuged for 15 min at 675 g, and the supernatant was analyzed using spectrophotometry (Agilent Cary 60 UV-Vis Spectrophotometer) and 1 cm 3 quartz cuvette. A 9:1 acetone: water blank was used to zero the instrument before absorbance measurements were taken. Absorbance values were measured at 750, 664, 647, and 630 nm and applied to the equation [chlorophyll-a (mg/L) = 11.85 × abs664 nm -1.54 × abs647 nm −0.08 × abs630 nm (Arar, 1997)]. Note, abs750 was subtracted from each of the absorbance measurements at the other three wavelengths for chlorophyll because abs750 is for suspended particles that may affect absorbance readings of chlorophyll pigments. Mass of chlorophyll within the extracted sample was determined by multiplying Chlorophyll-a concentration by volume of acetone: water used (10 mL), and that mass was normalized to both the mass of sample and the AFDW, as described above. These Chlorophyll-a values were not adjusted for phaeophytin levels.

Total Lipid and Fatty Acid Analysis
A subsample of each biofilm sample was used to determine total lipid and fatty acid content of mudflat biofilms. As with the chlorophyll sample sizes, each sample mass used equated to ∼20 mg of organic content when considering the AFDW (organic 1 | Explanatory variables for modeling spatial and temporal variation in measures of biofilm organic and fatty acid content.

Total Inorganic Nitrogen
Nitrogen is an essential nutrient in cell growth, homeostasis, and other cellular activities, whose availability can limit growth of marine diatoms.

Soluble Reactive Phosphorus
Phosphates are essential nutrients in cell growth, PUFA accumulation, homeostasis, and other cellular activities, whose availability can limit growth of marine diatoms.

Silica
Silica forms the basis of the frustule exoskeleton of diatoms, and its availability can limit growth of diatoms.

Particle Size Distribution Median
Particle size distribution of the substrate can affect the vertical migration of diatoms and is a proxy for tidal/wave energy, which affects biofilm presence. The median represented the most common particle size.

Salinity (24 h) 50 percentile
Salinity is associated with community composition of the microbial communities found in biofilm.
Temperature (24 h The influence of the Fraser River freshet may vary as a function of distance to the main outflow in the study area.

Distance to High Water Line (HWL)
Tidal/wave energy causes varying shear stresses, thereby affecting accumulation of biomass in mudflat sediments -areas closer to high water line should have lower tidal/wave energies. Nutrient inputs from terrestrial sources should also vary as a function from distance to shore, represented by the high water line.

Elevation
The elevation of the mudflat will determine the amount of time a sampling site is exposed to air, which will determine the amount of light and ambient conditions (i.e., temperature and CO 2 concentrations) the biofilms will be exposed to, thus affecting growth and fatty acid accumulation.
fraction) of the homogenized sample. The Folch extraction method (Folch et al., 1957) was used to extract the lipid from each of the samples. Tricosanoic acid (23:0; Nu-Chek Prep, Inc), which does not commonly occur naturally in this ecosystem, was used to quantify methylation efficiencies during the extraction and methylation procedures. The aqueous layer was discarded during the extraction because all fatty acids were dissolved in the chloroform layer. After extraction, the lipid-chloroform solution was evaporated under a nitrogen blanket before being re-suspended in 2.0 mL of hexane. Two-100 mL aliquots of the lipid-hexanes solution were taken to evaporate in a pre-weighed aluminum tin (∼0.25 mL). The mass difference after evaporation was used to gravimetrically determine the amount of total lipid in the biofilm samples. The remaining lipid-hexanes solution was methylated with a 1% solution of H 2 SO 4 in anhydrous methanol (2 mL), which was heat catalyzed for 90 min at 90 • C (VWR Digital 2 block heater).
After 90 min, a 1 mL aliquot of deionized water was added to samples to stop the reaction, and a 4 mL aliquot of hexanes was added to extract the fatty acids. The solution was then centrifuged at 2,500 rpm for 5 min to separate the hexanes and aqueous layers. The fatty acid-hexanes phase was separated from the aqueous phase with a Pasteur pipette, then dispensed into a clean centrifuge tube. Each sample was then re-extracted two more times with 4 mL hexane aliquots and combined into the sample in the centrifuge tube. Samples were once again evaporated until dry under a nitrogen blanket and re-suspended in 540 µL of hexanes. The final fatty acid-hexanes solution was pipetted into a Gas Chromatograph (GC) vial and stored at −80 • C until analysis.
A Shimadzu GC2010 Plus GC, equipped with an AOC-20i autosampler and a Flame Ionization Detector (FID), was used to identify and quantify the fatty acids in biofilm samples. The GC was fitted with a Supelco SP-2560, 100 m × 0.25 mm column used for fatty acid methyl esters (FAME). Samples were injected in splitless mode. Initial column temperature was 60 • C, which was then increased at 15 • C per minute until 180 • C, followed by 2 • C per minute to 240 • C. The column was held at 240 • C for 5 min to elute all fatty acids potentially remaining in the column. Helium was the carrier gas with a flow rate of 1.2 mL/min. The injector and detector temperatures were both at 250 • C through the duration of all sample runs. Fatty acids were identified by matching retention times from the FAME standard GLC-463 (Nu-Chek Prep, Inc) to the unknown peaks in the samples. Fatty acids were quantified by creating a series of dilutions of the FAME standard mix GLC-463 (Nu-Chek Prep, Inc), producing calibration curves for each of the fatty acids of interest after running them on the GC, then calculating the concentrations of the unknowns from the peak area of the fatty acid within the unknown. Methylation efficiencies were determined by comparing the theoretical amount of Tricosanoic Acid (23:0) added to each sample (10.2 µg) to the amount actually determined to be within each sample. All fatty acid mass fractions were then adjusted for methylation efficiencies, then normalized to both the mass of sample and the AFDW, as described above.

Particle Size Distribution Analysis
Particle size distribution analysis was conducted on samples collected in May, and assumed to have been constant throughout the study period. Analyses were conducted by the Department of Materials Science and Engineering, University of Toronto, based on several grams of lyophilized and homogenized samples. A Horiba Partica LA-950 Laser Diffraction Particle Size Analyzer was used, and summary statistics (median, mean, D10, D90) were derived for each location. D10 is the particle size diameter of the smallest grains (10th percentile of weight distribution), and D90 is the particle size diameter of the largest grains in the sample (90th percentile of weight distribution).

Water Sample Nutrient Analyses
Analysis for nutrients key to diatom growth were determined in the water column above each sampling site during tidal flood periods at the time of biofilm sampling. Before analyses, all samples were filtered with a 0.45 µm cellulose acetate membrane filter. Analyses for nitrite, nitrate and ammonia was colorimetric using Standard Method 4500-NO2-I (nitrite as nitrogen), Standard Method 4500-NO3-I (nitrate as nitrogen), and Standard Method 4500-NH3 H (ammonia as nitrogen), respectively, flow injection analysis, and Lachat Flow Injection System (APHA et al., 2017). Analysis for soluble reactive phosphorus was colorimetric using the Standard Method 4500-P F (phosphorus), a continuous flow analyzer method, and a Technicon AutoAnalyzer II System -segmented flow. Each analyte of interest within a sample was reacted with specific chemicals and analyzed at a wavelength specific to each chemical complex. The protocols for analysis, including the reactant chemicals and wavelengths used for colorimetric analysis are included in each method, which are cited in the References section. Total inorganic nitrogen (TIN) was calculated from a summation of the nitrate, nitrite, and ammonia data collected. Analyses for silica was conducted by Inductively Coupled Plasma-Optical Emission Spectrometry (ICP-OES) (Thermo Scientific; Model ICAP 6300 Duo) using standard solutions produced in the Environment and Climate Change Canada lab at the Canadian Centre for Inland Waters.

Data Selection and Biofilm Variables
We used a series of eight variables to represent the organic, total lipid, and fatty acid content of biofilm. Organic content was expressed as the mass of the organic fraction divided by the mass of the sample weight (g/g). Total lipid and fatty acid content varied both as a fraction of the organic matter within the sample, and by the total amount of organic matter enmeshed in the sediment. Therefore, our metric for total lipid and fatty acid content (that would be available to foraging birds) was normalized per gram of sample, such that: µg fatty acid/g sample = µg fatty acid/g organics × g organics/g sample We summed the fatty acid mass fractions into major groups based on their saturation levels, including total values for SFA, MUFA, and PUFA, including n-3, and n-6 fatty acids. Results are presented on a per gram basis, and summary statistics on measures normalized to AFDW are included in the Supplementary Material.
We considered two different combinations of summed fatty acids to test two hypotheses about (1) physiological condition of diatoms, and (2) community composition of biofilm. First, we recognized that algal cells accumulate TAG when they stop or slow down dividing, e.g., during periods of stress (Sharma et al., 2012). Since cells are dividing at a reduced rate, they are not building/accumulating new membranes, which are primarily composed of PUFA-rich phospholipids. Thus, the ratio of (SFA + MUFA)/PUFA would be elevated during these stress periods, as diatoms would be expected to produce proportionately more TAG relative to phospholipids. A ratio with a value >2.0 would indicate a preferential accumulation of total SFA and MUFA in excess of that which would be predicted if all fatty acid groups were present in the same amounts. Second, we assessed proportions of specific fatty acids relative to total fatty acids as chemotaxonomic biomarkers for microbial community composition (Dalsgaard et al., 2003;Kelly and Scheibling, 2012). We used the following markers: the sum of 16:1n-7 + EPA for diatoms (Shin et al., 2008), the sum of 22:6n-3 (DHA) + 18:1n-9 for dinoflagellates (Kelly and Scheibling, 2012), 18:2n-6 for cyanobacteria, and the sum of odd-chain length fatty acids for bacteria (Kunihiro et al., 2014). We acknowledge the inherent uncertainty in using such biomarkers (i.e., some fatty acids may be found in several taxa within microphytobenthos), and treat any results with the understanding that they only provide broad indices of species turnover.

Selection of Covariates
We examined correlations among three datasets derived from three separate sources (biofilm measures, salinity and temperature values from station devices, and nutrient analysis measures) by calculating Pearson correlation matrices among all variables within each data source (Supplementary Material). From these correlations, we identified a series of response and explanatory variables to identify important spatial and seasonal covariates of biofilm organic and fatty acid content related to known or hypothesized relationships between biofilm, diatoms, and environmental conditions ( Table 1).
To examine how the fatty acid content of the biofilm changed over time, we fitted a mixed effects model to each biofilm measure with sampling date as a categorical explanatory variable. Sampling station was included as a random effect to account for the repeated measures design, and controlled for consistent spatial differences. Biofilm measures were log e -transformed to normalize residuals. We tested for significance of the date term with a Type III Analysis of Variance (i.e., using Type III sums of squares) using the Satterthwaite's method [package lmerTest in R (Kuznetsova et al., 2017)], and a critical alpha level of P = 0.00625 to accommodate for multiple tests on correlated dependent variables. For each model, we checked the residuals for normality and heteroscedasticity. To differentiate whether biofilm measures varied among sampling dates, we conducted Tukey's all-pair comparisons using the glht function in package multcomp in R (Hothorn et al., 2008), and used a compact letter display to reveal groupings of data that were not statistically different.

Correlations With Covariates
We estimated the effect of each covariate individually using a mixed effects modeling approach. A mixed effects model was fitted with each biofilm measure as the response variable, including each covariate as a fixed effect in a model that also included date and sampling station as separate categorical random effects. All covariates were centered and scaled by subtracting the mean and dividing by the standard deviation (i.e., z-score), which allowed inferences about relative effect sizes by directly comparing parameter estimates. Because salinity and temperature data were only available for a subset of the stations their inclusion in a combined model would have resulted in loss of significant quantities of data related to other covariates, and therefore this approach using a series of tests allowed us to make best use of available data.
We tested for significance of each covariate with a t-test using the Satterthwaite's method [package lmerTest in R (Kuznetsova et al., 2017)]. For each model, we checked the residuals for normality and heteroscedasticity, and calculated pseudo-R 2 values (Nakagawa et al., 2017) using package piecewiseSEM (Lefcheck, 2016) in R. Two pseudo-R 2 values were derived for each model: the marginal R 2 considers only the variance explained by the fixed effects, and the conditional R 2 by both the fixed and random effects.

Temporal Trends in Biofilm Measures
The organic content of the sampled sediment ranged between 2 and 3% of the total sample mass, with a tendency toward higher values in August and September (Figure 2). Total lipid varied between a mean of 608 and 719 µg/g biofilm during April and May, respectively, and then declined to 450-237 µg/g biofilm in August and September, respectively. Median Chlorophylla values ranged from 32 to 47 µg/g biofilm with no clear seasonal pattern.
The five most abundant fatty acids found in biofilm during spring/summer sampling included 16:0 (palmitic acid), 16:1n-7 (palmitoleic acid) 20:5n-3 (EPA), 22:3n-3 (docosatrienoic acid), and 22:6n-3 (DHA) ( Table 2). All summed measures of fatty acid content varied widely by date, with highest values observed on April 24 and May 11 (Figure 2; F 4,325 > 15.2, p < 0.0001). For example, SFA had a mean value of 46 µg/g biofilm on 11 April, which doubled to 98 and 104 µg/g biofilm on 24 April and 11 May, respectively. Total SFA then declined to 64 and 44 µg/g biofilm during the August and September sample dates, similar to values observed in early spring (Figure 2). This same pattern of an approximate doubling in measures of fatty acid content from values observed in early spring to late April and May, followed by a decline to summer and early fall period (August/September), was apparent in all the summed measures of fatty acids, including SFA, MUFA, PUFA, n-3, and n-6 fatty acids (Figure 2).
The (SFA + MUFA)/PUFA ratio varied strongly by sampling date (F 4,325 = 25.8, p < 0.0001) ranging from 2.3 to 2.9 in spring to 1.9-2.1 in summer. Confidence intervals for summer values of the ratio overlapped 2.0, indicating that summed SFA and MUFA were not being preferentially accumulated during August and September.
The biomarkers for microbial community composition also showed strongly seasonal variation (Figure 3; F 4,325 = 4.8, p < 0.001). The marker for diatoms varied from 40 to 44% during the spring samples, and then dropped to ∼36% during summer samples. This change was accompanied by concurrent increases in the markers for dinoflagellates, cyanobacteria and bacteria, indicating a shift in the community composition of biofilm between spring and summer.

Correlations With Covariates
Covariates related to nutrients and water conditions varied among sampling dates (Figure 4). TIN was highest (0.25 mg/L) from 11 April to 11 May), after which values remained at ∼0.10 mg/L throughout summer. Similarly, phosphorus declined from 0.021 mg/L in the 11 April to 11 May period, to an average of 0.01 mg/L throughout summer. Both silica and salinity showed clear patterns related to the Fraser River freshet. Silica ranged between 3.2 and 5.1 mg/L and demonstrated little variance over time, except for a strong peak on 11 May, when silica increased to a mean value of 10.8 mg/L. Median salinity from the past 24 h, when averaged across locations, had a mean value of 13.7 PSU in April, which declined over the spring to its lowest mean value of 1.9 PSU during the freshet period in May, and then increased to range between 15−16 PSU during August and September. Median temperature from the past 24 h (i.e., before biofilm sampling) increased from 9.5 • C on 11 April to ∼20 • C in August and September. The 95th percentile of temperature ranged from a mean of 17 • C in April to >22 • C by September.
Biofilm measures were significantly correlated with a suite of spatial attributes and water conditions, of which elevation, salinity (24 h) 50th percentile, temperature (24 h) 50th percentile, and temperature (24 h) 95th percentile were the most influential, with details varying by each of the eight biofilm measures ( Table 3). Elevation also had a strong and positive effect on all biofilm measures (Table 3), indicating the important role played by local topography in determining biofilm organic and fatty acid content. The other spatial covariate, distance to the HWL, was also negatively correlated with organic content, chlorophyll-a, and SFA, indicating higher values closer to shore ( Table 3).
Soluble reactive phosphorus, TIN, and silica showed correlations with biofilm measures. Contrary to expectations, soluble reactive phosphorus and TIN had consistent and negative correlations with measures of fatty acid content ( Table 3). Silica had a negative correlation with organic content and a positive correlation with measures of fatty acid content ( Table 3).
The two measures of temperature had contrasting correlations with biofilm measures: the 95th percentile of temperature over the past 24 h had a negative association with organic content and chlorophyll-a, whereas median temperature was positively correlated with organic content (Table 3). Of the variables related to water conditions, median salinity over the past 24 h had negative correlations with SFA, n-3, and PUFA, and a positive correlation with total lipid ( Table 3).

DISCUSSION
Shorebirds often travel many thousands of kilometers on their annual migrations between wintering and breeding areas. These arduous journeys are only possible if the birds are able to refuel with energy and nutrients thereby maintaining themselves in top physiological condition. To achieve this, birds depend on key stopover sites to provide safe staging places to recover and restock the energy and essential nutrients required for the next leg of their migration (Warnock, 2010). Disruptions to habitat quantity or quality at stopover sites can lead to steep declines in shorebird numbers (Weber et al., 1999;Studds et al., 2017). However, it is not just the location of stopover sites that determines their suitability, there is also a confluence of physical and chemical factors that, operating synergistically, make them ideal places for the generally well-timed production (mainly in the biofilm) of accessible energy and essential nutrients required by shorebirds during their annual migrations (Mathot et al., 2018).
We demonstrate that intertidal mudflat biofilm in the Fraser River estuary, especially in spring when these communities are dominated by diatoms, is a rich source of lipid and essential fatty acids. These riches of energy and nutrients occur around the same time that Western Sandpipers and other shorebirds make their annual northward migration through the area in large congregations (Drever et al., 2014;Schnurr et al., 2019). These seasonal changes in fatty acid content may result from three (non-mutually exclusive) processes: (i) changes in the proportion of diatom biomass within the biofilm which results in more fatty acids; (ii) physiological changes within diatoms already present on the mudflats, resulting in an increase in per cell total lipid and individual fatty acid content, and (iii) seasonal species turnover within the biofilm community. Here, we describe the spatial and temporal confluence of triggers that stimulate lipid production in diatoms, and discuss how each of them may function to promote the observed temporal trajectory of essential fatty acids in these intertidal mudflat biofilm communities.

Spatial Variability
Spatially, we found measures of organic content, total lipid, chlorophyll-a, and fatty acids tended to be higher closer to the HWL and at sites with raised elevation ( Table 3). Elevated areas have longer times (i.e., several hours per tidal period) when the mudflat is exposed to ambient air conditions. These conditions likely allow algal cells to photosynthesize at higher rates since they are exposed to higher photon flux densities, warmer springtime temperatures (Schnurr et al., 2019), and higher rates of gas exchange (i.e., CO 2 into the cells and O 2 out of the biofilms). The sites closer to shore would have longer exposure times (to ambient air conditions) and may experience less shear stress associated with tidal and wave energy (Underwood and Paterson, 1993). Generally, low shear stress facilitates formation and growth of thicker algal biofilms (de Brouwer et al., 2005;Besemer et al., 2007). The positive correlation with elevation might also reflect the presence of intertidal hummocks, which appear to be associated with higher densities of epipelic diatoms (Beninger et al., 2018).

Differences in Diatom Biomass
Spring is a time of rapid change in the Fraser River estuary. Temperature, day length and photon flux densities increase (Schnurr et al., 2019), and the tidal cycles switch such that lowest tides occur during daylight hours rather than overnight (Thomson, 1981). These higher photon flux densities (Jensen and Revsbech, 1989;Schnurr and Allen, 2015;Schnurr et al., 2016b) and temperatures (Blanchard et al., 1997;Kudo et al., 2000;Jiang and Gao, 2004;Scholz and Liebezeit, 2013) result in enhanced rates of photosynthesis and growth in microphytobenthos. Since diatoms are the main primary producers in intertidal biofilm at our study site (Beninger et al., 2011;Schnurr et al., 2019), increases in photosynthetic activity likely cause increases in overall diatom biomass in the biofilm. Thus, we found a positive correlation between organic content and temperature (24 h) 50th percentile (Table 3). Temperatures in the 15−20 • C range foster ideal benthic diatom growth rates (Admiraal, 1976). The spring temperatures experienced on Roberts Bank in 2017 were often within these optimal temperatures (Figure 4), and outside temperature ranges that can result in growth inhibition and/or cell death (Thompson and Guo, 1992;Kudo et al., 2000;Jiang and Gao, 2004;Woelfel et al., 2014). Although not measured directly, the effect of photon flux density is encapsulated by the strong effect of date on the measures of fatty acid abundance. In addition, the diatom fatty acid biomarkers were higher in the spring compared to summer (Figure 3), consistent with the hypothesis that higher fatty acid content of biofilm results from an overall higher proportion of diatom biomass.
We found that fatty acid content was significantly reduced during summer, and this was likely a result of temperatures that were too high for diatoms to grow and reproduce. Indeed, we regularly observed maximal daytime temperatures on mudflat biofilms in the summertime to be greater than 33 • C, which is above optimal temperatures for diatom growth (Thompson and Guo, 1992;Kudo et al., 2000;Jiang and Gao, 2004;Woelfel et al., 2014). Higher temperatures during summer would have also resulted in a stronger "de-watering" effect, wherein changes in the sediment bulk density can result in lower measures of Chl a content of microphytobenthos (Perkins et al., 2003).

Physiological Changes
Differences in fatty acid content of biofilm between spring and summer may also occur from physiological changes within the diatom communities themselves. The existence of a lipid accumulation response during spring was supported by the larger (SFA + MUFA)/PUFA ratios that ranged from 2.3 to 2.9, indicating a preferential production of triacylglycerol relative to summer samples. Laboratory studies have shown lipid accumulation in microalgae occurs in two stages (Rodolfi et al., 2009). First, when all requirements for growth are met, microalgae divide rapidly and synthesize (mostly) phospholipiddominated membrane lipids, which are predominately composed of MUFA and PUFA (Solovchenko, 2012). Second, when any growth factor becomes limiting, cells enter a lipid-producing phase characterized by a slowdown or cessation of cell division accompanied by an accumulation of triacylglycerol, predominantly composed of SFA and MUFA (Sharma et al., 2012;Solovchenko, 2012). Fatty acid content during the spring sampling dates showed a wide variance, with some locations having very high values relative to the mean (Figure 2). Total fatty acid values had a double mode distribution, and these very high values could account for up to 85% of the total fatty acid content in the spring samples (Supplementary Material). Presently we have no way to differentiate the relative roles of increases in diatom biomass and lipid accumulation response in determining the total amount of fatty acids in intertidal biofilm, and we suggest that these local sites with very high values result from their combined action, i.e., very high local bursts of fatty acids become available when diatoms that are experiencing good growing conditions are suddenly triggered to accumulate lipid.
Sudden shifts in water chemistry such as salinity or pH can cause lipid accumulation responses in microalgae (Sharma et al., 2012). On Roberts Bank, the daily discharge of the Fraser River increases steadily during April and May on account of the spring freshet fueled by snow melt (Supplementary Material), and this influx of freshwater into the tidal estuary results in rapid changes in salinity and water chemistry. We consider rapid changes in salinity as the strongest driver in the lipid FIGURE 2 | Temporal changes in biofilm measures on Roberts Banks mudflats located on the Fraser River estuary, British Columbia, April-September 2017. Biofilm measures are expressed in units of per 1 g of sample, including inedible sediment. Note that y-axis extends from 0 to 0.90 percentile of each measure to avoid compression of scale from outliers. Box plots represent the distribution of observed values, where midline is the median, with the upper and lower limits of the box being 75th and 25th percentiles. Whiskers extend up to 1.5 × the interquartile range, and outliers are depicted as points. Blue circles indicate predicted means from linear mixed effects models, and bounds are 95% predictions intervals from fixed effects. Letters indicate groupings that are not statistically significantly different according to Tukey's pairwise comparisons.
FIGURE 3 | Temporal changes in the ratio of fatty acid groups and chemotaxonomic markers on Roberts Banks mudflats located on the Fraser River estuary, British Columbia, April-September 2017. Note that y-axis for biomarker plots extend from 0.025 to 0.975 percentiles of each measure to avoid compression of scale from outliers. Box plots represent the distribution of observed values, where midline is the median, with the upper and lower limits of the box being 75th and 25th percentiles. Whiskers extend up to 1.5 × the interquartile range, and outliers are depicted as points. Blue circles indicate predicted means from linear mixed effects models, and bounds are 95% prediction intervals from fixed effects. Letters indicate groupings that are not statistically significantly different according to Tukey's pairwise comparisons.
accumulation response because the rapid drop in salinity in May coincided with the peak abundance of fatty acids (Figure 4), and the covariate modeling ( Table 3) indicated salinity had a stronger role on fatty acid content than temperature. Further, seasonal changes are accompanied by wide daily fluctuations related to tidal cycles (i.e., salinity wedges), whereby dense saltwater flooding the mudflats causes a spike in salinity (Correll, 1978) until the tide peaks and then begins to recede, when diffusion and currents cause the salinity to decrease. This twicedaily event would contribute to a lipid accumulation response in response to the osmotic/salinity stress. In the spring, fatty acid content was highest during the May samples, which occurred during spring tides (Supplementary Material) that, in addition to resulting in longer exposure times, would have accentuated daily variance in salinity.
A lipid accumulation response can be also induced by factors such as changes in nutrients, light levels, and temperature (Sharma et al., 2012;Schnurr et al., 2016a). Diatoms can increase their total lipid (including fatty acids) contents (up to threefold) when "starved" of nitrogen and/or silica for relatively short periods of time (i.e., 4−120 h) (Opute, 1974;Shifrin and Chisholm, 1981;Roessler, 1988;Sriharan et al., 1991;FIGURE 4 | Temporal changes in nutrients and water conditions on Roberts Bank, Fraser River Delta, British Columbia, 2017. Blue circles indicate predicted means, and bounds are 95% prediction intervals from linear mixed effects models. Variables plotted are Total Inorganic Nitrogen, Soluble Reactive Phosphorus, Silica, Salinity (24 h) 50 percentile, Temperature (24 h) 50th percentile, and Temperature (24 h) 95th percentile. Beta coefficients are derived from mixed effects models (see text for details), and sorted by absolute values from largest to smallest. P-values shown as 0.00 are <0.001.
Marginal R 2 considers only the variance explained by the fixed effects, and the conditional R 2 by both the fixed and random effects (the total variance explained by the model). Only covariates that met a cut-off of p < 0.1 are included. Hu et al., 2008). Although the negative correlations between fatty acid content with N and P ( Table 3) are consistent with lipid accumulation response, these nutrients on Roberts Bank changed gradually over time (Figure 4), and therefore we interpret this correlation as a seasonal depletion, and that N and P are not limiting in a nutrient-rich marine environment. Silica had a negative correlation with organic content, and a positive correlation with measures of fatty acid content (Table 3). Silica (driven by the freshet) is likely fueling the diatom productivity selectively, wherein diatoms require this nutrient to form their exterior frustules, in contrast to other taxa in the microphytobenthos.

Community Composition
The composition of algal communities can have a pronounced effect on their fatty acid profiles (Galloway and Winder, 2015;Schnurr and Allen, 2015). Thus, differences in fatty acid content between spring and summer may have also resulted from changes in community composition of the intertidal biofilm. The Roberts Bank diatom community is dominated by Navicula spp. and Nitzschia spp. (Beninger et al., 2011;Schnurr et al., 2019). Based on fatty acid markers, we suggest that the overall diatom biomass within the biofilm declines from spring to summer (Figure 3), a change accompanied by an increase in bacteria, cyanobacteria, and dinoflagellates. The highest temperatures (>30 • C) observed in the summertime would inhibit diatom growth and reproduction (Renaud et al., 1994;Kudo et al., 2000;Jiang and Gao, 2004;Scholz and Liebezeit, 2013;Woelfel et al., 2014), but favor the growth of cyanobacteria, chlorophytes and bacteria (Schnurr and Allen, 2015). These changes tie the fatty acid content of biofilm to the presence of epipelic diatoms, with the caveat that the effect of species turnover on fatty acid content of intertidal biofilm needs to be better explored with more detailed taxonomic data.

Implications for Shorebird Migration
Hitherto, there has been only a limited understanding of why shorebirds choose where, when and how long to stop and forage during long-distance migration. For example, why Western Sandpipers stop at only a few large estuarine mudflats during their northward breeding migration along the Pacific Flyway (Mathot et al., 2018), but are dissipated over multiple smaller sites on their return to non-breeding grounds, has been enigmatic (O'Reilly and Wingfield, 1995). Conventional ornithological hypotheses to explain this dichotomy in migration strategies are based on either access to energy (Farmer and Wiens, 1999) or predation threats (Lank et al., 2003). Our finding that the northward arrival of Western Sandpipers and other shorebirds on Roberts Bank mirrors seasonal increases in total lipid and in the abundance of essential fatty acids found in intertidal biofilm provides a new explanation. The discovery complements previous studies demonstrating that shorebirds alter their foraging behaviors to select areas with high microphytobenthic biomass Jiménez et al., 2015) and the activities of shorebirds can significantly affect benthic diatom growth (Jauffrais et al., 2015). The inference is that migration success for shorebirds hinges on their presence at a stopover site that generally coincides with maximal total lipid and fatty acid production at the site. In the case of Roberts Bank, the main trigger for benthic diatoms switching to their lipid/fatty acid accumulation phase appears to be stress associated with high amplitude oscillations in salinity. We speculate that equivalent mechanisms triggering lipid/fatty acid accumulation responses in diatom-dominated biofilm assemblages are likely also affecting the quality of migratory stopover locations on mudflats elsewhere on the Pacific and along other international flyways. Investigating the existence of and, if so, elucidating the nature of stress triggers at all these sites should be a research priority. Of particular importance is discerning whether these triggers are driven by similar salinity-associated circumstances to Roberts Bank or an alternative mechanism. Meanwhile, maintaining the salinity trigger is an imperative for the effective conservation of migratory shorebirds in the Fraser River estuary and, if a pattern emerges on other stopover sites, this new understanding would galvanize a paradigm shift in shorebird conservation worldwide. Springtime accumulation of biofilm biomass and fatty acids represents a crucial "pulse" of fuel and essential nutrients into estuarine food webs. Invertebrate populations foraging on primary producers in the microphytobenthos can show rapid increases in abundance (Sahan et al., 2007), and thus accumulation of lipid and fatty acids up the food web (i.e., into fish, birds, etc.). Given the important role of essential fatty acids in the growth and reproduction of a wide variety of organisms within aquatic food webs (Kainz et al., 2004), understanding the physiological response of algae/diatoms to complex and interacting environmental factors in estuarine mudflats is fundamental to designing strategies to conserve and protect migrating birds and estuarine ecosystems in their entirety. We found that total lipid and fatty acid content in biofilm on mudflats varied between spring and summer and that these nutritionally important compounds were significantly affected by the outflow of freshwater from the Fraser River. To the north of our study site lies Sturgeon Bank, another large estuarine mudflat that has experienced extensive alterations to its coastal oceanography and sediment transport (Atkins et al., 2016). The Sturgeon Bank area now sees reduced use by shorebirds during the critical northward migration period relative to Roberts Bank (Jardine et al., 2015). Such differences underscores the need to protect these mudflat habitats in as pristine a state as possible so they may continue to furnish shorebirds with the energy and essential fatty acids required for long-distance migration.

DATA AVAILABILITY STATEMENT
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.