Transcriptional profiling of reproductive development, lipid storage and molting throughout the last juvenile stage of the marine copepod Calanus finmarchicus

Revealing the profound influence of diapause on gene expression: Insights from the annual transcriptome of the copepod Calanus finmarchicus

Annual rhythms are observed in living organisms with numerous ecological implications. In the zooplanktonic copepod Calanus finmarchicus, such rhythms are crucial regarding its phenology, body lipid accumulation, and global carbon storage. Climate change drives annual biological rhythms out of phase with the prevailing environmental conditions with yet unknown but potentially catastrophic consequences. However, the molecular dynamics underlying phenology are still poorly described. In a rhythmic analysis of C. finmarchicus annual gene expression, results reveal that more than 90% of the transcriptome shows significant annual rhythms, with abrupt and dramatic upheaval between the active and diapause life cycle states. This work explores the implication of the circadian clock in the annual timing, which may control epigenetic mechanisms to profoundly modulate gene expression in response to calendar time. Results also suggest an increased light sensitivity during diapause that would ensure the photoperiodic entrainment of the endogenous annual clock.

Effects of petrogenic pollutants on North Atlantic and Arctic Calanus copepods: From molecular mechanisms to population impacts

Oil and gas industries in the Northern Atlantic Ocean have gradually moved closer to the Arctic areas, a process expected to be further facilitated by sea ice withdrawal caused by global warming. Copepods of the genus Calanus hold a key position in these cold-water food webs, providing an important energetic link between primary production and higher trophic levels. Due to their ecological importance, there is a concern about how accidental oil spills and produced water discharges may impact cold-water copepods. In this review, we summarize the current knowledge of the toxicity of petroleum on North Atlantic and Arctic Calanus copepods. We also review how recent development of high-quality transcriptomes from RNA-sequencing of copepods have identified genes regulating key biological processes, like molting, diapause and reproduction in Calanus copepods, to suggest linkages between exposure, molecular mechanisms and effects on higher levels of biological organization. We found that the available ecotoxicity threshold data for these copepods provide valuable information about their sensitivity to acute petrogenic exposures; however, there is still insufficient knowledge regarding underlying mechanisms of toxicity and the potential for long-term implications of relevance for copepod ecology and phenology. Copepod transcriptomics has expanded our understanding of how key biological processes are regulated in cold-water copepods. These advances can improve our understanding of how pollutants affect biological processes, and thus provide the basis for new knowledge frameworks spanning the effect continuum from molecular initiating events to adverse effects of regulatory relevance. Such efforts, guided by concepts such as adverse outcome pathways (AOPs), enable standardized and transparent characterization and evaluation of knowledge and identifies research gaps and priorities. This review suggests enhancing mechanistic understanding of exposure-effect relationships to better understand and link biomarker responses to adverse effects to improve risk assessments assessing ecological effects of pollutant mixtures, like crude oil, in Arctic areas.

Unmasking microsatellite deceptiveness and debunking hybridization with SNPs in four marine copepod species of Calanus

Interspecific hybridization events are on the rise in natural systems due to climate change disrupting species barriers. Across taxa, microsatellites have long been the molecular markers of choice to identify admixed individuals. However, with the advent of high-throughput sequencing easing the generation of genome-wide datasets, incorrect reports of hybridization resulting from microsatellite technical artefacts have been uncovered in a growing number of taxa. In the marine zooplankton genus Calanus (Copepoda), whose species are used as climate change indicators, microsatellite markers have suggested hybridization between C. finmarchicus and C. glacialis, while other nuclear markers (InDels) never detected any admixed individuals, leaving the scientific community divided. Here, for the first time, we investigated the potential for hybridization among C. finmarchicus, C. glacialis, C. helgolandicus and C. hyperboreus using two large and independent SNP datasets. These were derived firstly from a protocol of target-capture applied to 179 individuals collected from 17 sites across the North Atlantic and Arctic Oceans, including sympatric areas, and second from published RNA sequences. All SNP-based analyses were congruent in showing that Calanus species are distinct and do not appear to hybridize. We then thoroughly re-assessed the microsatellites showing hybrids, with the support of published transcriptomes, and identified technical issues plaguing eight out of 10 microsatellites, including size homoplasy, paralogy, potential for null alleles and even two primer pairs targeting the same locus. Our study illustrates how deceptive microsatellites can be when applied to the investigation of hybridization.

Response of moulting genes and gut microbiome to nano-plastics and copper in juvenile horseshoe crab Tachypleus tridentatus

Nanoplastics (NPs) and heavy metals are typical marine pollutants, affecting the gut microbiota composition and molting rate of marine organisms. Currently, there is a lack of research on the toxicological effects of combined exposure to horseshoe crabs. In this study, we investigated the effects of NPs and copper on the expression of molt-related genes and gut microbiome in juvenile tri-spine horseshoe crabs Tachypleus tridentatus by exposing them to NPs (100 nm, 104 particles L-1) and/or Cu2+ (10 μgL-1) in seawater for 21 days. Compared with the control group, the relative mRNA expression of ecdysone receptor (EcR), retinoid x receptor (RXR), calmodulin-A-like isoform X1 (CaM X1), and heat shock 70 kDa protein (Hsp70) were significantly increased under the combined stress of NPs and Cu2+. There were no significant differences in the diversity and abundance indices of the gut microbial population of horseshoe crabs between the NPs and/or Cu2+ groups and the control group. According to linear discriminant analysis, Oleobacillus was the most abundant microorganism in the NPs and Cu2+ stress groups. These results indicate that exposure to either NPs stress alone or combined NPs and Cu2+ stress can promote the expression levels of juvenile molting genes. NPs exposure has a greater impact on the gut microbial community structure of juvenile horseshoe crabs compared to Cu2+ exposure. This study is helpful for predicting the growth and development of horseshoe crabs under complex environmental pollution.

Transcriptome responses in copepods Calanus finmarchicus, Calanus glacialis and Calanus hyperboreus exposed to phenanthrene and benzo[a]pyrene

Arctic and sub-arctic pelagic organisms can be exposed to effluents and spills from offshore petroleum-related activities and thus it is important to understand how they respond to crude oil related contaminants such as polycyclic aromatic hydrocarbons (PAHs). The copepod species Calanus finmarchicus, Calanus glacialis and Calanus hyperboreus represent key links in the arctic marine food web. We performed a transcriptome analysis of the three species exposed to phenanthrene (Phe) and benzo[a]pyrene (BaP) representing low and high molecular weight PAHs, respectively. Differential expression of several genes involved in many cellular pathways was observed after 72 h exposure to Phe (0.1 μM) and BaP (0.1 μM). In C. finmarchicus and C. glacialis, the exposure resulted in up-regulation of genes encoding enzymes in xenobiotic biotransformation, particularly the phase II cytosolic sulfonation system that include 3'-phosphoadenosine 5'-phosphosulfate synthase (PAPSS) and sulfotransferases (SULTs). The sulfonation pathway genes were more strongly induced by BaP than Phe in C. finmarchicus and C. glacialis but were not affected in C. hyperboreus. However, a larger number of genes and pathways were modulated in C. hyperboreus by the PAHs including genes encoding xenobiotic biotransformation and lipid metabolism enzymes, suggesting stronger responses in this species. The results suggest that the cytosolic sulfonation is a major phase II conjugation pathway for PAHs in C. finmarchicus and C. glacialis. Some of the biotransformation systems affected are known to be involved in metabolism of endogenous compounds such as ecdysteroids, which may suggest potential interference with physiological and developmental processes of the copepod species.

Annual transcriptome of a key zooplankton species, the copepod Calanus finmarchicus

The copepod Calanus finmarchicus (Crustacea, Copepoda) is a key zooplanktonic species with a crucial position in the North Atlantic food web and significant contributor to ocean carbon flux. Like many other high latitude animals, it has evolved a programmed arrested development called diapause to cope with long periods of limited food supply, while growth and reproduction are timed to take advantage of seasonal peaks in primary production. However, anthropogenic warming is inducing changes in the expected timing of phytoplankton blooms, suggesting phenological mismatches with negative consequences for the N. Atlantic ecosystem. While diapause mechanisms are mainly studied in terrestrial arthropods, specifically on laboratory model species, such as the fruit fly Drosophila, the molecular investigations of annual rhythms in wild marine species remain fragmentary. Here we performed a rigorous year‐long monthly sampling campaign of C. finmarchicus in a Scottish Loch (UK; 56.45°N, 5.18°W) to generate an annual transcriptome. The mRNA of 36 samples (monthly triplicate of 25 individuals) have been deeply sequenced with an average depth of 137 ± 4 million reads (mean ± SE) per sample, aligned to the reference transcriptome, and filtered. We detail the quality assessment of the datasets and provide a high‐quality resource for the investigation of wild annual transcriptomic rhythms (35,357 components) in a key diapausing zooplanktonic species.

Insights into the species evolution of Calanus copepods in the northern seas revealed by de novo transcriptome sequencing

Copepods of the zooplankton genus Calanus play a key role in marine ecosystems in the northern seas. Although being among the most studied organisms on Earth, due to their ecological importance, genomic resources for Calanus spp. remain scarce, mostly due to their large genome size (from 6 to 12 Gbps). As an alternative to whole-genome sequencing in Calanus spp., we sequenced and de novo assembled transcriptomes of five Calanus species: Calanus glacialis, C. hyperboreus, C. marshallae, C. pacificus, and C. helgolandicus. Functional assignment of protein families based on clusters of orthologous genes (COG) and gene ontology (GO) annotations showed analogous patterns of protein functions across species. Phylogenetic analyses using maximum likelihood (ML) of 191 protein-coding genes mined from RNA-seq data fully resolved evolutionary relationships among seven Calanus species investigated (five species sequenced for this study and two species with published datasets), with gene and site concordance factors showing that 109 out of 191 protein-coding genes support a separation between three groups: the C. finmarchicus group (including C. finmarchicus, C. glacialis, and C. marshallae), the C. helgolandicus group (including C. helgolandicus, C. sinicus, and C. pacificus) and the monophyletic C. hyperboreus group. The tree topology obtained in ML analyses was similar to a previously proposed phylogeny based on morphological criteria and cleared certain ambiguities from past studies on evolutionary relationships among Calanus species.

A novel approach to co-expression network analysis identifies modules and genes relevant for moulting and development in the Atlantic salmon louse (Lepeophtheirus salmonis)

Background

The salmon louse (Lepeophtheirus salmonis) is an obligate ectoparasitic copepod living on Atlantic salmon and other salmonids in the marine environment. Salmon lice cause a number of environmental problems and lead to large economical losses in aquaculture every year. In order to develop novel parasite control strategies, a better understanding of the mechanisms of moulting and development of the salmon louse at the transcriptional level is required.

Methods

Three weighted gene co-expression networks were constructed based on the pairwise correlations of salmon louse gene expression profiles at different life stages. Network-based approaches and gene annotation information were applied to identify genes that might be important for the moulting and development of the salmon louse. RNA interference was performed for validation. Regulatory impact factors were calculated for all the transcription factor genes by examining the changes in co-expression patterns between transcription factor genes and deferentially expressed genes in middle stages and moulting stages.

Results

Eight gene modules were predicted as important, and 10 genes from six of the eight modules have been found to show observable phenotypes in RNA interference experiments. We knocked down five hub genes from three modules and observed phenotypic consequences in all experiments. In the infection trial, no copepodids with a RAB1A-like gene knocked down were found on fish, while control samples developed to chalimus-1 larvae. Also, a FOXO-like transcription factor obtained highest scores in the regulatory impact factor calculation.

Conclusions

We propose a gene co-expression network-based approach to identify genes playing an important role in the moulting and development of salmon louse. The RNA interference experiments confirm the effectiveness of our approach and demonstrated the indispensable role of a RAB1A-like gene in the development of the salmon louse. We propose that our approach could be generalized to identify important genes associated with a phenotype of interest in other organisms.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-08054-7).

Effects of marine mine tailing exposure on the development, growth, and lipid accumulation in Calanus finmarchicus

Marine tailing disposal (MTD) is sometimes practiced as an alternative to traditional mine tailing deposition on land. Environmental challenges connected to MTD include spreading of fine particulate matter in the water column and the potential release of metals and processing chemicals. This study investigated if tailing exposure affects the marine copepod Calanus finmarchicus, and whether effects are related to exposure to mineral particles or the presence of metals and/or processing chemicals in the tailings. We investigated the impacts of three different tailing compositions: calcium carbonate particles with and without processing chemicals and fine-grained tailings from a copper ore. Early life stages of C. finmarchicus were exposed over several developmental stages to low and high suspension concentrations for 15 days, and their development, oxygen consumption and biometry determined. The data was fitted in a dynamic energy budget (DEB) model to determine mechanisms underlying responses and to understand the primary modes of action related to mine tailing exposure. Results show that copepods exposed to tailings generally exhibited slower growth and accumulated less lipids. The presence of metals and processing chemicals did not influence these responses, suggesting that uptake of mineral particles was responsible for the observed effects. This was further supported by the applied DEB model, confirming that ingestion of tailing particles while feeding can result in less energy being available for growth and development.

Diapause vs. reproductive programs: transcriptional phenotypes in a keystone copepod

Many arthropods undergo a seasonal dormancy termed "diapause" to optimize timing of reproduction in highly seasonal environments. In the North Atlantic, the copepod Calanus finmarchicus completes one to three generations annually with some individuals maturing into adults, while others interrupt their development to enter diapause. It is unknown which, why and when individuals enter the diapause program. Transcriptomic data from copepods on known programs were analyzed using dimensionality reduction of gene expression and functional analyses to identify program-specific genes and biological processes. These analyses elucidated physiological differences and established protocols that distinguish between programs. Differences in gene expression were associated with maturation of individuals on the reproductive program, while those on the diapause program showed little change over time. Only two of six filters effectively separated copepods by developmental program. The first one included all genes annotated to RNA metabolism and this was confirmed using differential gene expression analysis. The second filter identified 54 differentially expressed genes that were consistently up-regulated in individuals on the diapause program in comparison with those on the reproductive program. Annotated to oogenesis, RNA metabolism and fatty acid biosynthesis, these genes are both indicators for diapause preparation and good candidates for functional studies.

Comparison between cultured and wild Pacific white shrimp (Penaeus vannamei) vitellogenesis: next-generation sequencing and relative expression of genes directly and indirectly related to reproduction

Shrimp fisheries are among the most important fisheries worldwide, and shrimp culture has increased considerably in recent years. Most current studies on reproduction-related genes have been conducted on cultured shrimp. However, gene expression is intimately linked to physiological and environmental conditions, and therefore an organism’s growth environment has a great influence on reproduction. Thus, gene expression profiling, should be applied in fisheries studies. Here, we identified the expression patterns of 76 reproduction-related genes in P. vannamei via the analysis of pooled transcriptomes from a time-series experiment encompassing a full circadian cycle. The expression patterns of genes associated both directly (Vtg, ODP, and ProR) and indirectly (FAMet, CruA1, and CruC1) with reproduction were evaluated, as these genes could be used as molecular markers of previtellogenic and vitellogenic maturation stages. The evaluated genes were prominently upregulated during vitellogenic stages, with specific expression patterns depending on the organism’s environment, diet, and season. Vtg, ProR, ODP, and FaMet could serve as molecular markers for both wild and cultured organisms.

Transcriptomic analysis elucidates the molecular processes associated with hydrogen peroxide-induced diapause termination in Artemia-encysted embryos

Treatment with hydrogen peroxide (H₂O₂) raises the hatching rate through the development and diapause termination of Artemia cysts. To comprehend the upstream genetic regulation of diapause termination activated by exterior H₂O₂ elements, an Illumina RNA-seq analysis was performed to recognize and assess comparative transcript amounts to explore the genetic regulation of H₂O₂ in starting the diapause termination of cysts in Artemia salina. We examined three groupings treated with no H₂O₂ (control), 180 μM H₂O₂ (low) and 1800 μM H₂O₂ (high). The results showed a total of 114,057 unigenes were identified, 41.22% of which were functionally annotated in at least one particular database. When compared to control group, 34 and 98 differentially expressed genes (DEGs) were upregulated in 180 μM and 1800 μM H₂O₂ treatments, respectively. On the other hand, 162 and 30 DEGs were downregulated in the 180 μM and 1800 μM H₂O₂ treatments, respectively. Cluster analysis of DEGs demonstrated significant patterns among these types of 3 groups. GO and KEGG enrichment analysis showed the DEGs involved in the regulation of blood coagulation (GO: 0030193; GO: 0050818), regulation of wound healing (GO:0061041), regulation of hemostasis (GO: 1900046), antigen processing and presentation (KO04612), the Hippo signaling pathway (KO04391), as well as the MAPK signaling pathway (KO04010). This research helped to define the diapause-related transcriptomes of Artemia cysts using RNA-seq technology, which might fill up a gap in the prevailing body of knowledge.

De novo transcriptome assembly of the Southern Ocean copepod Rhincalanus gigas sheds light on developmental changes in gene expression

Copepods are small crustaceans that dominate most zooplankton communities in terms of both abundance and biomass. In the polar oceans, a subset of large lipid-storing copepods occupy central positions in the food web because of their important role in linking phytoplankton and microzooplankton with higher trophic levels. In this paper, we generated a high-quality de novo transcriptome for Rhincalanus gigas, the largest-and among the most abundant-of the Southern Ocean copepods. We then conducted transcriptional profiling to characterize the developmental transition between late-stage juveniles and adult females. We found that juvenile R. gigas substantially upregulate lipid synthesis and glycolysis pathways relative to females, as part of a developmental gene expression program that also implicates processes such as muscle growth, chitin formation, and ion transport. This study provides the first transcriptional profile of a developmental transition within Rhincalanus gigas or any endemic Southern Ocean copepod, thereby extending our understanding of copepod molecular physiology.

Lipid metabolism in Calanus finmarchicus is sensitive to variations in predation risk and food availability

Late developmental stages of the marine copepods in the genus Calanus can spend extended periods in a dormant stage (diapause) that is preceded by the accumulation of large lipid stores. We assessed how lipid metabolism during development from the C4 stage to adult is altered in response to predation risk and varying food availability, to ultimately understand more of the metabolic processes during development in Calanus copepods. We used RNA sequencing to assess if perceived predation risk in combination with varied food availability affects expression of genes associated with lipid metabolism and diapause preparation in C. finmarchicus. The lipid metabolism response to predation risk differed depending on food availability, time and life stage. Predation risk caused upregulation of lipid catabolism with high food, and downregulation with low food. Under low food conditions, predation risk disrupted lipid accumulation. The copepods showed no clear signs of diapause preparation, supporting earlier observations of the importance of multiple environmental cues in inducing diapause in C. finmarchicus. This study demonstrates that lipid metabolism is a sensitive endpoint for the interacting environmental effects of predation pressure and food availability. As diapause may be controlled by lipid accumulation, our findings may contribute towards understanding processes that can ultimately influence diapause timing.

Insights into the genetic basis of predator-induced response in Daphnia galeata

Phenotypic plastic responses allow organisms to rapidly adjust when facing environmental challenges-these responses comprise morphological, behavioral but also life-history changes. Alteration of life-history traits when exposed to predation risk have been reported often in the ecological and genomic model organism Daphnia. However, the molecular basis of this response is not well understood, especially in the context of fish predation. Here, we characterized the transcriptional profiles of two Daphnia galeata clonal lines with opposed life histories when exposed to fish kairomones. First, we conducted a differential gene expression, identifying a total of 125 candidate transcripts involved in the predator-induced response, uncovering substantial intraspecific variation. Second, we applied a gene coexpression network analysis to find clusters of tightly linked transcripts revealing the functional relations of transcripts underlying the predator-induced response. Our results showed that transcripts involved in remodeling of the cuticle, growth, and digestion correlated with the response to environmental change in D. galeata. Furthermore, we used an orthology-based approach to gain functional information for transcripts lacking gene ontology (GO) information, as well as insights into the evolutionary conservation of transcripts. We could show that our candidate transcripts have orthologs in other Daphnia species but almost none in other arthropods. The unique combination of methods allowed us to identify candidate transcripts, their putative functions, and evolutionary history associated with predator-induced responses in Daphnia. Our study opens up to the question as to whether the same molecular signature is associated with fish kairomones-mediated life-history changes in other Daphnia species.

RNA-Seq and differential gene expression analysis in Temora stylifera copepod females with contrasting non-feeding nauplii survival rates: an environmental transcriptomics study

Background

Copepods are fundamental components of pelagic food webs, but reports on how molecular responses link to reproductive success in natural populations are still scarce. We present a de novo transcriptome assembly and differential expression (DE) analysis in Temora stylifera females collected in the Gulf of Naples, Mediterranean Sea, where this copepod dominates the zooplankton community. High-Throughput RNA-Sequencing and DE analysis were performed from adult females collected on consecutive weeks (May 23rd and 30th 2017), because opposite naupliar survival rates were observed. We aimed at detecting key genes that may have influenced copepod reproductive potential in natural populations and whose expression was potentially affected by phytoplankton-derived oxylipins, lipoxygenase-derived products strongly impacting copepod naupliar survival.

Results

On the two sampling dates, temperature, salinity, pH and oxygen remained stable, while variations in phytoplankton cell concentration, oxylipin concentration and oxylipin-per-diatom-cell production were observed. T. stylifera naupliar survival was 25% on May 23rd and 93% on May 30th. De novo assembly generated 268,665 transcripts (isoforms) and 120,749 unique ‘Trinity predicted genes’ (unigenes), of which 50% were functionally annotated. Out of the 331 transcript isoforms differentially expressed between the two sampling dates, 119 sequences were functionally annotated (58 up- and 61 down-regulated). Among predicted genes (unigenes), 144 sequences were differentially expressed and 31 (6 up-regulated and 25 down-regulated) were functionally annotated. Most of the significantly down-regulated unigenes and isoforms were A5 Putative Odorant Binding Protein (Obp). Other differentially expressed sequences (isoforms and unigenes) related to developmental metabolic processes, protein ubiquitination, response to stress, oxidation-reduction reactions and hydrolase activities. DE analysis was validated through Real Time-quantitative PCR of 9 unigenes and 3 isoforms.

Conclusions

Differential expression of sequences involved in signal detection and transduction, cell differentiation and development offered a functional interpretation to the maternally-mediated low naupliar survival rates observed in samples collected on May 23rd. Down-regulation of A5 Obp along with higher quantities of oxylipins-per-litre and oxylipins-per-diatom-cell observed on May 23rd could suggest oxylipin-mediated impairment of naupliar survival in natural populations of T. stylifera. Our results may help identify biomarker genes explaining variations in copepod reproductive responses at a molecular level.

Predation risk alters life history strategies in an oceanic copepod

The ubiquitous oceanic copepod Calanus finmarchicus is the major link between primary producers and important fish stocks in the North Atlantic Ocean and adjacent seas. Despite over a century of research on growth and development of this key species, the effect of predation risk on these processes remains elusive. We tested how food level and chemical cues from a fish predator influence growth and development of C. finmarchicus, using a predator naïve laboratory population. Copepods reached adult stage earlier both in response to high food and to predator cues in our experiment. High food also increased growth and lipid accumulation. In contrast, perceived predation risk triggered reduced size and lipid fullness, indicating a decoupling of growth and development rates. Our results demonstrate that chemical predator cues can influence life history strategies in C. finmarchicus, and suggest that present and future patterns in oceanic zooplankton size and population dynamics may also reflect differences in predation risk.

De Novo Transcriptome Assembly and Gene Expression Profiling of the Copepod Calanus helgolandicus Feeding on the PUA-Producing Diatom Skeletonema marinoi

Diatoms are the dominant component of the marine phytoplankton. Several diatoms produce secondary metabolites, namely oxylipins, with teratogenic effects on their main predators, crustacean copepods. Our study reports the de novo assembled transcriptome of the calanoid copepod Calanus helgolandicus feeding on the oxylipin-producing diatom Skeletonema marinoi. Differential expression analysis was also performed between copepod females exposed to the diatom and the control flagellate Prorocentrum minimum, which does not produce oxylipins. Our results showed that transcripts involved in carbohydrate, amino acid, folate and methionine metabolism, embryogenesis, and response to stimulus were differentially expressed in the two conditions. Expression of 27 selected genes belonging to these functional categories was also analyzed by RT-qPCR in C. helgolandicus females exposed to a mixed solution of the oxylipins heptadienal and octadienal at the concentration of 10 µM, 15 µM, and 20 µM. The results confirmed differential expression analysis, with up-regulation of genes involved in stress response and down-regulation of genes associated with folate and methionine metabolism, embryogenesis, and signaling. Overall, we offer new insights on the mechanism of action of oxylipins on maternally-induced embryo abnormality. Our results may also help identify biomarker genes associated with diatom-related reproductive failure in the natural copepod population at sea.

The β-oxidation pathway is downregulated during diapause termination in Calanus copepods

Calanus copepods are keystone species in marine ecosystems, mainly due to their high lipid content, which is a nutritious food source for e.g. juvenile fish. Accumulated lipids are catabolized to meet energy requirements during dormancy (diapause), which occurs during the last copepodite stage (C5). The current knowledge of lipid degradation pathways during diapause termination is limited. We characterized changes in lipid fullness and generated transcriptional profiles in C5s during termination of diapause and progression towards adulthood. Lipid fullness of C5s declined linearly during developmental progression, but more β-oxidation genes were upregulated in early C5s compared to late C5s and adults. We identified four possible master regulators of energy metabolism, which all were generally upregulated in early C5s, compared to late C5s and adults. We discovered that one of two enzymes in the carnitine shuttle is absent from the calanoid copepod lineage. Based on the geographical location of the sampling site, the field-samples were initially presumed to consist of C. finmarchicus. However, the identification of C. glacialis in some samples underlines the need for performing molecular analyses to reliably identify Calanus species. Our findings contributes to a better understanding of molecular events occurring during diapause and diapause termination in calanoid copepods.

Diapause within the Context of Life-History Strategies in Calanid Copepods (Calanoida: Crustacea)

Post-embryonic diapause in copepods is an adaptation that allows species in the copepod family Calanidae to thrive in high-latitude environments by transforming a short spring phytoplankton bloom into large numbers of lipid-rich individuals capable of surviving a long period of starvation. The copepods, with their high-energy lipid reservoirs, are a critical food source for higher trophic levels, making the Calanidae a key component of high-latitude marine ecosystems. The physiological ecology of the developmental program remains poorly understood. However, new studies using high-throughput RNA sequencing approaches are giving detailed access to physiological status by generating gene expression profiles for both field-collected and laboratory-incubated individuals. These are beginning to characterize the diapause phenotype, elucidate the transcriptional and physiological progression through the diapause program, and illustrate the effects of organism-environment interactions. This paper reviews gene expression profiling studies on the life cycle and diapause program of Neocalanus flemingeri Miller (1988) that were conducted as part of a long-term observation program in the northern Gulf of Alaska. It summarizes recent findings and relates them to the ecology of this species and to that of other calanids.

A Crude Awakening: Effects of Crude Oil on Lipid Metabolism in Calanoid Copepods Terminating Diapause

Calanus finmarchicus and Calanus glacialis are keystone zooplankton species in North Atlantic and Arctic marine ecosystems because they form a link in the trophic transfer of nutritious lipids from phytoplankton to predators on higher trophic levels. These calanoid copepods spend several months of the year in deep waters in a dormant state called diapause, after which they emerge in surface waters to feed and reproduce during the spring phytoplankton bloom. Disruption of diapause timing could have dramatic consequences for marine ecosystems. In the present study, Calanus C5 copepodites were collected in a Norwegian fjord during diapause and were subsequently experimentally exposed to the water-soluble fraction of a naphthenic North Sea crude oil during diapause termination. The copepods were sampled repeatedly while progressing toward adulthood and were analyzed for utilization of lipid stores and for differential expression of genes involved in lipid metabolism. Our results indicate that water-soluble fraction exposure led to a temporary pause in lipid catabolism, suggested by (i) slower utilization of lipid stores in water-soluble fraction-exposed C5 copepodites and (ii) more genes in the β-oxidation pathway being downregulated in water-soluble fraction-exposed C5 copepodites than in the control C5 copepodites. Because lipid content and/or composition may be an important trigger for termination of diapause, our results imply that the timing of diapause termination and subsequent migration to the surface may be delayed if copepods are exposed to oil pollution during diapause or diapause termination. This delay could have detrimental effects on ecosystem dynamics.

Regional heterogeneity impacts gene expression in the subarctic zooplankter Neocalanus flemingeri in the northern Gulf of Alaska

Marine pelagic species are being increasingly challenged by environmental change. Their ability to persist will depend on their capacity for physiological acclimatization. Little is known about limits of physiological plasticity in key species at the base of the food web. Here we investigate the capacity for acclimatization in the copepod Neocalanus flemingeri, which inhabits the Gulf of Alaska, a heterogeneous and highly seasonal environment. RNA-Seq analysis of field-collected pre-adults identified large regional differences in expression of genes involved in metabolic and developmental processes and response to stressors. We found that lipid synthesis genes were up-regulated in individuals from Prince William Sound and down-regulated in the Gulf of Alaska. Up-regulation of lipid catabolic genes in offshore individuals suggests they are experiencing nutritional deficits. The expression differences demonstrate physiological plasticity in response to a steep gradient in food availability. Our transcriptional analysis reveals mechanisms of acclimatization that likely contribute to the observed resilience of this population.

Effects of Nylon Microplastic on Feeding, Lipid Accumulation, and Moulting in a Coldwater Copepod

Microplastic debris is a pervasive environmental contaminant that has the potential to impact the health of biota, although its modes of action remain somewhat unclear. The current study tested the hypothesis that exposure to fibrous and particulate microplastics would alter feeding, impacting on lipid accumulation, and normal development (e.g., growth, moulting) in an ecologically important coldwater copepod Calanus finmarchicus. Preadult copepods were incubated in seawater containing a mixed assemblage of cultured microalgae (control), with the addition of ∼50 microplastics mL^-1 of nylon microplastic granules (10-30 μm) or fibers (10 × 30 μm), which are similar in shape and size to the microalgal prey. The additive chemical profiles showed the presence of stabilizers, lubricants, monomer residues, and byproducts. Prey selectivity was significantly altered in copepods exposed to nylon fibers (ANOVA, P < 0.01) resulting in a nonsignificant 40% decrease in algal ingestion rates (ANOVA, P = 0.07), and copepods exposed to nylon granules showed nonsignificant lipid accumulation (ANOVA, P = 0.62). Both microplastics triggered premature moulting in juvenile copepods (Bernoulli GLM, P < 0.01). Our results emphasize that the shape and chemical profile of a microplastic can influence its bioavailability and toxicity, drawing attention to the importance of using environmentally relevant microplastics and chemically profiling plastics used in toxicity testing.

Towards population genomics in non-model species with large genomes: a case study of the marine zooplankton Calanus finmarchicus

Advances in next-generation sequencing technologies and the development of genome-reduced representation protocols have opened the way to genome-wide population studies in non-model species. However, species with large genomes remain challenging, hampering the development of genomic resources for a number of taxa including marine arthropods. Here, we developed a genome-reduced representation method for the ecologically important marine copepod Calanus finmarchicus (haploid genome size of 6.34 Gbp). We optimized a capture enrichment-based protocol based on 2656 single-copy genes, yielding a total of 154 087 high-quality SNPs in C. finmarchicus including 62 372 in common among the three locations tested. The set of capture probes was also successfully applied to the congeneric C. glacialis. Preliminary analyses of these markers revealed similar levels of genetic diversity between the two Calanus species, while populations of C. glacialis showed stronger genetic structure compared to C. finmarchicus. Using this powerful set of markers, we did not detect any evidence of hybridization between C. finmarchicus and C. glacialis. Finally, we propose a shortened version of our protocol, offering a promising solution for population genomics studies in non-model species with large genomes.

The crustacean ecdysone cassette: A gatekeeper for molt and metamorphosis

Arthropods have long been utilized as models to explore molecular function, and the findings derived from them can be applied throughout metazoa, including as a basis for medical research. This has led to the adoption of many representative insect models beyond Drosophila, as each lends its own unique perspective to questions in endocrinology and genetics. However, non-insect arthropods are yet to be realised for the potential insight they may provide in such studies. The Crustacea are among the most ancient arthropods from which insects descended, comprising a huge variety of life histories and ecological roles. Of the events in a typical crustacean development, metamorphosis is perhaps the most ubiquitous, challenging and highly studied. Despite this, our knowledge of the endocrinology which underpins metamorphosis is rudimentary at best; although several key molecules have been identified and studied in depth, the link between them is quite nebulous and leans heavily on well-explored insect models, which diverged from the Pancrustacea over 450 million years ago. As omics technologies become increasingly accessible, they bring the prospect of explorative molecular research which will allow us to uncover components and pathways unique to crustaceans. This review reconciles known components of crustacean metamorphosis and reflects on our findings in insects to outline a future search space, with focus given to the ecdysone cascade. To expand our knowledge of this ubiquitous endocrine system not only aids in our understanding of crustacean metamorphosis, but also provides a deeper insight into the adaptive capacity of arthropods throughout evolution.

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.

ArHsp40 and ArHsp40-2 contribute to stress tolerance and longevity in Artemia franciscana, but only ArHsp40 influences diapause entry

Embryos of the crustacean Artemia franciscana develop either ovoviviparously or oviparously, yielding swimming larvae (nauplii) or encysted gastrulae (cysts), respectively. Nauplii moult several times and become adults whereas cysts enter diapause, a state of dormancy characterized by exceptionally low metabolism and high stress tolerance. Synthesis of molecular chaperones such as the J-domain proteins ArHsp40 and ArHsp40-2 occurs during embryo development and post-diapause growth of A. franciscana and they influence development and stress tolerance. To further investigate J-domain protein function, ArHsp40 and ArHsp40-2 were each knocked down by RNA interference. Reductions in ArHsp40 and ArHsp40-2 had no effect on adult survival, time to release of cysts and nauplii from females and first-brood size. However, knockdown of both A. franciscana J-domain proteins reduced the longevity and heat tolerance of nauplii, with the loss of ArHsp40 having a greater effect. The knockdown of ArHsp40, but not of ArHsp40-2, caused approximately 50% of cysts to abort diapause entry and hatch without exposure to an exogenous signal such as low temperature and/or desiccation. Cysts lacking ArHsp40 that entered diapause exhibited decreased stress tolerance as did cysts with reduced ArHsp40-2, the latter to a lesser degree. The longevity of nauplii hatching prematurely from cysts was less than for nauplii arising by other means. The results expand our understanding of Hsp40 function in A. franciscana stress tolerance and development, especially during diapause, and they provide the first example of a molecular chaperone that influences diapause entry.

De novo transcriptome assembly and differential gene expression analysis of the calanoid copepod Acartia tonsa exposed to nickel nanoparticles

The calanoid copepod Acartia tonsa is a reference species in standardized ecotoxicology bioassay. Despite this interest, there is a lack of knowledge on molecular responses of A. tonsa to contaminants. We generated a de novo assembled transcriptome of A. tonsa exposed 4 days to 8.5 and 17 mg/L nickel nanoparticles (NiNPs), which have been shown to reduce egg hatching success and larval survival but had no effects on the adults. Aims of our study were to 1) improve the knowledge on the molecular responses of A. tonsa copepod and 2) increase the genomic resources of this copepod for further identification of potential biomarkers of NP exposure. The de novo assembled transcriptome of A. tonsa consisted of 53,619 unigenes, which were further annotated to nr, GO, KOG and KEGG databases. In particular, most unigenes were assigned to Metabolic and Cellular processes (34-45%) GO terms, and to Human disease (28%) and Organismal systems (23%) KEGG categories. Comparison among treatments showed that 373 unigenes were differentially expressed in A. tonsa exposed to NiNPs at 8.5 and 17 mg/L, with respect to control. Most of these genes were downregulated and took part in ribosome biogenesis, translation and protein turnover, thus suggesting that NiNPs could affect the copepod ribosome synthesis machinery and functioning. Overall, our study highlights the potential of toxicogenomic approach in gaining more mechanistic and functional information about the mode of action of emerging compounds on marine organisms, for biomarker discovering in crustaceans.

Variations in Copepod Proteome and Respiration Rate in Association with Diel Vertical Migration and Circadian Cycle

The diel vertical migration of zooplankton is a process during which individuals spend the night in surface waters and retreat to depth during the daytime, with substantial implications for carbon transport and the ecology of midwater ecosystems. The physiological consequences of this daily pattern have, however, been poorly studied beyond investigations of speed and the energetic cost of swimming. Many other processes are likely influenced, such as fuel use, energetic trade-offs, underlying diel (circadian) rhythms, and antioxidant responses. Using a new reference transcriptome, proteomic analyses were applied to compare the physiological state of a migratory copepod, Pleuromamma xiphias, immediately after arriving to the surface at night and six hours later. Oxygen consumption was monitored semi-continuously to explore underlying cyclical patterns in metabolic rate under dark-dark conditions. The proteomic analysis suggests a distinct shift in physiology that reflects migratory exertion and changes in metabolism. These proteomic analyses are supported by the respiration experiments, which show an underlying cycle in metabolic rate, with a peak at dawn. This project generates molecular tools (transcriptome and proteome) that will allow for more detailed understanding of the underlying physiological processes that influence and are influenced by diel vertical migration. Further, these studies suggest that P. xiphias is a tractable model for continuing investigations of circadian and diel vertical migration influences on plankton physiology. Previous studies did not account for this cyclic pattern of respiration and may therefore have unrepresented respiratory carbon fluxes from copepods by about 24%.

Inter-individual variability in copepod microbiomes reveals bacterial networks linked to host physiology

Copepods harbor diverse bacterial communities, which collectively carry out key biogeochemical transformations in the ocean. However, bulk copepod sampling averages over the variability in their associated bacterial communities, thereby limiting our understanding of the nature and specificity of copepod-bacteria associations. Here, we characterize the bacterial communities associated with nearly 200 individual Calanus finmarchicus copepods transitioning from active growth to diapause. We find that all individual copepods sampled share a small set of "core" operational taxonomic units (OTUs), a subset of which have also been found associated with other marine copepod species in different geographic locations. However, most OTUs are patchily distributed across individual copepods, thereby driving community differences across individuals. Among patchily distributed OTUs, we identified groups of OTUs correlated with common ecological drivers. For instance, a group of OTUs positively correlated with recent copepod feeding served to differentiate largely active growing copepods from those entering diapause. Together, our results underscore the power of individual-level sampling for understanding host-microbiome relationships.

Complementary mechanisms for neurotoxin resistance in a copepod

Toxin resistance is a recurring evolutionary response by predators feeding on toxic prey. These adaptations impact physiological interaction and community ecology. Mechanisms for resistance vary depending on the predator and the nature of the toxin. Potent neurotoxins like tetrodotoxin (TTX) and saxitoxin (STX) that are highly toxic to humans and other vertebrates, target conserved voltage-gated sodium channels (Na_V) of nerve and muscle, causing paralysis. The copepod Calanus finmarchicus consumes the STX-producing dinoflagellate, Alexandrium fundyense with no effect on survival. Using transcriptomic approaches to search for the mechanism that confers resistance in C. finmarchicus, we identified splice variants of Na_Vs that were predicted to be toxin resistant. These were co-expressed with putatively non-resistant form in all developmental stages. However its expression was unresponsive to toxin challenge nor was there any up-regulation of genes involved in multi-xenobiotic resistance (MXR) or detoxification (phases I or II). Instead, adults consistently regulated genes encoding digestive enzymes, possibly to complement channel resistance by limiting toxin assimilation via the digestive process. The nauplii, which were more susceptible to STX, did not regulate these enzymes. This study demonstrates how deep-sequencing technology can elucidate multiple mechanisms of toxin resistance concurrently, revealing the linkages between molecular/cellular adaptations and the ecology of an organism.

A deep transcriptomic resource for the copepod crustacean Labidocera madurae: A potential indicator species for assessing near shore ecosystem health

Coral reef ecosystems of many sub-tropical and tropical marine coastal environments have suffered significant degradation from anthropogenic sources. Research to inform management strategies that mitigate stressors and promote a healthy ecosystem has focused on the ecology and physiology of coral reefs and associated organisms. Few studies focus on the surrounding pelagic communities, which are equally important to ecosystem function. Zooplankton, often dominated by small crustaceans such as copepods, is an important food source for invertebrates and fishes, especially larval fishes. The reef-associated zooplankton includes a sub-neustonic copepod family that could serve as an indicator species for the community. Here, we describe the generation of a de novo transcriptome for one such copepod, Labidocera madurae, a pontellid from an intensively-studied coral reef ecosystem, Kāne'ohe Bay, Oahu, Hawai'i. The transcriptome was assembled using high-throughput sequence data obtained from whole organisms. It comprised 211,002 unique transcripts, including 72,391 with coding regions. It was assessed for quality and completeness using multiple workflows. Bench-marking-universal-single-copy-orthologs (BUSCO) analysis identified transcripts for 88% of expected eukaryotic core proteins. Targeted gene-discovery analyses included searches for transcripts coding full-length "giant" proteins (>4,000 amino acids), proteins and splice variants of voltage-gated sodium channels, and proteins involved in the circadian signaling pathway. Four different reference transcriptomes were generated and compared for the detection of differential gene expression between copepodites and adult females; 6,229 genes were consistently identified as differentially expressed between the two regardless of reference. Automated bioinformatics analyses and targeted manual gene curation suggest that the de novo assembled L. madurae transcriptome is of high quality and completeness. This transcriptome provides a new resource for assessing the global physiological status of a planktonic species inhabiting a coral reef ecosystem that is subjected to multiple anthropogenic stressors. The workflows provide a template for generating and assessing transcriptomes in other non-model species.

A transcriptomic resource for the northern krill Meganyctiphanes norvegica based on a short-term temperature exposure experiment

The northern krill, Meganyctiphanes norvegica, is an important component of the pelagic food web across the North Atlantic. Widespread from the Mediterranean to the Subarctic Atlantic, populations appear to be strongly adapted to local temperatures, and seem to have very little plasticity. The goal of this study was to create and annotate a de novo transcriptome assembly to allow for comparative and physiological studies and to explore the gene expression response of M. norvegica from the Gulf of Maine to two different temperature conditions. Our Trinity assembly produced 405,497 transcripts with ~16% annotation success versus nr with a stringent cutoff (>1e^-10), and substantial cross-annotation versus FlyBase and other published pelagic crustacean transcriptomes. There were 122 transcripts that were differentially expressed based on our 2-day 9 versus 12°C temperature exposure, and their annotation suggested changes in energetic metabolism and molting. These results generate a useful molecular resource for further more directed studies as well as provide initial insight into the physiological processes that may shape the temperature response of the northern krill.

Circadian signaling in the Northern krill Meganyctiphanes norvegica: In silico prediction of the protein components of a putative clock system using a publicly accessible transcriptome

The Northern krill Meganyctiphanes norvegica is a significant component of the zooplankton community in many regions of the North Atlantic Ocean. In the areas it inhabits, M. norvegica is of great importance ecologically, as it is both a major consumer of phytoplankton/small zooplankton and is a primary food source for higher-level consumers. One behavior of significance for both feeding and predator avoidance in Meganyctiphanes is diel vertical migration (DVM), i.e., a rising from depth at dusk and a return to depth at dawn. In this and other euphausiids, an endogenous circadian pacemaker is thought, at least in part, to control DVM. Currently, there is no information concerning the identity of the genes/proteins that comprise the M. norvegica circadian system. In fact, there is little information concerning the molecular underpinnings of circadian rhythmicity in crustaceans generally. Here, a publicly accessible transcriptome was used to identify the molecular components of a putative Meganyctiphanes circadian system. A complete set of core clock proteins was deduced from the M. norvegica transcriptome (clock, cryptochrome 2, cycle, period and timeless), as was a large suite of proteins that likely function as modulators of the core clock (e.g., doubletime), or serves as inputs to it (cryptochrome 1) or outputs from it (pigment dispersing hormone). This is the first description of a "complete" (core clock through putative output pathway signals) euphausiid clock system, and as such, provides a foundation for initiating molecular investigations of circadian signaling in M. norvegica and other krill species, including how clock systems may regulate DVM and other behaviors.

De novo transcriptome assembly of the calanoid copepod Neocalanus flemingeri: A new resource for emergence from diapause

Copepods, small planktonic crustaceans, are key links between primary producers and upper trophic levels, including many economically important fishes. In the subarctic North Pacific, the life cycle of copepods like Neocalanus flemingeri includes an ontogenetic migration to depth followed by a period of diapause (a type of dormancy) characterized by arrested development and low metabolic activity. The end of diapause is marked by the production of the first brood of eggs. Recent temperature anomalies in the North Pacific have raised concerns about potential negative effects on N. flemingeri. Since diapause is a developmental program, its progress can be tracked using through global gene expression. Thus, a reference transcriptome was developed as a first step towards physiological profiling of diapausing females using high-throughput Illumina sequencing. The de novo transcriptome, the first for this species was designed to investigate the diapause period. RNA-Seq reads were obtained for dormant to reproductive N. flemingeri females. A high quality de novo transcriptome was obtained by first assembling reads from each individual using Trinity software followed by clustering with CAP3 Assembly Program. This assembly consisted of 140,841transcripts (contigs). Bench-marking universal single-copy orthologs analysis identified 85% of core eukaryotic genes, with 79% predicted to be complete. Comparison with other calanoid transcriptomes confirmed its quality and degree of completeness. Trinity assembly of reads originating from multiple individuals led to fragmentation. Thus, the workflow applied here differed from the one recommended by Trinity, but was required to obtain a good assembly.

Mortality and transcriptional effects of inorganic mercury in the marine copepod Calanus finmarchicus

Inorganic mercury (Hg) is highly toxic to organisms including crustaceans and displays multiple toxic modes of action (MoA). The main aim of this investigation was to assess the acute and sublethal toxicity mediated by mercury chloride (HgCl₂) in the marine copepod Calanus finmarchicus. A combination of short-term static studies to determine acute toxicity and a transcriptional investigation to characterize the sublethal MoA of HgCl₂ were conducted with an in-house continuous culture of C. finmarchicus. Transcriptional changes were determined by a custom 6.6 k C. finmarchicus Agilent oligonucleotide microarray and quantitative RT-PCR analysis. Data demonstrate that HgCl₂ produced a concentration- and time-dependent reduction in survival (NOEC_{48 h} = 6.9 μg/L [Hg²⁺] and LC₅₀ of 279, 73, 48, and 34 µg/L [Hg²⁺] after 24, 48, 72, and 96 h, respectively) and that exposure to sublethal concentrations of HgCl₂ (5 μg/L [Hg2⁺]) induced differential expression of 98 features (probes) on the microarray. Gene ontology (GO) and toxicological pathway analyses suggested that the main MOA were (1) uncoupling of mitochondrial oxidative phosphorylation (OXPHOS) and ATP production, (2) oxidative stress and macromolecular damage, (3) inactivation of cellular enzymes, (4) induction of cellular apoptosis and autophagocytosis, (5) over-excitation of glutamate receptors (neurotoxicity), (6) disruption of calcium homeostasis and signaling, and (7) modulation of nuclear receptor activity involved in vitamin D receptor signaling. Quantitative RT-PCR analysis verified that oligoarray performed reliably in terms of specificity and response, thus demonstrating that Hg²⁺ exerts multiple potential MoA in C. finmarchicus.

Maternal polycyclic aromatic hydrocarbon (PAH) transfer and effects on offspring of copepods exposed to dispersed oil with and without oil droplets

Copepods of the genus Calanus have the potential for accumulating lipophilic oil components due to their high lipid content and found to filter and ingest oil droplets during exposure. As female copepods produce eggs at the expense of lipid storage, there is a concern for transfer of lipophilic contaminants to offspring. To assess the potential for maternal transfer of oil components, ovigerous female copepods (Calanus finmarchicus) were exposed to filtered and unfiltered oil dispersions for 4 days, collected and eggs maintained in clean seawater and hatching and gene expression examined in hatched nauplii. Oil droplet exposure contributed to polycyclic aromatic hydrocarbon (PAH) uptake in dispersion-treated adult copepods, as displayed through PAH body residue analyses and fluorescence microscopy. Applying the latter methodology, transfer of heavy PAH from copepod mothers to offspring were detected Subtle effects were observed in offspring as evidenced by a temporal reduction in hatching success appear to be occurring only when mothers were exposed to the unfiltered oil dispersions. Offspring reared in clean water through to late naupliar stages were collected for RNA extraction and preparation of libraries for high-throughput transcriptome sequencing. Differentially expressed genes were identified through pairwise comparisons between treatments. Among these, several expressed genes have known roles in responses to chemical stress including xenobiotic metabolism enzymes, antioxidants, chaperones, and components of the inflammatory response. While gene expression results suggest a transgenerational activation of stress responses, the increase in relatively small number of differentially expressed genes suggests a minor long-term effect on offspring following maternal exposure.

Molecular Characterization of Copepod Photoreception

Copepod crustaceans are an abundant and ecologically significant group whose basic biology is guided by numerous visually guided behaviors. These behaviors are driven by copepod eyes, including naupliar eyes and Gicklhorn's organs, which vary widely in structure and function among species. Yet little is known about the molecular aspects of copepod vision. In this study we present a general overview of the molecular aspects of copepod vision by identifying phototransduction genes from newly generated and publicly available RNA-sequencing data and assemblies from 12 taxonomically diverse copepod species. We identify a set of 10 expressed transcripts that serve as a set of target genes for future studies of copepod phototransduction. Our more detailed evolutionary analyses of the opsin gene responsible for forming visual pigments found that all of the copepod species investigated express two main groups of opsins: middle-wavelength-sensitive (MWS) opsins and pteropsins. Additionally, there is evidence from a few species (e.g., Calanus finmarchicus, Eurytemora affinis, Paracyclopina nana, and Lernaea cyprinacea) for the expression of two additional groups of opsins-the peropsins and rhodopsin 7 (Rh7) opsins-at low levels or distinct developmental stages. An ontogenetic analysis of opsin expression in Calanus finmarchicus found the expression of a single dominant MWS opsin, as well as evidence for differences in expression across development in some MWS, pteropsin, and Rh7 opsins, with expression peaking in early naupliar through early copepodite stages.

Spliced leader-based analyses reveal the effects of polycyclic aromatic hydrocarbons on gene expression in the copepod Pseudodiaptomus poplesia

Polycyclic aromatic hydrocarbons (PAHs) are a group of toxic and carcinogenic pollutants that can adversely affect the development, growth and reproduction of marine organisms including copepods. However, knowledge on the molecular mechanisms regulating the response to PAH exposure in marine planktonic copepods is limited. In this study, we investigated the survival and gene expression of the calanoid copepod Pseudodiaptomus poplesia upon exposure to two PAHs, 1, 2-dimethylnaphthalene (1, 2-NAPH) and pyrene. Acute toxicity responses resulted in 96-h LC₅₀ of 788.98μgL^-1 and 54.68μgL^-1 for 1, 2-NAPH and pyrene, respectively. Using the recently discovered copepod spliced leader as a primer, we constructed full-length cDNA libraries from copepods exposed to sublethal concentrations and revealed 289 unique genes of diverse functions, including stress response genes and novel genes previously undocumented for this species. Eighty-three gene families were specifically expressed in PAH exposure libraries. We further analyzed the expression of seven target genes by reverse transcription-quantitative PCR in a time-course test with three sublethal concentrations. These target genes have primary roles in detoxification, oxidative defense, and signal transduction, and include different forms of glutathione S-transferase (GST), glutathione peroxidases (GPX), peroxiredoxin (PRDX), methylmalonate-semialdehyde dehydrogenase (MSDH) and ras-related C3 botulinum toxin substrate (RAC1). Expression stability of seven candidate reference genes were evaluated and the two most stable ones (RPL15 and RPS20 for 1, 2-NAPH exposure, RPL15 and EF1D for pyrene exposure) were used to normalize the expression levels of the target genes. Significant upregulation was detected in GST-T, GST-DE, GPX4, PRDX6 and RAC1 upon 1, 2-NAPH exposure, and GST-DE and MSDH upon pyrene exposure. These results indicated that the oxidative stress was induced and that signal transduction might be affected by PAH exposure in P. poplesia. However, gene upregulation was followed by a reduction in expression level towards 96h, indicating a threshold value of exposure time that leads to depressed gene expression. Prolonged exposure may cause dysfunction of detoxification and antioxidant machinery in P. poplesia. The transcriptional responses of GST-T, GPX2 and GPX4 upon pyrene exposure were minimal. Our results reveal the different sensitivity of P. poplesia to two PAHs at both the individual and transcriptional levels. As the first attempt, this study proved that copepod spliced leader is useful for obtaining full-length cDNA in P. poplesia exposed to PAHs and provided a valuable gene resource for this non-model species. This approach can be applied to other calanoid copepods exposed to various stressors, particularly under field conditions.

Mechanisms of Action of Compounds That Enhance Storage Lipid Accumulation in Daphnia magna

Accumulation of storage lipids in the crustacean Daphnia magna can be altered by a number of exogenous and endogenous compounds, like 20-hydroxyecdysone (natural ligand of the ecdysone receptor, EcR), methyl farnesoate, pyrirproxyfen (agonists of the methyl farnesoate receptor, MfR), and tributyltin (agonist of the retinoid X acid receptor, RXR). This effect, analogous to the obesogenic disruption in mammals, alters Daphnia's growth and reproductive investment. Here we propose that storage lipid accumulation in droplets is regulated in Daphnia by the interaction between the nuclear receptor heterodimer EcR:RXR and MfR. The model was tested by determining changes in storage lipid accumulation and on gene transcription in animals exposed to different effectors of RXR, EcR, and MfR signaling pathways, either individually or in combination. RXR, EcR, and MfR agonists increased storage lipid accumulation, whereas fenarimol and testosterone (reported inhibitors of ecdysteroid synthesis and an EcR antagonist, respectively) decreased it. Joint effects of mixtures with fenarimol, testosterone, and ecdysone were antagonistic, mixtures of juvenoids showed additive effects following a concentration addition model, and combinations of tributyltin with juvenoids resulted in greater than additive effects. Co-exposures of ecdysone with juvenoids resulted in deregulation of ecdysone- and farnesoid-regulated genes, accordingly with the observed changes in lipid accumulation These results indicate the requirement of ecdysone binding to the EcR:RXR:MfR complex to regulate lipid storage and that an excess of ecdysone disrupts the whole process, probably by triggering negative feedback mechanisms.

Effects of Elevated Temperature and Food Supply on the Termination of Over-Summering and Subsequent Development of the Calanoid Copepod Calanus sinicus: Morphology, Physiology and Gene Expression

The copepod Calanus sinicus Brodsky dominates the zooplankton in the Yellow Sea, China, and undergoes over-summering within the Yellow Sea Cold Water Mass (YSCWM). Termination of over-summering and subsequent development are regarded as key processes in population recruitment, and are probably linked to environmental variations in the YSCWM. In this study, we examined the effects of temperature (9 and 18°C) and food conditions (0.1 μg C mL^-1 and unfed) on metabolic rates, morphological characteristics, and relative gene expressions of six genes involved in molting, gonad development, lipid catabolism, and stress tolerance processes of C. sinicus during termination of over-summering and subsequent development. Both elevated temperature and external food supply rapidly ended over-summering of C. sinicus, accompanied by up-regulation of the ecdysteroid receptor (EcR) gene expression and increased metabolic rates. These environmental conditions resulted in irreversible termination of over-summering and ensure the success of molting. During subsequent development, the lipid reserve in oil sacs could permit only early gonad development. The food supply might be a trigger to activate the final maturity of gonad by up-regulating expression of the vitellogenin receptor (VgR) gene. Thus, food played an indispensable role in population recruitment after termination of over-summering, whereas the elevated temperature accelerated these physiological processes. This study revealed the first dynamic profiles of physiological processes involved in over-summering termination and the subsequent development of C. sinicus using morphological, physiological and molecular methods simultaneously, confirmed the quiescent state of over-summering C5 copepodites, detected the effects of environmental changes on over-summering termination and subsequent development, and provided a foundation for future investigations of the mechanisms involved in over-summering in YSCWM.

Diversity of insulin-like peptide signaling system proteins in Calanus finmarchicus (Crustacea; Copepoda) - Possible contributors to seasonal pre-adult diapause

Calanus finmarchicus, an abundant calanoid copepod in the North Atlantic Ocean, is both a major grazer on phytoplankton and an important forage species for invertebrate and vertebrate predators. One component of the life history of C. finmarchicus is the overwintering dormancy of sub-adults, a feature key for the annual recruitment of this species in early spring. While little is known about the control of dormancy in C. finmarchicus, one hypothesis is that it is an insect-like diapause, where the endocrine system is a key regulator. One group of hormones implicated in the control of insect diapause is the insulin-like peptides (ILPs). Here, C. finmarchicus transcriptomic data were used to predict ILP signaling pathway proteins. Four ILP precursors were identified, each possessing a distinct A- and B-chain peptide; these peptides are predicted to form bioactive heterodimers via inter-chain disulfide bridging. Two ILP receptors, which likely represent splice variants of a common gene, were identified. Three insulin-degrading enzymes were also discovered, as were proteins encoding the transcription factor FOXO, a downstream target of ILP that has been implicated in the regulation of insect diapause, and insulin receptor substrate, a protein putatively linking the ILP receptor and FOXO. RNA-Seq data suggest that some C. finmarchicus insulin pathway transcripts are differentially expressed across development. As in insects, the ILP signaling system may be involved in controlling C. finmarchicus' organism-environment interactions (e.g., regulation of seasonal sub-adult diapause), a hypothesis that can now be investigated using these data.

Construction and Characterization of Two Novel Transcriptome Assemblies in the Congeneric Porcelain Crabs Petrolisthes cinctipes and P. manimaculis

Crustaceans have commonly been used as non-model systems in basic biological research, especially physiological regulation. With the recent and rapid adoption of functional genomic tools, crustaceans are increasingly becoming model systems for ecological investigations of development and evolution and for mechanistic examinations of genotype-phenotype interactions and molecular pathways of response to environmental stressors. Comparative transcriptomic approaches, however, remain constrained by a lack of sequence data in closely related crustacean taxa. We identify challenges in the use of functional genomics tools in comparative analysis among decapod crustacean in light of recent advances. We present RNA-seq data from two congeneric species of porcelain crabs (Petrolisthes cinctipes and P. manimaculis) used to construct two de novo transcriptome assemblies with ∼194K and ∼278K contigs, respectively. We characterize and contrast these assemblies and compare them to a previously generated EST sequence library for P. cinctipes We also discuss the potential use of these data as a case-study system in the broader context of crustacean comparative transcriptomics.

Massive Shift in Gene Expression during Transitions between Developmental Stages of the Gall Midge, Mayetiola Destructor

Mayetiola destructor is a destructive pest of wheat and has six developmental stages. Molecular mechanisms controlling the transition between developmental stages remain unknown. Here we analyzed genes that were expressed differentially between two successive developmental stages, including larvae at 1, 3, 5, and 7 days, pupae, and adults. A total of 17,344 genes were expressed during one or more of these studied stages. Among the expressed genes, 38-68% were differently expressed between two successive stages, with roughly equal percentages of up- and down-regulated genes. Analysis of the functions of the differentially expressed genes revealed that each developmental stage had some unique types of expressed genes that are characteristic of the physiology at that stage. This is the first genome-wide analysis of genes differentially expressed in different stages in a gall midge. The large dataset of up- and down-regulated genes in each stage of the insect shall be very useful for future research to elucidate mechanisms regulating insect development and other biological processes.

Acute toxicity of dispersed crude oil on the cold-water copepod Calanus finmarchicus: Elusive implications of lipid content

In this investigation, acute toxicity data were used from two previously reported studies where cold-water copepods were exposed to mechanically dispersed (MD) and chemically (CD) dispersed oil. In one of these studies, concentration-dependent mortality was observed, whereas no apparent relationship between exposure concentration and mortality was found in the other. The only marked difference between the studies is that copepods in the first experiment displayed a lower lipid sac volume (on average) than in the second one. In this study additional biometric data on lipid content were utilized and observed effects and toxicokinetics modeling applied in order to investigate whether differences in sensitivity between copepod cohorts might be explained by differences in lipid content. Results suggest that although a considerable lipid sac might retard toxicokinetics, the observed differences in lipid volume are not sufficient to explain differences in toxicity. Further, there are no apparent indications that acute toxic stress leads to lipid depletion, or that acute increased mortality rate selectively affects lipid-poor individuals. It is conceivable that other potential explanations exist, but the causal relationship between lipid content and increased mortality frequency remains elusive.