Spiral valve parasites of blue and common thresher sharks as indicators of shark feeding behaviour and ecology.

This study documented the parasite faunas of the spiral valves of blue sharks Prionace glauca (L. 1758) and common thresher sharks Alopias vulpinus (Bonnaterre, 1788) caught in the California Current Large Marine Ecosystem (CCLME) north of the Mexican border. The spiral valves of 18 blue and 19 thresher sharks caught in the CCLME from 2009 to 2013 were examined for parasites. Seven parasite taxa were found in blue sharks and nine in threshers. The tetraphyllidean cestode Anthobothrium sp. (78% prevalence) was the most common parasite in blue sharks, and the phyllobothriid cestode Paraorygmatobothrium sp. (90% prevalence) was the most common in threshers. An adult nematode of the genus Piscicapillaria was found in threshers for the first time and may be a new species. Adult individuals of Hysterothylacium sp. were found in both shark species. The adult acanthocephalan Rhadinorhynchus cololabis and remains of the parasitic copepod Pennella sp. - both parasites of Pacific saury, Cololabis saira - were found in the intestines of threshers, indicating recent feeding on saury. This study paves the way for a more comprehensive examination, including more samples and a wider variety of shark species, to provide a greater understanding of shark feeding behaviour and possibly provide information on shark population biology.

Evidence of positive selection associated with placental loss in tiger sharks.

BACKGROUND: All vertebrates initially feed their offspring using yolk reserves. In some live-bearing species these yolk reserves may be supplemented with extra nutrition via a placenta. Sharks belonging to the Carcharhinidae family are all live-bearing, and with the exception of the tiger shark (Galeocerdo cuvier), develop placental connections after exhausting yolk reserves. Phylogenetic relationships suggest the lack of placenta in tiger sharks is due to secondary loss. This represents a dramatic shift in reproductive strategy, and is likely to have left a molecular footprint of positive selection within the genome. RESULTS: We sequenced the transcriptome of the tiger shark and eight other live-bearing shark species. From this data we constructed a time-calibrated phylogenetic tree estimating the tiger shark lineage diverged from the placental carcharhinids approximately 94 million years ago. Along the tiger shark lineage, we identified five genes exhibiting a signature of positive selection. Four of these genes have functions likely associated with brain development (YWHAE and ARL6IP5) and sexual reproduction (VAMP4 and TCTEX1D2). CONCLUSIONS: Our results indicate the loss of placenta in tiger sharks may be associated with subsequent adaptive changes in brain development and sperm production.

Environmental context explains Lévy and Brownian movement patterns of marine predators.

An optimal search theory, the so-called Lévy-flight foraging hypothesis, predicts that predators should adopt search strategies known as Lévy flights where prey is sparse and distributed unpredictably, but that Brownian movement is sufficiently efficient for locating abundant prey. Empirical studies have generated controversy because the accuracy of statistical methods that have been used to identify Lévy behaviour has recently been questioned. Consequently, whether foragers exhibit Lévy flights in the wild remains unclear. Crucially, moreover, it has not been tested whether observed movement patterns across natural landscapes having different expected resource distributions conform to the theory's central predictions. Here we use maximum-likelihood methods to test for Lévy patterns in relation to environmental gradients in the largest animal movement data set assembled for this purpose. Strong support was found for Lévy search patterns across 14 species of open-ocean predatory fish (sharks, tuna, billfish and ocean sunfish), with some individuals switching between Lévy and Brownian movement as they traversed different habitat types. We tested the spatial occurrence of these two principal patterns and found Lévy behaviour to be associated with less productive waters (sparser prey) and Brownian movements to be associated with productive shelf or convergence-front habitats (abundant prey). These results are consistent with the Lévy-flight foraging hypothesis, supporting the contention that organism search strategies naturally evolved in such a way that they exploit optimal Lévy patterns.

Evaluating the suitability of close-kin mark-recapture as a demographic modelling tool for a critically endangered elasmobranch population.

Estimating the demographic parameters of contemporary populations is essential to the success of elasmobranch conservation programmes, and to understanding their recent evolutionary history. For benthic elasmobranchs such as skates, traditional fisheries-independent approaches are often unsuitable as the data may be subject to various sources of bias, whilst low recapture rates can render mark-recapture programmes ineffectual. Close-kin mark-recapture (CKMR), a novel demographic modelling approach based on the genetic identification of close relatives within a sample, represents a promising alternative approach as it does not require physical recaptures. We evaluated the suitability of CKMR as a demographic modelling tool for the critically endangered blue skate (Dipturus batis) in the Celtic Sea using samples collected during fisheries-dependent trammel-net surveys that ran from 2011 to 2017. We identified three full-sibling and 16 half-sibling pairs among 662 skates, which were genotyped across 6291 genome-wide single nucleotide polymorphisms, 15 of which were cross-cohort half-sibling pairs that were included in a CKMR model. Despite limitations owing to a lack of validated life-history trait parameters for the species, we produced the first estimates of adult breeding abundance, population growth rate, and annual adult survival rate for D. batis in the Celtic Sea. The results were compared to estimates of genetic diversity, effective population size (N e ), and to catch per unit effort estimates from the trammel-net survey. Although each method was characterized by wide uncertainty bounds, together they suggested a stable population size across the time-series. Recommendations for the implementation of CKMR as a conservation tool for data-limited elasmobranchs are discussed. In addition, the spatio-temporal distribution of the 19 sibling pairs revealed a pattern of site fidelity in D. batis, and supported field observations suggesting an area of critical habitat that could qualify for protection might occur near the Isles of Scilly.

Feeding ecology of broadbill swordfish (Xiphias gladius) in the California current.

The feeding ecology of broadbill swordfish (Xiphias gladius) in the California Current was described based on analysis of stomach contents collected by fishery observers aboard commercial drift gillnet boats from 2007 to 2014. Prey were identified to the lowest taxonomic level and diet composition was analyzed using univariate and multivariate methods. Of 299 swordfish sampled (74 to 245 cm eye-to-fork length), 292 non-empty stomachs contained remains from 60 prey taxa. Genetic analyses were used to identify prey that could not be identified visually. Diet consisted mainly of cephalopods but also included epipelagic and mesopelagic teleosts. Jumbo squid (Dosidicus gigas) and Gonatopsis borealis were the most important prey based on the geometric index of importance. Swordfish diet varied with body size, location and year. Jumbo squid, Gonatus spp. and Pacific hake (Merluccius productus) were more important for larger swordfish, reflecting the ability of larger specimens to catch large prey. Jumbo squid, Gonatus spp. and market squid (Doryteuthis opalescens) were more important in inshore waters, while G. borealis and Pacific hake predominated offshore. Jumbo squid was more important in 2007-2010 than in 2011-2014, with Pacific hake being the most important prey item in the latter period. Diet variation by area and year probably reflects differences in swordfish preference, prey availability, prey distribution, and prey abundance. The range expansion of jumbo squid that occurred during the first decade of this century may particularly explain their prominence in swordfish diet during 2007-2010. Some factors (swordfish size, area, time period, sea surface temperature) that may influence dietary variation in swordfish were identified. Standardizing methods could make future studies more comparable for conservation monitoring purposes.

Characteristics of the spiny dogfish (Squalus acanthias) nuclear genome.

Sequenced shark nuclear genomes are underrepresented, with reference genomes available for only four out of nine orders so far. Here, we present the nuclear genome, with annotations, of the spiny dogfish (Squalus acanthias), a shark of interest to biomedical and conservation efforts, and the first representative of the second largest order of sharks (Squaliformes) with nuclear genome annotations available. Using Pacific Biosciences Continuous Long Read data in combination with Illumina paired-end and Hi-C sequencing, we assembled the genome de novo, followed by RNA-Seq-supported annotation. The final chromosome-level assembly is 3.7 Gb in size, has a BUSCO completeness score of 91.6%, and an error rate of less than 0.02%. Annotation predicted 33,283 gene models in the spiny dogfish's genome, of which 31,979 are functionally annotated.

The complete mitochondrial genome of the Atlantic spiny lumpsucker Eumicrotremus spinosus (Fabricius, 1776).

The complete mitogenome of the Atlantic spiny lumpsucker (Eumicrotremus spinosus) was generated using the PacBio Sequel II HiFi sequencing platform. The mitogenome assembly has a length of 19,281 bp and contains 13 protein-coding sequences, 22 tRNA genes, 2 rRNA genes, one control region containing the D-loop (2383 bp) and a duplicate control region (1133 bp) Phylogenetic analysis using maximum likelihood revealed that E. spinosus is closely related to the Siberian lumpsucker (E. asperrimus). The mitogenome of the spiny lumpsucker will be useful in population genomics and systematic studies of Cyclopteridae, Liparidae, and Cottidae.

Common aquatic pollutants modify hemocyte immune responses in Biomphalaria glabrata.

Disruptions to reproductive health in wildlife species inhabiting polluted environments is often found to occur alongside compromised immunity. However, research on impacts of aquatic pollution on freshwater mollusc immune responses is limited despite their importance as vectors of disease (Schistosomiasis) in humans, cattle and wild mammals. We developed an in vitro 'tool-kit' of well-characterized quantitative immune tests using Biomphalaria glabrata hemocytes. We exposed hemocytes to environmentally-relevant concentrations of common aquatic pollutants (17ß-estradiol, Bisphenol-A and p,p'-DDE) and measured key innate immune responses including motility, phagocytosis and encapsulation. Additionally, we tested an extract of a typical domestic tertiary treated effluent as representative of a 'real-world' mixture of chemicals. Encapsulation responses were stimulated by p,p'-DDE at low doses but were suppressed at higher doses. Concentrations of BPA (above 200 ng/L) and p,p'-DDE (above 500 ng/L) significantly inhibited phagocytosis compared to controls, whilst hemocyte motility was reduced by all test chemicals and the effluent extract in a dose-dependent manner. All responses occurred at chemical concentrations considered to be below the cytotoxic thresholds of hemocytes. This is the first time a suite of in vitro tests has been developed specifically in B. glabrata with the purpose of investigating the impacts of chemical pollutants and an effluent extract on immunity. Our findings indicate that common aquatic pollutants alter innate immune responses in B. glabrata, suggesting that pollutants may be a critical, yet overlooked, factor impacting disease by modulating the dynamics of parasite transmission between molluscs and humans.

The complete mitochondrial genomes of the flapper skate Dipturus intermedius and the longnose skate Dipturus oxyrinchus.

We describe the complete mitochondrial genomes of the flapper skate Dipturus intermedius (Parnell 1837) and the longnose skate Dipturus oxyrinchus (Linnaeus 1758), which have been obtained by Sanger sequencing. We report the length of the sequences to be 16,906 and 16,911 bp, respectively. The length and structure of gene regions, containing 13 protein-coding regions, 22 tRNA genes, two rRNA genes, and two non-coding areas, resemble those of related skate species. Despite D. intermedius being considered a cryptic species with D. batis, the full mitogenomes confirm that D. intermedius and D. oxyrinchus are more genetically similar. In comparison to other Dipturus species, D. intermedius is missing a whole codon in its cytochrome oxidase subunit 2 gene. These mitogenomes will be a useful resource furthering investigation of the population genetic differences and evolutionary history of skate species.

Population and seascape genomics of a critically endangered benthic elasmobranch, the blue skate Dipturus batis.

The blue skate (Dipturus batis) has a patchy distribution across the North-East Atlantic Ocean, largely restricted to occidental seas around the British Isles following fisheries-induced population declines and extirpations. The viability of remnant populations remains uncertain and could be impacted by continued fishing and by-catch pressure, and the projected impacts of climate change. We genotyped 503 samples of D. batis, obtained opportunistically from the widest available geographic range, across 6 350 single nucleotide polymorphisms (SNPs) using a reduced-representation sequencing approach. Genotypes were used to assess the species' contemporary population structure, estimate effective population sizes and identify putative signals of selection in relation to environmental variables using a seascape genomics approach. We identified genetic discontinuities between inshore (British Isles) and offshore (Rockall and Faroe Island) populations, with differentiation most pronounced across the deep waters of the Rockall Trough. Effective population sizes were largest in the Celtic Sea and Rockall, but low enough to be of potential conservation concern among Scottish and Faroese sites. Among the 21 candidate SNPs under positive selection was one significantly correlated with environmental variables predicted to be affected by climate change, including bottom temperature, salinity and pH. The paucity of well-annotated elasmobranch genomes precluded us from identifying a putative function for this SNP. Nevertheless, our findings suggest that climate change could inflict a strong selective force upon remnant populations of D. batis, further constraining its already-restricted habitat. Furthermore, the results provide fundamental insights on the distribution, behaviour and evolutionary biology of D. batis in the North-East Atlantic that will be useful for the establishment of conservation actions for this and other critically endangered elasmobranchs.

Antipodean white sharks on a Mediterranean walkabout? Historical dispersal leads to genetic discontinuity and an endangered anomalous population.

The provenance of white sharks (Carcharodon carcharias) in the Mediterranean is both a conundrum and an important conservation issue. Considering this species's propensity for natal philopatry, any evidence that the Mediterranean stock has little or no contemporary immigration from the Atlantic would suggest that it is extraordinarily vulnerable. To address this issue we sequenced the mitochondrial control region of four rare Mediterranean white sharks. Unexpectedly, the juvenile sequences were identical although collected at different locations and times, showing little genetic differentiation from Indo-Pacific lineages, but strong separation from geographically closer Atlantic/western Indian Ocean haplotypes. Historical long-distance dispersal (probably a consequence of navigational error during past climatic oscillations) and potential founder effects are invoked to explain the anomalous relationships of this isolated 'sink' population, highlighting the present vulnerability of its nursery grounds.

Genotypic and phenotypic correlates with proliferative kidney disease-induced mortality in Atlantic salmon.

Heritable variation in resistance to pathogens has been reported in many fish species, but little is known about its genetic architecture. To extend understanding, an investigation was made of the association of resistance to proliferative kidney disease (PKD) in 4 second filial generation (F2) families of Atlantic salmon with molecular markers from different genetic linkage groups in the species' genome, following a natural disease outbreak. PKD causes serious mortality in cultured salmonids. In addition to mortality, associations with growth-related traits were also examined, as immune responses are energetically costly and have been observed to reduce growth. Associations were investigated for 34 microsatellite markers and 5 restriction fragment length polymorphism (RFLP) loci from 3 regions of the growth hormone 1 gene (GH1). The phenotypic and genotypic character of survivors was compared with unexposed fish derived from the same families. Mortality was not size-selective, but growth in the survivors was reduced, and fish had a lower condition factor than unexposed fish, suggesting an energetic cost to resistance. Five markers showed significant allele frequency differences between survivors and unexposed fish, albeit in single families. Prior to correction for multiple tests, 2 of these markers were also linked to variation in growth-related traits among survivors, along with a further 7 markers. Though sample sizes constrained the power of the analysis, the study points to regions of the salmon genome that may contain quantitative trait loci related to PKD resistance, on which further work on the genetic architecture of PKD resistance in this species could focus.

The complete mitochondrial genome of the blue skate Dipturus batis.

The complete mitochondrial genome of the blue skate Dipturus batis is described from shotgun sequencing on an Illumina next-generation sequencing platform. We report a 16,911 bp long sequence similar in size to other members of the genus, containing 13 protein-coding regions, 22 tRNA genes, 2 rRNA genes, and 2 non-coding areas. Phylogenetic analysis was performed using the complete mitochondrial genomes of 17 related species, placing D. batis within the Rajini tribe of the Rajidae family, consistent with current taxonomy. The new resource adds to a growing database of rajid mitogenomes which will help resolve phylogenetic relationships within the family.

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.

Of mice and (Viking?) men: phylogeography of British and Irish house mice.

The west European subspecies of house mouse (Mus musculus domesticus) has gained much of its current widespread distribution through commensalism with humans. This means that the phylogeography of M. m. domesticus should reflect patterns of human movements. We studied restriction fragment length polymorphism (RFLP) and DNA sequence variations in mouse mitochondrial (mt) DNA throughout the British Isles (328 mice from 105 localities, including previously published data). There is a major mtDNA lineage revealed by both RFLP and sequence analyses, which is restricted to the northern and western peripheries of the British Isles, and also occurs in Norway. This distribution of the 'Orkney' lineage fits well with the sphere of influence of the Norwegian Vikings and was probably generated through inadvertent transport by them. To form viable populations, house mice would have required large human settlements such as the Norwegian Vikings founded. The other parts of the British Isles (essentially most of mainland Britain) are characterized by house mice with different mtDNA sequences, some of which are also found in Germany, and which probably reflect both Iron Age movements of people and mice and earlier development of large human settlements. MtDNA studies on house mice have the potential to reveal novel aspects of human history.

Distance decay of similarity among parasite communities of three marine invertebrate hosts.

The similarity in species composition between two communities generally decays as a function of increasing distance between them. Parasite communities in vertebrate definitive hosts follow this pattern but the respective relationship in intermediate invertebrate hosts of parasites with complex life cycles is unknown. In intermediate hosts, parasite communities are affected not only by the varying vagility of their definitive hosts (dispersing infective propagules) but also by the necessary coincidence of all their hosts in environmentally suitable localities. As intermediate hosts often hardly move they do not contribute to parasite dispersal. Hence, their parasite assemblages may decrease faster in similarity with increasing distance than those in highly mobile vertebrate definitive hosts. We use published field survey data to investigate distance decay of similarity in trematode communities from three prominent coastal molluscs of the Eastern North-Atlantic: the gastropods Littorina littorea and Hydrobia ulvae, and the bivalve Cerastoderma edule. We found that the similarity of trematode communities in all three hosts decayed with distance, independently of local sampling effort, and whether or not the parasites used the mollusc as first or second intermediate host in their life cycle. In H. ulvae, the halving distance (i.e. the distance that halves the similarity from its initial similarity at 1 km distance) for the trematode species using birds as definitive hosts was approximately two to three times larger than for species using fish. The initial similarities (estimated at 1 km distance) among trematode communities were relatively higher, whereas mean halving distances were lower, compared to published values for parasite communities in vertebrate hosts. We conclude that the vagility of definitive hosts accounts for a high similarity at the local scale, while the strong decay of similarity across regions is a consequence of the low probability that all necessary hosts and suitable environmental conditions coincide on a large scale.

Biomphalaria glabrata transcriptome: cDNA microarray profiling identifies resistant- and susceptible-specific gene expression in haemocytes from snail strains exposed to Schistosoma mansoni.

BACKGROUND: Biomphalaria glabrata is an intermediate snail host for Schistosoma mansoni, one of the important schistosomes infecting man. B. glabrata/S. mansoni provides a useful model system for investigating the intimate interactions between host and parasite. Examining differential gene expression between S. mansoni-exposed schistosome-resistant and susceptible snail lines will identify genes and pathways that may be involved in snail defences. RESULTS: We have developed a 2053 element cDNA microarray for B. glabrata containing clones from ORESTES (Open Reading frame ESTs) libraries, suppression subtractive hybridization (SSH) libraries and clones identified in previous expression studies. Snail haemocyte RNA, extracted from parasite-challenged resistant and susceptible snails, 2 to 24 h post-exposure to S. mansoni, was hybridized to the custom made cDNA microarray and 98 differentially expressed genes or gene clusters were identified, 94 resistant-associated and 4 susceptible-associated. Quantitative PCR analysis verified the cDNA microarray results for representative transcripts. Differentially expressed genes were annotated and clustered using gene ontology (GO) terminology and Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analysis. 61% of the identified differentially expressed genes have no known function including the 4 susceptible strain-specific transcripts. Resistant strain-specific expression of genes implicated in innate immunity of invertebrates was identified, including hydrolytic enzymes such as cathepsin L, a cysteine proteinase involved in lysis of phagocytosed particles; metabolic enzymes such as ornithine decarboxylase, the rate-limiting enzyme in the production of polyamines, important in inflammation and infection processes, as well as scavenging damaging free radicals produced during production of reactive oxygen species; stress response genes such as HSP70; proteins involved in signalling, such as importin 7 and copine 1, cytoplasmic intermediate filament (IF) protein and transcription enzymes such as elongation factor 1alpha and EF-2. CONCLUSION: Production of the first cDNA microarray for profiling gene expression in B. glabrata provides a foundation for expanding our understanding of pathways and genes involved in the snail internal defence system (IDS). We demonstrate resistant strain-specific expression of genes potentially associated with the snail IDS, ranging from signalling and inflammation responses through to lysis of proteinacous products (encapsulated sporocysts or phagocytosed parasite components) and processing/degradation of these targeted products by ubiquitination.

Identification of genes involved in interactions between Biomphalaria glabrata and Schistosoma mansoni by suppression subtractive hybridization.

Biomphalaria glabrata is an intermediate snail host for Schistosoma mansoni, a medically important schistosome. In order to identify transcripts involved in snail-schistosome interactions, subtractive cDNA libraries were prepared, using suppression subtractive hybridization (SSH) between a parasite-exposed schistosome-resistant and a susceptible strain of B. glabrata, and also between schistosome-exposed and unexposed snails from the resistant snail line. Separate libraries were made from both haemocytes and the haemopoietic organ. Subtraction was performed in both directions enriching for cDNAs differentially expressed between parasite-exposed resistant and susceptible samples and up or down-regulated in the resistant line after challenge. The resulting eight libraries were screened and eight genes, differentially expressed between the haemocytes of resistant and susceptible snail strains, were identified and confirmed with reverse transcriptase PCR, including two transcripts expected to be involved in the stress response mechanism for regulating the damaging oxidative burst pathways involved in cytotoxic killing of the parasite: the iron-storage and immunoregulatory molecule, ferritin, and HtrA2, a serine protease involved in the cellular stress response. Transcripts with elevated levels in the resistant strain, had the same expression patterns in the subtracted libraries and unsubtracted controls; higher levels in exposed resistant snails compared to susceptible ones and down-regulated in exposed compared with unexposed resistant snails. Differential expression of two of the transcripts with no known function from the susceptible strain, was independently confirmed in a repeat exposure experiment.

Biomphalaria glabrata transcriptome: identification of cell-signalling, transcriptional control and immune-related genes from open reading frame expressed sequence tags (ORESTES).

Biomphalaria glabrata is the major intermediate snail host for Schistosoma mansoni, one of the important schistosomes infecting man. Much remains to be discovered concerning specific molecules mediating the defence events in these intermediate hosts, triggered by invading schistosomes. An expressed sequence tag (EST) gene discovery strategy known as ORESTES has been employed to identify transcripts that might be involved in snail-schistosome interactions in order to examine gene expression patterns in infected B. glabrata. Over 3930 ESTs were sequenced from cDNA libraries made from both schistosome-exposed and unexposed snails using different tissue types, producing a database of 1843 non-redundant clones. The non-redundant set has been assessed for gene ontology and KEGG pathway assignments. This approach has revealed a number of signalling, antioxidant and immune-related gene homologues that, based on current understanding of molluscan and other comparative systems, might play an important role in the molluscan defence response towards infection.

World without borders-genetic population structure of a highly migratory marine predator, the blue shark (Prionace glauca).

Highly migratory, cosmopolitan oceanic sharks often exhibit complex movement patterns influenced by ontogeny, reproduction, and feeding. These elusive species are particularly challenging to population genetic studies, as representative samples suitable for inferring genetic structure are difficult to obtain. Our study provides insights into the genetic population structure one of the most abundant and wide-ranging oceanic shark species, the blue shark Prionace glauca, by sampling the least mobile component of the populations, i.e., young-of-year and small juveniles (<2 year; N = 348 individuals), at three reported nursery areas, namely, western Iberia, Azores, and South Africa. Samples were collected in two different time periods (2002-2008 and 2012-2015) and were screened at 12 nuclear microsatellites and at a 899-bp fragment of the mitochondrial control region. Our results show temporally stable genetic homogeneity among the three Atlantic nurseries at both nuclear and mitochondrial markers, suggesting basin-wide panmixia. In addition, comparison of mtDNA CR sequences from Atlantic and Indo-Pacific locations also indicated genetic homogeneity and unrestricted female-mediated gene flow between ocean basins. These results are discussed in light of the species' life history and ecology, but suggest that blue shark populations may be connected by gene flow at the global scale. The implications of the present findings to the management of this important fisheries resource are also discussed.