Acessibilidade / Reportar erro

First molecular evidence for two new associate copepods of genus Clausidium Kossmann, 1874 (Copepoda: Cyclopoida: Clausidiidae) from the Persian Gulf and Gulf of Oman

Abstract

Clausidium Kossman, 1874 is a genus of copepods that is found in subtropical to temperate coastal areas. All species of the genus occur on the bodies of mud shrimp of the families Callianassidae and Upogebiidae. Based on morphological data from light scanning and confocal laser scanning microscopy, there are four species of Clausidium copepod in Iran. In this study we address Clausidium iranensis Sepahvand, Kihara and Boxshall, 2019 and Clausidium persiaensis Sepahvand and Kihara, 2017 that were reported on the body of the burrowing shrimps Neocallichirus jousseaumei (Nobili, 1904) and Callianidea typa Milne Edwards, 1837, respectively. We undertook analyses of mitochondrial DNA gene sequences (CO1) to evaluate taxonomic status and taxonomic relationships of the Clausidium species. The result demonstrates that two major clades, with strong support, can be identified within the Clausidium copepods in the southern waters of Iran, representing distinct taxonomic entities at the species rank. Our data indicate that CO1 can be a powerful tool for species identification and delimitation. In the case of Clausidium copepods, the general utility of CO1 for taxonomic relationship inferences within a genus or a family is still under investigated. Our study adds the first genetic data from these copepods from the Persian Gulf and Gulf of Oman.

Keywords
CO1; Copepoda; DNA barcoding; north-west Indian Ocean

INTRODUCTION

Biological research, including biodiversity and ocean monitoring, require accurate species records and continuous observation (DeSalle and Amato, 2004DeSalle, R. and Amato, G. 2004. The expansion of conservation genetics. Nature Reviews Genetics, 5: 702-712. ). Traditionally, the biodiversity of marine habitats has been examined using morphological identification, needing specialized taxonomic expertise. Furthermore, morphological approaches lack the resolution to identify cryptic taxa or deal with morphological plasticity (Hebert et al., 2004Hebert, P.D.; Penton, E.H.; Burns, J.M.; Janzen, D.H. and Hallwachs, W. 2004. Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator. Proceedings of the National Academy of Sciences of the United States of America, 101: 14812-14817.; Pfenninger et al., 2006Pfenninger, M.; Cordellier, M. and Streit, B. 2006. Comparing the efficacy of morphologic and DNA-based taxonomy in the freshwater gastropod genus Radix (Basommatophora, Pulmonata). BMC Evolutionary Biology, 6: 1-14.). Integrating data from multiple lines of evidence will support a more reliable species delineation as the baseline for subsequent biodiversity evaluations (Dayrat, 2005Dayrat, B. 2005. Towards integrative taxonomy. Biological Journal of the Linnean Society, 85: 407-415.). Molecular approaches have been broadly used to provide a faster and more precise species identification for many taxa (Tautz et al., 2003Tautz, D.; Arctander, P.; Minelli, A.; Thomas, R.H. and Vogler, A.P. 2003. A plea for DNA taxonomy. Trends in Ecology and Evolution, 18: 70-74.; Vogler and Monaghan 2007Vogler, A.P and Monaghan, M.T. 2007. Recent advances in DNA taxonomy. Journal of Zoological Systematics and Evolutionary Research, 45: 1-10.; Pereira et al., 2008Pereira, F.; Carneiro, J. and Amorim, A. 2008. Identification of species with DNA-based technology: current progress and challenges. Recent Patents on DNA and Gene Sequences, 2: 187-199.).

In spite of the recent progress in other diverse genetic and genomic indicators, the barcode region of the mitochondrial cytochrome c oxidase subunit I (CO1) gene remains a useful, and in some cases unique, diagnostic character for species-level identification of copepods (Blanco-Bercial et al., 2014Blanco-Bercial, L.; Cornils, A.; Copley N. and Bucklin A. 2014. DNA barcoding of marine copepods: assessment of analytical approaches to species identification. PLoS Currents Tree of Life, 1. doi: 10.1371/currents.tol.cdf8b74881f87e3b01d56b43791626d2.
https://doi.org/10.1371/currents.tol.cdf...
).

Ghost shrimps or burrowing shrimps, representatives of the infraorders Axiidea and Gebiidea, are among the most common benthic macro-invertebrates in coastal regions of the Persian Gulf and Gulf of Oman (Sepahvand et al., 2013Sepahvand, V.; Sari, A.; Salehi, H.; Nabavi, S.M.B. and Ghorbanzadeh, S.G. 2013. Littoral mud shrimps (Decapoda: Gebiidea & Axiidea) of the Persian Gulf and Gulf of Oman, Iran. Journal of the Marine Biological Association of the United Kingdom, 93: 999-100.). These cryptic shrimps are adapted to a burrowing lifestyle and their burrows can also be occupied by a variety of organisms, including copepods (Dworschak et al., 2012Dworschak, P.C.; Felder D.F.; and Tudge, C.C. 2012. Infraorders Axiidea de Saint Laurent 1979 and Gebiidea de Saint Laurent 1979 (formerly known collectively as Thalassinidea). p. 109-219. In: F.R. Schram and J.C. von Vaupel Klein (eds), Treatise on Zoology - Anatomy, Taxonomy, Biology. The Crustacea. Complementary to the volumes translated from the French of the Traite´ de Zoologie [founded by P.-P. Grasse´]. 9 Part B. Eucarida: Decapoda: Astacidea p.p. (Enoplometopoidea, Nephropoidea), Glypheidea, Axiidea, Gebiidea, and Anomura. Brill, Leiden. ). Callianidea typa Milne Edwards, 1837Milne Edwards, H. 1834-1840. Histoire Naturelle des Crustacés, Comprenant l´Anatomie, la Physiologie et la Classification de ces Animaux. Encyclopédique Roret, Paris. Vol. III, 638 pp., plates 1-42. and Neocallichirus jousseaumei (Nobili, 1904Nobili, G. 1904. Diagnoses preliminaires de vingt-huit especes nouvelles de Stomatopodes et Decapodes Macroures de la mer Rouge. Bulletin du Muséum National d’Histoire Naturelle, Paris, 10: 230-238.) are two widely distributed burrowing shrimps in the southern waters of Iran that are reported as hosts for Clausidium Kossman, 1874 copepods (Sepahvand et al., 2017Sepahvand, V.; Rastegar-Pouyani, N.; Kihara, T.C. and Momtazi, F. 2017. A new species of Clausidium Kossmann, 1874 (Crustacea, Copepoda, Cyclopoida, Clausidiidae) associated with ghost shrimps from Iran. Nauplius, 25: 1-16. ; 2019Sepahvand, V.; Kihara, T.C. and Boxshall, G.A. 2019. A new species of Clausidium Kossmann, 1874 (Copepoda: Cyclopoida) associated with ghost shrimps from the Persian Gulf, including female-male interlocking mechanisms and remarks on host specificity. Systematic Parasitology, 96: 171-189.). Copepods of the genus Clausidium are external associates of burrowing decapods of the families Callianassidae and Upogebiidae (Kihara and Rocha, 2013Kihara, T.C. and Rocha, C.E.F. 2013. First record of Clausidium (Copepoda, Clausidiidae) from Brazil: a new species associated with ghost shrimps Neocallichirus grandimana (Gibbes, 1850) (Decapoda, Callianassidae). ZooKeys, 335: 47-67.).

Most clausidiid copepods live in loose association with marine invertebrate hosts (Huys and Boxshall, 1991Huys, R. and Boxshall, G.A. 1991. Copepod Evolution. London, Ray Society, 468p.), and species of Clausidium are recorded exclusively living in association with ghost shrimps (Boxshall and Halsey, 2004Boxshall, G.A and Halsey, S.H. 2004. An Introduction to Copepod Diversity. London, The Ray Society, 966p.). Although it has been suggested that members of Clausidium are parasitic on their hosts (Wilson, 1935Wilson, C.B. 1935. Parasitic copepods from the Pacific Coast. American Midland Naturalist, 16: 776-797.; Pillai, 1959Pillai, N.K. 1959. On two new species of Clausidium (Copepoda: Cyclopoida) parasitic on the shrimp Callianassa. Journal of the Marine Biological Association of India, 1: 57-65.), this relationship has yet to be conclusively demonstrated (Hayes, 1976Hayes, H.J. 1976. Biology of Clausidium dissimile an epizoic copepod. Pensacola, University of West Florida, Master thesis, 103p.). There is very scarce available documentation on the biology of these copepods, or their interactions with their host, or with the environment. Although Clausidium species are rarely recorded because of the cryptic lifestyle of their hosts, a total of 18 species of Clausidium have been described to date and it is hypothesized that each species shows a preference for a specific host Sepahvand et al. (2019Sepahvand, V.; Kihara, T.C. and Boxshall, G.A. 2019. A new species of Clausidium Kossmann, 1874 (Copepoda: Cyclopoida) associated with ghost shrimps from the Persian Gulf, including female-male interlocking mechanisms and remarks on host specificity. Systematic Parasitology, 96: 171-189.).

The only available molecular study on Clausidium was carried out by Huys et al. (2012 Huys, R.; Fatih, F.; Ohtsuka, S. and Llewellyn-Hughes, J. 2012. Evolution of the bomolochiform superfamily complex (Copepoda: Cyclopoida): New insights from ssrDNA and morphology, and origin of umazuracolids from polychaete-infesting ancestors rejected. International Journal for Parasitology, 42: 71-92.). In that study, the 18S rRNA gene, in combination with morphological features, were investigated in an integrative approach to study a Clausidiiform complex (Huys et al., 2012 Huys, R.; Fatih, F.; Ohtsuka, S. and Llewellyn-Hughes, J. 2012. Evolution of the bomolochiform superfamily complex (Copepoda: Cyclopoida): New insights from ssrDNA and morphology, and origin of umazuracolids from polychaete-infesting ancestors rejected. International Journal for Parasitology, 42: 71-92.). Here, we use CO1 gene sequences to delineate species boundaries and find cryptic diversity within Clausidium in the Persian Gulf and the Gulf of Oman.

MATERIALS AND METHODS

Sample collection

Copepods specimens were obtained from the body of N. jousseaumei and C. typa (Fig. 1) from the Iranian coast of the Persian Gulf and Gulf of Oman (Fig. 2). Specimens were preserved in 96 % ethanol at the sampling site and after 24 h, these specimens were transferred to fresh 96 % ethanol at the Iranian National Institute for Oceanography and Atmospheric Science (INIOAS). Clausidium copepods were later sorted using a microscope-mounted camera at the German Center for Marine Biodiversity Research (DZMB) in Wilhelmshaven, Germany. Photographs of 96 selected specimens in ethanol were obtained using a camera for further morphological identification. Individuals were preserved in 96% ethanol and stored at -20 °C, as morphological vouchers for future reference. Specimens were identified to species level using diagnostic morphological characters based on the identification key of Kihara and Rocha (2013Kihara, T.C. and Rocha, C.E.F. 2013. First record of Clausidium (Copepoda, Clausidiidae) from Brazil: a new species associated with ghost shrimps Neocallichirus grandimana (Gibbes, 1850) (Decapoda, Callianassidae). ZooKeys, 335: 47-67.).

Figure 1.
A, Callianidea typa, lateral view of cephalothorax with copepods; B, Clausidium persiaensis, habitus, ventral view, female; C, C. persiaensis, confocal laser scanning microscopy maximum projections; D, Neocallichirus jousseaumei, habitus, dorsal view; E, Clausidium persiaensis, habitus, ventral view; F, C. iranensis, confocal laser scanning microscopy maximum projections, dorsal view. Scale bars: A, 1 mm; B, 0.5 mm; C, 100 µm; D, 1 cm; E, 1 mm; F, 100 µm.

Figure 2.
Map of sampling stations in the Persian Gulf and the Gulf of Oman.

DNA extraction and gene amplification

DNA extractions from 63 selected morphologically identified specimens stored at -20 °C before use, were carried out using the 30-40 μl Chelex (InstaGene Matrix, Bio−Rad) protocol. The primers LCO1490: 5' -GGTCAACAAATCATAAAGATATTGG-3' and HCO2198: 5'-TAAACTTCAGGGTGACCAAAAAATCA-3' (Folmer et al.,1994Folmer, O.; Black M.; Hoeh, W.; Lutz, R. and Vrijenhoek, R. 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Molecular marine biology and biotechnology, 3: 294-299.) were used to amplify the CO1 gene. PCR cycles consisted of an initial denaturation at 95 °C for 5 min, followed by a denaturation at 95 °C for 30 s, annealing at 45 °C for 1 min, and extension at 72 °C for 1 min, for 60 cycles and a final elongation 72 °C for 7 min. The number of cycles was decreased if the concentration of PCR product reached an optimum. The PCR was performed using IllustraPuReTaq Ready−To−Go PCR Beads (GE Healthcare) in 25 μL volume containing 22 μL H2O, 0.5 μL of each primer (10 pmol/μL) and 2 μL of DNA templates. All PCR products were checked by electrophoresis on a 1 % agarose/TBE gel containing 1 % GelRed. PCR product purifications and sequencing was carried out by Macrogen (Amsterdam, Netherlands). All the COI sequences are deposited in GenBank and accession numbers are provided in Tab. 1.

Table 1.
Sequence accession numbers in NCBI, host of Clausidium copepods, and locality of sample collection.

Genetic analyses

Sequences of CO1 were aligned using Clustal W, as implemented in the Bioedit software version 7.0.5.3 (Hall, 1999Hall, T.A. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium, Series 41: 95-98.). The DNA substitution model of TIM1+I+G (AIC = 6054.89, −lnL = 2967.44, k = 60, p-inv = 0.43, gamma shape = 0.73) was estimated using the Akaike information criterion run in jModelTest, version 0.1.1 (Posada, 2008Posada, D. 2008. jModelTest: phylogenetic model averaging. Molecular Biology and Evolution, 25: 1253-1256.). Taxonomic relationships were reconstructed using Maximum Parsimony method in MEGAX (Kumar et al., 2018Kumar, S.; Stecher, G.; Li, M.; Knyaz, C. and Tamura, K. 2018. MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Molecular Biology and Evolution, 35: 1547-1549.), the Maximum Likelihood method implemented in PHYML version 2.4.4 (Guindon and Gascuel, 2003Guindon, S. and Gascuel, O. 2003. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic Biology, 52: 696-704.) and a Bayesian Inference tree in MrBayes version 3.1.2 (Huelsenbeck and Ronquist, 2001Huelsenbeck, J.P and Ronquist, F. 2001. MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics, 17: 754-755.; Ronquist and Huelsenbeck, 2003Ronquist, F. and Huelsenbeck, J.P. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics, 19: 1572-1574.). Bayesian Inference was performed with two simultaneous runs and four search chains within each run (three heated chains and one cold chain) for 10,000,000 generations, sampling trees every 1000 generations using the Markov chain Monte Carlo method. Reliability of nodes was assessed using 1000 bootstrap replications for all methods (Felsenstein, 1985Felsenstein, J. 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution, 39: 783-791.). Genetic distances were calculated to quantify sequence divergence between species using Kimura’s (Kimura, 1980Kimura, M. 1980. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences.Journal of Molecular Evolution, 16: 111-120.) two-parameter (K2P) model by MEGAX (Kumar et al., 2018Kumar, S.; Stecher, G.; Li, M.; Knyaz, C. and Tamura, K. 2018. MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Molecular Biology and Evolution, 35: 1547-1549.). Genetic diversity was measured for each species based on haplotype diversity (Hd) and nucleotide diversity (π). Values for the number of polymorphic sites, parsimony informative sites, haplotype frequencies, and the average number of nucleotide differences between sequences were estimated. These genetic diversity values were computed using the software DnaSP v5.0 (Librado and Rozas, 2009Librado, P. and Rozas J. 2009. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics, 25: 1451-1452.). Taxonomic relationship networks were calculated using the software PopArt (Leigh and Bryant, 2015Leigh, J.W. and Bryant D. 2015. POPART: full-feature software for haplotype network construction. Methods in Ecology and Evolution, 6: 1110-1116.; http://popart.otago.ac.nz).

RESULTS

Maximum Likelihood (ML), Maximum Parsimony (MP) and Bayesian Inference (BI) based on CO1 sequences obtained trees with similar topologies (with minor changes), each confirming the taxonomic status of two recently described morphospecies of C. persiaensis and C. iranensis as distinct lineages. The species were represented by monophyletic clades on all the taxonomic relationship analyses. Moreover, in all trees, Clausidium grouped into two major clades: one clade included sequences from Qeshm Island in the Persian Gulf as C. iranensis, while the other included samples from the Oman Gulf as C. persiaensis. Monophyly of these clades was strongly supported by BI posterior probability, MP, and ML bootstrap values (1/100/100 respectively for all the Clausidium clades) (Fig. 3).

Figure 3.
Maximum Likelihood (ML) gene tree based on partial CO1 gene sequences of the Clausidium species from the Persian Gulf and the Oman Gulf. The average branch lengths are proportional to the number of substitutions per site. Numbers above the nodes indicate Bayesian posterior probabilities and numbers below the nodes represent bootstrap support values for MP/ML (1000 replicates).

The rate of variation between sites was modeled with a gamma distribution (shape parameter = 1). The analysis used 27 nucleotide sequences. Codon positions included were 1st+2nd+3rd+Noncoding and all positions containing gaps and missing data were eliminated. There was a total of 645 positions in the final dataset. The mean genetic divergence (%) in the CO1 gene sequences between the species of the genera Hemicyclops Boeck,1872, Clausidium, and Pseudobradya Sars G.O., 1904 are shown in Tab. 2. The Kimura’s (1980Kimura, M. 1980. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences.Journal of Molecular Evolution, 16: 111-120.) two-parameter (K2P) genetic distance between C. persiaensis and C. iranensis was 24.6 % and between these two species and outgroups ranged from 26 to 43.5 %. The mean genetic p- distance between the two species was 18.2 %. The mean genetic divergence (%) within C. persiaensis and C. iranensis were 0.44 % and 0.13 %, respectively. K2P interspecies genetic distances within genus Clausidium (between C. persiaensis and C. iranensis) were 99.2 times higher than the mean intraspecific genetic divergence.

Table 2.
The percentage of mean genetic divergence in the CO1 gene between the species of the genera Hemicyclops, Clausidium, and Pseudobradya. Numbers in upper right (Bold) are P-distances and in lower left are Kimura’s two-parameter (K2P) distances.

Based on the 648 bp of CO1 examined in C. persiaensis (n = 15), 12 sites were variable (polymorphic) and 8 sites were parsimony informative, resulting in the identification of 9 haplotypes (Fig. 4). Haplotype diversity (Hd) and nucleotide diversity (π) were 0.914 and 0.0044 for this species, respectively. Furthermore, based on the 648 bp of COI examined in C. iranensis (n = 7), 3 sites were variable (polymorphic) and 2 sites were parsimony informative, resulting in the identification of 3 haplotypes. Haplotype diversity (Hd) and nucleotide diversity (π) were 0.714 and 0.0013 for C. iranensis, respectively.

Figure 4.
95 % minimum spanning haplotype network of COI haplotypes of Clausidium persiaensis from the Persian Gulf, Qeshm Island. The size of the circle is proportional to the frequency of that haplotype.

DISCUSSION

This study is the first to assess the species delimitation of burrowing shrimp-associated copepods via molecular markers in the Persian Gulf and Gulf of Oman. The results extracted by molecular analysis in the present study absolutely confirmed morphospecies: C. persiaensis and C. iranensis.

The morphological species concept is most commonly applied in Clausidium copepod taxonomy. The molecular sequence data available for Clausidium copepods is still very scarce and the relationships within the genera of the Clausidiidae remain elusive, because of the wide host range they utilize and the different morphologies of this group. Copepods of the family Clausidiidae represent an early offshoot of the Poecilostome lineage within the order Cyclopoida Burmeister 1834 (see Khodami et al., 2017Khodami, S.; McArthur, J.; Blanco-Bercial, L. and Martinez Arbizu, P. 2017. Molecular phylogeny and revision of copepod orders (Crustacea: Copepoda). Scientific Reports, 7: 9164.). Most clausidiids live in loose association with marine invertebrate hosts (Huys and Boxshall, 1991Huys, R. and Boxshall, G.A. 1991. Copepod Evolution. London, Ray Society, 468p.) and species of Clausidium are recorded exclusively living in association with ghost or burrowing shrimps (Boxshall and Halsey, 2004Boxshall, G.A and Halsey, S.H. 2004. An Introduction to Copepod Diversity. London, The Ray Society, 966p.).

Clausidium persiaensis and C. iranensis can be readily identified based on their morphological characteristics (Sepahvand et al., 2017Sepahvand, V.; Rastegar-Pouyani, N.; Kihara, T.C. and Momtazi, F. 2017. A new species of Clausidium Kossmann, 1874 (Crustacea, Copepoda, Cyclopoida, Clausidiidae) associated with ghost shrimps from Iran. Nauplius, 25: 1-16. ; 2019Sepahvand, V.; Kihara, T.C. and Boxshall, G.A. 2019. A new species of Clausidium Kossmann, 1874 (Copepoda: Cyclopoida) associated with ghost shrimps from the Persian Gulf, including female-male interlocking mechanisms and remarks on host specificity. Systematic Parasitology, 96: 171-189.). Clausidium iranensis shares the armature formula of swimming legs 2 to 4 with C. persiaensis but can be easily distinguished by unique characteristics of the females: the prominent spine on endopodal segment 1 of the antenna, the armature of the maxilliped, and the elongated basis of the swimming legs (Sepahvand et al., 2019Sepahvand, V.; Kihara, T.C. and Boxshall, G.A. 2019. A new species of Clausidium Kossmann, 1874 (Copepoda: Cyclopoida) associated with ghost shrimps from the Persian Gulf, including female-male interlocking mechanisms and remarks on host specificity. Systematic Parasitology, 96: 171-189.). In this study, we were able to show that they are genetically distinct lineages based on the CO1 gene. Moreover, the mean genetic distance between C. persiaensis and C. iranensis was 99.2-fold higher than the mean intraspecific genetic variation for each species. Our findings, furthermore, emphasize the effectiveness of molecular data, such as CO1 gene sequences, in species delimitation and identification for marine metazoans. However, some drawbacks of using CO1, and mtDNA in general, for species identification, include the possible co-amplification of nuclear mitochondrial pseudogenes (numts) (Song et al., 2008Song, H.; Buhay, J.E.; Whiting, M.F. and Crandall, K.A. 2008. Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. Proceedings of the National Academy of Sciences of the United States of America, 105: 13486-13491.; Hazkani-Covo, 2010Hazkani-Covo, E.; Zeller, R. M. and Martin, W. 2010. Molecular Poltergeists: Mitochondrial DNA Copies numts in Sequenced Nuclear Genomes. PLoS Genet, 6(2): e1000834.), introgression through hybridization or incomplete lineage sorting, and heteroplasmy (Hoeh et al., 1991Hoeh, W.R.; Blakley K.H. and Brown, W.M. 1991. Heteroplasmy suggests limited biparental inheritance of Mytilus mitochondrial DNA. Science, 251: 1488-1490.).

Our taxonomic relationship analyses demonstrated that, two major clades strongly supported by BI posterior probability, MP, and ML bootstrap values (1/100/100 respectively for all Clausidium clades) can be identified within Clausidium copepods in the southern waters of Iran. Each clade comprises individuals of one of the two morphospecies with substantial genetic divergence between the clades. The degree of genetic divergence between C. persiaensis and C. iranensis was 24.8 % and between these two species and outgroups ranged from 26 to 43.5 %. This large genetic divergence taxonomically supports placing the two distinct taxa at the rank of species. Here, we found congruency between morphological and molecular species delimitation of the genus Clausidium in Iran.

Despite tremendous effort, we were only able to obtain sequences from two species of four recognized species of Clausidium copepods recorded in the Persian Gulf and the Gulf of Oman. This highlights the challenges of collecting Clausidium copepods that are only associated with ghost or burrowing shrimps because they have such a cryptic lifestyle and are hard to find. Sampling these copepods is usually very time consuming and expensive. Moreover, extracting high quality DNA from the collected specimens can be daunting.

The lack of a complete DNA barcode library is the most limiting parameter for precise and trustworthy discrimination and documentation of species of copepods. In fact, DNA barcodes are currently available for only ~ 400 copepod species, including many parasitic and freshwater taxa (Blanco-Bercial et al., 2014Blanco-Bercial, L.; Cornils, A.; Copley N. and Bucklin A. 2014. DNA barcoding of marine copepods: assessment of analytical approaches to species identification. PLoS Currents Tree of Life, 1. doi: 10.1371/currents.tol.cdf8b74881f87e3b01d56b43791626d2.
https://doi.org/10.1371/currents.tol.cdf...
). In addition, extensive coverage of species diversity is especially critical for efficient resolution in large datasets using automated methods.

In conclusion, our data indicate that CO1 may be a powerful tool for species identification and delimitation. In the case of Clausidium copepods, the general utility of CO1 for taxonomic relationship inferences within a genus, or a family, is still under investigation. Hence, our study adds a few, but highly important, data points for future comparative studies.

ACKNOWLEDGEMENTS

We are very grateful to Dr. Terue Cristina Kihara and Dr. Sahar Khodami from Senckenberg am Meer, German Center for Marine Biodiversity Research, Wilhelmshaven, for help and encouragement during this study. The authors also express their gratitude to Dr. Abdolvahab Maghsoudlou from the Iranian National Institute for Oceanography and Atmospheric Science for his cooperation in this study.

REFERENCES

  • Boxshall, G.A and Halsey, S.H. 2004. An Introduction to Copepod Diversity. London, The Ray Society, 966p.
  • Blanco-Bercial, L.; Cornils, A.; Copley N. and Bucklin A. 2014. DNA barcoding of marine copepods: assessment of analytical approaches to species identification. PLoS Currents Tree of Life, 1. doi: 10.1371/currents.tol.cdf8b74881f87e3b01d56b43791626d2.
    » https://doi.org/10.1371/currents.tol.cdf8b74881f87e3b01d56b43791626d2
  • Dayrat, B. 2005. Towards integrative taxonomy. Biological Journal of the Linnean Society, 85: 407-415.
  • DeSalle, R. and Amato, G. 2004. The expansion of conservation genetics. Nature Reviews Genetics, 5: 702-712.
  • Dworschak, P.C.; Felder D.F.; and Tudge, C.C. 2012. Infraorders Axiidea de Saint Laurent 1979 and Gebiidea de Saint Laurent 1979 (formerly known collectively as Thalassinidea). p. 109-219. In: F.R. Schram and J.C. von Vaupel Klein (eds), Treatise on Zoology - Anatomy, Taxonomy, Biology. The Crustacea. Complementary to the volumes translated from the French of the Traite´ de Zoologie [founded by P.-P. Grasse´]. 9 Part B. Eucarida: Decapoda: Astacidea p.p. (Enoplometopoidea, Nephropoidea), Glypheidea, Axiidea, Gebiidea, and Anomura. Brill, Leiden.
  • Felsenstein, J. 1985. Confidence limits on phylogenies: an approach using the bootstrap. Evolution, 39: 783-791.
  • Folmer, O.; Black M.; Hoeh, W.; Lutz, R. and Vrijenhoek, R. 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Molecular marine biology and biotechnology, 3: 294-299.
  • Kihara, T.C. and Rocha, C.E.F. 2013. First record of Clausidium (Copepoda, Clausidiidae) from Brazil: a new species associated with ghost shrimps Neocallichirus grandimana (Gibbes, 1850) (Decapoda, Callianassidae). ZooKeys, 335: 47-67.
  • Khodami, S.; McArthur, J.; Blanco-Bercial, L. and Martinez Arbizu, P. 2017. Molecular phylogeny and revision of copepod orders (Crustacea: Copepoda). Scientific Reports, 7: 9164.
  • Kimura, M. 1980. A simple method for estimating evolutionary rate of base substitutions through comparative studies of nucleotide sequences.Journal of Molecular Evolution, 16: 111-120.
  • Kumar, S.; Stecher, G.; Li, M.; Knyaz, C. and Tamura, K. 2018. MEGA X: Molecular Evolutionary Genetics Analysis across computing platforms. Molecular Biology and Evolution, 35: 1547-1549.
  • Guindon, S. and Gascuel, O. 2003. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Systematic Biology, 52: 696-704.
  • Hall, T.A. 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symposium, Series 41: 95-98.
  • Hayes, H.J. 1976. Biology of Clausidium dissimile an epizoic copepod. Pensacola, University of West Florida, Master thesis, 103p.
  • Hazkani-Covo, E.; Zeller, R. M. and Martin, W. 2010. Molecular Poltergeists: Mitochondrial DNA Copies numts in Sequenced Nuclear Genomes. PLoS Genet, 6(2): e1000834.
  • Huelsenbeck, J.P and Ronquist, F. 2001. MrBayes: Bayesian inference of phylogenetic trees. Bioinformatics, 17: 754-755.
  • Huys, R. and Boxshall, G.A. 1991. Copepod Evolution. London, Ray Society, 468p.
  • Huys, R.; Fatih, F.; Ohtsuka, S. and Llewellyn-Hughes, J. 2012. Evolution of the bomolochiform superfamily complex (Copepoda: Cyclopoida): New insights from ssrDNA and morphology, and origin of umazuracolids from polychaete-infesting ancestors rejected. International Journal for Parasitology, 42: 71-92.
  • Hebert, P.D.; Penton, E.H.; Burns, J.M.; Janzen, D.H. and Hallwachs, W. 2004. Ten species in one: DNA barcoding reveals cryptic species in the neotropical skipper butterfly Astraptes fulgerator Proceedings of the National Academy of Sciences of the United States of America, 101: 14812-14817.
  • Hoeh, W.R.; Blakley K.H. and Brown, W.M. 1991. Heteroplasmy suggests limited biparental inheritance of Mytilus mitochondrial DNA. Science, 251: 1488-1490.
  • Librado, P. and Rozas J. 2009. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics, 25: 1451-1452.
  • Leigh, J.W. and Bryant D. 2015. POPART: full-feature software for haplotype network construction. Methods in Ecology and Evolution, 6: 1110-1116.
  • Pereira, F.; Carneiro, J. and Amorim, A. 2008. Identification of species with DNA-based technology: current progress and challenges. Recent Patents on DNA and Gene Sequences, 2: 187-199.
  • Pillai, N.K. 1959. On two new species of Clausidium (Copepoda: Cyclopoida) parasitic on the shrimp Callianassa Journal of the Marine Biological Association of India, 1: 57-65.
  • Milne Edwards, H. 1834-1840. Histoire Naturelle des Crustacés, Comprenant l´Anatomie, la Physiologie et la Classification de ces Animaux. Encyclopédique Roret, Paris. Vol. III, 638 pp., plates 1-42.
  • Nobili, G. 1904. Diagnoses preliminaires de vingt-huit especes nouvelles de Stomatopodes et Decapodes Macroures de la mer Rouge. Bulletin du Muséum National d’Histoire Naturelle, Paris, 10: 230-238.
  • Pfenninger, M.; Cordellier, M. and Streit, B. 2006. Comparing the efficacy of morphologic and DNA-based taxonomy in the freshwater gastropod genus Radix (Basommatophora, Pulmonata). BMC Evolutionary Biology, 6: 1-14.
  • Posada, D. 2008. jModelTest: phylogenetic model averaging. Molecular Biology and Evolution, 25: 1253-1256.
  • Ronquist, F. and Huelsenbeck, J.P. 2003. MrBayes 3: Bayesian phylogenetic inference under mixed models. Bioinformatics, 19: 1572-1574.
  • Sepahvand, V.; Sari, A.; Salehi, H.; Nabavi, S.M.B. and Ghorbanzadeh, S.G. 2013. Littoral mud shrimps (Decapoda: Gebiidea & Axiidea) of the Persian Gulf and Gulf of Oman, Iran. Journal of the Marine Biological Association of the United Kingdom, 93: 999-100.
  • Sepahvand, V.; Rastegar-Pouyani, N.; Kihara, T.C. and Momtazi, F. 2017. A new species of Clausidium Kossmann, 1874 (Crustacea, Copepoda, Cyclopoida, Clausidiidae) associated with ghost shrimps from Iran. Nauplius, 25: 1-16.
  • Sepahvand, V.; Kihara, T.C. and Boxshall, G.A. 2019. A new species of Clausidium Kossmann, 1874 (Copepoda: Cyclopoida) associated with ghost shrimps from the Persian Gulf, including female-male interlocking mechanisms and remarks on host specificity. Systematic Parasitology, 96: 171-189.
  • Song, H.; Buhay, J.E.; Whiting, M.F. and Crandall, K.A. 2008. Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified. Proceedings of the National Academy of Sciences of the United States of America, 105: 13486-13491.
  • Tautz, D.; Arctander, P.; Minelli, A.; Thomas, R.H. and Vogler, A.P. 2003. A plea for DNA taxonomy. Trends in Ecology and Evolution, 18: 70-74.
  • Vogler, A.P and Monaghan, M.T. 2007. Recent advances in DNA taxonomy. Journal of Zoological Systematics and Evolutionary Research, 45: 1-10.
  • Wilson, C.B. 1935. Parasitic copepods from the Pacific Coast. American Midland Naturalist, 16: 776-797.

Publication Dates

  • Publication in this collection
    16 Apr 2021
  • Date of issue
    2021

History

  • Received
    16 Mar 2020
  • Accepted
    05 Feb 2021
Sociedade Brasileira de Carcinologia Instituto de Biociências, UNESP, Campus Botucatu, Rua Professor Doutor Antônio Celso Wagner Zanin, 250 , Botucatu, SP, 18618-689 - Botucatu - SP - Brazil
E-mail: editor.nauplius@gmail.com