Transcriptomic responses of the calanoid copepod Calanus finmarchicus to the saxitoxin producing dinoflagellate Alexandrium fundyense

A short-term exposure to saxitoxin triggers a multitude of deleterious effects in Daphnia magna at levels deemed safe for human health

Harmful algal blooms and the toxins produced during these events are a human and environmental health concern worldwide. Saxitoxin and its derivatives are potent natural aquatic neurotoxins produced by certain freshwater cyanobacteria and marine algae species during these bloom events. Saxitoxins effects on human health are well studied, however its effects on aquatic biota are still largely unexplored. This work aims at evaluating the effects of a pulse acute exposure (24 h) of the model cladoceran Daphnia magna to 30 μg saxitoxin L-1, which corresponds to the safety guideline established by the World Health Organization (WHO) for these toxins in recreational freshwaters. Saxitoxin effects were assessed through a comprehensive array of biochemical (antioxidant enzymes activity and lipid peroxidation), genotoxicity (alkaline comet assay), neurotoxicity (total cholinesterases activity), behavioral (swimming patterns), physiological (feeding rate and heart rate), and epigenetic (total 5-mC DNA methylation) biomarkers. Exposure resulted in decreased feeding rate, heart rate, total cholinesterases activity and catalase activity. Contrarily, other antioxidant enzymes, namely glutathione-S-transferases and selenium-dependent Glutathione peroxidase had their activity increased, together with lipid peroxidation levels. The enhancement of the antioxidant enzymes was not sufficient to prevent oxidative damage, as underpinned by lipid peroxidation enhancement. Accordingly, average DNA damage level was significantly increased in STX-exposed daphnids. Total DNA 5-mC level was significantly decreased in exposed organisms. Results showed that even a short-term exposure to saxitoxin causes significant effects on critical molecular and cellular pathways and modulates swimming patterns in D. magna individuals. This study highlights sub-lethal effects caused by saxitoxin in D. magna, suggesting that these toxins may represent a marked challenge to their thriving even at a concentration deemed safe for humans by the WHO.

Transcriptome response of a marine copepod in response to environmentally-relevant concentrations of saxitoxin

Paralytic shellfish toxins (PSTs) can pose a serious threat to human health. Among them, saxitoxin (STX) is one of the most potent natural neurotoxins. Here, the copepod Tigriopus japonicus, was exposed to environmentally relevant concentrations (2.5 and 25 μg/L) STX for 48 h. Although no lethal effects were observed at both concentrations, the transcriptome was significantly altered, and displayed a concentration-dependent response. STX exposure decreased the copepod's metabolism and compromised immune defense and detoxification. Additionally, STX disturbed signal transduction, which might affect other cellular processes. STX exposure could inhibit the copepod's chitin metabolism, disrupting its molting process. Also, the processes related to damage repair and protection were up-regulated to fight against high concentration exposure. Collectively, this study has provided an early warning of PSTs for coastal ecosystem not only because of their potent toxicity effect but also their bioaccumulation that can transfer up the food chain after ingestion by copepods.

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.

A review on aquatic toxins - Do we really know it all regarding the environmental risk posed by phytoplankton neurotoxins?

Aquatic toxins are potent natural toxins produced by certain cyanobacteria and marine algae species during harmful cyanobacterial and algal blooms (CyanoHABs and HABs, respectively). These harmful bloom events and the toxins produced during these events are a human and environmental health concern worldwide, with occurrence, frequency and severity of CyanoHABs and HABs being predicted to keep increasing due to ongoing climate change scenarios. These contexts, as well as human health consequences of some toxins produced during bloom events have been thoroughly reviewed before. Conversely, the wider picture that includes the non-human biota in the assessment of noxious effects of toxins is much less covered in the literature and barely covered by review works. Despite direct human exposure to aquatic toxins and related deleterious effects being responsible for the majority of the public attention to the blooms' problematic, it constitutes a very limited fraction of the real environmental risk posed by these toxins. The disruption of ecological and trophic interactions caused by these toxins in the aquatic biota building on deleterious effects they may induce in different species is paramount as a modulator of the overall magnitude of the environmental risk potentially involved, thus necessarily constraining the quality and efficiency of the management strategies that should be placed. In this way, this review aims at updating and consolidating current knowledge regarding the adverse effects of aquatic toxins, attempting to going beyond their main toxicity pathways in human and related models' health, i.e., also focusing on ecologically relevant model organisms. For conciseness and considering the severity in terms of documented human health risks as a reference, we restricted the detailed revision work to neurotoxic cyanotoxins and marine toxins. This comprehensive revision of the systemic effects of aquatic neurotoxins provides a broad overview of the exposure and the hazard that these compounds pose to human and environmental health. Regulatory approaches they are given worldwide, as well as (eco)toxicity data available were hence thoroughly reviewed. Critical research gaps were identified particularly regarding (i) the toxic effects other than those typical of the recognized disease/disorder each toxin causes following acute exposure in humans and also in other biota; and (ii) alternative detection tools capable of being early-warning signals for aquatic toxins occurrence and therefore provide better human and environmental safety insurance. Future directions on aquatic toxins research are discussed in face of the existent knowledge, with particular emphasis on the much-needed development and implementation of effective alternative (eco)toxicological biomarkers for these toxins. The wide-spanning approach followed herein will hopefully stimulate future research more broadly addressing the environmental hazardous potential of aquatic toxins.

De novo transcriptomes of six calanoid copepods (Crustacea): a resource for the discovery of novel genes

This study presents eight new high-quality de novo transcriptomes from six co-occurring species of calanoid copepods, the first published for Neocalanus plumchrus, N. cristatus, Eucalanus bungii and Metridia pacifica and additional ones for N. flemingeri and Calanus marshallae. They are ecologically-important members of sub-arctic North Pacific marine zooplankton communities. 'Omics data for this diverse and numerous taxonomic group are sparse and difficult to obtain. Total RNA from single individuals was used to construct gene libraries that were sequenced on an Illumina Next-Seq platform. Quality filtered reads were assembled with Trinity software and validated using multiple criteria. The study's primary purpose is to provide a resource for gene expression studies. The integrated database can be used for quantitative inter- and intra-species comparisons of gene expression patterns across biological processes. An example of an additional use is provided for discovering novel and evolutionarily-significant proteins within the Calanoida. A workflow was designed to find and characterize unannotated transcripts with homologies across de novo assemblies that have also been shown to be eco-responsive.

Chemosensory-Related Genes in Marine Copepods

Living organisms deeply rely on the acquisition of chemical signals in any aspect of their life, from searching for food, mating and defending themselves from stressors. Copepods, the most abundant and ubiquitous metazoans on Earth, possess diversified and highly specified chemoreceptive structures along their body. The detection of chemical stimuli activates specific pathways, although this process has so far been analyzed only on a relatively limited number of species. Here, in silico mining of 18 publicly available transcriptomes is performed to delve into the copepod chemosensory genes, improving current knowledge on the diversity of this multigene family and on possible physiological mechanisms involved in the detection and analysis of chemical cues. Our study identifies the presence of ionotropic receptors, chemosensory proteins and gustatory receptors in copepods belonging to the Calanoida, Cyclopoida and Harpacticoida orders. We also confirm the absence in these copepods of odorant receptors and odorant-binding proteins agreeing with their insect specificity. Copepods have evolved several mechanisms to survive in the harsh marine environment such as producing proteins to respond to external stimulii. Overall, the results of our study open new possibilities for the use of the chemosensory genes as biomarkers in chemical ecology studies on copepods and possibly also in other marine holozooplankters.

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.

Marine Neurotoxins' Effects on Environmental and Human Health: An OMICS Overview

Harmful algal blooms (HAB), and the consequent release of toxic metabolites, can be responsible for seafood poisoning outbreaks. Marine wildlife can accumulate these toxins throughout the food chain, which presents a threat to consumers’ health. Some of these toxins, such as saxitoxin (STX), domoic acid (DA), ciguatoxin (CTX), brevetoxin (BTX), tetrodotoxin (TTX), and β-N-methylamino-L-alanine (BMAA), cause severe neurological symptoms in humans. Considerable information is missing, however, notably the consequences of toxin exposures on changes in gene expression, protein profile, and metabolic pathways. This information could lead to understanding the consequence of marine neurotoxin exposure in aquatic organisms and humans. Nevertheless, recent contributions to the knowledge of neurotoxins arise from OMICS-based research, such as genomics, transcriptomics, proteomics, and metabolomics. This review presents a comprehensive overview of the most recent research and of the available solutions to explore OMICS datasets in order to identify new features in terms of ecotoxicology, food safety, and human health. In addition, future perspectives in OMICS studies are discussed.

First Report of OvoA Gene in Marine Arthropods: A New Candidate Stress Biomarker in Copepods

Ovothiol is one of the most powerful antioxidants acting in marine organisms as a defense against oxidative stress during development and in response to environmental cues. The gene involved in the ovothiol biosynthesis, OvoA, is found in almost all metazoans, but open questions existed on its presence among arthropods. Here, using an in silico workflow, we report a single OvoA gene in marine arthropods including copepods, decapods, and amphipods. Phylogenetic analyses indicated that OvoA from marine arthropods separated from the other marine phyla (e.g., Porifera, Mollusca) and divided into two separate branches, suggesting a possible divergence through evolution. In the copepod Calanus finmarchicus, we suggest that OvoA has a defense role in oxidative stress as shown by its high expression in response to a toxic diet and during the copepodite stage, a developmental stage that includes significant morphological changes. Overall, the results of our study open possibilities for the use of OvoA as a biomarker of stress in copepods and possibly also for other marine holozooplankters. The finding of OvoA in copepods is also promising for the drug discovery field, suggesting the possibility of using copepods as a new source of bioactive compounds to be tested in the marine biotechnological sector.

The genome of the European estuarine calanoid copepod Eurytemora affinis: Potential use in molecular ecotoxicology

In this study, we sequenced and assembled the genome of a European estuarine calanoid copepod using Oxford Nanopore PromethION and Illumina HiSeq 2500 platforms. The length of the assembled genome was 776.1 Mb with N50 = 474.9 kb (BUSCO 85.9%), and the genome consisted of 2473 contigs. A total of 18,014 genes were annotated and orthologous gene clusters were analyzed in comparison to other copepods. In addition, genome-wide identification of cytochrome P450s, glutathione S-transferases, and ATP-binding cassette transporters in E. affinis was performed to determine gene repertoire of these detoxification-related gene families. Results revealed the presence of species-specific gene inventories, indicating that these gene families have evolved through species-specific gene loss/expansion processes, possibly due to adaptation to different environmental stressors. Our study provides a new inventory of the European estuarine calanoid copepod E. affinis genome with emphasis on phase I, II, and III detoxification systems.

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.

Differential gene expression profile of male and female copepods in response to cadmium exposure

In this study, the whole transcriptome and sex-specific differential gene expression of the copepod Pseudodiaptomus annandalei exposed to cadmium (Cd) were investigated. P. annandalei were exposed to 40 μg/L Cd from the naupliar stage to male and female adults. High-throughput transcriptome sequencing (RNA-seq) was performed with copepod samples using an Illumina Hiseq™ 2000 platform. TransDecoder analysis found 32,625 putative open reading frame contigs. At p-values of <0.001, a total of 4756 differentially expressed genes (DEGs) (2216 up-regulated and 2540 down-regulated genes) were found in male copepods. Whereas a total of 2879 DEGs (2007 up-regulated and 872 down-regulated genes) were found in female copepods. A few selected cellular stress response genes, involved in xenobiotic metabolism, energy metabolism, growth, and development as a result of Cd exposure in the copepods were discussed. The study showed that most of these processes were changed in a sex-specific manner, accounting for the different sensitivities of male and female copepods. Results suggest and reinforce that sex is an important factor to be considered in ecotoxicogenomics.

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.

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.

Transcriptomic responses of Artemia salina exposed to an environmentally relevant dose of Alexandrium minutum cells or Gonyautoxin2/3

Toxicities of the marine algae Alexandrium minutum and its excreted gonyautoxins (GTXs) to the marine crustacean Artemia salina were investigated. Mortality was observed for neither larvae nor adult A. salina exposed to A. minutum at a density of 5000 cells/mL or 0.5 μM GTX2/3. After exposure, the full transcriptome of adult A. salina was assembled and functionally annotated. A total of 599,286 transcripts were obtained, which were clustered into 515,196 unigenes. Results of the transcriptional effect level index revealed that direct exposure to the toxic algae A. minutum caused greater alterations in the transcriptome than did exposure to the extracellular product GTX2/3. Mechanisms of effects were different between exposure of A. salina to A. minutum cells or GTX2/3. Exposure to A. minutum modulated formation of the ribonucleoprotein complex and metabolism of amino acids and lipids in A. salina. Exposure to GTX2/3 exposure inhibited expression of genes related to metabolism of chitin, which might result in disruption of molting process or disturbed sheath morphogenesis. Overall, effects on transcription observed in this study represent the first report based on application of next generation sequencing techniques to investigate the transcriptomic response of A. salina exposed to an environmentally realistic level of A. minutum or GTX2/3.

Diverse expression regulation of Hsp70 genes in scallops after exposure to toxic Alexandrium dinoflagellates

Heat shock proteins 70KD (Hsp70s) are highly conserved molecular chaperones with essential roles against biotic and abiotic stressors. Marine bivalves inhabit highly complex environments and could accumulate paralytic shellfish toxins (PSTs), the well-noted neurotoxins generated during harmful algal blooms. Here, we systematically analyzed Hsp70 genes (CfHsp70s) in Zhikong scallop (Chlamys farreri), an important aquaculture mollusk in China. Sixty-five CfHsp70s from eight sub-families were identified, and 47 of these genes showed expansion in the Hspa12 sub-family. After exposure to different PST-producing dinoflagellates, Alexandrium minutum and Alexandrium catenella, diverse CfHsp70s regulation presented in scallop hepatopancreas, mainly accumulating incoming PSTs, and kidneys, transforming PSTs into higher toxic analogs. All the up-regulated CfHsp70s were from CfHsp70B2, CfHspa12, and CfHspa5 sub-families. CfHsp70B2 sub-family was mainly induced in the hepatopancreas, and CfHspa12 sub-family was highly induced in the kidneys. CfHsp70s up-regulation under two dinoflagellates exposure was stronger in the kidneys (log₂FC: 19.5 and 18.6) than that in hepatopancreas (log₂FC: 4.3 and 6.1). Exposure to different Alexandrium species had varying effects, that in hepatopancreas, CfHsp70B2s were chronically induced only after A. catenella exposure, whereas in kidney, CfHspa12s were more acutely induced after exposure of A. minutum than A. caenella. Moreover, in Yesso scallops (Patinopecten yessoensis), only Hspa12s were up-regulated in hepatopancreas after A. catenella exposure, and all the Hsp70B2s were down-regulated. These organ-, toxin-, and species-dependent Hsp70 regulation suggested the functional diversity of duplicated Hsp70s in response to the stress by PST-producing algae. Our findings provide insights into the evolution and functional characteristics of Hsp70s in scallops.

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.

Molecular mechanism of glucose-6-phosphate utilization in the dinoflagellate Karenia mikimotoi

Phosphorus (P) is an essential nutrient for marine phytoplankton as for other living organisms, and the preferred form, dissolved inorganic phosphate (DIP), is often quickly depleted in the sunlit layer of the ocean. Phytoplankton have developed mechanisms to utilize organic forms of P (DOP). Hydrolysis of DOP to release DIP by alkaline phosphatase is believed to be the most common mechanism of DOP utilization. Little effort has been made, however, to understand other potential molecular mechanisms of utilizing different types of DOP. This study investigated the bioavailability of glucose-6-phosphate (G6P) and its underlying molecular mechanism in the dinoflagellate Karenia mikimotoi. Suppression Subtraction Hybridization (SSH) was used to identify genes up- and down-regulated during G6P utilization compared to DIP condition. The results showed that G6P supported the growth and yield of K. mikimotoi as efficiently as DIP. Neither DIP release nor AP activity was detected in the cultures grown in G6P medium, however, suggesting direct uptake of G6P. SSH analysis and RT-qPCR results showed evidence of metabolic modifications, particularly that mitochondrial ATP synthase f1gamma subunit and thioredoxin reductase were up-regulated while diphosphatase and pyrophosphatase were down-regulated in the G6P cultures. All the results indicate that K. mikimotoi has developed a mechanism other than alkaline phosphatase to utilize G6P.

Different transcriptomic responses of two marine copepods, Tigriopus japonicus and Pseudodiaptomus annandalei, to a low dose of mercury chloride (HgCl2)

Mercury (Hg) pollution is a ubiquitous and serious concern in marine environments, but the response mechanisms of marine animals to Hg pollution (i.e., toxicity/tolerance) are poorly understood. To compare the global responses of two marine copepods (Tigriopus japonicus and Pseudodiaptomus annandalei), we analyzed whole transcriptomes using RNA-seq technology in response to Hg treatment (a nominal 10μg/L HgCl₂ in seawater) for 5h. Hg was strikingly accumulated in both copepods under treatment. The Hg concentration in P. annandalei was higher under metal exposure by approximately 1.4-fold compared with treated T. japonicus. Among transcriptomic data, 101 genes in T. japonicus and 18 genes in P. annandalei were differentially regulated in response to Hg exposure. The up-regulated genes in T. japonicus were concerned with stress, growth, and development, while the down-regulated ones were mainly related to immune response. In P. annandalei, most of the differentially expressed genes were up-regulated, and all were involved in stress response. Our work indicated that Hg exhibits endocrine-disrupting potential at the transcriptomic level in marine copepods. Overall, our study demonstrates the species-specific molecular responses of these two copepods to Hg pollution.

Glutathione S-Transferase Regulation in Calanus finmarchicus Feeding on the Toxic Dinoflagellate Alexandrium fundyense

The effect of the dinoflagellate, Alexandrium fundyense, on relative expression of glutathione S-transferase (GST) transcripts was examined in the copepod Calanus finmarchicus. Adult females were fed for 5-days on one of three experimental diets: control (100% Rhodomonas spp.), low dose of A. fundyense (25% by volume, 75% Rhodomonas spp.), and high dose (100% A. fundyense). Relative expression of three GST genes was measured using RT-qPCR on days 0.5, 1, 2 and 5 in two independent experiments. Differential regulation was found for the Delta and the Sigma GSTs between 0.5 to 2 days, but not on day 5 in both experiments. The third GST, a microsomal, was not differentially expressed in either treatment or day. RT-qPCR results from the two experiments were similar, even though experimental females were collected from the Gulf of Maine on different dates and their reproductive output differed. In the second experiment, expression of 39 GSTs was determined on days 2 and 5 using RNA-Seq. Global gene expression analyses agreed with the RT-qPCR results. Furthermore, the RNA-Seq measurements indicated that only four GSTs were differentially expressed under the experimental conditions, and the response was small in amplitude. In summary, the A. fundyense diet led to a rapid and transient response in C. finmarchicus in three cytosolic GSTs, while a fourth GST (Omega I) was significantly up-regulated on day 5. Although there was some regulation of GSTs in response the toxic dinoflagellate, the tolerance to A. fundyense by C. finmarchicus is not dependent on the long-term up-regulation of specific GSTs.