High similarity of genetic population structure in the false clown anemonefish (<Emphasis Type="Italic">Amphiprion ocellaris</Emphasis>) found in microsatellite and mitochondrial control region analysis

Many studies, using various marker systems, have been conducted on the genetic population structure of marine organisms to reveal connectivity among locations and dispersal capabilities. Although mitochondrial sequence markers are widely used, their accuracy is controversially discussed in the context of small scale population genetic discrimination. In the present study, the genetic population structure of the False Clown Anemonefish (Amphiprion ocellaris) in the Indo-Malay Archipelago was revealed by screening six microsatellite loci. Results were congruent to previous mitochondrial control region results, with three major genetic breaks within the Indo-Malay Archipelago. Similar to the mitochondrial DNA (mtDNA) analysis, microsatellite data showed a correlation of genetic structure to historical ocean basin separation during Pleistocene sea level low stands, geographic distance, and dominant current patterns. However, microsatellite divergences are not as deep as the mtDNA divergence, suggesting either that admixture of mtDNA lineages is slower than that of nuclear microsatellites, providing a rather historic picture of separation, or the stronger differentiation signal is due to lower effective population sizes presented by mtDNA. As well, the microsatellite analysis did not give a better resolution on the small scale as expected. This study showed that depending on the genetic markers used, different stages of population separation might be illuminated.     

Geographic Structure of Mitochondrial and Nuclear Gene Polymorphisms in Australian Green Turtle Populations and Male-Biased Gene Flow

The genetic structure of green turtle (Chelonia mydas) rookeries located around the Australian coast was assessed by (1) comparing the structure found within and among geographic regions, (2) comparing microsatellite loci vs. restriction fragment length polymorphism analyses of anonymous single copy nuclear DNA (ascnDNA) loci, and (3) comparing the structure found at nuclear DNA markers to that of previously analyzed mitochondrial (mtDNA) control region sequences. Significant genetic structure was observed over all regions at both sets of nuclear markers, though the microsatellite data provided greater resolution in identifying significant genetic differences in pairwise tests between regions. Inferences about population structure and migration rates from the microsatellite data varied depending on whether statistics were based on the stepwise mutation or infinite allele model, with the latter being more congruent with geography. Estimated rates of gene flow were generally higher than expected for nuclear DNA (nDNA) in comparison to mtDNA, and this difference was most pronounced in comparisons between the northern and southern Great Barrier Reef (GBR). The genetic data combined with results from physical tagging studies indicate that the lack of nuclear gene divergence through the GBR is likely due to the migration of sGBR turtles through the courtship area of the nGBR population, rather than male-biased dispersal. This example highlights the value of combining comparative studies of molecular variation with ecological data to infer population processes.     

Genetic evidence for female-biased dispersal and gene flow in a polygynous primate

Many models of sex-biased dispersal predict that the direction of sex-bias depends upon a species' mating system. In agreement with this, almost all polygynous mammals show male-biased dispersal whereas largely monogamous birds show female-biased dispersal (FBD). The hamadryas baboon (Papio hamadryas hamadryas) is polygynous and so dispersal is predicted to be male biased, as is found in all other baboon subspecies, but there are conflicting field data showing both female and male dispersal. Using 19 autosomal genetic markers genotyped in baboons from four Saudi Arabian populations, we found strong evidence for FBD in post-dispersal adults but not, as expected, in pre-dispersal infants and young juveniles, when we compared male and female: population structure (F(st)), inbreeding (F(is)), relatedness (r), and the mean assignment index (mAIc). Furthermore, we found evidence for female-biased gene flow as population genetic structure (F(st)), was about four times higher for the paternally inherited Y, than for either autosomal markers or for maternally inherited mtDNA. These results contradict the direction of sex-bias predicted by the mating system and show that FBD has evolved recently from an ancestral state of male-biased dispersal. We suggest that the cost-benefit balance of dispersal to males and females is tightly linked to the unique hierarchical social structure of hamadryas baboons and that dispersal and social organization have coevolved.     

Phylogeography, intraspecific structure and sex-biased dispersal of Dall's porpoise, Phocoenoides dalli, revealed by mitochondrial and microsatellite DNA analyses

Mitochondrial DNA (mtDNA) control-region sequences and microsatellite loci length polymorphisms were used to estimate phylogeographical patterns (historical patterns underlying contemporary distribution), intraspecific population structure and gender-biased dispersal of Phocoenoides dalli dalli across its entire range. One-hundred and thirteen animals from several geographical strata were sequenced over 379 bp of mtDNA, resulting in 58 mtDNA haplotypes. Analysis using F(ST) values (based on haplotype frequencies) and phi(ST) values (based on frequencies and genetic distances between haplotypes) yielded statistically significant separation (bootstrap values P < 0.05) among most of the stocks currently used for management purposes. A minimum spanning network of haplotypes showed two very distinctive clusters, differentially occupied by western and eastern populations, with some common widespread haplotypes. This suggests some degree of phyletic radiation from west to east, superimposed on gene flow. Highly male-biased migration was detected for several population comparisons. Nuclear microsatellite DNA markers (119 individuals and six loci) provided additional support for population subdivision and gender-biased dispersal detected in the mtDNA sequences. Analysis using F(ST) values (based on allelic frequencies) yielded statistically significant separation between some, but not all, populations distinguished by mtDNA analysis. R(ST) values (based on frequencies of and genetic distance between alleles) showed no statistically significant subdivision. Again, highly male-biased dispersal was detected for all population comparisons, suggesting, together with morphological and reproductive data, the existence of sexual selection. Our molecular results argue for nine distinct dalli-type populations that should be treated as separate units for management purposes.     

Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes

Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N(e) estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N(e)) that are very large, making small selection intensities significant evolutionary forces in this population.     

Phylogeography of the Temminck’s Stint (Calidris temminckii): historical vicariance but little present genetic structure in a regionally endangered Palearctic wader

Aim We study the population differentiation and phylogeography of the Temminck’s Stint (Calidris temminckii). Specifically, we seek signs of past and present population size changes and dispersal events and evaluate management and conservation unit status of the populations. We also study the possibility of introgression as the origin of two mitochondrial DNA (mtDNA) lineages found and estimate the divergence time of the lineages. Location Northern Eurasia. Methods We analysed 583 bp of mtDNA control region domains I and II and 11 microsatellite loci from 13 localities throughout the breeding range. In addition, we used mitochondrial cytochrome c oxidase subunit I (COI), a barcoding gene, to search for signs of introgression. Results More population differentiation was found from microsatellites than from mtDNA, although differentiation was weak in both markers. Signs of past population growth were observed, in addition to more recent decline in some areas. Both control region and COI sequences revealed two maternal lineages coexisting in Fennoscandia and in north‐west Siberia. No signs of introgression were detected. Lineage divergence time was estimated to have occurred during the glacial periods of Pleistocene. Main conclusions Slight differences in mtDNA and microsatellite differentiation and diversity may reflect different features – such as the mutation rate and effective population size – of the markers used, or female‐biased dispersal pattern and high male site‐fidelity of the species. The coexistence of the two mitochondrial lineages is most likely a consequence of post‐glacial mixing of two refugial Pleistocene populations. Based on genetic information alone, global conservation concerns are not imminent. However, fast decline of a marginal Bothnian Bay population and the smallness and remoteness of a Central Yakutian population warrant conservation actions.     

Phylogeography and population structure of European sea bass in the north‐east Atlantic

The European sea bass Dicentrarchus labrax represents a historically and commercially valuable species in the north‐east Atlantic, although the demographic history and the patterns of geographical structure of the species in the north‐east Atlantic remain poorly understood. The present study investigates the population genetic structure of sea bass in north‐western European waters, employing different genetic markers [a portion of the mitochondrial (mt)DNA control region and 13 nuclear microsatellites] aiming to unravel demographic history and population connectivity. The results obtained show a previously unrecognized pattern of population divergence at mtDNA, with three strikingly different lineages identified. Extant sea bass populations, including the Mediterranean lineage, derive from an Atlantic ancestor. A much increased number of nuclear microsatellite loci (comparatively to previous studies) still fail to detect biologically meaningful patterns of spatial genetic structuring in the North Atlantic. Past Pleistocene glacial and interglacial events and some degree of female philopatry might be at the basis of the current geographical separation of the Atlantic lineages that has been identified. Signatures of sudden demographic expansions are more evident in the most recent mitochondrial lineages, and their slight, yet significant, geographical segregation leads to the hypothesis that present‐day spawning grounds for European sea bass may still to some extent be linked to their most recent glacial refugia. © 2011 The Linnean Society of London, Biological Journal of the Linnean Society, 2011, 104 , 364–377.     

Long-Range Gene Flow and the Effects of Climatic and Ecological Factors on Genetic Structuring in a Large,Solitary Carnivore: The Eurasian Lynx

Due to their high mobility, large terrestrial predators are potentially capable of maintaining high connectivity, and therefore low genetic differentiation among populations. However, previous molecular studies have provided contradictory findings in relation to this. To elucidate patterns of genetic structure in large carnivores, we studied the genetic variability of the Eurasian lynx, Lynx lynx throughout north-eastern Europe using microsatellite, mitochondrial DNA control region and Y chromosome-linked markers. Using SAMOVA we found analogous patterns of genetic structure based on both mtDNA and microsatellites, which coincided with a relatively little evidence for male-biased dispersal. No polymorphism for the cytochrome b and ATP6 mtDNA genes and Y chromosome-linked markers were found. Lynx inhabiting a large area encompassing Finland, the Baltic countries and western Russia formed a single genetic unit, while some marginal populations were clearly divergent from others. The existence of a migration corridor was suggested to correspond with distribution of continuous forest cover. The lowest variability (in both markers) was found in lynx from Norway and Białowieża Primeval Forest (BPF), which coincided with a recent demographic bottleneck (Norway) or high habitat fragmentation (BPF). The Carpathian population, being monomorphic for the control region, showed relatively high microsatellite diversity, suggesting the effect of a past bottleneck (e.g. during Last Glacial Maximum) on its present genetic composition. Genetic structuring for the mtDNA control region was best explained by latitude and snow cover depth. Microsatellite structuring correlated with the lynx''s main prey, especially the proportion of red deer (Cervus elaphus) in its diet. Eurasian lynx are capable of maintaining panmictic populations across eastern Europe unless they are severely limited by habitat continuity or a reduction in numbers. Different correlations of mtDNA and microsatellite population divergence patterns with climatic and ecological factors may suggest separate selective pressures acting on males and females in this solitary carnivore.     

High gene flow on a continental scale in the polyandrous Kentish plover Charadrius alexandrinus

Gene flow promotes genetic homogeneity of species in time and space. Gene flow can be modulated by sex‐biased dispersal that links population genetics to mating systems. We investigated the phylogeography of the widely distributed Kentish plover Charadrius alexandrinus. This small shorebird has a large breeding range spanning from Western Europe to Japan and exhibits an unusually flexible mating system with high female breeding dispersal. We analysed genetic structure and gene flow using a 427‐bp fragment of the mitochondrial (mtDNA) control region, 21 autosomal microsatellite markers and a Z microsatellite marker in 397 unrelated individuals from 21 locations. We found no structure or isolation‐by‐distance over the continental range. However, island populations had low genetic diversity and were moderately differentiated from mainland locations. Genetic differentiation based on autosomal markers was positively correlated with distance between mainland and each island. Comparisons of uniparentally and biparentally inherited markers were consistent with female‐biased gene flow. Maternally inherited mtDNA was less structured, whereas the Z‐chromosomal marker was more structured than autosomal microsatellites. Adult males were more related than females within genetic clusters. Taken together, our results suggest a prominent role for polyandrous females in maintaining genetic homogeneity across large geographic distances.     

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.     

Inferring sex-biased dispersal from population genetic tools: a review

Sex-biased dispersal, where individuals of one sex stay or return to their natal site (or group) to breed while individuals of the other sex are prone to disperse, is a wide-spread pattern in vertebrate organisms. In general, mammals exhibit male-biased dispersal whereas birds exhibit female-bias. Dispersal estimates are often difficult to obtain from direct field observations. Here we describe different methods for inferring sex-specific dispersal using population genetic tools and discuss the problems they can raise. We distinguish two types of methods: those based on bi-parental markers (eg comparison of male/female relatedness, F(st) and assignment probabilities) and those relying on the comparison between markers with different modes of inheritance (eg mtDNA markers and microsatellites). Finally, we discuss statistical problems that are encountered with these different methods (eg pseudoreplication, problems due to the comparison of distinct markers). While the genetic methods to detect sex-biased dispersal are now relatively well developed, their interpretation can prove problematic due to the confounding effects of factors such as the mating system of the species. Moreover, the relative power of these methods is not well known and requires further investigation.     

Genetic structure and mating system of Euterpe edulis Mart. Populations: a comparative analysis using microsatellite and allozyme markers

A comparative study between microsatellite and allozyme markers was conducted on the genetic structure and mating system in natural populations of Euterpe edulis Mart. Three cohorts, including seedlings, saplings, and adults, were examined in 4 populations using 10 allozyme loci and 10 microsatellite loci. As expected, microsatellite markers had a much higher degree of polymorphism than allozymes, but estimates of multilocus outcrossing rate ( = 1.00), as well as estimates of genetic structure (F(IS), G(ST)), were similar for the 2 sets of markers. Estimates of R(ST), for microsatellites, were higher than those of G(ST), but results of both statistics revealed a close agreement for the genetic structure of the species. This study provides support for the important conclusion that allozymes are still useful and reliable markers to estimate population genetic parameters. Effects of sample size on estimates from hypervariable loci are also discussed in this paper.     

Sea snakes rarely venture far from home

The extent to which populations are connected by dispersal influences all aspects of their biology and informs the spatial scale of optimal conservation strategies. Obtaining direct estimates of dispersal is challenging, particularly in marine systems, with studies typically relying on indirect approaches to evaluate connectivity. To overcome this challenge, we combine information from an eight-year mark-recapture study with high-resolution genetic data to demonstrate extremely low dispersal and restricted gene flow at small spatial scales for a large, potentially mobile marine vertebrate, the turtleheaded sea snake (Emydocephalus annulatus). Our mark-recapture study indicated that adjacent bays in New Caledonia (<1.15 km apart) contain virtually separate sea snake populations. Sea snakes could easily swim between bays but rarely do so. Of 817 recaptures of marked snakes, only two snakes had moved between bays. We genotyped 136 snakes for 11 polymorphic microsatellite loci and found statistically significant genetic divergence between the two bays (F(ST)= 0.008, P < 0.01). Bayesian clustering analyses detected low mixed ancestry within bays and genetic relatedness coefficients were higher, on average, within than between bays. Our results indicate that turtleheaded sea snakes rarely venture far from home, which has strong implications for their ecology, evolution, and conservation.     

Scale-specific sex-biased dispersal in the Valais shrew unveiled by genetic variation on the Y chromosome, autosomes, and mitochondrial DNA

We investigated sex specificities in the evolutionary processes shaping Y chromosome, autosomes, and mitochondrial DNA patterns of genetic structure in the Valais shrew (Sorex antinorii), a mountain dwelling species with a hierarchical distribution. Both hierarchical analyses of variance and isolation-by-distance analyses revealed patterns of population structure that were not consistent across maternal, paternal, and biparentally inherited markers. Differentiation on a Y microsatellite was lower than expected from the comparison with autosomal microsatellites and mtDNA, and it was mostly due to genetic variance among populations within valleys, whereas the opposite was observed on other markers. In addition, there was no pattern of isolation by distance for the Y, whereas there was strong isolation by distance on mtDNA and autosomes. We use a hierarchical island model of coancestry dynamics to discuss the relative roles of the microevolutionary forces that may induce such patterns. We conclude that sex-biased dispersal is the most important driver of the observed genetic structure, but with an intriguing twist: it seems that dispersal is strongly male biased at large spatial scale, whereas it is mildly biased in favor of females at local scale. These results add to recent reports of scale-specific sex-biased dispersal patterns, and emphasize the usefulness of the Y chromosome in conjunction with mtDNA and autosomes to infer sex specificities.     

What Maintains the Central North Pacific Genetic Discontinuity in Pacific Herring?

Pacific herring show an abrupt genetic discontinuity in the central North Pacific that represents secondary contact between refuge populations previously isolated during Pleistocene glaciations. Paradoxically, high levels of gene flow produce genetic homogeneity among ocean-type populations within each group. Here, we surveyed variability in mtDNA control-region sequences (463 bp) and nine microsatellite loci in Pacific herring from sites across the North Pacific to further explore the nature of the genetic discontinuity around the Alaska Peninsula. Consistent with previous studies, little divergence (Φ_ST  = 0.011) was detected between ocean-type populations of Pacific herring in the North West Pacific, except for a population in the Yellow Sea (Φ_ST  = 0.065). A moderate reduction in genetic diversity for both mtDNA and microsatellites in the Yellow Sea likely reflects founder effects during the last colonization of this sea. Reciprocal monophyly between divergent mtDNA lineages (Φ_ST  = 0.391) across the Alaska Peninsula defines the discontinuity across the North Pacific. However, microsatellites did not show a strong break, as eastern Bering Sea (EBS) herring were more closely related to NE Pacific than to NW Pacific herring. This discordance between mtDNA and microsatellites may be due to microsatellite allelic convergence or to sex-biased dispersal across the secondary contact zone. The sharp discontinuity between Pacific herring populations may be maintained by high-density blocking, competitive exclusion or hybrid inferiority.     

Colonization, dispersal, and hybridization influence phylogeography of North Atlantic sea urchins (Strongylocentrotus droebachiensis)

We used frequency-based and coalescent-based phylogeographic analysis of sea urchin (Strongylocentrotus droebachiensis) mitochondrial DNA (mtDNA) sequences and previously published microsatellite data to understand the relative influence of colonization and gene flow from older (north Pacific) and younger (northeast Atlantic) sea urchin populations on genetic variation in the northwest Atlantic. We found strong evidence of survival of northwestern Atlantic populations in local Pleistocene glacial refugia: most haplotypes were the same as or closely related to Pacific haplotypes, with deep gene genealogies that reflect divergence times within the northwestern Atlantic that are much older than the last glacial maximum. We detected gene flow across the North Atlantic in the form of haplotypes shared with or recently descended from European populations. We also found evidence of significant introgression of haplotypes from a closely related species (S. pallidus). The relative magnitude of gene flow estimated by coalescent methods (and the effective population size differences among oceanic regions) depended on the genetic marker used. In general, we found very small effective population size in the northeastern Atlantic and high trans-Arctic gene flow between the Pacific and northwestern Atlantic. Both analyses suggested significant back-migration to the Pacific. However, microsatellites more strongly reflected older Pacific migration (with similar effective population sizes across the Arctic), whereas mtDNA sequences appeared to be more sensitive to recent trans- Atlantic dispersal (with larger differences in effective population size). These differences across marker types might have several biological or methodological causes, and they suggest caution in interpretation of the results from a single locus or class of markers.     

When east meets west: population structure of a high-latitude resident species,the boreal chickadee (Poecile hudsonicus)

The population genetic structure of northern boreal species has been strongly influenced both by the Quaternary glaciations and the presence of contemporary barriers, such as mountain ranges and rivers. We used a combination of mitochondrial DNA (mtDNA), nuclear microsatellites and spatial distribution modelling to study the population genetic structure of the boreal chickadee (Poecile hudsonicus), a resident passerine, and to investigate whether historical or contemporary barriers have influenced this northern species. MtDNA data showed evidence of eastern and western groups, with secondary admixture occurring in central Canada. This suggests that the boreal chickadee probably persisted in multiple glacial refugia, one in Beringia and at least one in the east. Palaeo-distribution modelling identified suitable habitat in Beringia (Alaska), Atlantic Canada and the southern United States, and correspond to divergence dates of 60–96 kya. Pairwise F_ST values for both mtDNA and microsatellites were significant for all comparisons involving Newfoundland, though mtDNA data suggest a more recent separation. Furthermore, unlike mtDNA data, nuclear data support population connectivity among the continental populations, possibly due to male-biased dispersal. Although both are significant, the isolation-by-distance signal is much stronger for mtDNA (r²=0.51) than for microsatellites (r²=0.05), supporting the hypothesis of male-biased dispersal. The population structure of the boreal chickadee was influenced by isolation in multiple refugia and contemporary barriers. In addition to geographical distance, physical barriers such as the Strait of Belle Isle and northern mountains in Alaska are restricting gene flow, whereas the Rocky Mountains in the west are a porous barrier.     

Population genetic structure and male-biased dispersal in the queenless ant Diacamma cyaneiventre

In this study we investigated the population genetic structure of the queenless ant Diacamma cyaneiventre. This species, lacking winged queens, is likely to have a restricted female dispersal. We used both mitochondrial and microsatellite markers to assess the consequence of such restricted female dispersal at three geographical scales: within a given locality (< 1 km), between localities within a given region (< 10 km) and between regions (> 36 km). Within a locality, a strong population structure was observed for mitochondrial DNA (mtDNA) whereas weak or nonexistent population genetic structure was observed for the microsatellites (around 5% of the value for mtDNA). Male gene flow was estimated to be about 20-30 times higher than female gene flow at this scale. At a larger spatial scale, very strong genetic differentiation for both markers was observed between localities - even within a single region. Female dispersal is nonexistent at these scales and male dispersal is very restricted, especially between regions. The phylogeographical structure of the mtDNA haplotypes as well as the very low genetic diversity of mtDNA within localities indicate that new sites are colonized by a single migration event from adjacent localities, followed by successive colony fissions. These patterns of genetic variability and differentiation agree with what is theoretically expected when colonization events are kin-structured and when, following colonization, dispersion is mainly performed by males.     

Sex-biased gene flow and colonization in the Formosan lesser horseshoe bat: inference from nuclear and mitochondrial markers

Sex-biased behaviours are expected to play an important role in partitioning genetic variance in animal populations. Comparing genetic structure at markers with different modes of inheritance provides a means of detecting these behaviours and their consequences for population genetic structure. In colonially breeding mammals, the common combination of female philopatry and male vagility can promote contrasting patterns of genetic differentiation between the sexes, both via their effects on recurrent gene flow and on colonization. We examined sex differences in gene flow and structure by comparing maternally inherited mitochondrial DNA (mtDNA) and biparentally inherited autosomal loci in the Formosan lesser horseshoe bat Rhinolophus monoceros . We found that genetic partitioning was higher at mtDNA than autosomal markers in both sexes, indicative of female-biased philopatry and male-biased dispersal. Across Taiwan, isolation-by-distance was detected for all sex/marker combinations but was steeper for mtDNA than for nuclear markers. We suggest that isolation-by-distance shown from mtDNA at large scales is likely to reflect the stepwise founding of new breeding colonies by females during colonization. In contrast, no isolation-by-distance was found at smaller distances of up to 100 km, indicating that gene flow and/or recent shared ancestry homogenises genetic structure among nearby sites. Our results highlight the value of an indirect genetic approach to understanding sex-biased behaviours and their consequences in a little-studied species.     

Multilocus phylogeography of the sea snake Hydrophis curtus reveals historical vicariance and cryptic lineage diversity

The Indo‐Australian archipelago (IAA) supports the world's highest diversity of marine fish, invertebrates and reptiles. Many of the marine fish and invertebrates show congruent phylogeographic patterns, supporting a view that the region's complex geo‐climatic history has played an important role in generating its exceptional biodiversity. Here, we examine population genetic structure of the viviparous sea snake, Hydrophis curtus, to assess how past and present barriers to gene flow in the IAA have contributed to genetic and species diversity in a fully marine reptile. Mitochondrial and anonymous nuclear sequences and ten microsatellite loci were used to identify patterns of historical genetic structure and population expansion, reconstruct dated genealogies and assess levels of recent gene flow. These markers revealed strong concordant geographic structure within H. curtus with a prominent genetic break between populations broadly distributed in the Indian Ocean and the West Pacific. These populations were estimated to have diverged in the late Pliocene or early Pleistocene, and microsatellite admixture analyses suggested limited recent gene flow between them despite the current lack of barriers to dispersal, indicating possible cryptic species. Subsequent divergence in the mid–late Pleistocene was detected within the West Pacific clade among the populations in the Phuket‐Thailand region, South‐East Asia and Australia, and two of these populations also showed genetic signals of recent range expansions. Our results show that climatic fluctuations during the Plio‐Pleistocene generated high levels of cryptic genetic diversity in H. curtus, and add to similar findings for diverse other marine groups in the IAA.