Genomic DNA from the fin clip was extracted and sequenced on an Illumina novaseq 6000 (Illumina, USA) high-throughput sequencing platform in accordance with the manufacturer’s instructions. After assembly and annotation, haplotype diversity, TCS network analysis, AMOVA analysis, population pairwise genetic distances, and UPGMA tree construction were conducted based on the concatenated sequences of 18S and ITS2.
In total, eleven nrDNA haplotypes were detected based on the concatenated sequences of 18S and ITS2. Amongst, three haplotypes were the main haplotypes, as representatives of three corresponding Clusters. There were two major Clusters in China, however, the Cluster in Vietnam was significantly divergent from the other two in China, likely due to the lack of river connection between China and Vietnam. Interestingly, based on low FST value, we found that gene flow occurred between the isolated island, Hainan Province, and the mainland China of Guangxi Province, probably as a result of exposed continental shelf connected them during glacial periods. In general, combing our data and literature data, genetic diversity and differentiation of