Patterns of population genetic structure and diversity across bumble bee communities in the Pacific Northwest

Patterns of genetic structure and diversity are largely mediated by a species’ ecological niche and sensitivity to climate variation. Some species with narrow ecological niches have been found to exhibit increased population differentiation, limited gene flow across populations, and reduced population genetic diversity. In this study, we examine patterns of population genetic structure and diversity of four bumble bee species that are broadly sympatric, but do not necessarily inhabit the same ecological niche in the Pacific Northwest of the United States. Testing for the effect of isolation by geographic distance (IBD) with linearized F _st and D _est found that Bombus sylvicola and B. mixtus exhibited significant IBD across populations. In contrast, both B. melanopygus and B. flavifrons, two species that are distributed across a broad elevation gradient, exhibited no IBD, a result further corroborated by Bayesian a priori population assignment tests. Furthermore, we discovered that B. sylvicola populations distributed on the Olympic Peninsula have significantly less average allelic diversity than populations distributed in the Cascade Mountains. Our results suggest that populations distributed in the Olympic Mountains represent a distinct genetic cluster relative to the Cascade Mountains, with B. sylvicola and B. mixtus likely experiencing the greatest degree of population genetic differentiation relative to B. flavifrons and B. melanopygus. While bumble bees are known to co-exist across a diversity of habitats, our results demonstrate that underlying population genetic structure and diversity may not necessarily be similar across species, and are largely governed by their respective niches.     

Genetic structure of the populations of <Emphasis Type="Italic">Dactylorhiza ochroleuca</Emphasis> and <Emphasis Type="Italic">D. incarnata</Emphasis> (Orchidaceae) in the area of their joint growth in Russia and Belarus

We carried out an allozyme analysis to investigate polymorphism and genetic structure of the populations of D. incarnata and D. ochroleuca in regions of their joint growth in Russia and Belarus. We found that D. ochroleuca individuals in the populations of the Urals and Siberia, which are distant fragments from the main range of the species, do not differ significantly from individuals within the main part of the area (Belarus) on the basis of the allelic composition of eight gene loci. We revealed that D. ochroleuca and D. incarnata are differentiated by different alleles of the GDH locus. Thus, we established a genetic marker suitable to distinguish these closely related taxa. In addition to the GDH locus, D. ochroleuca and D. incarnata in the places of their joint growth, differ in the allelic structure of the PGI and NADHD loci. D. incarnata from the Urals and Siberia were polymorphic for both loci, and individuals from Belarus were polymorphic for one locus (PGI). In contrast, all D. ochroleuca individuals growing in sympatric populations with polymorphic D. incarnata were homozygous for the same alleles. Thus, comparison of the genetic structure of D. ochroleuca and D. incarnata points to the existence of a genetic isolation and a functioning isolation mechanism even under conditions of their joint growth. We found that the GDH locus in D. incarnata is polymorphic only in populations which grow together with D. ochroleuca, with exception a few examples. Thus, we conclude that variability of the GDH locus in D. incarnata is associated with hybridization with D. ochroleuca.     

Assessment of spatial–temporal variation in natural populations of <Emphasis Type="Italic">Brassica incana</Emphasis> in south Italy: implications for conservation

Brassica incana is a secondary-gene pool wild relative of Brassica oleracea. Twenty-two B. incana populations are recorded in Italy, where the species has recently been pointed out as in priority need of conservation. While data on the spatial and temporal variation of B. incana Italian populations are completely lacking, this information is useful in conservation planning for this species. Three populations from the Sorrento peninsula and from the islands of Ischia and Capri, collected in 1984 and 2012, were characterised for 12 morpho-phenological and 21 genetic traits to assess their spatial–temporal variation. The populations were quite different for morpho-phenological and genetic traits. Spatial differentiation was high and easily explained by the isolation. Temporal differentiation between the Sorrento and Ischia populations was high and explained by a reduction in the population census across time, while it was not significant between the two Capri accessions as such pointing to a major effect of genetic drift. Numerical dimension is extremely relevant in evaluating conservation priorities since it has a major impact on population dynamics over time. The Sorrento and Ischia populations are under threat and urgently need conservation actions, suggesting an alarming scenario for the survival of other crop wild relative populations which are similar in census. Our data also show that, in an allogamous and self-incompatible species like B. incana, populations of 100–200 individuals maintain high allelic diversity. According to obtained results, natural populations of species with similar reproductive system and census can be considered at low risk of genetic erosion.     

Two new species of linebellies (Genus Pseudoscopelus) with a revised key of species (Perciformes: Chiasmodontidae)

The swallowerfishes previously identified as Pseudoscopelus altipinnis and P. cf. altipinnis (Prokofiev, Kukuev, 2005) are described as a complex of two new, closely-related, allopatric species having antitropical distribution. P. astronesthidens—superspecies is most closely related to P. altipinnis due to its dark orobranchial cavity and the details of its photophore arrangement, but is unique within the genus in the disconnection between the mxf and apf series of photophores and in the details of jaws dentition. P. astronesthidens sp. n. (North Atlantic) differs from P. australis sp. n. (Southern Hemisphere) mainly in the long upper jaw and pectoral fins and in the presence of the ppf series of photophores. P. altipinnis Parr, 1933 is considered as a senior synonym of P. microps Fowler, 1934. A revised key for identification of the linebellies of the genus Pseudoscopelus is given.     

Surviving in isolation: genetic variation,bottlenecks and reproductive strategies in the Canarian endemic <Emphasis Type="Italic">Limonium macrophyllum</Emphasis> (Plumbaginaceae)

Oceanic archipelagos are typically rich in endemic taxa, because they offer ideal conditions for diversification and speciation in isolation. One of the most remarkable evolutionary radiations on the Canary Islands comprises the 16 species included in Limonium subsection Nobiles, all of which are subject to diverse threats, and legally protected. Since many of them are single-island endemics limited to one or a few populations, there exists a risk that a loss of genetic variation might limit their long-term survival. In this study, we used eight newly developed microsatellite markers to characterize the levels of genetic variation and inbreeding in L. macrophyllum, a species endemic to the North-east of Tenerife that belongs to Limonium subsection Nobiles. We detected generally low levels of genetic variation over all populations (H _T = 0.363), and substantial differentiation among populations (F _ST = 0.188; R _ST = 0.186) coupled with a negligible degree of inbreeding (F?=?0.042). Obligate outcrossing may have maintained L. macrophyllum relatively unaffected by inbreeding despite the species’ limited dispersal ability and the genetic bottlenecks likely caused by a prolonged history of grazing. Although several factors still constitute a risk for the conservation of L. macrophyllum, the lack of inbreeding and the recent positive demographic trends observed in the populations of this species are factors that favour its future persistence.     

Reproductive Isolation Among <Emphasis Type="Italic">Drosophila arizonae</Emphasis> from Geographically Isolated Regions of North America

A long-standing goal of speciation research is to describe how reproductive isolating barriers develop, when they arise along the ‘speciation continuum’, and to measure the strength with which they restrict gene flow. Drosophila arizonae and D. mojavensis are a recently diverged sister species pair distributed from the southwestern United States through southern Mexico. While incipient speciation in D. mojavensis has been studied for decades, relatively little attention has been directed toward D. arizonae, despite the fact that previous studies have revealed evidence for significant genetic differentiation among populations separated by geographic barriers. Here, we examine the potential for both pre- and post-mating reproductive isolation in D. arizonae from geographically isolated parts of North America. We find evidence for strong premating isolation between flies from northern mainland Mexico and southern mainland Mexico, but no evidence for postmating isolation in any cross. This study highlights the utility of the D. arizonae system for further investigation into the early evolution of premating isolation, and reinforces the potential of the D. arizonae/D. mojavensis system as a whole for studying the evolution of reproductive isolation at a range of divergence times.     

Genetic diversity and population structure of <Emphasis Type="Italic">Taxus cuspidata</Emphasis> Sieb. et Zucc. ex Endl. (Taxaceae) in Russia according to data of the nucleotide polymorphism of intergenic spacers of the chloroplast genome

The genetic diversity and population structure of Taxus cuspidata Sieb. et Zucc. ex Endl. on the Russian part of its range was studied for the first time on the basis of the nucleotide polymorphism of intergenic spacers psbA–trnH, trnL–trnF, and trnS–trnfM of chloroplast DNA. A high level of gene (h = 0.807) and nucleotide (π = 0.0227) diversity was revealed. The data of AMOVA showed that the interpopulation component accounted for 12% of genetic variability (F _ST = 0.12044, P = 0.0000). We revealed 15 haplotypes, four of which were shown to be unique. The presence of common haplotypes in the majority of populations, the absence of phylogenetic structure, and low values or even the absence of nucleotide divergence show that the T. cuspidata populations studied are fragments of the once common ancestor population. Geographical isolation, which resulted from climatic changes in the Pleistocene–Holocene, as well as from human activities, did not produce a significant effect on the genetic structure of the species.     

Genetic relationship and biological status of the industrially important yeast <Emphasis Type="Italic">Saccharomyces eubayanus</Emphasis> Sampaio et al.

The genomes of the recently discovered yeast Saccharomyces eubayanus and traditional S. cerevisiae are known to be found in the yeast S. pastorianus (syn. S. carlsbergensis), which are essential for brewing. The cryotolerant yeast S. bayanus var. uvarum is of great importance for production of some wines. Based on ascospore viability and meiotic recombination of the control parental markers in hybrids, we have shown that there is no complete interspecies post-zygotic isolation between the yeasts S. eubayanus, S. bayanus var. bayanus and S. bayanus var. uvarum. The genetic data presented indicate that all of the three taxa belong to the same species.     

Population genetic analysis of Euro-Arctic polar cod <Emphasis Type="Italic">Boreogadus saida</Emphasis> suggests fjord and oceanic structuring

Polar cod, Boreogadus saida, is a key species in Arctic marine ecosystems; however, its genetic population structure is largely undescribed. The population genetic structure of 472 B. saida specimens among nine locations in the north-east Atlantic was revealed using 12 microsatellite loci. Pairwise F _ST comparisons showed significant population differentiation between B. saida sampled inside fjords in Svalbard and north-east Greenland, as compared to B. saida from the shelf. The observed genetic variation was not a function of isolation by distance, and it is speculated that B. saida populations inhabiting fjords may have become reproductively isolated from shelf-dwelling B. saida during the last post-glacial recolonization.     

<Emphasis Type="Italic">Pseudoscopelus sagamianus</Emphasis> Tanaka, 1908 (Chiasmodontidae): Designation of the neotype,morphological characteristics of the North Pacific population,and species variation within its whole area

North Pacific Pseudoscopelus sagamianus (from the Kuroshio zone) are characterized morphologically in detail. The neotype is established. The North Pacific population and the group of populations from the Atlantic and Indian oceans are compared in detail, significant differences are found between them in coloration of the orobranchial cavity and in the pattern of photophores of series trf and prcf. In the author’s opinion, these differences are of the subspecies level. Potential areas of subspecies are discussed, and the area of the whole species is verified. A possibility of validity of the name P. pierbartus Spitz et al., 2007, for the group of populations from the Atlantic and Indian oceans is considered. Additions and changes to methods of investigations of Pseudoscopelus suggested by Melo et al. (2007) are critically discussed.     

No evidence of intrinsic reproductive isolation between two reciprocally non-monophyletic,ecologically differentiated mountain plants at an early stage of speciation

Adaptation to dissimilar habitats can trigger phenotypic and genetic differences between populations, which may, in the absence of gene flow, ultimately lead to ecological speciation. Reproductive isolation of diverging populations is a critical step at the onset of speciation. An excellent example for exploring the extent of reproductive isolation at early stages of speciation is provided by Heliosperma pusillum and H. veselskyi (Caryophyllaceae), two reciprocally non-monophyletic, ecologically differentiated species from the Alps. Interspecific gene flow—as revealed by recent genetic studies—is rare even between geographically close populations. Cross pollinations and fitness experiments revealed no evidence of intrinsic reproductive barriers, since fitness parameters measured under uniform conditions were not lower in inter- than in intraspecific crosses. Further, morphometric analyses of the offspring clearly showed that the differentiation of parental species is heritable. As parental phenotypes are likely adaptive, the intermediate morphology of hybrids may lead to reduced hybrid fitness in parental habitats. Altogether, H. pusillum and H. veselskyi provide an increasingly well characterised model system offering exciting insights into early stages of ecological speciation.     

Non-congruent geographic patterns of genetic divergence across European species of <Emphasis Type="Italic">Branchinecta</Emphasis> (Anostraca: Branchinectidae)

Three of the five European species of Branchinecta have a disjunct distribution. In this study, we analyze populations of B. ferox and B. orientalis for mitochondrial (cox1) and nuclear (ITS2) molecular markers. We compare intraspecific genetic divergences between geographically distant populations of B. orientalis, from its only known Spanish population (originally described as B. cervantesi) and from a Hungarian population (assigned to B. orientalis since its discovery), with data from two relatively close Iberian populations of B. ferox. Results indicate that isolation between B. ferox and B. orientalis clades is ancient, and that the clade including the two Iberian populations of B. ferox is geographically structured. Conversely, Iberian and Hungarian populations of B. orientalis do not show geographical structure for the mitochondrial fragment. Lack of geographic structure coupled with very low genetic distances indicates that current Iberian and Hungarian populations of B. orientalis originated from a common population stock, and that the time elapsed since their separation has not been long enough to render the clades reciprocally monophyletic. We hypothesize that colonization of the Iberian Peninsula by B. orientalis is probably the consequence of a single recent dispersal event, and consequently we confirm the synonymy between B. cervantesi and B. orientalis.     

Genetic diversity of the house mouse <Emphasis Type="Italic">Mus musculus</Emphasis> and geographic distribution of its subspecies-specific RAPD markers on the territory of Russia

Genetic diversity and geographic distribution of taxon-specific RAPD markers was examined in ten local populations of the house mouse Mus musculus (n = 42). The house mice were generally characterized by moderate genetic variation: polymorphism P ₉₉ = 60%, P ₉₅ = 32.57%; heterozygosity H = 0.12; the observed allele number n _a = 1.6; the effective allele number n _e = 1.18; the within-population differentiation ?_s = 0.388; and Shannon index I = 0.19. The degree of genetic isolation of individual local populations was greatly variable. The genetic subdivision index G _st varied from 0.162 to 0.770 at the gene flow of Nm = 2.58?0.149, while the among-population distances D _N varied from 0.026 to 0.178. The largest part of the genetic diversity was found among the populations (H _T = 0.125), while the within-population diversity was twice lower (H _S = 0.06). The samples examined were well discriminated relative to the sets of RAPD markers. The character distribution pattern provided conditional subdivision of the mice into the “western” and the “eastern” groups with the putative boarder along the Baikal Lake. The first group was characterized by the prevalence of the markers typical of M. m. musculus and M. m. domesticus. The second group was characterized by the prevalence of the markers typical of M. m. musculus, M. m. gansuensis, M. m. castaneus, M. m. domesticus, and M. m. wagneri. The genotype of the nominative subspecies M. m. musculus was background for all populations. In the populations examined some of earlier described subspecies-specific molecular markers were found at different frequencies, pointing to the involvement of several subspecies of M. musculus in the process of hybridization.     

Evaluation and comparison of the genetic structure of <Emphasis Type="Italic">Bunias orientalis</Emphasis> populations in their native range and two non-native ranges

We studied the invasive warty cabbage Bunias orientalis (Brassicaceae) in three geographically distinct areas. Using inter-simple sequence repeat fingerprinting, we analyzed warty cabbages, including non-native populations, from the eastern Baltic and western Siberian regions and native populations from southwestern Russia. The eastern Baltic region and western Siberia represent the two opposite directions of B. orientalis spread in climatically different zones. The genetic structures of the native and non-native B. orientalis populations were assessed through analysis of molecular variance (AMOVA) and the Bayesian clustering method and by determining the main measures of genetic diversity. AMOVA revealed considerable population differentiation in both the native and invasive ranges. Our results did not indicate a decrease in genetic diversity in the non-native populations of B. orientalis. Similar measures of genetic diversity and genetic structure were determined in the invasive populations in two geographically and ecologically distinct, non-native regions located in Europe and Asia. In both of these regions, higher genetic diversity was detected in the non-native populations than in the native region populations, which may be due to multiple introductions. However, Bayesian clustering analysis revealed slightly different sources of invasive populations in the two non-native regions. Genetic diversity patterns revealed the lack of isolation by distance between populations and confirmed the influence of anthropogenic factors on the spread of B. orientalis. The significance of native populations as germplasm resources for breeding is discussed.     

Existence of two strains of <Emphasis Type="Italic">Habrobracon hebetor</Emphasis> (Hymenoptera: Braconidae): a complex in Thailand and Japan

Habrobracon hebetor Say (Hymenoptera: Braconidae) is a cosmopolitan gregarious ectoparasitoid that attacks larvae of several species of Lepidoptera. Although there are two genetically different strains within H. hebetor, distribution of the strains has been poorly understood. In 2010, in Thailand, where H. hebetor has been known as a parasitoid of stored grain pests, it was found that H. hebetor attacked Opisina arenosella Walker (Lepidoptera: Oecophoridae), which is an invasive pest of coconut palm. For correct identification of this H. hebetor, we conducted DNA analysis and cross tests using populations collected from O. arenosella and stored grain pests in Thailand and populations in Japan known as H. hebetor. We obtained 413 bp of mitochondrial cytochrome oxidase I (COI) sequences and 414 bp of 16S rRNA gene sequences, and both indicated that there are two distinct clades within H. hebetor: one contains insects from Thailand, Spain, India, and Barbados; the other contains insects from Japan and the USA. There were no genetic differences or sexual isolation between Thai populations from different hosts. Our results also showed that populations in Thailand were sexually isolated from a H. hebetor population in Japan.     

Genetic differentiation among natural populations of the lizard complex <Emphasis Type="Italic">Darevskia raddei</Emphasis> as inferred from genome microsatellite marking

The article presents the genetic parameters of the populations of lizards of the Darevskia raddei complex (D. raddei nairensis and D. raddei raddei) and the populations of D. valentini calculated on the basis of the analysis of variability of 50 allelic variants of the three nuclear genome microsatellite-containing loci of 83 individuals. It was demonstrated that the F_st genetic distances between the populations of D. raddei nairensis and D. raddei raddei were not statistically significantly different from the F_st genetic distances between the populations of different species, D. raddei and D. valentini. At the same time, these distances were statistically significantly higher than the F_st distances between the populations belonging to one species within the genus Darevskia. These data suggest deep divergence between the populations of D. raddei raddei and D. raddei nairensis of the D. raddei complex and there arises the question on considering them as separate species.     

Molecular phylogeography and paleodistribution modeling of the boreal tree species <Emphasis Type="Italic">Ulmus lamellosa</Emphasis> (T.Wang et S. L. Chang) (Ulmaceae) in China

The uplift of mountains and climatic oscillations are important for understanding of the demographic history and genetic structure of species. We investigated the biogeographic history of the boreal tree species Ulmus lamellosa (Ulmaceae) in China, by using a combined phylogeographic and paleodistribution modeling approach. In this study, 14 populations of endangered U. lamellosa were analyzed by using chloroplast DNA (cpDNA) sequences. A high level of genetic differentiation (Φ _ST = 86.22%) among populations with a significant phylogeographic pattern (N _ST > G _ST, P < 0.05) was found in U. lamellosa. Ten haplotypes were detected by combining chloroplast DNA data, and haplotype 3 (H3) was found to be common and widespread. The intraspecific divergence of all U. lamellosa cpDNA haplotypes (9.27 Ma; 95% HPD 5.17–13.33 Ma) most probably began in the late Miocene. The pairwise difference among haplotypes and neutrality tests (Tajima’s D and Fu’s Fs statistic) indicated that populations of U. lamellosa, except group I, have not experienced recent sudden expansions. Multiple refuge areas were identified across the entire distribution ranges of U. lamellosa. The low level of gene flow (Nm = 0.14) among populations may have resulted from isolation resulting from distance and complex topography during climatic oscillations; this isolation was probably the major process that shaped the present distribution of haplotypes. These results support the hypothesis that U. lamellosa persisted in situ during glaciations and occupied multiple localized glacial refugia, contrary to the hypotheses of large-scale range contraction and long-distance southward migration.     

Are microbial symbionts involved in the speciation of the gall-inducing aphid, <Emphasis Type="Italic">Slavum wertheimae</Emphasis>?

Microbial symbionts have come to be recognized as agents in the speciation of their eukaryote hosts. In this study, we asked if bacterial symbionts are, or were in the past, involved in the speciation of the gall-inducing aphid Slavum wertheimae (Hemiptera: Aphididae). This aphid is specific to the tree Pistacia atlantica, which has a fragmented distribution among mesic and xeric habitats, leading to corresponding fragmentation of the aphid population. Previous studies revealed genetic differentiation among populations of the gall-inducing aphid, suggesting cryptic allopatric speciation. Pistacia atlantica trees show no such variation. By means of diagnostic PCR, we screened several populations of S. wertheimae from mesic and xeric sites in Israel for the presence of nine known aphid symbionts: Arsenophonus, Hamiltonella, Regiella, Rickettsia, Rickettsiella, Serratia, Spiroplasma, Wolbachia, and X-type, as well as Cardinium, known to be a reproductive manipulator. Only one symbiont, Wolbachia, was detected in S. wertheimae. Wolbachia was found in all the aphids of the mesic populations, compared to 26% in the aphids from the xeric populations. Multilocus Sequence typing of Wolbachia revealed new haplotypes in the fbpA and coxA genes in both the mesic and xeric populations. Phylogenetic analysis showed that Wolbachia of S. wertheimae is closely related to Wolbachia strains from assorted hosts, mostly lepidopterans, but only distantly related to Wolbachia strains from other aphid species. We conclude that the cryptic speciation of mesic and xeric populations of S. wertheimae was likely driven by geographical isolation rather than by Wolbachia.     

Genetic structure of populations and natural hybridization between <Emphasis Type="Italic">Dactylorhiza salina</Emphasis> and <Emphasis Type="Italic">D. incarnata</Emphasis> (Orchidaceae)

The results of studying the polymorphism and genetic structure of populations of D. salina and D. incarnata growing in Zabaykalsky krai and Buryatia are represented according to the data of allozyme analysis of eight genetic loci (PGI, NADHD, SKDH, GDH, PGM, DIA, ADH, and IDH). The specificity of the allelic structure of loci SKDH, PGM, and IDH is established, for which D. salina and D. incarnata reliably differ from each other. It is shown that interspecies introgressive hybrid complexes with different genetic structures were formed in Transbaikalia. Places of mass growth of D. incarnata were observed to have single plants of D. salina, the interspecies hybrids of the first and subsequent generations. Places of mass growth of D. salina were observed to contain only the hybrids that are not hybrids of the first generation. They were heterozygous not for three loci with differentiating alleles of both parents, SKDH, PGM, and IDH, but for only one of them. The degree of genetic differentiation among five populations of D. salina was on average 7.5% and that of D. incarnata was 7.1%, which in accordance with Wright’s estimation relates to mean values. The average value of FST for all studied populations of the two related species of the genus Dactylorhiza was 0.478, indicating a very high degree of genetic differentiation between D. salina and D. incarnata growing in Transbaikalia. The greatest differences between the species are for the allelic structure of loci SKDH, PGM, and IDH (FST was equal to 0.705, 0.976, and 0.762, respectively). Analysis of molecular variance (AMOVA) showed that populations of D. salina and D. incarnata in the zone where their ranges in Zabaykalsky krai and Buryatya overlap have essential differences both for the variation of alleles frequencies of eight loci (71%, d.f. = 9) and for the variability of genotypes (61%, d.f. = 9). Despite the fact that D. salina and D. incarnata explicitly share a gene flow as a result of interspecies hybridization, the genetic differentiation of populations of these related species remains at a high level.     

Pelagic larval dispersal habits influence the population genetic structure of clam <Emphasis Type="Italic">Gomphina aequilatera</Emphasis> in China

Pelagic larval dispersal habits influence the population genetic structure of marine mollusk organisms via gene flow. The genetic information of the clam Gomphina aequilatera (short larval stage, 10 days) which is ecologically and economically important in the China coast is unknown. To determine the influence of planktonic larval duration on the genetic structure of G. aequilatera. Mitochondrial markers, cytochrome oxidase subunit i (COI) and 12S ribosomal RNA (12S rRNA), were used to investigate the population structure of wild G. aequilatera specimens from four China Sea coastal locations (Zhoushan, Nanji Island, Zhangpu and Beihai). Partial COI (685 bp) and 12S rRNA (350 bp) sequences were determined. High level and significant F_ST values were obtained among the different localities, based on either COI (F_ST?=?0.100–0.444, P?<?0.05) or 12S rRNA (F_ST?=?0.193–0.742, P?<?0.05), indicating a high degree of genetic differentiation among the populations. The pairwise N_m between Beihai and Zhoushan for COI was 0.626 and the other four pairwise N_m values were >?1, indicating extensive gene flow among them. The 12S rRNA showed the same pattern. AMOVA test results for COI and 12S rRNA indicated major genetic variation within the populations: 77.96% within and 22.04% among the populations for COI, 55.73% within and 44.27% among the populations for 12S rRNA. A median-joining network suggested obvious genetic differentiation between the Zhoushan and Beihai populations. This study revealed the extant population genetic structure of G. aequilatera and showed a strong population structure in a species with a short planktonic larval stage.