ICES Journal of Marine Science: Journal du Conseil Advance Access originally published online on November 28, 2006
ICES Journal of Marine Science: Journal du Conseil 2007 64(2):328-331; doi:10.1093/icesjms/fsl023
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Identification of Girella punctata and G. leonina by PCR-RFLP analysis
Department of Marine Science and Resources, Nihon University, Fujisawa, Kanagawa 252-8510, Japan
Correspondence to H. Sugita: tel/fax: +81 466 84 3679; e-mail: sugita{at}brs.nihon-u.ac.jp
Itoi, S., Saito, T., Shimojo, M., Washio, S., and Sugita, H. 2007. Identification of Girella punctata and G. leonina by PCR-RFLP analysis. ICES Journal of Marine Science, 64: 328331.Two Girella species, Girella punctata and G. leonina, are sympatric sister species with an extensive overlap in their respective distributions on shallow rocky reefs from Hong Kong to the south of the Japanese Islands. Juveniles of the two species cannot be discriminated easily on the basis of external characters. In this study, after morphological identification of the species, sequencing analysis was carried out for the partial 16S ribosomal RNA gene and for the D-loop region in mitochondrial DNA. A total of 109 specimens was examined. Restriction site mapping of the sequences suggested that the electrophoretic patterns of polymerase chain reaction (PCR)-based restriction fragment length polymorphism (RFLP) analysis of the gene products would produce a species-specific banding pattern. Subsequently, the PCR-RFLP analysis showed that the method was as effective for separating the two morphologically similar species of the genus Girella as the sequencing analysis.
Keywords: D-loop, Girella leonina, Girella punctata, mtDNA, PCRRFLP, 16S rRNA
Received 3 February 2006; accepted 5 October 2006; advance access publication 28 November 2006.
| Introduction |
|---|
|
|
|---|
Marine teleosts belonging to the genus Girella are widely distributed in the Indo-Pan Pacific and eastern Atlantic Ocean (Yagishita and Nakabo, 2003). Among the 15 known species of Girella, three species, i.e. Girella punctata, G. leonina, and G. mezina, coexist in East Asia between Hong Kong and the southern Pacific Ocean off the Japanese Islands. The three species are closely related, although their food habits are slightly different (Yagishita and Nakabo, 2000, 2003). Kanda and Yamaoka (1995) suggested that G. mezina is herbivorous and G. leonine is omnivorous, whereas G. punctata overlaps the other two species in respect of its food preferences. Recently, we estimated that the divergence time between G. mezina and G. punctata + G. leonina was about 7.3 million years ago, whereas that between G. punctata and G. leonina was about 6.0 million years ago (Itoi et al., 2006). We also studied the phylogenetic relationship among the species of the genus Girella and demonstrated that G. punctata and G. leonina in East Asia are sister species.
Identification of G. punctata and G. leonina juveniles on the basis of morphological characters alone is often confusing, whereas the morphology of G. mezina can easily be separated from the other two species (Okuno, 1962, 1963, 1971). Therefore, the phylogenetic relationship between G. punctata and G. leonina has always been a subject of considerable interest with regard to finding distinct features that could discriminate the two species at the juvenile stage. Generally, the 16S ribosomal RNA (rRNA) gene is the region used to investigate the relationship between intergeneric and interfamily phylogenies (Nishida et al., 1998), whereas the D-loop region is adequate for resolving the relationships among closely related taxa, such as local races, subspecies, and sibling species (Moritz et al., 1987).
Here, we establish a simple and highly sensitive method using polymerase chain reaction (PCR)-based restriction fragment length polymorphism (RFLP) analysis for G. punctata and G. leonina, to serve as an important tool for ecological research.
| Material and methods |
|---|
|
|
|---|
A total of 95 specimens of G. punctata were collected from the seashore at Joetsu, Niigata (n = 28, 168478 g); Kashima, Ibaraki (n = 5, 77291 g); Kawasaki, Kanagawa (n = 5, 58139 g); Shimoda, Shizuoka (n = 12, 74116 g); Muroto, Kochi (n = 29, 5501100 g); Izumi, Kagoshima (n = 9, 269448 g); and Goto, Nagasaki (n = 7, 401800 g) in Japan (Figure 1). In addition, 14 specimens of G. leonina were collected from the seashore at Shimoda, Shizuoka (n = 9, 24168 g); Goto, Nagasaki (n = 4, 388900 g); and Izumi, Kagoshima (n = 1, 394 g) in Japan (Figure 1). Immediately after collection, the fish were transported to the laboratory and kept at 20°C until analysis.
|
The fish were examined for morphological characters, including the numbers of lateral line scales with pores, the numbers of spines and soft rays in the dorsal fin and soft rays in the anal fin, the depth of emargination of the caudal fin, and the colour of the opercular flap, as described by Yagishita and Nakabo (2000). Differences in morphology between G. punctata and G. leonina were as follows. The number of lateral line scales with pores in G. punctata is 5056, whereas in G. leonina, it is 6065. The pigmentation of the opercular flap is blackish in G. leonina, whereas in G. punctata, it is not black. The shape of the caudal fin in G. punctata is shallowly emarginated, whereas in G. leonina, it is deeply emarginated, with acute upper and lower lobes.
DNA extraction, PCR amplification, and PCR product sequencing
Total genomic DNA was extracted from the skeletal muscle, liver, or caudal fin of all 109 G. punctata and G. leonina, using the method of Sezaki et al. (1999). DNA fragments corresponding to the D-loop region and the 16S rRNA gene in mitochondrial DNA (mtDNA) were amplified by PCR. Primers 16SAR-L (5'-CGCCTGTTTATCAAAAACAT-3') and 16SBR-H (5'-CCGGTCTGAACTCAGATCACGT-3') used to amplify the partial 16S rRNA gene fragment (Figure 2) were obtained from Palumbi et al. (1991). Primers fDloop_F (5'-TTCCTGGCATTTGGTTCCTACTTCAG-3') (Itoi et al., 2006) and ftRPhe_R (5'-CCATCTTAACATCTTCAGTGTTATGC-3') for amplifying the partial D-loop region flanking by a part of transfer RNA (tRNA)Phe gene (Figure 2) were designed by referring to the corresponding regions in Pagrus major (GenBank database accession number NC_003196
[GenBank]
), Takifugu rubripes (AJ421455
[GenBank]
), and Tranchurus tranchurus (AB108498
[GenBank]
). PCR amplification was performed using a reaction mixture containing genomic DNA as a template, 1 units of Taq DNA polymerase (Promega, Madison, WI, USA), 2 µl of 10 x Taq DNA polymerase buffer (Promega), 0.8 µl of 10 µM primers, 1.6 µl of 2.5 mM dNTPs, and 1.6 µl of 25 mM MgCl2, and the total volume was brought to 20 µl with sterile water. The thermal cycling profile of PCR consisted of initial denaturation at 95°C for 1 min, followed by 40 cycles of denaturation at 95°C for 10 s, annealing at 55°C for 10 s, and extension at 72°C for 45 s, with a final extension step at 72°C for 2 min. PCR products were cloned into the TA site of pGEMT Easy vector (Promega) according to Marchuk et al. (1991), using the Escherischia coli strain DH5
as a host bacterium.
|
Sequencing of PCR products was performed for both strands with an ABI PRISM 3100 Genetic Analyser (Applied Biosystems, Foster City, CA, USA) using a BigDye Terminator v3.1 Cycle Sequencing Ready Reaction Kit (Applied Biosystems). Alignment of nucleotide sequences partially coding for 16S rRNA and D-loop regions in the mtDNA of G. punctata and G. leonina was carried out using CLUSTAL W (Thompson et al., 1994).
PCR-RFLP analysis
The PCR conditions employed for RFLP analysis were the same as described previously. The RFLP analysis of the gene products was performed by digesting 5 µl of each amplified product with 10 units of restriction enzyme; HinfI (Nippongene, Tokyo, Japan) for 16S rRNA and DdeI (Toyobo, Osaka, Japan) or XbaI (Nippongene) for the D-loop region. The reaction was carried out at 37°C for 2 h in the reaction buffer supplied with the kit. The digested samples were subjected to electrophoresis on a 3% agarose gel and stained with ethidium bromide.
| Results |
|---|
|
|
|---|
Sequencing analysis
Partial regions of the 16S rRNA gene and the D-loop region in the mtDNA from 14 G. punctata and 5 G. leonina were amplified by PCR with appropriately designed oligonucleotide primers. A BLAST search of
600 bp nucleotide sequence of the 16S rRNA gene amplified with primers 16SAR-L and 16SBR-H for each of the two Girella species matched exactly the corresponding sequences deposited in the databases (G. punctata, AB208648
[GenBank]
, AB233479
[GenBank]
AB233491
[GenBank]
; G. leonina, AB208649
[GenBank]
, AB233475
[GenBank]
AB233478
[GenBank]
; data not shown). Similarly,
500 bp fragment of the D-loop region amplified with primers fDloop_F and ftRPhe_R also corresponded to the species-specific sequences in the databases (data not shown). Comparison of the partial mtDNA sequences from G. punctata and G. leonina revealed that these two Girella spp. could be distinguished from each other using the RFLP pattern of the 16S rRNA gene generated with HinfI and that of the D-loop region with DdeI or XbaI.
Species identification
A unique restriction site for HinfI was found in the 613 bp 16S rRNA gene product amplified from G. punctata (Figure 2), which resulted in two restriction fragments measuring 490 and 123 bp (Figure 3). In the case of the 497 bp D-loop region amplified with primers fDloop_F and ftRPhe_R, one restriction site for XbaI was found in G. punctata, but not in G. leonina. In contrast, the D-loop region of G. punctata and G. leonina contained one and two sites recognizable by the restriction enzyme DdeI, respectively (Figure 2). As a result, RFLP with DdeI produced fragments of 308 and 189 bp in G. punctata and of 307, 123, and 66 bp in G. leonina (Figure 3). Therefore, the two Girella species, G. punctata and G. leonina, can easily be discriminated by the PCR-RFLP analysis with these three restriction enzymes (Figure 3).
|
| Discussion |
|---|
|
|
|---|
Girella punctata and G. leonina are widely distributed, living on coastal rocky reefs from southern Japan to Hong Kong, and exhibit considerable intraspecific variation with respect to their morphology. Many authors, including Jordan and Thompson (1912), Matsubara (1955), and Araga (1984), described G. punctata as characteristically having 15 dorsal fin spines, one more than found in G. leonina. However, Nakabo (1993) found considerable intraspecific variation in the number of dorsal fin spines in both species. In a recent study, Yagishita and Nakabo (2000) found that the opercular flap of G. punctata was blackish dorsally in some specimens, a feature not described in the previous studies. Owing to such confusion, Girella has always been a subject of controversy with regard to the taxonomic status of its species, particularly those found in East Asia. Specimens were classified variously as Crenidens melanichthys (Richardson, 1846) or G. melanichthys; this name was placed in varying junior synonym with G. punctata or C. leoninus (the original combination for G. leonina; Richardson, 1846; Yagishita and Nakabo, 2000), and G. melanichthys was incorrectly used as the valid name for G. leonina (Yagishita and Nakabo, 2000). Yagishita and Nakabo (2000) recently revised the classification of the genus Girella in East Asia.
In this study, after identification of G. punctata and G. leonina by morphological analysis and nucleotide sequencing, samples were subjected to PCR-RFLP analysis. The results using partial mtDNA sequences (16S rRNA gene and D-loop region) showed that the technique was as effective for separating the constituent species of the genus Girella as the sequencing analysis (Figure 3). The results demonstrate that the phylogeny derived using mtDNA corresponds with the grouping obtained using the morphological characters of G. punctata and G. leonina.
Analytical methods to distinguish two Girella species, G. punctata and G. leonina, from tissues with no morphological determination are of great methodological importance for ecological research on coastal fish. We have shown that the DNA-based method that discriminates between the two Girella species without sequencing their genes can be applied easily in the field. Although it is not possible to discriminate the juveniles of the two species on the basis of their morphological characters alone, the same would be even more difficult at the larval stage. For example, the larvae of Girella species found at Tanoura inlet, Shimoda, Shizuoka, Japan were thought to be G. punctata (Yoshihara, 1998; Yoshihara et al., 1998, 1999). However, after rearing in tanks for several months, the larvae appeared to be G. punctata and G. leonina (unpublished data). Therefore, the PCR-RFLP analysis reported in this study would potentially prove to be a useful tool for identifying the larvae of the two species.
In conclusion, we have been able to distinguish between G. punctata and G. leonina not only on the basis of differences in their morphological characters, but also on the basis of their nucleotide sequences of both partial 16S rRNA gene and the D-loop region in mtDNA. The present study established a molecular biological method for distinguishing between the two Girella species, G. punctata and G. leonina, by PCR-RFLP analysis based on the mtDNA sequences.
| Acknowledgements |
|---|
This study was supported, in part, by the Open Research Center Project of Japan's Ministry of Education, Culture, Science, Sports, and Technology.
| References |
|---|
|
|
|---|
-
Araga C. (1984) Family Girellidae. In Masuda H., Amaoka K., Araga C., Uyeno T., Yoshino T. (Eds.). The Fishes of the Japanese Archipelago(Tokai University Press, Tokyo) pp. 165166 (in Japanese).
Itoi S., Saito T., Washio S., Shimojo M., Takai N., Yoshihara K., Sugita H. (2006) Speciation of two sympatric coastal fish species, Girella punctata and Girella leonina (Perciforms, Kyphosidae). Organisms, Diversity and Evolution doi:10.1016/j.ode.2006.01.002.
Jordan D. S. and Thompson W. F. (1912) A review of the Sparidae and related families of perch-like fishes found in the waters of Japan. Proceedings of the United States National Museum 41:521601.
Kanda M. and Yamaoka K. (1995) Tooth and gut morphology in relation to feeding in three girellid species (Perciformes: Girellidae) from southern Japan. Netherlands Journal of Zoology 45:495512.
Marchuk D., Drumm M., Saulono A., Collins F. S. (1991) Construction of T-vectors, a rapid and general system for direct cloning of unmodified PCR products. Nucleic Acids Research 19:1154.
Matsubara K. (1955) Fish Morphology and Hierarchy. 1.(Ishizakishoten, Tokyo) pp. xi + 789 (in Japanese).
Moritz C., Dowling T. E., Brown W. M. (1987) Evolution of animal mitochondrial DNA: relevance for population biology and systematics. Annual Review of Ecology and Systematics 18:269292.[CrossRef][Web of Science]
Nakabo T. (1993) Girellidae. In Nakabo T. (Ed.). Fishes of Japan with Pictorial Keys to the Species(Tokai University Press, Tokyo) pp. 771 (in Japanese).
Nishida M., Ohkawa T., Iwata H. (1998) Methods of analysis of genetic population structure with mitochondrial DNA markers. Fish Genetics and Breeding Science 26:81100 (in Japanese with English abstract).
Okuno R. (1962) Distribution of youngs of two reef fishes, Girella punctata Gray and G. melanichthys (Richardson), in Tanabe Bay and the relationship found between their schooling behaviors. Publications of the Seto Marine Biological Laboratory 10:293306.
Okuno R. (1963) Observations and discussions on the social behaviors of marine fishes. Publications of the Seto Marine Biological Laboratory 11:281336.
Okuno R. (1971) Ecology of Rocky Fish. , Sogensha, Osaka (in Japanese).
Palumbi S., Martin A., Romano S., McMillan W. O., Stice L., Grabowski G. (1991) The Simple Fool's Guide to PCR, Version 2(Department of Zoology and Kewalo Marine Laboratory, University of Hawaii, Honolulu).
Richardson J. (1846) Report on the ichthyology of the seas of China and Japan (from the Report of the British Association for the Advancement of Science for 1845). (Richard and John E. Taylor, London) pp. 187320.
Sezaki K., Begum R. A., Wongrat P., Srivastava M. P., SriKantha S., Kikuchi K., Ishihara H., et al. (1999) Molecular phylogeny of Asian freshwater and marine stingrays based on the DNA nucleotide and deduced amino acid sequences of the cytochrome b gene. Fisheries Science 65:563570.
Thompson J. D., Higgins D. G., Gibson T. J. (1994) CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Research 22:46734680.
Yagishita N. and Nakabo T. (2000) Revision of the genus Girella (Girellidae) from East Asia. Ichthyological Research 47:119135.
Yagishita N. and Nakabo T. (2003) Evolutionary trend in feeding habits of Girella (Perciformes: Girellidae). Ichthyological Research 50:358366.[CrossRef]
Yoshihara K. (1998) Population dynamics of juvenile nibbler, Girella punctata, at coastal waters of Izu Peninsula. Suisanzoshoku 46:371372.
Yoshihara K., Kurakata S., Fujita C., Yanagihara M., Ikegami T., Wada T. (1998) Movement and growth of tagged juvenile nibbler at the coastal waters. Suisanzoshoku 46:177182.
Yoshihara K., Okamoto H., Kataoka D. (1999) Lengthweight relationship of juvenile nibbler, Girella punctata, at the coastal waters of Izu Peninsula. Suisanzoshoku 47:343348.
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||


