Relative contribution of multiple stressors on copepod density and diversity dynamics in the Belgian part of the North Sea

Contribution of combined stressors on density and gene expression dynamics of the copepod Temora longicornis in the North Sea

The impact of multiple environmental and anthropogenic stressors on the marine environment remains poorly understood. Therefore, we studied the contribution of environmental variables to the densities and gene expression of the dominant zooplankton species in the Belgian part of the North Sea, the calanoid copepod Temora longicornis. We observed a reduced density of copepods, which were also smaller in size, in samples taken from nearshore locations when compared to those obtained from offshore stations. To assess the factors influencing the population dynamics of this species, we applied generalised additive models. These models allowed us to quantify the relative contribution of temperature, nutrient levels, salinity, turbidity, concentrations of photosynthetic pigments, as well as chemical pollutants such as polychlorinated biphenyls and polycyclic aromatic hydrocarbons (PAHs), on copepod density. Temperature and Secchi depth, a proxy for turbidity, were the most important environmental variables predicting the densities of T. longicornis, followed by summed PAH and chlorophyll concentrations. Analysing gene expression in field-collected adults, we observed significant variation in metabolic and stress-response genes. Temperature correlated significantly with genes involved in proteolytic activities, and encoding heat shock proteins. Yet, concentrations of anthropogenic chemicals did not induce significant differences in the gene expression of genes involved in the copepod's fatty acid metabolism or well-known stress-related genes, such as glutathione transferases or cytochrome P450. Our study highlights the potential of gene expression biomonitoring and underscores the significance of a changing environment in future studies.

Environmentally realistic concentrations of chlorinated, brominated, and fluorinated persistent organic pollutants induce the unfolded protein response as a shared stress pathway in the liver of Atlantic cod (Gadus morhua)

In the North Sea and North Atlantic coastal areas, fish experience relatively high background levels of persistent organic pollutants. This study aimed to compare the mode of action of environmentally relevant concentrations of mixtures of halogenated compounds in Atlantic cod. Juvenile male cod with mean weight of 840 g were exposed by gavage to dietary mixtures of chlorinated (PCBs, DDT analogs, chlordane, lindane, and toxaphene), brominated (PBDEs), and fluorinated (PFOS) compounds for 4 weeks. One group received a combined mixture of all three compound groups. The results showed that the accumulated levels of chemicals in cod liver after 4 weeks of exposure reflected concentrations found in wild fish in this region. Pathway analysis revealed that the treatment effects by each of the three groups of chemicals (chlorinated, brominated, and fluorinated) converged on activation of the unfolded protein response (UPR). Upstream regulator analysis predicted that almost all the key transcription factors (XBP1, ERN1, ATF4, EIF2AK3, and NFE2L2) regulating the UPR were significantly activated. No additive effect was observed in cod co-treated with all three compound groups. In conclusion, the genome-wide transcriptomic study suggests that the UPR pathway is a sensitive common target of halogenated organic environmental pollutants in fish.

Decadal decline of dominant copepod species in the North Sea is associated with ocean warming: Importance of marine heatwaves

Time-series are crucial to understand the status of zooplankton communities and to anticipate changes that might affect the entire food web. Long-term time series allow us to understand impacts of multiple environmental and anthropogenic stressors, such as chemical pollution and ocean warming, on the marine ecosystems. Here, a recent time series (2018-2022) of abundance data of four dominant calanoid and one harpacticoid copepod species from the Belgian Part of the North Sea was combined with previously collected (2009-2010, 2015-2016) datasets for the same study area. The time series reveals a significant decrease (up to two orders of magnitude) in calanoid copepod abundance (Temora longicornis, Acartia clausi, Centropages spp., Calanus helgolandicus), while this was not the case for the harpacticoid Euterpina acutifrons. We applied generalized additive models to quantify the relative contribution of temperature, nutrients, salinity, primary production, turbidity and pollution (anthropogenic chemicals, i.e., polychlorinated biphenyls and polycyclic aromatic hydrocarbons) to the population dynamics of these species. Temperature, turbidity and chlorophyll a concentrations were the only variables consistently showing a relative high contribution in all models predicting the abundances of the selected species. The observed heat waves which occurred during the summer periods of the investigated years coincided with population collapses (versus population densities in non-heatwave years) and are considered the most likely cause for the observed copepod abundance decreases. Moreover, the recorded water temperatures during these heatwaves correspond to the physiological thermal limit of some of the studied species. As far as we know, this is the first study to observe ocean warming and marine heat waves having such a dramatic impact (population collapse) on the dominant zooplankton species in shallow coastal areas.

Wild canids and the ecological traps facing the climate change and deforestation in the Amazon Forest

Ecological traps occur when species choose to settle in lower-quality habitats, even if this reduces their survival or productivity. This happens in situations of drastic environmental changes, resulting from anthropogenic pressures. In long term, this could mean the extinction of the species. We investigated the dynamics of occurrence and distribution of three canid species (Atelocynus microtis, Cerdocyon thous, and Spheotos venaticus) considering human threats to their habitats in the Amazon Rainforest. We analyzed the environmental thresholds for the occurrence of these species and related to the future projections of climatic niches for each one. All three species will be negatively affected by climate change in the future, with losses of up to 91% of the suitable area of occurrence in the Brazilian Amazon. A. microtis appear to be more forest-dependent and must rely on the goodwill of decision-makers to be maintained in the future. For C. thous and S. venaticus, climatic variables and those associated with anthropogenic disturbances that modulate their niches today may not act the same way in the future. Even though C. thous is least dependent on the Amazon Forest; this species may be affected in the future due to the ecological traps. S. venaticus, can also undergo the same process, but perhaps more drastically due to the lower ecological plasticity of this species compared to C. thous. Our results suggest that the ecological traps may put these two species at risk in the future. Using the canid species as a model, we had the opportunity to investigate these ecological effects that can affect a large part of the Amazonian fauna in the current scenario. Considering the high degree of environmental degradation and deforestation in the Amazon Rainforest, the theory of ecological traps must be discussed at the same level as the habitat loss, considering the strategies for preserving the Amazon biodiversity.

Copepod Assemblages in A Large Arctic Coastal Area: A Baseline Summer Study

To provide a baseline description of copepod assemblages in the Pechora Sea, an estuarine area with great economical and ecological importance, we conducted a survey during the summer season. A total of 24 copepod taxa were identified in the study, with Acartia longiremis, Calanus finmarchicus, Centropages hamatus, Copepoda nauplii, Eurytemora affinis, Oithona similis, Pseudocalanus spp., and Temora longicornis being the most numerous. The high diversity (Shannon index = 2.51 ± 0.06), density (18,720 ± 3376 individuals m−3) and biomass (89 ± 18 mg dry mass m−3) of copepods were revealed. Populations of common small copepod taxa were dominated by the young stages, indicating spawning, while older copepodites prevailed among medium- and large-sized species, showing that their reproduction occurred before our survey. Cluster analysis indicated three groups of stations that mainly differed in the abundance of particular species. There were clear associations between copepod assemblages and environmental variables. Statistical analyses showed significant correlations between copepod abundance and water temperature or sampling depth, while other factors had a lesser influence. Our results suggest a strong effect of local circulation and currents on the spatial pattern of the copepod assemblages in the study area. This study may be useful for future biomonitoring in the south-eastern Barents Sea.

Neonicotinoid Insecticides from a Marine Perspective: Acute and Chronic Copepod Testing and Derivation of Environmental Quality Standards

Neonicotinoid insecticides have become of global concern for the aquatic environment. Harpacticoid copepods are among the organisms most sensitive to neonicotinoids. We exposed the brackish copepod Nitocra spinipes to 4 neonicotinoid insecticides (clothianidin, imidacloprid, thiacloprid, and thiamethoxam) to investigate acute toxicity on adults (96-h exposure) and effects on larval development (7-d exposure). We used these results in combination with publicly available ecotoxicity data to derive environmental quality standards (EQS). These EQS were ultimately used in a single-substance and mixture risk assessment for the Belgian part of the North Sea. Acute toxicity testing revealed that immobilization is a more sensitive endpoint than mortality, with 96-h median effect concentration (EC50) values of 6.9, 7.2, 25, and 120 µg L^-1 for clothianidin, thiacloprid, imidacloprid, and thiamethoxam, respectively. In addition, the larval development tests resulted in 7-d no-observed-effect concentrations (NOECs) of 2.5, 2.7, 4.2, and >99 µg L^-1 for clothianidin, thiacloprid, imidacloprid, and thiamethoxam, respectively. The derived saltwater annual average (AA-)EQS were 0.05, 0.0048, 0.002, and 0.016 µg L^-1 for clothianidin, thiacloprid, imidacloprid, and thiamethoxam, respectively. Finally, the risk characterization revealed some exceedances of the AA-EQS in Belgian harbors for imidacloprid (number of exceedances, n = 2/4), for thiacloprid (n = 1/4), for thiamethoxam (n = 1/4), and for the mixture of the 4 neonicotinoids (n = 4/4), but not at the open sea. At the open sea site, the toxic unit sums relative to the AA-EQS were 0.72 and 0.22, suggesting no mixture risk, albeit with a relatively small margin of safety. Including short-term EC10 (96-h) values of N. spinipes for the AA-EQS derivation led to a refinement of the AA-EQS for clothianidin and thiamethoxam, suggesting their use for the AA-EQS derivation because one of the overarching goals of the definition of EQS is to protect species at the population level. Environ Toxicol Chem 2021;40:1353-1367. © 2021 SETAC.

Climate change and oil pollution: A dangerous cocktail for tropical zooplankton

Climate change and oil pollution pose a major threat to tropical marine ecosystems and to the coastal communities relying on their resources. The Gulf of Guinea is severely affected by multiple human induced stressors, but the potential impacts of these on marine productivity remain unknown. We investigated the combined effects of heatwaves (climate stressor) and the polycyclic aromatic hydrocarbon pyrene (proxy for oil) on the copepod Centropages velificatus. We quantified survival, reproduction and fecal pellet production of females exposed to concentrations of 0, 10, 100 and 100+ nM (saturated) pyrene under simulated heatwaves of different thermal intensity (+3 °C and +5 °C above control treatment temperature). Thermal stress due to both moderate and intensive heatwaves resulted in reduced survival and egg production. The negative effects of pyrene were only measurable at the high pyrene concentrations. However, thermal stress increased the sensitivity of C. velificatus to pyrene, indicating a synergistic interaction between the two stressors. We document that the interaction of multiple stressors can result in cumulative impacts that are stronger than expected based on single stressor studies. Further research is urgently needed to evaluate the combined impact of climatic and anthropogenic stressors on the productivity of coastal ecosystems, particularly in the tropical areas.

Spatio-temporal patterns in the gene expression of the calanoid copepod Temora longicornis in the Belgian part of the North Sea

Marine zooplankton are increasingly being affected by recent environmental changes, such as climate change, and respond with profound spatial relocations and shifts in phenology and physiology. In order to predict whether populations are able to persist or adapt to such new conditions, it is essential to understand the molecular basis of such adaptations, which ultimately get translated into these physiological responses. To explore variation in population gene expression across time and space, we investigated transcriptome-level profiles of the calanoid copepod Temora longicornis, that were collected at four different locations in the Belgian Part of the North Sea (BPNS) on three different time points (April, June, October) in 2018. RNA-seq analysis of field collected adults identified large seasonal differences in gene expression, mainly between spring-summer and autumn samples. The largest log-fold changes occurred in a set of genes encoding for ribosomal and myosin (heavy chain) transcripts. Enrichment analysis revealed a strong seasonal pattern in vitellogenin, cuticle and glycolytic gene expression as well. We also found a positive correlation between vitellogenin expression and densities of T. longicornis. No clear spatial variation in expression patterns was found in the BPNS. This study underlines the potential of field gene expression studies for biomonitoring purposes and the significance of considering seasonal variation in future studies.

Evaluating the potential of direct RNA nanopore sequencing: Metatranscriptomics highlights possible seasonal differences in a marine pelagic crustacean zooplankton community

The implementation of cost-effective monitoring programs for zooplankton remains challenging due to the requirements of taxonomical expertise and the high costs of sampling and species identification. To reduce costs, molecular methods have been proposed as alternatives to morphology-based monitoring. Metatranscriptomics can contribute to promote both cost-effectiveness and accuracy of biological assessments of aquatic ecosystems. Here, we describe and evaluate the construction of a metatranscriptome dataset from a pelagic crustacean zooplankton community. We sampled zooplankton in one marine station, named LW02, in the North Sea, in both winter and summer, and generated transcripts using Oxford Nanopore Technology (ONT), a third-generation nanopore-based sequencing technology. ONT is, uniquely, capable of sequencing RNA directly, rather than depending on reverse transcription and PCR, and applicable to be used directly in the field. We found that metatranscriptomics is capable of species detection, including screening for the presence of endoparasites, hence competing with morphological identification. Taxonomic analysis based on ribosomal 18S transcripts identified calanoid copepods, particularly Temora longicornis and Acartia clausi, as the most abundant community members. Moreover, up to 40.4% and 50.5% of all sequences could be assigned to predicted genes in the winter and summer sample, respectively. The most abundant mRNA transcripts with known function coded for essential metabolic processes. GO term annotation revealed that genes involved in glycolytic and translation-related processes were most expressed in the community. Although small in scale, our study provides the basis for future efforts to characterize the metatranscriptome of marine zooplankton communities and its application in biomonitoring programs.

Seasonal and spatial fatty acid profiling of the calanoid copepods Temora longicornis and Acartia clausi linked to environmental stressors in the North Sea

The Belgian part of the North Sea (BPNS) is subjected to multiple environmental stressors. The impact of these stressors includes the modulation of fatty acid (FA) composition of the zooplankton. This study recorded temporal and spatial patterns of the FA profiles of two dominant calanoid copepods within the BPNS: Temora longicornis (Müller, 1785) and Acartia clausi (Giesbrecht, 1889). By means of distance-based linear modelling and by applying multi model inference to generalized additive models, environmental stressors were linked to patterns of the FA profiles of these species. The FA profiles of A. clausi and T. longicornis showed distinct intraspecific, spatial and temporal differences within the BPNS. Temperature and algal food quality (marked by the ratio of silicate concentration to dissolved inorganic nitrogen concentration, SiO₄/DIN) were the most important drivers of seasonal fluctuations in the DHA/EPA ratio of both species. DHA/EPA ratio can be used as marker for stress in copepods in the BPNS in order to have a quick indication of food quality changes at the basis of the food web.

The transcriptome of the marine calanoid copepod Temora longicornis under heat stress and recovery

Understanding the impacts of global change in zooplankton communities is crucial, as alterations in the zooplankton communities can affect entire marine ecosystems. Despite the economic and ecological importance of the calanoid copepod Temora longicornis in the Belgian part of the North Sea, molecular data is still very limited for this species. Using HiSeq Illumina sequencing, we sequenced the whole transcriptome of T. longicornis, after being exposed to realistic temperatures of 14 and 17 °C. After both an acute (1 day) and a more sustained (5 days) thermal exposure to 17 °C, we investigated gene expression differences with animals exposed to 14 °C, which may be critical for the thermal acclimation and resilience of this copepod species. We also studied the possibility of a short term stress recovery of a heat shock. A total of 179,569 transcripts were yielded, of which 44,985 putative ORF transcripts were identified. These transcripts were subsequently annotated into roughly 22,000 genes based on known sequences using Gene Ontology (GO) and KEGG databases. Temora only showed a mild response to both the temperature and the duration of the exposure. We found that the expression of 27 transcripts varied significantly with an increase in temperature of 3 °C, of which eight transcripts were differentially expressed after acute exposure only. Gene set enrichment analysis revealed that, overall, T. longicornis was more impacted by a sustained thermal exposure, rather than an immediate (acute) exposure, with two times as many enriched GO terms in the sustained treatment. We also identified several general stress responses independent of exposure time, such as modified protein synthesis, energy mobilisation, cuticle and chaperone proteins. Finally, we highlighted candidate genes of a possible recovery from heat exposure, identifying similar terms as those enriched in the heat treatments, i.e. related to for example energy metabolism, cuticle genes and extracellular matrix. The data presented in this study provides the first transcriptome available for T. longicornis which can be used for future genomic studies.