Redescription of Taeniacanthus aulacocephali (Copepoda: Cyclopoida: Taeniacanthidae) parasitic on Uranoscopus japonicus Houttuyn (Uranoscopidae) from Pacific coast, Japan, with adaptation of the non-destructive DNA technique.

Taeniacanthus aulacocephali Izawa, 2021 (Copepoda: Cyclopoida: Taeniacanthidae) was redescribed from the branchial cavity and gill filaments of Uranoscopus japonicus Houttuyn (Perciformes: Uranoscopidae) collected from the Pacific coast of the Kochi and Wakayama prefectures, Japan. This is the second record of the copepod, and the finding from U. japonicus represents the new host record. The species is characterized by several distinguishing features: 1) a decrease in the width of the habitus between the second and fourth pedigerous segments; 2) the ratio of prosome/body length; 3) the presence of eight setae on the exopodal terminal segment of leg 2; 4) an un-bifurcated maxilliped claw surrounded by 14-28 transverse ridges; and 5) the presence of an inner coxal seta on legs 2 and 3. The newly collected specimens were subjected to a modified non-destructive DNA extraction method and morphological description based on the same copepod individual, while preserving a morphologically describable specimen. Sequences of 18S rDNA, 28S rDNA and the mitochondrial cytochrome c oxidase subunit 1 mitochondrial gene (cox1) were obtained.

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

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

Elevated temperature as a dominant driver to aggravate cadmium toxicity: Investigations through toxicokinetics and omics.

Despite the great interest in the consequences of global change stressors on marine organisms, their interactive effects on cadmium (Cd) bioaccumulation/biotoxicity are very poorly explored, particularly in combination with the toxicokinetic model and molecular mechanism. According to the projections for 2100, this study investigated the impact of elevated pCO2 and increased temperature (isolated or joint) on Cd uptake dynamics and transcriptomic response in the marine copepod Tigriopus japonicus. Toxicokinetic results showed significantly higher Cd uptake in copepods under increased temperature and its combination with elevated pCO2 relative to the ambient condition, linking to enhanced Cd bioaccumulation. Transcriptome analysis revealed that, under increased temperature and its combination with elevated pCO2, up-regulated expression of Cd uptake-related genes but down-regulation of Cd exclusion-related genes might cause increased cellular Cd level, which not only activated detoxification and stress response but also induced oxidative stress and concomitant apoptosis, demonstrating aggravated Cd biotoxicity. However, these were less pronouncedly affected by elevated pCO2 exposure. Therefore, temperature seems to be a primary factor in increasing Cd accumulation and its toxicity in the future ocean. Our findings suggest that we should refocus the interactive effects between climate change stressors and Cd pollution, especially considering temperature as a dominant driver.

Multigenerational effects of glyphosate-based herbicide and emamectin benzoate insecticide on the reproduction and gene expression of the copepod Pseudodiaptomus annandalei (Sewell, 1919).

This study investigates the effects of glyphosate-based herbicide (GLY) and pure emamectin benzoate (EB) insecticide on the brackish copepod Pseudodiaptomus annandalei. The 96h median lethal concentration (96 h LC50) was higher in the GLY exposure (male: 3420.96 ± 394.67 µg/L; female: 3093.46 ± 240.67 µg/L) than in the EB (male: 79.10 ± 7.30 µg/L; female: 6.38 ± 0.72 µg/L). Based on the result of 96h LC50, we further examined the effects of GLY and EB exposures at sub-lethal concentrations on the naupliar production of P. annandalei. Subsequently, a multigenerational experiment was conducted to assess the long-term impact of GLY and EB at concentrations 375 µg/L, and 0.025 µg/L respectively determined by sub-lethal exposure testing. During four consecutive generations, population growth, clutch size, prosome length and width, and sex ratio were measured. The copepods exposed to GLY and EB showed lower population growth but higher clutch size than the control group in most generations. Gene expression analysis indicated that GLY and EB exposures resulted in the downregulation of reproduction-related (vitellogenin) and growth-related (myosin heavy chain) genes, whereas a stress-related gene (heat shock protein 70) was upregulated after multigenerational exposure. The results of the toxicity test after post-multigenerational exposure indicated that the long-term GLY-exposed P. annandalei displayed greater vulnerability towards GLY toxicity compared to newly-exposed individuals. Whereas, the tolerance of EB was significantly higher in the long-term exposed copepod than in newly-exposed individuals. This suggests that P. annandalei might have greater adaptability towards EB toxicity than towards GLY toxicity. This study reports for the first time the impacts of common pesticides on the copepod P. annandalei, which have implications for environmental risk assessment and contributes to a better understanding of copepod physiological responses towards pesticide contaminations.

Molecular evidence for a new endemic species of Acartia (Copepoda, Calanoida) from the Southeast Pacific coast.

The loss of biodiversity in marine populations is one of the consequences of the increased events of extreme environmental conditions in the oceans, which can condition the persistence of populations to future scenarios of climate change. Therefore, it is extremely necessary to explore and monitor the genetic diversity of natural populations. In the Southeast Pacific Ocean (SEPO), specifically on the coast of Chile, the presence of the copepod Acartia tonsa has been indicated solely using morphological evidence, due to the absence of genetic information. In the present work, the genetic diversity, population structure and phylogenetic position within the genus Acartia, of populations identified morphologically as A. tonsa, was evaluated by amplification of the mitochondrial cytochrome c oxidase subunit I and nuclear marker 18 s. Our results showed that the populations identified as A. tonsa correspond to a new monophyletic group endemic to SEPO (GMYC = 1.00; PTP = 0.95). The populations showed moderate to high genetic diversity with an incipient structuring between populations and biogeographic zones. Our results suggest that despite the homogenizing effect of the Humboldt Current, isolation by distance and contrasting environmental conditions at different geographic scales have an important influence on the genetic diversity of zooplankton in the SEPO region.

Insights from Hi-C data regarding the Pacific salmon louse (Lepeophtheirus salmonis) sex chromosomes.

Salmon lice, Lepeophtheirus salmonis (family Caligidae), are ectoparasites that have negatively impacted the salmon aquaculture industry and vulnerable wild salmon populations. Researchers have studied salmon lice to better understand their biology to develop effective control strategies. In this study, we updated the chromosome-level reference genome assembly of the Pacific subspecies of L. salmonis using Hi-C data. The previous version placed contigs/scaffolds using an Atlantic salmon louse genetic map. By utilizing Hi-C data from Pacific salmon lice, we were able to properly assign locations to contigs/scaffolds previously unplaced or misplaced. This resulted in a more accurate genome assembly and a more comprehensive characterization of the sex chromosome unique to females (W). We found evidence that the same ZW-ZZ mechanism is common in both Atlantic and Pacific subspecies of salmon lice using PCR assays. The W chromosome was approximately 800 kb in size, which is ∼30 times smaller than the Z chromosome (24 Mb). The W chromosome contained 61 annotated genes, including 32 protein-coding genes, 27 long noncoding RNA (lncRNA) genes, and 2 pseudogenes. Among these 61 genes, 39 genes shared homology to genes found on other chromosomes, while 20 were unique to the W chromosome. Two genes of interest on the W chromosome, prohibitin-2 and kinase suppressor of ras-2, were previously identified as potential sex-linked markers in the salmon louse. However, we prioritized the 20 unique genes on the W chromosome as sex-determining candidates. This information furthers our understanding of the biology of this ectoparasite and may help in the development of more effective management strategies.

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.

Molecular characterization of Lernathropus kroyeri from intensive aquaculture and pathophysiology of infested sea bass.

The copepod Lernathropus kroyeri constitutes one of the major parasites for the Mediterranean aquaculture, infesting the sea bass Dicentrarchus labrax causing thus disruptions of growth performance and occasionally mortalities. Despite the large spread and the high frequency of this parasite in mariculture farms of Eastern Mediterranean, L. kroyeri genetic profile from aquaculture as well as the pathophysiological response of D. labrax have not been studied so far. Keeping this in mind, in the present study we investigated the L. kroyeri infestation on D. labrax from two farms in Greece, examining both healthy and heavy parasitized individuals. Assays included histopathology, phylogenetic reconstruction of the parasite and physiological response of the fish by the means of antioxidant, inflammatory metabolic and stress related gene expression analysis at both mRNA and protein levels. Genetic analysis indicated that L. kroyeri composes a monophyletic group, highly phylogenetically distant from other congeneric groups. Heavy infested D. labrax witnessed a significantly increased immune response that further led to oxidative stress and metabolic alterations. Overall, our results demonstrate the, seasonally independent, high infestation of this parasitic copepods, which continue to affect Mediterranean intensive aquaculture systems.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.