Submit or Track your Manuscript LOG-IN

Genetic Structure of Octopus minor (Sasaki, 1920) (Cephalopoda: Octopoda) from Chinese Waters using Mitochondrial ATPase 6 Gene




Genetic Structure of Octopus minor (Sasaki, 1920) (Cephalopoda: Octopoda) from Chinese Waters using Mitochondrial ATPase 6 Gene

Faiz Muhammad1,2, Canfeng Dou1, Zhen-ming Lü1,*, Li Gong1, Xun Du1 and Muhammad Shafi3

1National Engineering Research Center of Marine Facilities Aquaculture, College of Marine Sciences and Technology, Zhejiang Ocean University, Zhoushan, Zhejiang, PR China

2Center of Excellence in Marine Biology, University of Karachi, Karachi, Pakistan

3Lasbella University of Agriculture, Water and Marine Sciences, Baluchistan, Pakistan


Genetic variability and structure in Octopus minor was investigated from 154 individuals using mtDNA ATPase 6, sampled from Dalian, Dongshan, Nantong, Qingdao, Shanghai, Wenzhou, Xiamen, and Zhoushan China. Based on 672 bp fragment of Mitochondrial DNA subunit ATPase 6, thirty haplotypes and 140 polymorphic sites were noted, whereas (0.740, 0.048) high and low haplotype and nucleotide diversity observed. The data was further tested with AMOVA that validates 91.24% variation among populations and 8.76% within population. Genetic fixation index was also calculated as 0.91243*. The neighbor-joining phylogenetic analysis revealed three lineages. The Dongshan population distinctly has shown separate lineage. The haplotype networking has shown that Hap1 contributed 48.051% and appeared in six populations whereas, Hap6, Hap11, Hap 16 and Hap17 were also shared by some populations, rest of the haplotypes were unique. This study solely provide species specific information on its genetic structure and population variability in Chinese waters.

Article Information

Received 03 June 2017

Revised 30 August 2017

Accepted 12 December 2017

Available online 06 December 2018

Authors’ Contributions

ZL designed and supervised the work. FM performed experimental work and wrote the paper. XD collected samples and extracted DNA. CD helped in data analysis. LG helped in experiments. MS helped in paper writing.

Key words

Cephalopods, Octopus minor, ATPase 6, Genetic structure, Haplotype networking.


* Corresponding author:

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

Copyright 2019 Zoological Society of Pakistan

Octopus minor (Sasaki, 1920) has profound commercial and medicinal importance (Qian et al., 2010). It is distributed in Chinese waters from Bohai Sea to South China Sea (Dong, 1988, 1991), including Korean Peninsula and Japanese archipelago (Okutani et al., 1987; Lu et al., 2012). It has also synonymized as O. variabilis (Yamamoto, 1942; Dong, 1988; Lu et al., 2012, 2013). Moreover, propagation efficiency of marine organisms immensely influence the genetic structure (Daniels et al., 2002; Liu et al., 2007). Increased distribution capabilities causes no regional divergence in gene frequency (Gyllensten, 1985; Liu et al., 2007), whereas concise dispersal competence isolate gene flow (Riginos and Victor, 2001; Han et al., 2008). O. minor has little dispersal capacity and weak migratory potential; male adults are slow mover or remain sessile crawling away on the seabed or burrowing deep in mud. The juvenile grow spontaneously without a planktonic larval stage leading to low gene flow and genetic differentiation (Dong, 1988, 1991).

An extensive range of investigations has been conducted in population genetics of marine organisms including cephalopod (Shaw, 2003; Gao et al., 2016; Chang et al., 2010; Kang et al., 2012). The genetic studies play a fundamental role in the development and management of resources and contribute immensely in searching population history, genetic diversity and geographical partitioning throughout natural dispersal range of targeted organisms. Due to the recent advancement in population genetics these attributes can directly be analyzed by using genetic markers (Hutchings, 2000). ATPase 6 (mtDNA marker) gene has been considered as fast evolving and useful genetic marker for assessment of population structure (Ovender and Street, 2003). The present investigation was designed to probe genetic variability and structure in O. minor using ATPase 6 gene from eight sampling locations of China.


Materials and methods

Specimens of Octopus minor were collected from commercial trawl catches from eight different geographical locations, Dalian, Donghsan, Nantong, Qingdao, Shanghai, Wenzhou, Xiamen, and Zhoushan (Supplementary Fig. S1). Samples were preserved following standard procedure at the site, and were transported to the laboratory in 95% ethyl alcohol and placed in refrigerator until further investigation. Total genomic DNA was isolated from muscle tissues (Liu et al., 2017), extracted DNA was electrophoresed and concentrated DNA was further diluted with DEPC water to reach working concentration (20ng/ µl) for PCR reactions.

The complete ATPase 6 gene was amplified by PCR, using primers ATP6F2 5´ACATGCCACAATTATCACCA 3´ andATP6R 5´ AGCTTAATAAGAGCGTAACG 3´. The pair of primers were designed from mtDNA gene sequence taken from NCBI accession number AB158363. The PCR mixture consisted of 0.2µM of each primer, 5.0µL of 10XTaq Plus polymerase buffer, 0.2 mM of dNTPs, 2 units of Taq Plus DNA polymerase, and 1 µL of DNA template. The PCR conditions were as follows: denaturation at 94 ◦C for 5 min, 35 cycles at 94 ◦C for 30 s, annealing at 55 ◦C for 30 s, and extension at 72 ◦C for 30s, and the final extension at 72 C for 10 min. The PCR products were electrophoresed on a 1.5% agarose gel to check integrity, and visualized, using Molecular Imager Gel Doc XR system (Bio- Rad USA). The PCR products were sequenced, using ABI 3730 automated sequencer for which the same PCR primers were used. The DNA sequences were aligned using CLUSTALW and sequence composition was estimated using software MEGA 6 (Tamura et al., 2013). Analysis of genetic differentiation, AMOVA, molecular diversity indices, genetic differentiation values and Fst values were determined with software ARLEQUIN (Excoffier and Lischer, 2010). Haplotype and nucleotide diversity was estimated, using DnaSP (Librado and Rozas, 2009). The neighbor joining tree was constructed to evaluate genetic relationship between the populations using MEGA 6 (Tamura et al., 2013). The haplotype networking was created, using NETWORK software version (Bandelt et al., 1999).



ATPase 6 gene was sequenced from 154 individuals of eight populations. Only 672 bp fragment of ATPase 6 of O. minor was used for data analyses. The average nucleotide composition of T, C, A, and G were 42.5%, 18.4%, 31.9% and 7.2%, respectively. A total of 30 haplotypes, 140 polymorphic sites were identified based on 154 individuals from eight populations. The haplotype diversity was higher (0.740) whereas nucleotide diversity was lower (0.04845) in the investigated populations of O. minor. The average number of nucleotide variations, “K” was 32.028. The details of each population are given in Table I. The high haplotype diversity was observed in Zhoushan following Xiamen, Nantong, Qingdao, Wenzhou, Dongsha, Dalian, and Shanghai, respectively.


Table I.- Nucleotide and haplotype diversities for different populations.



P. sites

Hap. (h)

Hap. div.(Hd)

Nucl. div.(π)

Avg. diff. (K)

























































n, sample number; P. sites, number of polymorphic sites; Hap., Haplotype; Hap. div., Haplotype diversity; Nucl. div., Nucleotide diversity; Avg. diff., average number of differences; DL, Dalian; DS, Dongshan; NT, Nantong; Q, Qingdao; SH, Shanghai; W, Wenzhou; X, Xiamen; Z, Zhoushan.


The AMOVA test validates 8.76% within population and 91.24% among populations (Supplementary Table I). The population structure illustrated high and significant FST value of 0.912* at P < 0.05. The Tajima’s D test defined significant negative values ranging (-0.684 - -2.728). The detailed description of Tajima’s D, Fu’s Fs, corresponding P values and mismatch distribution parameter estimates are shown in Supplementary Table II, whereas the pairwise FST between populations of O. minor are described in Table II. The phylogenetic analysis revealed three lineages. The lineage one clustered five populations. Second lineage included the Xiamen and Wenzhou populations while Lineage three distinctly distinguishes Dongshan population (Fig. 1).


Table II.- Pairwise FST between eight populations of O. minor.





















































For abbreviations, see Table I.


The phylogenetic relationship illustrated, using neighbor-joining (NJ) method which exhibit three lineages. The lineage one is extended group that describes populations from six location. However Dongshan population was distinctly different from rest of populations. The Lineage II is supported with a bootstrap value of 64% whereas the Lineage III supported with a bootstrap value of 100%. Lineage II included the individuals of Xiamen and Wenzhou populations, and Lineage III distinctly describes the Donghshan population (Fig. 2).

The haplotype network (Fig. 2) revealed that haplotype 1 (Hap1) contributed 48.051% of total O. minor individuals and appeared in six populations. Hap6 and Hap 11 shared by Nantong and Xiamen populations, whereas the Hap16 and Hap 17 was shared by Wenzhou and individuals of Xiamen population. Rest of the haplotypes are unique in their respective populations.



The genetic structure and relationship of O. minor populations in the present study was analyzed using ATPase 6 gene which is essential for successful recovery and long-term conservation strategies of important cephalopod species. The genetic diversity ensures evolutionary potential and even a small decrease or increase in genetic diversity causes significant threat to wild populations inhabiting variable environment (Yang et al., 2008). Therefore, such studies play major role in collection information on the conservation genetics of a species (Crandall et al., 1999; Divya et al., 2015).

In this investigation, one hundred fifty four individuals from eight different localities of Chinese coastal waters were investigated. The results revealed 140 polymorphic sites, 30 haplotypes and high Hd values and low nucleotide diversity. AMOVA analysis of ATPase 6 sequences revealed genetic variation in O. minor. These variations are similar to findings reported by earlier workers in a number of Cephalopod species (Chang et al., 2010; Guo et al., 2011; Li et al., 2013; et al., 2011).

Wright (1978) described relationship between FST values and degree of differentiation: FST 0-0.05 denotes little differentiation among populations; FST 0.05-0.15 indicate moderate differentiation, whereas FST 0.15-0.25 means high distinction and above 0.25 demonstrate very high differentiation (Gao et al., 2016).

The pairwise FST values shows highest FST value in Wenzhou population then Donghsan population, other populations have comparatively lower FST values. The phylogenetic tree and pairwise distance results may suggest a sub-species of Donghsan population. The morphological analysis also suspected this population as a sub-species.

The genetic diversity among populations or marine organisms influenced by certain factors, such as geographic segregation, habitat differences, and human activities (Hedrick, 1999; Chang et al., 2010). The lifestyle differences make it hard to carry out suitable gene exchange; as a result, variations revealed different geographical clades. Such general rules are also applied to target species because long armed octopus O. minor also has meager dispersal potential and it is adapted to benthic environment. Adult O. minor are sessile and juveniles develop directly, benthic eggs without a planktonic larval stage (Dong, 1988, 1991). Therefore, O. minor fishery require special care for exploitation and management of future resource along the coast of al. (2013) reported that long armed octopus fishery needs a special care for exploitation when many populations involved over a wide range of geographic distance.

This study would solely provide species specific information on its genetic structure and population variability in Chinese waters.



This research was supported by the National Natural Science Foundation of China (NSFC) (41576131) and Talented Yong Scientist Program (PAK-15-012).


Supplementary material

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


Statement of conflict of interest

Authors declared no conflict of interest.



Bandelt, H.J., Forster, P. and Rohl, A., 1999. Mol. Biol. Evol., 16: 37-48.

Chang, K.M., Li, H., Lü, Z.M. and Chi, C.F., 2010. Oceanol. Limnol. Sin., 41: 307-314.

Crandall, K.A., Posada, D. and Vasco, D., 1999. Anim .Conserv., 2: 317-319.

Daniels, S.R., Stewart, B.A. and Cook, P.A., 2002. Hydrobiologia, 468: 171-179.

Dong, Z.Z., 1988. Fauna sinica: Phylum Mollusca (Class Cephalopoda). Science Press, Beijing, China.

Dong, Z.Z., 1991. Biology of the economic species of cephalopods in the World Oceans. Shandong Science and Technology Pres, Jinan, China.

Divya, P.R., Gopalakrishnan, A., Basheer, V.S., Raja Swaminathan, Mohitha, C., Linu J., Kumar R., Manoj, P. and Jena, J.K., 2015. Mitochond. DNA, 26: 189-194.

Excoffier, L. and Lischer, H.E.L., 2010. Mol. Ecol. Res., 10: 564-567.

Guo, B.Y., Zhou, C., Lu, ZM, Li, J.J. and Wu, C.W., 2011. Oceanol. Limnol. Sin., 42: 868-873.

Gao, X., Zheng, X., Bo, Q. and Li, Q., 2016. Biochem. Syst. Ecol., 66: 129-136.

Gyllensten, U.B., 1985. J. Fish Biol., 26: 691-699.

Han, Z.Q., Gao, T.X., Yanagimoto, T. and Sakurai, Y., 2008. Fish. Sci., 74: 770-780.

Hedrick, P.W., 1999. Evolution, 53: 313-318.

Hutchings, J.A., 2000. Nature, 406: 882-885.

Kang, J.H., Kim, Y.K., Park, J.Y., An, C.M. and Jun, J.C., 2012. Mol. Biol. Rep., 39: 8277-8286.

Librado, P. and Rozas, J., 2009. Bioinformatics, 25: 1451-1452.

Li, H.M., Lü, Z.M., Liu, L.Q., Wu, C.W. and Zhang, J.S., 2013. Ocean Limnosin.,44: 626-631.

Liu, J.X., Gao, T.X., Wu, S.F. and Zhang, Y.P., 2007. Mol. Ecol., 16: 275-288.

Liu, Z., Song, N., Yanagimoto, T., Han, Z., Shui, B. and Gao, T., 2017. Pakitan J. Zool., 49: 111-120.

Lu, C.C., Zhang, X.D. and Lin, X.Z., 2012. In: Proceedings of the 1st Mainland and Taiwan symposium of marine biodiversity studies (eds. M. Lin and C.G. Wang). Ocean Press, Beijing, pp. 76-87.

Lü, Z.M., Li, H., Wu.,C.W., Fan, Z.J. and Zhang, J.S., 2011. J. Fish. Sci. China, 18: 29-37.

Lü, Z.M., Liu, L.Q., Li, H., Wu, C.N. and Zhang, J.S., 2013. Biochem. Syst. Ecol., 51: 224-231.

Okutani, T., Tagawa, M. and Horikawa, H., 1987. Cephalopods from continental shelf and slope around Japan: the intensive research of unexploited fishery resources on continental slopes. Jap. Fish. Res. Cons. Soc., Tokyo, pp. 164-165.

Ovender, R.J. and Street, R., 2003. Mar. Freshw. Res., 54: 127-137.

Qian, Y.S., Zheng, X.D., Wang, P. and Li, Q., 2010. Mar. Sci., 34: 14-18.

Riginos, C. and Victor, B.C., 2001. Proc. R. Soc. London Ser. B., 268: 1931-1936.

Shaw, P.W., 2003. Conserv. Genet., 4: 533-535.

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

Wright, S., 1978. Evolution and the genetics of population, variability within and among natural populations. The University of Chicago Press, Chicago.

Yang, B., Chen, X.Y. and Yang, J.X., 2008. Zool. Res., 29: 379-385.

Yamamoto, T., 1942. Japan J. Malacol., 12: 9-20.

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