Molecular Systematics of Mahseers (Cyprinidae) in Malaysia Inferred from Sequencing of a Mitochondrial Cytochrome C Oxidase I (COI) Gene
Total Page:16
File Type:pdf, Size:1020Kb
Pertanika J. Trop. Agric. Sci. 31(2): 263 - 269 (2008) ISSN: 1511-3701 ©Universiti Putra Malaysia Press Molecular Systematics of Mahseers (Cyprinidae) in Malaysia Inferred from Sequencing of a Mitochondrial Cytochrome C Oxidase I (COI) Gene Yuzine B. Esa1,2*, Siti Shapor Siraj2, Siti Khalijah Daud2, Japning Jeffrine Rovie Ryan3, Khairul Adha A. Rahim1 and Soon Guan, Tan2 1Faculty of Resource Science and Technology, Universiti Malaysia Sarawak, 94300 Kota Samarahan, Sarawak, Malaysia 2Department of Biology, Faculty of Science, Universiti Putra Malaysia, 43400 UPM, Serdang, Selangor, Malaysia 3Institute of Biodiversity, Bukit Rengit, 28500 Lanchang, Pahang, Malaysia *E-mail: [email protected] ABSTRACT This study examined the molecular systematics among three Mahseers (Tor douronensis, Tor tambroides and Neolissochilus stracheyi) using partial sequencing of a Cytochrome C Oxidase I (COI) mitochondrial DNA segment (466bp). The phylogenetic results using the Neighbour-Joining (NJ) method supported the monophyletic status (hence the taxonomic status) among the three putative Mahseer species. The close genetic relationships (0.1- 0.4%) found between T. tambroides samples from Peninsular Malaysia (kelah fish) and those from Sarawak (empurau fish) also supported their classification as belonging to the same species. The phylogenetic analysis also showed that the T. douronensis mtDNA consisted of three highly distinct lineages supported by high bootstrap values, with the Sabah samples forming its own cluster. Thus, this phylogenetic study, although based on a limited number of samples and only a single mtDNA gene managed to provide useful insights into the systematic status of the Mahseers found in Malaysia. Keywords: Freshwater fish, Mahseers, COI, molecular systematics INTRODUCTION decreasing rapidly due to environmental Freshwater fishes of the genus Tor Gray, degradation and increased fishing pressure (Ng, commonly known as the Mahseers, belong to 2004). the family Cyprinidae (subfamily Cyprininae) So far, very little taxonomic work to (Inger and Chin, 1962; Mohsin and Ambak, systematically sort out Malaysian Mahseer has 1983; Roberts, 1989; Kottelat et al., 1993). There been documented. The most cited work was by are currently 17 described species under the Mohsin and Ambak (1983) who described Tor genus Tor from all across Asia (Ng, 2004) but tambroides and Tor soro as two valid Mahseers only three species were reported in Malaysia: Tor found in Peninsular Malaysia while a more recent tambroides Bleeker, Tor tambra Valenciennes, and view by Ng (2004) suggested the occurrence of Tor douronensis Valenciennes (Kottelat et al., 1993; three species; T. tambroides, T. tambra and T. Kottelat and Whitten, 1996; Ng, 2004). The douronensis. Other taxonomic works recognized taxonomic status of Tor soro Valenciennes had T. douronensis and T. tambroides as two valid species been revised and it is currently re-classified as (Roberts, 1989; Kottelat et al., 1993; Rainboth, Neolissochilus stracheyi (Rainboth, 1996). Mahseers 1996; Zhou and Chu, 1996), although Roberts are important as food fish as well as ornamental (1999) classified them to be a single species, and and recreational fishes. However, recently the a junior synonym to T. tambra. The presence of population sizes of their natural stocks are the median lobe has been characterized as a Received: 12 July 2007 Accepted: 30 June 2008 * Corresponding Author 14. 25/2008-Yuzine 263 24/3/05, 9:56 AM Yuzine B. Esa et al. diagnostic morphological character distinguishing the present study also aimed to clarify the the genus Tor from the genus Neolissochillus phylogenetic relationships among Mahseer fishes (Rainboth, 1996; Ng, 2004), though it cannot be in Malaysia but with a few different approaches. used consistently as a marker to discriminate First, we utilized direct sequencing of the between fishes of the genus Tor. Thus, the Cytochrome C Oxidase I (COI) mtDNA gene region, application of molecular techniques (such as DNA a gene with a faster evolutionary rate compared sequencing) can provide new and better insights to the 16S rRNA (Simon et al., 1994), and thus into the unresolved taxonomy and phylogenetic capable of providing a better resolution at the relationships of all the putative Mahseers in interspecific level. Secondly, we included Malaysia (Nguyen et al., 2006). additional Tor samples from Peninsular Malaysia Nguyen et al. (2006) recently produced the (kelah fish) and Sabah (pelian fish), to compare first molecular work on Mahseers in Malaysia by with the T. douronensis (semah fish) and T. examining the genetic diversity and phylogenetic tambroides (empurau fish) of Sarawak. Thirdly, relationships of broodstocks of T. douronensis N. stracheyi representing the genus Neolissochilus and T. tambroides cultured in Sarawak (Borneo) were included in the phylogenetic study to through sequencing analysis of the mitochondrial quantify the genetic differences between the two DNA (mtDNA) 16S rRNA gene region. Thus, genera. Fig. 1: Map showing sampling locations and sample size (N) in Peninsular Malaysia, Sarawak and Sabah. TD= T. douronensis, TT= T. tambroides, NS= N. stracheyi 264 Pertanika J. Trop. Agric. Sci. Vol. 31(2) 2008 14. 25/2008-Yuzine 264 24/3/05, 9:56 AM Molecular Systematics of Mahseers in Malaysia Inferred from Sequencing of a Mitochondrial Cytochrome COI Gene MATERIALS AND METHODS between populations was calculated using the A total of 111 individuals of the three putative Tamura-Nei distance (Tamura and Nei, 1993), mahseers (T. douronensis, T. tambroides and N. based on unequal base frequencies and unequal stracheyi) were collected from several locations in ratios of transition to transversion (Ti:Tv) Peninsular Malaysia, Sarawak and Sabah (Fig. 1). implemented in MEGA. Total DNA was extracted using a modified CTAB method (Grewe et al., 1993) in the presence of RESULTS AND DISCUSSION Proteinase K. The isolated genomic DNA was The sequence analysis of the partial COI gene used for the mtDNA analysis. (466 base pairs) revealed a total of 24 haplotypes A 500 bp segment of the COI gene was in the nucleotide data set: 14 haplotypes amplified with the oligonucleotide primers COIf belonging to T. douronensis, six haplotypes (5’ CCTGCAGGAGGAGGAGAYCC 3’, forward) belonging to T. tambroides and four haplotypes and COIe (5’ CCAGAGATTAGA GGGAATC AG belonging to N. stracheyi (Fig. 2). The sequence TG 3’, reverse) (Palumbi et al. 1991). of each of the haplotypes was deposited in the Approximately, 50-100 ng of the template DNA GenBank (GeneBank Reference Numbers: was amplified in a 25 ml reaction mixture DQ532824-DQ532827 and EF192444-192463). In containing 50 mM 10X buffer, 2 mM MgCl2, 0.2 total, 74 (15.9%) variable sites were found, of mM of each dNTP (Promega), 0.1 mM of each which 56 (12.0%) were parsimony informative primer, and 0.5 units of Taq DNA Polymerase sites, while 392 (84.1%) were monomorphic sites. (Promega). The cycle parameters consisted of Transitional changes occurred more frequently 35 cycles of denaturation (95oC, 30 seconds), than transversional changes as is typical of animal annealing (45oC, 30 seconds), and extension mitochondrial genomes (Briolay et al., 1998). (72oC, 60 seconds). The amplified products were The phylogenetic results obtained by using visualized on 1% agarose gel containing ethidium the NJ method supported the monophyletic bromide, run for approximately 30 min at 90 V status among the three mahseers (Fig. 3), and photographed under UV light. The purified although the bootstrap support between T. PCR products were directly sequenced using the tambroides and N. stracheyi was low (58%). The BigDye® Terminator v3.0 Cycle Sequencing kit high genetic divergence separating T. douronensis (ACGT) on an ABI 377 automated sequencer and T. tambroides confirmed their status as distinct (PE Applied Biosystem) using only the forward species (Table 1). Likewise, the high genetic primer (COIf). Sequencing reaction using the divergence separating the N. stracheyi lineage reverse primer (COIe) was subsequently carried from the Tor lineages (7.7-8.7%) also supported out on some of the samples (haplotypes) to its recent reclassification from the genus Tor verify the polymorphism in the DNA sequence into the genus Neolissochilus (Rainboth, 1996). initially detected using the forward primer. The close genetic relationships (0.1-0.4%) The CHROMAS (Version 1.45) program was found between T. tambroides samples from used to display the fluorescence-based DNA Peninsular Malaysia (kelah fish) and those from sequencing analysis results. The multiple Sarawak (empurau fish) supported their sequence alignments were done using the taxonomic status as belonging to the same species CLUSTAL X program version 1.81 (Thompson (Table 1). However, the overall very low level of et al., 1997), and subsequently aligned by eye. genetic differentiation within and among T. The Molecular Evolutionary Genetics Analysis tambroides populations may have resulted from (MEGA) version 3.1 (Kumar et al., 2004) program the limited number of samples (2-17) analysed was used to construct a neighbour-joining (NJ) for each population, but was consistent with (Saitou and Nei, 1987) tree using two indigenous those found by Nguyen et al. (2006) in the cyprinids (Barbonymus gonionotus (Genbank Sarawak populations. accession number: DQ532806) and Barbonymus The phylogenetic analysis also revealed that schwanenfeldii (Genbank accession number: the T. douronensis mtDNA consisted of three DQ532805)) as outgroup