- Short report
- Open Access
First insights on the mitochondrial genetic variability of Lightiella magdalenina (Crustacea), the sole Mediterranean cephalocarid species
Journal of Biological Research-Thessaloniki volume 21, Article number: 5 (2014)
Here we report the first insight into the mitochondrial (Cytochrome c Oxidase subunit I - COI and Cytochrome b - Cyt b) genetic variation of the only Mediterranean cephalocarid Lightiella magdalenina.
COI sequences provide a scenario of low intraspecific variability, while significant genetic divergence occurs between L. magdalenina and L. incisa. Interestingly, Cyt b sequences reveal a higher degree of intraspecific variability, with no shared haplotypes between the sites considered.
In the future, COI and Cyt b molecular markers could be used as valuable tools to shed new light into the extant species within the genus Lightiella thus providing molecular support to the taxonomical identifications carried out on a morphological basis.
Cephalocarids are a class of small benthic crustaceans distributed from the intertidal zone to a depth of approximately 1550 m . They have been recently considered to be phylogenetically related to the class Remipedia (Crustacea), using both nuclear  and mitochondrial  markers. Currently, only one order (Brachypoda) has been described with one family (Hutchinsoniellidae) containing twelve species belonging to five genera (Hutchinsoniella, Lightiella, Sandersiella, Chiltoniella, Hampsonellus) .
The genus Lightiella harbours the greatest number of species (five) and the largest distribution (for a review see ), with L. serendipita Jones, 1961 reported from San Francisco Bay, California, L. incisa Gooding, 1963 from the Gulf of Mexico and Caribbean Sea, L. monniotae Cals & Delamare Deboutteville, 1970 from New Caledonia, L. floridana McLaughlin, 1976 from the western coast of Florida, and L. magdalenina Carcupino, Floris, Addis, Curini-Galletti, Castelli, 2006 from Sardinia. However, each species of this genus presents a very restricted distribution, with the exception of L. incisa, which is reported to populate the entire Caribbean Sea and the Gulf of Mexico.
The five species belonging to the genus Lightiella show limited morphological differentiation [4, 5]. Given the limited number of specimens used for taxonomical classification, their detailed external morphology is often incomplete e.g. . In light of these considerations, Olesen et al.  advise that the genus needs an urgent overhaul. In particular, these authors suggest that the species richness within the whole Cephalocarida class may have been overestimated. A molecular taxonomic approach could address this problem although molecular data are so far available for only three cephalocarid species: Hutchinsoniella macracantha, L. magdalenina, and L. incisa. Mitochondrial and nuclear genes sequences (see [7–14] for details) are available for H. macracantha. Sequences of mitochondrial genes are known for L. magdalenina (see  for details), while mitochondrial and nuclear gene sequences are known for L. incisa (see  for details).
The species Lightiella magdalenina is the only cephalocarid reported in the Mediterranean and is only found so far from its type locality in La Maddalena Archipelago (Sardinia, Italy) (see  for details). Cladistic analysis based on external morphology  suggest that L. magdalenina is more closely related to the Caribbean species L. incisa and L. floridana than to the Pacific Ocean species L. serendipita (California, Pacific Ocean).
In this study, we provide the first preliminary insight into the genetic variation of the isolated Mediterranean cephalocarid L. magdalenina using sequences of the mitochondrial DNA Cytochrome c Oxidase subunit I (COI) and Cytochrome b (Cyt b) coding genes from the two close Mediterranean sites where this species is found.
We also provide new specific primers for the COI gene. Results obtained were compared with data from other species of cephalocarids and remipeds.
Eighteen specimens (10 adults and 8 larvae) of Lightiella magdalenina were collected at a depth of 15–25 m from muddy sand of the southern shore of La Maddalena Archipelago (Sardinia, Italy) (Figure 1). Sampling was carried out, using SCUBA, during the 2007–2011 spring/summer seasons. We analysed specimens from the type locality (Punta San Giorgio) and from an adjacent site (Isolotto Roma) (approximately 600 m farther), both characterised by bottoms with Posidonia oceanica and Caulerpa racemosa (see Figure 1 for details). No individuals of L. magdalenina were found when sampling campaigns were carried out in other eight sites from five western Mediterranean localities in northern Sardinia and southern France (see Figure 1 for details).
DNA was extracted from whole specimens by means of the QIAGEN® DNeasy Tissue kit (Hilden, Germany). DNA concentration was estimated through fluorometric measurements (average value per sample: 10 ng/μL).
Since COI universal primers for marine invertebrates  did not result in good quality sequences, the COI region was amplified using specific primers designed in the present study (L: 5′-gttactctggggggattcgga-3′ and H: 5′-ggccaaaatagaagaaaccccagcta-3′). Universal primers reported by Boore and Brown  (L: 5′-ggwtaygtwytwccwtgrggwcarat-3′ and H: 5′-gcrtawgcraawarraartaycaytcwgg-3′) were used for the Cyt b gene. Attempts to design new specific primers for this region, without critically affecting the quality of the amplified fragment, failed. Each 25 μL PCR mixture contained approximately 50 ng of total genomic DNA, 0.4 μM of each primer, and 2.5 U of EuroTaq DNA Polymerase (Euroclone®, Italy) in a reaction mix prepared according to the manufacturer’s instructions. The MgCl2 concentration was set at 2.5 mM.
PCR amplification was performed with an initial denaturation step of 2 min at 94°C, followed by 35 cycles of 1 min at 94°C, 1 min and 30 sec at 49°C, and 1 min and 30 sec at 72°C. The final elongation step was for 5 min at 72°C with a final cooling at 4°C. In all experiments, negative controls and replicates were included. Electrophoretic runs were carried out at 4 V/cm for 20 min in 2% agarose gels prepared using 0.5 × TBE buffer. The gels were then stained with ethidium bromide (10 mg/mL).
PCR products were purified with ExoSAP-IT (USB Corporation®, USA) and sequenced using an external sequencing core service (Macrogen Inc.®, Europe). Sequences were aligned using Clustal W , implemented in the BioEdit 188.8.131.52 software package , and deposited in GenBank [COI:JX013538-JX013555; Cyt b:JX013556-JX013565].
In order to characterise genetic variation among individuals in both populations, the number of polymorphic sites (S), of haplotypes (H), the estimates of haplotype diversity (h), and nucleotide diversity (π), were computed using the software DnaSP 5.10 .
Genetic relationships among haplotypes were investigated by constructing a median-joining network using the software Network 4.6 [http://www.fluxus-engineering.com/]. In the absence of random recurrent polymorphisms across the whole dataset, an identical weight (10) was given to all mutations to preserve the original genetic information. We ran the Maximum Parsimony (MP) calculation post-processing option to delete all superfluous median vectors and links that were not contained in the shortest trees of the network. For the COI gene, the analysis was performed by adding also the sequences of two other cephalocarids, Lightiella incisa [GenBank:GQ328968] and Hutchinsoniella macracantha [GenBank:AF370852]. In addition, two sequences, representative of the two families of the class Remipedia, Pleomothra apletocheles (Godzilliidae) [GenBank:GU067682] and Speleonectes tulumensis (Speleonectidae) [GenBank:NC_005938], were included in the analysis as outgroups. For the Cyt b gene, the analysis was carried out including the cephalocarid Hutchinsoniella macracantha [GenBank:AY456189] plus the remipede Speleonectes tulumensis [GenBank:NC_005938] as outgroup.
Eighteen individuals were sequenced and a 216 base pair (bp) alignment was obtained. A previously published sequence [GenBank:EU530536]  belonging to an individual from Punta San Giorgio (see Table 1 for details) was added to the dataset. All 19 specimens shared the same haplotype, except one specimen from Isolotto Roma [GenBank:JX013545], bearing a single non-synonymous transversion (G to C). Indices of genetic variation are reported in Table 1.
Results obtained by the median-joining network analysis (Figure 2a) provided support for strong genetic differentiation between L. magdalenina and L. incisa. Sixteen (for L. magdalenina) and 12 (for L. incisa) point mutations represent the apomorphic polymorphisms of the two species. Twenty-seven point mutations separated H. macracantha from the group including L. magdalenina and L. incisa.
Cyt b analysis
As a possible consequence of the limited homology between the universal primers  and the annealing region in our specimens, scorable Cyt b sequences were obtained for only ten individuals (three from Punta San Giorgio and seven from Isolotto Roma). The already published Cyt b sequence [GenBank:EU530537] was not included in the dataset since a BLAST search revealed highest similarity with Saprolegnia ferax (a fungus belonging to the family Saprolegniaceae) (81% identity, coverage 87%) and indicated a possible fungal contamination of this sample. A 348 bp long alignment was obtained. In total, four haplotypes were detected. Values of genetic diversity were higher than those obtained for the COI gene (Table 1). No haplotypes were shared between Punta San Giorgio and Isolotto Roma samples. All of the ten polymorphic sites retrieved were biallelic with seven synonymous substitutions. Eight transversions and two transitions occurred.
The median-joining network analysis (Figure 2b) revealed a well-defined divergence, supported by 85 to 125 point mutations (this pattern contains a cycle with three median vectors separated from each other by 39 mutational steps) between L. magdalenina and H. macracantha. Within the L. magdalenina group, the most common haplotype was exclusive of individuals from Isolotto Roma, while two haplotypes belonging to Punta San Giorgio and one from Isolotto Roma diverged by two to five point mutations from the root sequence.
Five decades after the first description of Cephalocarida  our knowledge of this taxon is still limited, especially for species other than Hutchinsoniella macracantha. The present work evidences that if COI universal primers  fail, even a 216 bp long internal portion (amplified by the new primers provided here) of the widely used 710 bp long Folmer  region can be used for the taxonomic identification of the species belonging to the genus Lightiella. Consistently, previous studies [22, 23] showed that partial fragments of the Folmer region of about or less than two hundred bp may correctly assess phylogenetic/phylogeographic traits of the species. In such a context, our preliminary analyses provide molecular support to the morphology-based distinction reported between L. magdalenina and L. incisa .
Sequencing of the Cyt b gene, performed on L. magdalenina specimens, showed higher values of intraspecific genetic variability than the COI gene. Indeed, the absence of shared Cyt b haplotypes between L. magdalenina specimens from the two neighbouring sites sampled in this study is consistent with the possible occurrence of two distinct populations. In spite of the close geographical proximity of the two sites (only 600 m apart), the isolation can be explained by restricted population size  together with a low dispersion capacity of the species . In fact, L. magdalenina, lacks a high-dispersal phase, given that both larval stages and adults show a strictly benthic lifestyle [1, 25]. However, it should be taken into consideration that this pattern could be affected by the small number of samples from Punta San Giorgio (three individuals), thus preventing supported inferences.
In conclusion, this study reports the first data on the genetic variation of a cephalocarid species, providing a preliminary survey on the mitochondrial variability of the only Mediterranean species, L. magdalenina.
The high level of interspecific genetic divergence reported between L. magdalenina and L. incisa, using the COI gene, offers molecular support to previous morphological analysis , allowing to consider this marker as a valuable tool for molecular taxonomic assays within the genus Lightiella. On the other hand, if confirmed by further analysis on a larger number of samples, the estimates of intraspecific variation based on the Cyt b gene for L. magdalenina suggest the possible use of this marker for assessing the phylogeographic patterns of the species.
These results highlight the importance of performing joint (COI and Cyt b) molecular analyses in cephalocarids for further taxonomic revision of the genus Lightiella and population genetic surveys.
Addis A, Biagi F, Floris A, Puddu E, Carcupino M: Larval development of Lightiella magdalenina (Crustacea, Cephalocarida). Mar Biol 2007, 152: 733–744. 10.1007/s00227-007-0735-8
Regier JC, Shultz JW, Zwick A, Hussey A, Ball B, Wetzer R, Martin JW, Cunningham CW: Arthropod relationships revealed by phylogenomic analysis of nuclear protein-coding sequences. Nature 2010, 463: 1079–1083. 10.1038/nature08742
Koenemann S, Jenner RA, Hoenemann M, Stemme T, von Reumont BM: Arthropod phylogeny revisited, with a focus on crustacean relationships. Arthropod Struct Dev 2010, 39: 88–110. 10.1016/j.asd.2009.10.003
Olesen J, Haug JT, Maas A, Waloszek D: External morphology of Lightiella monniotae (Crustacea, Cephalocarida) in the light of Cambrian ‘Orsten’ crustaceans. Arthropod Struct Dev 2011, 40: 449–478. 10.1016/j.asd.2011.04.002
Carcupino M, Floris A, Addis A, CA , Curini-Galletti M: A new species of the genus Lightiella : the first record of Cephalocarida (Crustacea) in Europe. Zool J Linn Soc-Lond 2006, 148: 209–220. 10.1111/j.1096-3642.2006.00237.x
Cals P, Delamare Deboutteville C: Une nouvelle espèce de Crustacé Céphalocaride de l’hémisphère austral. CR Acad Sci III-Vie 1970, 270: 2444–2447.
Spears T, Abele LG: Crustacean phylogeny inferred from 18S rDNA. In Arthropod Relationships. 1st edition. Edited by: Fortey RA, Thomas RH. London: Chapman and Hall; 1998:169–187.
Colgan DJ, McLauchlan A, Wilson GDF, Livingston SP, Edgecombe GD, Macaranas J, Cassis G, Gray MR: Histone H3 and U2 snRNA DNA sequences and arthropod molecular evolution. Aust J Zool 1998, 46: 419–437. 10.1071/ZO98048
Regier JC, Shultz JW: Molecular phylogeny of arthropods and the significance of the Cambrian "explosion" for molecular systematics. Am Zool 1998, 38: 918–928.
Regier JC, Shultz JW: Elongation factor-2: A useful gene for arthropod phylogenetics. Mol Phylogenet Evol 2001, 20: 136–148. 10.1006/mpev.2001.0956
Shultz JW, Regier JC: Phylogenetic analysis of arthropods using two nuclear protein-encoding genes supports a crustacean + hexapod clade. P Roy Soc B-Biol Sci 2000, 267: 1011–1019. 10.1098/rspb.2000.1104
Giribet G, Edgecombe GD, Wheeler WC: Arthropod phylogeny based on eight molecular loci and morphology. Nature 2001, 413: 157–161. 10.1038/35093097
Lavrov DV, Brown WM, Boore JL: Phylogenetic position of the Pentastomida and (pan)crustacean relationships. P Roy Soc B-Biol Sci 2004, 271: 537–544. 10.1098/rspb.2003.2631
Richter S, Olesen J, Wheeler WC: Phylogeny of Branchiopoda (Crustacea) based on a combined analysis of morphological data and six molecular loci. Cladistics 2007, 23: 301–336. 10.1111/j.1096-0031.2007.00148.x
Sanna D, Addis A, Carcupino M, Francalacci P: Molecular data on two mitochondrial genes of a newly discovered crustacean species ( Lightiella magdalenina , Cephalocarida). J Mar Biol Assoc UK 2010, 90: 827–829. 10.1017/S0025315409991160
Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R: DNA primers for amplification of mitochondrial cytochrome c oxidase subunit I from diverse metazoan invertebrates. Mol Mar Biol Biotechnol 1994, 3: 294–299.
Boore JL, Brown WM: Mitochondrial genomes of Galathealinum , Helobdella , and Platynereis : sequence and gene arrangement comparisons indicate that Pogonophora is not a phylum and Annelida and Arthropoda are not sister taxa. Mol Biol Evol 2000, 17: 87–106. 10.1093/oxfordjournals.molbev.a026241
Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res 1994, 22: 4673–4680. 10.1093/nar/22.22.4673
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser 1999, 41: 95–98.
Rozas J, Rozas R: DnaSP version 3: an integrated program for molecular population genetics and molecular evolution analysis. Bioinformatics 1999, 15: 174–175. 10.1093/bioinformatics/15.2.174
Sanders HL: The Cephalocarida, a new subclass of Crustacea from Long Island Sound. Proc Natl Acad Sci USA 1955, 41: 61–66. 10.1073/pnas.41.1.61
Hajibabaei M, Smith MA, Janzen DH, Rodriguez JJ, Whitfield JB, Hebert PDN: A minimalist barcode can identify a specimen whose DNA is degraded. Mol Ecol Notes 2006, 6: 959–964. 10.1111/j.1471-8286.2006.01470.x
Sanna D, Lai T, Francalacci P, Curini-Galletti M, Casu M: Population structure of the Monocelis lineata (Proseriata, Monocelididae) species complex assessed by phylogenetic analysis of the mitochondrial Cytochrome c Oxidase subunit I (COI) gene. Genet Mol Biol 2009, 32: 864–867. 10.1590/S1415-47572009005000076
Addis A, Fabiano F, Delogu V, Carcupino M: Reproduction system morphology of Lightiella magdalenina (Crustacea, Cephalocarida): functional and adaptive implications. Invertebr Reprod Dev 2013, 57: 142–155. 10.1080/07924259.2012.704408
Sanders HL, Hessler RR: The larval development of Lightiella incisa Gooding (Cephalocarida). Crustaceana 1964, 7: 81–97. 10.1163/156854064X00326
We thank La Maddalena Archipelago National Park and especially Mr. Mirko Ugo for logistical and technical support during the collection of specimens. We also thank the two anonymous referees for their comments and criticisms which greatly helped in improving the MS. This work was supported by the Regione Autonoma della Sardegna with funds from Operative Program FSE 2007–2013 L.R.7/2007 “Promozione della ricerca scientifica e dell’innovazione tecnologica in Sardegna” (grant to AA), by MURST (Ministero per l’Università e la Ricerca Scientifica e Tecnologica) ex-60% (grant to PF), and by Fondazione Banco di Sardegna, 2013 (grant to MC).
The authors declare that they have no competing interests.
DS designed the study, conceived and designed the experiments, performed the molecular and statistical analysis and drafted the manuscript. AA participated in the collection of specimens and drafting of the manuscript. FS participated in the molecular analysis and revision of the manuscript. FF participated in collection of specimens and revision of the manuscript. MC and PF conceived the study and participated in its design, coordination, and drafting; also contributed reagents/materials. All authors read and approved the final manuscript.