Submit or Track your Manuscript LOG-IN

Genetic Diversity and Geographic Differentiation in Northern Snakehead (Channa argus) Based on Mitochondrial Cyt b Gene




Genetic Diversity and Geographic Differentiation in Northern Snakehead (Channa argus) Based on Mitochondrial Cyt b Gene

Aiguo Zhou1,3,4, Shaolin Xie1,4, Zhenlu Wang1, Lanfen Fan1, Yanfeng Chen2, QiaoYe1, Fang Zeng1 and Jixing Zou1,*

1College of Marine Sciences, South China Agricultural University, Guangzhou 510642, Guangdong, China

2School of Life Science and Engineering, Foshan University, Foshan, Guangdong 528231, China

3Qingyuan North River Fishery Science Institute, Qingyuan 511510, Guangdong, China

4Qingyuan Fisheries Technology Limited Company, Qingyuan 511510, Guangdong, China


The great northern snakehead, white type Channa argus is one of the most important economic and conservational fish in China, but an uncertain classification has existed for many years. In present study, mitochondrial Cyt b gene was used to analyze the relationships and genetic diversity of two color morphs of northern snakehead. For the Chnannidae family, the average genetic distances of Cyt b gene was 0.171 (ranged 0.000 to 0.248) with the inter-species genetic distances. The intra-species genetic distance ranged from 0.000 to 0.002. For the haplotypes of two color morphs of Channa argus, the mean pair-wise genetic distances were estimated as 0.002, which were within the intra-species genetic distance interval for Channa species, indicating that they belonged to the same species, but not subspecies at the molecular level. Moreover, take Lateolabrax maculates and Epinephelus coioides as outgroups, molecular phylogenetic tree showed that all the haploids gather together as a branch and crossing each other. These indicated that the white type Channa argus should be regard as an albino of biocolor type.

Article Information

Received 23 June 2017

Revised 24 July 217

Accepted 09 August 2017

Available online 28 November 2018

Authors’ Contributions

JZ and AZ designed the study and drafted the paper. SX and ZW collected the samples. LF and YC pretreated the samples. FZ and QY analysed the data.

Key words

Channa argus, Cyt b gene, Genetic diversity, Geographic differentiation.


* Corresponding author:

0030-9923/2019/0001-0359 $ 9.00/0

Copyright 2019 Zoological Society of Pakistan

The great northern snakeheads, C. argus (Perciformes, Channoidei, Channidae) is widely distributed in Asia and Africa, which was generally dark green with large black blotches (Berra, 2007; Bhat et al., 2014). While the white type C. argus was only found in the Jialing River in Sichuan (105.05E, 29.58N) of China, which is white without any blotches. Kimura (1934) and Shih (1936) regarded the two color morphs as distinct species. However, some studies suggest that the two color morphs as color variations of one single species (Wang et al., 1992, 1993), and similar findings were obtained via our previous studies (Zhou et al., 2016, 2017). In order to further elucidate the genetic relationship between the two color morphs as well as clarify the genetic relatedness among the family Channidae, mitochondrial Cyt b gene was used to determine the level of sequence divergence between them.

The evolution speed of Cyt b gene was relatively moderate, which is suitable for the analysis of population and interspecific differences. It has been widely used for the analysis of population genetic structure and phylogenetic aspects in recent years, such as birds, voles, insects, shrimp, crabs, shellfish and fish (Awan et al., 2017; Gao et al., 2017; Braby and Zwick, 2015; Kartavtsev et al., 2016; Kim et al., 2015). In this study, we used mitochondrial Cyt b gene sequences as molecular marker to investigate the relationships and genetic diversity of two color morphs of northern snakehead and other Chnannidae family species. Based on the results, the germplasm resources and taxonomic position of white type C. argus were discussed. Our aims were to provide crucial genetic information for fisheries management and conservation in white type C. argus.

Materials and methods

All samples were identified according to morphological characters (Cheng and Zheng, 1987; Courtenay and Williams, 2004). Tissue samples were taken as fin clips and preserved in 95% ethanol. Basic information of sampling sites and size (n) of Channa species were given in Table I.


Table I.- Basic information of sampling sites and size (n) of Channa species.



Location name

GPS coordinate

Altitude (m)

Date (size)


Biocolor type C. argus

Neijiang city, Sichuan province, China

29°35'3.16"N, 105° 2'49.95"E





White type

C. argus

Neijiang city, Sichuan province, China

29°34'14.64"N, 105° 4'1.95"E



Ziyang city, Sichuan province, China

30° 6'28.77"N, 104°38'39.46"E




C. maculata

San Jiaozhen, Zhongshan city, Guangdong province, China

22°39'14.31"N, 113°25'43.95"E




C. maculata x

C. argus

Guangzhou city, Guangdong province, China

23° 5'43.82"N, 113°14'19.60"E




C. asiatica

Guangzhou city, Guangdong province, China

23° 5'48.06"N, 113°14'5.13"E




C. micropeltes

Guangzhou city, Guangdong province, China

23° 5'28.31"N, 113°14'1.34"E




C. striata

Vientiane city, Laos

17°57'52.27"N, 102°35'17.97"E




C. lucius

Ha noi City, Vietnam

21° 1'41.45"N, 105°51'48.94"E




Total genomic DNA was extracted from the caudal fin using a standard extraction kit (DNeasy tissue kit, Baitaike Biotech Co., Ltd., China). The amplified primers were given in follows, Cyt b-F: AAAACCACCGTTGTTATTC; Cyt b-R: GCTCTGACGCTGAGCTAC C. The sizes of the amplified fragment length was 1200 bp (Zhou et al., 2016). The amplification reaction was carried out in 50 μl volumes consisting of 25 μl of 2× PCR mix buffer and 0.5 μl of 2.5U/μl Taq DNA polymerase (Takara Bio Inc., Japan), 2 μl of 100 ng/μl DNA template, 2 μl of 10 mM each primer and 18.5 μl ddH2O. Thermal cycling condition consisted of 94°C for 5 min, followed by 35 cycles of 94°C for 30 sec, annealing at 55°C for 30 sec, and an extension temperature of 72°C for 60 sec, and then followed by a final extension of 72°C for 10 min. The PCR product was purified using the QIAquick PCR purification kit (Qiagen, Hilden, Germany). Sequencing was conducted on an ABI377 automatic sequencer with both forward and reverse primers.

The nucleotide sequences of this study (KU852444-KU852447, KU852458-KU852460, KU852459, KT358955-KT358956) and GenBank database (LT577206, KC823607, KC774689, AP006042), were aligned using Clustal-X (Thompson et al., 1997) and then was edited and connected using BioEdit (Hall, 1999). DNAstar was used to estimate the number of haplotypes and its frequencies (Excoffier et al., 2005). The dendrogram of 12 populations was constructed using Kimura 2-parameter model (K2P) in MEGA 6.0 based on the pairwise genetic distance (Tamura et al., 2013).


Results and discussion

A total of 1026 bp of fragments of Cyt b gene was sequenced for 75 individuals. All the sequences were aligned unambiguously with 807bp based on the sequences of other Chnannidae family species in the NCBI database. The variation sites (310), parsimony informative polymorphic sites (229) and conversion/transversion ratio (1.70) were obtained, which indicated that the nucleotide substitution did not reach saturation, these characteristics suggested that the mtDNA Cyt b gene could be used as an effective molecular marker to identify different Channidae species.


Table II.- Haplotype information of Cyt b gene among white and biocolor type C. argus.


Populations (number of the haplotypes and its frequencies)

Common haplotype

White type

C. argus


Biocolor type

C. argus


Cyt b




WtcaHB2 and BtcaHB1










Note: WtcaH and BtcaH represent the haplotypes of white and biocolor type C. argus. B is Cyt b. The same as below.


Three haplotypes (WtcaHB1–B3) of white type C. argus and one haplotype (BtcaHB1) of biocolor type C. argus, were defined from 60 nucleotide sequences (Table II). The nucleotide sequences of all the test Channa species and haplotypes were submitted in GenBank. The sequences of both the haplotypes WtcaHB2 and BtcaHB1, had the largest distribution.

Considering Lateolabrax maculates and Epinephelus coioides as outgroups, sequences alignment via MEGA6.0, showed that all the haploids gather together as a branch and cross each other. Genetic distance was calculated based on the Kimura 2-parameter model. The result showed that the interspecific genetic distance was 0.000-0.248 among 12 kinds of Channidae species, while the maximum was 0.002 and the intraspecific genetic distance was 0.000 (Supplementary Table SI). Interspecific sequence difference is a prerequisite for accurate identification of species for different species of generic and species level (Peng et al., 2009). Previous research shows that the genetic distances within species are generally less than 0.020, and most are less than 0.010 in 11 phylums 13320 species of animals (Hebert et al., 2003). In the present study, the intraspecific average genetic distance was 0.000 and the maximum genetic distances between different haplotypes was 0.002, which was less than 0.010, respectively. These indicate that the white and biocolor type C. argus have not yet reached the level of subspecies differentiation, and thus it may be able to correct the previous classification for the white type C. argus (Shih, 1936). Previous study also have found that the synonym phenomenon is extremely common in the subfamily species of Cynoglossidae, the genetic distance between Cynoglossus lighti and C. joyneri is merely 0.004, which is considered as same species of different morphological types (Liu et al., 2010).


Comparative analysis of the construction of NJ tree, which showed that all haplotypes were clustered into a single group with a high confidence value between the white and biocolor type C. argus (Fig. 1). As a gene sequence information, phylogenetic information can estimate the biological related group, genetic relationship or effective and reflect the size of the development information of a gene, and estimate the accuracy when constructing a phylogenetic tree (Takashima et al., 2004). From molecular phylogenetic trees, we can see that all the selected Channidae species are clustered together, which was separated from the outgroups, and this was mainly reflected in the species level and above. Meanwhile, all haplotypes were clustered into one and cross each other. These results indicate that the white type should be served as an albino of biocolor type C. argus, which was consistent with the result of our previous studies (Zhou et al., 2016, 2017), and the molecular marker has a certain application value in phylogenetic and germplasm analysis. In order to valid the present study and obtain more supportive data, molecular markers are recommended to be used for further comprehensive analysis.



This work was supported by the Science and Technology Planning Project of Guangdong Province (2017A020225035); Produce-learn-research Project of Guangdong Province (2011B090400270); Fund Fostering Talents for Young Scholars of South China Agricultural University (201707N025); Talent introduction special funds of South China Agricultural University and Scientific Research Staring Foundation for Young Scholars of College of Marine Sciences. We also wish to express our appreciation to our anonymous reviewers for providing valuable comments on the manuscript.


Supplementary material

There is supplementary material associated with this article. Access the material online at:


Statement of conflict of interest

The authors report no conflicts of interest. The authors alone are responsible for the content and writing of this article.



Awan, A.R., Firyal, S., Tayyab, M., Rukh, L., Zia ul Haq, M., Saeed, S. and Wasim, M., 2017. Pakistan J. Zool., 49: 283- 288.

Berra, T.M., 2007. Freshwater fish distribution. University of Chicago Press, Chicago.

Bhat, A.A., Haniffa, M.A., Milton, M.J., Paray, B.A., Divya, P.R. and Gopalakrishnan, A., 2014. Int. J. Biodiv. Conserv., 6: 363-372.

Braby, M.F. and Zwick, A., 2015. Invertebr. System., 29: 487-509.

Cheng, Q.T. and Zheng, B.S., 1987. Systematic synopsis of Chinese fishes. Science Press, Beijing.

Chen, S.J., He, C.B., Mu, Y.L., Liu, W.D., Zhou, Z.C., Gao, X.G. and Cong, L.L., 2008. J. Fish. Sci. China, 15: 12-21.

Courtenay, W.R. and Williams, J.D., 2004. Snakeheads (Pisces, Channidae): A biological synopsis and risk assessment. US Geological Survey.

Ding, R.H., 1994. The fishes of Sichuan, China, 1st edition. Sichuan Science and Technology Press, Chengdu, pp. 554.

Excoffier, L., Laval, G. and Schneider, S., 2005. Evolut. Bioinform., 1: 47.

Gao, J., Yue, L.L., Jiang, X., Ni, L., Ashraf, M.A., Zhou, Y., Li, K. and Xiao, J., 2017. Pakistan J. Zool., 49: 1185-1195.

Hall, T.A., 1999. Nucl. Acids Symp. Ser., 41: 95-98.

Hebert, P.D.N., Cywinska, A. and Ball, S.L., 2003. Proc. R. Soc. London B: Biol. Sci., 270: 313-321.

Kartavtsev, Y.P., Sharina, S.N., Saitoh, K., Imoto, J.M., Hanzawa, N. and Redin, A.D., 2016. Mitochondrial DNA, 27: 667-678.

Kim, S.J., Pak, S.J. and Ju, S.J., 2015. Mitochondrial DNA, 26: 127-128.

Kimura, S., 1934. J. Shanghai Sci. Inst., 3: 11-247.

Liu, S.F., Liu, J.X., Zhuang, Z.M., Gao, T.X., Han, Z.Q. and Chen, D.G., 2010. Biodiv. Sci., 18: 275-282.

Peng, J.L., Wang, X.Z., Wang, D. and He, S.P., 2009. Acta Hydrobiol. Sin., 2: 271-276.

Shih, H.J., 1936. Bull. Fan. Mem. Inst. Biol. Peiping (Zool.), 7: 81-82.

Takashima, S., Ise, H., Zhao, P., Akaike, T. and Nikaido, T., 2004. Cell Struct. Funct., 29: 73-84.

Tamura, K., Stecher, G., Peterson, D., Filipski, A. and Kumar, S., 2013. Mol. Biol. Evolut., 30: 2725-2729.

Thompson, J.D., Gibson, T.J., Plewniak, F., Jeanmougin, F. and Higgins, D.G., 1997. Nucl. Acids Res., 25: 4876-4882.

Wang, J.X., Zhao, X.F., Zhuo, C.W. and Liao, Z.G., 1992. Trans. Oceanol. Limnol., 2: 51-57.

Wang, J.X., Liao, Z.G., Zhang, Z.G. and Zhao, X.F., 1993. J. Southw. China Norm. Univ. (Nat. Sci.), 2: 168-172.

Zhou, A.G., Chen, J.T., Xie S.L., Chen, Y.F. and Zou, J.X., 2016. Mitochondrial DNA Part A, 27: 1419-1420.

Zhou, A.G., Wang, C., Jiang, W.Z., Li, Z.G., Chen, Y.F., Xie, S.L., Luo, J.Z. and Zou, J.X., 2016. Mitochondrial DNA Part A, 6: 1-3.

Zhou, A.G., Xie, S.L., Wang, Z.L., Fan, L.F., Wang, C., Ye, Q., Chen, Y.F. and Zou, J.X., 2017. Mitochondrial DNA Part B, 2: 283-286.

To share on other social networks, click on P-share. What are these?

Pakistan Journal of Zoology


Vol. 51, Iss. 5, Pages 1599-1997


Click here for more

Subscribe Today

Receive free updates on new articles, opportunities and benefits

Subscribe Unsubscribe