Genetic structure and demographic history of Lymantria dispar (Linnaeus, 1758) (Lepidoptera: Erebidae) in its area of origin and adjacent areas

We analyzed the population genetic structure and demographic history of 20 Lymantria dispar populations from Far East Asia using microsatellite loci and mitochondrial genes. In the microsatellite analysis, the genetic distances based on pairwise F_ST values ranged from 0.0087 to 0.1171. A NeighborNet network based on pairwise F_ST genetic distances showed that the 20 regional populations were divided into five groups. Bayesian clustering analysis (K = 3) demonstrated the same groupings. The populations in the Korean Peninsula and adjacent regions, in particular, showed a mixed genetic pattern. In the mitochondrial genetic analysis based on 98 haplotypes, the median‐joining network exhibited a star shape that was focused on three high‐frequency haplotypes (Haplotype 1: central Korea and adjacent regions, Group 1; Haplotype 37: southern Korea, Group 2; and Haplotype 90: Hokkaido area, Group 3) connected by low‐frequency haplotypes. The mismatch distribution dividing the three groups was unimodal. In the neutral test, Tajima's D and Fu's FS tests were negative. We can thus infer that the Far East Asian populations of L. dispar underwent a sudden population expansion. Based on the age expansion parameter, the expansion time was inferred to be approximately 53,652 years before present (ybp) for Group 1, approximately 65,043 ybp for Group 2, and approximately 76,086 ybp for Group 3. We propose that the mixed genetic pattern of the inland populations of Far East Asia is due to these expansions and that the inland populations of the region should be treated as valid subspecies that are distinguishable from other subspecies by genetic traits.     

Intraspecific genetic variation in the common midwife toad (Alytes obstetricans): subspecies assignment using mitochondrial and microsatellite markers

The common midwife toad (Alytes obstetricans), widely distributed in the northern half of the Iberian Peninsula and part of Western Europe, is currently subdivided into four subspecies: A. o. obstetricans, A. o. boscai, A. o. pertinax and A. o. almogavarii. However, the delimitation of these subspecies and their ranges are still under discussion because strong discordances have been found between morphological and molecular data, and especially among different genetic markers. Here, we screen a set of novel microsatellite loci and mtDNA sequences of A. obstetricans populations representative of all currently recognized subspecies to investigate the correspondence between genetic groupings inferred from clustering analysis of microsatellite genotypes and the described subspecies and test whether patterns of mtDNA variation are concordant with those genetic clusters. Our results confirm previous expectations of extremely high intraspecific diversity in A. obstetricans in Iberia. Analyses of microsatellite and mtDNA data were concordant in recovering five well‐defined groups, of which three correspond to previously defined subspecies, while the two additional clusters correspond to populations of subspecies A. o. boscai separated by the Douro River. Our results suggest the occurrence of two distinct genetic units within A. o. boscai that likely result from a long independent evolutionary history, thus deserving special attention from a conservation point of view.     

Genetic diversity loss in a biodiversity hotspot: ancient DNA quantifies genetic decline and former connectivity in a critically endangered marsupial

The extent of genetic diversity loss and former connectivity between fragmented populations are often unknown factors when studying endangered species. While genetic techniques are commonly applied in extant populations to assess temporal and spatial demographic changes, it is no substitute for directly measuring past diversity using ancient DNA (aDNA). We analysed both mitochondrial DNA (mtDNA) and nuclear microsatellite loci from 64 historical fossil and skin samples of the critically endangered Western Australian woylie (Bettongia penicillata ogilbyi), and compared them with 231 (n = 152 for mtDNA) modern samples. In modern woylie populations 15 mitochondrial control region (CR) haplotypes were identified. Interestingly, mtDNA CR data from only 29 historical samples demonstrated 15 previously unknown haplotypes and detected an extinct divergent clade. Through modelling, we estimated the loss of CR mtDNA diversity to be between 46% and 91% and estimated this to have occurred in the past 2000–4000 years in association with a dramatic population decline. In addition, we obtained near‐complete 11‐loci microsatellite profiles from 21 historical samples. In agreement with the mtDNA data, a number of ‘new’ microsatellite alleles was only detected in the historical populations despite extensive modern sampling, indicating a nuclear genetic diversity loss >20%. Calculations of genetic diversity (heterozygosity and allelic rarefaction) showed that these were significantly higher in the past and that there was a high degree of gene flow across the woylie's historical range. These findings have an immediate impact on how the extant populations are managed and we recommend the implementation of an assisted migration programme to prevent further loss of genetic diversity. Our study demonstrates the value of integrating aDNA data into current‐day conservation strategies.     

Genetic differentiation of the pine processionary moth at the southern edge of its range: contrasting patterns between mitochondrial and nuclear markers

The pine processionary moth (Thaumetopoea pityocampa) is an important pest of coniferous forests at the southern edge of its range in Maghreb. Based on mitochondrial markers, a strong genetic differentiation was previously found in this species between western (pityocampa clade) and eastern Maghreb populations (ENA clade), with the contact zone between the clades located in Algeria. We focused on the moth range in Algeria, using both mitochondrial (a 648 bp fragment of the tRNA‐cox2) and nuclear (11 microsatellite loci) markers. A further analysis using a shorter mtDNA fragment and the same microsatellite loci was carried out on a transect in the contact zone between the mitochondrial clades. Mitochondrial diversity showed a strong geographical structure and a well‐defined contact zone between the two clades. In particular, in the pityocampa clade, two inner subclades were found whereas ENA did not show any further structure. Microsatellite analysis outlined a different pattern of differentiation, with two main groups not overlapping with the mitochondrial clades. The inconsistency between mitochondrial and nuclear markers is probably explained by sex‐biased dispersal and recent afforestation efforts that have bridged isolated populations.     

Weak large‐scale population genetic structure in a philopatric seabird,the European Shag Phalacrocorax aristotelis

Quantifying population genetic structure is fundamental to testing hypotheses regarding gene flow, population divergence and dynamics across large spatial scales. In species with highly mobile life‐history stages, where it is unclear whether such movements translate into effective dispersal among discrete philopatric breeding populations, this approach can be particularly effective. We used seven nuclear microsatellite loci and mitochondrial DNA (ND2) markers to quantify population genetic structure and variation across 20 populations (447 individuals) of one such species, the European Shag, spanning a large geographical range. Despite high breeding philopatry, rare cross‐sea movements and recognized subspecies, population genetic structure was weak across both microsatellites and mitochondrial markers. Furthermore, although isolation‐by‐distance was detected, microsatellite variation provided no evidence that open sea formed a complete barrier to effective dispersal. These data suggest that occasional long‐distance, cross‐sea movements translate into gene flow across a large spatial scale. Historical factors may also have shaped contemporary genetic structure: cluster analyses of microsatellite data identified three groups, comprising colonies at southern, mid‐ and northern latitudes, and similar structure was observed at mitochondrial loci. Only one private mitochondrial haplotype was found among subspecies, suggesting that this current taxonomic subdivision may not be mirrored by genetic isolation.     

Genetic isolation in an endemic African habitat specialist

The Chestnut‐banded Plover Charadrius pallidus is a Near‐Threatened shorebird species endemic to mainland Africa. We examined levels of genetic differentiation between its two morphologically and geographically distinct subspecies, C. p. pallidus in southern Africa (population size 11 000–16 000) and C. p. venustus in eastern Africa (population size 6500). In contrast to other plover species that maintain genetic connectivity over thousands of kilometres across continental Africa, we found profound genetic differences between remote sampling sites. Phylogenetic network analysis based on four nuclear and two mitochondrial gene regions, and population genetic structure analyses based on 11 microsatellite loci, indicated strong genetic divergence, with 2.36% mitochondrial sequence divergence between individuals sampled in Namibia (southern Africa) and those of Kenya and Tanzania (eastern Africa). This distinction between southern and eastern African populations was also supported by highly distinct genetic clusters based on microsatellite markers (global F_ST = 0.309,  = 0.510, D = 0.182). Behavioural factors that may promote genetic differentiation in this species include habitat specialization, monogamous mating behaviour and sedentariness. Reliance on an extremely small number of saline lakes for breeding and limited dispersal between populations are likely to promote reproductive and genetic isolation between eastern and southern Africa. We suggest that the two Chestnut‐banded Plover subspecies may warrant elevation to full species status. To assess this distinction fully, additional sample collection will be needed, with analysis of genetic and phenotypic traits from across the species’ entire breeding range.     

Global genetic diversity of Aedes aegypti

Mosquitoes, especially Aedes aegypti, are becoming important models for studying invasion biology. We characterized genetic variation at 12 microsatellite loci in 79 populations of Ae. aegypti from 30 countries in six continents, and used them to infer historical and modern patterns of invasion. Our results support the two subspecies Ae. aegypti formosus and Ae. aegypti aegypti as genetically distinct units. Ae. aegypti aegypti populations outside Africa are derived from ancestral African populations and are monophyletic. The two subspecies co‐occur in both East Africa (Kenya) and West Africa (Senegal). In rural/forest settings (Rabai District of Kenya), the two subspecies remain genetically distinct, whereas in urban settings, they introgress freely. Populations outside Africa are highly genetically structured likely due to a combination of recent founder effects, discrete discontinuous habitats and low migration rates. Ancestral populations in sub‐Saharan Africa are less genetically structured, as are the populations in Asia. Introduction of Ae. aegypti to the New World coinciding with trans‐Atlantic shipping in the 16th to 18th centuries was followed by its introduction to Asia in the late 19th century from the New World or from now extinct populations in the Mediterranean Basin. Aedes mascarensis is a genetically distinct sister species to Ae. aegypti s.l. This study provides a reference database of genetic diversity that can be used to determine the likely origin of new introductions that occur regularly for this invasive species. The genetic uniqueness of many populations and regions has important implications for attempts to control Ae. aegypti, especially for the methods using genetic modification of populations.     

Genetic structure of the critically endangered Red‐headed Wood Pigeon Columba janthina nitens and its implications for the management of threatened island populations

The Red‐headed Wood Pigeon Columba janthina nitens is endemic to the Ogasawara Islands, an oceanic island chain located 1000 km south of the main islands of Japan. The subspecies is at high risk of extinction because of its small population size and restricted habitat range. We undertook genetic analyses of this pigeon using sequences of a portion of the mitochondrial control region and five microsatellite markers to estimate the genetic characteristics of two wild populations from the Bonin and Volcano Islands, as well as one captive breeding population. The genetic diversity of the wild individuals was exceptionally low in both the mitochondria (nucleotide diversity = 0.00105) and at the microsatellite (3.2 alleles per locus and H_E = 0.12) loci. Higher numbers of microsatellite genotypes were observed in the Volcano Islands population than in the Bonin Islands population, which may be because of the relatively low impact of human disturbance. The most common mitochondrial haplotypes and microsatellite alleles observed in the two wild populations were completely fixed in the captive population. Our results suggest that the genetic diversity of the captive population needs to be increased. However, introduction of a wild individual into a captive population can lead to a decreased genetic diversity in the wild population and therefore should be done with caution. The genetic differentiation between the Bonin and the Volcano island groups was low, and the populations of the two island groups should be regarded as a single evolutionarily significant unit. However, special consideration is required for habitat conservation in the Volcano Islands, which may be functioning as a sanctuary for the Red‐headed Wood Pigeon. For the long‐term conservation of threatened bird species that live on remote oceanic islands, determination of management units considering gene flow caused by their flying capacity and maintenance of genetically suitable wild and captive populations are essential.     

Extremely low nucleotide polymorphism in Pinus krempfii Lecomte,a unique flat needle pine endemic to Vietnam

Pinus krempfii Lecomte is a morphologically and ecologically unique pine, endemic to Vietnam. It is regarded as vulnerable species with distribution limited to just two provinces: Khanh Hoa and Lam Dong. Although a few phylogenetic studies have included this species, almost nothing is known about its genetic features. In particular, there are no studies addressing the levels and patterns of genetic variation in natural populations of P. krempfii. In this study, we sampled 57 individuals from six natural populations of P. krempfii and analyzed their sequence variation in ten nuclear gene regions (approximately 9 kb) and 14 mitochondrial (mt) DNA regions (approximately 10 kb). We also analyzed variation at seven chloroplast (cp) microsatellite (SSR) loci. We found very low haplotype and nucleotide diversity at nuclear loci compared with other pine species. Furthermore, all investigated populations were monomorphic across all mitochondrial DNA (mtDNA) regions included in our study, which are polymorphic in other pine species. Population differentiation at nuclear loci was low (5.2%) but significant. However, structure analysis of nuclear loci did not detect genetically differentiated groups of populations. Approximate Bayesian computation (ABC) using nuclear sequence data and mismatch distribution analysis for cpSSR loci suggested recent expansion of the species. The implications of these findings for the management and conservation of P. krempfii genetic resources were discussed.     

Population genetic structure and taxonomy of the common dolphin (Delphinus sp.) at its southernmost range limit: New Zealand waters

New Zealand is the southernmost limit of the common dolphin's (genus Delphinus) distribution in the Pacific Ocean. In this area, common dolphins occur in both coastal and oceanic habitats, exhibit seasonal and resident occurrence, and present high morphological variability. Here we investigated the population structure and the taxonomic identity of common dolphins (Delphinus sp.) within New Zealand waters using 14 microsatellite loci, 577 bp of the mtDNA control region, and 1,120 bp of the mtDNA cytochrome b gene across 90 individuals. We found high genetic variability and evidence of population expansion. Phylogenetic analyses conducted to clarify the taxonomic status of New Zealand common dolphins did not show any clustering reflecting geographic origin or morphotypes. The microsatellite analysis showed genetic differentiation between Coastal and Oceanic putative populations, while mtDNA revealed significant genetic differentiation only between the Hauraki Gulf and other putative groups. Our results suggest that differences in habitat choice and possible female site fidelity may play a role in shaping population structure of New Zealand common dolphins.     

Genetic variation in Tertiary relics: The case of eastern‐Mediterranean Abies (Pinaceae)

The eastern‐Mediterranean Abies taxa, which include both widely distributed species and taxa with minuscule ranges, represent a good model to study the impacts of range size and fragmentation on the levels of genetic diversity and differentiation. To assess the patterns of genetic diversity and phylogenetic relationships among eastern‐Mediterranean Abies taxa, genetic variation was assessed by eight nuclear microsatellite loci in 52 populations of Abies taxa with a focus on those distributed in Turkey and the Caucasus. Both at the population and the taxon level, the subspecies or regional populations of Abies nordmanniana s.l. exhibited generally higher allelic richness, private allelic richness, and expected heterozygosity compared with Abies cilicica s.l. Results of both the Structure analysis and distance‐based approaches showed a strong differentiation of the two A. cilicica subspecies from the rest as well as from each other, whereas the subspecies of A. nordmanniana were distinct but less differentiated. ABC simulations were run for a set of scenarios of phylogeny and past demographic changes. For A. ×olcayana, the simulation gave a poor support for the hypothesis of being a taxon resulting from a past hybridization, the same is true for Abies equi‐trojani: both they represent evolutionary branches of Abies bornmuelleriana.     

Genetic differentiation and intraspecific structure of Eastern Tropical Pacific spotted dolphins,Stenella attenuata,revealed by DNA analyses

Mitochondrial DNA (mtDNA) control region sequences and microsatellite loci length polymorphisms were used to investigate genetic differentiation in spotted dolphins (Stenella attenuata) in the Eastern Tropical Pacific and to examine the intraspecific structure of the coastal subspecies (Stenella attenuata graffmani). One-hundred and thirty-five animals from several coastal areas and 90 offshore animals were sequenced for 455 bp of the mitochondrial control region, resulting in 112 mtDNA haplotypes. Phylogenetic analyses and the existence of shared haplotypes between the two subspecies suggest recent and/or current gene flow. Analyses using χ², F _ST (based on haplotype frequencies) and Φ_ST values (based on frequencies and genetic distances between haplotypes) yielded statistically significant separation (randomized permutation values P<0.05) among four different coastal populations and between all but one of these and the offshore subspecies (overall F _ST=0.0691). Ninety-one coastal animals from these four geographic populations and 50 offshore animals were genotyped for seven nuclear microsatellite loci. Analysis using F _ST values (based on allelic frequencies) yielded statistically significant separation between most coastal populations and offshore animals, although no coastal populations were distinguished. These results argue for the existence of some genetic isolation between offshore and inshore populations and among some inshore populations, suggesting that these should be treated as separate units for management purposes.     

Microsatellites from the vairone Leuciscus souffia (Pisces: Cyprinidae) and their application to closely related species

Eleven microsatellites were isolated from the vairone Leuciscus souffia (Risso 1826), an endangered fish that inhabits river systems in and around the Alps in Europe. The level of genetic diversity was assessed in 29 individuals of the subspecies L. s. souffia, and their variability was further estimated in seven individuals of a different subspecies, L. s. muticellus. Eight of these microsatellite loci were also applied to seven closely related cyprinid species. Availability of the reported microsatellite loci will facilitate the investigation of population genetic structure of these species with applications for the development of conservation strategies and phylogeographical approaches.     

Population genetics of Wolbachia‐infected,parthenogenetic and uninfected,sexual populations of Tetrastichus coeruleus (Hymenoptera: Eulophidae)

Wolbachia are endosymbiotic bacteria known to manipulate the reproduction of their hosts. These manipulations are expected to have consequences on the population genetics of the host, such as heterozygosity levels, genetic diversity and gene flow. The parasitoid wasp Tetrastichus coeruleus has populations that are infected with parthenogenesis‐inducing Wolbachia and populations that are not infected. We studied the population genetics of T. coeruleus between and within Wolbachia‐infected and uninfected populations, using nuclear microsatellites and mitochondrial DNA. We expected reduced genetic diversity in both DNA types in infected populations. However, migration and gene flow could introduce new DNA variants into populations. We therefore paid special attention to individuals with unexpected (genetic) characteristics. Based on nuclear and mitochondrial DNA, two genetic clusters were evident: a thelytokous cluster containing all Wolbachia‐infected, parthenogenetic populations and an arrhenotokous cluster containing all uninfected, sexual populations. Nuclear and mitochondrial DNA did not exhibit concordant patterns of variation, although there was reduced genetic diversity in infected populations for both DNA types. Within the thelytokous cluster, there was nuclear DNA variation, but no mitochondrial DNA variation. This nuclear DNA variation may be explained by occasional sex between infected females and males, by horizontal transmission of Wolbachia, and/or by novel mutations. Several females from thelytokous populations were uninfected and/or heterozygous for microsatellite loci. These unexpected characteristics may be explained by migration, by inefficient transmission of Wolbachia, by horizontal transmission of Wolbachia, and/or by novel mutations. However, migration has not prevented the build‐up of considerable genetic differentiation between thelytokous and arrhenotokous populations.     

Unexpected cryptic species diversity in the widespread coral Seriatopora hystrix masks spatial‐genetic patterns of connectivity

Mounting evidence of cryptic species in a wide range of taxa highlights the need for careful analyses of population genetic data sets to unravel within‐species diversity from potential interspecies relationships. Here, we use microsatellite loci and hierarchical clustering analysis to investigate cryptic diversity in sympatric and allopatric (separated by 450 km) populations of the widespread coral Seriatopora hystrix on the Great Barrier Reef. Structure analyses delimited unique genetic clusters that were confirmed by phylogenetic and extensive population‐level analyses. Each of four sympatric yet distinct genetic clusters detected within S. hystrix demonstrated greater genetic cohesion across regional scales than between genetic clusters within regions (<10 km). Moreover, the magnitude of genetic differentiation between different clusters (>0.620 G”_ST) was similar to the difference between S. hystrix clusters and the congener S. caliendrum (mean G”_ST 0.720). Multiple lines of evidence, including differences in habitat specificity, mitochondrial identity, Symbiodinium associations and morphology, corroborate the nuclear genetic evidence that these distinct clusters constitute different species. Hierarchical clustering analysis combined with more traditional population genetic methods provides a powerful approach for delimiting species and should be regularly applied to ensure that ecological and evolutionary patterns interpreted for single species are not confounded by the presence of cryptic species.     

Different genetic structures revealed resident populations of a specialist parasitoid wasp in contrast to its migratory host

Genetic comparisons of parasitoids and their hosts are expected to reflect ecological and evolutionary processes that influence the interactions between species. The parasitoid wasp, Cotesia vestalis, and its host diamondback moth (DBM), Plutella xylostella, provide opportunities to test whether the specialist natural enemy migrates seasonally with its host or occurs as resident population. We genotyped 17 microsatellite loci and two mitochondrial genes for 158 female adults of C. vestalis collected from 12 geographical populations, as well as nine microsatellite loci for 127 DBM larvae from six separate sites. The samplings covered both the likely source (southern) and immigrant (northern) areas of DBM from China. Populations of C. vestalis fell into three groups, pointing to isolation in northwestern and southwestern China and strong genetic differentiation of these populations from others in central and eastern China. In contrast, DBM showed much weaker genetic differentiation and high rates of gene flow. TESS analysis identified the immigrant populations of DBM as showing admixture in northern China. Genetic disconnect between C. vestalis and its host suggests that the parasitoid did not migrate yearly with its host but likely consisted of resident populations in places where its host could not survive in winter.     

A genetic demographic analysis of Lake Malawi rock‐dwelling cichlids using spatio‐temporal sampling

We estimated the effective population sizes (N_e) and tested for short‐term temporal demographic stability of populations of two Lake Malawi cichlids: Maylandia benetos, a micro‐endemic, and Maylandia zebra, a widespread species found across the lake. We sampled a total of 351 individuals, genotyped them at 13 microsatellite loci and sequenced their mitochondrial D‐loop to estimate genetic diversity, population structure, demographic history and effective population sizes. At the microsatellite loci, genetic diversity was high in all populations. Yet, genetic diversity was relatively low for the sequence data. Microsatellites yielded mean N_e estimates of 481 individuals (±99 SD) for M. benetos and between 597 (±106.3 SD) and 1524 (±483.9 SD) individuals for local populations of M. zebra. The microsatellite data indicated no deviations from mutation–drift equilibrium. Maylandia zebra was further found to be in migration–drift equilibrium. Temporal fluctuations in allele frequencies were limited across the sampling period for both species. Bayesian Skyline analyses suggested a recent expansion of M. zebra populations in line with lake‐level fluctuations, whereas the demographic history of M. benetos could only be estimated for the very recent past. Divergence time estimates placed the origin of M. benetos within the last 100 ka after the refilling of the lake and suggested that it split off the sympatric M. zebra population. Overall, our data indicate that micro‐endemics and populations in less favourable habitats have smaller N_e, indicating that drift may play an important role driving their divergence. Yet, despite small population sizes, high genetic variation can be maintained.     

Fine‐scale population structure of Collichtys lucidus populations inferred from microsatellite markers

The spinyhead croaker Collichthys lucidus (Richardson) is a small sciaenid species distributed along the inshore waters of northwestern Pacific Ocean, and now has been listed as Key Protected Commercial Sources of Aquatic Animals and Plants in China. To delineate stock boundaries and inform conservation policy for its management, samples were collected from eight locations across the Chinese coastal waters and analyzed at nine microsatellite loci. C. lucidus populations showed low genetic diversity (expected heterozygosity = 0.445–0.542; observed heterozygosity = 0.392–0.539; Polymorphism Information Content = 0.268–0.684). Strong genetic fdifferentiation (F_st = 0.065–0.510, all significant after Bonferroni correction) among all populations and high levels of self‐recruitment (89.2%–91.5%) were observed, which suggested limited genetic exchange for this species. Clustering results of discriminant analysis of principal components and STRUCTURE found strong support for obvious genetic clusters (populations FZ, XM and SZ vs. populations SH, YRE, ZS, WZ and ND). The results of the present study not only supported the phylogeographic pattern of north‐south differentiation, but also suggested that C. lucidus populations may be predominantly sustained by self‐replenishment rather than by recruitment from distant populations.     

Population genetic isolation and limited connectivity in the purple finch (Haemorhous purpureus)

Using a combination of mitochondrial and z‐linked sequences, microsatellite data, and spatio‐geographic modeling, we examined historical and contemporary factors influencing the population genetic structure of the purple finch (Haemorhous purpureus). Mitochondrial DNA data show the presence of two distinct groups corresponding to the two subspecies, H. p. purpureus and H. p. californicus. The two subspecies likely survived in separate refugia during the last glacial maximum, one on the Pacific Coast and one east of the Rocky Mountains, and now remain distinct lineages with little evidence of gene flow between them. Southwestern British Columbia is a notable exception, as subspecies mixing between central British Columbia and Vancouver Island populations suggests a possible contact zone in this region. Z‐linked data support two mitochondrial groups; however, Coastal Oregon and central British Columbia sites show evidence of mixing. Contemporary population structure based on microsatellite data identified at least six genetic clusters: three H. p. purpureus clusters, two H. p. californicus clusters, and one mixed cluster, which likely resulted from high site fidelity and isolation by distance, combined with sexual selection on morphological characters reinforcing subspecies differences.     

Fine‐scale patterns of genetic divergence within and between morphologically variable subspecies of the sea urchin Heliocidaris erythrogramma (Echinometridae)

The spatial scale over which genetic divergences occur between populations and the extent that they are paralleled by morphological differences can vary greatly among marine species. In the present study, we use a hierarchical spatial design to investigate genetic structure in Heliocidaris erythrogramma occurring on near shore limestone reefs in Western Australia. These reefs are inhabited by two distinct subspecies: the thick‐spined Heliocidaris erythrogramma armigera and the thin‐spined Heliocidaris erythrogramma erythrogramma, each of which also have distinct colour patterns. In addition to pronounced morphological variation, H. erythrogramma exhibits a relatively short (3–4 days) planktonic phase before settlement and metamorphosis, which limits their capacity for dispersal. We used microsatellite markers to determine whether patterns of genetic structure were influenced more by morphological or life history limitations to dispersal. Both individual and population‐level analyses found significant genetic differentiation between subspecies, which was independent of geographical distance. Genetic diversity was considerably lower within H. e. erythrogramma than within H. e. armigera and genetic divergence was four‐fold greater between subspecies than among populations within subspecies. This pattern was consistent even at fine spatial scales (< 5 km). We did detect some evidence of gene flow between the subspecies; however, it appears to be highly restricted. Within subspecies, genetic structure was more clearly driven by dispersal capacity, although weak patterns of isolation‐by‐distance suggest that there may be other factors limiting gene exchange between populations. Our results show that spatial patterns of genetic structure in Western Australian H. erythrogramma is influenced by a range of factors but is primarily correlated with the distribution of morphologically distinct subspecies. This suggests the presence of reproductive barriers to gene exchange between them and demonstrates that morphological variation can be a good predictor of genetic divergence. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 103 , 578–592.