Highly structured populations of copepods at risk to deep-sea mining: Integration of genomic data with demogenetic and biophysical modelling

Copepoda is the most abundant taxon in deep-sea hydrothermal vents, where hard substrate is available. Despite the increasing interest in seafloor massive sulphides exploitation, there have been no population genomic studies conducted on vent meiofauna, which are known to contribute over 50% to metazoan biodiversity at vents. To bridge this knowledge gap, restriction-site-associated DNA sequencing, specifically 2b-RADseq, was used to retrieve thousands of genome-wide single-nucleotide polymorphisms (SNPs) from abundant populations of the vent-obligate copepod Stygiopontius lauensis from the Lau Basin. SNPs were used to investigate population structure, demographic histories and genotype-environment associations at a basin scale. Genetic analyses also helped to evaluate the suitability of tailored larval dispersal models and the parameterization of life-history traits that better fit the population patterns observed in the genomic dataset for the target organism. Highly structured populations were observed on both spatial and temporal scales, with divergence of populations between the north, mid, and south of the basin estimated to have occurred after the creation of the major transform fault dividing the Australian and the Niuafo'ou tectonic plate (350 kya), with relatively recent secondary contact events (<20 kya). Larval dispersal models were able to predict the high levels of structure and the highly asymmetric northward low-level gene flow observed in the genomic data. These results differ from most studies conducted on megafauna in the region, elucidating the need to incorporate smaller size when considering site prospecting for deep-sea exploitation of seafloor massive sulphides, and the creation of area-based management tools to protect areas at risk of local extinction, should mining occur.

Characterization of the complete mitochondrial genome of Ergasilus anchoratus Markevich, 1946 (Ergasilidae) and phylogeny of Copepoda

The mitochondrial (mt) genome can provide data for phylogenetic analyses and evolutionary biology. Herein, we sequenced and annotated the complete mt genome of Ergasilus anchoratus. This mt genome was 13852 bp long and comprised 13 protein-coding genes (PCGs), 22 tRNAs and 2 rRNAs. All PCGs used the standard ATN start codons and complete TAA/TAG termination codons. A majority of tRNA genes exhibited standard cloverleaf secondary structures, with the exception of one tRNA that lacked the TψC arm (trnC), and three tRNAs that lacked the DHU arm (trnR, trnS1 and trnS2). Phylogenetic analyses conducted using Bayesian inference (BI) and maximum likelihood (ML) methods both supported Ergasilidae as a monophyletic family forming a sister group to Lernaea cyprinacea and Paracyclopina nana. It also supported the monophyly of orders Calanoida, Cyclopoida, and Siphonostomatoida; and the monophyly of families Harpacticidae, Ergasilidae, Diaptomidae, and Calanidae. The gene orders of E. anchoratus and Sinergasilus undulatus were identical, which represents the first instance of two identical gene orders in copepods. More mt genomes are needed to better understand the phylogenetic relationships within Copepoda in the future.

Establishing an astaxanthin-rich live feed strain of Pseudodiaptomus annandalei

This study aimed to establish an astaxanthin-rich strain of the calanoid copepod Pseudodiaptomus annandalei, through selective breeding based on RGB (red, green and blue) value, a parameter indicating color intensity. We evaluated the RGB value frequency distributions of the copepod populations, and selected individuals with the highest 10% and the lowest 10% RGB value over six generations. The RGB value, nauplii production, clutch interval and clutch number were assessed, and the genetic gain was calculated across generations (G0-G5). Two strains of copepods were selected and defined as dark body copepod strain (DBS) and light body copepod strain (LBS) at the end of experiment. Results revealed significantly lower RGB values (male: 121.5 ± 14.1; female: 108.8 ± 15) in the G5 DBS population compared to the G0 (male: 163.9 ± 13.1; female: 162.2 ± 14.6), with higher genetic gains of RGB values during G0 to G2. While DBS females exhibited longer clutch intervals in the G3 and G4, there was no significant difference in nauplii production between the two strains across all generations. Significantly higher astaxanthin content was found in the DBS copepods (0.04 μg/ ind.) compared to the LBS copepods (0.01 μg/ ind.) and the non-selective copepods (0.02 μg/ ind.) 20 months post selective breeding, validating the stability of the desired trait in the DBS strain. This study successfully established an astaxanthin-rich strain of P. annandalei, which provides implications for enhancing marine and brackish larviculture production.

Sex-Biased Transcription Expression of Vitellogenins Reveals Fusion Gene and MicroRNA Regulation in the Sea Louse Caligus rogercresseyi

The caligid ectoparasite, Caligus rogercresseyi, is one of the main concerns in the Chilean salmon industry. The molecular mechanisms displayed by the parasite during the reproductive process represent an opportunity for developing novel control strategies. Vitellogenin is a multifunctional protein recognized as a critical player in several crustaceans' biological processes, including reproduction, embryonic development, and immune response. This study aimed to characterize the C. rogercresseyi vitellogenins, including discovering novel transcripts and regulatory mechanisms associated with microRNAs. Herein, vitellogenin genes were identified by homology analysis using the reference sea louse genome, transcriptome database, and arthropods vitellogenin-protein database. The validation of expression transcripts was conducted by RNA nanopore sequencing technology. Moreover, fusion gene profiling, miRNA target analysis, and functional validation were performed using luciferase assay. Six putative vitellogenin genes were identified in the C. rogercresseyi genome with high homology with other copepods vitellogenins. Furthermore, miR-996 showed a putative role in regulating the Cr_Vitellogenin1 gene, which is highly expressed in females. Moreover, vitellogenin-fusion genes were identified in adult stages and highly regulated in males, demonstrating sex-related expression patterns. In females, the identified fusion genes merged with several non-vitellogenin genes involved in biological processes of ribosome assembly, BMP signaling pathway, and biosynthetic processes. This study reports the genome array of vitellogenins in C. rogercresseyi for the first time, revealing the putative role of fusion genes and miRNA regulation in sea lice biology.

Effects of ocean warming on the fatty acid and epigenetic profile of Acartia tonsa: A multigenerational approach

The effects of climate change are becoming more prevalent, and it is important to know how copepods, the most abundant class in zooplankton, will react to changing temperatures as they are the main food source for secondary consumers. They act as key transferers of nutrients from primary producers to organisms higher up the food chain. Little is known about the effects of temperature changes on copepods on the long term, i.e., over several generations. Especially the epigenetic domain seems to be understudied and the question remains whether the nutritional value of copepods will permanently change with rising water temperatures. In this research, the effects of temperature on the fatty acid and epigenetic profiles of the abundant planktonic copepod Acartia tonsa were investigated, since we expect to see a link between these two. Indeed, changing methylation patterns helped copepods to deal with higher temperatures, which is in line with the relative abundance of the most important fatty acids, e.g., DHA. However, this pattern was only observed when temperature increased slowly. A sudden increase in temperature showed the opposite effect; Acartia tonsa did not show deviant methylation patterns and the relative abundance of DHA and other important fatty acids dropped significantly after several generations. These results suggest that local fluctuations in temperature have a greater effect on Acartia tonsa than an elevation of the global mean.

Molecular quantification of parasitic sea louse larvae depends on species and life stage

Sea lice cause substantial economic and environmental harm to Norway's aquaculture industry and wild salmonid populations. Rapid, accurate quantification of lice larval densities in coastal waters remains the greatest bottleneck for providing empirical data on infestation risk within wild salmon habitats and aquaculture production regions. We evaluated the capability of droplet digital PCR (ddPCR) as an absolute quantification method for the planktonic stages of two parasitic louse species, Lepeophtheirus salmonis (Krøyer) and Caligus elongatus (von Nordman). Results demonstrated linear relationships between the DNA quantity measured and the number of spiked larvae for both species and life stages. However, L. salmonis contained a significantly greater number of DNA copies than C. elongatus individuals and for C. elongatus, nauplii displayed a significantly higher number of DNA copies than copepodids. Our results suggest that ddPCR can effectively enumerate louse larvae, but interpreting ddPCR results differ between the two louse species. Obtaining larval abundance estimates from marine plankton samples will depend on the nauplii to copepodid ratio for C. elongatus, but not for L. salmonis.

New species of Dermoergasilus Ho & Do, 1982 (Copepoda: Cyclopoida: Ergasilidae) parasitizing endemic cichlid Paretroplus polyactis (Bleeker) in Madagascar

Dermoergasilus madagascarensis n. sp. is described from the gills of Paretroplus polyactis, an endemic cichlid fish in Madagascar, using a combined morphological (light microscopy and SEM) and molecular approach (partial 18S rDNA, 28S rDNA, and COI sequences). The new species is characterized mainly by possessing: (i) roughly pentagonal cephalosome; (ii) antennal endopodal segments covered with slightly inflated membrane; (iii) maxillule bearing 2 equally long outer setae and a minute inner seta; (iv) interpodal sternites of swimming legs ornamented with 3–4 rows of spinules; (v) genital segment and first abdominal somite both barrel-shaped; and (vi) a caudal ramus projecting into a digitiform process with inconspicuous terminal seta and bearing 3 terminal setae. The obtained DNA sequences of Malagasy species represent the first molecular data for species of Dermoergasilus. The 28S rDNA phylogeny showed the affiliation of D. madagascarensis n. sp. to Ergasilidae and its sister relationship with cosmopolitan Ergasilus sieboldi von Nordmann, 1832. The first checklist for all species of Dermoergasilus is provided.

Diel metabolic patterns revealed by in situ transcriptome and proteome in a vertically migratory copepod

Zooplankton undergo a diel vertical migration (DVM) which exposes them to gradients of light, temperature, oxygen, and food availability on a predictable daily schedule. Disentangling the co-varying and potentially synergistic interactions on metabolic rates has proven difficult, despite the importance of this migration for the delivery of metabolic waste products to the distinctly different daytime (deep) and nighttime (surface) habitats. This study examines the transcriptomic and proteomic profiles of the circumglobal migratory copepod, Pleuromamma xiphias, over the diel cycle. The transcriptome showed that 96% of differentially expressed genes were upregulated during the middle of the day - the period often considered to be of lowest zooplankton activity. The changes in protein abundance were more spread out over time, peaking (42% of comparisons) in the early evening. Between 9:00 and 15:00, both the transcriptome and proteome datasets showed increased expression related to chitin synthesis and degradation. Additionally, at 09:00 and 22:00, there were increases in myosin and vitellogenin proteins, potentially linked to the stress of migration and/or reproductive investment. Based on protein abundances detected, there is an inferred switch in broad metabolic processes, shifting from electron transport system in the day to glycolysis and glycogen mobilization in the afternoon/evening. These observations provide evidence of the diel impact of DVM on transcriptomic and proteomic pathways that likely influence metabolic processes and subsequent excretion products, and clarify how this behaviour results in the direct rapid transport of waste metabolites from the surface to the deep ocean.

Comprehensive molecular characterisation of the complete mitogenome of Ergasilus tumidus and phylogenetic relationships of Copepoda inferred from mitogenomes

Although parasitic copepods of the genus Ergasilus von Nordmann, 1832 are globally distributed parasites of fish, their phylogenetic relationships with other Copepoda are not clear, and the characteristics of their mitochondrial genomes (mitogenomes) are not thoroughly understood. The objective of this study was to address these knowledge gaps by sequencing the complete mitogenome of Ergasilus tumidus Markevich, 1940. The complete mitogenome (GenBank Acc. No. OQ596537) was 14,431 bp long and it comprised 13 protein-coding genes (PCGs), 22 tRNAs, two tRNAs, and two control regions (CRs). Phylogenetic analyses, conducted using concatenated nucleotide and amino acid sequences of 13 protein-coding genes, produced two partially incongruent topologies. While the order Calanoida was consistently resolved as the sister lineage to the other three orders, topological instability was observed in the relationships of the orders Cyclopoida, Siphonostomatoida and Harpacticoida. Siphonostomatoida clustered with Cyclopoida in the nucleotide-based phylogeny, but with Harpacticoida in the amino acid-based phylogeny. The latter topology conforms to the widely accepted relationships, but we speculate that the former topology is more likely to be the correct one. Our study provides a complete mitogenome sequence of E. tumidus, which helps us better understand the molecular evolution of the genus Ergasilus. Additionally, we suggest a different perspective on the controversial phylogenetic relationships among Siphonostomatoida, Cyclopoida and Harpacticoida, diverging from previously accepted views.

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.

Analysis of the transcriptional pathways associated with the induction of quiescent embryonic arrest in the calanoid copepod Acartia tonsa

The copepod species Acartia tonsa (Dana)(Crustacea) have the unique ability to induce quiescent embryonic dormancy if adverse environmental conditions occur; a characteristic shared by 41 other species belonging to the superfamily Centropagoida in the Calanoida order. However, the transcriptional changes characterizing this process are not known. Here, we compare the transcriptome of embryos in arrested quiescence with the normal development to identify pathways and differentially regulated transcripts involved in quiescent embryogenesis. Quiescence was induced by incubating eggs at 4 °C with anoxia for 26 h(hr), while eggs undergoing normal immediate development were incubated at 16.9 °C in normoxia for 7 h (where gastrulation occurs) or 14 h (where organogenesis occurs) before collecting for RNA extraction and analysis by RNA-sequencing. Results indicate that the expression profile of the quiescent embryo is not as different from the normal embryonic gastrulation as initially expected: None of the mapped transcripts is uniquely expressed in quiescence. Moreover, in quiescence a large proportion of the annotated transcripts display expression values halfway in-between the normal, immediate developmental stages of gastrulation and organogenesis. In depth comparison between the organogenesis stage and quiescent samples, reveal a high degree of divergence, confirming that a developmental arrest has been induced through quiescence. Specifically: Stress response transcripts are prominent in the quiescent phase with a transcript like the mammalian autophagy gene Sequestosome-1/p62 (SQSTM) being upregulated. The present analysis provides a better understanding of the molecular mechanisms characterizing the quiescent embryonic state of A. tonsa.

Disproportionate role of nuclear-encoded proteins in organismal and mitochondrial thermal performance in a copepod

Determining the mechanisms by which organisms evolve thermal tolerance is crucial to predicting how populations may respond to changes in local temperature regimes. Although evidence of relationships between mitochondrial background and thermal adaptation have been found, the presence of both nuclear-encoded and mitochondrial DNA (mtDNA)-encoded proteins warrants experiments aimed at parsing out the relative role of each genome in thermal adaptation. We investigated the relative role of mtDNA-encoded products in thermal tolerance between two divergent populations of Tigriopus californicus using first-generation (F1) hybrids that vary in maternally inherited mtDNA but are heterozygous for population-specific alleles across nuclear loci. We tested two measures of thermal tolerance, (1) survivorship to acute thermal stress and (2) thermal stability of mitochondrial performance in Complex I-fueled ATP synthesis, both across a range of increasing temperatures. We found that the southern population (San Diego, CA, USA) outperformed the northern population (Strawberry Hill, OR, USA) in survivorship, and that both reciprocal F1 hybrid crosses had intermediate survival. Mitochondria from the San Diego population displayed greater stability in ATP synthesis with increasing temperatures compared with those from Strawberry Hill. Interestingly, hybrids from both cross directions had synthesis profiles that were very similar to that of Strawberry Hill. Taken together, these results suggest that the relative role of the mtDNA in these phenotypes is negligible compared with that of elements encoded by nuclear DNA in this system.

Highly structured populations of deep-sea copepods associated with hydrothermal vents across the Southwest Pacific, despite contrasting life history traits

Hydrothermal vents are extreme environments, where abundant communities of copepods with contrasting life history traits co-exist along hydrothermal gradients. Here, we discuss how these traits may contribute to the observed differences in molecular diversity and population genetic structure. Samples were collected from vent locations across the globe including active ridges and back-arc basins and compared to existing deep-sea hydrothermal vent and shallow water data, covering a total of 22 vents and 3 non-vent sites. A total of 806 sequences of mtDNA from the Cox1 gene were used to reconstruct the phylogeny, haplotypic relationship and demography within vent endemic copepods (Dirivultidae, Stygiopontius spp.) and non-vent-endemic copepods (Ameiridae, Miraciidae and Laophontidae). A species complex within Stygiopontius lauensis was studied across five pacific back-arc basins at eight hydrothermal vent fields, with cryptic species being restricted to the basins they were sampled from. Copepod populations from the Lau, North Fiji and Woodlark basins are undergoing demographic expansion, possibly linked to an increase in hydrothermal activity in the last 10 kya. Highly structured populations of Amphiascus aff. varians 2 were also observed from the Lau to the Woodlark basins with populations also undergoing expansion. Less abundant harpacticoids exhibit little to no population structure and stable populations. This study suggests that similarities in genetic structure and demography may arise in vent-associated copepods despite having different life history traits. As structured meta-populations may be at risk of local extinction should major anthropogenic impacts, such as deep-sea mining, occur, we highlight the importance of incorporating a trait-based approach to investigate patterns of genetic connectivity and demography, particularly regarding area-based management tools and environmental management plans.

Transcriptomic analysis reveals responses to a polluted sediment in the Mediterranean copepod Acartia clausi

Marine sediments are regarded as sinks for several classes of contaminants. Characterization and effects of sediments on marine biota now require a multidisciplinary approach, which includes chemical and ecotoxicological analyses and molecular biomarkers. Here, a gene expression study was performed to measure the response of adult females of the Mediterranean copepod Acartia clausi to elutriates of polluted sediments (containing high concentrations of polycyclic aromatic hydrocarbons, PAHs, and heavy metals) from an industrial area in the Southern Tyrrhenian Sea (Bagnoli-Coroglio). Functional annotation of the A. clausi transcriptome generated as reference here, showed a good quality of the assembly and great homology with other copepod and crustacean sequences in public databases. This is one of the few available transcriptomic resources for this widespread copepod species of great ecological relevance in temperate coastal areas. Differential expression analysis between females exposed to the elutriate and those in control seawater identified 1000 differentially expressed genes, of which 743 up- and 257 down-regulated. Within the up-regulated genes, the most represented functions were related to proteolysis (lysosomal protease, peptidase, cathepsin), response to stress and detoxification (heat-shock protein, superoxide dismutase, glutathione-S-transferase, cytochrome P450), and cytoskeleton structure (α- and β-tubulin). Down-regulated genes were mostly involved with ribosome structure (ribosomal proteins) and DNA binding (histone proteins, transcription factors). Overall, these results suggest that processes such as transcription, translation, protein degradation, metabolism of biomolecules, reproduction, and xenobiotic detoxification were altered in the copepod in response to polluted elutriates. In conclusion, our results contribute to gaining information on the transcriptomic responses of copepods to polluted sediments. They will also prompt the selection of genes of interest to be used as biomarkers of exposure to PAHs and heavy metals in molecular toxicology studies on copepods, and in general, in comparative functional genomic studies on marine zooplankton.

Genetic incompatibilities in reciprocal hybrids between populations of Tigriopus californicus with low to moderate mitochondrial sequence divergence

All mitochondrial-encoded proteins and RNAs function through interactions with nuclear-encoded proteins, which are critical for mitochondrial performance and eukaryotic fitness. Coevolution maintains inter-genomic (i.e., mitonuclear) compatibility within a taxon, but hybridization can disrupt coevolved interactions, resulting in hybrid breakdown. Thus, mitonuclear incompatibilities may be important mechanisms underlying reproductive isolation and, potentially, speciation. Here we utilize Pool-seq to assess the effects of mitochondrial genotype on nuclear allele frequencies in fast- and slow-developing reciprocal inter-population F2 hybrids between relatively low-divergence populations of the intertidal copepod Tigriopus californicus. We show that mitonuclear interactions lead to elevated frequencies of coevolved (i.e., maternal) nuclear alleles on two chromosomes in crosses between populations with 1.5% or 9.6% fixed differences in mitochondrial DNA nucleotide sequence. However, we also find evidence of excess mismatched (i.e., noncoevolved) alleles on three or four chromosomes per cross, respectively, and of allele frequency differences consistent with effects involving only nuclear loci (i.e., unaffected by mitochondrial genotype). Thus, our results for low-divergence crosses suggest an underlying role for mitonuclear interactions in variation in hybrid developmental rate, but despite substantial effects of mitonuclear coevolution on individual chromosomes, no clear bias favoring coevolved interactions overall.

Silencing of the Vasa gene by RNA Interference Affects Embryonic Development and Reproductive Output in the Sea Louse Caligus rogercresseyi

The sea louse Caligus rogercresseyi is a major ectoparasitic copepod that causes significant economic losses in the salmon farming industry. Despite recent advancements, the mechanisms underlying germline and embryo development in this species remain poorly understood. The Vasa gene encodes a highly conserved DEAD box helicase that is required for germ cell formation and function in many species. In this study, the Vasa gene was characterized in C. rogercresseyi, and its expression and function were analyzed. Phylogenetic analysis showed that the Cr-Vasa gene product formed clusters in clades with Vasa proteins from closely related species of crustaceans. Cr-Vasa gene expression patterns were assessed by qPCR, and the results showed a significantly higher relative expression level in adult females compared to copepodid, chalimus, and adult male stages. Tissue-specific localization of Cr-Vasa mRNA in C. rogercresseyi was determined using chromogenic in situ hybridization, and strong positive signal was observed in male testes, but also in the intestine and cuticle, while in females, it was observed in the ovaries, oocytes, cuticle, intestine, and egg strings. RNAi-mediated gene silencing of Cr-Vasa impacted embryonic development and reproductive output in adult female lice. Females from the dsVasa-treated group displayed unusual phenotypes, including shorter egg strings with numerous extra-embryonic inclusions, irregularly shaped abnormal embryos, and aborted egg strings. This study provides insights into the role of the Vasa gene in C. rogercresseyi embryonic development and reproductive output, which may have implications for the control of this parasitic copepod in the salmon farming industry.

Joint effects of temperature and copper exposure on developmental and gene-expression responses of the marine copepod Tigriopus japonicus

There is growing contamination of copper (Cu) in the marine environment, particularly after the ban of organotin compounds and the increase of the use of Cu-based antifouling paints. Although there are increasing research interests in temperature-dependent chemical toxicity to aquatic organisms, most existing studies focused on acute impacts of chemicals at high concentrations. This study aimed to investigate the interacting effect of temperature and copper exposure at environmentally relevant concentrations on survival and development in the marine copepod Tigriopus japonicus with a partial life-cycle toxicity test. Expressions of five stress response genes in the copepod, namely two glutathione S-transferases (GST-S and GST-O), two heat shock proteins (HSP70 and HSP90), and glutathione reductase (GR) were also investigated. The copepod's survival was significantly impaired at 15 °C after development to adult stage, while its developmental time reduced significantly with increasing temperature. Copper at the two environmentally relevant test concentrations had no significant impacts on these apical endpoints whereas the interaction between Cu and temperature was more significant in modulating gene expressions. GST-S, GST-O and HSP90 genes in copepods exposed to 100 µg Cu L^-1 were significantly upregulated at 20 °C. At 32 °C, most genes were either insignificantly expressed or down-regulated, compared to the control, likely suggesting that thermal stress inhibited the copepod's antioxidative defense system. Overall, the results revealed that the joint Cu and thermal stresses have significantly elicited antioxidative system in the copepods. It clearly demonstrated the need for more fundamental studies about potential impacts of different environmental factors such as temperature on chemical toxicity under realistic scenario of marine pollution.

Mining Lepeophtheirus salmonis RNA-Seq data for qPCR reference genes and their application in Caligus elongatus

Lepeophtheirus salmonis and Caligus elongatus are two parasitic copepod species posing a significant threat to salmonid aquaculture. Consequently, several gene expression studies are executed each year to gain new knowledge and treatment strategies. Though, to enable accurate gene expression measurements by quantitative real time PCR, stable reference genes are needed. Previous studies have mainly focused on a few genes selected based on their function as housekeeping genes, as these are often stably expressed in various cells and tissues. In the present study, however, RNA-sequencing data from 127 L. salmonis samples from different life stages and diverse environmental conditions were used to identify new candidate reference genes displaying low variation. From this, six genes were selected, and the stability validated by qPCR on samples from different life stages. Since neither a genome nor comprehensive RNA sequencing data are available for C. elongatus, homologous genes to those identified for L. salmonis were identified within a C. elongatus transcriptome assembly and validated by qPCR in different life stages. Overall, the genes eukaryotic translation initiation factor 1A (EIF1A) and serine/threonine-protein phosphatase 1 (PP1) displayed the highest stability in L. salmonis, while the combination of PP1 and ribosomal protein S13 (RPS13) was found to have the highest stability in C. elongatus. These genes are well-suited reference genes for qPCR applications which allow for accurate normalization of target genes.

Embryonic and post-embryonic development in the parasitic copepod Ive ptychoderae (Copepoda: Iviidae): Insights into its phylogenetic position

Parasitic copepods are frequently discovered in many marine animals, and they exhibit great species diversity with remarkable morphological adaptations to their parasitic lifestyle. Similar to their free-living relatives, parasitic copepods usually develop through complex life cycle, but they eventually transform into a modified adult form with reduced appendages. Although the life cycle and distinct larval stages have been described in a few species of parasitic copepods, particularly those infecting commercially valuable marine animals (such as fishes, oysters, and lobsters), very little is known about the developmental process of the species that transformed into extremely simplified adult body plan. This paucity also causes some difficulties when investigating the taxonomy and phylogeny of this kind of parasitic copepods. Here we describe the embryonic development and a series of sequential larval stages of a parasitic copepod, Ive ptychoderae, which is a vermiform endoparasite living inside the hemichordate acorn worms. We devised laboratory regimes that enable us raising large quantity of embryos and free living larvae, and obtaining post-infested I. ptychoderae samples from the host tissues. Using defined morphological features, the embryonic development of I. ptychoderae can be categorized into eight stages (1-, 2-, 4-, 8-, 16- cell stages, blastula, gastrula, and limb bud stages) and the post-embryonic development comprises six larval stages (2 naupliar and 4 copepodid stages). Based on the comparisons of morphological characters in the nauplius stage, our results provide evidence to support that the Ive-group is more closely related to the Cyclopoida, which represents one of the two major clades that contain many highly transformed parasitic copepods. Thus, our results help to resolve the problematic phylogenetic position of the Ive-group in previous study based on analysis using 18S rDNA sequences. Combining with more molecular data, future comparative analyses on the morphological features of copepodid stages will further refine our understanding of the phylogenetic relationships of parasitic copepods.

Morphological and mitochondrial approaches of Hatschekia sargi (Copepoda: Hatschekiidae) as a parasite of Epinephelus chlorostigma

Little information is available until now about the copepods infecting different fish species. Therefore, this study aimed to provide light on siphonostomatoids infecting Epinephelus chlorostigma. Twenty fish specimens were taken from the Red Sea coast (Jeddah, Saudi Arabia), and ectoparasitic copepods were investigated. Light microscopy and molecular tools were used to examine copepods isolated from fish. Parasitological indexes were calculated and showed that 60% of the examined fish were infected with a mean intensity of 12 parasite/fish. Morphological examination revealed that this copepod species is characterized by all unique features of the genus Hatschekia with special reference to Hatschekia sargi. The taxonomic position of the recovered species in the Hatschekiidae family within Siphonostomatoida was confirmed using phylogenetic analysis based on partial mitochondrial cytochrome c oxidase subunit I (mt COI) gene sequences. The mt COI gene query revealed that the recovered Hatschekia species is closely related to Hatschekia maculatus (gb| JQ664005.1). This study discovers a new host for Hatschekia species isolated from Saudi Arabia and conducts the first genomic investigation of the mt COI gene.

Transcriptional and toxic responses to saxitoxin exposure in the marine copepod Tigriopus japonicus

Saxitoxin (STX) is a highly toxic marine neurotoxin produced by phytoplankton and a growing threat to ecosystems worldwide due to the spread of toxic algae. Although STX is an established sodium channel blocker, the overall profile of transcriptional levels in STX-exposed organisms has yet to be described. Here, we describe a toxicity assay and transcriptome analysis of the copepod Tigriopus japonicus exposed to STX. The half-maximal lethal concentration of STX was 12.35 μM, and a rapid mortality slope was evident at concentrations between 12 and 13 μM. STX induced changes in swimming behavior among the copepods after 10 min of exposure. In transcriptome analysis, gene ontology revealed that the genes involved in nervous system and gene expression were highly enriched. In addition, the congenital neurological disorder and nuclear factor erythroid 2-related factor 2-mediated oxidative stress pathways were identified to be the most significant in network analysis and toxicity pathway analysis, respectively. This study provides valuable information about the effects of STX and related transcriptional responses in T. japonicus.

A ferritin gene in the marine copepod Acartia tonsa as a highly sensitive biomonitor for nano-contamination

Toxicology is not only for eco-risk assessments, but also for the real-time environmental monitoring based on the quick response of specific biomarkers. Ferritin gene (ftn) is a potential biomarker involving in crucial protective responses in biota. However, little information is available concerning the ftn in marine copepod Acartia tonsa (A. tonsa), a model organism widely applied in toxicology assessments. Our study for the first time identified and characterized the ftn in A. tonsa, along with its time-dependent transcriptional response to the reproductive toxicity of two newly emerged nanomaterials. The full-length cDNA of ftn contains a 114-bp 5'-untranslated region (UTR), a 236-bp 3'-untranslated region, and a 510-bp open reading frame which encodes an 18.51 kDa polypeptide composed of 169 amino acids. The ftn sequence has an iron binding signature and a potential phosphorylation site, which is closely-related to the ftn of Calanus sinicus and Pseudodiaptomus annandalei genes at the phylogenetical level. The ftn showed a quick and highly sensitive response to nanomaterial exposures, even at no observed effect concentrations. In detail, after exposure to nickel nanomaterials (up to 17.0 mg/L), the ftn was significantly upregulated immediately at 0.5 h and peaked at 9.5-fold in adults within 48 h, along with a significant reduction of egg hatching rate. When exposed to CdSe/ZnS quantum dots (up to 135 mg/L), no significant change in egg productions or hatching rates was observed, while the expression of ftn still significantly increased to over 3.0-fold in the initial 48 h. After that, the upregulation of ftn induced by CdSe/ZnS quantum dots or nickel nanoparticles both gradually returned back within 96 h. These findings demonstrate the highly sensitive response of this new cloned ftn to nanomaterial exposures, and highlight the suitability of ftn in A. tonsa as a promising biomonitor for nano-contamination in marine environments.

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.

RNA-seq transcriptome analysis and identification of the theromacin antimicrobial peptide of the copepod Apocyclops royi

Copepods, including Apocyclops royi, are small aquatic crustaceans and one of the important foods for fish and shellfish larvae. However, studies of the host-pathogen interactions and understanding of infectious disease in copepods are still very limited, yet they are likely to be a significant factor in the sustainable development of copepod aquaculture. In the present study, we performed de novo RNA sequence analysis of A. royi-TH (a Thai isolate of A. royi), which yielded 4.80 Gb bases of clean data and a total of 29,786 unigenes. Annotation was then performed by comparison against seven functional databases, yielding 17,617 (NR: 59.15%), 2,969 (NT: 9.97%), 15,023 (SwissProt: 50.44%), 14,543 (KOG: 48.82%), 15,077 (KEGG: 50.62%), 6,763(GO: 22.71%), and 15,841 (InterPro: 53.18%) unigenes. In comparison to the components of the shrimp Toll pathway, LGBP, Spätzle, Toll receptors, MyD88, Pelle, TRAF6, Dorsal, and Cactus homologs were successfully identified in A. royi-TH. Additionally, a novel antimicrobial peptide (Theromacin-like) was characterized in A. royi (ArTM-like). The ArTM-like ORF was 279 bp and predicted to encode for 92 amino acid residues, with a mature peptide of 75 amino acids and a molecular mass of 8.56 kDa. The genomic organization of the ArTM-like gene consisted of three exons and two introns. Expression analysis indicated that ArTM-like mRNA was abundantly expressed in copepodid and adult stages as an immune responsive gene after infection with the pathogenic Vibrio parahaemolyticus-(AHPND)-causing strain. Altogether, the knowledge obtained in this study will provide a basis for future functional studies of the molecular mechanisms in copepod immunity that may eventually be applied for disease prevention in copepod aquaculture.

Sea lice (Lepeophtherius salmonis) detection and quantification around aquaculture installations using environmental DNA

The naturally occurring ectoparasite salmon lice (Lepeophtherirus salmonis) poses a great challenge for the salmon farming industry, as well as for wild salmonids in the Northern hemisphere. To better control the infestation pressure and protect the production, there is a need to provide fish farmers with sensitive and efficient tools for rapid early detection and monitoring of the parasitic load. This can be achieved by targeting L. salmonis DNA in environmental samples. Here, we developed and tested a new L. salmonis specific DNA-based assay (qPCR assay) for detection and quantification from seawater samples using an analytical pipeline compatible with the Environmental Sample Processor (ESP) for autonomous water sample analysis of gene targets. Specificity of the L. salmonis qPCR assay was demonstrated through in-silico DNA analyses covering sequences of different L. salmonis isolates. Seawater was spiked with known numbers of nauplii and copepodite free-swimming (planktonic) stages of L. salmonis to investigate the relationship with the number of marker gene copies (MGC). Finally, field samples collected at different times of the year in the vicinity of a salmon production farm in Western Norway were analyzed for L. salmonis detection and quantification. The assay specificity was high and a high correlation between MGC and planktonic stages of L. salmonis was established in the laboratory conditions. In the field, L. salmonis DNA was consequently detected, but with MGC number below that expected for one copepodite or nauplii. We concluded that only L. salmonis tissue or eDNA residues were detected. This novel study opens for a fully automatized L. salmonis DNA quantification using ESP robotic to monitor the parasitic load, but challenges remain to exactly transfer information about eDNA quantities to decisions by the farmers and possible interventions.

Identification and characterization of homeobox gene clusters in harpacticoid and calanoid copepods

In this study, we have identified the entire complement of typical homeobox (Hox) genes (Lab, Pb, Dfd, Scr, Antp, Ubx, Abd-A, and Abd-B) in harpacticoid and calanoid copepods and compared them with the cyclopoid copepod Paracyclopina nana. The harpacticoid copepods Tigriopus japonicus and Tigriopus kingsejongensis have seven Hox genes (Lab, Dfd, Scr, Antp, Ubx, Abd-A, and Abd-B) and the Pb and Ftz genes are also present in the cyclopoid copepod P. nana. In the Hox gene cluster of the calanoid copepod Eurytemora affinis, all the Hox genes were present linearly in the genome but the Antp gene was duplicated. Of the three representative copepods, the P. nana Hox gene cluster was the most compact due to its small genome size. The Hox gene expression profile patterns in the three representative copepods were stage-specific. The Lab, Dfd, Scr, Pb, Ftz, and Hox3 genes showed a high expression in early developmental stages but Antp, Ubx, Abd-A, and Abd-B genes were mostly expressed in later developmental stages, implying that these Hox genes may be closely associated with the development of segment identity during early development.

An update of the salmon louse (Lepeophtheirus salmonis) reference genome assembly

Salmon lice have plagued the salmon farming industry and have negatively impacted salmon populations in the wild. In response, researchers have generated high density genetic maps, genome assemblies, transcriptomes, and whole-genome resequencing data to better understand this parasite. In this study, we used long-read sequencing technology to update the previous genome assemblies of Atlantic Ocean salmon lice with a more contiguous assembly and a more comprehensive gene catalog of Pacific Ocean salmon lice. We were also able to further characterize genomic features previously identified from other studies by using published resequenced genomes of 25 Atlantic and 15 Pacific salmon lice. One example was further characterizing the ZW sex chromosomes. For both the Atlantic and Pacific Ocean salmon lice subspecies, we found that the female W-chromosome is only a small fraction of the Z-chromosome and that the vast majority of the W and Z-chromosome do not contain conserved regions (i.e. pseudoautosomal regions). However, conserved orthologous protein sequences can still be identified between the W- and Z-chromosomes.

Frequent Paternal Mitochondrial Inheritance and Rapid Haplotype Frequency Shifts in Copepod Hybrids

Mitochondria are assumed to be maternally inherited in most animal species, and this foundational concept has fostered advances in phylogenetics, conservation, and population genetics. Like other animals, mitochondria were thought to be solely maternally inherited in the marine copepod Tigriopus californicus, which has served as a useful model for studying mitonuclear interactions, hybrid breakdown, and environmental tolerance. However, we present PCR, Sanger sequencing, and Illumina Nextera sequencing evidence that extensive paternal mitochondrial DNA (mtDNA) transmission is occurring in inter-population hybrids of T. californicus. PCR on four types of crosses between three populations (total sample size of 376 F1 individuals) with 20% genome-wide mitochondrial divergence showed 2% to 59% of F1 hybrids with both paternal and maternal mtDNA, where low and high paternal leakage values were found in different cross directions of the same population pairs. Sequencing methods further verified nucleotide similarities between F1 mtDNA and paternal mtDNA sequences. Interestingly, the paternal mtDNA in F1s from some crosses inherited haplotypes that were uncommon in the paternal population. Compared to some previous research on paternal leakage, we employed more rigorous methods to rule out contamination and false detection of paternal mtDNA due to non-functional nuclear mitochondrial DNA fragments. Our results raise the potential that other animal systems thought to only inherit maternal mitochondria may also have paternal leakage, which would then affect the interpretation of past and future population genetics or phylogenetic studies that rely on mitochondria as uniparental markers.

The salmon louse genome may be much larger than sequencing suggests

The genome size of organisms impacts their evolution and biology and is often assumed to be characteristic of a species. Here we present the first published estimates of genome size of the ecologically and economically important ectoparasite, Lepeophtheirus salmonis (Copepoda, Caligidae). Four independent L. salmonis genome assemblies of the North Atlantic subspecies Lepeophtheirus salmonis salmonis, including two chromosome level assemblies, yield assemblies ranging from 665 to 790 Mbps. These genome assemblies are congruent in their findings, and appear very complete with Benchmarking Universal Single-Copy Orthologs analyses finding > 92% of expected genes and transcriptome datasets routinely mapping > 90% of reads. However, two cytometric techniques, flow cytometry and Feulgen image analysis densitometry, yield measurements of 1.3-1.6 Gb in the haploid genome. Interestingly, earlier cytometric measurements reported genome sizes of 939 and 567 Mbps in L. salmonis salmonis samples from Bay of Fundy and Norway, respectively. Available data thus suggest that the genome sizes of salmon lice are variable. Current understanding of eukaryotic genome dynamics suggests that the most likely explanation for such variability involves repetitive DNA, which for L. salmonis makes up ≈ 60% of the genome assemblies.

First transcriptome of the copepod Gladioferens pectinatus subjected to chronic contaminant exposures

Contaminants are often at low concentrations in ecosystems and their effects on exposed organisms can occur over long periods of time and across multiple generations. Alterations to subcellular mechanistic pathways in response to exposure to contaminants can provide insights into mechanisms of toxicity that methods measuring higher levels of biological may miss. Analysis of the whole transcriptome can identify novel mechanisms of action leading to impacts in exposed biota. The aim of this study was to characterise how exposures to copper, benzophenone and diclofenac across multiple generations altered molecular expression pathways in the marine copepod Gladioferens pectinatus. Results of the study demonstrated differential gene expression was observed in cultures exposure to diclofenac (569), copper (449) and benzophenone (59). Pathways linked to stress, growth, cellular and metabolic processes were altered by exposure to all three contaminants with genes associated with oxidative stress and xenobiotic regulation also impacted. Protein kinase functioning, cytochrome P450, transcription, skeletal muscle contraction/relaxation, mitochondrial phosphate translocator, protein synthesis and mitochondrial methylation were all differentially expressed with all three chemicals. The results of the study also suggested that using dimethyl sulfoxide as a dispersant influenced the transcriptome and future research may want to investigate it's use in molecular studies. Data generated in this study provides a first look at transcriptomic response of G. pectinatus exposed to contaminants across multiple generations, future research is needed to validate the identified biomarkers and link these results to apical responses such as population growth to demonstrate the predictive capacity of molecular tools.

Production of Copepod Luciferases via Baculovirus Expression System

Secreted copepod luciferases (CopLucs) represent highly homologous enzymes which catalyze the oxidation of a low molecular weight substrate, coelenterazine, with the emission of blue light (λ_max = 485-488 nm), that is called bioluminescence (BL). The well-studied Gaussia (GLuc) and Metridia (MLuc) luciferases originally cloned from the marine copepods Gaussia princeps and Metridia longa belong to the group of the smallest natural luciferases. Their minimal molecular weight, high luminescent activity, cofactor-independent BL, and the ability to be secreted due to the own signal peptide open up the horizons for genetic engineering of CopLuc-based sensitive biosensors for in vivo imaging and in vitro analytical applications. The "standard" soluble bacterial expression of the recombinant CopLucs and luciferase-based hybrid proteins is hampered by the presence of high amounts of intramolecular disulfide bonds (up to 5 per molecule). Here, we describe the universal protocol for highly effective secreted expression of disulfide-rich CopLucs using their own signal peptide in insect cells and their purification from serum-free culture medium. The suggested protocol allows obtaining high-purity CopLucs folded in their native form with the yield of up to 5 mg per liter.

Evidence for hybrid breakdown in production of red carotenoids in the marine invertebrate Tigriopus californicus

The marine copepod, Tigriopus californicus, produces the red carotenoid pigment astaxanthin from yellow dietary precursors. This ‘bioconversion’ of yellow carotenoids to red is hypothesized to be linked to individual condition, possibly through shared metabolic pathways with mitochondrial oxidative phosphorylation. Experimental inter-population crosses of lab-reared T. californicus typically produces low-fitness hybrids is due in large part to the disruption of coadapted sets nuclear and mitochondrial genes within the parental populations. These hybrid incompatibilities can increase variability in life history traits and energy production among hybrid lines. Here, we tested if production of astaxanthin was compromised in hybrid copepods and if it was linked to mitochondrial metabolism and offspring development. We observed no clear mitonuclear dysfunction in hybrids fed a limited, carotenoid-deficient diet of nutritional yeast. However, when yellow carotenoids were restored to their diet, hybrid lines produced less astaxanthin than parental lines. We observed that lines fed a yeast diet produced less ATP and had slower offspring development compared to lines fed a more complete diet of algae, suggesting the yeast-only diet may have obscured effects of mitonuclear dysfunction. Astaxanthin production was not significantly associated with development among lines fed a yeast diet but was negatively related to development in early generation hybrids fed an algal diet. In lines fed yeast, astaxanthin was negatively related to ATP synthesis, but in lines fed algae, the relationship was reversed. Although the effects of the yeast diet may have obscured evidence of hybrid dysfunction, these results suggest that astaxanthin bioconversion may still be related to mitochondrial performance and reproductive success.

Transcriptome and morphological analysis in Caligus rogercresseyi uncover the effects of Atlantic salmon vaccination with IPath®

It is known that iron transporter proteins and their regulation can modulate the fish's immune system, suggesting these proteins as a potential candidate for fish vaccines. Previous studies have evidenced the effects of Atlantic salmon immunized with the chimeric iron-related protein named IPath® against bacterial and ectoparasitic infections. The present study aimed to explore the transcriptome modulation and the morphology of the sea louse Caligus rogercresseyi in response to Atlantic salmon injected with IPath®. Herein, Atlantic salmon were injected with IPath® and challenged to sea lice in controlled laboratory conditions. Then, female adults were collected after 25 days post-infection for molecular and morphological evaluation. Transcriptome analysis conducted in lice collected from immunized fish revealed high modulation of transcripts compared with the control groups. Notably, the low number of up/downregulated transcripts was mainly found in lice exposed to the IPath® fish group. Among the top-25 differentially expressed genes, Vitellogenin, Cytochrome oxidases, and proteases genes were strongly downregulated, suggesting that IPath® can alter lipid transport, hydrogen ion transmembrane transport, and proteolysis. The morphological analysis in lice collected from IPath® fish revealed abnormal embryogenesis and inflammatory processes of the genital segment. Furthermore, head kidney, spleen, and skin were also analyzed in immunized fish to evaluate the transcription expression of immune and iron homeostasis-related genes. The results showed downregulation of TLR22, MCHII, IL-1β, ALAs, HO, BLVr, GSHPx, and Ferritin genes in head kidney and skin tissues; meanwhile, those genes did not show significant differences in spleen tissue. Overall, our findings suggest that IPath® can be used to enhance the fish immune response, showing a promissory commercial application against lice infections.

Genomic analysis of the carboxylesterase family in the salmon louse (Lepeophtheirus salmonis)

The pyrethroid deltamethrin and the macrocyclic lactone emamectin benzoate (EMB) are used to treat infestations of farmed salmon by parasitic salmon lice, Lepeophtheirus salmonis. While the efficacy of both compounds against Atlantic populations of the parasite has decreased as a result of the evolution of resistance, the molecular mechanisms of drug resistance in L. salmonis are currently not fully understood. The functionally diverse carboxylesterases (CaE) family includes members involved in pesticide resistance phenotypes of terrestrial arthropods. The present study had the objective to characterize the CaE family in L. salmonis and assess its role in drug resistance. L. salmonis CaE homologues were identified by homology searches in the parasite's transcriptome and genome. The transcript expression of CaEs predicted to be catalytically competent was studied using quantitative reverse-transcription PCR in drug susceptible and multi-resistant L. salmonis. The above strategy led to the identification of 21 CaEs genes/pseudogenes. Phylogenetic analyses assigned 13 CaEs to clades involved in neurodevelopmental signaling and cell adhesion, while three sequences were predicted to encode secreted enzymes. Ten CaEs were identified as being potentially catalytically competent. Transcript expression of acetylcholinesterase (ace1b) was significantly increased in multi-resistant lice compared to drug-susceptible L. salmonis, with transcript abundance further increased in preadult-II females following EMB exposure. In summary, results from the present study demonstrate that L. salmonis possesses fewer CaE gene family members than most arthropods characterized so far. Drug resistance in L. salmonis was associated with overexpression of ace1b.

iTRAQ-based proteomic profiling, pathway analyses, and apoptotic mechanism in the Antarctic copepod Tigriopus kingsejongensis in response to ultraviolet B radiation

iTRAQ proteomic profiling was conducted to examine the proteomic responses of the Antarctic copepod Tigriopus kingsejongensis under ultraviolet B (UVB) exposure. Of the 5507 proteins identified, 3479 proteins were annotated and classified into 25 groups using clusters of orthologous genes analysis. After exposing the T. kingsejongensis to 12 kJ/m² UVB radiation, 77 biological processes were modulated over different time periods (0, 6, 12, 24, and 48 h) compared with the control. A Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis showed that UVB exposure in T. kingsejongensis downregulated ribosome and glyoxylate and dicarboxylate metabolism at all time points. Furthermore, antioxidant and chaperone proteins were highly downregulated in response to UVB exposure, causing protein damage and activating apoptotic processes in the 48 h UVB exposure group. These proteomic changes show the mechanisms that underlie the detrimental effects of UVB on the cellular defense systems of the Antarctic copepod T. kingsejongensis.

Genome-wide identification of heat shock proteins in harpacticoid, cyclopoid, and calanoid copepods: Potential application in marine ecotoxicology

Constant evolution of omics-technologies has provided access to identification of various important gene families. Recently, genome assemblies on widely used ecotoxicological model species, including rotifers and copepods have been completed and representative detoxification-related gene families have been discovered for biomarker genes. However, despite ubiquitous presence of stress-response proteins, limited information on full genome-wide report on heat shock proteins (Hsps) is available. Various studies have demonstrated multiple cellular functions of Hsps in living organisms as an important biomarker in response to abiotic and biotic stressors, however, full genome-wide identification of Hsps, particularly in aquatic invertebrates, has not been reported. This is the first study to report the entire Hsps and basal gene expression levels in three regional-specific copepods: Tigriopus japonicus and kingsejongensis, Paracyclopina nana, and Eurytemora affnis, and how each Hsp family gene is regulated at a basal level.

Hydrogen peroxide treatment and its impacts on Lepeophtheirus salmonis originating from the Bay of Fundy, Canada

Hydrogen peroxide (H₂ O₂ ) is used to treat sea lice infections of farmed salmonids in the Atlantic and Pacific Oceans and issues with resistance to this treatment, and others are a major threat to the sustainability of the industry. The objectives of this study were to determine how H₂ O₂ exposure affects survival and antioxidant-related gene expression in salmon lice (Lepeophtheirus salmonis) collected from the Bay of Fundy, New Brunswick. The maximum recommended dose of H₂ O₂ is 1,800 mg/L, while the EC₅₀ values (with 95% CI) for the population tested were 1,486 (457, 2,515) mg/L for males and 2,126 (984, 3,268) mg/L for females. Neither temperature nor pretreatment with emamectin benzoate (EMB) impacted survival after H₂ O₂ exposure. RT-qPCR was performed on pre-adult sea lice exposed to H₂ O₂ and showed that four genes classically involved in the response to oxidative stress were unchanged between treated and control groups. Seven genes were found to be significantly upregulated in males and one in females. This is the first report on the efficacy and molecular responses of Atlantic Canada sea lice to H₂ O₂ treatment.

Lineage-level divergence of copepod glycerol transporters and the emergence of isoform-specific trafficking regulation

Transmembrane conductance of small uncharged solutes such as glycerol typically occurs through aquaglyceroporins (Glps), which are commonly encoded by multiple genes in metazoan organisms. To date, however, little is known concerning the evolution of Glps in Crustacea or what forces might underly such apparent gene redundancy. Here, we show that Glp evolution in Crustacea is highly divergent, ranging from single copy genes in species of pedunculate barnacles, tadpole shrimps, isopods, amphipods and decapods to up to 10 copies in diplostracan water fleas although with monophyletic origins in each lineage. By contrast the evolution of Glps in Copepoda appears to be polyphyletic, with surprisingly high rates of gene duplication occurring in a genera- and species-specific manner. Based upon functional experiments on the Glps from a parasitic copepod (Lepeophtheirus salmonis), we show that such lineage-level gene duplication and splice variation is coupled with a high rate of neofunctionalization. In the case of L. salmonis, splice variation of a given gene resulted in tissue- or sex-specific expression of the channels, with each variant evolving unique sites for protein kinase C (PKC)- or protein kinase A (PKA)-regulation of intracellular membrane trafficking. The combined data sets thus reveal that mutations favouring a high fidelity control of intracellular trafficking regulation can be a selection force for the evolution and retention of multiple Glps in copepods.

Genetic homogeneity coupled with morphometric variability suggests high phenotypic plasticity in the sea louse Caligus rogercresseyi (Boxshall and Bravo, 2000), infecting farmed salmon (Salmo salar) along a wide latitudinal range in southern Chile

The sea louse Caligus rogercresseyi is the most important pathogen causing "caligidosis" in the Chilean salmon industry. In this study, using cox1 gene, we evaluate the genetic variation of C. rogercresseyi from farmed Salmo salar along a latitudinal range (40°-52°S) in south Chile to determine whether morphological differences are explained by genetic or environmental factors. Female parasites were randomly collected from S. salar at five farms. Body variation was examined using multivariate analyses and genetic heterogeneity was explored with AMOVA. C. rogercresseyi exhibited significant morphometric variability among sites and parasites collected from >54°S were the longest ones. Parasites did not show genetic structure among farms. Thus, C. rogercresseyi infesting salmons is panmictic along an extensive latitudinal range in south Chile. The same genetic pattern can be explained by the frequent movement of parasitized S. salar among farms in that region. Phenotypic plasticity in parasites could be explained by natural or aquaculture-mediated environment variability. C. rogercreseyi from 54°S could favor the local spread of this disease, suggesting an immediate health risk for the recent salmon industry in that region. Further research is required to confirm genetic homogeneity of this parasite along its geographical distribution using more powerful markers (e.g. SNPs).

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.

An Experimental Test of Adaptive Introgression in Locally Adapted Populations of Splash Pool Copepods

As species struggle to keep pace with the rapidly warming climate, adaptive introgression of beneficial alleles from closely related species or populations provides a possible avenue for rapid adaptation. We investigate the potential for adaptive introgression in the copepod, Tigriopus californicus, by hybridizing two populations with divergent heat tolerance limits. We subjected hybrids to strong heat selection for 15 generations followed by whole-genome resequencing. Utilizing a hybridize evolve and resequence (HER) technique, we can identify loci responding to heat selection via a change in allele frequency. We successfully increased the heat tolerance (measured as LT50) in selected lines, which was coupled with higher frequencies of alleles from the southern (heat tolerant) population. These repeatable changes in allele frequencies occurred on all 12 chromosomes across all independent selected lines, providing evidence that heat tolerance is polygenic. These loci contained genes with lower protein-coding sequence divergence than the genome-wide average, indicating that these loci are highly conserved between the two populations. In addition, these loci were enriched in genes that changed expression patterns between selected and control lines in response to a nonlethal heat shock. Therefore, we hypothesize that the mechanism of heat tolerance divergence is explained by differential gene expression of highly conserved genes. The HER approach offers a unique solution to identifying genetic variants contributing to polygenic traits, especially variants that might be missed through other population genomic approaches.

Highlighting the promise of qPCR-based environmental monitoring: response of the ribosomal RNA:DNA ratio of calanoid copepods to toxic cyanobacteria

Calanoid copepods are integral to aquatic food webs and may drive the bioaccumulation of toxins and heavy metals, spread of infectious diseases, and occurrence of toxic cyanobacterial harmful algal blooms (HABs) in freshwater aquatic systems. However, interrelationships between copepod and cyanobacterial population dynamics and ecophysiology remain unclear. Insights into these relationships are important to aquatic resource management, as they may help guide mitigation efforts. We developed a calanoid copepod qPCR assay to investigate how copepod abundance and physiological status relate to the abundance of cyanobacteria and the concentration of total microcystin in a HAB-prone freshwater multi-use eutrophic lake. Through in silico and in vitro validation of primers and analyses of time series, we demonstrate that our assay can be used as a reliable tool for environmental monitoring. Importantly, copepod RNA:DNA ratios on and shortly after the day when microcystin concentration was at its highest within the lake were not significantly lower (or higher) than before or after this period, suggesting that copepods may have been tolerant of microcystin levels observed and capable of perpetuating bloom events by consuming competitors of toxic cyanobacteria.

Chromosome-scale genome assembly of the sea louse Caligus rogercresseyi by SMRT sequencing and Hi-C analysis

Caligus rogercresseyi, commonly known as sea louse, is an ectoparasite copepod that impacts the salmon aquaculture in Chile, causing losses of hundreds of million dollars per year. In this study, we report a chromosome-scale assembly of the sea louse (C. rogercresseyi) genome based on single-molecule real-time sequencing (SMRT) and proximity ligation (Hi-C) analysis. Coding RNAs and non-coding RNAs, and specifically long non-coding RNAs (lncRNAs) and microRNAs (miRNAs) were identified through whole transcriptome sequencing from different life stages. A total of 23,686 protein-coding genes and 12,558 non-coding RNAs were annotated. In addition, 6,308 lncRNAs and 5,774 miRNAs were found to be transcriptionally active from larvae to adult stages. Taken together, this genomic resource for C. rogercresseyi represents a valuable tool to develop sustainable control strategies in the salmon aquaculture industry.

Investigation of deltamethrin resistance in salmon lice (Lepeophtheirus salmonis) provides no evidence for roles of mutations in voltage-gated sodium channels

BACKGROUND

The pyrethroid deltamethrin is used to treat infestations of farmed salmon by parasitic salmon lice, Lepeophtheirus salmonis (Krøyer). However, the efficacy of deltamethrin for salmon delousing is threatened by resistance development. In terrestrial arthropods, knockdown resistance (kdr) mutations of the voltage-gated sodium channel (Na_v ), the molecular target for pyrethroids, can cause deltamethrin resistance. A putative kdr mutation of an L. salmonis sodium channel homologue (LsNa_v 1.3 I936V) has been identified previously. At the same time, deltamethrin resistance of L. salmonis has been shown to be inherited maternally and to be associated with mitochondrial DNA (mtDNA) mutations. This study assessed potential roles of the above putative kdr mutation as a determinant of deltamethrin resistance in laboratory strains and field populations of L. salmonis.

RESULTS

The deltamethrin-resistant L. salmonis strain IoA-02 expresses the LsNa_v 1.3 I936V mutation but was susceptible to the non-ester pyrethroid etofenprox, a compound against which pyrethroid-resistant arthropods are usually cross-resistant if resistance is caused by Na_v mutations. In a family derived from a cross between an IoA-02 male and a drug-susceptible female lacking the kdr mutation, deltamethrin resistance was not associated with the genotype at the LsNa_v 1.3 locus (P > 0.05). Similarly, in Scottish field populations of L. salmonis, LsNa_v 1.3 I936V showed no association with deltamethrin resistance. By contrast, genotypes at the mtDNA loci A14013G and A9030G were significantly associated with deltamethrin resistance (P < 0.001).

CONCLUSION

In the studied L. salmonis isolates, deltamethrin resistance was unrelated to the LsNa_v 1.3 I936V mutation, but showed close association with mtDNA mutations.

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.

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.

Combined effects of exposure to engineered silver nanoparticles and the water-soluble fraction of crude oil in the marine copepod Calanus finmarchicus

While it is likely that ENPs may occur together with other contaminants in nature, the combined effects of exposure to both ENPs and environmental contaminants are not studied sufficiently. In this study, we investigated the acute and sublethal toxicity of PVP coated silver nanoparticles (AgNP) and ionic silver (Ag⁺; administered as AgNO₃) to the marine copepod Calanus finmarchicus. We further studied effects of single exposures to AgNPs (nominal concentrations: low 15 μg L^-1 NPL, high 150 μg L^-1 NPH) or Ag⁺ (60 μg L^-1), and effects of co-exposure to AgNPs, Ag⁺ and the water-soluble fraction (WSF; 100 μg L^-1) of a crude oil (AgNP + WSF; Ag⁺+WSF). The gene expression and the activity of antioxidant defense enzymes SOD, CAT and GST, as well as the gene expression of HSP90 and CYP330A1 were determined as sublethal endpoints. Results show that Ag⁺ was more acutely toxic compared to AgNPs, with 96 h LC₅₀ concentrations of 403 μg L^-1 for AgNPs, and 147 μg L^-1 for Ag⁺. Organismal uptake of Ag following exposure was similar for AgNP and Ag⁺, and was not significantly different when co-exposed to WSF. Exposure to AgNPs alone caused increases in gene expressions of GST and SOD, whereas WSF exposure caused an induction in SOD. Responses in enzyme activities were generally low, with significant effects observed only on SOD activity in NPL and WSF exposures and on GST activity in NPL and NPH exposures. Combined AgNP and WSF exposures caused slightly altered responses in expression of SOD, GST and CYP330A1 genes compared to the single exposures of either AgNPs or WSF. However, there was no clear pattern of cumulative effects caused by co-exposures of AgNPs and WSF. The present study indicates that the exposure to AgNPs, Ag⁺, and to a lesser degree WSF cause an oxidative stress response in C. finmarchicus, which was slightly, but mostly not significantly altered in combined exposures. This indicated that the combined effects between Ag and WSF are relatively limited, at least with regard to oxidative stress.

Warmer temperature increases mercury toxicity in a marine copepod

Marine biota have been co-challenged with ocean warming and mercury (Hg) pollution over many generations because of human activities; however, the molecular mechanisms to explain their combined effects are not well understood. In this study, a marine planktonic copepod Pseudodiaptomus annandalei was acutely exposed to different temperature (22 and 25 °C) and Hg (0 and 118 μg/L) treatments in a 24-h cross-factored experiment. Hg accumulation and its subcellular fractions were determined in the copepods after exposure. The expression of the genes of superoxide dismutase (SOD), glutathione peroxidase (GPx), metallothionein1 (mt1), heat shock protein 70 (hsp70), hsp90, hexokinase (hk), and pyruvate kinase (pk) was also analyzed. Both the Hg treatment alone and the combined exposure of warmer temperature plus Hg pollution remarkably facilitated Hg bioaccumulation in the exposed copepods. Compared with the Hg treatment alone, the combined exposure increased total Hg accumulation and also the amount of Hg stored in the metal-sensitive fractions (MSF), suggesting elevated Hg toxicity in P. annandalei under a warmer environment, given that the MSF is directly related to metal toxicity. The warmer temperature significantly up-regulated the mRNA levels of mt1, hsp70, hsp90, and hk, indicating the copepods suffered from thermal stress. With exposure to Hg, the mRNA level of SOD increased strikingly but the transcript levels of hsp90, hk, and pk decreased significantly, indicating that Hg induced toxic events (e.g., oxidative damage and energy depletion). Particularly, in contrast to the Hg treatment alone, the combined exposure significantly down-regulated the mRNA levels of SOD and GPx but up-regulated the mRNA levels of mt1, hsp70, hsp90, hk, and pk. Collectively, the results of this study indicate that ocean warming will potentially boost Hg toxicity in the marine copepod P. annandalei, which is information that will increase the accuracy of the projections of marine ecosystem responses to the joint effects of climate change stressors and metal pollution on the future ocean.

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

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

A pond-side test for Guinea worm: Development of a loop-mediated isothermal amplification (LAMP) assay for detection of Dracunculus medinensis

Guinea worm Dracunculus medinensis causes debilitating disease in people and is subject to an ongoing global eradication programme. Research and controls are constrained by a lack of diagnostic tools. We developed a specific and sensitive LAMP method for detecting D. medinensis larval DNA in copepod vectors. We were able to detect a single larva in a background of field-collected copepods. This method could form the basis of a “pond-side test” for detecting potential sources of Guinea worm infection in the environment, in copepods, including in the guts of fish as potential transport hosts, enabling research, surveillance and targeting of control measures. The key constraint on the utility of this assay as a field diagnostic, is a lack of knowledge of variation in the temporal and spatial distribution of D. medinensis larvae in copepods in water bodies in the affected areas and how best to sample copepods to obtain a reliable diagnostic sample. These fundamental knowledge gaps could readily be addressed with field collections of samples across areas experiencing a range of worm infection frequencies, coupled with field and laboratory analyses using LAMP and PCR.

•

LAMP tests were developed for Dracunculus medinensis and D. insignis and were shown to be specific and sensitive.

•

A LAMP test was developed to amplify DNA from copepods to use as an internal amplification control during testing.

•

Samples of copepods taken from ponds could be tested using the LAMP tests and Dracunculus medinensis could be detected.

•

Results are achieved in less than 30 min using just the Genie III instrument and no other laboratory equipment is necessary.