Microsatellite and mitochondrial DNA analyses of the genetic structure of blacktip shark (Carcharhinus limbatus) nurseries in the northwestern Atlantic, Gulf of Mexico, and Caribbean Sea

Phylogeography and population structure of the red stingray, Dasyatis akajei inferred by mitochondrial control region

The red stingray Dasyatis akajei is distributed in both marine and freshwater, but little is known about its phylogeography and population structure. We sampled 107 individuals from one freshwater region and 6 coastal localities within the distribution range of D. akajei. Analyses of the first hypervariable region of mitochondrial DNA control region of 474 bp revealed only 17 polymorphism sites that defined 28 haplotypes, with no unique haplotype for the freshwater population. A high level of haplotype diversity and low nucleotide diversity were observed in both marine (h = 0.9393 ± 0.0104, π = 0.0069 ± 0.0040) and freshwater populations (h = 0.8333 ± 0.2224, π = 0.0084 ± 0.0063). Significant level of genetic structure was detected between four marine populations (TZ, WZ, ND and ZZ) via both hierarchical molecular variance analysis (AMOVA) and pairwise FST (with two exceptions), which is unusual for elasmobranchs detected previously over such short geographical distance. However, limited sampling suggested that the freshwater population was not particularly distinct (p > 0.05), but additional samples would be needed to confirm it. Demersal and slow-moving characters likely have contributed to the genetically heterogeneous population structure. The demographic history of D. akajei examined by mismatch distribution analyses, neutrality tests and Bayesian skyline analyses suggested a sudden population expansion dating to upper Pleistocene. The information on genetic diversity and genetic structure will have implications for the management of fisheries and conservation efforts.

Population expansion and genetic structure in Carcharhinus brevipinna in the southern Indo-Pacific

BACKGROUND

Quantifying genetic diversity and metapopulation structure provides insights into the evolutionary history of a species and helps develop appropriate management strategies. We provide the first assessment of genetic structure in spinner sharks (Carcharhinus brevipinna), a large cosmopolitan carcharhinid, sampled from eastern and northern Australia and South Africa.

METHODS AND FINDINGS

Sequencing of the mitochondrial DNA NADH dehydrogenase subunit 4 gene for 430 individuals revealed 37 haplotypes and moderately high haplotype diversity (h = 0.6770 ±0.025). While two metrics of genetic divergence (ΦST and F ST) revealed somewhat different results, subdivision was detected between South Africa and all Australian locations (pairwise ΦST, range 0.02717-0.03508, p values ≤ 0.0013; pairwise F ST South Africa vs New South Wales = 0.04056, p = 0.0008). Evidence for fine-scale genetic structuring was also detected along Australia's east coast (pairwise ΦST = 0.01328, p < 0.015), and between south-eastern and northern locations (pairwise ΦST = 0.00669, p < 0.04).

CONCLUSIONS

The Indian Ocean represents a robust barrier to contemporary gene flow in C. brevipinna between Australia and South Africa. Gene flow also appears restricted along a continuous continental margin in this species, with data tentatively suggesting the delineation of two management units within Australian waters. Further sampling, however, is required for a more robust evaluation of the latter finding. Evidence indicates that all sampled populations were shaped by a substantial demographic expansion event, with the resultant high genetic diversity being cause for optimism when considering conservation of this commercially-targeted species in the southern Indo-Pacific.

Phylogeography of sea oats (Uniola paniculata), a dune-building coastal grass in southeastern North America

We tested the hypothesis that Uniola paniculata populations are divided into eastern and western lineages, with the primary geographic break at the southern tip of Florida, as observed in codistributed animal taxa. We asked whether the geographic distribution of chloroplast DNA (cpDNA) variation in U. paniculata corresponds to 1) genetic structure in nuclear variation reported in previous studies, and 2) the geographic distribution of morphological adaptive traits reported in previous studies. We sampled 66 populations and performed phylogeographic analyses using sequence variations in maternally inherited cpDNA. We reconstructed the intraspecific phylogenetic network with TCS software and identified phylogeographic breaks in the species using Monmonier's algorithm. Analyses identified 6 cpDNA haplotypes and 2 major lineages: eastern (Atlantic) and western (Gulf), with a phylogeographic break at the southern tip of Florida. The data suggest U. paniculata survived the last glacial maximum (LGM) in southern refugia. Following the LGM, differential leading-edge recolonization explains the current distribution of haplotypes into 2 lineages. Populations containing a haplotype from outside its native range are likely due to human-mediated transplantation. The genetic structure of cpDNA variation has weak correlation with nuclear DNA variation, and there is partial concordance between the geographic distribution of cpDNA and morphological variation.

Microsatellite analyses of blacktip reef sharks (Carcharhinus melanopterus) in a fragmented environment show structured clusters

The population dynamics of shark species are generally poorly described because highly mobile marine life is challenging to investigate. Here we investigate the genetic population structure of the blacktip reef shark (Carcharhinus melanopterus) in French Polynesia. Five demes were sampled from five islands with different inter-island distances (50–1500 km). Whether dispersal occurs between islands frequently enough to prevent moderate genetic structure is unknown. We used 11 microsatellites loci from 165 individuals and a strong genetic structure was found among demes with both F-statistics and Bayesian approaches. This differentiation is correlated with the geographic distance between islands. It is likely that the genetic structure seen is the result of all or some combination of the following: low gene flow, time since divergence, small effective population sizes, and the standard issues with the extent to which mutation models actually fit reality. We suggest low levels of gene flow as at least a partial explanation of the level of genetic structure seen among the sampled blacktip demes. This explanation is consistent with the ecological traits of blacktip reef sharks, and that the suitable habitat for blacktips in French Polynesia is highly fragmented. Evidence for spatial genetic structure of the blacktip demes we studied highlights that similar species may have populations with as yet undetected or underestimated structure. Shark biology and the market for their fins make them highly vulnerable and many species are in rapid decline. Our results add weight to the case that total bans on shark fishing are a better conservation approach for sharks than marine protected area networks.

Ancient divergence in the trans-oceanic deep-sea shark Centroscymnus crepidater

Unravelling the genetic structure and phylogeographic patterns of deep-sea sharks is particularly challenging given the inherent difficulty in obtaining samples. The deep-sea shark Centroscymnus crepidater is a medium-sized benthopelagic species that exhibits a circumglobal distribution occurring both in the Atlantic and Indo-Pacific Oceans. Contrary to the wealth of phylogeographic studies focused on coastal sharks, the genetic structure of bathyal species remains largely unexplored. We used a fragment of the mitochondrial DNA control region, and microsatellite data, to examine genetic structure in C. crepidater collected from the Atlantic Ocean, Tasman Sea, and southern Pacific Ocean (Chatham Rise). Two deeply divergent (3.1%) mtDNA clades were recovered, with one clade including both Atlantic and Pacific specimens, and the other composed of Atlantic samples with a single specimen from the Pacific (Chatham Rise). Bayesian analyses estimated this splitting in the Miocene at about 15 million years ago. The ancestral C. crepidater lineage was probably widely distributed in the Atlantic and Indo-Pacific Oceans. The oceanic cooling observed during the Miocene due to an Antarctic glaciation and the Tethys closure caused changes in environmental conditions that presumably restricted gene flow between basins. Fluctuations in food resources in the Southern Ocean might have promoted the dispersal of C. crepidater throughout the northern Atlantic where habitat conditions were more suitable during the Miocene. The significant genetic structure revealed by microsatellite data suggests the existence of present-day barriers to gene flow between the Atlantic and Pacific populations most likely due to the influence of the Agulhas Current retroflection on prey movements.

Direct genetic evidence for reproductive philopatry and associated fine-scale migrations in female blacktip reef sharks (Carcharhinus melanopterus) in French Polynesia

Conservation of top predators has been emphasized as essential in an ecosystem due to their role in trophic chain regulation. Optimizing conservation strategies for these endangered marine top predators requires direct estimates of breeding patterns and connectivity as these are essential to understanding the population dynamics. There have been some attempts to investigate breeding patterns of reef sharks from litter reconstruction using molecular analyses. However, direct fine-scale migrations of female sharks for parturition as well as connectivity at a medium scale like between islands remain mostly unknown. We used microsatellite DNA markers and a likelihood-based parentage analysis to determine breeding patterns of female blacktip reef sharks in Moorea (Society Islands, French Polynesia). Most females gave birth at their home island but some migrated to specific nursery areas outside the area they are attached to, sometimes going to another island 50 km away across deep ocean. Our analysis also revealed that females migrated to the same nursery for every birthing event. Many offspring showed a high level of inbreeding indicating an overall reduced population size, restricted movements and dispersal, or specific mating behaviour. Females represent the vectors that transport the genes at nursery grounds, and their fidelity should thus define reproductive units. As females seem to be philopatric, males could be the ones dispersing genes between populations. These results highlight the need to conserve coastal zones where female reef sharks seem to exhibit philopatry during the breeding season.

Inclusion of South American samples reveals new population structuring of the blacktip shark (Carcharhinus limbatus) in the western Atlantic

Carcharhinus limbatus has a cosmopolitan distribution and marked genetic structuring, mainly because of its philopatric behavior. However, analysis of this structuring has not previously included South American populations. In the present study, we analyzed a sample of adult individuals collected on the northern coast of Brazil and compared the sequences of the mitochondrial control region with those of populations already genotyped. Relatively high haplotype diversity (12 haplotypes, genetic diversity of 0.796) was observed, similar to that in other populations but with a much larger number of private alleles. In contrast to populations studied previously, which were represented by neonates, the pronounced allelic variability found in the South American individuals may have resulted from migrations from other populations in the region that have yet to be genotyped. This population was also genetically distinct from the other Atlantic populations (F_st > 0.8), probably because of female philopatry, and apparently separated from the northwestern Atlantic group 1.39 million years ago. These findings indicate that the C. limbatus population from northern Brazil is genetically distinct from all other populations and should be considered as a different management unit for the protection of stocks.

Common misconceptions in molecular ecology: echoes of the modern synthesis

The field of molecular ecology has burgeoned into a large discipline spurred on by technical innovations that facilitate the rapid acquisition of large amounts of genotypic data, by the continuing development of theory to interpret results, and by the availability of computer programs to analyse data sets. As the discipline grows, however, misconceptions have become enshrined in the literature and are perpetuated by routine citations to other articles in molecular ecology. These misconceptions hamper a better understanding of the processes that influence genetic variation in natural populations and sometimes lead to erroneous conclusions. Here, we consider eight misconceptions commonly appearing in the literature: (i) some molecular markers are inherently better than other markers; (ii) mtDNA produces higher F(ST) values than nDNA; (iii) estimated population coalescences are real; (iv) more data are always better; (v) one needs to do a Bayesian analysis; (vi) selective sweeps influence mtDNA data; (vii) equilibrium conditions are critical for estimating population parameters; and (viii) having better technology makes us smarter than our predecessors. This is clearly not an exhaustive list and many others can be added. It is, however, sufficient to illustrate why we all need to be more critical of our own understanding of molecular ecology and to be suspicious of self-evident truths.

Evidence for reproductive philopatry in the bull shark Carcharhinus leucas

Reproductive philopatry in bull sharks Carcharhinus leucas was investigated by comparing mitochondrial (NADH dehydrogenase subunit 4, 797 base pairs and control region genes 837 base pairs) and nuclear (three microsatellite loci) DNA of juveniles sampled from 13 river systems across northern Australia. High mitochondrial and low microsatellite genetic diversity among juveniles sampled from different rivers (mitochondrial φ(ST) = 0·0767, P < 0·05; microsatellite F(ST) = -0·0022, P > 0·05) supported female reproductive philopatry. Genetic structure was not further influenced by geographic distance (P > 0·05) or long-shore barriers to movement (P > 0·05). Additionally, results suggest that C. leucas in northern Australia has a long-term effective population size of 11 000-13 000 females and has undergone population bottlenecks and expansions that coincide with the timing of the last ice-ages.

Molecular markers: progress and prospects for understanding reproductive ecology in elasmobranchs

Application of modern molecular tools is expanding the understanding of elasmobranch reproductive ecology. High-resolution molecular markers provide information at scales ranging from the identification of reproductively isolated populations in sympatry (i.e. cryptic species) to the relationships among parents, offspring and siblings. This avenue of study has not only augmented the current understanding of the reproductive biology of elasmobranchs but has also provided novel insights that could not be obtained through experimental or observational techniques. Sharing of genetic polymorphisms across ocean basins indicates that for some species there may be gene flow on global scales. The presence, however, of morphologically similar but genetically distinct entities in sympatry suggests that reproductive isolation can occur with minimal morphological differentiation. This review discusses the recent findings in elasmobranch reproductive biology like philopatry, hybridization and polyandry while highlighting important molecular and analytical techniques. Furthermore, the review examines gaps in current knowledge and discusses how new technologies may be applied to further the understanding of elasmobranch reproductive ecology.

A review of the application of molecular genetics for fisheries management and conservation of sharks and rays

Since the first investigation 25 years ago, the application of genetic tools to address ecological and evolutionary questions in elasmobranch studies has greatly expanded. Major developments in genetic theory as well as in the availability, cost effectiveness and resolution of genetic markers were instrumental for particularly rapid progress over the last 10 years. Genetic studies of elasmobranchs are of direct importance and have application to fisheries management and conservation issues such as the definition of management units and identification of species from fins. In the future, increased application of the most recent and emerging technologies will enable accelerated genetic data production and the development of new markers at reduced costs, paving the way for a paradigm shift from gene to genome-scale research, and more focus on adaptive rather than just neutral variation. Current literature is reviewed in six fields of elasmobranch molecular genetics relevant to fisheries and conservation management (species identification, phylogeography, philopatry, genetic effective population size, molecular evolutionary rate and emerging methods). Where possible, examples from the Indo-Pacific region, which has been underrepresented in previous reviews, are emphasized within a global perspective.

Habitat use and spatial segregation of adult spottail sharks Carcharhinus sorrah in tropical nearshore waters

An array of acoustic receivers deployed in Cleveland Bay, north Queensland, Australia, passively tracked 20 adult spottail sharks Carcharhinus sorrah over 2 years (2009-2010) to define patterns in movement and habitat use. Individuals were present in the study site for long periods, ranging from 8 to 408 days (mean = 185). Size and location of home ranges did not vary over time. A high level of segregation occurred among C. sorrah, with individuals using different types of habitat and showing strong attachment to specific regions. The depth of habitat individuals used varied between sexes. Males tended to use a narrow range of habitat depths within the study site (2·8-6·0 m), whereas females used shallower habitats (1·4-6·2 m) and displayed a seasonal shift in the depth of habitat used. Mean monthly habitat depth used varied by as much as 2 m for females, with individuals using shallower habitats during the winter months. Long-term presence and consistent home ranges suggest that Cleveland Bay provides important habitat for C. sorrah. By defining patterns in the use of nearshore habitats for C. sorrah, this study improves the understanding of the movement and habitat use of smaller-bodied coastal sharks and may help provide guidance for the management of their populations.

Population structure and phylogeography of the short-tailed stingray, Dasyatis brevicaudata (Hutton 1875), in the Southern Hemisphere

There is accumulating evidence that the degree of vagility explains little of the extent of population subdivision found within elasmobranch species. Instead, patterns of gene flow in elasmobranchs appear more closely linked to the presence of dispersal barriers, either physical or biological. Here, we investigate the potential role of some of these isolating mechanisms in shaping the population structure of a member of the stingray family Dasyatidae (Dasyatis brevicaudata) at various scales (southern hemisphere vs. coastal New Zealand). Analyses of the mitochondrial DNA control region from 176 individuals revealed significant genetic structure between South Africa, Australia, and New Zealand populations (analysis of molecular variance [AMOVA], overall Ф(ST) = 0.67, P < 0.001), although New Zealand and Australia shared some haplotypes. Surprisingly, significant population differentiation was found among several coastal New Zealand locations (AMOVA, overall Ф(ST) = 0.05, P < 0.05). However, data did not support the genetic differentiation between individuals from an offshore breeding area and mainland individuals. Comparisons suggest that these stingrays exhibit similar levels of population differentiation as other coastal elasmobranchs, with high divergence across oceanic basins and lower differentiation along continuous coastal habitats. Differences in coastal population structuring in elasmobranch species studied to date may be attributed to species-specific preferences for coastal habitats, which may be linked to life history functions (e.g., feeding and pupping).

Genetic diversity despite population collapse in a critically endangered marine fish: the smalltooth sawfish (Pristis pectinata)

Sawfish (family Pristidae) are among the most critically endangered marine fish in the world, yet very little is known about how genetic bottlenecks, genetic drift, and inbreeding depression may be affecting these elasmobranchs. In the US Atlantic, the smalltooth sawfish (Pristis pectinata) has declined to 1-5% of its abundance in the 1900s, and its core distribution has contracted to southwest Florida. We used 8 polymorphic microsatellite markers to show that this remnant population still exhibits high genetic diversity in terms of average allelic richness (18.23), average alleles per locus (18.75, standard deviation [SD] 6.6) and observed heterozygosity (0.43-0.98). Inbreeding is rare (mean individual internal relatedness = -0.02, SD 0.14; F(IS) = -0.011, 95% confidence interval [CI] = -0.039 to 0.011), even though the estimated effective population size (N(e)) is modest (250-350, 95% CI = 142-955). Simulations suggest that the remnant smalltooth sawfish population will probably retain >90% of its current genetic diversity over the next century even at the lower estimate of N(e). There is no evidence of a genetic bottleneck accompanying last century's demographic bottleneck, and we discuss hypotheses that could explain this. We also discuss features of elasmobranch life history and population biology that could make them less vulnerable than other large marine vertebrates to genetic change associated with reduced population size.

High connectivity on a global scale in the pelagic wahoo, Acanthocybium solandri (tuna family Scombridae)

The population genetic structure and phylogeography of wahoo, Acanthocybium solandri, were investigated on a global scale with intron six of lactate dehydrogenase-A (ldhA6, 8 locations, N = 213) and mtDNA cytochrome b (Cytb, 10 locations, N = 322). Results show extensive sharing of haplotypes across the wahoo's entire global range, and analyses were unable to detect significant structure (nuclear F(ST) = 0.0125, P = 0.106; mtDNA Phi(ST) < 0.0001, P = 0.634). Power analyses indicated 95% confidence in detecting nuclear F(ST) > or = 0.0389 and mtDNA Phi(ST) > or = 0.0148. These findings appear unique, as most other tunas, billfishes, and oceanic sharks exhibit significant population structure on the scale of East-West Atlantic, Atlantic vs. Indian-Pacific, or East-West Pacific. Overall nuclear heterozygosity (H = 0.714) and mtDNA haplotype diversity (h = 0.918) are both high in wahoo, while overall mtDNA nucleotide diversity (pi = 0.006) and nuclear nucleotide diversity (pi = 0.004) are uniformly low, indicating a recent increase in population size. Coalescence analyses yield an estimate of effective female population size (NeF) at approximately 816,000, and a population bottleneck approximately 690,000 years ago. However, conclusions about population history from our Cytb data set are not concordant with a control region survey, a finding that will require further investigation. This is the first example of a vertebrate with a single globally distributed population, a finding we attribute to extensive dispersal at all life stages. The indications of a worldwide stock for wahoo reinforce the mandate for international cooperation on fisheries issues.

Populations analysis of the Brazilian Sharpnose Shark Rhizoprionodon lalandii (Chondrichthyes: Carcharhinidae) on the São Paulo coast, Southern Brazil: inferences from mt DNA sequences

Sharks of the genus Rhizoprionodon can be considered some of the most important predators along the trophic coastal marine ecosystems and represent an important economic resource for the small-scale fisheries, especially on the Brazilian coastline. In order to analyze the population structure of the shark Rhizoprionodon lalandii of São Paulo, Southeastern coast of Brazil, levels of genetic diversity were identified by nucleotide sequence analyses of the mitochondrial DNA control region. The results obtained from this study present moderate values of haplotype diversity and low nucleotide diversity. Although the AMOVA tests (ΦST = 0.08394, P < 0.01) had shown slightly differences among the studied samples, evidence for the occurrence of population structuring was not found, which may be a general feature of sharks living in coastal areas.

Low genetic differentiation across three major ocean populations of the whale shark, Rhincodon typus

BACKGROUND

Whale sharks are a declining species for which little biological data is available. While these animals are protected in many parts of their range, they are fished legally and illegally in some countries. Baseline biological and ecological data are needed to allow the formulation of an effective conservation plan for whale sharks. It is not known, for example, whether the whale shark is represented by a single worldwide panmictic population or by numerous, reproductively isolated populations. Genetic analysis of population structure is one essential component of the baseline data required for whale shark conservation.

METHODOLOGY/PRINCIPAL FINDINGS

We have identified 8 polymorphic microsatellites in the whale shark and used these markers to assess genetic variation and population structure in a panel of whale sharks covering a broad geographic region. This is the first record of microsatellite loci in the whale shark, which displayed an average of 9 alleles per locus and mean H(o) = 0.66 and H(e) = 0.69. All but one of the eight loci meet the expectations of Hardy-Weinberg equilibrium. Analysis of these loci in whale sharks representing three major portions of their range, the Pacific (P), Caribbean (C), and Indian (I) Oceans, determined that there is little population differentiation between animals sampled in different geographic regions, indicating historical gene flow between populations. F(ST) values for inter-ocean comparisons were low (PxC = 0.0387, CxI = 0.0296 and PxI = -0.0022), and only CxI approached statistical significance (p = 0.0495).

CONCLUSIONS/SIGNIFICANCE

We have shown only low levels of genetic differentiation between geographically distinct whale shark populations. Existing satellite tracking data have revealed both regional and long-range migration of whale sharks throughout their range, which supports the finding of gene flow between populations. Whale sharks traverse geographic and political boundaries during their life history and interbreed with animals from distant populations; conservation efforts must therefore target international protection for this species.

IUCN classification zones concord with, but underestimate, the population genetic structure of the zebra shark Stegostoma fasciatum in the Indo-West Pacific

The Indo-West Pacific (IWP), from South Africa in the western Indian Ocean to the western Pacific Ocean, contains some of the most biologically diverse marine habitats on earth, including the greatest biodiversity of chondrichthyan fishes. The region encompasses various densities of human habitation leading to contrasts in the levels of exploitation experienced by chondrichthyans, which are targeted for local consumption and export. The demersal chondrichthyan, the zebra shark, Stegostoma fasciatum, is endemic to the IWP and has two current regional International Union for the Conservation of Nature (IUCN) Red List classifications that reflect differing levels of exploitation: 'Least Concern' and 'Vulnerable'. In this study, we employed mitochondrial ND4 sequence data and 13 microsatellite loci to investigate the population genetic structure of 180 zebra sharks from 13 locations throughout the IWP to test the concordance of IUCN zones with demographic units that have conservation value. Mitochondrial and microsatellite data sets from samples collected throughout northern Australia and Southeast Asia concord with the regional IUCN classifications. However, we found evidence of genetic subdivision within these regions, including subdivision between locations connected by habitat suitable for migration. Furthermore, parametric F(ST) analyses and Bayesian clustering analyses indicated that the primary genetic break within the IWP is not represented by the IUCN classifications but rather is congruent with the Indonesian throughflow current. Our findings indicate that recruitment to areas of high exploitation from nearby healthy populations in zebra sharks is likely to be minimal, and that severe localized depletions are predicted to occur in zebra shark populations throughout the IWP region.

A genetic assessment of polyandry and breeding-site fidelity in lemon sharks

We here employ 11 microsatellite markers and recently developed litter reconstruction methods to infer mating system parameters (i.e. polyandry and breeding-site fidelity) at a lemon shark nursery site in Marquesas Key, Florida. Four hundred and eight juvenile or subadult sharks were genotyped over eight complete breeding seasons. Using this information, we were able to infer family structure, as well as fully or partially reconstruct genotypes of 46 mothers and 163 fathers. Multiple litter reconstruction methods were used, and novel simulations helped define apparent bias and precision of at least some mating system parameters. For Marquesas Key, we find that adult female lemon sharks display high levels of polyandry (81% of all litters sampled) and stronger fidelity to the nursery site than do males. Indeed, few male sharks sired offspring from more than one litter during the course of the study. These findings were quite similar to previous results from another lemon shark nursery site (Bimini, Bahamas), suggesting conserved mating system parameters despite significant variation in early life-history traits (i.e. body size and growth) among sites. The finding of at least some site fidelity in females also supports the need for careful conservation of each nursery.

Population genetic structure of Earth's largest fish, the whale shark (Rhincodon typus)

Large pelagic vertebrates pose special conservation challenges because their movements generally exceed the boundaries of any single jurisdiction. To assess the population structure of whale sharks (Rhincodon typus), we sequenced complete mitochondrial DNA control regions from individuals collected across a global distribution. We observed 51 single site polymorphisms and 8 regions with indels comprising 44 haplotypes in 70 individuals, with high haplotype (h = 0.974 +/- 0.008) and nucleotide diversity (pi = 0.011 +/- 0.006). The control region has the largest length variation yet reported for an elasmobranch (1143-1332 bp). Phylogenetic analyses reveal no geographical clustering of lineages and the most common haplotype was distributed globally. The absence of population structure across the Indian and Pacific basins indicates that oceanic expanses and land barriers in Southeast Asia are not impediments to whale shark dispersal. We did, however, find significant haplotype frequency differences (AMOVA, Phi(ST) = 0.107, P < 0.001) principally between the Atlantic and Indo-Pacific populations. In contrast to other recent surveys of globally distributed sharks, we find much less population subdivision and no evidence for cryptic evolutionary partitions. Discovery of the mating and pupping areas of whale sharks is key to further population genetic studies. The global pattern of shared haplotypes in whale sharks provides a compelling argument for development of broad international approaches for management and conservation of Earth's largest fish.

Swimming in the deep end of the gene pool: global population structure of an oceanic giant

Despite the impression held by some that few biological mysteries remain, even evocative species such as humpback whales (Megaptera novaeangliae), white sharks (Carcharodon carcharias) and green turtles (Chelonia mydas) have poorly documented movement patterns, reproductive strategies and population dynamics despite years of dedicated research. This is largely due to the difficulty of observing wide-ranging marine species over the majority of their life cycle. The advent of powerful tracking devices has certainly improved our understanding, but it is usually only with molecular tools that the nature of population structure becomes apparent. In this issue of Molecular Ecology, Castro and colleagues have provided the first global-scale assessment of population structure for the largest fish--whale sharks (Rhincodon typus). Whale sharks can reach lengths > 12 m and are a popular tourist attraction at places where they aggregate, yet for most of their life cycle, we know little indeed of where they go and how they interact with other populations. Previous tracking studies imply a high dispersal capacity, but only now have Castro and colleagues demonstrated high gene flow and haplotype diversity among the major ocean basins where they are found.

Isolation and genetic diversity of endangered grey nurse shark (Carcharias taurus) populations

Anthropogenic impacts are believed to be the primary threats to the eastern Australian population of grey nurse sharks (Carcharias taurus), which is listed as critically endangered, and the most threatened population globally. Analyses of 235 polymorphic amplified fragment length polymorphisms (AFLP) loci and 700 base pairs of mitochondrial DNA control region provide the first account of genetic variation and geographical partitioning (east and west coasts of Australia, South Africa) in C. taurus. Assignment tests, analysis of relatedness and Fst values all indicate that the Australian populations are isolated from South Africa, with negligible migration between the east and west Australian coasts. There are significant differences in levels of genetic variation among regions. Australian C. taurus, particularly the eastern population, has significantly less AFLP variation than the other sampling localities. Further, the eastern Australian sharks possess only a single mitochondrial haplotype, also suggesting a small number of founding individuals. Therefore, historical, rather than anthropogenic processes most likely account for their depauperate genetic variation. These findings have implications for the viability of the eastern Australian population of grey nurse sharks.

Low worldwide genetic diversity in the basking shark (Cetorhinus maximus)

The basking shark (Cetorhinus maximus) is found in temperate waters throughout the world's oceans, and has been subjected to extensive exploitation in some regions. However, little is known about its current abundance and genetic status. Here, we investigate the diversity of the mitochondrial DNA control region among samples from the western North Atlantic, eastern North Atlantic, Mediterranean Sea, Indian Ocean and western Pacific. We find just six haplotypes defined by five variable sites, a comparatively low genetic diversity of pi=0.0013 and no significant differentiation between ocean basins. We provide evidence for a bottleneck event within the Holocene, estimate an effective population size (Ne) that is low for a globally distributed species, and discuss the implications.

Worldwide phylogeography of the blacktip shark (Carcharhinus limbatus) inferred from mitochondrial DNA reveals isolation of western Atlantic populations coupled with recent Pacific dispersal

Although many coastal shark species have widespread distributions, the genetic relatedness of worldwide populations has been examined for few species. The blacktip shark, (Carcharhinus limbatus), inhabits tropical and subtropical coastal waters throughout the world. In this study, we examined the genetic relationships of blacktip shark populations (n = 364 sharks) throughout the majority of the species' range using the entire mitochondrial control region (1067-1070 nucleotides). Two geographically distinct maternal lineages (western Atlantic, Gulf of Mexico, and Caribbean Sea clades, and eastern Atlantic, Indian, and Pacific Ocean clades) were identified and shallow population structure was detected throughout their geographic ranges. These findings indicate that a major population subdivision exists across the Atlantic Ocean, but not the Pacific Ocean. The historical dispersal of this widespread, coastal species may have been interrupted by the rise of the Isthmus of Panama. This scenario implies historical dispersal across the Pacific Ocean (supported by the recovery of the same common haplotype from the Philippines, Hawaii, and the Gulf of California reflecting recent/contemporary dispersal abilities) and an oceanic barrier to recent migration across the Atlantic. Genetic structure within the eastern Atlantic/Indo-Pacific (Phi(ST) = 0.612, P < 0.001) supports maternal philopatry throughout this area, expanding previous western Atlantic findings. Eastern Atlantic/Indo-Pacific C. limbatus control region haplotypes were paraphyletic to Carcharhinus tilstoni haplotypes in our maximum-parsimony analysis. The greater divergence of western Atlantic C. limbatus than C. tilstoni from eastern Atlantic/Indo-Pacific C. limbatus reflects the taxonomic uncertainty of western Atlantic C. limbatus.

Comparative phylogeography of unglaciated eastern North America

Regional phylogeographical studies involving co-distributed animal and plant species have been conducted for several areas, most notably for Europe and the Pacific Northwest of North America. Until recently, phylogeographical studies in unglaciated eastern North America have been largely limited to animals. As more studies emerge for diverse lineages (including plants), it seems timely to assess the phylogeography across this region: (i) comparing and contrasting the patterns seen in plants and animals; (ii) assessing the extent of pseudocongruence; and (iii) discussing the potential applications of regional phylogeography to issues in ecology, such as response to climatic change. Unglaciated eastern North America is a large, geologically and topographically complex area with the species examined having diverse distributions. Nonetheless, some recurrent patterns emerge: (i) maritime - Atlantic vs. Gulf Coast; (ii) Apalachicola River discontinuity; (iii) Tombigbee River discontinuity; (iv) the Appalachian Mountain discontinuity; (v) the Mississippi River discontinuity; and (vi) the Apalachicola River and Mississippi River discontinuities. Although initially documented in animals, most of these patterns are also apparent in plants, providing support for phylogeographical generalizations. These patterns may generally be attributable to isolation and differentiation during Pleistocene glaciation, but in some cases may be older (Pliocene). Molecular studies sometimes agree with longstanding hypotheses of glacial refugia, but also suggest additional possible refugia, such as the southern Appalachian Mountains and areas close to the Laurentide Ice Sheet. Many species exhibit distinct patterns that reflect the unique, rather than the shared, aspects of species' phylogeographical histories. Furthermore, similar modern phylogeographical patterns can result from different underlying causal factors operating at different times (i.e. pseudocongruence). One underemphasized component of pseudocongruence may result from the efforts of researchers to categorize patterns visually - similar patterns may, in fact, not fully coincide, and inferring agreement may obscure the actual patterns and lead to erroneous conclusions. Our modelling analyses indicate no clear spatial patterning and support the hypothesis that phylogeographical structure in diverse temperate taxa is complex and was not shaped by just a few barriers.

Phylogeography and population structure of thornback rays (Raja clavata L., Rajidae)

The phylogeography of thornback rays (Raja clavata) was assessed from European waters, using five nuclear microsatellite loci and mitochondrial cytochome b sequences. Strong regional differentiation was found between the Mediterranean basin, the Azores and the European continental shelf. Allelic and haplotype diversities were high in Portuguese populations, consistent with the existence of a refugium along the Iberian Peninsula. Unexpectedly, high diversity was also found in the English Channel/North Sea area. The lowest genetic diversity was found in the Black Sea. Populations sampled from the Mediterranean, Adriatic and Black Seas were characterized by a single mitochondrial haplotype. This haplotype was also the most ancestral and widespread outside of the Mediterranean basin except for the Azores. Populations from the Azores were dominated by a second ancestral haplotype which was shared with British populations. Results from multidimensional scaling, amova and nested clade analysis indicate that British waters are a secondary contact zone recolonized from at least two refugia--one around the Iberian Peninsula and one possibly in the Azores. Links to a potential refugium known as the Hurd Deep, between Cornwall and Brittany, are discussed. Finally, a historical demographic analysis indicates that thornback ray populations started to expand between 580,000 and 362,000 years ago, which suggests that the Last Glacial Maximum (20,000 years ago) had mainly affected the distribution of populations rather than population size.

Global phylogeography of the scalloped hammerhead shark (Sphyrna lewini)

Large marine fishes typically have little population genetic structure. The exceptions are associated with sedentary behaviour, disjunct distributions, or reproductive philopatry. Scalloped hammerhead sharks (Sphyrna lewini) incorporate the contrasting traits of oceanic habitat (usually associated with high dispersal) and possible fidelity to nursery grounds (for reproductive females). To evaluate the expectations of these contrasting behaviours, we examined the global genetic structure of S. lewini based on collections (n = 271 individuals) from 20 nursery areas. A 548-bp fragment of mitochondrial DNA control region revealed 22 polymorphic sites, 24 haplotypes, and three lineages distinguished by 2.56-3.77% sequence divergence. Coalescence analyses based on a provisional molecular clock indicate an origin in the Indo-West Pacific with late Pleistocene radiations into the central Pacific (Hawaii) and eastern Pacific (Central America), as well as recent interchange between oceans via southern Africa. Population subdivisions are strong (overall Phi(ST) = 0.749, P < 0.0001 and among oceans Phi(ST) = 0.598, P < 0.0098). Genetic discontinuity within oceans (Phi(ST) = 0.519, P < 0.0001) is primarily associated with oceanic barriers (migration across oceans M approximately 0), with much less structure along continental margins (M > 10). We conclude that nursery populations linked by continuous coastline have high connectivity, but that oceanic dispersal by females is rare. Although we cannot rule out philopatry to natal nurseries, oceanic barriers appear to have a much stronger influence on the genetic architecture of this species and may indicate a mechanism for recent evolutionary radiations in the genus Sphyrna.

Minding the gap: frequency of indels in mtDNA control region sequence data and influence on population genetic analyses

Insertions and deletions (indels) result in sequences of various lengths when homologous gene regions are compared among individuals or species. Although indels are typically phylogenetically informative, occurrence and incorporation of these characters as gaps in intraspecific population genetic data sets are rarely discussed. Moreover, the impact of gaps on estimates of fixation indices, such as F(ST), has not been reviewed. Here, I summarize the occurrence and population genetic signal of indels among 60 published studies that involved alignments of multiple sequences from the mitochondrial DNA (mtDNA) control region of vertebrate taxa. Among 30 studies observing indels, an average of 12% of both variable and parsimony-informative sites were composed of these sites. There was no consistent trend between levels of population differentiation and the number of gap characters in a data block. Across all studies, the average influence on estimates of PhiST was small, explaining only an additional 1.8% of among population variance (range 0.0-8.0%). Studies most likely to observe an increase in PhiST with the inclusion of gap characters were those with < 20 variable sites, but a near equal number of studies with few variable sites did not show an increase. In contrast to studies at interspecific levels, the influence of indels for intraspecific population genetic analyses of control region DNA appears small, dependent upon total number of variable sites in the data block, and related to species-specific characteristics and the spatial distribution of mtDNA lineages that contain indels.