Global collision-risk hotspots of marine traffic and the world's largest fish, the whale shark.

Marine traffic is increasing globally yet collisions with endangered megafauna such as whales, sea turtles, and planktivorous sharks go largely undetected or unreported. Collisions leading to mortality can have population-level consequences for endangered species. Hence, identifying simultaneous space use of megafauna and shipping throughout ranges may reveal as-yet-unknown spatial targets requiring conservation. However, global studies tracking megafauna and shipping occurrences are lacking. Here we combine satellite-tracked movements of the whale shark, Rhincodon typus, and vessel activity to show that 92% of sharks' horizontal space use and nearly 50% of vertical space use overlap with persistent large vessel (>300 gross tons) traffic. Collision-risk estimates correlated with reported whale shark mortality from ship strikes, indicating higher mortality in areas with greatest overlap. Hotspots of potential collision risk were evident in all major oceans, predominantly from overlap with cargo and tanker vessels, and were concentrated in gulf regions, where dense traffic co-occurred with seasonal shark movements. Nearly a third of whale shark hotspots overlapped with the highest collision-risk areas, with the last known locations of tracked sharks coinciding with busier shipping routes more often than expected. Depth-recording tags provided evidence for sinking, likely dead, whale sharks, suggesting substantial "cryptic" lethal ship strikes are possible, which could explain why whale shark population declines continue despite international protection and low fishing-induced mortality. Mitigation measures to reduce ship-strike risk should be considered to conserve this species and other ocean giants that are likely experiencing similar impacts from growing global vessel traffic.

Seasonally mediated niche partitioning in a vertically compressed pelagic predator guild.

Niche partitioning among closely related, sympatric species is a fundamental concept in ecology, and its mechanisms are of broad interest for understanding ecosystem functioning and predicting the impacts of human-driven environmental change. However, identifying mechanisms by which top marine predators partition available resources has been especially challenging given the difficulty of quantifying resource use of large pelagic animals. In the eastern tropical Pacific (ETP), three large, highly mobile and ecologically similar pelagic predators (blue marlin (Makaira nigricans), black marlin (Istiompax indica) and sailfish (Istiophorus platypterus)) coexist in a vertically compressed habitat. To evaluate each species' ecological niche, we leveraged a decade of recreational fisheries data, multi-year satellite tracking with high-resolution dive data, and stable isotope analysis. Fishery interaction and telemetry-based three-dimensional seasonal utilization distributions suggested high spatial and temporal overlap among species; however, seasonal and diel variability in diving behaviour produced spatial partitioning, leading to low trophic overlap among species. Expanding oxygen minimum zones will reduce the available vertical habitat within predator guilds, likely leading to increases in interspecific competition. Thus, understanding the mechanisms of habitat partitioning among predators in the vertically compressed ETP can provide insight into how predators in other ocean regions may respond to vertically limited habitats.

Patrolling the border: Billfish exploit the hypoxic boundary created by the world's largest oxygen minimum zone.

Pelagic predators must contend with low prey densities that are irregularly distributed and dynamic in space and time. Based on satellite imagery and telemetry data, many pelagic predators will concentrate horizontal movements on ephemeral surface fronts-gradients between water masses-because of enhanced local productivity and increased forage fish densities. Vertical fronts (e.g. thermoclines, oxyclines) can be spatially and temporally persistent, and aggregate lower trophic level and diel vertically migrating organisms due to sharp changes in temperature, water density or available oxygen. Thus, vertical fronts represent a stable and potentially energy rich habitat feature for diving pelagic predators but remain little explored in their capacity to enhance foraging opportunities. Here, we use a novel suite of high-resolution biologging data, including in situ derived oxygen saturation and video, to document how two top predators in the pelagic ecosystem exploit the vertical fronts created by the oxygen minimum zone of the eastern tropical Pacific. Prey search behaviour was dependent on dive shape, and significantly increased near the thermocline and hypoxic boundary for blue marlin Makaira nigricans and sailfish Istiophorus platypterus, respectively. Further, we identify a behaviour not yet reported for pelagic predators, whereby the predator repeatedly dives below the thermocline and hypoxic boundary (and by extension, below the prey). We hypothesize this behaviour is used to ambush prey concentrated at the boundaries from below. We describe how habitat fronts created by low oxygen environments can influence pelagic ecosystems, which will become increasingly important to understand in the context of global change and expanding oxygen minimum zones. We anticipate that our findings are shared among many pelagic predators where strong vertical fronts occur, and additional high-resolution tagging is warranted to confirm this.

Elucidating the role of competition in driving spatial and trophic niche patterns in sympatric juvenile sharks.

The coexistence of ecologically and morphologically similar species is often facilitated by the partitioning of ecological niches. While subordinate species can reduce competition with dominant competitors through spatial and/or trophic segregation, empirical support from wild settings, particularly those involving large-bodied taxa in marine ecosystems, are rare. Shark nursery areas provide an opportunity to investigate the mechanisms of coexistence. We used experimental and field studies of sympatric juvenile sharks (blacktip reef shark, Carcharhinus melanopterus; sicklefin lemon shark, Negaprion acutidens) to investigate how competitive ability influenced realized niches at St. Joseph Atoll, Seychelles. Captive trials revealed that sicklefin lemon sharks were dominant over blacktip reef sharks, consistently taking food rewards. In the field, blacktip reef sharks were captured over a broader area than sicklefin lemon sharks, but daily space use of actively tracked sharks showed a high degree of overlap across microhabitats. While stomach contents analysis revealed that blacktip reef shark diets included a broader range of prey items, stable isotope analysis demonstrated significantly higher mean δ13C values for sicklefin lemon sharks, suggesting diverging dietary preferences. Overall, our results matched theoretical predictions of subordinate competitors using a greater range of habitats and displaying broader feeding niches than competitively dominant species. While separating the realized and fundamental niche of marine predators is complicated, we provide evidence that resource partitioning is at least partially driven by interspecific competition.

White shark genome reveals ancient elasmobranch adaptations associated with wound healing and the maintenance of genome stability.

The white shark (Carcharodon carcharias; Chondrichthyes, Elasmobranchii) is one of the most publicly recognized marine animals. Here we report the genome sequence of the white shark and comparative evolutionary genomic analyses to the chondrichthyans, whale shark (Elasmobranchii) and elephant shark (Holocephali), as well as various vertebrates. The 4.63-Gbp white shark genome contains 24,520 predicted genes, and has a repeat content of 58.5%. We provide evidence for a history of positive selection and gene-content enrichments regarding important genome stability-related genes and functional categories, particularly so for the two elasmobranchs. We hypothesize that the molecular adaptive emphasis on genome stability in white and whale sharks may reflect the combined selective pressure of large genome sizes, high repeat content, high long-interspersed element retrotransposon representation, large body size, and long lifespans, represented across these two species. Molecular adaptation for wound healing was also evident, with positive selection in key genes involved in the wound-healing process, as well as Gene Ontology enrichments in fundamental wound-healing pathways. Sharks, particularly apex predators such as the white shark, are believed to have an acute sense of smell. However, we found very few olfactory receptor genes, very few trace amine-associated receptors, and extremely low numbers of G protein-coupled receptors. We did however, identify 13 copies of vomeronasal type 2 (V2R) genes in white shark and 10 in whale shark; this, combined with the over 30 V2Rs reported previously for elephant shark, suggests this gene family may underlie the keen odorant reception of chondrichthyans.

Genomic Assessment of Global Population Structure in a Highly Migratory and Habitat Versatile Apex Predator, the Tiger Shark (Galeocerdo cuvier).

Understanding the population dynamics of highly mobile, widely distributed, oceanic sharks, many of which are overexploited, is necessary to aid their conservation management. We investigated the global population genomics of tiger sharks (Galeocerdo cuvier), a circumglobally distributed, apex predator displaying remarkable behavioral versatility in its diet, habitat use (near coastal, coral reef, pelagic), and individual movement patterns (spatially resident to long-distance migrations). We genotyped 242 tiger sharks from 10 globally distributed locations at more than 2000 single nucleotide polymorphisms. Although this species often conducts massive distance migrations, the data show strong genetic differentiation at both neutral (FST = 0.125-0.144) and candidate outlier loci (FST = 0.570-0.761) between western Atlantic and Indo-Pacific sharks, suggesting the potential for adaptation to the environments specific to these oceanic regions. Within these regions, there was mixed support for population differentiation between northern and southern hemispheres in the western Atlantic, and none for structure within the Indian Ocean. Notably, the results demonstrate a low level of population differentiation of tiger sharks from the remote Hawaiian archipelago compared with sharks from the Indian Ocean (FST = 0.003-0.005, P < 0.01). Given concerns about biodiversity loss and marine ecosystem impacts caused by overfishing of oceanic sharks in the midst of rapid environmental change, our results suggest it imperative that international fishery management prioritize conservation of the evolutionary potential of the highly genetically differentiated Atlantic and Indo-Pacific populations of this unique apex predator. Furthermore, we suggest targeted management attention to tiger sharks in the Hawaiian archipelago based on a precautionary biodiversity conservation perspective.

Transcriptome-Derived Microsatellites Demonstrate Strong Genetic Differentiation in Pacific White Sharks.

Recent advances in genome-scale sequencing technology have allowed the development of high resolution genetic markers for the study of nonmodel taxa. In particular, transcriptome sequencing has proven to be highly useful in generating genomic markers for use in population genetic studies, allowing for insight into species connectivity, as well as local adaptive processes as many transcriptome-derived markers are found within or associated with functional genes. Herein, we developed a set of 30 microsatellite markers from a heart transcriptome for the white shark (Carcharodon carcharias), a widely distributed and globally vulnerable marine predator. Using these markers as well as 10 published anonymous genomic microsatellite loci, we provide 1) the first nuclear genetic assessment of the cross-Pacific connectivity of white sharks, and 2) a comparison of the levels of inferred differentiation across microsatellite marker sets (i.e., transcriptome vs. anonymous) to assess their respective utility to elucidate the population genetic dynamics of white sharks. Significant (FST = 0.083, P = 0.05; G″ST = 0.200; P = 0.001) genetic differentiation was found between Southwestern Pacific (n = 19) and Northeastern Pacific (n = 20) white sharks, indicating restricted, cross Pacific gene flow in this species. Transcriptome-derived microsatellite marker sets identified much higher (up to 2×) levels of genetic differentiation than anonymous genomic markers, underscoring potential utility of transcriptome markers in identifying subtle population genetic differences within highly vagile, globally distributed marine species.Subject areas: Population structure and phylogeography; Conservation genetics and biodiversity.

Comparative transcriptomics of elasmobranchs and teleosts highlight important processes in adaptive immunity and regional endothermy.

BACKGROUND: Comparative genomic and/or transcriptomic analyses involving elasmobranchs remain limited, with genome level comparisons of the elasmobranch immune system to that of higher vertebrates, non-existent. This paper reports a comparative RNA-seq analysis of heart tissue from seven species, including four elasmobranchs and three teleosts, focusing on immunity, but concomitantly seeking to identify genetic similarities shared by the two lamnid sharks and the single billfish in our study, which could be linked to convergent evolution of regional endothermy. RESULTS: Across seven species, we identified an average of 10,877 Swiss-Prot annotated genes from an average of 32,474 open reading frames within each species' heart transcriptome. About half of these genes were shared between all species while the remainder included functional differences between our groups of interest (elasmobranch vs. teleost and endotherms vs. ectotherms) as revealed by Gene Ontology (GO) and selection analyses. A repeatedly represented functional category, in both the uniquely expressed elasmobranch genes (total of 259) and the elasmobranch GO enrichment results, involved antibody-mediated immunity, either in the recruitment of immune cells (Fc receptors) or in antigen presentation, including such terms as "antigen processing and presentation of exogenous peptide antigen via MHC class II", and such genes as MHC class II, HLA-DPB1. Molecular adaptation analyses identified three genes in elasmobranchs with a history of positive selection, including legumain (LGMN), a gene with roles in both innate and adaptive immunity including producing antigens for presentation by MHC class II. Comparisons between the endothermic and ectothermic species revealed an enrichment of GO terms associated with cardiac muscle contraction in endotherms, with 19 genes expressed solely in endotherms, several of which have significant roles in lipid and fat metabolism. CONCLUSIONS: This collective comparative evidence provides the first multi-taxa transcriptomic-based perspective on differences between elasmobranchs and teleosts, and suggests various unique features associated with the adaptive immune system of elasmobranchs, pointing in particular to the potential importance of MHC Class II. This in turn suggests that expanded comparative work involving additional tissues, as well as genome sequencing of multiple elasmobranch species would be productive in elucidating the regulatory and genome architectural hallmarks of elasmobranchs.

Global population genetic dynamics of a highly migratory, apex predator shark.

Knowledge of genetic connectivity dynamics in the world's large-bodied, highly migratory, apex predator sharks across their global ranges is limited. One such species, the tiger shark (Galeocerdo cuvier), occurs worldwide in warm temperate and tropical waters, uses remarkably diverse habitats (nearshore to pelagic) and possesses a generalist diet that can structure marine ecosystems through top-down processes. We investigated the phylogeography and the global population structure of this exploited, phylogenetically enigmatic shark by using 10 nuclear microsatellites (n = 380) and sequences from the mitochondrial control region (CR, n = 340) and cytochrome oxidase I gene (n = 100). All three marker classes showed the genetic differentiation between tiger sharks from the western Atlantic and Indo-Pacific ocean basins (microsatellite FST  > 0.129; CR ΦST  > 0.497), the presence of North vs. southwestern Atlantic differentiation and the isolation of tiger sharks sampled from Hawaii from other surveyed locations. Furthermore, mitochondrial DNA revealed high levels of intraocean basin matrilineal population structure, suggesting female philopatry and sex-biased gene flow. Coalescent- and genetic distance-based estimates of divergence from CR sequences were largely congruent (dcorr  = 0.0015-0.0050), indicating a separation of Indo-Pacific and western Atlantic tiger sharks <1 million years ago. Mitochondrial haplotype relationships suggested that the western South Atlantic Ocean was likely a historical connection for interocean basin linkages via the dispersal around South Africa. Together, the results reveal unexpectedly high levels of population structure in a highly migratory, behaviourally generalist, cosmopolitan ocean predator, calling for management and conservation on smaller-than-anticipated spatial scales.

Genetic Diversity of White Sharks, Carcharodon carcharias, in the Northwest Atlantic and Southern Africa.

The white shark, Carcharodon carcharias, is both one of the largest apex predators in the world and among the most heavily protected marine fish. Population genetic diversity is in part shaped by recent demographic history and can thus provide information complementary to more traditional population assessments, which are difficult to obtain for white sharks and have at times been controversial. Here, we use the mitochondrial control region and 14 nuclear-encoded microsatellite loci to assess white shark genetic diversity in 2 regions: the Northwest Atlantic (NWA, N = 35) and southern Africa (SA, N = 131). We find that these 2 regions harbor genetically distinct white shark populations (Φ ST = 0.10, P < 0.00001; microsatellite F ST = 0.1057, P < 0.021). M-ratios were low and indicative of a genetic bottleneck in the NWA (M-ratio = 0.71, P < 0.004) but not SA (M-ratio = 0.85, P = 0.39). This is consistent with other evidence showing a steep population decline occurring in the mid to late 20th century in the NWA, whereas the SA population appears to have been relatively stable. Estimates of effective population size ranged from 22.6 to 66.3 (NWA) and 188 to 1998.3 (SA) and evidence of inbreeding was found (primarily in NWA). Overall, our findings indicate that white population dynamics within NWA and SA are determined more by intrinsic reproduction than immigration and there is genetic evidence of a population decline in the NWA, further justifying the strong domestic protective measures that have been taken for this species in this region. Our study also highlights how assessment of genetic diversity can complement other sources of information to better understand the status of threatened marine fish populations.

Comparative population genetics and evolutionary history of two commonly misidentified billfishes of management and conservation concern.

BACKGROUND: Misidentifications between exploited species may lead to inaccuracies in population assessments, with potentially irreversible conservation ramifications if overexploitation of either species is occurring. A notable showcase is provided by the realization that the roundscale spearfish (Tetrapturus georgii), a recently validated species, has been historically misidentified as the morphologically very similar and severely overfished white marlin (Kajikia albida) (IUCN listing: Vulnerable). In effect, no information exists on the population status and evolutionary history of the enigmatic roundscale spearfish, a large, highly vagile and broadly distributed pelagic species. We provide the first population genetic evaluation of the roundscale spearfish, utilizing nuclear microsatellite and mitochondrial DNA sequence markers. Furthermore, we re-evaluated existing white marlin mitochondrial genetic data and present our findings in a comparative context to the roundscale spearfish. RESULTS: Microsatellite and mitochondrial (control region) DNA markers provided mixed evidence for roundscale spearfish population differentiation between the western north and south Atlantic regions, depending on marker-statistical analysis combination used. Mitochondrial DNA analyses provided strong signals of historical population growth for both white marlin and roundscale spearfish, but higher genetic diversity and effective female population size (1.5-1.9X) for white marlin. CONCLUSIONS: The equivocal indications of roundscale spearfish population structure, combined with a smaller effective female population size compared to the white marlin, already a species of concern, suggests that a species-specific and precautionary management strategy recognizing two management units is prudent for this newly validated billfish.

Characterization of the heart transcriptome of the white shark (Carcharodon carcharias).

BACKGROUND: The white shark (Carcharodon carcharias) is a globally distributed, apex predator possessing physical, physiological, and behavioral traits that have garnered it significant public attention. In addition to interest in the genetic basis of its form and function, as a representative of the oldest extant jawed vertebrate lineage, white sharks are also of conservation concern due to their small population size and threat from overfishing. Despite this, surprisingly little is known about the biology of white sharks, and genomic resources are unavailable. To address this deficit, we combined Roche-454 and Illumina sequencing technologies to characterize the first transciptome of any tissue for this species. RESULTS: From white shark heart cDNA we generated 665,399 Roche 454 reads (median length 387-bp) that were assembled into 141,626 contigs (mean length 503-bp). We also generated 78,566,588 Illumina reads, which we aligned to the 454 contigs producing 105,014 454/Illumina consensus sequences. To these, we added 3,432 non-singleton 454 contigs. By comparing these sequences to the UniProtKB/Swiss-Prot database we were able to annotate 21,019 translated open reading frames (ORFs) of ≥ 20 amino acids. Of these, 19,277 were additionally assigned Gene Ontology (GO) functional annotations. While acknowledging the limitations of our single tissue transcriptome, Fisher tests showed the white shark transcriptome to be significantly enriched for numerous metabolic GO terms compared to the zebra fish and human transcriptomes, with white shark showing more similarity to human than to zebra fish (i.e. fewer terms were significantly different). We also compared the transcriptome to other available elasmobranch sequences, for signatures of positive selection and identified several genes of putative adaptive significance on the white shark lineage. The white shark transcriptome also contained 8,404 microsatellites (dinucleotide, trinucleotide, or tetranucleotide motifs ≥ five perfect repeats). Detailed characterization of these microsatellites showed that ORFs with trinucleotide repeats, were significantly enriched for transcription regulatory roles and that trinucleotide frequency within ORFs was lower than for a wide range of taxonomic groups including other vertebrates. CONCLUSION: The white shark heart transcriptome represents a valuable resource for future elasmobranch functional and comparative genomic studies, as well as for population and other biological studies vital for effective conservation of this globally vulnerable species.

The behavioural and genetic mating system of the sand tiger shark, Carcharias taurus, an intrauterine cannibal.

Sand tiger sharks (Carcharias taurus) have an unusual mode of reproduction, whereby the first embryos in each of the paired uteri to reach a certain size ('hatchlings') consume all of their smaller siblings during gestation ('embryonic cannibalism' or EC). If females commonly mate with multiple males ('behavioural polyandry') then litters could initially have multiple sires. It is possible, however, that EC could exclude of all but one of these sires from producing offspring thus influencing the species genetic mating system ('genetic monogamy'). Here, we use microsatellite DNA profiling of mothers and their litters (n = 15, from two to nine embryos per litter) to quantify the frequency of behavioural and genetic polyandry in this system. We conservatively estimate that nine of the females we examined (60%) were behaviourally polyandrous. The genetic mating system was characterized by assessing sibling relationships between hatchlings and revealed only 40 per cent genetic polyandry (i.e. hatchlings were full siblings in 60% of litters). The discrepancy stemmed from three females that were initially fertilized by multiple males but only produced hatchlings with one of them. This reveals that males can be excluded even after fertilizing ova and that some instances of genetic monogamy in this population arise from the reduction in litter size by EC. More research is needed on how cryptic post-copulatory and post-zygotic processes contribute to determining paternity and bridging the behavioural and genetic mating systems of viviparous species.

Hunting behavior of a solitary sailfish Istiophorus platypterus and estimated energy gain after prey capture.

Foraging behavior and interaction with prey is an integral component of the ecological niche of predators but is inherently difficult to observe for highly mobile animals in the marine environment. Billfishes have been described as energy speculators, expending a large amount of energy foraging, expecting to offset high costs with periodic high energetic gain. Surface-based group feeding of sailfish, Istiophorus platypterus, is commonly observed, yet sailfish are believed to be largely solitary roaming predators with high metabolic requirements, suggesting that individual foraging also represents a major component of predator-prey interactions. Here, we use biologging data and video to examine daily activity levels and foraging behavior, estimate metabolic costs, and document a solitary predation event for a 40 kg sailfish. We estimate a median active metabolic rate of 218.9 ± 70.5 mgO2 kg-1 h-1 which increased to 518.8 ± 586.3 mgO2 kg-1 h-1 during prey pursuit. Assuming a successful predation, we estimate a daily net energy gain of 2.4 MJ (5.1 MJ acquired, 2.7 MJ expended), supporting the energy speculator model. While group hunting may be a common activity used by sailfish to acquire energy, our calculations indicate that opportunistic individual foraging events offer a net energy return that contributes to the fitness of these highly mobile predators.

Genomes of endangered great hammerhead and shortfin mako sharks reveal historic population declines and high levels of inbreeding in great hammerhead.

Despite increasing threats of extinction to Elasmobranchii (sharks and rays), whole genome-based conservation insights are lacking. Here, we present chromosome-level genome assemblies for the Critically Endangered great hammerhead (Sphyrna mokarran) and the Endangered shortfin mako (Isurus oxyrinchus) sharks, with genetic diversity and historical demographic comparisons to other shark species. The great hammerhead exhibited low genetic variation, with 8.7% of the 2.77 Gbp genome in runs of homozygosity (ROH) > 1 Mbp and 74.4% in ROH >100 kbp. The 4.98 Gbp shortfin mako genome had considerably greater diversity and <1% in ROH > 1 Mbp. Both these sharks experienced precipitous declines in effective population size (Ne) over the last 250 thousand years. While shortfin mako exhibited a large historical Ne that may have enabled the retention of higher genetic variation, the genomic data suggest a possibly more concerning picture for the great hammerhead, and a need for evaluation with additional individuals.

Three complete mitochondrial genomes of shortfin mako sharks, Isurus oxyrinchus, from the Atlantic and Pacific Oceans.

We present complete mitogenome sequences of three shortfin mako sharks (Isurus oxyrinchus) sampled from the western Pacific, and eastern and western Atlantic oceans. Mitogenome sequence lengths ranged between 16,699 bp and 16,702 bp, and all three mitogenomes contained one non-coding control region, two rRNA genes, 22 tRNA genes, and 13 protein-coding genes. Comparative assessment of five mitogenomes from globally distributed shortfin makos (the current three and two previously published mitogenomes) yielded 98.4% identity, with the protein-coding genes ATP8, ATP6, and ND5 as the most variable regions (sequence identities of 96.4%, 96.5%, and 97.6%, respectively). These mitogenome sequences contribute resources for assessing the genetic population dynamics of this endangered oceanic apex predator.

Genetic population dynamics of the critically endangered scalloped hammerhead shark (Sphyrna lewini) in the Eastern Tropical Pacific.

The scalloped hammerhead shark, Sphyrna lewini, is a Critically Endangered, migratory species known for its tendency to form iconic and visually spectacular large aggregations. Herein, we investigated the population genetic dynamics of the scalloped hammerhead across much of its distribution in the Eastern Tropical Pacific (ETP), ranging from Costa Rica to Ecuador, focusing on young-of-year animals from putative coastal nursery areas and adult females from seasonal aggregations that form in the northern Galápagos Islands. Nuclear microsatellites and partial mitochondrial control region sequences showed little evidence of population structure suggesting that scalloped hammerheads in this ETP region comprise a single genetic stock. Galápagos aggregations of adults were not comprised of related individuals, suggesting that kinship does not play a role in the formation of the repeated, annual gatherings at these remote offshore locations. Despite high levels of fisheries exploitation of this species in the ETP, the adult scalloped hammerheads here showed greater genetic diversity compared with adult conspecifics from other parts of the species' global distribution. A phylogeographic analysis of available, globally sourced, mitochondrial control region sequence data (n = 1818 sequences) revealed that scalloped hammerheads comprise three distinct matrilines corresponding to the three major world ocean basins, highlighting the need for conservation of these evolutionarily unique lineages. This study provides the first view of the genetic properties of a scalloped hammerhead aggregation, and the largest sample size-based investigation of population structure and phylogeography of this species in the ETP to date.

Shark virgin birth produces multiple, viable offspring.

Facultative automictic parthenogenesis has only recently been confirmed in the most ancient jawed vertebrates, the chondrichthyan fishes (sharks, batoids, and chimeras). To date, however, in both documented cases, the females have only produced a single parthenogen offspring, and none of these have lived for more than 3 days. We present a genetically verified case of automictic parthenogenesis by a white-spotted bamboo shark (Chiloscyllium plagiosum), in which at least 2 parthenogens were produced and survived for 5 years or more. These findings demonstrate that some female sharks are capable of producing, multiple, viable offspring through parthenogenesis.

The complete mitochondrial genome of a gray reef shark, Carcharhinus amblyrhynchos (Carcharhiniformes: Carcharhinidae), from the Western Indian Ocean.

We present the mitochondrial genome sequence of a gray reef shark, Carcharhinus amblyrhynchos (Bleeker 1856), a coral reef associated species. This is the first mitogenome for this species from the western Indian Ocean. The mitogenome is 16,705 bp in length, has 13 protein-coding genes, 22 tRNA genes, 2 rRNA genes and a non-coding control region, and demonstrates a gene arrangement congruent with other shark and most vertebrate species. This mitogenome provides a genomic resource for assisting with population, evolutionary and conservation studies for the gray reef shark, which is increasingly under threat from fisheries.

Spatio-temporal genetic tagging of a cosmopolitan planktivorous shark provides insight to gene flow, temporal variation and site-specific re-encounters.

Migratory movements in response to seasonal resources often influence population structure and dynamics. Yet in mobile marine predators, population genetic consequences of such repetitious behaviour remain inaccessible without comprehensive sampling strategies. Temporal genetic sampling of seasonally recurring aggregations of planktivorous basking sharks, Cetorhinus maximus, in the Northeast Atlantic (NEA) affords an opportunity to resolve individual re-encounters at key sites with population connectivity and patterns of relatedness. Genetic tagging (19 microsatellites) revealed 18% of re-sampled individuals in the NEA demonstrated inter/multi-annual site-specific re-encounters. High genetic connectivity and migration between aggregation sites indicate the Irish Sea as an important movement corridor, with a contemporary effective population estimate (Ne) of 382 (CI = 241-830). We contrast the prevailing view of high gene flow across oceanic regions with evidence of population structure within the NEA, with early-season sharks off southwest Ireland possibly representing genetically distinct migrants. Finally, we found basking sharks surfacing together in the NEA are on average more related than expected by chance, suggesting a genetic consequence of, or a potential mechanism maintaining, site-specific re-encounters. Long-term temporal genetic monitoring is paramount in determining future viability of cosmopolitan marine species, identifying genetic units for conservation management, and for understanding aggregation structure and dynamics.