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

Characterization of the Complete Mitochondrial Genome of the Elongate Loach and Its Phylogenetic Implications in Cobitidae

The elongate loach is an endemic fish in China. Previous studies have provided some insights into the mitochondrial genome composition and the phylogenetic relationships of the elongate loach inferred using protein-coding genes (PCGs), yet detailed information about it remains limited. Therefore, in this study we sequenced the complete mitochondrial genome of the elongate loach and analyzed its structural characteristics. The PCGs and mitochondrial genome were used for selective stress analysis and genomic comparative analysis. The complete mitochondrial genome of the elongate loach, together with those of 35 Cyprinidae species, was used to infer the phylogenetic relationships of the Cobitidae family through maximum likelihood (ML) reconstruction. The results showed that the genome sequence has a full length of 16,591 bp, which includes 13 PCGs, 22 transfer RNA genes (tRNA), 2 ribosomal RNA genes (rRNA), and 2 non-coding regions (CR D-loop and light chain sub-chain replication origin OL). Overall, the elongate loach shared the same gene arrangement and composition of the mitochondrial genes with other teleost fishes. The Ka/Ks ratios of all mitochondrial PCGs were less than 1, indicating that all of the PCGs were evolving under purifying selection. Genome comparison analyses showed a significant sequence homology of species of Leptobotia. A significant identity between L. elongata and the other five Leptobotia species was observed in the visualization result, except for L. mantschurica, which lacked the tRNA-Arg gene and had a shorter tRNA-Asp gene. The phylogenetic tree revealed that the Cobitidae species examined here can be grouped into two clades, with the elongate loach forming a sister relationship with L. microphthalma. This study could provide additional inferences for a better understanding of the phylogenetic relationships among Cobitidae species.

Genetic stock structure of the silky shark Carcharhinus falciformis in the Indo-Pacific Ocean

The silky shark, Carcharhinus falciformis, is a cosmopolitan species commonly caught as a bycatch for longline fisheries. However, the genetic stock structure for the Indo-Pacific Ocean is not well-defined yet. Here, we used eight microsatellite loci to examine the genetic stock structure and effective population size of 307 silky sharks across 5 Indo-Pacific sampling locations. A major genetic break was found between Aceh and the remaining locations (FST = 0.0505-0.0828, p = 0.001). The Indian Ocean displayed a slightly lower effective population estimate (Ne) compared to the Pacific Ocean, potentially due to the higher fishing pressure in the Indian Ocean region. The lowest Ne was found in the Aceh population (Ne = 2.3), suggesting it might be a small and endemic population. These findings offer valuable information for the conservation and management of the silky shark. We suggest that the population around Aceh waters constitutes a distinct stock and should be managed independently. Further investigations into migratory and movement patterns are needed to define the boundaries of different stocks, ensuring effective management the silky shark across the Indo-Pacific region.

Phylogeography of sharks and rays: a global review based on life history traits and biogeographic partitions

Considerable research exists on the life history traits, evolutionary history, and environmental factors that shape the population genetic structure of marine organisms, including sharks and rays. Conservation concerns are particularly strong for this group as they are highly susceptible to anthropogenic stressors due to a combination of life history traits including late maturity and low fecundity. Here, we provide a review and synthesis of the global phylogeography of sharks and rays. We examined existing data for 40 species of sharks belonging to 17 genera and 19 species of rays belonging to 11 genera. Median joining haplotype networks were constructed for each species for the mtDNA cytochrome C oxidase subunit I (COI), and an Analysis of Molecular Variance (AMOVA) was conducted to understand patterns of genetic diversity and structure across the three major ocean basins-the Indian, Atlantic and Pacific Oceans. Haplotype networks showed very shallow coalescence in most species, a finding previously reported for marine teleosts. Star topologies were predominant among sharks while complex mutational topologies predominated among rays, a finding we attribute to extremely limited dispersal in the early life history of rays. Population structuring varied amongst species groups, apparently due to differences in life history traits including reproductive philopatry, site fidelity, pelagic habitat, migratory habits, and dispersal ability. In comparison to reef-associated and demersal species, pelagic and semi pelagic species showed lower levels of structure between and within ocean basins. As expected, there is variation between taxa and groups, but there are also some broad patterns that can guide management and conservation strategies.

Characterization of mitochondrial genome of Indian Ocean blue-spotted maskray, Neotrygon indica and its phylogenetic relationship within Dasyatidae Family

The present study characterized complete mitochondrial genome of Blue-spotted maskray, Neotrygon indica and studied the evolutionary relationship of the species within the Dasyatidae family. The total length of the mitogenome was 17,974 bp including 37 genes and a non-coding control region. The average frequency of nucleotides in protein-coding genes was A: 29.1 %, T: 30.2 %, G: 13.0 % and C: 27.7 % with AT content of 59.3 %. The values of AT and GC skewness were -0.018 and -0.338, respectively. Comparative analyses showed a large number of average synonymous substitutions per synonymous site (Ks) in gene NADH4 (5.07) followed by NADH5 (4.72). High values of average number of non-synonymous substitutions per non-synonymous site (Ka) were observed in genes ATPase8 (0.54) and NADH2 (0.44). Genes NADH4L and NADH2 showed high interspecific genetic distance values of 0.224 ± 0.001 and 0.213 ± 0.002, respectively. Heat map analysis showed variation in codon usage among different species of the Dasyatidae family. The phylogenetic tree showed a sister relationship between the Dasyatinae and the Neotrygoninae subfamilies. Neotrygon indica formed as a sister species to the clade consisting of N. varidens and N. orientalis. Based on the present results, Neotrygon indica could have diverged from the common ancestor of the two latter in the Plio-Pleistocene. The present study showed distinct characteristics of N. indica from its congeners through comparative mitogenomics.

Mitochondrial marker implies fishery separate management units for spotted sardinella, Amblygaster sirm (Walbaum, 1792) populations in the South China Sea and the Andaman Sea

The spotted sardinella, Amblygaster sirm (Walbaum, 1792), is a commercial sardine commonly caught in Malaysia. Lack of management of these marine species in Malaysian waters could lead to overfishing and potentially declining fish stock populations. Therefore, sustainable management of this species is of paramount importance to ensure its longevity. As such, molecular information is vital in determining the A. sirm population structure and management strategy. In the present study, mitochondrial DNA Cytochrome b was sequenced from 10 A. sirm populations: the Andaman Sea (AS) (two), South China Sea (SCS) (six), Sulu Sea (SS) (one), and Celebes Sea (CS) (one). Accordingly, the intra-population haplotype diversity (Hd) was high (0.91-1.00), and nucleotide diversity (π) was low (0.002-0.009), which suggests a population bottleneck followed by rapid population growth. Based on the phylogenetic trees, minimum spanning network (MSN), population pairwise comparison, and F _ST,and supported by analysis of molecular variance (AMOVA) and spatial analysis of molecular variance (SAMOVA) tests, distinct genetic structures were observed (7.2% to 7.6% genetic divergence) between populations in the SCS and its neighboring waters, versus those in the AS. Furthermore, the results defined A. sirm stock boundaries and evolutionary between the west and east coast (which shares the same waters as western Borneo) of Peninsular Malaysia. In addition, genetic homogeneity was revealed throughout the SCS, SS, and CS based on the non-significant F _STpairwise comparisons. Based on the molecular evidence, separate management strategies may be required for A. sirm of the AS and the SCS, including its neighboring waters.

Little Sharks in a Big World: Mitochondrial DNA Reveals Small-scale Population Structure in the California Horn Shark (Heterodontus francisci)

The California horn shark (Heterodontus francisci) is a small demersal species distributed from southern California and the Channel Islands to Baja California and the Gulf of California. These nocturnal reef predators maintain small home-ranges as adults, and lay auger-shaped egg cases that become wedged into the substrate. While population trends are not well documented, this species is subject to fishing pressure through portions of its range and has been identified as vulnerable to overexploitation. Here we present a survey of 318 specimens from across the range, using mtDNA control region sequences to provide the first genetic assessment of H. francisci. Overall population structure (ΦST = 0.266, P < 0.001) is consistent with limited dispersal as indicated by life history, with two distinct features. Population structure along the continuous coastline is low, with no discernable breaks from Santa Barbara, CA to Bahia Tortugas (Baja California Sur, Mexico); however, there is a notable partition at Punta Eugenia (BCS), a well-known biogeographic break between tropical and subtropical marine faunas. In contrast, population structure is much higher (max ΦST = 0.601, P < 0.05) between the coast and adjacent Channel Islands, a minimum distance of 19 km, indicating that horn sharks rarely disperse across deep habitat and open water. Population structure in most elasmobranchs is measured on a scale of hundreds to thousands of kilometers, but the California Horn Shark has population partitions on an unprecedented small scale, indicating a need for localized management strategies which ensure adequate protection of distinct stocks.

Genomic insights into the historical and contemporary demographics of the grey reef shark

Analyses of genetic diversity can shed light on both the origins of biodiversity hotspots, as well as the conservation status of species that are impacted by human activities. With these objectives, we assembled a genomic dataset of 14,935 single nucleotide polymorphisms from 513 grey reef sharks (Carcharhinus amblyrhynchos) sampled across 17 locations in the tropical Indo-Pacific. We analysed geographic variation in genetic diversity, estimated ancient and contemporary effective population size (N_e) across sampling locations (using coalescent and linkage disequilibrium methods) and modelled the history of gene flow between the Coral Triangle and the Coral Sea. Genetic diversity decreased with distance away from the Coral Triangle and north-western Australia, implying that C. amblyrhynchos may have originated in this region. Increases in N_e were detected across almost all sampling locations 40,000-90,000 generations ago (approximately 0.6-1.5 mya, given an estimated generation time of 16.4 years), suggesting a range expansion around this time. More recent, secondary increases in N_e were inferred for the Misool and North Great Barrier Reef sampling locations, but joint modelling did not clarify whether these were due to population growth, migration, or both. Despite the greater genetic diversity and ancient N_e observed at sites around Australia and the Coral Triangle, remote reefs around north-western New Caledonia had the highest contemporary N_e, demonstrating the importance of using multiple population size assessment methods. This study provides insight into both the past and present demographics of C. amblyrhynchos and contributes to our understanding of evolution in marine biodiversity hotspots.

Seascape genomics of coastal bottlenose dolphins along strong gradients of temperature and salinity

Heterogeneous seascapes and strong environmental gradients in coastal waters are expected to influence adaptive divergence, particularly in species with large population sizes where selection is expected to be highly efficient. However, these influences might also extend to species characterized by strong social structure, natal philopatry and small home ranges. We implemented a seascape genomic study to test this hypothesis in Indo-Pacific bottlenose dolphins (Tursiops aduncus) distributed along the environmentally heterogeneous coast of southern Australia. The datasets included oceanographic and environmental variables thought to be good predictors of local adaptation in dolphins and 8,081 filtered single nucleotide polymorphisms (SNPs) genotyped for individuals sampled from seven different bioregions. From a neutral perspective, population structure and connectivity of the dolphins were generally influenced by habitat type and social structuring. Genotype-environment association analysis identified 241 candidate adaptive loci and revealed that sea surface temperature and salinity gradients influenced adaptive divergence in these animals at both large- (1,000s km) and fine-scales (<100 km). Enrichment analysis and annotation of candidate genes revealed functions related to sodium-activated ion transport, kidney development, adipogenesis and thermogenesis. The findings of spatial adaptive divergence and inferences of putative physiological adaptations challenge previous suggestions that marine megafauna is most likely to be affected by environmental and climatic changes via indirect, trophic effects. Our work contributes to conservation management of coastal bottlenose dolphins subjected to anthropogenic disturbance and to efforts of clarifying how seascape heterogeneity influences adaptive diversity and evolution in small cetaceans.

Brown banded bamboo shark (Chiloscyllium punctatum) shows high genetic diversity and differentiation in Malaysian waters

The demersal brown banded bamboo shark Chiloscyllium punctatum is a major component of sharks landed in Malaysia. However, little is known about their population structure and the effect of high fishing pressure on these weak swimming sharks. Both mitochondrial DNA control region (1072 bp) and NADH dehydrogenase subunit 2 (1044 bp) were used to elucidate the genetic structure and connectivity of C. punctatum among five major areas within the Sundaland region. Our findings revealed (i) strong genetic structure with little present day mixing between the major areas, (ii) high intra-population genetic diversity with unique haplotypes, (iii) significant correlation between genetic differentiation and geographical distance coupled with detectable presence of fine scale geographical barriers (i.e. the South China Sea), (iv) historical directional gene flow from the east coast of Peninsular Malaysia towards the west coast and Borneo, and (v) no detectable genetic differentiation along the coastline of east Peninsular Malaysia. Genetic patterns inferred from the mitochondrial DNA loci were consistent with the strong coastal shelf association in this species, the presence of contemporary barriers shaped by benthic features, and limited current-driven egg dispersal. Fine scale population structure of C. punctatum highlights the need to improve genetic understanding for fishery management and conservation of other small-sized sharks.

Population structure, connectivity, and demographic history of an apex marine predator, the bull shark Carcharhinus leucas

Knowledge of population structure, connectivity, and effective population size remains limited for many marine apex predators, including the bull shark Carcharhinus leucas. This large-bodied coastal shark is distributed worldwide in warm temperate and tropical waters, and uses estuaries and rivers as nurseries. As an apex predator, the bull shark likely plays a vital ecological role within marine food webs, but is at risk due to inshore habitat degradation and various fishing pressures. We investigated the bull shark's global population structure and demographic history by analyzing the genetic diversity of 370 individuals from 11 different locations using 25 microsatellite loci and three mitochondrial genes (CR, nd4, and cytb). Both types of markers revealed clustering between sharks from the Western Atlantic and those from the Western Pacific and the Western Indian Ocean, with no contemporary gene flow. Microsatellite data suggested low differentiation between the Western Indian Ocean and the Western Pacific, but substantial differentiation was found using mitochondrial DNA. Integrating information from both types of markers and using Bayesian computation with a random forest procedure (ABC-RF), this discordance was found to be due to a complete lack of contemporary gene flow. High genetic connectivity was found both within the Western Indian Ocean and within the Western Pacific. In conclusion, these results suggest important structuring of bull shark populations globally with important gene flow occurring along coastlines, highlighting the need for management and conservation plans on regional scales rather than oceanic basin scale.

Review of Current Conservation Genetic Analyses of Northeast Pacific Sharks

Conservation genetics is an applied science that utilizes molecular tools to help solve problems in species conservation and management. It is an interdisciplinary specialty in which scientists apply the study of genetics in conjunction with traditional ecological fieldwork and other techniques to explore molecular variation, population boundaries, and evolutionary relationships with the goal of enabling resource managers to better protect biodiversity and identify unique populations. Several shark species in the northeast Pacific (NEP) have been studied using conservation genetics techniques, which are discussed here. The primary methods employed to study population genetics of sharks have historically been nuclear microsatellites and mitochondrial (mt) DNA. These markers have been used to assess genetic diversity, mating systems, parentage, relatedness, and genetically distinct populations to inform management decisions. Novel approaches in conservation genetics, including next-generation DNA and RNA sequencing, environmental DNA (eDNA), and epigenetics are just beginning to be applied to elasmobranch evolution, physiology, and ecology. Here, we review the methods and results of past studies, explore future directions for shark conservation genetics, and discuss the implications of molecular research and techniques for the long-term management of shark populations in the NEP.

Switch from sexual to parthenogenetic reproduction in a zebra shark

Parthenogenesis is a natural form of asexual reproduction in which embryos develop in the absence of fertilisation. Most commonly found in plants and invertebrate organisms, an increasing number of vertebrate species have recently been reported employing this reproductive strategy. Here we use DNA genotyping to report the first demonstration of an intra-individual switch from sexual to parthenogenetic reproduction in a shark species, the zebra shark Stegostoma fasciatum. A co-housed, sexually produced daughter zebra shark also commenced parthenogenetic reproduction at the onset of maturity without any prior mating. The demonstration of parthenogenesis in these two conspecific individuals with different sexual histories provides further support that elasmobranch fishes may flexibly adapt their reproductive strategy to environmental circumstances.

Complete mitochondrial DNA genome of bonnethead shark, Sphyrna tiburo, and phylogenetic relationships among main superorders of modern elasmobranchs

Elasmobranchs are one of the most diverse groups in the marine realm represented by 18 orders, 55 families and about 1200 species reported, but also one of the most vulnerable to exploitation and to climate change. Phylogenetic relationships among main orders have been controversial since the emergence of the Hypnosqualean hypothesis by Shirai (1992) that considered batoids as a sister group of sharks. The use of the complete mitochondrial DNA (mtDNA) may shed light to further validate this hypothesis by increasing the number of informative characters. We report the mtDNA genome of the bonnethead shark Sphyrna tiburo, and compare it with mitogenomes of other 48 species to assess phylogenetic relationships. The mtDNA genome of S. tiburo, is quite similar in size to that of congeneric species but also similar to the reported mtDNA genome of other Carcharhinidae species. Like most vertebrate mitochondrial genomes, it contained 13 protein coding genes, two rRNA genes and 22 tRNA genes and the control region of 1086 bp (D-loop). The Bayesian analysis of the 49 mitogenomes supported the view that sharks and batoids are separate groups.

Reef sharks: recent advances in ecological understanding to inform conservation

Sharks are increasingly being recognized as important members of coral-reef communities, but their overall conservation status remains uncertain. Nine of the 29 reef-shark species are designated as data deficient in the IUCN Red List, and three-fourths of reef sharks had unknown population trends at the time of their assessment. Fortunately, reef-shark research is on the rise. This new body of research demonstrates reef sharks' high site restriction, fidelity and residency on coral reefs, their broad trophic roles connecting reef communities and their high population genetic structure, all information that should be useful for their management and conservation. Importantly, recent studies on the abundance and population trends of the three classic carcharhinid reef sharks (grey reef shark Carcharhinus amblyrhynchos, blacktip reef shark Carcharhinus melanopterus and whitetip reef shark Triaenodon obesus) may contribute to reassessments identifying them as more vulnerable than currently realized. Because over half of the research effort has focused on only these three reef sharks and the nurse shark Ginglymostoma cirratum in only a few locales, there remain large taxonomic and geographic gaps in reef-shark knowledge. As such, a large portion of reef-shark biodiversity remains uncharacterized despite needs for targeted research identified in their red list assessments. A research agenda for the future should integrate abundance, life history, trophic ecology, genetics, habitat use and movement studies, and expand the breadth of such research to understudied species and localities, in order to better understand the conservation requirements of these species and to motivate effective conservation solutions.

Genetic structure, population demography and seasonal occurrence of blacktip shark Carcharhinus limbatus in Bimini, the Bahamas

A longline survey was conducted from 2004 to 2014 to investigate the demographic population structure and seasonal abundance of the blacktip shark Carcharhinus limbatus in the Bimini Islands, the Bahamas. All individuals sampled (n = 242) were sub-adult or adults [70·1-145·1 cm pre-caudal length (LPC) range] with no neonates or YOY recorded in Bimini. Carcharhinus limbatus abundance peaked in September, coincident with the largest ratio of female to male sharks and a peak in fresh mating wounds on females. Mitochondrial control region (mtCR) DNA sequences were obtained from C. limbatus at Bimini to test whether Bimini C. limbatus are most closely related to geographically proximate populations sampled on the south-eastern coast of the U.S.A., the closest known nursery areas for this species. Nine mtCR haplotypes were observed in 32 individuals sampled at Bimini [haplotype diversity (h) = 0·821, nucleotide diversity (π) = 0·0015]. Four haplotypes observed from Bimini matched those previously found in the northern Yucatan (Mexico)-Belize and two matched a haplotype previously found in the U.S.A. Four haplotypes were novel but were closely related to the northern Yucatan-Belizean haplotypes. Pair-wise ΦST analysis showed that Bimini was significantly differentiated from all of the populations previously sampled (U.S.A. Atlantic, U.S.A. Gulf of Mexico, northern Yucatan, Belize and Brazil). This indicates that C. limbatus sampled from Bimini are unlikely from the described, proximate U.S.A. nurseries.

Population genetics of four heavily exploited shark species around the Arabian Peninsula

The northwestern Indian Ocean harbors a number of larger marine vertebrate taxa that warrant the investigation of genetic population structure given remarkable spatial heterogeneity in biological characteristics such as distribution, behavior, and morphology. Here, we investigate the genetic population structure of four commercially exploited shark species with different biological characteristics (Carcharhinus limbatus, Carcharhinus sorrah, Rhizoprionodon acutus, and Sphyrna lewini) between the Red Sea and all other water bodies surrounding the Arabian Peninsula. To assess intraspecific patterns of connectivity, we constructed statistical parsimony networks among haplotypes and estimated (1) population structure; and (2) time of most recent population expansion, based on mitochondrial control region DNA and a total of 20 microsatellites. Our analysis indicates that, even in smaller, less vagile shark species, there are no contemporary barriers to gene flow across the study region, while historical events, for example, Pleistocene glacial cycles, may have affected connectivity in C. sorrah and R. acutus. A parsimony network analysis provided evidence that Arabian S. lewini may represent a population segment that is distinct from other known stocks in the Indian Ocean, raising a new layer of conservation concern. Our results call for urgent regional cooperation to ensure the sustainable exploitation of sharks in the Arabian region.

Species composition and aspects of the biology of Orectolobiformes from Indonesian waters

The biological aspects of members of the order Orectolobiformes exploited in Indonesian waters are given. Seven species belonging to five families were recorded in the catches at various landing sites in southern Indonesia. Of these, Chiloscyllium punctatum was the most abundant species landed, contributing >50% of the number of orectolobiforms recorded. The biological data obtained varied greatly between the species listed. The total length (L_T ) at maturity for some species, e.g. C. punctatum and Nebrius ferrugineus, varied from that which has been previously recorded for those species in other regions. This highlights the need for regional-specific biological data for fisheries managers and conservation assessors.

Mitogenomics of the Speartooth Shark challenges ten years of control region sequencing

Background

Mitochondrial DNA markers have long been used to identify population boundaries and are now a standard tool in conservation biology. In elasmobranchs, evolutionary rates of mitochondrial genes are low and variation between distinct populations can be hard to detect with commonly used control region sequencing or other single gene approaches. In this study we sequenced the whole mitogenome of 93 Critically Endangered Speartooth Shark Glyphis glyphis from the last three river drainages they inhabit in northern Australia.

Results

Genetic diversity was extremely low (π =0.00019) but sufficient to demonstrate the existence of barriers to gene flow among river drainages (AMOVA Φ_ST =0.28283, P <0.00001). Surprisingly, the comparison with single gene sub-datasets revealed that ND5 and 12S were the only ones carrying enough information to detect similar levels of genetic structure. The control region exhibited only one mutation, which was not sufficient to detect any structure among river drainages.

Conclusions

This study strongly supports the use of single river drainages as discrete management units for the conservation of G. glyphis. Furthermore when genetic diversity is low, as is often the case in elasmobranchs, our results demonstrate a clear advantage of using the whole mitogenome to inform population structure compared to single gene approaches. More specifically, this study questions the extensive use of the control region as the preferential marker for elasmobranch population genetic studies and whole mitogenome sequencing will probably uncover a large amount of cryptic population structure in future studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-014-0232-x) contains supplementary material, which is available to authorized users.

The use of positive reinforcement in training zebra sharks (Stegostoma fasciatum)

Positive reinforcement training (PRT) was used on 4 adult zebra sharks, Stegostoma fasciatum, housed at the Downtown Aquarium, Denver, to determine the ability of zebra sharks to become desensitized to various stimuli associated with veterinary procedures. One male and 3 female sharks were trained for 12 weeks. As a result of PRT, all 4 zebra sharks were desensitized to staying within a closed holding tank off of the main exhibit, the presence of multiple trainers in the closed holding tank, and tactile stimulation. One of the 4 zebra sharks was also successfully desensitized to the presence of a stretcher being brought into the holding tank. All of these procedures are common in veterinary examinations, and it is hoped that desensitization to these stimuli will reduce the stress associated with examinations. The training accomplished has allowed for easier maintenance of the zebra sharks by the aquarium staff and an improvement in the care of the sharks.

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 the Indo-West Pacific maskrays (Dasyatidae, Neotrygon): a complex example of chondrichthyan radiation in the Cenozoic

Maskrays of the genus Neotrygon (Dasyatidae) have dispersed widely in the Indo-West Pacific being represented largely by an assemblage of narrow-ranging coastal endemics. Phylogenetic reconstruction methods reproduced nearly identical and statistically robust topologies supporting the monophyly of the genus Neotrygon within the family Dasyatidae, the genus Taeniura being consistently basal to Neotrygon, and Dasyatis being polyphyletic to the genera Taeniurops and Pteroplatytrygon. The Neotrygon kuhlii complex, once considered to be an assemblage of color variants of the same biological species, is the most derived and widely dispersed subgroup of the genus. Mitochondrial (COI, 16S) and nuclear (RAG1) phylogenies used in synergy with molecular dating identified paleoclimatic fluctuations responsible for periods of vicariance and dispersal promoting population fragmentation and speciation in Neotrygon. Signatures of population differentiation exist in N. ningalooensis and N. annotata, yet a large-scale geological event, such as the collision between the Australian and Eurasian Plates, coupled with subsequent sea-level falls, appears to have separated a once homogeneous population of the ancestral form of N. kuhlii into southern Indian Ocean and northern Pacific taxa some 4–16 million years ago. Repeated climatic oscillations, and the subsequent establishment of land and shallow sea connections within and between Australia and parts of the Indo-Malay Archipelago, have both promoted speciation and established zones of secondary contact within the Indian and Pacific Ocean basins.

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.

Demographic processes underlying subtle patterns of population structure in the scalloped hammerhead shark, Sphyrna lewini

Genetic diversity (θ), effective population size (N(e)), and contemporary levels of gene flow are important parameters to estimate for species of conservation concern, such as the globally endangered scalloped hammerhead shark, Sphyrna lewini. Therefore, we have reconstructed the demographic history of S. lewini across its Eastern Pacific (EP) range by applying classical and coalescent population genetic methods to a combination of 15 microsatellite loci and mtDNA control region sequences. In addition to significant population genetic structure and isolation-by-distance among seven coastal sites between central Mexico and Ecuador, the analyses revealed that all populations have experienced a bottleneck and that all current values of θ are at least an order of magnitude smaller than ancestral θ, indicating large decreases in N(e) (θ = 4N(e)μ), where μ is the mutation rate. Application of the isolation-with-migration (IM) model showed modest but significant genetic connectivity between most sampled sites (point estimates of Nm = 0.1-16.7), with divergence times (t) among all populations significantly greater than zero. Using a conservative (i.e., slow) fossil-based taxon-specific phylogenetic calibration for mtDNA mutation rates, posterior probability distributions (PPDs) for the onset of the decline in N(e) predate modern fishing in this region. The cause of decline over the last several thousand years is unknown but is highly atypical as a post-glacial demographic history. Regardless of the cause, our data and analyses suggest that S. lewini was far more abundant throughout the EP in the past than at present.

Spatially explicit Bayesian clustering models in population genetics

This article reviews recent developments in Bayesian algorithms that explicitly include geographical information in the inference of population structure. Current models substantially differ in their prior distributions and background assumptions, falling into two broad categories: models with or without admixture. To aid users of this new generation of spatially explicit programs, we clarify the assumptions underlying the models, and we test these models in situations where their assumptions are not met. We show that models without admixture are not robust to the inclusion of admixed individuals in the sample, thus providing an incorrect assessment of population genetic structure in many cases. In contrast, admixture models are robust to an absence of admixture in the sample. We also give statistical and conceptual reasons why data should be explored using spatially explicit models that include admixture.

Global population structure of the spiny dogfish Squalus acanthias, a temperate shark with an antitropical distribution

The spiny dogfish (Squalus acanthias) is a temperate, coastal squaloid shark with an antitropical distribution in the Atlantic and Pacific oceans. The global population structure of this species is poorly understood, although individuals are known to undergo extensive migrations within coastal waters and across ocean basins. In this study, an analysis of the global population structure of the spiny dogfish was conducted using eight polymorphic nuclear microsatellite markers and a 566-bp fragment of the mitochondrial ND2 gene region. A low level of genetic divergence was found among collections from the Atlantic and South Pacific basins, whereas a high level of genetic divergence was found among Pacific Ocean collections. Two genetically distinct groups were recovered by both marker classes: one exclusive to North Pacific collections, and one including collections from the South Pacific and Atlantic locations. The strong genetic break across the equatorial Pacific coincides with major regional differences in the life-history characters of spiny dogfish, suggesting that spiny dogfish in areas on either side of the Pacific equator have been evolving independently for a considerable time. Phylogeographic analyses indicate that spiny dogfish populations had a Pacific origin, and that the North Atlantic was colonized as a result of a recent range expansion from the South American coast. Finally, the available data strongly argue for the taxonomic separation of the North Pacific spiny dogfish from S. acanthias and a re-evaluation of the specific status of S. acanthias is warranted.

The extent of population genetic subdivision differs among four co-distributed shark species in the Indo-Australian archipelago

Background

The territorial fishing zones of Australia and Indonesia are contiguous to the north of Australia in the Timor and Arafura Seas and in the Indian Ocean to the north of Christmas Island. The area surrounding the shared boundary consists of a variety of bio-diverse marine habitats including shallow continental shelf waters, oceanic trenches and numerous offshore islands. Both countries exploit a variety of fisheries species, including whaler (Carcharhinus spp.) and hammerhead sharks (Sphyrna spp.). Despite their differences in social and financial arrangements, the two countries are motivated to develop complementary co-management practices to achieve resource sustainability. An essential starting point is knowledge of the degree of population subdivision, and hence fisheries stock status, in exploited species.

Results

Populations of four commercially harvested shark species (Carcharhinus obscurus, Carcharhinus sorrah, Prionace glauca, Sphyrna lewini) were sampled from northern Australia and central Indonesia. Neutral genetic markers (mitochondrial DNA control region sequence and allelic variation at co-dominant microsatellite loci) revealed genetic subdivision between Australian and Indonesian populations of C. sorrah. Further research is needed to address the possibility of genetic subdivision among C. obscurus populations. There was no evidence of genetic subdivision for P. glauca and S. lewini populations, but the sampling represented a relatively small part of their distributional range. For these species, more detailed analyses of population genetic structure is recommended in the future.

Conclusion

Cooperative management between Australia and Indonesia is the best option at present for P. glauca and S. lewini, while C. sorrah and C. obscurus should be managed independently. On-going research on these and other exploited shark and ray species is strongly recommended. Biological and ecological similarity between species may not be a predictor of population genetic structure, so species-specific studies are recommended to provide new data to assist with sustainable fisheries management.