Spatial and temporal variability in Scotian Shelf zooplankton communities

Zooplankton are indicators of the response of marine ecosystems to environmental variability. The relationships between zooplankton communities on the Scotian Shelf and hydrographic and geographic regions of the Scotian Shelf in the 1990s and 2000s were described using complementary data sets, each resolving different space and time scales. The Atlantic Zone Monitoring Program (AZMP) sampled Scotian Shelf zooplankton from the whole water column twice per year at stations along three cross-shelf transects and semi-monthly at a fixed station on the inshore central shelf, while Continuous Plankton Recorder (CPR) samples were collected from near-surface waters approximately monthly on an along-shelf transect. Variability patterns were compared among these three data sets to identify robust spatial and interannual trends. Stations were clustered based on taxonomic composition, and spatial clusters were compared to hydrographic boundaries and bathymetry to determine whether temporal changes in community composition were driven by changes in water mass distributions on the shelf. This project identifies zooplankton community and abundance shifts that may affect fish recruitment in the northwest Atlantic and contributes to development of ecosystem-based fisheries management on the Scotian Shelf.

Spatio-Temporal Variability of the North Sea Cod Recruitment in Relation to Temperature and Zooplankton

The North Sea cod (Gadus morhua, L.) stock has continuously declined over the past four decades linked with overfishing and climate change. Changes in stock structure due to overfishing have made the stock largely dependent on its recruitment success, which greatly relies on environmental conditions. Here we focus on the spatio-temporal variability of cod recruitment in an effort to detect changes during the critical early life stages. Using International Bottom Trawl Survey (IBTS) data from 1974 to 2011, a major spatio-temporal change in the distribution of cod recruits was identified in the late 1990s, characterized by a pronounced decrease in the central and southeastern North Sea stock. Other minor spatial changes were also recorded in the mid-1980s and early 1990s. We tested whether the observed changes in recruits distribution could be related with direct (i.e. temperature) and/or indirect (i.e. changes in the quantity and quality of zooplankton prey) effects of climate variability. The analyses were based on spatially-resolved time series, i.e. sea surface temperature (SST) from the Hadley Center and zooplankton records from the Continuous Plankton Recorder Survey. We showed that spring SST increase was the main driver for the most recent decrease in cod recruitment. The late 1990s were also characterized by relatively low total zooplankton biomass, particularly of energy-rich zooplankton such as the copepod Calanus finmarchicus, which have further contributed to the decline of North Sea cod recruitment. Long-term spatially-resolved observations were used to produce regional distribution models that could further be used to predict the abundance of North Sea cod recruits based on temperature and zooplankton food availability.

Temporal change in UK marine communities: trends or regime shifts?

A regime shift is a large, sudden, and long-lasting change in the dynamics of an ecosystem, affecting multiple trophic levels. There are a growing number of papers that report regime shifts in marine ecosystems. However, the evidence for regime shifts is equivocal, because the methods used to detect them are not yet well developed. We have collated over 300 biological time series from seven marine regions around the UK, covering the ecosystem from phytoplankton to marine mammals. Each time series consists of annual measures of abundance for a single group of organisms over several decades. We summarised the data for each region using the first principal component, weighting either each time series or each biological component (e.g. plankton, fish, benthos) equally. We then searched for regime shifts using Rodionov’s regime shift detection (RSD) method, which found regime shifts in the first principal component for all seven marine regions. However, there are consistent temporal trends in the data for six of the seven regions. Such trends violate the assumptions of RSD. Thus, the regime shifts detected by RSD in six of the seven regions are likely to be artefacts caused by temporal trends. We are therefore developing more appropriate time series models for both single populations and whole communities that will explicitly model temporal trends and should increase our ability to detect true regime shift events.

The emergence of automated high-frequency flow cytometry: revealing temporal and spatial phytoplankton variability

Phytoplankton observation is the product of a number of trade-offs related to sampling processes, required level of diversity and size spectrum analysis capabilities of the techniques involved. Instruments combining the morphological and high-frequency analysis for phytoplankton cells are now available. This paper presents an application of the automated high-resolution flow cytometer Cytosub as a tool for analysing phytoplanktonic cells in their natural environment. High resolution data from a temporal study in the Bay of Marseille (analysis every 30 min over 1 month) and a spatial study in the Southern Indian Ocean (analysis every 5 min at 10 knots over 5 days) are presented to illustrate the capabilities and limitations of the instrument. Automated high-frequency flow cytometry revealed the spatial and temporal variability of phytoplankton in the size range 1−∼50 μm that could not be resolved otherwise. Due to some limitations (instrumental memory, volume analysed per sample), recorded counts could be statistically too low. By combining high-frequency consecutive samples, it is possible to decrease the counting error, following Poisson’s law, and to retain the main features of phytoplankton variability. With this technique, the analysis of phytoplankton variability combines adequate sampling frequency and effective monitoring of community changes.

Temporal and spatial variation in the Nazaré Canyon (Western Iberian margin): Inter-annual and canyon heterogeneity effects on meiofauna biomass and diversity

The Nazaré Canyon on the Portuguese Margin (NE Atlantic) was sampled during spring-summer for three consecutive years (2005–2007), permitting the first inter-annual study of the meiofaunal communities at the Iberian Margin at two abyssal depths (~3500 m and ~4400 m). Using new and already published data, the meiofauna standing stocks (abundance and biomass) and nematode structural and functional diversity were investigated in relation to the sediment biogeochemistry (e.g. organic carbon, nitrogen, chlorophyll a, phaeopigments) and grain size. A conspicuous increase in sand content from 2005 to 2006 and decrease of phytodetritus at both sites, suggested the occurrence of one or more physical disturbance events. Nematode standing stocks and trophic diversity decreased after these events, seemingly followed by a recovery/recolonisation period in 2007, which was strongly correlated with an increase in the quantity and bioavailability of phytodetrital organic matter supplied. Changes in meiofauna assemblages, however, also differed between stations, likely because of the contrasting hydrodynamic and food supply conditions. Higher meiofauna and nematode abundances, biomass and trophic complexity were found at the shallowest canyon station, where the quantity, quality and bioavailability of food material were higher than at the deeper site. The present results suggest that even though inter-annual variations in the sedimentary environment can regulate the meiofauna in the abyssal Nazaré Canyon, heterogeneity between sampling locations in the canyon were more pronounced.

Impacts of bioturbation on temporal variation in bacterial and archaeal nitrogen-cycling gene abundance in coastal sediments

In marine environments, macrofauna living in or on the sediment surface may alter the structure, diversity and function of benthic microbial communities. In particular, microbial nitrogen (N)-cycling processes may be enhanced by the activity of large bioturbating organisms. Here, we study the effect of the burrowing mud shrimp Upogebia deltaura upon temporal variation in the abundance of genes representing key N-cycling functional guilds. The abundance of bacterial genes representing different N-cycling guilds displayed different temporal patterns in burrow sediments in comparison with surface sediments, suggesting that the burrow provides a unique environment where bacterial gene abundances are influenced directly by macrofaunal activity. In contrast, the abundances of archaeal ammonia oxidizers varied temporally but were not affected by bioturbation, indicating differential responses between bacterial and archaeal ammonia oxidizers to environmental physicochemical controls. This study highlights the importance of bioturbation as a control over the temporal variation in nitrogen-cycling microbial community dynamics within coastal sediments.

On the temporal consistency of chlorophyll products derived from three ocean-colour sensors

Satellite ocean-colour sensors have life spans lasting typically five-to-ten years. Detection of long-term trends in chlorophyll-a concentration (Chl-a) using satellite ocean colour thus requires the combination of different ocean-colour missions with sufficient overlap to allow for cross-calibration. A further requirement is that the different sensors perform at a sufficient standard to capture seasonal and inter-annual fluctuations in ocean colour. For over eight years, the SeaWiFS, MODIS-Aqua and MERIS ocean-colour sensors operated in parallel. In this paper, we evaluate the temporal consistency in the monthly Chl-a time-series and in monthly inter-annual variations in Chl-a among these three sensors over the 2002–2010 time period. By subsampling the monthly Chl-a data from the three sensors consistently, we found that the Chl-a time-series and Chl-a anomalies among sensors were significantly correlated for >90% of the global ocean. These correlations were also relatively insensitive to the choice of three Chl-a algorithms and two atmospheric-correction algorithms. Furthermore, on the subsampled time-series, correlations between Chl-a and time, and correlations between Chl-a and physical variables (sea-surface temperature and sea-surface height) were not significantly different for >92% of the global ocean. The correlations in Chl-a and physical variables observed for all three sensors also reflect previous theories on coupling between physical processes and phytoplankton biomass. The results support the combining of Chl-a data from SeaWiFS, MODIS-Aqua and MERIS sensors, for use in long-term Chl-a trend analysis, and highlight the importance of accounting for differences in spatial sampling among sensors when combining ocean-colour observations.

How many Coccolithovirus genotypes does it take to terminate an Emiliania huxleyi bloom?

Giant viruses are known to be significant mortality agents of phytoplankton, often being implicated in the terminations of large Emiliania huxleyi blooms. We have previously shown the high temporal variability of E. huxleyi-infecting coccolithoviruses (EhVs) within a Norwegian fjord mesocosm. In the current study we investigated EhV dynamics within a naturally-occurring E. huxleyi bloom in the Western English Channel. Using denaturing gradient gel electrophoresis and marker gene sequencing, we uncovered a spatially highly dynamic Coccolithovirus population that was associated with a genetically stable E. huxleyi population as revealed by the major capsid protein gene (mcp) and coccolith morphology motif (CMM), respectively. Coccolithoviruses within the bloom were found to be variable with depth and unique virus populations were detected at different stations sampled indicating a complex network of EhV-host infections. This ultimately will have significant implications to the internal structure and longevity of ecologically important E. huxleyi blooms.

Spatio-temporal variability in ammonia oxidation and ammonia oxidising bacteria and archaea in coastal sediments of the Western English Channel

The abundance of ammonia-oxidising bacterial (AOB) and ammonia-oxidising archaeal (AOA) (amoA) genes and ammonia oxidation rates were compared bimonthly from July 2008 to May 2011 in 4 contrasting coastal sediments in the western English Channel. Despite a higher abundance of AOA amoA genes within all sediments and at all time-points, rates of ammonia oxidation correlated with AOB and not AOA amoA gene abundance. Sediment type was a major factor in determining both AOB amoA gene abundance and AOB community structure, possibly due to deeper oxygen penetration into the sandier sediments, increasing the area available for ammonia oxidation. Decreases in AOB amoA gene abundance were evident during summer and autumn, with maximum abundance and ammonia oxidation rates occurring in winter and early spring. PCR-DGGE of AOB amoA genes indicated that no seasonal changes to community composition occurred; however, a gradual movement in community composition occurred at 3 of the sites studied. The lack of correlation between AOA amoA gene abundance and ammonium oxidation rates, or any other environmental variable measured, may be related to the higher spatial variation amongst measurements, obscuring temporal trends, or the bimonthly sampling, which may have been too infrequent to capture temporal variability in the deposition of fresh organic matter. Alternatively, AOA may respond to changing substrate concentrations by an increase or decrease in transcript rather than gene abundance.

Temporal changes in total and size-fractioned chlorophyll-a in surface waters of three provinces in the Atlantic Ocean (September to November) between 2003 and 2010

Phytoplankton total chlorophyll concentration (TCHLa) and phytoplankton size structure are two important ecological indicators in biological oceanography. Using high performance liquid chromatography (HPLC) pigment data, collected from surface waters along the Atlantic Meridional Transect (AMT), we examine temporal changes in TCHLa and phytoplankton size class (PSC: micro-, nano- and pico-phytoplankton) between 2003 and 2010 (September to November cruises only), in three ecological provinces of the Atlantic Ocean. The HPLC data indicate no significant change in TCHLa in northern and equatorial provinces, and an increase in the southern province. These changes were not significantly different to changes in TCHLa derived using satellite ocean-colour data over the same study period. Despite no change in AMT TCHLa in northern and equatorial provinces, significant differences in PSC were observed, related to changes in key diagnostic pigments (fucoxanthin, peridinin, 19′-hexanoyloxyfucoxanthin and zeaxanthin), with an increase in small cells (nano- and pico-phytoplankton) and a decrease in larger cells (micro-phytoplankton). When fitting a three-component model of phytoplankton size structure — designed to quantify the relationship between PSC and TCHLa to each AMT cruise, model parameters varied over the study period. Changes in the relationship between PSC and TCHLa have wide implications in ecology and marine biogeochemistry, and provide key information for the development and use of empirical ocean-colour algorithms. Results illustrate the importance of maintaining a time-series of in-situ observations in remote regions of the ocean, such as that acquired in the AMT programme.

Temporal variability in total, micro- and nano-phytoplankton primary production at a coastal site in the Western English Channel

Primary productivity and subsequent carbon cycling in the coastal zone have a significant impact on the global carbon budget. It is currently unclear how anthropogenic activity could alter these budgets but long term coastal time series of hydrological, biogeochemical and biological measurements represent a key means to better understand past drivers, and hence to predicting future seasonal and inter-annual variability in carbon fixation in coastal ecosystems. An 8-year time series of primary production from 2003 to 2010, estimated using a recently developed absorption-based algorithm, was used to determine the nature and extent of change in primary production at a coastal station (L4) in the Western English Channel (WEC). Analysis of the seasonal and inter-annual variability in production demonstrated that on average, nano- and pico-phytoplankton account for 48% of the total carbon fixation and micro-phytoplankton for 52%. A recent decline in the primary production of nano- and pico-phytoplankton from 2005 to 2010 was observed, corresponding with a decrease in winter nutrient concentrations and a decrease in the biomass of Phaeocystis sp. Micro-phytoplankton primary production (PPM) remained relatively constant over the time series and was enhanced in summer during periods of high precipitation. Increases in sea surface temperature, and decreases in wind speeds and salinity were associated with later spring maxima in PPM. Together these trends indicate that predicted increases in temperature and decrease in wind speeds in future would drive later spring production whilst predicted increases in precipitation would also continue these blooms throughout the summer at this site.

Next Generation Sequencing Reveals the Hidden Diversity of Zooplankton Assemblages

Zooplankton play an important role in our oceans, in biogeochemical cycling and providing a food source for commercially important fish larvae. However, difficulties in correctly identifying zooplankton hinder our understanding of their roles in marine ecosystem functioning, and can prevent detection of long term changes in their community structure. The advent of massively parallel next generation sequencing technology allows DNA sequence data to be recovered directly from whole community samples. Here we assess the ability of such sequencing to quantify richness and diversity of a mixed zooplankton assemblage from a productive time series site in the Western English Channel. Methodology/Principle Findings Plankton net hauls (200 µm) were taken at the Western Channel Observatory station L4 in September 2010 and January 2011. These samples were analysed by microscopy and metagenetic analysis of the 18S nuclear small subunit ribosomal RNA gene using the 454 pyrosequencing platform. Following quality control a total of 419,041 sequences were obtained for all samples. The sequences clustered into 205 operational taxonomic units using a 97% similarity cut-off. Allocation of taxonomy by comparison with the National Centre for Biotechnology Information database identified 135 OTUs to species level, 11 to genus level and 1 to order, <2.5% of sequences were classified as unknowns. By comparison a skilled microscopic analyst was able to routinely enumerate only 58 taxonomic groups. Conclusions Metagenetics reveals a previously hidden taxonomic richness, especially for Copepoda and hard-to-identify meroplankton such as Bivalvia, Gastropoda and Polychaeta. It also reveals rare species and parasites. We conclude that Next Generation Sequencing of 18S amplicons is a powerful tool for elucidating the true diversity and species richness of zooplankton communities. While this approach allows for broad diversity assessments of plankton it may become increasingly attractive in future if sequence reference libraries of accurately identified individuals are better populated.

High-quality RNA extraction from copepods for Next Generation Sequencing: A comparative study

Despite the ecological importance of copepods, few Next Generation Sequencing studies (NGS) have been performed on small crustaceans, and a standard method for RNA extraction is lacking. In this study, we compared three commonly-used methods: TRIzol®, Aurum Total RNA Mini Kit and Qiagen RNeasy Micro Kit, in combination with preservation reagents TRIzol® or RNAlater®, to obtain high-quality and quantity of RNA from copepods for NGS. Total RNA was extracted from the copepods Calanus helgolandicus, Centropages typicus and Temora stylifera and its quantity and quality were evaluated using NanoDrop, agarose gel electrophoresis and Agilent Bioanalyzer. Our results demonstrate that preservation of copepods in RNAlater® and extraction with Qiagen RNeasy Micro Kit were the optimal isolation method for high-quality and quantity of RNA for NGS studies of C. helgolandicus. Intriguingly, C. helgolandicus 28S rRNA is formed by two subunits that separate after heat-denaturation and migrate along with 18S rRNA. This unique property of protostome RNA has never been reported in copepods. Overall, our comparative study on RNA extraction protocols will help increase gene expression studies on copepods using high-throughput applications, such as RNA-Seq and microarrays.