Deep genome-wide divergences among species in White Cloud Mountain minnow Tanichthys albonubes (Cypriniformes: Tanichthyidae) complex: Conservation and species management implications.

Identifying cryptic species is important for the assessments of biodiversity. Further, untangling mechanisms underlying the origins of cryptic species can facilitate our understanding of evolutionary processes. Advancements in genomic approaches for non-model systems have offered unprecedented opportunities to investigate these areas. The White Cloud Mountain minnow (Tanichthys albonubes) is a popular freshwater pet fish worldwide but its wild populations in China are critically endangered. Recent research based on a few molecular markers suggested that this species in fact comprised seven cryptic species, of which six were previously unknown. Here, we tested six of these cryptic species and quantified genomic interspecific divergences between species in the T. albonubes complex by analyzing genome-wide restriction site-associated DNA sequencing (RADseq) data generated from 189 individuals sampled from seven populations (including an outgroup congeneric species, T. micagemmae). We found that six cryptic species previously suggested were well supported by RADseq data. The genetic diversity of each species in the T. albonubes complex was low compared with T. micagemmae and the contemporary effective population sizes (Ne) of each cryptic species were small. Phylogenetic analysis showed seven clades with high support values confirmed with Neighbor-Net trees. The pairwise divergences between species in T. albonubes complex were deep and the highly differentiated loci were evenly distributed across the genome. We proposed that the divergence level of T. albonubes complex is at a late stage of cryptic speciation and lacking gene flow. Our findings provide new insights into cryptic speciation and have important implications for conservation and species management of T. albonubes complex.

Cryptic species in White Cloud Mountain minnow, Tanichthys albonubes: Taxonomic and conservation implications.

Cryptic species describe two or more species that had mistakenly been considered to be a single species, a phenomenon that has been found throughout the tree of life. Recognizing cryptic species is key to estimating the real biodiversity of the world and understanding evolutionary processes. Molecular methods present an unprecedented opportunity for biologists to question whether morphologically similar populations are actually cryptic species. The minnow Tanichthys albonubes is a critically endangered freshwater fish and was classified as a second-class state-protected animal in China. Previous studies have revealed highly divergent lineages with similar morphological characters in this species. Herein, we tested for cryptic species across the ranges of all known wild populations of this minnow. Using multilocus molecular (one mitochondrial gene, two nuclear genes and 13 microsatellite loci) and morphological data for 230 individuals from eight populations, we found deep genetic divergence among these populations with subtle morphological disparity. Morphological examination found variance among these populations in the number of branched anal-fin rays. Based on genetic data, we inferred eight monophyletic groups that were well supported by haplotype network and population clustering analyses. Species delimitation methods suggested eight putative species in the T. albonubes complex. Molecular dating suggested that these cryptic species diverged in the period from the Pliocene to the Pleistocene. Based on these findings, we propose the existence of seven cryptic species in the T. albonubes complex. Our results highlight the need for a taxonomic revision of Tanichthys. What is more, the conservation status of and conservation strategies for the T. albonubes complex should be reassessed as soon as possible.

Geographical and temporal origins of Neocaridina species (Decapoda: Caridea: Atyidae) in Taiwan.

BACKGROUND: The freshwater species on Taiwan Island have been documented to have originated from mainland China and the Japanese islands from multiple events and by multiple colonization routes. Moreover, the sequences from the mitochondrial DNA cytochrome c oxidase subunit I (COI) have been used for DNA barcoding to identify the species. This study used the COI sequences to identify Neocaridina species in Taiwan and to examine their geographical and temporal origins. RESULTS: In total, 479 specimens were collected from 35 localities, which covered almost all rivers in Taiwan. In addition, some sequences were downloaded from GenBank. The maximum likelihood (ML) tree displayed that all sequences were sorted into 13 taxa (clades), and all sequences in Taiwan were sorted into four clades. The Bayesian skyline plots revealed that these four Neocaridina species have declined recently in Taiwan. CONCLUSIONS: All results support that (1) there are four Neocaridina species in Taiwan, which are N. davidi, N. saccam, N. ketagalan and an undescribed Neocaridina species (N. sp.); (2) these four species colonized Taiwan Island in four colonization events; (3) N. sp. colonized Taiwan first; (4) after the island reached its shape, N. ketagalan and N. saccam colonized Taiwan from the Japanese islands and mainland China, respectively; (5) N. davidi colonized northern Taiwan last; and (6) the cyclic glacial and landform changes in East Asia shaped the colonization events and population structures of the Neocaridina species.

Mitochondrial Genetic Diversity of the Freshwater Snail Melanoides tuberculata.

According to geological history, Peninsular Malaysia and Borneo formed at different times and were once connected during Quaternary glaciations. To determine how this history has influenced phylogeography, our study examined the population genetic structure of the tropical freshwater gastropod Melanoides tuberculata across Peninsular Malaysia and Borneo using the sequences from mitochondrial DNA 16S rRNA and cytochrome oxidase subunit I genes (1168 bp). In total, 104 specimens were collected from seventeen populations. All mtDNA haplotypes were identified as belonging to two highly divergent lineages, and these lineages were almost allopatric in their distributions. Our study found that the freshwater fauna in Malaysia might be divided into four regions: northeast Peninsular Malaysia, northwest Peninsular Malaysia, south Peninsular Malaysia, and Borneo. The phylogeography of M. tuberculata in Malaysia was shaped by the landforms of Peninsular Malaysia and by the paleo-river systems in the Sunda continental shelf. In addition, our study found that these two lineages in Malaysia have invaded the globe. These results suggest that Malaysia is located in important shipping lanes throughout the world, and the populations of M. tuberculate might be widely distributed throughout the world by shipping.

The population history of Garra orientalis (Teleostei: Cyprinidae) using mitochondrial DNA and microsatellite data with approximate Bayesian computation.

BACKGROUND: The South China landmass has been characterized by a complex geological history, including mountain lifting, climate changes, and river capture/reversal events. To determine how this complexity has influenced the landmass's phylogeography, our study examined the phylogeography of Garra orientalis, a cyprinid widely distributed in South China, using sequences from the mitochondrial DNA control region and cytochrome b gene (1887 bp) and polymorphisms of thirteen microsatellite loci. RESULTS: In total, 157 specimens were collected from eight populations. All 88 mtDNA haplotypes were identified as belonging to three major lineages, and these lineages were almost allopatric in their distributions. The results of a statistical dispersal-vicariance analysis suggested that the ancestral populations of G. orientalis were distributed south of the Yunkai Mountains, including on Hainan Island. The mtDNA data revealed a strong relationship between phylogeny and geography. In the microsatellite analysis, a total of 339 alleles with an average of 26 alleles per locus were observed across thirteen microsatellite loci. A clustering algorithm for microsatellite data revealed an admixture-like genetic structure. Although the mtDNA and microsatellite data sets displayed a discordant population structure, the results of an approximate Bayesian computation approach showed that these two markers revealed congruent historical signals. The population history of G. orientalis reflects vicariance events and dispersal related to the complex geological history of South China. CONCLUSION: Our results (i) found that the discordances between mtDNA and microsatellite markers were accounted for by admixtures; (ii) showed that the Wuzhishan and Yinggeling mountain ranges and Qiongzhou Strait were important barriers limiting gene exchange between populations on both sides; (iii) indicated that during glaciation and inter-glacial periods, the strait and continental shelves were exposed and sank, which contributed with the dispersion and differentiation of populations; and (iv) displayed that the admixtures between lineages took place in coastal populations and then colonized the tributaries of the Pearl River.

Phylogeography of the inshore fish, Bostrychus sinensis, along the Pacific coastline of China.

This study assesses the phylogeography of the Chinese four-eyed sleeper (Bostrychus sinensis) with one mitochondrial and one nuclear genes and two morphological characters. Phylogenetic and population genetic analyses of the sequences reveals two phylogeographic lineages from the East and South China Seas, which are corroborated by the morphological data. The vicariance of the two lineages is attributed to the Pleistocene Ice Age exposure of the Taiwan Strait and consequent connection of Taiwan to the mainland, which thereby introduced an ecological barrier to gene flow between populations in the East and South China Seas. The distributions of the two lineages now overlap in the East China Sea and this secondary contact is attributed to biased northward migration along the two main currents of the Taiwan Strait following its interglacial re-flooding. In conclusion, this study reinforces the importance of "vicariance, then secondary contact" due to Late Pliocene and Pleistocene sea-level changes to the phylogeography of marine species. Specifically, it corroborates the importance of Pleistocene sea-level changes in the Taiwan Strait to the phylogeography of Chinese inshore species.

Mitochondrial DNA phylogeography of Semisulcospira libertina (Gastropoda: Cerithioidea: Pleuroceridae): implications the history of landform changes in Taiwan.

The mitochondrial DNA cytochrome c oxidase subunit I sequences from 95 specimens of Semisulcospira libertina in Taiwan were identified as two major phylogroups, exhibiting a southern and northern distribution, north of Formosa Bank and south of Miaoli Plateau. The genetic distance between these two phylogroups was 12.20%, and the distances within-phylogroups were 4.97 and 5.56%. According to a molecular clock of 1.56% per lineage per million years, the divergence time between these two major phylogroups was estimated at 4.94 million years ago (mya), with the two phylogroups forming at 3.64 and 3.75 mya, respectively. Moreover, the geological events have suggested that Taiwan Island emerged above sea level at 4-5 mya, and became its present shape at 2 mya. These results suggested that these two phylogroups might originate from two independent ancestral populations or divergent before colonizing Taiwan. Within South phylogroup, the initial colonization was hypothesized to be in Kaoping River (WT), followed by its northward. The high divergence between south- and north of WT River was influenced by the formation of the Kaoping foreland basins. Within North phylogroup, the colonization was from central sub-region through paleo-Miaoli Plateau to northern and northeastern sub-regions. This study showed that the landform changes might have shaped the genetic structure of S. libertina in concert. Apparently, two cryptic species or five different genetic stocks of S. libertina could be identified; these results are useful for the evaluation and conservation of S. libertina in Taiwan.

Phylogeography of the Chinese beard eel, Cirrhimuraena chinensis Kaup, inferred from mitochondrial DNA: a range expansion after the last glacial maximum.

The Chinese beard eel (Cirrhimuraena chinensis Kaup) is an intertidal fish and a model organism for the study of impacts caused by topological fluctuations during the Pleistocene and current intricate hydrological conditions on fauna living in the coastal areas of China. In this study, we examined the phylogeographical pattern, population genetic profile and demographical history of C. chinensis using mitochondrial DNA (cytochrome b (cyt b) and control region (CR)) from 266 individuals sampled in seven localities across the coastal area of southeastern China. The combined data indicated high levels of haplotype diversity and low levels of nucleotide diversity. Analyses of molecular variance (AMOVA) and FST statistics suggested the absence of a significant population structure across the Chinese coast. Neutrality tests, mismatch distributions and Bayesian skyline plots uniformly indicated a recent population expansion. The phylogeographic structure of C. chinensis may be attributed to past population expansion and long-distance pelagic larval dispersal facilitated by present-day ocean currents.

Phylogeography and genetic differentiation among populations of the Moon Turban Snail Lunella granulata Gmelin, 1791 (Gastropoda: Turbinidae).

We examined the genetic variation and phylogeographic relationships among 10 populations of Lunella granulata from mainland China, Penghu Archipelago, Taiwan Island, and Japan using mitochondrial COI and 16S markers. A total of 45 haplotypes were obtained in 112 specimens, and relatively high levels of haplotype diversity (h = 0.903) and low levels of nucleotide diversity (π = 0.0046) were detected. Four major phylogenetic lineage clusters were revealed and were concordant with their geographic distribution, agreeing with the haplotype network. These results suggested that geographic barrier isolating effects were occurring among the populations. This hypothesis was also supported by a significant genetic differentiation index (FST = 0.709) and by a spatial analysis of molecular variance (SAMOVA) analysis. A mismatch distribution analysis, neutrality tests and Bayesian skyline plots found a single significant population expansion. This expansion occurred on the coast of mainland China before 20-17 ka. Consequently, although the dispersal ability of the planktonic stage and the circulation of ocean currents generally promote genetic exchanges among populations, L. granulata has tended to maintain distinct genetic groups that reflect the respective geographic origins of the constituent lineages. Although the circulation of ocean currents, in principle, may still play a role in determining the genetic composition of populations, long-distance migration between regions is difficult even at the planktonic stage.

DNA barcoding the ichthyofauna of the Beibu Gulf: Implications for fisheries management in a seafood market hub.

The Beibu Gulf in China is situated in the tropics, in the western Pacific Ocean. It is an emblematic region combining proximity to a marine biodiversity hotspot and a major seafood hub. Intensification of marine fishing and ocean warming led to a drastic decline in fish populations in the Beibu Gulf during the last decades. This situation urges the development of molecular resources of the Beibu Gulf fish fauna in order to enable automated molecular identifications at the species level for next-generation monitoring. With this objective, we present the results of a large-scale campaign to DNA barcode fishes of the Beibu Gulf. We successfully generated 789 new DNA barcodes corresponding to 263 species which, together with 291 sequences mined from Genbank and BOLD, resulted in a reference library of 1080 sequences from 285 species. Based on the use of four DNA-based species delimitation methods (BIN, ASAP, mPTP, mGMYC), a total of 285 Molecular Operational Taxonomical Units (MOTUs). A single case of cryptic diversity was detected in Scomberomorus guttatus and a single species pair was not captured by delimitation methods. Intraspecific K2P genetic distances averaged 0.36% among sequences within species, whereas K2P genetic distances among species within genera averaged 6.96%. The most speciose families in open water trawling differ from those at fish market, and discrepancies with historical data are discussed in the light of recently documented stock collapses.

Strict consistency between genetic and topographic landscapes of the brown tree frog (Buergeria robusta) in Taiwan.

Taiwan presents an excellent opportunity to build a phylogeographic paradigm for fine-scaled differentiation occurring within short distances on an single island. Due to the limitation of habitat availability on the island, demographic histories of species in Taiwan were strongly influenced by glacial-interglacial cycles. Nevertheless, there are relatively few studies demonstrating such phylogeographic patterns for islands, especially in subtropical Asia. In this study, we aim to construct the genetic landscape of a philopatric stream frog Buergeria robusta by an intense and fine-scaled collection throughout the island. The deduced genetic landscape of B. robusta presented extremely high congruence with the actual topography of Taiwan. Two major lineages were found on the eastern and the western sides of Taiwan with a non-overlapping distribution, indicating the importance of the Central Mountain Range as the major biogeographic barrier. Both clades showed a strong and congruent tendency of demographic or distributional expansion in recent history based on different analyses. Population expansion of such a subtropical lowland species might be a result from a release of available habitat in post-glacial periods.

Isolation and characterization of polymorphic microsatellite loci from Metapenaeopsis barbata using PCR-Based Isolation of Microsatellite Arrays (PIMA).

The red-spot prawn, Metapenaeopsis barbata, is a commercially important, widely distributed demersal species in the Indo-West Pacific Ocean. Overfishing has made its populations decline in the past decade. To study conservation genetics, eight polymorphic microsatellite loci were isolated. Genetic characteristics of the SSR (simple sequence repeat) fingerprints were estimated in 61 individuals from adjacent seas of Taiwan and China. The number of alleles, ranging from 2 to 4, as well as observed and expected heterozygosities in populations, ranging from 0.048 to 0.538, and 0.048 and 0.654, respectively, were detected. No deviation from Hardy-Weinberg expectations was detected at either locus. No significant linkage disequilibrium was detected in locus pairs. The polymorphic microsatellite loci will be useful for investigations of the genetic variation, population structure, and conservation genetics of this species.

Phylogeographical analysis on Squalidus argentatus recapitulates historical landscapes and drainage evolution on the island of Taiwan and mainland China.

Phylogeographical analyses on Squalidus argentatus samples from thirteen localities within mainland China and Taiwan were conducted for biogeographic studies, as their dispersal strictly depends on geological evolution of the landmasses. A total of 95 haplotypes were genotyped for mtDNA cyt b gene in 160 specimens from nine river systems. Relatively high levels of haplotype diversity (h = 0.984) and low levels of nucleotide diversity (π = 0.020) were detected in S. argentatus. Two major phylogenetic haplotype groups, A and B, were revealed via phylogenetic analysis. The degree of intergroup divergence (3.96%) indicates that these groups diverged about 4.55 myr (million years) ago. Haplotype network and population analyses indicated significant genetic structure (F(ST) = 0.775), largely concordant with the geographical location of the populations. According to SAMOVA analysis, we divided these populations into four units: Yangtze-Pearl, Qiantang-Minjiang, Jiulong-Beijiang and Taiwan groups. Mismatch distribution analysis, neutrality tests and Bayesian skyline plots indicated a significant population expansion for lineage A and B, approximately dated 0.35 and 0.04 myr ago, respectively. We found strong geographical organization of the haplotype clades across different geographic scales that can be explained by episodes of dispersal and population expansion followed by population fragmentation and restricted gene flow.

Multilocus Phylogeography and Population Genetic Analyses of Opsariichthys hainanensis Reveal Pleistocene Isolation Followed by High Gene Flow around the Gulf of Tonkin.

The ichthyofauna of continental islands is characterized by immigration through a land bridge due to fluctuating sea levels. Hainan Island is adjacent to the southern margin of mainland China and provides opportunities for understanding the origin and diversification of freshwater fishes. The aim of our study was to evaluate the level of genetic variation and phylogeographic structure of Opsariichthys hainanensis on Hainan Island and mainland China, using mtDNA cyt b gene (1140 bp) and D-loop (926 bp), nuclear RAG1 gene (1506 bp), and 12 microsatellite loci. Mitochondrial phylogenetic analysis identified five major lineages according to the geographical distribution from different populations. We suggested that two dispersal events occurred: the population in the Changhua River migrated to the Red River (Lineage B), and the populations in the South Hainan region moved northwards to the North Hainan region. However, populations in Northwest Hainan Island dispersed to the populations around the Gulf of Tonkin (Lineage A1) and populations in Northeast Hainan Island dispersed to the populations in mainland China (Lineage A2). Our results indicated that the populations of O. hainanensis suffered a bottleneck event followed by a recent population expansion supported by the ABC analysis. We suggest that O. hainanensis populations were found mostly in the lowlands and a lack of suitable freshwater habitat in southern mainland China and Hainan during the Last Interglacial period, and then expansion occurred during the Last Glacial Maximum.

Complete mitogenomes of four Trichiurus species: A taxonomic review of the T.lepturus species complex.

Four Trichiurus species, T.japonicus, T.lepturus, T.nanhaiensis, and T.brevis, from the coasts of the China Seas, have been identified and their entire mitochondrial genomes (mitogenomes) have been sequenced by next-generation sequencing technology. A comparative analysis of five mitogenomes was conducted, including the mitogenome of T.gangeticus. The mitogenomes contained 16.568-16.840 bp and encoded 36 typical mitochondrial genes (13 protein-coding, 2 ribosomal RNA-coding, and 21 transfer RNA-coding genes) and two typical noncoding control regions. Although tRNAPro is absent from Trichiurus mitogenomes, when compared with the 22 tRNAs reported in other vertebrates, the gene arrangements in the mitogenomes of the studied species are consistent with those in most teleost mitogenomes. The full-length sequences and protein-coding genes (PCGs) in the mitogenomes of the five species had obvious AT biases and negative GC skew values. Our study indicate that the specimens in the Indian Ocean are neither T.lepturus nor T.nanhaiensis but they are T.gangeticus; the Trichiurus species composition in the Indian Ocean is totally different from that in Pacific and Atlantic oceans; there are at least two Trichiurus species in Indian Ocean; and the worldwide systematics and diversity of the genus Trichiurus need to be reviewed.

Genetic pattern and demographic history of cutlassfish (Trichiurus nanhaiensis) in South China Sea by the influence of Pleistocene climatic oscillations.

Trichiurus nanhaiensis is one of the most important commercial fish species in the South China Sea. This study aimed to investigate the level of genetic variation and population genetic structure of T. nanhaiensis in the South China Sea for the first time, using 281 individuals collected from seven locations along the coast of mainland China, Taiwan, and Hainan Island. A high level of haplotype diversity and low nucleotide diversity were detected in the mitochondrial DNA cyt b gene and nuDNA RYR 3 gene. The overall expected heterozygosity (He = 0.693) among the seven populations ranged from 0.681 to 0.706 in microsatellite DNA data, which revealed high levels of genetic diversity. Significant genetic differentiation was found in Taidong populations in Taiwan, revealing the prevention of gene flow caused by the Kuroshio Current. Two major lineages based on the cyt b gene suggested that the Taiwan Strait acted as a geographic barrier for T. nanhaiensis during the glacier periods in the late Pleistocene. The Bayesian skyline plot also revealed that population demographic expansion of T. nanhaiensis was estimated to have occurred in 0.1 Mya. Our results indicated that all populations of T. nanhaiensis had experienced a recent genetic bottleneck following recent expansion based on ABC analysis.

Population genetic diversity and structure of Rhinogobius candidianus (Gobiidae) in Taiwan: Translocation and release.

Rhinogobius candidianus is a freshwater goby distributed in north, northwest, west, and south Taiwan, but this species has been introduced to east Taiwan and became dominant. To investigate its native population genetic diversity and structure and evaluate the sources and diversity of translocated populations, the mitochondrial DNA control region and cytochrome b gene (1981 bp) from 220 specimens were analyzed. These results indicated that (1) the east populations originated from two sources in west Taiwan; (2) translocated populations exist in east Taiwan and south Taiwan; (3) many populations have likely been moved secondarily by human intervention; (4) the effective size of the populations had declined greatly; (5) within the native populations, the ancestral populations colonized Taiwan during the land bridge phase in the Pleistocene through north Taiwan; (6) the landform changes in Taiwan shaped the population structure; and (7) the landforms of the coastline during glaciation also shaped the native range. The low-level genetic diversity, high population differentiation, and population decline greatly suggest the need for resource management and conservation interventions. Four clades (α-δ) should be managed as four distinct evolutionarily significant units, while the translocated populations should be managed as separate management units. Moreover, the translocated populations in east Taiwan should be evaluated and monitored carefully.

Genetic population structure of the alpine species Rhododendron pseudochrysanthum sensu lato (Ericaceae) inferred from chloroplast and nuclear DNA.

BACKGROUND: A complex of incipient species with different degrees of morphological or ecological differentiation provides an ideal model for studying species divergence. We examined the phylogeography and the evolutionary history of the Rhododendron pseudochrysanthum s. l. RESULTS: Systematic inconsistency was detected between gene genealogies of the cpDNA and nrDNA. Rooted at R. hyperythrum and R. formosana, both trees lacked reciprocal monophyly for all members of the complex. For R. pseudochrysanthum s.l., the spatial distribution of the cpDNA had a noteworthy pattern showing high genetic differentiation (FST=0.56-0.72) between populations in the Yushan Mountain Range and populations of the other mountain ranges. CONCLUSION: Both incomplete lineage sorting and interspecific hybridization/introgression may have contributed to the lack of monophyly among R. hyperythrum, R. formosana and R. pseudochrysanthum s.l. Independent colonizations, plus low capabilities of seed dispersal in current environments, may have resulted in the genetic differentiation between populations of different mountain ranges. At the population level, the populations of Central, and Sheishan Mountains may have undergone postglacial demographic expansion, while populations of the Yushan Mountain Range are likely to have remained stable ever since the colonization. In contrast, the single population of the Alishan Mountain Range with a fixed cpDNA haplotype may have experienced bottleneck/founder's events.

Population structure in the endangered cyprinid fish Pararasbora moltrechti in Taiwan, based on mitochondrial and microsatellite DNAs.

The genetic structure of four populations of Pararasbora moltrechti, an endemic species of the Cyprinidae in Taiwan, was investigated based on the genetic variation of mtDNA Cyt b gene and five microsatellite loci. High haplotype diversity (h = 0.92) but low nucleotide diversity (0.004) in mtDNA was detected in this endangered species. In total, 33 haplotypes and four clusters were identified in its mtDNA. Nevertheless, low correspondence was found between geographical division and mtDNA clusters. In contrast, Bayesian cluster analysis of the microsatellite data identified four genetic groups and revealed highly structured populations. Significantly negative Tajima's D statistics and mismatch distribution analyses suggest that P. moltrechti populations may have experienced a demographic expansion. In light of the results of a nested clade analysis of mtDNA haplotypes, we conclude that recent population fluctuations and restricted gene flow played major roles in shaping the spatial genetic structure of P. moltrechti populations.

Development and characterization of 20 microsatellite mrkers for Chinese back sleeper, Bostrychus sinensis.

Twenty microsatellite markers were isolated and characterized from the Chinese black sleeper, Bostrychus sinensis. Loci were screened in 30 individuals from Taiwan. For each locus, the number of alleles varied from 4 to 22 with mean expected and observed heterozygosity of 0.79 and 0.66, respectively. One locus significantly deviated from Hardy-Weinberg equilibrium after Bonferroni correction and no significant linkage disequilibrium was detected. This set of microsatellites will provide a suitable tool for population genetic studies of Chinese black sleeper.