The complete mitochondrial genome of the cyclopoid copepod Paracyclopina nana: a highly divergent genome with novel gene order and atypical gene numbers

Nanoplastics induce epigenetic signatures of transgenerational impairments associated with reproduction in copepods under ocean acidification

Ocean acidification (OA) is one of many major global climate changes that pose a variety of risks to marine ecosystems in different ways. Meanwhile, there is growing concern about how nanoplastics (NPs) affect marine ecosystems. Combined exposure of marine organisms to OA and NPs is inevitable, but their interactive effects remain poorly understood. In this study, we investigated the multi- and transgenerational toxicity of NPs on copepods under OA conditions for ten generations. The findings revealed that OA and NPs have a synergistic negative effect on copepod reproduction across generations. In particular, the transgenerational groups showed reproductive impairments in the F1 and F2 generations (F1T and F2T), even though they were never exposed to NPs. Moreover, our epigenetic examinations demonstrated that the observed intergenerational reproductive impairments are associated with differential methylation patterns of specific genes, suggesting that the interaction of OA and NPs can pose a significant threat to the sustainability of copepod populations through epigenetic modifications. Overall, our findings provide valuable insight into the intergenerational toxicity and underlying molecular mechanisms of responses to NPs under OA conditions.

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.

Geographic Variation of Phyllodiaptomus tunguidus Mitogenomes: Genetic Differentiation and Phylogeny

Phyllodiaptomus tunguidus (Copepoda: Calanoida) is largely endemic to and widespread in freshwater in southern China, where it inhabits a complex landscape from lowland to highland across an elevation gradient of 2000m. A deep genetic differentiation can be expected between its most distant geographic populations. Here, we sequenced nine mitogenomes from diverse populations. All mitogenomes contained 37 genes, including 13 protein-coding genes (PCG), two rRNA genes, 22 tRNA genes and one control region. Their base composition, genetic distance and tRNA structure indeed revealed a wide differentiation between mitogenomes. Two P. tunguidus from Guangxi near Vietnam differed from the other seven by up to 10.1%. Their tRNA-Arg had a complete clover-leaf structure, whereas that of the others did not contain an entire dihydrouridine arm. The nine mitogenomes also differed in the length of rRNA. NJ, ML, and Bayesian analyses all split them into two clades, viz. the two P. tunguidus from Guangxi (Clade 1), and the other seven (Clade 2). Both the structure and phylogeny of the mitogenomes suggest that P. tunguidus has complex geographic origin, and its populations in Clade 1 have long lived in isolation from those in Clade 2. They currently reach the level of subspecies or cryptic species. An extensive phylogenetic analysis of Copepoda further verified that Diaptomidae is the most recently diverging family in Calanoida and that P. tunguidus is at the evolutionary apex of the family.

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.

Interspecific biotransformation and detoxification of arsenic compounds in marine rotifer and copepod

The toxicity of arsenic (As) has been reported to be different depending on their chemical forms. However, its toxicity mechanisms largely remain unknown. In this study, to investigate toxicity mechanism of As in marine zooplanktons, namely, the rotifer Brachionus plicatilis and the copepod Paracyclopina nana, metabolites of As were analyzed by using a high-performance liquid chromatography coupled with inductively coupled plasma mass spectrometry with in vivo toxicity and antioxidant responses in response to inorganic As, including arsenate (As^V) and arsenite (As^III). While As^III was more toxic than As^V in both organisms, the rotifer B. plicatilis exhibited stronger tolerance, compared to the copepod P. nana. The As speciation analysis revealed differences in biotransformation processes in two species with B. plicatilis having a more simplified process than P. nana, contributing to a better tolerance against As in the rotifer B. plicatilis compared to P. nana. Moreover, the levels of GSH content and the regulation of omega class glutathione S-transferases were different in response to oxidative stress between B. plicatilis and P. nana. These results suggest that the rotifer B. plicatilis has a unique survival strategy with more efficient biotransformation and antioxidant responses, compared to P. nana, conferring higher tolerance to As.

Effects of different routes of exposure to metals on bioaccumulation and population growth of the cyclopoid copepod Paracyclopina nana

We examined effects of the three metals cadmium (Cd), copper (Cu), and nickel (Ni) on two subpopulations of the cyclopoid copepod Paracyclopina nana. We sought to investigate the effects of metal exposure on population growth and structure of P. nana and to understand the parameters affecting the metal bioaccumulation in copepods. A first experiment tested the hypothesis of competition between these metals in a mixture using a P. nana mass culture in 10 L beakers with the sublethal concentrations (1/3 of LC50) as determined for E. affinis. A second experiment pursued the same with a P. nana population which was adapted to a higher Cu concentration for several generations (226.9 ± 15.9 μg g^-1 dw Cu in copepods) and using the proper sublethal concentrations for P. nana. After 96 h of exposure, results from the first experiment showed a decreasing population growth and instead of an increasing metal accumulation in copepods. Cd also appeared to be more accumulated when it was alone, confirming the hypothesis of metal competition in mixture. Results from the second experiment revealed less marked effects. When metal concentrations increased in the treatment it decreased in copepods, indicating depuration activity in the population already adapted to metal exposure. This paper is the first one investigating the parameters affecting the bioaccumulation capacity of P. nana in response to metals. It offers a better understanding of copepod responses to metal contamination in a complex aquatic environment.

The complete mitochondrial genomes of two vent squat lobsters, Munidopsis lauensis and M. verrilli: Novel gene arrangements and phylogenetic implications

Hydrothermal vents are considered as one of the most extremely harsh environments on the Earth. In this study, the complete mitogenomes of hydrothermal vent squat lobsters, Munidopsis lauensis and M. verrilli, were determined through Illumina sequencing and compared with other available mitogenomes of anomurans. The mitogenomes of M. lauensis (17,483 bp) and M. verrilli (17,636 bp) are the largest among all Anomura mitogenomes, while the A+T contents of M. lauensis (62.40%) and M. verrilli (63.99%) are the lowest. The mitogenomes of M. lauensis and M. verrilli display novel gene arrangements, which might be the result of three tandem duplication-random loss (tdrl) events from the ancestral pancrustacean pattern. The mitochondrial gene orders of M. lauensis and M. verrilli shared the most similarities with S. crosnieri. The phylogenetic analyses based on both gene order data and nucleotide sequences (PCGs and rRNAs) revealed that the two species were closely related to Shinkaia crosnieri. Positive selection analysis revealed that eighteen residues in seven genes (atp8, Cytb, nad3, nad4, nad4l, nad5, and nad6) of the hydrothermal vent anomurans were positively selected sites.

Genome-wide identification and expression of the entire 52 glutathione S-transferase (GST) subfamily genes in the Cu2+-exposed marine copepods Tigriopus japonicus and Paracyclopina nana

In this study, the entire glutathione S-transferases (GSTs), the major phase II detoxification enzyme, were identified in two marine copepod species Tigriopus japonicus and Paracyclopina nana. The genome-wide identification of GSTs in T. japonicus and P. nana resulted in 32 and 20 GSTs in total, respectively. Among the identified GSTs, two specific classes of GSTs, specifically sigma and delta/epsilon GSTs were the dominant form of cytosolic GSTs in T. japonicus, while delta/epsilon and mu classes were dominant cytosolic GSTs in P. nana. In addition, Membrane-Associated Proteins in Eicosanoid and Glutathione metabolism (MAPEG) family were found in relatively higher proportion compared to other classes. Moreover, sigma, delta/epsilon, and microsomal GSTs have shown to expand through tandem duplication. To validate the detoxification function of the identified GSTs, both copepods were exposed to copper (Cu²⁺) and the reactive oxygen species (ROS) level and GST activity were measured. With integration of phylogenetic analysis and xenobiotic-mediated GST mRNA expression patterns along with previous enzymatic activities, the functional divergence among species-specific GST genes was clearly observed. This study covers full identification of GST classes in two marine copepod species and their important role in marine environmental ecotoxicology.

Adverse effects of a synthetic pyrethroid insecticide cypermethrin on life parameters and antioxidant responses in the marine copepods Paracyclopina nana and Tigriopus japonicus

To find the adverse effects induced by cypermethrin, the ecotoxicological model copepods Tigriopus japonicus and Paracyclopina nana were exposed under cypermethrin, which is a widely used type-II pyrethroid insecticide in agriculture. Despite its large-scale application as insecticide in agriculture, little information is available on its adverse effects on aquatic invertebrates. In this study, the toxicity of cypermethrin on two copepods was assessed based on life parameters (growth rate and reproduction), oxidative stress and consequent antioxidant enzymatic activities, and gene expression profiles of antioxidants. The acute toxicity alone demonstrated that P. nana is less tolerant and sensitive against cypermethrin, compared to T. japonicus. However, under chronic exposure, life parameters of both P. nana and T. japonicus were severely affected by cypermethrin. Among antioxidant enzymatic activities, superoxide dismutase (SOD) and glutathione S-transferase (GST), in particular, demonstrated significant increase in response to cypermethrin. Furthermore, temporal-mRNA expression profile showed modulations in antioxidant related genes in response to cypermethrin. Our results provide the underlying toxic mechanism of cypermethrin and the oxidative stress defense mechanism and species-specific tolerance against cypermethrin in two model copepods species.

Complete mitochondrial genome of the first deep-sea spongicolid shrimp Spongiocaris panglao (Decapoda: Stenopodidea): Novel gene arrangement and the phylogenetic position and origin of Stenopodidea

Stenopodidea Claus, 1872 (Crustacea: Decapoda) is one of the major groups of decapods crustaceans. Hitherto, only one complete mitochondrial genome (mitogenome) from the family Stenopodidae is available for the infraorder Stenopodidea. Here, we determined the complete mitogenome of Spongiocaris panglao de Grave and Saito, 2016 using Illumina sequencing, representing the first species from the family Spongicolidae. The 15,909 bp genome is a circular molecule and consists of 13 protein-coding genes (PCGs), 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes and one control region. Although the overall genome organization is typical for metazoans, the mitogenome of S. panglao shows some derived characters. A + T content of 77.42% in S. pamglao mitogenome is second-highest among the dacapods described to date. The trnR gene exhibit modified secondary structure with the TψC loop completely missing, which might be a putative autapomorphy of S. pamglao mitogenome. Compared with the shallow-water stenopodidean species S. hispidus, the control region of S. pamglao exhibits three characteristics: larger size, higher A + T content, and more tandem repeat sequences. The gene order exhibited difference from the ancestral mitogenome pattern of the Pancrustacea, with 5 tRNA genes rearrangement. The result from BI was agreed with most morphological characters and molecular evidences, revealing that Stenopodidea and Reptantia had the closest relationship, as the sister group of Caridea. Still, the alternative hypothesis supported from ML topology cannot be completely rejected based on the current data. Estimated times revealed that the two stenopodideans families Stenopodidae and Spongicolidae diverged from each other around 122 Mya. The divergence time of spongicolid shrimp is in good agreement with the origin of their hexactinellid hosts (78-144 Mya).

Mitochondrial genomes of the key zooplankton copepods Arctic Calanus glacialis and North Atlantic Calanus finmarchicus with the longest crustacean non-coding regions

We determined the nearly complete mitochondrial genomes of the Arctic Calanus glacialis and its North Atlantic sibling Calanus finmarchicus, which are key zooplankton components in marine ecosystems. The sequenced part of C. glacialis mitogenome is 27,342 bp long and consists of two contigs, while for C. finmarchicus it is 29,462 bp and six contigs, what makes them the longest reported copepod mitogenomes. The typical set of metazoan mitochondrial genes is present in these mitogenomes, although the non-coding regions (NCRs) are unusually long and complex. The mitogenomes of the closest species C. glacialis and C. finmarchicus, followed by the North Pacific C. sinicus, are structurally similar and differ from the much more typical of deep-water, Arctic C. hyperboreus. This evolutionary trend for the expansion of NCRs within the Calanus mitogenomes increases mitochondrial DNA density, what resulted in its similar density to the nuclear genome. Given large differences in the length and structure of C. glacialis and C. finmarchicus mitogenomes, we conclude that the species are genetically distinct and thus cannot hybridize. The molecular resources presented here: the mitogenomic and rDNA sequences, and the database of repetitive elements should facilitate the development of genetic markers suitable in pursuing evolutionary research in copepods.

Adverse effects of methylmercury (MeHg) on life parameters, antioxidant systems, and MAPK signaling pathways in the rotifer Brachionus koreanus and the copepod Paracyclopina nana

To evaluate the adverse effects of MeHg on the rotifer Brachionus koreanus and the copepod Paracyclopina nana, we assessed the effects of MeHg toxicity on life parameters (e.g. growth retardation and fecundity), antioxidant systems, and mitogen-activated protein kinase (MAPK) signaling pathways at various concentrations (1ng/L, 10ng/L, 100ng/L, 500ng/L, and 1000ng/L). MeHg exposure resulted in the growth retardation with the increased ROS levels but decreased glutathione (GSH) levels in a dose-dependent manner in both B. koreanus and P. nana. Antioxidant enzymatic activities (e.g. glutathione S-transferase [GST], glutathione reductase [GR], and glutathione peroxidase [GPx]) in B. koreanus showed more positive responses compared the control but in P. nana, those antioxidant enzymatic activities showed subtle changes due to different no observed effect concentration (NOEC) values among the two species. Expression of antioxidant genes (e.g. superoxide dismutase [SOD], GSTs, glutathione peroxidase [GPx], and catalase [CAT]) also demonstrated similar effects as shown in antioxidant enzymatic activities. In B. koreanus, the level of p-ERK was decreased in the presence of 1000ng/L MeHg, while the levels of p-ERK and p-p38 in P. nana were reduced in the presence of 10ng/L MeHg. However, p-JNK levels were not altered by MeHg in B. koreanus and P. nana, compared to the corresponding controls. In summary, life parameters (e.g. reduced fecundity and survival rate) were closely associated with effects on the antioxidant system in response to MeHg. These observations provide a better understanding on the adverse effects of MeHg on in vivo life parameters and molecular defense mechanisms in B. koreanus and P. nana.

Assessing the identity and expression level of the cytochrome P450 20A1 (CYP20A1) gene in the BPA-, BDE-47, and WAF-exposed copepods Tigriopus japonicus and Paracyclopina nana

CYP20A1 is a member of the cytochrome P450 (CYP) superfamily, identified as an orphan P450 without any assigned biological function; hence, its continued status as an "orphan" gene. In order to address this shortcoming in our understanding of this superfamily, we sought to characterize the CYP20A1 gene in the copepods Tigriopus japonicus (Tj-CYP20A1) and Paracyclopina nana (Pn-CYP20A1) at their mRNA transcriptional level. We assessed the response of this gene's expression in various developmental stages and in response to treatment with bisphenol A (BPA), 2, 2', 4, 4'-tetrabromodiphenyl ether (BDE-47), and water accommodated fractions (WAFs) of crude oil. As shown in the vertebrate CYP20A1, both Tj-CYP20A1 and Pn-CYP20A1 contained characteristic conserved motifs and domain regions (I helix, K helix and heme-binding motifs) with unusual amino acid sequences apparent in their gene structure. Also molecular characterization of the putative responsive elements in the promoter regions was performed. We observed transcriptional up-regulation of these genes during post-embryonic developmental stages including sex-specific up-regulation in adults. In addition, concentration- and time-dependent mRNA transcripts in response to xenobiotics (BPA, BDE-47, and WAFs) were seen. This study focuses on the molecular elucidation of CYP20A1 genes and their interactions with xenobiotics in the copepods T. japonicus and P. nana that provides important insight into the biological importance of CYP20A1 in invertebrates.

Adverse effects of methylmercury (MeHg) on life parameters, antioxidant systems, and MAPK signaling pathways in the copepod Tigriopus japonicus

Methylmercury (MeHg) is a concerning environmental pollutant that bioaccumulates and biomagnifies in the aquatic food web. However, the effects of MeHg on marine zooplankton are poorly understood even though zooplankton are considered key mediators of the bioaccumulation and biomagnification of MeHg in high-trophic marine organisms. Here, the toxicity of MeHg in the benthic copepod Tigriopus japonicus was assessed, and its adverse effects on growth rate and reproduction were demonstrated. Antioxidant enzymatic activities were increased in the presence of MeHg, indicating that these enzymes play an important role in the defense response to MeHg, which is regulated by a complex mechanism. Subsequent activation of different patterns of mitogen-activated protein kinase (MAPK) pathways was demonstrated, providing a mechanistic approach to understand the signaling pathways involved in the effects of MeHg. Our results provide valuable information for understanding the toxicity of MeHg and the underlying defense mechanism in response to MeHg exposure in marine zooplankton.

Ecdysone receptor (EcR) and ultraspiracle (USP) genes from the cyclopoid copepod Paracyclopina nana: Identification and expression in response to water accommodated fractions (WAFs)

Ecdysteroid hormones are pivotal in the development, growth, and molting of arthropods, and the hormone pathway is triggered by binding ecdysteroid to a heterodimer of the two nuclear receptors; ecdysone receptors (EcR) and ultraspiracle (USP). We have characterized EcR and USP genes, and their 5'-untranslated region (5'-UTR) from the copepod Paracyclopina nana, and studied mRNA transcription levels in post-embryonic stages and in response to water accommodated fractions (WAFs) of crude oil. The open reading frames (ORF) of EcR and USP were 1470 and 1287bp that encoded 490 and 429 amino acids with molecular weight of 121.18 and 105.03kDa, respectively. Also, a well conserved DNA-binding domain (DBD) and ligand-binding domain (LBD) were identified which confirmed by phylogenetic analysis. Messenger RNA transcriptional levels of EcR and USP were developmental stage-specific in early post-embryonic stages (N3-4). However, an evoked expression of USP was observed throughout copepodid stage and in adult females. WAFs (40 and 80%) were acted as an ecdysone agonist in P. nana, and elicited the mRNA transcription levels in adults. Developmental stage-specific transcriptional activation of EcR and USP in response to WAFs was observed. USP gene was down-regulated in the nauplius in response to WAF, whereas up-regulation of USP was observed in the adults. This study represents the first data of molecular elucidation of EcR and USP genes and their regulatory elements from P. nana and the developmental stage specific expression in response to WAFs, which can be used as potential biomarkers for environmental stressors with ecotoxicological evaluations in copepods.

Genome-wide identification of ATP-binding cassette (ABC) transporters and their roles in response to polycyclic aromatic hydrocarbons (PAHs) in the copepod Paracyclopina nana

The ATP-binding cassette (ABC) protein superfamily is one of the largest gene families and is highly conserved in all domains. The ABC proteins play roles in several biological processes, including multi-xenobiotic resistance (MXR), by functioning as transporters in the cellular membrane. They also mediate the cellular efflux of a wide range of substrates against concentration gradients. In this study, 37 ABC genes belonging to eight distinct subfamilies were identified in the marine copepod Paracyclopina nana and annotated based on a phylogenetic analysis. Also, the functions of P-glycoproteins (P-gp) and multidrug resistance-associated proteins (MRPs), conferring MXR, were verified using fluorescent substrates and specific inhibitors. The activities of MXR-mediated ABC proteins and their transcriptional level were examined in response to polyaromatic hydrocarbons (PAHs), main components of the water-accommodated fraction. This study increases the understanding of the protective role of MXR in response to PAHs over the comparative evolution of ABC gene families.

Adverse effects of microplastics and oxidative stress-induced MAPK/Nrf2 pathway-mediated defense mechanisms in the marine copepod Paracyclopina nana

Microplastic pollution causes a major concern in the marine environment due to their worldwide distribution, persistence, and adverse effects of these pollutants in the marine ecosystem. Despite its global presence, there is still a lack of information on the effect of microplastics on marine organisms at the molecular level. Herein we demonstrated ingestion and egestion of nano- (0.05 μm) and micro-sized (0.5 and 6 μm) polystyrene microbeads in the marine copepod Paracyclopina nana, and examined molecular responses to exposure to microbeads with in vivo endpoints such as growth rate and fecundity. Also, we proposed an adverse outcome pathway for microplastic exposure that covers molecular and individual levels. This study provides the first insight into the mode of action in terms of microplastic-induced oxidative stress and related signaling pathways in P. nana.

BDE-47 induces oxidative stress, activates MAPK signaling pathway, and elevates de novo lipogenesis in the copepod Paracyclopina nana

Brominated flame retardant, 2, 2', 4, 4'-tetrabromodiphenyl ether (BDE-47), has received grave concerns as a persistent organic pollutant, which is toxic to marine organisms, and a suspected link to endocrine abnormalities. Despite the wide distribution in the marine ecosystem, very little is known about the toxic impairments on marine organisms, particularly on invertebrates. Thus, we examined the adverse effects of BDE-47 on life history trait (development), oxidative markers, fatty acid composition, and lipid accumulation in response to BDE-47-induced stress in the marine copepod Paracyclopina nana. Also, activation level of mitogen-activated protein kinase (MAPK) signaling pathways along with the gene expression profile of de novo lipogenesis (DNL) pathways were addressed. As a result, BDE-47 induced oxidative stress (e.g. reactive oxygen species, ROS) mediated activation of extracellular signal-regulated kinase (ERK) and c-Jun-N-terminal kinase (JNK) signaling cascades in MAPK pathways. Activated MAPK pathways, in turn, induced signal molecules that bind to the transcription factors (TFs) responsible for lipogenesis to EcR, SREBP, ChREBP promoters. Also, the stress stimulated the conversion of saturated fatty acids (SFAs) to polyunsaturated fatty acids (PUFAs), a preparedness of the organism to adapt the observed stress, which could be correlated with the elongase and desaturase gene (e.g. ELO3, Δ5-DES, Δ9-DES) expressions, and then extended to the delayed early post-embryonic development and increased accumulation of lipid droplets in P. nana. This study will provide a better understanding of how BDE-47 effects on marine invertebrates particularly on the copepods, an important link in the marine food chain.

Potential of the small cyclopoid copepod Paracyclopina nana as an invertebrate model for ecotoxicity testing

Aquatic invertebrates contribute significantly to environmental impact assessment of contaminants in aquatic ecosystems. Much effort has been made to identify viable and ecologically relevant invertebrate test organisms to meet rigorous regulatory requirements. Copepods, which are ecologically important and widely distributed in aquatic organisms, offer a huge opportunity as test organisms for aquatic toxicity testing. They have a major role not only in the transfer of energy in aquatic food chains, but also as a medium of transfer of aquatic pollutants across the tropic levels. In this regard, a supratidal and benthic harpacticoid copepod Tigriopus japonicus Mori (order Harpacticoida) has shown promising characteristics as a test organism in the field of ecotoxicology. Because there is a need to standardize a battery of test organisms from species in different phylogenetic and critical ecosystem positions, it is important to identify another unrelated planktonic species for wider application and comparison. In this regard, the cyclopoid copepod Paracyclopina nana Smirnov (order Cyclopoida) has emerged as a potential test organism to meet such requirements. Like T. japonicus, it has a number of features that make it a candidate worth consideration in such efforts. Recently, the genomics of P. nana has been unraveled. Data on biochemical and molecular responses of P. nana against exposure to environmental chemicals and other stressors have been collected. Recently, sequences and expression profiles of a number of genes in P. nana encoding for heat shock proteins, xenobiotic-metabolizing enzymes, and antioxidants have been reported. These genes serve as potential biomarkers in biomonitoring of environmental pollutants. Moreover, the application of gene expression techniques and the use of its whole transcriptome have allowed evaluation of transcriptional changes in P. nana with the ultimate aim of understanding the mechanisms of action of environmental stressors. Whole-animal bioassays and gene expression studies indicate that P. nana may serve as an excellent tool to evaluate the impact of diverse disturbances in the marine environment. With a better understanding of toxicological mechanisms, ecotoxicologists will be able to understand defense mechanisms against toxicants in copepods. In this review, we illustrate the potential of P. nana as an alternative as well as a complementary invertebrate model organism for risk assessment of aquatic pollutants.

Adverse effects of MWCNTs on life parameters, antioxidant systems, and activation of MAPK signaling pathways in the copepod Paracyclopina nana

Engineered multi-walled carbon nanotubes (MWCNTs) have received widespread applications in a broad variety of commercial products due to low production cost. Despite their significant commercial applications, CNTs are being discharged to aquatic ecosystem, leading a threat to aquatic life. Thus, we investigated the adverse effect of CNTs on the marine copepod Paracyclopina nana. Additional to the study on the uptake of CNTs and acute toxicity, adverse effects on life parameters (e.g. growth, fecundity, and size) were analyzed in response to various concentrations of CNTs. Also, as a measurement of cellular damage, oxidative stress-related markers were examined in a time-dependent manner. Moreover, activation of redox-sensitive mitogen-activated protein kinase (MAPK) signaling pathways along with the phosphorylation pattern of extracellular signal-regulated kinase (ERK), p38, and c-Jun-N-terminal kinases (JNK) were analyzed to obtain a better understanding of molecular mechanism of oxidative stress-induced toxicity in the copepod P. nana. As a result, significant inhibition on life parameters and evoked antioxidant systems were observed without ROS induction. In addition, CNTs activated MAPK signaling pathway via ERK, suggesting that phosphorylated ERK (p-ERK)-mediated adverse effects are the primary cause of in vitro and in vivo endpoints in response to CNTs exposure. Moreover, ROS-independent activation of MAPK signaling pathway was observed. These findings will provide a better understanding of the mode of action of CNTs on the copepod P. nana at cellular and molecular level and insight on possible ecotoxicological implications in the marine environment.

Identification of Hox genes and rearrangements within the single homeobox (Hox) cluster (192.8 kb) of the cyclopoid copepod (Paracyclopina nana)

We report the first identification of the entire complement of the eight typical homeobox (hox) genes (lab, pb, Dfd, scr, antp, ubx, Abd-A, and Abd-B) and the ftz gene in a 192.8 kb region in the cyclopoid copepod Paracyclopina nana. A Hox3 gene ortholog was not present in the P. nana hox gene cluster, while the P. nana Dfd gene was transcribed in the opposite direction to the Daphnia pulex Dfd gene, but in the same direction as the Dfd genes of the fruit fly Drosophila melanogaster and red flour beetle Tribolium castaneum. The location of the lab and pb genes was switched in the P. nana hox cluster, while the order of the remaining hox genes was generally conserved with those of other arthropods. J. Exp. Zool. (Mol. Dev. Evol.) 9999B:XX-XX, 2016. © 2016 Wiley Periodicals, Inc.

In vivo effects of UV radiation on multiple endpoints and expression profiles of DNA repair and heat shock protein (Hsp) genes in the cycloid copepod Paracyclopina nana

To evaluate the effects of ultraviolet (UV) radiation on energy acquisition and consumption, the copepod Paracyclopina nana was irradiated with several doses (0-3kJ/m(2)) of UV. After UV radiation, we measured the re-brooding success, growth pattern of newly hatched nauplii, ingestion rate, and assimilation of diet. In addition, we checked the modulated patterns of DNA repair and heat shock protein (hsp) chaperoning genes of P. nana. UV-B radiation induced a significant reduction (7-87%) of the re-brooding rate of ovigerous females, indicating that UV-induced egg sac damage is closely correlated with a reduction in the hatching rate of UV-irradiated ovigerous female offspring. Using chlorophyll a and stable carbon isotope incubation experiments, we found a dose-dependent decrease (P<0.05) in food ingestion and the rate of assimilation to the body in response to UV radiation, implying that P. nana has an underlying ability to shift its balanced-energy status from growth and reproduction to DNA repair and adaptation. Also, expression of P. nana base excision repair (BER)-associated genes and hsp chaperoning genes was significantly increased in response to UV radiation in P. nana. These findings indicate that even 1kJ/m(2) of UV radiation induces a reduction in reproduction and growth patterns, alters the physiological balance and inhibits the ability to cope with UV-induced damage in P. nana.

Effects of UV radiation on hatching, lipid peroxidation, and fatty acid composition in the copepod Paracyclopina nana

To evaluate the effects of UV radiation on the reproductive physiology and macromolecules in marine zooplankton, several doses of UV radiation were used to treat the copepod Paracyclopina nana, and we analyzed in vivo endpoints of their life cycle such as mortality and reproductive parameters with in vitro biochemical biomarkers such as reactive oxygen species (ROS), the modulated enzyme activity of glutathione S-transferase (GST) and superoxide dismutase (SOD), and the production of a byproduct of peroxidation (e.g. malonedialdehyde, MDA). After UV radiation, the survival rate of P. nana was significantly reduced. Also, egg sac damage and a reduction in the hatching rate of offspring were observed in UV-irradiated ovigerous females. According to the assessed biochemical parameters, we found dose-dependent increases in ROS levels and high levels of the lipid peroxidation decomposition product by 2 kJ m(-2), implying that P. nana was under off-balanced status by oxidative stress-mediated cellular damage. Antioxidant enzyme activities of GST and SOD increased over different doses of UV radiation. To measure UV-induced lipid peroxidation, we found a slight reduction in the composition of essential fatty acids such as eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA). These findings indicate that UV radiation can induce oxidative stress-triggered lipid peroxidation with modulation of antioxidant enzyme activity, leading to a significant effect on mortality and reproductive physiology (e.g. fecundity). These results demonstrate the involvement of UV radiation on essential fatty acids and its susceptibility to UV radiation in the copepod P. nana compared to other species.

A description of the complete mitochondrial genomes of Amphiporus formidabilis, Prosadenoporus spectaculum and Nipponnemertes punctatula (Nemertea: Hoplonemertea: Monostilifera)

We sequenced the complete mitochondrial genomes (mitogenomes) of three Hoplonemertea species, Amphiporus formidabilis, Prosadenoporus spectaculum and Nipponnemertes punctatula, which are 14,616, 14,655 and 15,354 bp in length, respectively. Each of the three circular mitogenomes consists of 37 typical genes and some non-coding regions. The nucleotide composition of the coding strand is biased toward T, almost a half of total nucleotides in these mitogenomes. There are many poly-T tracts across these mitogenomes, which exhibit T-number variation within different clones of protein-coding genes, mainly resulting from false PCR amplification. The major non-coding regions have tandem repeat motifs and hairpin-like structures that may be associated with the initiation of replication or transcription. Data published to date for nemerteans show that Palaeonemertea species usually bear the largest mitogenomes, while representatives in the more recently derived Distromatonemertea clade bear the smallest ones; and that the gene arrangement of mitogenomes seems to be variable within the phylum Nemertea, but stable within either of Heteronemertea and Hoplonemertea.

Gamma radiation induces growth retardation, impaired egg production, and oxidative stress in the marine copepod Paracyclopina nana

Accidental nuclear radioisotope release into the ocean from nuclear power plants is of concern due to ecological and health risks. In this study, we used the marine copepod Paracyclopina nana to examine the effects of radioisotopes on marine organisms upon gamma radiation, and to measure the effects on growth and fecundity, which affect population and community structure. Upon gamma radiation, mortality (LD50 - 96 h=172 Gy) in P. nana was significantly increased in a dose-dependent manner in ovigerous P. nana females. For developmental impairment of gamma-irradiated nauplii, we observed growth retardation; in over 30 Gy-irradiated groups, offspring did not grow to adults. Particularly, over 50 Gy-irradiated ovigerous P. nana females did not have normal bilateral egg sacs, and their offspring did not develop normally to adulthood. Additionally, at over 30 Gy, we found dose-dependent increases in oxidative levels with elevated antioxidant enzyme activities and DNA repair activities. These findings indicate that gamma radiation can induce oxidative stress and DNA damage with growth retardation and impaired reproduction.

The complete mitochondrial genome of the black mud crab, Scylla serrata (Crustacea: Brachyura: Portunidae) and its phylogenetic position among (pan)crustaceans

The black mud crab, Scylla serrata (Forskål 1775), is the most economically important edible crab in South-East Asia. In the present study, the complete mitochondrial genome of black mud crab, S. serrata, was determined with the sequential polymerase chain reaction and primer walking sequencing. The complete mitochondrial genome was 15,721 bp in length with an A+T content of 69.2 % and contained 37 mitochondrial genes (13 protein coding genes (PCGs), 2 ribosomal RNA genes and 22 transfer RNA genes) and a control region (CR). The analysis of the CR sequence shows that it contains a multitude of repetitive fragments which can fold into hairpin-like or secondary structures and conserved elements as in other arthropods. The gene order of S. serrata mainly retains as the pancrustacean ground pattern, except for a single translocation of trnH. The gene arrangement of S. serrata appears to be a typical feature of portunid crabs. Phylogenetic analyses with concatenated amino acid sequences of 12 PCGs establishes that S. serrata in a well-supported monophyletic Portunidae and is consistent with previous morphological classification. Moreover, the phylogenomic results strongly support monophyletic Pancrustacea (Hexapoda plus "Crustaceans"). Within Pancrustacea, this study identifies Malacostraca + Entomostraca and Branchiopoda as the sister group to Hexapoda, which confirms that "Crustacea" is not monophyletic. Cirripedia + Remipedia appear to be a basal lineage of Pancrustacea. The present study also provides considerable data for the application of both population and phylogenetic studies of other crab species.

The complete mitochondrial genome of Galba pervia (Gastropoda: Mollusca), an intermediate host snail of Fasciola spp

Complete mitochondrial (mt) genomes and the gene rearrangements are increasingly used as molecular markers for investigating phylogenetic relationships. Contributing to the complete mt genomes of Gastropoda, especially Pulmonata, we determined the mt genome of the freshwater snail Galba pervia, which is an important intermediate host for Fasciola spp. in China. The complete mt genome of G. pervia is 13,768 bp in length. Its genome is circular, and consists of 37 genes, including 13 genes for proteins, 2 genes for rRNA, 22 genes for tRNA. The mt gene order of G. pervia showed novel arrangement (tRNA-His, tRNA-Gly and tRNA-Tyr change positions and directions) when compared with mt genomes of Pulmonata species sequenced to date, indicating divergence among different species within the Pulmonata. A total of 3655 amino acids were deduced to encode 13 protein genes. The most frequently used amino acid is Leu (15.05%), followed by Phe (11.24%), Ser (10.76%) and IIe (8.346%). Phylogenetic analyses using the concatenated amino acid sequences of the 13 protein-coding genes, with three different computational algorithms (maximum parsimony, maximum likelihood and Bayesian analysis), all revealed that the families Lymnaeidae and Planorbidae are closely related two snail families, consistent with previous classifications based on morphological and molecular studies. The complete mt genome sequence of G. pervia showed a novel gene arrangement and it represents the first sequenced high quality mt genome of the family Lymnaeidae. These novel mtDNA data provide additional genetic markers for studying the epidemiology, population genetics and phylogeographics of freshwater snails, as well as for understanding interplay between the intermediate snail hosts and the intra-mollusca stages of Fasciola spp..

Multiple rearrangements in mitochondrial genomes of Isopoda and phylogenetic implications

In this study, we analyse the evolutionary dynamics and phylogenetic implications of gene order rearrangements in five newly sequenced mitochondrial (mt) genomes and four published mt genomes of isopod crustaceans. The sequence coverage is nearly complete for four of the five newly sequenced species, with only the control region and some tRNA genes missing, while in Janira maculosa only two thirds of the genome could be determined. Mitochondrial gene order in isopods seems to be more plastic than that in other crustacean lineages, making all nine known mt gene orders different. Especially the asellote Janira is characterized by many autapomorphies. The following inferred ancestral isopod mt gene order exists slightly modified in modern isopods: nad1, tnrL1, rrnS, control region, trnS1, cob, trnT, nad5, trnF. We consider the inferred gene translocation events leading to gene rearrangements as valuable characters in phylogenetic analyses. In this first study covering major isopod lineages, potential apomorphies were identified, e.g., a shared relative position of trnR in Valvifera. We also report one of the first findings of homoplasy in mitochondrial gene order, namely a shared relative position of trnV in unrelated isopod lineages. In addition to increased taxon sampling secondary structure, modification in tRNAs and GC-skew inversion may be potentially fruitful subjects for future mt genome studies in a phylogenetic context.

Genomic resources for sea lice: analysis of ESTs and mitochondrial genomes

Sea lice are common parasites of both farmed and wild salmon. Salmon farming constitutes an important economic market in North America, South America, and Northern Europe. Infections with sea lice can result in significant production losses. A compilation of genomic information on different genera of sea lice is an important resource for understanding their biology as well as for the study of population genetics and control strategies. We report on over 150,000 expressed sequence tags (ESTs) from five different species (Pacific Lepeophtheirus salmonis (49,672 new ESTs in addition to 14,994 previously reported ESTs), Atlantic L. salmonis (57,349 ESTs), Caligus clemensi (14,821 ESTs), Caligus rogercresseyi (32,135 ESTs), and Lernaeocera branchialis (16,441 ESTs)). For each species, ESTs were assembled into complete or partial genes and annotated by comparisons to known proteins in public databases. In addition, whole mitochondrial (mt) genome sequences of C. clemensi (13,440 bp) and C. rogercresseyi (13,468 bp) were determined and compared to L. salmonis. Both nuclear and mtDNA genes show very high levels of sequence divergence between these ectoparastic copepods suggesting that the different species of sea lice have been in existence for 37-113 million years and that parasitic association with salmonids is also quite ancient. Our ESTs and mtDNA data provide a novel resource for the study of sea louse biology, population genetics, and control strategies. This genomic information provides the material basis for the development of a 38K sea louse microarray that can be used in conjunction with our existing 44K salmon microarray to study host-parasite interactions at the molecular level. This report represents the largest genomic resource for any copepod species to date.

Observing copepods through a genomic lens

Background

Copepods outnumber every other multicellular animal group. They are critical components of the world's freshwater and marine ecosystems, sensitive indicators of local and global climate change, key ecosystem service providers, parasites and predators of economically important aquatic animals and potential vectors of waterborne disease. Copepods sustain the world fisheries that nourish and support human populations. Although genomic tools have transformed many areas of biological and biomedical research, their power to elucidate aspects of the biology, behavior and ecology of copepods has only recently begun to be exploited.

Discussion

The extraordinary biological and ecological diversity of the subclass Copepoda provides both unique advantages for addressing key problems in aquatic systems and formidable challenges for developing a focused genomics strategy. This article provides an overview of genomic studies of copepods and discusses strategies for using genomics tools to address key questions at levels extending from individuals to ecosystems. Genomics can, for instance, help to decipher patterns of genome evolution such as those that occur during transitions from free living to symbiotic and parasitic lifestyles and can assist in the identification of genetic mechanisms and accompanying physiological changes associated with adaptation to new or physiologically challenging environments. The adaptive significance of the diversity in genome size and unique mechanisms of genome reorganization during development could similarly be explored. Genome-wide and EST studies of parasitic copepods of salmon and large EST studies of selected free-living copepods have demonstrated the potential utility of modern genomics approaches for the study of copepods and have generated resources such as EST libraries, shotgun genome sequences, BAC libraries, genome maps and inbred lines that will be invaluable in assisting further efforts to provide genomics tools for copepods.

Summary

Genomics research on copepods is needed to extend our exploration and characterization of their fundamental biological traits, so that we can better understand how copepods function and interact in diverse environments. Availability of large scale genomics resources will also open doors to a wide range of systems biology type studies that view the organism as the fundamental system in which to address key questions in ecology and evolution.

The complete mitochondrial genome of Macrobrachium nipponense

The complete mitochondrial (mt) genome sequence plays an important role in the accurate determination of phylogenetic relationships among metazoans. Herein, we determined the complete mt genome sequence, structure and organization of Macrobrachium nipponense (M. nipponense) (GenBank ID: NC_015073.1) and compared it to that of Macrobrachium lanchesteri (M. lanchesteri) and Macrobrachium rosenbergii (M. rosenbergii). The 15,806 base pair (bp) M. nipponense mt genome, which is comprised of 37 genes, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs) and 2 ribosomal RNAs (rRNAs), is slightly larger than that of M. lanchesteri (15,694bp, GenBank ID: NC_012217.1) and M. rosenbergii (15,772bp, GenBank ID: NC_006880.1). The M. nipponense genome contains a high AT content (66.0%), which is a common feature among metazoan mt genomes. Compared with M. lanchesteri and M. rosenbergii, we found a peculiar non-coding region of 950bp with a microsatellite-like (TA)(6) element and many hairpin structures. The 13 PCGs are comprised of a total of 3707 codons, excluding incomplete termination codons, and the most frequently used amino acid is Leu (16.0%). The predicted start codons in the M. nipponense mt genome include ATG, ATC and ATA. Seven PCGs use TAA as a stop codon, whereas two use TAG, three use T and only one uses TA. Twenty-three of the genes are encoded on the L strand, and ND1, ND4, ND5, ND4L, 12S rRNA, 16S rRNA, tRNA(His), tRNA(Pro), tRNA(Phe), tRNA(Val), tRNA(Gln), tRNA(Cys), tRNA(Tyr) and a tRNA(Leu) are encoded on the H strand. The two rRNAs of M. nipponense and M. rosenbergii are encoded on the H strand, whereas the M. lanchesteri rRNAs are encoded on the L stand.

Complete mitochondrial genome sequence of the spinyhead croaker Collichthys lucidus (Perciformes, Sciaenidae) with phylogenetic considerations

The complete mitochondrial genome of the spinyhead croaker Collichthys lucidus was determined using long-PCR and primer walking methods. It is a circular molecule of 16,451 bp in length with a standard set of 22 tRNAs, 2 rRNAs, 13 protein-coding genes as well as a non-coding control region in the same order as those of the other bony fishes. C. lucidus mitogenome exhibited a clear strand-specific bias in nucleotide composition, as evidence by a GC- skew of the H-strand of -0.319. The strand-specific bias was also reflected in the codon usage of genes oriented in opposite directions. All tRNA genes except for tRNA( Ser(AGY)) harbored the typical cloverleaf secondary structures and possessed anticodons that matched the vertebrate mitochondrial genetic code. Contrary to the typical structure of control region consistig of TAS, central, and CSB domains, there were no central conserved blocks available in C. lucidus mitogenome. Despite extensive studies based on both morphology and molecules, phylogenetic position of C. lucidus with Sciaenidae is still controversial. Our phylogenetic results provided more evidence to support previous morphological studies and consistently placed C. lucidus as a sister taxon to Collichthys niveatus, with both of these taxa forming the monophyletic Collichthys.

Distinctive mitochondrial genome of Calanoid copepod Calanus sinicus with multiple large non-coding regions and reshuffled gene order: useful molecular markers for phylogenetic and population studies

Background

Copepods are highly diverse and abundant, resulting in extensive ecological radiation in marine ecosystems. Calanus sinicus dominates continental shelf waters in the northwest Pacific Ocean and plays an important role in the local ecosystem by linking primary production to higher trophic levels. A lack of effective molecular markers has hindered phylogenetic and population genetic studies concerning copepods. As they are genome-level informative, mitochondrial DNA sequences can be used as markers for population genetic studies and phylogenetic studies.

Results

The mitochondrial genome of C. sinicus is distinct from other arthropods owing to the concurrence of multiple non-coding regions and a reshuffled gene arrangement. Further particularities in the mitogenome of C. sinicus include low A + T-content, symmetrical nucleotide composition between strands, abbreviated stop codons for several PCGs and extended lengths of the genes atp6 and atp8 relative to other copepods. The monophyletic Copepoda should be placed within the Vericrustacea. The close affinity between Cyclopoida and Poecilostomatoida suggests reassigning the latter as subordinate to the former. Monophyly of Maxillopoda is rejected. Within the alignment of 11 C. sinicus mitogenomes, there are 397 variable sites harbouring three 'hotspot' variable sites and three microsatellite loci.

Conclusion

The occurrence of the circular subgenomic fragment during laboratory assays suggests that special caution should be taken when sequencing mitogenomes using long PCR. Such a phenomenon may provide additional evidence of mitochondrial DNA recombination, which appears to have been a prerequisite for shaping the present mitochondrial profile of C. sinicus during its evolution. The lack of synapomorphic gene arrangements among copepods has cast doubt on the utility of gene order as a useful molecular marker for deep phylogenetic analysis. However, mitochondrial genomic sequences have been valuable markers for resolving phylogenetic issues concerning copepods. The variable site maps of C. sinicus mitogenomes provide a solid foundation for population genetic studies.

Molecular characterization and expression of vitellogenin (Vg) genes from the cyclopoid copepod, Paracyclopina nana exposed to heavy metals

Induction of vitellogenin (Vg) has been used as a biomarker of exposure to heavy metals and endocrine-disruption chemicals (EDCs) in aquatic organisms. Here, we identified the full-length Vg1 and Vg2 sequence from the brackish water copepod, Paracyclopina nana. Vg1 gene contained 5718bp of the open reading frame (ORF) that encoded the putative protein of 1905 amino acids residues, while Vg2 gene consisted of 5442bp of ORF, encoding the putative protein of 1813 amino acids residues. P. nana Vgs showed highly conserved domains in the N-terminal region. The phylogenetic analysis revealed that P. nana Vgs are distinct from other arthropods, such as insects and decapods, as it formed a clade with other copepods, Tigriopus japonicus and salmon louse (Lepeophtheirus salmonis). The expression of Vg transcripts was detectable after the copepodid stages 4-5. Female copepods expressed over 83 times and 223 times more Vg1 and Vg2 transcripts, respectively, than males. When copepods were exposed to heavy metals (0.1mg/L Cd, 0.4mg/L Cu, and 2mg/L AsIII) for 24, 48, 72, and 96h, P. nana Vg transcripts were highly induced in a time-dependent manner. Interestingly, Vg2 gene was more susceptible than Vg1 to trace heavy metal exposure. This finding indicates that P. nana Vgs provide a potential indicator for assessing the toxic effect of heavy metals. In addition, we suggest P. nana as a potential model species for risk-assessment to environmental pollutants in brackish water.

Mitochondrial genome of Thais clavigera (Mollusca: Gastropoda): affirmation of the conserved, ancestral gene pattern within the mollusks

Class Gastropoda includes a large number of described species, many with extensively rearranged mitochondrial genomes. We sequenced the mitogenome of the rock shell, Thais clavigera (Gastropoda: Muricidae), an intertidal snail, using long PCR with primers designed on the basis of expressed sequence tags. The mitogenome of T. clavigera consists of 2 rRNAs, 22 tRNAs, and 13 protein-coding genes, but no control region. Structural comparisons revealed that the order Sorbeoconcha, including T. clavigera, have nearly identical mitochondrial gene patterns. However, they have an inversion between a tRNA(Phe)-tRNA(Glu) cluster that comprises 21 genes, but most of the remaining structure is similar to the putative mollusk ground pattern. These findings will provide a better insight into mitochondrial gene rearrangement over the course of gastropod evolution.

The complete mitogenome of the hydrothermal vent crab Xenograpsus testudinatus (Decapoda, Brachyura) and comparison with brachyuran crabs

In this study, we analyzed the complete mitochondrial (mt) genome of a hydrothermal vent crab Xenograpsus testudinatus (Decapoda: Brachyura) obtained from the hydrothermal vents off Kueishantao Island, Taiwan, which extend from the deep sea Okinawa Trench. The mitogenome of X. testudinatus was 15,796 bp in length and contained the same 37 genes (e.g. 2 rRNAs, 22 tRNAs, and 13 PCGs) found in other metazoan mitogenomes. Analysis of the structural mt gene order in X. testudinatus revealed that the 13 PCGs, excluding a translocation of ND6-Cyt b cluster, were similarly ordered when compared to the pancrustacean ground pattern; however the tRNAs were severely rearranged. Phylogenetic analysis of decapod mitogenomes showed that the molecular taxonomy of the vent crab was in accordance with its morphological systematics. Together, these findings suggest that the vent crab studied here has little mitochondrial genetic variation when compared with morphologically defined conspecifics from other marine habitats.