Mitochondrial DNA Markers for PCR-Based Phylogenetic Analysis of Ark Shells ()
1. Introduction
Ark shells are a group of bivalves, in which red blood cells circulate through open blood-vascular systems, and hemoglobin pigment in both blood and tissue cells colors the body red [1]. Approximately 200 species of ark shells have been classified into the family Arcidae [2] and widely distributed globally in shallow tropical and temperate seas [3]. Arcidae is characterized by several derived characters including the duplivincular ligament, taxodont dentition, and shell microstructures consisting of outer crossed lamellar and inner complex crossed lamellar layers [4]. Many Arcidae species have considerable commercial value. In particular, Scapharca, Anadara, and Tegillarca have been actively collected for food throughout Asia [5]. In Japan, Scapharca broughtonii, S. satowi, S. kagoshimensis, S. inaequivalvis, S. globosa ursus, Anadara antiquate, and Tegillarca granosa are consumed by humans (Figure 1). Ark shells vary in size and are mostly elongate or squarish in shape. Shell surfaces are sculptured by radial ribs and covered with thick velvety periostracum [6]. These morphological features are important indicator of species identification, although the radial ribs are often similar between species (Table 1). Hence, morphologically identification of sympatrically related species is difficult, because their radial ribs are similar in number and change with maturity. Hence, genetic markers are required for accurate ark shells identification.
The use of molecular techniques to analyze relationships between species and populations has become widespread, and nucleotide sequence variations in the mitochondrial DNA (mtDNA) offer a powerful tool in molecular phylogenetics of marine organisms [7]. MtDNA markers have been widely used for discrimination of closely related species [8] and also for analyses of population genetic structure [9], however, there are few reports on the genetic studies of Arcidae. Such genetic data are of particular interest because Arcidae may reveal molecular evidence of gene flow or genetic isolation, which is undetectable using traditional morphological studies. In this study, genetic differences in the mtDNA cytochrome c
Figure 1. Photos of left valves of 7 Arcidae species.
Table 1. Morphological and geographical characterization of 7 Arcidae species based on Okutani [6] and sample profiles in this study.
oxidase subunit I (COI), 16S ribosomal RNA (16S rRNA), cytochrome b (Cyt b), and 12S ribosomal RNA (12S rRNA) genes between Scapharca, Anadara, and Tegillarca species were investigated using PCR amplification and subsequent sequence analysis to develop markers that discriminate between Arcidae species in Japan.
2. Materials and Methods
2.1. Samples
A total of 32 specimens were collected from May 2008 to November 2009, and their shell characteristics were recorded (Table 1). Adductor muscle was obtained from each specimen and immediately stored at −20˚C until DNA extraction.
2.2. DNA Extraction
High quality total genomic DNA was prepared from small scraps of frozen adductor muscle according to the modified urea-SDS-proteinase K method [10-12]. Samples were incubated in the extraction buffer (10 mM TrisHCl, pH 7.5, 20 mM EDTA, pH 8.0, 1% SDS, and 4 M urea) containing 25 µg proteinase K at 55˚C and 5 M NaCl was then added and mixed. DNA was isolated with phenol-chloroform-isoamyl alcohol and subsequent chloroform-isoamyl alcohol followed by precipitation with ethanol. DNA pellets were washed with ethanol, dried, and resuspended in 10T0.1E (10 mM Tris-HCl, pH 7.5, 0.1 mM EDTA, pH 8.0).
2.3. PCR Amplification
PCR amplification of the mtDNA regions was performed in GoTaq Green PCR Master Mix (Promega) containing 2 mM MgCl2, 0.5 µM each primer listed in Table 2, and template DNA in a Techgene thermal cycler (Techne). PCR protocol consisted of an initial denaturation at 95˚C
for 2 min, followed by 40 cycles of 10 sec at 95˚C, 20 sec at annealing temperature listed in Table 2, and 40 sec at 72˚C, and a final extension at 72˚C for 5 min. PCR products were analyzed using a DNA-1000 Reagent Kit (Shimadzu) containing a SYBR Gold Nucleic Acid Gel Stain (Invitrogen) in a MCE-202 MultiNA microchip electrophoresis system (Shimadzu).
2.4. Phylogenetic Analysis
Nucleotide sequencing of double strands of PCR product was accomplished using a BigDye Terminator version 3.1 Cycle Sequencing Kit (Applied Biosystems) with an automated 3730xl DNA Analyzer (Applied Biosystems). Nucleotide sequences of PCR products were aligned using CLUSTAL W [19] in MEGA version 5.5 [20] and deposited in GenBank (Scapharca broughtonii, AB690346; S. satowi, AB690347; S. kagoshimensis, AB690348; S. inaequivalvis, AB690349; S. globosa ursus, AB690350; Anadara antiquata, AB690351; and Tegillarca granosa, AB690352). A neighbor-joining tree with bootstrap analysis was constructed from 10,000 replicates based on the Kimura 2-parameter distance [21] using Cucullaea labiata (AB050892) and Barbatia fusca (AB050899) as outgroups.
3. Results
PCR tests were performed using 8 primer pairs to verify useful regions of the mtDNA for the phylogenetic analysis of 7 Arcidae species (Table 3).
Firstly, a universal primer pair COI-1 developed for PCR amplification of the mtDNA COI gene was examined. This primer pair amplified the mtDNA COI gene from 4 Arcidae species such as S. broughtonii, S. globosa ursus, A. antiquata, and T. granosa. Secondly, primer pairs COI-2 and COI-3 developed for PCR amplification of the bivalve the mtDNA COI gene were examined. Neither of these primer pairs amplified the mtDNA COI gene from 7 Arcidae species. Thirdly, a unique primer pair COI-4 developed for PCR amplification of Scapharca mtDNA COI gene was examined. The forward primer of COI-4 was designed as a nested PCR primer inside of the forward primer of the COI-1 [2]. This primer pair successfully amplified the mtDNA COI gene from 7 Arcidae species (Figure 2). Finally, universal primer pairs developed for PCR amplification of the mtDNA 16S rRNA, 12S rRNA, and Cyt b genes were examined. The primer pairs of 16S-1, 16S-2, and Cyt-1 did not consistently amplify the respective genes from 7 Arcidae species, though a small amount of the PCR product was obtained from S. kagoshimensis by Cyt-1. A fragment of approximately 630 bp was amplified from the mtDNA COI gene in all 7 Arcidae species using the unique primer pair COI-4, and a 481 bp portion of complementary nucleotide sequence was obtained after correction and alignment. The average base composition of these sequences was calculated to be A—25.10%, T— 37.84%, C—15.95%, and G—21.12%. Consistent with the known composition of mtDNA, this base composition was biased towards adenine and thymine (62.94%). Alignment of the nucleotide sequences of the mtDNA COI gene of 7 Arcidae species revealed a relatively high level of variability, and 111 variable sites were determined from a total of 481 sites. Of 111 variable sites, there were 73 transitions and 45 transversions, and neither insertion nor deletion. In addition, 39 variable sites at the first or second codon positions occurred 31 nonsynonymous amino acid substitutions (Table 4). The lowest and highest degrees of nucleotide sequence substitution were calculated to be 0.4% between S. broughtonii and T. granosa and 20.2% between S. satowi and A. antiquata, respectively (Table 5). Phylogenetic analysis using Cucullaea labiata and Barbatia fusca as outgroups
Table 3. PCR amplification of 7 Arcidae species.
Figure 2. Electropherograms for PCR amplification of the mtDNA COI gene from 7 Arcidae species. L, lower marker; U, upper marker; arrow, PCR product.
Table 4. Variable sites of the mtDNA COI gene of 7 Arcidae species.
Table 5. Pairwise sequence divergences of the mtDNA COI gene of 7 Arcidae species.
suggested that 7 Arcidae species formed a monophyletic clade, which was supported by higher bootstrap values than 69. These analyses additionally indicated that A. antiquata and T. granosa could be the sister group to Scapharca species (Figure 3).
4. Discussion
In recent years, bivalve resources have been threatened by overfishing or ecological changes, resulting in decreases in wild stocks in Japan [22]. In an effort to promote stock restoration, bivalve seeds were imported and released into coastal waters for farming, and seeds of many mollusk species including S. broughtonii and S. subcrenata were mixed with clam seeds and released into fishing areas of Japan [23].
Despite the commercial importance of Arcidae in Japan, to date little is known about its genetic characteristics and phylogenetic relationships in the native range [24]. Recent developments of DNA markers and genetic identification keys have helped to resolve questions surrounding bivalve taxonomic profiles. Molecular diagnostic markers were specifically developed to identify S. broughtonii in Korea and China [25,26], but these remain unavailable in Japan. In addition, molecular diagnostic markers for discriminating sympatrically distributed, closely related species belonging to Arcidae was also unavailable. In this study, we developed a reliable molecular diagnostic marker, which genetically discriminates between 7 Arcidae species belonging to Scapharca, Anadara, and Tegillarca based on the species-specific polymorphic segments of mtDNA COI gene.
Lee and Kim [2] reported genetic differentiation among S. broughtonii, S. satowi, and S. subcrenata in Korea based on a 599 bp portion of the mtDNA COI gene. The number and degree of nucleotide sequence variations among these 3 Scapharca species ranged from 27% to 48% and from 4.7% to 12.4%, respectively. Those among the same Scapharca species in Japan, however, ranged from 49% to 70% and from 11.2% to 16.5%, respectively, despite the shorter 481 bp portion of the corresponding
Figure 3. Neighbor-joining tree of 7 Arcidae species with Barbatia fusca and Cucullaea labiata as outgroup inferred from the mtDNA COI gene. Bootstrap values higher than 50 are shown at nodes.
the mtDNA COI gene (Table 5). These results suggest that genetic differentiation among S. broughtonii, S. satowi, and S. subcrenata may have progressed in Japan more than in Korea.
Matsumoto and Hayami [24] reported genetic relationship among S. broughtonii, S. satowi, S. subcrenata, Diluvarca tricenicosta, and Tegillarca granosa in Japan, A. antiquata in Philippine and T. nodifera in China based on deduced 302 amino acid sequences of the mtDNA COI gene. The neighbor-joining tree formed a monophyletic clade with Scapharca, Diluvarca, Anadara, and Tegillarca with high bootstrap support.
In the present study, Figure 3 showed a monophyletic clade including Scapharca, Anadara, and Tegillarca, using a neighbor-joining tree. These nucleotide and amino acid sequences of the mtDNA COI gene indicated a monophyletic clade of closely related genera that belong to Anadarinae. In conclusion, the mtDNA COI gene was verified as a sensitive DNA marker for species identifi- cation of Arcidae species. Although DNA sequencing is undoubtedly the most powerful approach, it is costly and time-consuming. Single nucleotide polymorphisms (SNPs) genotyping has been extensively used for species identification, even in bivalves [27]. We accordingly determined several sites containing species-specific SNPs, and developed a simple and easily accessible SNPs genotyping method to identify Arcidae species using speciesspecific PCR amplification of the mtDNA COI gene.
5. Acknowledgements
We thank to Dr. Ayako Yamakawa from Okinawa International University, Dr. Keiko Yamaguchi from Shimane University and Ms. Kumiko Tsuru from Oita Prefectural Agriculture, Forestry and Fisheries Research Center for providing the samples. This research was supported in part by Grant-in-Aid for JSPS Fellows Grant Number 24-2801 to Tomomi Tanaka as JSPS Research Fellow DC2, and by grants from the Agriculture, Forestry and Fisheries Council of Japan.