Evidence of geographical structuring in the Malaysian Snakehead, Channa striata based on partial segment of the CO1 gene

Evidence of geographical structuring in the Malaysian Snakehead, Channa striata based on partial segment of the CO1 gene

Abu Bakar Siti-Balkhis; Amirul Firdaus Jamaluddin Jamsari; Tan Shau Hwai; Zulfigar Yasin; Mohd Nor Siti-Azizah

School of Biological Sciences, Universiti Sains Malaysia, Penang, Malaysia

^{Send correspondence to}

ABSTRACT

Channa striata, locally known as "haruan", is economically important in fisheries and aquaculture industries in several Asian countries. DNA sequencing, based on a partial segment of the Cytochrome oxidase c subunit 1 (CO1) gene, was used to determine genetic variation in C. striata samples from four different populations on the west coast of Peninsular Malaysia. The highest nucleotide and haplotype diversities were observed in the Linggi population (π = 0.0067, h = 0.835), and the lowest in the Timah Tasoh population (π = 0.0008, h = 0.286). Apart from Kajang-Linggi, which was insignificant, F_ST values were significant (p < 0.05) in all pairwise-population comparisons. Consequently, it is inferred that genetic structuring C. striata populations in this region was largely shaped by a common origin, with secondary influences from geographical factors and isolation.

Key words:Channa striata, snakehead, mtDNA, CO1, isolation by distance.

Channa striata, locally known as "haruan", is a carnivorous freshwater air-breather (Jianguang and Fast, 1997). A popular food-fish in Malaysia (Ali, 1999; Mat Jais, 2007), Thailand and Indo-China, it is also cultivated in India, Pakistan and Thailand (Hossain et al., 2008). Very hardy, it can stay alive out of water for a long time, whereby is the normal marketing procedure (Hossain et al., 2008). It is also widely used in traditional medicine for wound-healing (Zakaria et al., 2004). Studies by Ng et al. (2004) showed that a C. striata extract is useful as an alternative treatment in, osteoarthritis, a joint-disorder involving the softening and disintegration of articular cartilage, with subsequent changes in the underlying bones.

Due to its higher mutation rate of base substitution compared to nuclear DNA (Qiongying et al., 2006), maternally inherited mitochondrial DNA is an efficient genetic marker in genetic-differentiation studies. Mitochondrial cytochrome c oxidase 1 (CO1) being the prime DNA bar-coding region for taxonomically identifying animals, such as fishes, insects, amphibians and sponges (Erpenbeck et al., 2005; Seifert et al,. 2007; Alessandrini et al., 2008; Rock et al., 2008; Smith et al., 2008), is sometimes useful for investigating phylogeographic groups within a single species (Yu et al., 2009). The objective herein was to determine the genetic variation in C. striata populations from four different locations on the west coast of Peninsular Malaysia, through DNA sequencing of CO1 of mtDNA.

Forty five samples were collected from two northern sites, Timah Tasoh (Perlis) and Jeniang (Kedah) and two southern Linggi (Negeri Sembilan) and Kajang (Selangor) (Figure 1). With the exception of Kajang, a cultured species of unknown origin, but believed to be from the southern region itself, all were wild populations. Prior to DNA extraction, fin tissues were preserved in TNES-Urea (100 mM Tris-HCL pH 7.5, 125 mM NaCl, 10 mM EDTA pH 7.5, 1% SDS, 3 M Urea), modified according to Valles-Jimenez et al. (2004).

DNA extraction was with an Aquagenomic Kit (Bio-Syntech, Salt Lake City, Utah, USA), according to manufacturer's instructions. A segment of the CO1 mtDNA gene was amplified using a pair of primers - forward primer L6154 (5'-AYC ARC AYY TRT TYT GRT TCT-3') and reverse primer H6556 (5' TGR AAR TGI CGI ACW ACR TA-3') (Telechea et al., 2006). PCR was done in a Peltier thermal cycler (MJ Research Waltham, MA, USA), with the following profile: pre-denaturation at 95 °C for 5 min, 30 cycles of denaturation at 94 °C, annealing at 50 °C and extension at 70 °C for 1 min each, followed by final extension at 72 °C for 5 min. The PCR products were then purified using Wizard^® SV Gel and a PCR Clean-Up System by Promega (Promega Madison, WI) and sequenced on an ABI3730XL Genetic Analyzer (Applied Biosystems, Foster City, CA, USA).

Sequences were edited using MEGA 4.0 (Tamura et al., 2007) and aligned with Clustal W 1.6, implemented in same software. Arlequin version 3.1 (Schneider et al., 2000) was used for calculating nucleotide and haplotype diversities, haplotype frequency and F_ST values. Bonferroni correction was applied, with global significance level at 0.05 to correct for multiple comparison. The relationships among all haplotypes were constructed using the Network program (Bandelt et al., 1999).

Aligned sequences of 364 bp in the mtDNA CO1 gene were obtained. Twelve unique haplotypes were identified from the four populations of 45 individuals. The sequences were all submitted to the Genbank under accession numbers GQ204314, GQ204319, GQ204321, GQ204323, GQ204326, GQ204330, HM192913 and HM192914. The highest nucleotide and haplotype diversities were observed in Linggi (π = 0.007, h = 0.835) and the lowest in Timah Tasoh (π = 0.0008, h = 0.286). Linggi presented 13 polymorphic sites, Jeniang and Kajang 7 each, and Timah Tasoh only one (Table 1).

The minimum spanning-network relationships among all the haplotypes are represented in Figure 1. Except for Timah Tasoh, Hap1, the most common haplotype, was found in the remaining three populations. From this common haplotype, two clades emerged, the northern and the southern. The spanning network showed that geographically close populations, i.e. Timah Tasoh and Jeniang (north) and Kajang and Linggi (south), were also genetically close.

The highest inter-population genetic distance occurred between Linggi and Jeniang, at 0.011, and the lowest between Timah Tasoh and Jeniang, at 0.004 (Table 2). Genetic differentiation, as revealed by pairwise F_ST analysis, was insignificant between the Kajang and Linggi populations, and significant between the remaining two (Table 2). The scatterplot revealed a highly positive correlation of genetic isolation versus geographical distance (Figure 2) thereby fitting the 'isolation by distance' model (Sun et al., 2004). Nevertheless, caution is called for, due to the small sample sizes. Due to geological events, such as river formation (Tzeng et al., 2006), the distribution of freshwater fishes is usually restricted, thereby lowering the chances of free-gene-flow.

Thus, fairly low intra-population genetic diversity (π = 0.0008-0.0047, h = 0.286-0.483) was the case, except in the Linggi population (π = 0.0067, h = 0.835). Esa et al. (2008), also through CO1-gene analysis, reported even lower variability indices (π = 0-0.001 and h = 0.186-0.450) in an endangered fish, Tor tambroides, from Sarawak (Malaysian Borneo). Furthermore, although low (0.4%-1.1%), the inter-population distances in C. striata were higher (0.1%-0.3%) than in Tor tambroides. Thus, the indications of prevailing genetic variability advocates initiating management procedures. Nonetheless, as a consequence of marker characteristics, an underestimation of the true levels of genetic variability is possible. Several recent studies employing the bar-coding approach have shown that the CO1 gene is a more useful marker for unambiguous-species identification, due to its genetically conserved within-species characteristics (Dasmahapatra and Mallet, 2006; Lara et al., 2009; Valdez-Moreno et al., 2009). Notwithstanding, other markers for population studies, such as microsatellites, or more variable mtDNA markers, such as the D-loop, should also be investigated.

The absence of the common haplotype Hap1 in the Timah Tasoh population is the evidence of its isolation. The presence of a dam could be a likely factor, through its construction possibly changing the original riverine environment into a lacustrine. Notwithstanding, although Hap1 is missing, the main Timah Tasoh haplotype (Hap3) only differed from it by two bases. Thus, the most probable scenario for this population would be a shared origin which evolved independently, due to its isolation. Wang et al. (2004) described the common haplotype as having immense potentiality for producing mutational derivatives. Through being more recently derived, the other haplotypes were unique.

Based on this preliminary analysis, it can be inferred that the genetic structuring of the Peninsular Malaysia C. striata population was largely shaped by a common origin, with secondary influences from geographical factors and isolation. Further studies, with more efficient markers and larger populations, especially from the northern and southern regions, are required to corroborate the findings.

Acknowledgments

We would like to thank all the colleagues involved in sample collection and laboratory work. This study was funded by a Universiti Sains Malaysia Research University Grant (1001/PBiologi/8150123) received by the corresponding author, and a Universiti Sains Malaysia Fellowship by the first.

Received: July 7, 2010; Accepted: March 23, 2011.

Associate Editor: Louis Bernard Klaczko

License information: This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

• Alessandrini F, Mazzanti M, Onofri V, Turchi C and Tagliabracci A (2008) MtDNA analysis for genetic identification of forensically important insects. Forensic Sci Int Genet Suppl Ser 1:584-585.
• Ali A (1999) Aspects of the reproductive biology of female snakehead (Channa striata Bloch) obtained from irrigated rice agroecosystem, Malaysia. Hydrobiologia 411:71-77.
• Bandelt HJ, Forster P and Röhl A (1999) Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol 16:37-48.
• Dasmahapatra KK and Mallet J (2006) DNA barcodes: Recent successes and future prospects. Heredity 97:254-255.
• Erpenbeck D, Hooper JNA and Worheide G (2005) CO1 phylogenies in diploblasts and the 'Barcoding of Life'-are we sequencing a suboptimal partition? Mol Ecol Notes 6:550-553.
• Esa Y, Siraj SS, Daud SK, Rahim KAA, Rovie J, Japning R and Guan SG (2008) Mitochondrial DNA diversity of Tor tambroides Valenciennes (Cyprinidae) from five natural populations in Malaysia. Zool Stud 47:360-367.
• Hossain MK, Latifa GA and Rahman MM (2008) Observations on induced breeding of snakehead murrel, Channa striatus (Bloch, 1793). Int J Sustain Crop Prod 3:65-68.
• Jianguang XH and Fast AW (1997) A bioenergetics model for an air-breathing fish, Channa striatus Environ Biol Fishes 50:309-318.
• Lara A, Ponce De Leo NJL, Rodríguez R, Casane D, Côté G, Bernatchez L and García-Machado E (2009) DNA barcoding of Cuban freshwater fishes: Evidence for cryptic species and taxonomic conflicts. Mol Ecol Resour 10:421-430.
• Mat Jais AM (2007) Molecular size of the bio-active component from haruan Channa striatus extract. J Appl Sci 7:2198-2199.
• Ng YTM, Shanthi G and Loqman MY (2004) Effect of orally administered Channa striatus extract against experimentally-induced osteoarthritis in rabbits. Intern J Appl Res Vet Med 2:171-175.
• Qiongying T, Huanzhang L, Mayden R and Bangxi X (2006) Comparison of evolutionary rates in the mitochondrial DNA cytochrome b gene and control region and their implications for phylogeny of the Cobitoidea (Teleostei, Cypriniformes). Mol Phylogenet Evol 39:347-357.
• Rock J, Costa FO, Walker DI, North AW, Hutchinson WF and Carvalho GR (2008) DNA barcodes of fish of the Scotia Sea, Antarctica, indicate priority groups for taxonomic and systematics. Antarct Sci 20:253-262.
• Schneider S, Roessli D and Excoffier L (2000) Arlequin v. 2000: A software for population genetic data analysis. Genetics and Biometry Laboratory, University of Geneva, Switzerland.
• Seifert KA, Samson RM, deWaard JR, Houbraken J, Lévesque CA, Moncalvo GLS and Hebert PDN (2007) Prospects for fungus identification using CO1 DNA barcodes, with Penicillium as a test case. Proc Natl Acad Sci USA 104:3901-3906.
• Smith MA, Poyarkov JNA and Hebert PDN (2008) CO1 DNA barcoding amphibians: Take the chance, meet the challenge. Mol Ecol Resour 8:235-246.
• Sun Y, Liu S, Zhao G, He S, Wu Q, Taniguchi N and Yu Q (2004) Genetic structure of Chinese sucker population Myxocyprinus asiaticus in Yangtze River based on mitochondrial DNA marker. Fish Sc 70:412-420.
• Tamura K, Dudley J, Nei M and Kumar S (2007) MEGA4: Molecular evolutionary genetics analysis (MEGA) software v. 4.0. Mol Biol Evol 24:1596-1599.
• Telechea F, Laudet V and Hanni C (2006) Phylogeny of the Gadidae (sensu Svetovidov, 1948) based on their morphology and two mitochondrial genes. Mol Phylogenet Evol 38:189-199.
• Tzeng CS, Lin YS, Lin SM, Wang TY and Feng FY (2006) The phylogeography and population demographics of selected freshwater fishes in Taiwan. Zool Stud 45:285-297.
• Wang JP, Lin HD, Huang S, Pan CH, Chen XL and Chiang TY (2004) Phylogeography of Varocorhinus barbatulus (Cyprinidae) in Taiwan based on nucleotide variation of mtDNA and allozymes. Mol Phylogenet Evol 31:1143-1156.
• Valdez-Moreno M, Ivanova NV, Elías -Gutierez M, Contreras-Balderas S and Hebert PDN (2009) Probing diversity in freshwater fishes from Mexico and Guatemala with DNA barcodes. J Fish Biol 74:377-402.
• Valles-Jimenez R, Cruz P and Perez-Enriquez R (2004) Population genetic structure of Pacific White Shrimp (Litopenaeus vannamei) from Mexico to Panama: Microsatellite DNA variation. Mar Biotechnol 6:475-484.
• Yu LL, Xiao YK, Zi NY, Jie K, Shen M and Li MC (2009) Genetic diversity and historical demography of Chinese shrimp Feneropenaeus chinensis in Yellow Sea and Bohai Sea based on mitochondrial DNA Analysis. Afr J Biotechnol 8:1193-1202.
• Zakaria ZA, Somchit MN, Sulaiman MR and Mat Jais AM (2004) Preliminary investigation properties of Haruan (Channa striatus) fillet extracted with various solvent systems. Pakistan J Biol Sc 7:1706-1710.

Publication Dates

• Publication in this collection
  17 June 2011
• Date of issue
  2011

History

• Received
  07 July 2010
• Accepted
  23 Mar 2011