Discordant patterns of genetic variation between mitochondrial and microsatellite markers in Acanthogobius ommaturus across the coastal areas of China

Acanthogobius ommaturus, which belongs to the family Gobiidae, is a euryhaline and demersal fish that is widely distributed in the coastal areas, harbors, and estuaries of China, D. P. R. Korea and Japan. In this study, the genetic diversity and genetic structure of five geographical populations of A. ommaturus was assessed using the mitochondrial hypervariable region gene and microsatellite markers. The results of the two genetic markers indicated that the A. ommaturus populations had a high level of genetic diversity. The mitochondrial marker detected weak genetic differentiation among populations, and the Neighbor-Joining tree showed that there was no obvious pedigree branches and geographic structure as well. However, population of Zhoushan showed significant genetic differentiation with other populations by microsatellite markers. The population of A.ommaturus has not experienced bottleneck effect recently. We speculated that the Pleistocene climate change and juvenile fish dispersal played an important role in the population differentiation of A. ommaturus.

Development and Characterization of Microsatellite Markers for Harpadon nehereus Based on High-Throughput Sequencing and Cross-Species Amplification in Three Myctophiformes Fishes

Harpadon nehereus is a widespread economical fish found in the coastal seas of China and has important ecological value in the marine ecosystem.H_(o)wever,its germplasm resources have been seriously degraded due to natural factors and anthropogenic activities.In this study,high-throughput sequencing was applied to search for microsatellite loci in H.nehereus transcriptome to provide references for its resource conservation and utilization.Polymorphic loci were developed by non-denaturing polyacrylamide gel electrophoresis,and their cross-species amplified ability was detected in three related species.A total of 5652 microsatellites were identified from 16974320 unigenes.Among the primer pairs designed for 100 SSRs for PCR amplification,80%were successfully amplified,and 26 loci were polymorphic with a high number of alleles from 3 to 11 each.The expected(H_(e))and observed(H_(o))heterozygosities were 0.355–0.885 and 0.375–0.958,respectively.Most of the loci were highly polymorphic(polymorphism information content:0.316–0.852;mean:0.713),and these markers can be applied in the population genetic diversity research of H.nehereus.H_(o)wever,the transferability of these primers was low,probably because of the close relation of the collected species.In follow-up work,simple sequence repeats will be excavated with genome-based technologies,and related species will be gathered to address the present inadequacies.

Mapping QTLs Affecting Economic Traits on BTA3 in Chinese Holstein with Microsatellite Markers

It had been demonstrated that the strong and highly significant quantitative trait locus（QTL） can affect protein percentage on Bos Taurus Autosome 3（BTA3） at the position 52 cM, near the microsatellite DIK4353, with the 95% confidence interval spanning from 25 to 57 cM in Chinese Holstein population using QTL-express, MQREML, and GRIDQTL softwares. This study herein focused on such region of fine mapping QTLs for milk production and functional traits with 16 microsatellite markers with coverage of 33 cM between the markers BMS2904 and MB099 on BTA3 in a daughter-designed Chinese Holstein population. A total of 1 298 Holstein cows and 7 sires were genotyped for 16 microsatellites with ABI 3700 DNA sequencer. The variance components QTL linkage analysis（LA） and linkage-disequilibrium（LD） analysis（LA/LD） was performed to map QTLs for 7 traits, i.e., 305-d milk yield, fat yield, protein yield, fat percentage, protein percentage, somatic cell score and persistency of milk yield. Four strong and highly significant QTLs were detected for fat yield, fat percentage, protein percentage and somatic cell score at the position 40, 30, 27 and 26 cM, respectively. Two minor QTLs for milk yield and persistency of milk yield were identified at 42 and 46 cM, respectively. These findings provided a general idea for the fine mapping of the causal mutation for milk production and functional traits on BTA3 in the future.

Development and characterization of novel polymorphic microsatellite markers for the Korean freshwater snail Semisulcospira coreana and cross-species amplification using next-generation sequencing

Korean freshwater snails of the genus Semisulcospira are widely distributed across East Asia.It has been a very popular nutritional food in Korea,and is an ecologically important water quality indicator because it lives only in clean water.However,no microsatellite markers have been generated to study the population genetic diversity of this genus.In the present study,we developed and characterized 18 novel microsatellite loci from Semisulcospira coreana genomic DNA.The microsatellites were isolated using 454 GS-FLX titanium sequencing and 18 markers were used for genotyping in S.coreana.In addition,we also tested the cross-species transferability of the microsatellite markers in four additional Semisulcospira spp.We identified 18 polymorphic loci and the number of alleles per loci,and their polymorphism information content values ranged from 2 to 17 and 0.203 to 0.902,respectively.The observed and expected heterozygosities of the loci ranged from 0.063 to 0.924 and 0.226 to 0.924,respectively.According to the analysis of the cross-species transferability of these markers,four species,S.forticosta,S.gottschei,S.tegulata,and S.libertina,showed a very high transferability(80%–85%).These results show that this set of nuclear markers could be useful for population genetics studies of this species and closely related species.

Intra-breed Genetic Variation of Fragrance Pigs Detected by Microsatellite Markers

[Objective] The genetic background of four types of fragrance pigs was studied using microsatellite molecular markers,in order to fully understand the genetic resources of miniature pigs in China.[Method] Using 27 pairs of microsatellite loci jointly recommended by Food and Agricultural Organization（FAO）and International Society for Animal Genetics（ISAG）,we detected the genotypes of 200 fragrance individuals belonging to four types（Jiuyang fragrance pig,Jianbai fragrance pig,Congjiang fragrance pig and Huanjiang fragrance pig）,and analyzed their Inter-and intra-breed genetic variations.[Result]The 23 loci detected in the test were high polymorphic;the mean heterozygosity（H） of Jiuyang fragrance pig,Jianbai fragrance pig,Congjiang fragrance pig and Huanjiang fragrance pig were 0.683 6,0.667 9,0.697 3 and 0.702 2,and their mean polymorphism information contents（PIC） were 0.6263,0.6063,0.6420 and 0.6415,respectively.[Conclusion]Four types of fragrance pigs detected in the test all had high intra-breed genetic variability.

Development of microsatellite markers and their utilization in genetic diversity analysis of cultivated and wild populations of the mud carp (Cirrhina molitorella)

Microsatellite markers have been increasingly used in genetic studies on fishery species because of their high applicability in selective breeding programs. Here we reported the development of microsatellite markers and their utilization in mud carp （Cirrhina molitorella）. An （CA）15 enriched library has been constructed for mud carp, using the magnetic beads enrichment procedure. Sequence analysis of 60 randomly picked positive colonies indicate that 56 （93.3%） of the colonies contain microsatellites. Microsatellite polymorphism was assessed using 10 mud carp individuals, and 12 microsatellite loci turned out to be polymorphic. We utilized these loci to study the genetic diversity of a wild population （WM） and a cultured population （CM） of the mud carp. A total of 109 alleles were detected with an average of 9.08 alleles per locus. The mean value of the observed heterozygosity of WM and CM was 0.6361 and 0.6417, respectively, and significant decrease of genetic diversity in CM was not observed. The genetic distance between the two populations was 0.1546 and the value of Gsr was 0.0473. This showed that there existed a slight genetic differentiation between WM and CM.

Development and Characterization of Microsatellite Markers in Brassica rapa ssp.chinensis and Transferability Among Related Species

Simple sequence repeat （SSR） or microsatellite marker is a valuable tool for several purposes, such as mapping, fingerprinting, and breeding. In the present study, an intersimple sequence repeat （ISSR）-PCR technique was applied for developing SSR markers in non-heading Chinese cabbage （Brassica rapa）. A total of 190 SSRs were obtained. Among these, AG or CT （54.7%） was the most frequent repeat, followed by AC or GT （31.6%） of the microsatellites. The average number of the SSRs length array was 16 and 10 times, respectively. Based on the determined SSR sequences, 143 SSR primer pairs were designed to evaluate their transferabilities among the related species of Brassica. The number of alleles produced per marker averaged 2.91, and the polymorphism information content （PIC） value ranged from 0 to 0.863 with an average of 0.540. Monomorphism was observed in 16 primer pairs. The transferability percentage in CC genome was higher than in BB genome. More loci occurred in the BBCC genome. This result supported the hypothesis that BB genome was divergent from A and C genomes, and AA and CC genomes were relatively close. The polymorphic primers can be exploited for further evolution, fingerprinting, and variety identification.

Germplasm of breeding Pseudosciaena crocea as revealed by microsatellite markers

The germplasm of breeding large yellow croaker （Pseudosciaena crocea Richardson） was revealed using 12 microsatellite markers. The results showed that the genetic diversities were on a mediated level in the bred Daiqu and Min - Yue stocks and two hybrid groups, as represented by 4.83 of the mean number of alleles and 0.561 of the average observed heterozygosity. The value of pair- wise differentiation coefficient（ FST ） was only 13. 1% between Daiqu and Min -Yue stocks, demonstrating the low level of differentiation and a close relationship. However, STRUCTURE simulations and phylogenetic tree based on the UPGMA method supported that they are geographically different populations of the same species with distinct genetic structures. Examinations of individual admixture showed that Min - Yue stock had been contaminated by alien individuals. Moreover, the genetic structures of the two hybrid groups resembled those of their parents, especially affected more by their female parents. Finally, the values of average observed heterozygosity between parents and their ascendants were compared and tested, as a result of no detectable differences （P 〉 0. 05）.

Genetic variability of half-smooth tongue sole Cynoglossus semilaevis populations using microsatellite markers

Ten highly variable microsatellite loci were performed in order to evaluate the genetic variability of wild and hatchery samples of half-smooth tongue sole Cynoglossu＇s semilaevis. A group of 200 genotypes belonging to four wild samples, Laizhou （LZ）, Weihai （WH）, Qingdao （QD）, Rizhao （RZ） and one hatchery sample, Mingbo （MB）, were screened. All of the： ten microsatellite loci screened in this study showed marked polymorphism. A total of 70 different alleles were observed over all loci. The number of alleles per locus ranged from 4.2 to 1006 the average of expected and observed heterozygosity ranged from 0.67 to 0.82, and from 0：77 to0.87 respectively. A total of 10, 16, 10, 11 and 5 unique alleles each were found in LZ, WH~ QD：RZ and MB populations. The effective number of alleles varied from 2.87 for HSTS_7 to 6.83 for HSTS_h. The number of average genotypes ranged from 6.0 for HSTS a to 15.4 for HSTS_h： As compared with the wild populations, the hatchery population, MB, showed significant genetic changes such as fewer alleles per locus （P〈0.05）, a smaller number of low frequency alleles （P〈0.05）, a small number of unique alleles and a small number of genotypes （P〈0.05）, all indicative of a reduction in genetic diversity.

Novel microsatellite markers suggest the mechanism of parthenogenesis in Extatosoma tiaratum is automixis with terminal fusion

Parthenogenetic reproduction is taxonomically widespread and occurs through various cytological mechanisms, which have different impact on the genetic variation of the offspring. Extatosoma tiaratum is a facultatively parthenogenetic Australian insect （Phasmatodea）, in which females oviposit continuously throughout their adult lifespan irre- spective of mating. Fertilized eggs produce sons and daughters through sexual reproduction and unfertilized eggs produce female offspring via parthenogenesis. Here, we developed novel microsatellite markers for E. tiaratum and characterized them by genotyping indi- viduals from a natural population. We then used the microsatellite markers to infer the cytological mechanism of parthenogenesis in this species. We found evidence suggesting parthenogenesis in E. tiaratum occurs through automixis with terminal fusion, resulting in substantial loss of microsatellite heterozygosity in the offspring. Loss of microsatel- lite heterozygosity may be associated with loss of heterozygosity in fitness related loci. The mechanism of parthenogenetic reproduction can therefore affect fitness outcomes and needs to be considered when comparing costs and benefits of sex versus parthenogenesis.

Morphometric analysis and fluorescent microsatellite markers to evaluate the genetic diversity of five populations of Penaeus japonicus in China

To understand the current status of germplasm resources of Penaeus japonicus,in this study,genetic diversity was analyzed based on morphometric analysis and fluorescent microsatellite markers.For morphometric analysis,the cumulative ratios of the PC1 and PC2(85.90%)were greater than 75%,which can explain the morphological variation among populations.For microsatellite analysis,the allele number in 10 microsatellite loci ranged from 27 to 53,and the effective allele number ranged from 21 to 42.The observed heterozygosity ranged from 0.463 to 0.983,whereas expected heterozygosity ranged from 0.932 to 0.966,and the average polymorphism information content was 0.947,which indicated a high level of genetic diversity.The unbiased genetic distance among populations ranged from 0.230 to 0.431,and the genetic similarity coefficient ranged from 0.654 to 0.782.Analysis of molecular variance(AMOVA)showed that the genetic variation within individuals was the main source of variation,accounting for 88.55%,and only 1.49%of the genetic variation came from among populations.

Microsatellite Analysis of Expressed Sequence Tags and Development of EST-SSR Markers for Melampsora spp.

In this study,64 498 ESTs of Melampsoraspp.was assembled into 1 998 contigs,and 604SSR loci were detected on these contigs,with 736.6 bp containing one SSR on average.Among these SSRs,trinucleotide repeats were the most abundant repeats（44.70%）.As for the composition of microsatellites, AC and AT repeats were the richest motif in dinucleotide repeats,and AGT and AAG repeats were the most frequent motifs in trinucleotide repeats,whereas（AAAN） n and（AAAAN） n repeats were dominant in tetranucleotide and pentanucleotide repeats,respectively.All the dominant repeat motifs of different types of SSRs were rich in A and T.In EST sequences of Melampsoraspp.genome,microsatellites longer than 20 bp accounted for about 15.07%.It was noticeable that microsatellites were highly rich in the expressed sequences of Melampsoraspp.genome,which implied that SSRs played a significant role in triggering the gene mutation in Melampsoraspp.genome.A total of 455 SSR primers were designed according to the detected microsatellites using Primer 5.0 and Oligo 6. 0,and 30 primer pairs were randomly selected for amplification test.Among these primer pairs,27 primer pairs succeed in amplification,with a successful rate of 90%.Eight primer pairs generated polymorphic fingerprints in Melampsoraspp.collected from different poplar genotypes,accounting for 26.7% of the total primer pairs.The EST-SSR markers developed fromMelampsoraspp.EST sequences provided important marker resources for studying Melampsoraspp.from the aspects of pathogen identification and survey of genetic variation.

Isolation and characterization of 49 polymorphic microsatellite loci for Decapterus maruadsi using SLAF-seq, and cross-amplifi cation to related species

Decapterus maruadsi is a commercially important species in China, but has been heavily exploited in some areas. There is a growing need to develop microsatellites promoting its genetic research for the adequate management of this ?shery resources. The recently developed speci?c-locus ampli?ed fragment sequencing(SLAF-seq) is an effcient and high-resolution method for genome-wide microsatellite markers discovery. In this study, 28 905 microsatellites(mono-to hexa-nucleotide repeats) were identi?ed using SLAF-seq technology, of which di-nucleotide was the most frequent(13 590, 47.02%), followed by mono-nucleotide(8 138, 28.15%), tri-nucleotide(5 727, 19.81%), tetra-nucleotide(1 104, 3.82%), pentanucleotide(234, 0.81%), and hexa-nucleotide(112, 0.39%). One hundred and thirty-two microsatellite loci(di-and tri-nucleotide) were randomly selected for ampli?cation and polymorphism, of which 49 were highly polymorphic and well-resolved. The average number of alleles per locus was 13.63, ranging from 4 to 25, and allele sizes varied between 110 bp and 309 bp. The observed heterozygosity( Ho) and expected heterozygosity( He) ranged from 0.233 to 1.000 and from 0.374 to 0.959, with mean values of 0.738 and0.836, respectively. The polymorphism information content(PIC) ranged from 0.341 to 0.941(mean=0.806).However, 12 loci deviated from Hardy-Weinberg equilibrium. Furthermore, transferability tests were also successful in validating the utility of the developed markers in ?ve phylogenetically related species of family Carangidae. A total of 48 microsatellite markers were successfully cross-ampli?ed in Decapterus macarellus, Decapterus macrosoma, Decapterus kurroides, Trachurus japonicus, and Selaroides leptolepis.The present microsatellites provided the ?rst known set of microsatellite DNA markers for D. maruadsi,D. macarellus, D. kurroides, and D. macrosoma, and would be useful for further population genetic and molecular phylogeny studies as well as help with the ?sheries management formulation and implementation of the understudied species.

Genetic Structure and Diversity of Parental Cultivars Involved in China Mainland Sugarcane Breeding Programs as Inferred from DNA Microsatellites

To understand genetic structure and diversity of parental cultivars involved in China Mainland sugarcane breeding programs, 92 elite parents and 4 wild relatives were genotyped with 18 microsatellite DNA markers. The genetic similarity （GS） values among the cultivars ranged from 0.346 to 0.960 with an average of 0.533. Among the introduced cultivars, India accessions had the closest genetic distance to China Mainland accessions （0.447）, while Australia accessions have the furthest distance （0.503）. A comparison of allelic diversity among geographical origins showed that there were 22 China Mainland specific alleles, of which 28% were derived from native S. spontaneaum germplasm in China. Model-based genetic structure, clustering, and principal components analyses consistently revealed there were five groups within the 96 accessions. Groups 1, 2, 4, and 5 consisted of all cultivars and group 3 only contained wild germplasm. Group 2 was characterized as the Introduction group with 46 cultivars predominantly introduced from Australia, Taiwan of China, India, and USA. Groups 1, 4, and 5 consisted of cultivars mostly originated from China Mainland, defined as the Complex group, Yacheng lines group, and F134/CP72-1210 group, respectively, upon their pedigree. By understanding the genetic relationships among the parental cultivars, breeders can gain a rational basis for expanding the gene pool and select the best parental accessions for crossing.

Correlation of Individual Heterozygosity of Microsatellite Marker Loci with Heterosis of Growth Traits in Pig Populations

To investigate the correlation of individual heterozygosity and heterosis of three traits in crossbred F1 pig populations, the F1 populations were built by random mating Yorkshire x Meishan （YM, n = 82）, and its reciprocal （MY, n =47） and two straightbred populations （Yorkshire = 34, Meishan = 55） were used as control groups. The heterosis of birth weight （BWT）, average daily gain （ADG）, and feed conversion ratio （FCR） were acquired as well. In the research, the significant marker loci for the heterosis of the three traits were observed by one-way ANOVA （P〈0.01） in a total of 39 marker loci on SSC4, SSC6, SSC7, SSC8, and SSC13, and the numbers of the significant marker loci were 12 （BWT）, 18 （ADG）, and 17 （FCR）, respectively, based on which the general heterozygosity （GH） was divided into significant marker loci heterozygosity （SH） and insignificant marker loci heterozygosity （IH）. Furthermore, the trends of alteration in heterosis with the stepwise increase in heterozygosity by 0.05 were explored. This was done by the regression analysis of the three kinds of heterozygosity against heterosis of the three traits. The results showed that, for BWT, the heterosis increased with the increase in GH （r=0.9337, P=0.0021） and SH （r=0.9165, P=0.0102）; for ADG, the heterosis increased with the increase in IH （r=0.7012, P=0.0353） and GH （r=0.7470, P=0.0537, near significant）; for FCR, the heterosis of feed efficiency increased with the increase in IH （r=0.8721, P=0.0022）. The results indicated that the correlation was not always higher or more significant for SH with heterosis than it was for IH or GH with heterosis, and it might be because of the reciprocal cancellation of the positive effect and negative effect of QTL linked to the significant marker loci.

应用微卫星标记对海南和广西恒河猴遗传多样性的研究

本研究利用11个微卫星位点,对海南和广西恒河猴进行了遗传检测,并通过POPGEN32软件计算各个微卫星座位的等位基因频率、有效等位基因数目(Ne)、多态信息含量(PIC)和遗传杂合度(H)。结果表明,所选择的11个微卫星位点均存在高度的遗传多态性,H为0.6848~0.河河猴的Ne、PIC和H的平均值均高于海南猴,分别为4.2583、0.7090、0.7706和4.2054、0.7025、0.7656,这种差别可能与地域来源有关。这些研究为微卫星标记分析恒河猴遗传多样性提供理论基础。

Comparative population genetic analyses suggest hybrid origin of Rhododendron pubicostatum,an endangered plant species with extremely small populations endemic to Yunnan,China

Gene flow between sympatric congeneric plants is thought to be very common and may pose serious threats to endangered species.In the present study,we evaluate the genetic diversity and divergence of three sympatric Rhododendron species in Jiaozi Mountain using newly developed microsatellites through the Illumina MiSeq sequencing approach.Genetic diversity of all three Rhododendron species studied was moderate in comparison to genetic parameters previously reported from species of this genus.Interestingly,genetic structure analysis of the three species identified a possible hybrid origin of the threatened Rh.pubicostatum.This sympatry should be considered a unimodal hybrid zone,since Rh.pubicostatum is predominant here.Unimodal hybrid zones are uncommon in Rhododendron,despite the fact that hybridization frequently occurs in the genus.Issues pertaining to the conservation of Rh.pubicostatum resulting from admixture of genetic material from its parental species are discussed.

Genetic structure,gene flow pattern,and association analysis of superior germplasm resources in domesticated upland cotton(Gossypium hirsutum L.)

Gene flow patterns and the genetic structure of domesticated crops like cotton are not well understood.Furthermore,marker-assisted breeding of cotton has lagged far behind that of other major crops because the loci associated with cotton traits such as fiber yield and quality have scarcely been identified.In this study,we used 19 microsatellites to first determine the population genetic structure and patterns of gene flow of superior germplasm resources in upland cotton.We then used association analysis to identify which markers were associated with 15 agronomic traits(including ten yield and five fiber quality traits).The results showed that the upland cotton accessions have low levels of genetic diversity(polymorphism information content=0.427),although extensive gene flow occurred among different ecological and geographic regions.Bayesian clustering analysis indicated that the cotton resources used in this study did not belong to obvious geographic populations,which may be the consequence of a single source of domestication followed by frequent genetic introgression mediated by human transference.A total of 82 maker-trait associations were examined in association analysis and the related ratios for phenotypic variations ranged from 3.04% to 47.14%.Interestingly,nine SSR markers were detected in more than one environmental condition.In addition,14 SSR markers were co-associated with two or more different traits.It was noteworthy that NAU4860 and NAU5077 markers detected at least in two environments were simultaneously associated with three fiber quality traits(uniformity index,specific breaking strength and micronaire value).In conclusion,these findings provide new insights into the population structure and genetic exchange pattern of cultivated cotton accessions.The quantitative trait loci of domesticated cotton identified will also be very useful for improvement of yield and fiber quality of cotton in molecular breeding programs.

Morphological characteristics and genetic differentiation of Lutraria maxima in coast waters off southeast China

To explore genetic diversity and estimate the genetic differences among populations of Lutraria maxima in the coastal waters off south to southeast China,the morphology of the species of five different geographical populations(Beihai,Weizhou Island,Zhanjiang,Xiamen,and Fuzhou)in Guangxi,Guangdong,and Fujian provinces was studied statistically in combination with the microsatellite markers.As revealed by morphological principal component analysis(PCA),the cumulative contribution rate of the first three principal components was 72.596%.The discrimination accuracy ranged from 47.5%to 80.0%,and the scatter plots of principal component and discriminant analysis were consistent in overall,showing that the Xiamen and Fuzhou populations were overlapped obviously.For microsatellite markers,10 pairs of polymorphic primers were obtained by high-throughput transcriptome sequencing,and used for genetic diversity analysis.It was showed that the average number of alleles and eff ective alleles observed in each population ranged from 8.100 to 10.900,and from 3.497 to 4.228,respectively.The average observed heterozygosity(H_(o))and expected heterozygosity(H_(e))in the five populations ranged from 0.541 to 0.615,and from 0.642 to 0.733,respectively.The genetic distance(DA)ranged from 0.078 to 0.523,and the population genetic differentiation index(F_(ST))ranged from 0.027 to 0.139.The unweighted pair-population method with arithmetic means(UPGMA)and structure analysis showed that the five populations could be divided into two main clusters,the Beibu Gulf group(Beihai and Weizhou Island)and the Southeast China Sea group(Zhanjiang,Xiamen,and Fuzhou),suggesting that L.maxima has been separated geographically by the barrier of the Leizhou Peninsula into two groups in evolution,which provided us with a scientific clue to better protect the bioresource and establish an appropriate fishery management stocks for L.maxima populations in south China.

Population Genetic Analysis of the Rice Stem Borer, Chilo suppressalis, in the South China

Genetic variation and patterns of genetic differentiation of the rice stem borer, Chilo suppressalis （Lepidoptera： Pyralidae）, from the South China were analyzed using 6 microsatellite markers and two partial mtDNA （cox1 and cox2） regions. All of the 6 microsatellite loci were polymorphic in the studied seven populations. The allelic richness per population ranged between 5.67 and 14.00, and average H E and H O values were 0.6246-0.8329 and 0.2634-0.6061, respectively. As the mitochondrial genome is a single genetic locus, we only present results for the concatenated data set （cox1 plus cox2 gene sequences, 513 bp）. The concatenated data showed high level of genetic diversity and there are 23 variable polymorphic sites among the 513 sites in concatenated data. Nearly all of （20 of 21） pairwise F ST comparisons among populations showed genetic differentiation with moderate to high pairwise F ST values based on microsatellite markers. However, for the mtDNA data, most of the seven populations did not show significant differentiation with other populations. The differences of population differentiation obtained with the two different genetic markers could be mainly attributed to the different mutation rates of microsatellite and mtDNA. There was not genetic structure existed in these studied populations based on microsatellite loci and mtDNA data. The analysis based on network, mismatch distribution, Tajima＇s D and F S indicated that the studied populations were from the recent same ancestor or the same refuge and followed by a sudden demographic expansion condition.