Allochthony, fatty acid and mercury trends in muscle of Eurasian perch (Perca fluviatilis) along boreal environmental gradients.

Environmental change, including joint effects of increasing dissolved organic carbon (DOC) and total phosphorus (TP) in boreal northern lakes may affect food web energy sources and the biochemical composition of organisms. These environmental stressors are enhanced by anthropogenic land-use and can decrease the quality of polyunsaturated fatty acids (PUFAs) in seston and zooplankton, and therefore, possibly cascading up to fish. In contrast, the content of mercury in fish increases with lake browning potentially amplified by intensive forestry practises. However, there is little evidence on how these environmental stressors simultaneously impact beneficial omega-3 fatty acid (n3-FA) and total mercury (THg) content of fish muscle for human consumption. A space-for-time substitution study was conducted to assess whether environmental stressors affect Eurasian perch (Perca fluviatilis) allochthony and muscle nutritional quality [PUFA, THg, and their derivative, the hazard quotient (HQ)]. Perch samples were collected from 31 Finnish lakes along pronounced lake size (0.03-107.5 km2), DOC (5.0-24.3 mg L-1), TP (5-118 µg L-1) and land-use gradients (forest: 50.7-96.4%, agriculture: 0-32.6%). These environmental gradients were combined using principal component analysis (PCA). Allochthony for individual perch was modelled using source and consumer δ2H values. Perch allochthony increased with decreasing lake pH and increasing forest coverage (PC1), but no correlation between lake DOC and perch allochthony was found. Perch muscle THg and omega-6 fatty acid (n6-FA) content increased with PC1 parallel with allochthony. Perch muscle DHA (22:6n3) content decreased, and ALA (18:3n3) increased towards shallower murkier lakes (PC2). Perch allochthony was positively correlated with muscle THg and n6-FA content, but did not correlate with n3-FA content. Hence, the quality of perch muscle for human consumption decreases (increase in HQ) with increasing forest coverage and decreasing pH, potentially mediated by increasing fish allochthony.

Role of potentially toxic cyanobacteria in crustacean zooplankton diet in a eutrophic lake.

The coexistence of potentially toxic bloom-forming cyanobacteria (CY) and generally smaller-sized grazer communities has raised the question of zooplankton (ZP) ability to control harmful cyanobacterial blooms and highlighted the need for species-specific research on ZP-CY trophic interactions in naturally occurring communities. A combination of HPLC, molecular and stable isotope analyses was used to assess in situ the importance of CY as a food source for dominant crustacean ZP species and to quantify the grazing on potentially toxic strains of Microcystis during bloom formation in large eutrophic Lake Peipsi (Estonia). Aphanizomenon, Dolichospermum, Gloeotrichia and Microcystis dominated bloom-forming CY, while Microcystis was the major genus producing cyanotoxins all over the lake. Grazing studies showed that CY, and especially colonial CY, formed a significant, and also preferred component of algae ingested by the cladocerans Bosmina spp. and Daphnia spp. while this was not the case for the more selective calanoid copepod Eudiaptomus gracilis. Molecular analyses confirmed the presence of CY, including Microcystis, in ZP guts. Further analyses using qPCR targeting cyanobacterial genus-specific mcyE synthase genes indicated that potentially toxic strains of Microcystis can be ingested directly or indirectly by all the dominant crustacean grazers. However, stable isotope analyses indicated that little, if any, assimilation from ingested bloom-forming CY occurred. The study suggests that CY, and particularly Microcystis with both potentially toxic and non-toxic strains, can be widely ingested by cladoceran grazers during a bloom event with implications for control of CY abundance and for transfer of CY toxins through the food web.

Influence of littoral periphyton on whole-lake metabolism relates to littoral vegetation in humic lakes.

The role of littoral habitats in lake metabolism has been underrated, especially in humic lakes, based on an assumption of low benthic primary production (PP) due to low light penetration into water. This assumption has been challenged by recent recognition of littoral epiphyton dominance of whole-lake PP in a small highly humic lake and of epiphyton as an important basal food source for humic lake biota. However, as these studies have mostly concerned single lakes, there is a need to test their wider generality. We studied the whole-lake PP and community respiration (CR) in eight small humic lakes in southern Finland during July 2015 using 14 C incorporation to measure pelagic PP and the changes in dissolved inorganic carbon in light and dark in situ incubations to measure CR and littoral PP by epiphyton. Changes in O2 concentration in both pelagic and littoral surface water were measured periodically from each lake and, additionally, continuously with a data logger from one lake during the study period. The results revealed that the littoral dominated whole-lake net primary production (NPP) in five of the eight lakes, which was supported by observed O2 supersaturation in the littoral surface water in most of the lakes. Calculated pelagic:littoral ratios by area correlated negatively with both littoral NPP and littoral contribution to whole-lake NPP. Moreover, there was a significant positive relationship between littoral proportion of whole-lake NPP and the fraction of lake surface area covered by littoral aquatic vegetation. This demonstrates that increased aquatic littoral vegetation cover increases the overall importance of the littoral to whole-lake PP in highly humic lakes. Littoral NPP also correlated strongly with littoral O2 saturation, and the continuously measured O2 revealed substantial temporal variation in O2 saturation, particularly in the littoral zone. Whole-lake gross primary production:community respiration (GPP:CR) ratios revealed that accounting for littoral metabolism produced a marked shift towards lake metabolic balance, although all the eight lakes remained net heterotrophic. This study emphasizes that littoral metabolism needs to be accounted for when estimating whole-lake C fluxes in all lakes, even in highly colored humic waters.

Oxygen and carbon isoscapes for the Baltic Sea: Testing their applicability in fish migration studies.

Conventional tags applied to individuals have been used to investigate animal movement, but these methods require tagged individuals be recaptured. Maps of regional isotopic variability known as "isoscapes" offer potential for various applications in migration research without tagging wherein isotope values of tissues are compared to environmental isotope values. In this study, we present the spatial variability in oxygen (δ18OH2O) and dissolved inorganic carbon (δ13 CDIC) isotope values of Baltic Sea water. We also provide an example of how these isoscapes can reveal locations of individual animal via spatial probability surface maps, using the high-resolution salmon otolith isotope data from salmon during their sea-feeding phase in the Baltic Sea. A clear latitudinal and vertical gradient was found for both δ18OH2O and δ13 CDIC values. The difference between summer and winter in the Baltic Sea δ18OH2O values was only slight, whereas δ13 CDIC values exhibited substantial seasonal variability related to algal productivity. Salmon otolith δ18Ooto and δ13Coto values showed clear differences between feeding areas and seasons. Our example demonstrates that dual isotope approach offers great potential for estimating probable fish habitats once issues in model parameterization have been resolved.

Land Use Affects Carbon Sources to the Pelagic Food Web in a Small Boreal Lake.

Small humic forest lakes often have high contributions of methane-derived carbon in their food webs but little is known about the temporal stability of this carbon pathway and how it responds to environmental changes on longer time scales. We reconstructed past variations in the contribution of methanogenic carbon in the pelagic food web of a small boreal lake in Finland by analyzing the stable carbon isotopic composition (δ13C values) of chitinous fossils of planktivorous invertebrates in sediments from the lake. The δ13C values of zooplankton remains show several marked shifts (approx. 10 ‰), consistent with changes in the proportional contribution of carbon from methane-oxidizing bacteria in zooplankton diets. The results indicate that the lake only recently (1950s) obtained its present state with a high contribution of methanogenic carbon to the pelagic food web. A comparison with historical and palaeobotanical evidence indicates that this most recent shift coincided with agricultural land-use changes and forestation of the lake catchment and implies that earlier shifts may also have been related to changes in forest and land use. Our study demonstrates the sensitivity of the carbon cycle in small forest lakes to external forcing and that the effects of past changes in local land use on lacustrine carbon cycling have to be taken into account when defining environmental and ecological reference conditions in boreal headwater lakes.

Top consumer abundance influences lake methane efflux.

Lakes are important habitats for biogeochemical cycling of carbon. The organization and structure of aquatic communities influences the biogeochemical interactions between lakes and the atmosphere. Understanding how trophic structure regulates ecosystem functions and influences greenhouse gas efflux from lakes is critical to understanding global carbon cycling and climate change. With a whole-lake experiment in which a previously fishless lake was divided into two treatment basins where fish abundance was manipulated, we show how a trophic cascade from fish to microbes affects methane efflux to the atmosphere. Here, fish exert high grazing pressure and remove nearly all zooplankton. This reduction in zooplankton density increases the abundance of methanotrophic bacteria, which in turn reduce CH4 efflux rates by roughly 10 times. Given that globally there are millions of lakes emitting methane, an important greenhouse gas, our findings that aquatic trophic interactions significantly influence the biogeochemical cycle of methane has important implications.

Inferring Phytoplankton, Terrestrial Plant and Bacteria Bulk δ¹³C Values from Compound Specific Analyses of Lipids and Fatty Acids.

Stable isotope mixing models in aquatic ecology require δ13C values for food web end members such as phytoplankton and bacteria, however it is rarely possible to measure these directly. Hence there is a critical need for improved methods for estimating the δ13C ratios of phytoplankton, bacteria and terrestrial detritus from within mixed seston. We determined the δ13C values of lipids, phospholipids and biomarker fatty acids and used these to calculate isotopic differences compared to the whole-cell δ13C values for eight phytoplankton classes, five bacterial taxa, and three types of terrestrial organic matter (two trees and one grass). The lipid content was higher amongst the phytoplankton (9.5±4.0%) than bacteria (7.3±0.8%) or terrestrial matter (3.9±1.7%). Our measurements revealed that the δ13C values of lipids followed phylogenetic classification among phytoplankton (78.2% of variance was explained by class), bacteria and terrestrial matter, and there was a strong correlation between the δ13C values of total lipids, phospholipids and individual fatty acids. Amongst the phytoplankton, the isotopic difference between biomarker fatty acids and bulk biomass averaged -10.7±1.1‰ for Chlorophyceae and Cyanophyceae, and -6.1±1.7‰ for Cryptophyceae, Chrysophyceae and Diatomophyceae. For heterotrophic bacteria and for type I and type II methane-oxidizing bacteria our results showed a -1.3±1.3‰, -8.0±4.4‰, and -3.4±1.4‰ δ13C difference, respectively, between biomarker fatty acids and bulk biomass. For terrestrial matter the isotopic difference averaged -6.6±1.2‰. Based on these results, the δ13C values of total lipids and biomarker fatty acids can be used to determine the δ13C values of bulk phytoplankton, bacteria or terrestrial matter with ± 1.4‰ uncertainty (i.e., the pooled SD of the isotopic difference for all samples). We conclude that when compound-specific stable isotope analyses become more widely available, the determination of δ13C values for selected biomarker fatty acids coupled with established isotopic differences, offers a promising way to determine taxa-specific bulk δ13C values for the phytoplankton, bacteria, and terrestrial detritus embedded within mixed seston.

Lake size and fish diversity determine resource use and trophic position of a top predator in high-latitude lakes.

Prey preference of top predators and energy flow across habitat boundaries are of fundamental importance for structure and function of aquatic and terrestrial ecosystems, as they may have strong effects on production, species diversity, and food-web stability. In lakes, littoral and pelagic food-web compartments are typically coupled and controlled by generalist fish top predators. However, the extent and determinants of such coupling remains a topical area of ecological research and is largely unknown in oligotrophic high-latitude lakes. We analyzed food-web structure and resource use by a generalist top predator, the Arctic charr Salvelinus alpinus (L.), in 17 oligotrophic subarctic lakes covering a marked gradient in size (0.5-1084 km(2)) and fish species richness (2-13 species). We expected top predators to shift from littoral to pelagic energy sources with increasing lake size, as the availability of pelagic prey resources and the competition for littoral prey are both likely to be higher in large lakes with multispecies fish communities. We also expected top predators to occupy a higher trophic position in lakes with greater fish species richness due to potential substitution of intermediate consumers (prey fish) and increased piscivory by top predators. Based on stable carbon and nitrogen isotope analyses, the mean reliance of Arctic charr on littoral energy sources showed a significant negative relationship with lake surface area, whereas the mean trophic position of Arctic charr, reflecting the lake food-chain length, increased with fish species richness. These results were supported by stomach contents data demonstrating a shift of Arctic charr from an invertebrate-dominated diet to piscivory on pelagic fish. Our study highlights that, because they determine the main energy source (littoral vs. pelagic) and the trophic position of generalist top predators, ecosystem size and fish diversity are particularly important factors influencing function and structure of food webs in high-latitude lakes.

Resistant microbial cooccurrence patterns inferred by network topology.

Although complex cooccurrence patterns have been described for microbes in natural communities, these patterns have scarcely been interpreted in the context of ecosystem functioning and stability. Here we constructed networks from species cooccurrences between pairs of microorganisms which were extracted from five individual aquatic time series, including a dystrophic and a eutrophic lake as well as an open ocean site. The resulting networks exhibited higher clustering coefficients, shorter path lengths, and higher average node degrees and levels of betweenness than those of random networks. Moreover, simulations demonstrated that taxa with a large number of cooccurrences and placement at convergence positions in the network, so-called "hubs" and "bottlenecks," confer resistance against random removal of "taxa." Accordingly, we refer to cooccurrences at convergence positions as system-relevant interdependencies, as they, like hubs and bottlenecks, determine network topology. These topology features of the cooccurrence networks point toward microbial community dynamics being resistant over time and thus could provide indicators for the state of ecosystem stability.

An empirical evaluation of the utility of convex hull and standard ellipse areas for assessing population niche widths from stable isotope data.

Stable isotope analyses are increasingly employed to characterise population niche widths. The convex hull area (TA) in a δ(13)C-δ(15)N biplot has been used as a measure of isotopic niche width, but concerns exist over its dependence on sample size and associated difficulties in among-population comparisons. Recently a more robust method was proposed for estimating and comparing isotopic niche widths using standard ellipse areas (SEA), but this approach has yet to be tested with empirical stable isotope data. The two methods measure different kind of isotopic niche areas, but both are now widely used to characterise isotopic niche widths of populations. We used simulated data and an extensive empirical dataset from two fish populations to test the influence of sample size on the observed isotopic niche widths (TA and SEA). We resampled the original datasets to generate 5000 new samples for different numbers of observations from 5 to 80 to examine the statistical distributions of niche area estimates for increasing sample size. Our results illustrate how increasing sample size increased the observed TA; even sample sizes much higher than n = 30 did not improve the precision for the TA method. SEA was less sensitive to sample size, but the natural variation in our empirical fish δ(13)C and δ(15)N data still resulted in considerable uncertainty around the mean estimates of niche width, reducing the precision particularly with sample sizes n<30. These results confirm that the TA method is less appropriate for estimating population isotopic niche areas using small samples, especially when considerable population level isotope variation is expected. The results also indicate a need for caution when using SEA as a measure of trophic niche widths for consumers, particularly with low sample sizes and when the distribution and range for population isotope values are not known.

Bacterial and phytoplankton responses to nutrient amendments in a boreal lake differ according to season and to taxonomic resolution.

Nutrient limitation and resource competition in bacterial and phytoplankton communities may appear different when considering different levels of taxonomic resolution. Nutrient amendment experiments conducted in a boreal lake on three occasions during one open water season revealed complex responses in overall bacterioplankton and phytoplankton abundance and biovolume. In general, bacteria were dominant in spring, while phytoplankton was clearly the predominant group in autumn. Seasonal differences in the community composition of bacteria and phytoplankton were mainly related to changes in observed taxa, while the differences across nutrient treatments within an experiment were due to changes in relative contributions of certain higher- and lower-level phylogenetic groups. Of the main bacterioplankton phyla, only Actinobacteria had a treatment response that was visible even at the phylum level throughout the season. With increasing resolution (from 75 to 99% sequence similarity) major responses to nutrient amendments appeared using 454 pyrosequencing data of 16S rRNA amplicons. This further revealed that OTUs (defined by 97% sequence similarity) annotated to the same highly resolved freshwater groups appeared to occur during different seasons and were showing treatment-dependent differentiation, indicating that OTUs within these groups were not ecologically coherent. Similarly, phytoplankton species from the same genera responded differently to nutrient amendments even though biovolumes of the majority of taxa increased when both nitrogen and phosphorus were added simultaneously. The bacterioplankton and phytoplankton community compositions showed concurrent trajectories that could be seen in synchronous succession patterns over the season. Overall, our data revealed that the response of both communities to nutrient changes was highly dependent on season and that contradictory results may be obtained when using different taxonomic resolutions.

Distinct and diverse anaerobic bacterial communities in boreal lakes dominated by candidate division OD1.

Lakes have a central role in the carbon cycle of the boreal landscape. These systems typically stratify in summer and their hypolimnetic microbial communities influence burial of biogenic organic matter in sediments. The composition of bacterial communities in these suboxic habitats was studied by pyrosequencing of 16S rRNA amplicons from five lakes with variable dissolved organic carbon (DOC) concentrations. Bacterioplankton communities in the hypolimnetic waters were clearly different from the surface layer with candidate division OD1, Chlorobi and Bacteroidetes as dominant community members. Several operational taxonomic units (OTUs) affiliated with candidate division OD1 were abundant and consistently present in the suboxic hypolimnion in these boreal lakes. The overall representation of this group was positively correlated with DOC and methane concentrations. Network analysis of time-series data revealed contrasting temporal patterns but suggested similar ecological roles among the abundant OTUs affiliated with candidate division OD1. Together, stable isotope data and taxonomic classification point to methane oxidation and autotrophic denitrification as important processes in the suboxic zone of boreal lakes. Our data revealed that while hypolimnetic bacterial communities are less dynamic, they appear to be more diverse than communities from the oxic surface layer. An appreciable proportion of the hypolimnetic bacteria belong to poorly described phyla.

Widespread contribution of methane-cycle bacteria to the diets of lake profundal chironomid larvae.

Reports of unexpectedly 13C-depleted chironomid larvae in lakes have led to an hypothesis that significant transfer of detrital organic matter to chironomid larvae may occur via methane-cycle bacteria. However, to date little is known of how such transfer might vary across species and lakes. We gathered data from 87 lakes to determine how widespread this phenomenon might be and to define boundaries for its likely magnitude. Carbon stable isotope values of chironomid larvae varied greatly between taxa. Very marked 13C-depletion was evident only in certain taxa, especially Chironomus plumosus, C. anthracinus, and C. tenuistylus, all characteristic of eutrophic or dystrophic lakes and known to be tolerant of low oxygen conditions. Furthermore, marked 13C-depletion was only found in larvae from lakes in which late-summer hypolimnetic oxygen depletion near the sediment surface was below an apparent threshold concentration of 2-4 mg O2/L. Similarly, application of a two-source mixing model suggested that methanotrophic bacteria made the greatest contribution to profundal chironomid growth (0-70% of larval carbon) when the late-summer oxygen concentration dropped below approximately 2 mg O2/L. Our study demonstrates that methane-derived carbon is an important, but often neglected, contribution to the flux of carbon through the food webs of many productive or dystrophic lakes.

Whole-lake dissolved inorganic 13C additions reveal seasonal shifts in zooplankton diet.

Sustained whole-lake additions of 13C-enriched dissolved inorganic carbon (DIC), intended to increase experimentally the delta13C of DIC in the epilimnion of a small lake with high dissolved organic carbon (DOC), were made during three seasonal periods (spring, summer, and autumn). Coupled with carbon and nitrogen stable isotope analysis of zooplankton and several of their putative food sources, these additions were used to investigate seasonal changes in the relative contributions of different food sources to zooplankton diet in the lake. Four main potential food sources were considered: phytoplankton, heterotrophic bacteria (HB), methanotrophic bacteria (MOB), and green sulfur bacteria (GSB). Because the number of potential food sources exceeded the number of isotopes analyzed, a computer program (IsoSource) was used to estimate the range of possible contributions of the various food sources. During all three periods the added inorganic 13C quickly increased the epilimnetic DIC delta13C by between 18 per thousand and 21 per thousand above the initial value of approximately -21 per thousand. This 13C enrichment of DIC was rapidly transmitted to the particulate organic matter (POM), which included photosynthetic phytoplankton. In spring and summer, delta13C of both adult and juvenile Daphnia increased by approximately 10 per thousand, indicating that Daphnia utilized autochthonous carbon. However, this 13C labeling of Daphnia was not so obvious during the autumn period, when their delta13C generally decreased. According to the IsoSource model outputs based on both delta13C and delta15N values, Daphnia utilized all four potential food source types during spring, summer, and autumn, but in different proportions. The possible contribution of phytoplankton to Daphnia diet was substantial (25-71%) in all seasons. The possible contributions of the bacterial food sources were more variable. The possible contribution of GSB was minor (0-20%) at all times and negligible in autumn. The possible contribution of HB was higher but very variable. Methanotrophic bacteria always made a significant contribution to Daphnia diet and were likely the single most important food source in autumn. Since both HB and MOB in this high-DOC lake probably depend largely on allochthonous organic carbon, our results highlight the seasonal variability in the potential importance of ecosystem subsidies in lake food webs.

Human dietary intake of organochlorines from Baltic herring: implications of individual fish variability and fisheries management.

This study examines the extent to which Finnish human dietary intake of organochlorines (PCDD/Fs and PCBs) originating from Northern Baltic herring can be influenced by fisheries management. This was investigated by estimation of human intake using versatile modeling tools (e.g., a herring population model and a bioenergetics model). We used a probabilistic approach to account for the variation in human intake of organochlorines originating from the variation among herring individuals. Our estimates were compared with present precautionary limits and recommendation for use. The results show that present consumption levels and frequencies of herring give a high probability of exceeding recommended intake limits of PCDD/Fs and PCBs. Furthermore, our results clearly demonstrate that in the risk management of dioxinlike organochlorines, regulating fishing (in this case increasing fishing pressure) is a far less effective way to decrease the risk than regulating the consumption of herring. Increased fishing would only slightly decrease organochlorine concentrations of herring in the Finnish fish market.

Terrestrial carbon and intraspecific size-variation shape lake ecosystems.

Conceptual models of lake ecosystem structure and function have generally assumed that energy in pelagic systems is derived from in situ photosynthesis and that its use by higher trophic levels depends on the average properties of individuals in consumer populations. These views are challenged by evidence that allochthonous subsidies of organic carbon greatly influence energy mobilization and transfer and the trophic structure of pelagic food webs, and that size variation within consumer species has major ramifications for lake community dynamics and structure. These discoveries represent conceptual shifts that have yet to be integrated into current views on lake ecosystems. Here, we assess key aspects of energy mobilization and size-structured community dynamics, and show how these processes are intertwined in pelagic food webs.

The contribution of insect prey to the total nitrogen content of sundews (Drosera spp.) determined in situ by stable isotope analysis.

• The contribution of insect prey to total N in the carnivorous plants, Drosera rotundifolia and D. intermedia, was quantified in situ and without any experimental manipulation using natural abundance stable isotope analysis. • Samples of D. rotundifolia and D. intermedia, insects and noncarnivorous reference plants were collected from three contrasting locations across Britain. The proportion of Drosera nitrogen obtained from insect prey was calculated by a mixing model using δ15 N values from the different plant groups. • The mean proportion of Drosera N derived from prey was 50%. There were significant differences in this proportion between sites, and significant differences within sites. There were significant differences between plant tissues and a significant negative relationship between the proportion of N derived from prey and the C : N ratio of Drosera tissues. • There was little evidence of differences in prey capture/utilisation in response to N availability, possibly due to a limited range in available N between the sites. However, evidence of a positive benefit of prey capture was apparent through the decrease in C : N ratio with increasing prey N concentrations in the plants.