Phylogeography and conservation of the bull shark (Carcharhinus leucas) inferred from mitochondrial and microsatellite DNA

Seasonal occurrence and individual variability of bull sharks, Carcharhinus leucas, in a marine reserve of the southwestern Gulf of California

Background

Studying how the bull sharks aggregate and how they can be driven by life history traits such as reproduction, prey availability, predator avoidance and social interaction in a National Park such as Cabo Pulmo, is key to understand and protect the species.

Methods

The occurrence variability of 32 bull sharks tracked with passive acoustic telemetry were investigated via a hierarchical logistic regression model, with inference conducted in a Bayesian framework, comparing sex, and their response to temperature and chlorophyll.

Results

Based on the fitted model, occurrence probability varied by sex and length. Juvenile females had the highest values, whereas adult males the lowest. A strong seasonality or day of the year was recorded, where sharks were generally absent during September-November. However, some sharks did not show the common pattern, being detected just for a short period. This is one of the first studies where the Bayesian framework is used to study passive acoustic telemetry proving the potential to be used in further studies.

Origins of juvenile green sea turtles (Chelonia mydas) in the Bahamas: A comparison of recent and historical rookery contributions

Conservation of green sea turtles (Chelonia mydas) benefits from knowledge of population connectivity across life stages. Green turtles are managed at the level of genetically discrete rookeries, yet individuals from different rookeries mix at foraging grounds; therefore, rookeries may be impacted by processes at foraging grounds. Bimini, Bahamas, hosts an important foraging assemblage, but rookery contributions to this assemblage have never been resolved. We generated mitochondrial DNA sequences for 96 foraging green turtles from Bimini and used Mixed Stock Analysis to determine rookery contributions to this population using 817 and 490 base pair (bp) rookery baseline data. The MSA conducted with 817 bp data indicated that Quintana Roo, Mexico, and Central Eastern Florida contributed most to the Bimini population. The MSA conducted with 490 bp data indicated that Southwest Cuba and Central Eastern Florida contributed the most to Bimini. The results of the second MSA differ from a previous study undertaken with 490 bp data, conducted in Great Inagua, Bahamas, which suggested that Tortuguero, Costa Rica, contributed the most to that foraging assemblage. Large credible intervals in our results do not permit explicit interpretation of individual rookery contributions, but our results do indicate substantial relative differences in rookery contributions to two Bahamian foraging assemblages which may be driven by oceanic currents, rookery sizes, and possibly juvenile natal homing. Our findings may implicate a shift in contributions to the Bahamas over two decades, highlighting the importance of regularly monitoring rookery contributions and resolving regional recruitment patterns to inform conservation.

High-order correlations in species interactions lead to complex diversity-stability relationships for ecosystems

How ecosystems maintain stability is an active area of research. Inspired by applications of random matrix theory in nuclear physics, May showed decades ago that in an ecosystem model with many randomly interacting species, increasing species diversity decreases the stability of the ecosystem. There have since been many additions to May's efforts, one being an improved understanding the effect of mutualistic, competitive, or predator-prey-like correlations between pairs of species. Here we extend a random matrix technique developed in the context of spin-glass theory to study the effect of high-order correlations among species interactions. The resulting analytically solvable models include next-to-nearest-neighbor correlations in the species interaction network, such as the enemy of my enemy is my friend, as well as higher-order correlations. We find qualitative differences from May and others' models, including nonmonotonic diversity-stability relationships. Furthermore, inclusion of particular next-to-nearest-neighbor correlations in predator-prey as opposed to mutualist-competitive networks causes the former to transition to being more stable at higher species diversity. We discuss potential applicability of our results to microbiota engineering and to the ecology of interpredator interactions, such as cub predation between lions and hyenas as well as companionship between humans and dogs.

Mitochondrial DNA genome evidence for the existence of a third divergent lineage in the western Atlantic Ocean for the bull shark (Carcharhinus leucas)

We report for the first time a highly divergent lineage in the Caribbean Sea for the bull shark (Carcharhinus leucas) based on the analysis of 51 mitochondrial DNA genomes of individuals collected in the western North Atlantic. When comparing the mtDNA control region obtained from the mitogenomes to sequences reported previously for Brazil, the Caribbean lineage remained highly divergent. These results support the existence of a discrete population in Central America due to a phylogeographic break separating the Caribbean Sea from the western North Atlantic, Gulf of Mexico and South America.

Development and characterization of thirteen novel microsatellite markers for use in Greenland sharks (Somniosus microcephalus), with cross-amplification in Pacific sleeper sharks (Somniosus pacificus)

Objective

The objectives of this work are to isolate, develop, and characterize polymorphic microsatellite markers for use in Greenland sharks (Somniosus microcephalus).

Results

Thirteen microsatellite loci were successfully amplified and yielded multi-locus genotypes for 36 S. microcephalus individuals from Grise Fjord (n = 16) and Svalbard (n = 20). Each locus yielded between 2 and 9 alleles and observed heterozygosity ranged from 0.11 to 0.70 when estimated across both sites. One locus and three loci deviated from HWE following Bonferroni correction, for individuals sampled from Grise Fjord and Svalbard, respectively. Cross-amplification was successful at every locus for five of the ten S. pacificus individuals.

Genomics versus mtDNA for resolving stock structure in the silky shark (Carcharhinus falciformis)

Conservation genetic approaches for elasmobranchs have focused on regions of the mitochondrial genome or a handful of nuclear microsatellites. High-throughput sequencing offers a powerful alternative for examining population structure using many loci distributed across the nuclear and mitochondrial genomes. These single nucleotide polymorphisms are expected to provide finer scale and more accurate population level data; however, there have been few genomic studies applied to elasmobranch species. The desire to apply next-generation sequencing approaches is often tempered by the costs, which can be offset by pooling specimens prior to sequencing (pool-seq). In this study, we assess the utility of pool-seq by applying this method to the same individual silky sharks, Carcharhinus falciformis, previously surveyed with the mtDNA control region in the Atlantic and Indian Oceans. Pool-seq methods were able to recover the entire mitochondrial genome as well as thousands of nuclear markers. This volume of sequence data enabled the detection of population structure between regions of the Atlantic Ocean populations, undetected in the previous study (inter-Atlantic mitochondrial SNPs FST values comparison ranging from 0.029 to 0.135 and nuclear SNPs from 0.015 to 0.025). Our results reinforce the conclusion that sampling the mitochondrial control region alone may fail to detect fine-scale population structure, and additional sampling across the genome may increase resolution for some species. Additionally, this study shows that the costs of analyzing 4,988 loci using pool-seq methods are equivalent to the standard Sanger-sequenced markers and become less expensive when large numbers of individuals (>300) are analyzed.

Identifying mechanisms of genetic differentiation among populations in vagile species: historical factors dominate genetic differentiation in seabirds

Elucidating the factors underlying the origin and maintenance of genetic variation among populations is crucial for our understanding of their ecology and evolution, and also to help identify conservation priorities. While intrinsic movement has been hypothesized as the major determinant of population genetic structuring in abundant vagile species, growing evidence indicates that vagility does not always predict genetic differentiation. However, identifying the determinants of genetic structuring can be challenging, and these are largely unknown for most vagile species. Although, in principle, levels of gene flow can be inferred from neutral allele frequency divergence among populations, underlying assumptions may be unrealistic. Moreover, molecular studies have suggested that contemporary gene flow has often not overridden historical influences on population genetic structure, which indicates potential inadequacies of any interpretations that fail to consider the influence of history in shaping that structure. This exhaustive review of the theoretical and empirical literature investigates the determinants of population genetic differentiation using seabirds as a model system for vagile taxa. Seabirds provide a tractable group within which to identify the determinants of genetic differentiation, given their widespread distribution in marine habitats and an abundance of ecological and genetic studies conducted on this group. Herein we evaluate mitochondrial DNA (mtDNA) variation in 73 seabird species. Lack of mutation-drift equilibrium observed in 19% of species coincided with lower estimates of genetic differentiation, suggesting that dynamic demographic histories can often lead to erroneous interpretations of contemporary gene flow, even in vagile species. Presence of land across the species sampling range, or sampling of breeding colonies representing ice-free Pleistocene refuge zones, appear to be associated with genetic differentiation in Tropical and Southern Temperate species, respectively, indicating that long-term barriers and persistence of populations are important for their genetic structuring. Conversely, biotic factors commonly considered to influence population genetic structure, such as spatial segregation during foraging, were inconsistently associated with population genetic differentiation. In light of these results, we recommend that genetic studies should consider potential historical events when identifying determinants of genetic differentiation among populations to avoid overestimating the role of contemporary factors, even for highly vagile taxa.

Population-level inferences from environmental DNA-Current status and future perspectives

Environmental DNA (eDNA) extracted from water samples has recently shown potential as a valuable source of population genetic information for aquatic macroorganisms. This approach offers several potential advantages compared with conventional tissue-based methods, including the fact that eDNA sampling is noninvasive and generally more cost-efficient. Currently, eDNA approaches have been limited to single-marker studies of mitochondrial DNA (mtDNA), and the relationship between eDNA haplotype composition and true haplotype composition still needs to be thoroughly verified. This will require testing of bioinformatic and statistical software to correct for erroneous sequences, as well as biases and random variation in relative sequence abundances. However, eDNA-based population genetic methods have far-reaching potential for both basic and applied research. In this paper, we present a brief overview of the achievements of eDNA-based population genetics to date, and outline the prospects for future developments in the field, including the estimation of nuclear DNA (nuDNA) variation and epigenetic information. We discuss the challenges associated with eDNA samples as opposed to those of individual tissue samples and assess whether eDNA might offer additional types of information unobtainable with tissue samples. Lastly, we provide recommendations for determining whether an eDNA approach would be a useful and suitable choice in different research settings. We limit our discussion largely to contemporary aquatic systems, but the advantages, challenges, and perspectives can to a large degree be generalized to eDNA studies with a different spatial and temporal focus.

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.

Reproductive biology, multiple paternity and polyandry of the bull shark Carcharhinus leucas

To improve understanding of bull shark Carcharhinus leucas reproductive biology, we analysed reproductive traits from 118 bull sharks caught along Reunion Island coasts (Western Indian Ocean), including 16 gravid females. Specific microsatellite loci were used to investigate the frequency of multiple paternity. Males and females reached maturity at c. 234 cm and 257 cm total length (L_T ), respectively, and litter sizes ranged from 5 to 14 embryos. Analysis of the 16 litters collected in various months of the year indicated that parturition occurs between October and December, with a size at birth c. 60-80 cm L_T and that the gestation period is probably c. 12 months. Assuming a 1 year resting period and a period of sperm storage (4-5 months) between mating (in June-September) and fertilisation, the reproductive cycle of bull sharks at Reunion Island would be biennial. At least 56.25% of the litters investigated were polyandrous, sired by 2-5 males. Several males that each sired several litters conceived during the same or distinct mating seasons were detected, suggesting both a seasonal aggregation of sharks to mate and some male fidelity to mating site. Altogether, these findings provide valuable information for both shark risk management and conservation of the species in the Western Indian Ocean.

Genetic population structure and demography of an apex predator, the tiger shark Galeocerdo cuvier

Population genetics has been increasingly applied to study large sharks over the last decade. Whilst large shark species are often difficult to study with direct methods, improved knowledge is needed for both population management and conservation, especially for species vulnerable to anthropogenic and climatic impacts. The tiger shark, Galeocerdo cuvier, is an apex predator known to play important direct and indirect roles in tropical and subtropical marine ecosystems. While the global and Indo-West Pacific population genetic structure of this species has recently been investigated, questions remain over population structure and demographic history within the western Indian (WIO) and within the western Pacific Oceans (WPO). To address the knowledge gap in tiger shark regional population structures, the genetic diversity of 286 individuals sampled in seven localities was investigated using 27 microsatellite loci and three mitochondrial genes (CR,COI, and cytb). A weak genetic differentiation was observed between the WIO and the WPO, suggesting high genetic connectivity. This result agrees with previous studies and highlights the importance of the pelagic behavior of this species to ensure gene flow. Using approximate Bayesian computation to couple information from both nuclear and mitochondrial markers, evidence of a recent bottleneck in the Holocene (2,000-3,000 years ago) was found, which is the most probable cause for the low genetic diversity observed. A contemporary effective population size as low as 111 [43,369] was estimated during the bottleneck. Together, these results indicate low genetic diversity that may reflect a vulnerable population sensitive to regional pressures. Conservation measures are thus needed to protect a species that is classified as Near Threatened.

Taking Advantage of the Genomics Revolution for Monitoring and Conservation of Chondrichthyan Populations

Chondrichthyes (sharks, rays, skates and chimaeras) are among the oldest extant predators and are vital to top-down regulation of oceanic ecosystems. They are an ecologically diverse group occupying a wide range of habitats and are thus, exploited by coastal, pelagic and deep-water fishing industries. Chondrichthyes are among the most data deficient vertebrate species groups making design and implementation of regulatory and conservation measures challenging. High-throughput sequencing technologies have significantly propelled ecological investigations and understanding of marine and terrestrial species’ populations, but there remains a paucity of NGS based research on chondrichthyan populations. We present a brief review of current methods to access genomic and metagenomic data from Chondrichthyes and discuss applications of these datasets to increase our understanding of chondrichthyan taxonomy, evolution, ecology and population structures. Last, we consider opportunities and challenges offered by genomic studies for conservation and management of chondrichthyan populations.

Trends in global shark attacks

Shark attacks are a global phenomenon that attracts widespread attention and publicity, often with negative outcomes for shark populations. Despite the widespread perceptions of shark attacks, trends in human water activities and shark populations are both dynamic, resulting in variable rates of shark attacks over space and time. Understanding variable trends in shark attacks may contribute to a better understanding of risk, and a more tempered response in the wake of an attack. We found that global shark attack rates are low, yet variable across global regions and over decades. Countries with low populations were found to have the highest rates of attack, while countries with high populations (U.S.A., Australia, South Africa) tended to have overall low attack rates, but also much more interannual variability. From the 1960s to the present, those countries with the highest populations also tended to be the places where attack rates have increased. Ultimately, shark attack risk is also driven by local conditions (e.g., time of day, species present); however, a global scale understanding of attack rates helps place risk into perspective and may contribute to a more scientifically-grounded discussion of sharks, and their management and conservation.

Strong trans-Pacific break and local conservation units in the Galapagos shark (Carcharhinus galapagensis) revealed by genome-wide cytonuclear markers

The application of genome-wide cytonuclear molecular data to identify management and adaptive units at various spatio-temporal levels is particularly important for overharvested large predatory organisms, often characterized by smaller, localized populations. Despite being "near threatened", current understanding of habitat use and population structure of Carcharhinus galapagensis is limited to specific areas within its distribution. We evaluated population structure and connectivity across the Pacific Ocean using genome-wide single-nucleotide polymorphisms (~7200 SNPs) and mitochondrial control region sequences (945 bp) for 229 individuals. Neutral SNPs defined at least two genetically discrete geographic groups: an East Tropical Pacific (Mexico, east and west Galapagos Islands), and another central-west Pacific (Lord Howe Island, Middleton Reef, Norfolk Island, Elizabeth Reef, Kermadec, Hawaii and Southern Africa). More fine-grade population structure was suggested using outlier SNPs: west Pacific, Hawaii, Mexico, and Galapagos. Consistently, mtDNA pairwise ΦST defined three regional stocks: east, central and west Pacific. Compared to neutral SNPs (FST = 0.023-0.035), mtDNA exhibited more divergence (ΦST = 0.258-0.539) and high overall genetic diversity (h = 0.794 ± 0.014; π = 0.004 ± 0.000), consistent with the longstanding eastern Pacific barrier between the east and central-west Pacific. Hawaiian and Southern African populations group within the west Pacific cluster. Effective population sizes were moderate/high for east/west populations (738 and 3421, respectively). Insights into the biology, connectivity, genetic diversity, and population demographics informs for improved conservation of this species, by delineating three to four conservation units across their Pacific distribution. Implementing such conservation management may be challenging, but is necessary to achieve long-term population resilience at basin and regional scales.

Fisheries-independent surveys identify critical habitats for young scalloped hammerhead sharks (Sphyrna lewini) in the Rewa Delta, Fiji

Sharp declines in numerous shark populations around the world have generated considerable interest in better understanding and characterising their biology, ecology and critical habitats. The scalloped hammerhead shark (SHS, Sphyrna lewini) is subject to a multitude of natural and anthropogenic threats that are often exacerbated within the coastal embayments and estuaries used during SHS early life stages. In this study, we describe the temporal and spatial distribution, age class composition, and reproductive biology of SHS in the Rewa Delta (RD), Fiji. A total of 1054 SHS (including 796 tagged individuals; 101 of which were recaptured) were captured from September 2014 to March 2016 in the RD. A majority of the captures in this area were neonates and young-of-the-year (YOY) (99.8%). Significant seasonality in patterns of occurrence of both neonates and YOY individuals suggests a defined parturition period during the austral summer. Between the seven sampling sites in the RD we also found significant differences in SHS neonate catch per unit of effort, and average total length of individuals. According to the data, the RD is likely to represent an important nursery area for SHS up to one year of age.

Bioecology and movements of bull sharks, Carcharhinus leucas , caught in a long-term longline survey off northeastern Brazil

ABSTRACT A robust understanding of habitat usage by coastal shark species, and how it overlaps with human presence in densely-populated regions is needed to inform the development of efficient conservation strategies for these important top predators. An intensive longline survey conducted in nearshore waters off northeastern Brazil from 2004 through 2014 caught a total of 18 bull sharks (Carcharhinus leucas) (male-female ratio = 0.63:1), which can be dangerous to humans. Although most sharks were sexually mature, there was no evidence that this region could be used as a parturition or nursery area. Prey items identified in the guts of the sharks comprised teleosts, mollusks and elasmobranchs. Additionally, one satellite-tagged bull shark covered a great distance (> 3,000 km) in 75 days at liberty, making most use of shallow waters (< 20 m depth) and presumably also entering an estuarine area. Although bull sharks are not an important fishery resource in this region, such a reduced abundance coupled with its affinity for coastal and inshore habitats highlights the potential vulnerability of C. leucas to deleterious anthropic interferences off northeastern Brazil.

Population structure and cryptic speciation in bonnethead sharks Sphyrna tiburo in the south-eastern U.S.A. and Caribbean

Population structure and lineage diversification within a small, non-dispersive hammerhead shark species, the bonnethead shark Sphyrna tiburo, was assessed. Sphyrna tiburo is currently described as one continuously distributed species along the Atlantic continental margins of North, Central and South America, but recent genetic analysis of an insular population (Trinidad) suggests the possibility of cryptic speciation. To address this issue S. tiburo were sampled at six sites along c. 6200 km of continuous, continental coastline and from one island location (Grand Bahama) across a discontinuity in their distribution (the Straits of Florida), in order to test if they constitute a single lineage over this distribution. A total of 1030 bp of the mitochondrial control region (CR) was obtained for 239 S. tiburo, revealing 73 distinct haplotypes, high nucleotide diversity (0·01089) and a pair of highly divergent lineages estimated to have separated 3·61-5·62 million years ago. Mitochondrial cytochrome oxidase I and nuclear internal transcribed spacer loci show the same pattern. Divergence is similar within S. tiburo to that observed between established elasmobranch sister species, providing further evidence of cryptic speciation. A global AMOVA based on CR confirms that genetic diversity is primarily partitioned among populations (Φ_ST = 0·828, P < 0·001) because the divergent lineages are almost perfectly segregated between Belize and North America-The Bahamas. An AMOVA consisting solely of the North American and Bahamian samples is also significantly different from zero (Φ_ST = 0·088, P < 0·001) and pairwise F_ST is significantly different between all sites. These findings suggest that S. tiburo comprises a species complex and supports previous research indicating fine population structure, which has implications for fisheries management and biodiversity conservation.

Structure and Genetic Variability of the Oceanic Whitetip Shark, Carcharhinus longimanus, Determined Using Mitochondrial DNA

Information regarding population structure and genetic connectivity is an important contribution when establishing conservation strategies to manage threatened species. The oceanic whitetip shark, Carcharhinus longimanus, is a highly migratory, large-bodied, pelagic shark listed by the IUCN (International Union for Conservation of Nature) Red List as "vulnerable" throughout its range and "critically endangered" in the western north Atlantic. In 2014, the species was protected globally under Appendix II of CITES (Convention on International Trade in Endangered Species), limiting and regulating trade. This study used partial sequences of mitochondrial DNA (mtDNA) control region to determine the population genetic structure of oceanic whitetip sharks across the Atlantic and Indian Oceans. 724 base pairs were obtained from 215 individuals that identifed nine polymorphic sites and defined 12 distinct haplotypes. Total nucleotide diversity (π) was 0.0013 and haplotype diversity (h) was 0.5953. The Analysis of Molecular Variance (AMOVA) evidenced moderate levels of population structure (ɸST = 0.1039) with restricted gene flow between the western and eastern Atlantic Ocean, and a strong relationship between the latter region and the Indian Ocean. Even though the oceanic whitetip is a highly migratory animal the results presented here show that their genetic variability is slightly below average of other pelagic sharks. Additionally, this study recommends that at least two populations in the Atlantic Ocean should be considered distinct (eastern and western Atlantic) and conservation efforts should be focused in areas with the greatest genetic diversity by environmental managers.

Long-term demographic decline and late glacial divergence in a Californian paleoendemic: Sequoiadendron giganteum (giant sequoia)

Mediterranean ecosystems comprise a high proportion of endemic taxa whose response to climate change will depend on their evolutionary origins. In the California flora, relatively little attention has been given to the evolutionary history of paleoendemics from a molecular perspective, yet they number among some of the world's most iconic plant species. Here, we address questions of demographic change in Sequoiadendron giganteum (giant sequoia) that is restricted to a narrow belt of groves in the Sierra Nevada Mountains. We ask whether the current distribution is a result of northward colonization since the last glacial maximum (LGM), restriction of a broader range in the recent past (LGM) or independent colonizations in the deeper past. Genetic diversity at eleven microsatellite loci decreased with increasing latitude, but partial regressions suggested this was a function of smaller population sizes in the north. Disjunct populations north of the Kings River were divergent from those south of the Kings River that formed a single cluster in Bayesian assignment tests. Demographic inferences supported a demographic contraction just prior to the LGM as the most likely scenario for the current disjunct range of the species. This contraction appeared to be superimposed upon a long‐term decline in giant sequoia over the last 2 million years, associated with increasing aridity due to the Mediterranean climate. Overall, low genetic diversity, together with competition in an environment to which giant sequoia is likely already poorly adapted, will pose major constraints on its success in the face of increasing aridity.

The Pillars of Hercules as a bathymetric barrier to gene flow promoting isolation in a global deep-sea shark (Centroscymnus coelolepis)

Knowledge of the mechanisms limiting connectivity and gene flow in deep-sea ecosystems is scarce, especially for deep-sea sharks. The Portuguese dogfish (Centroscymnus coelolepis) is a globally distributed and near threatened deep-sea shark. C. coelolepis population structure was studied using 11 nuclear microsatellite markers and a 497-bp fragment from the mtDNA control region. High levels of genetic homogeneity across the Atlantic (Φ(ST) = -0.0091, F(ST) = 0.0024, P > 0.05) were found suggesting one large population unit at this basin. The low levels of genetic divergence between Atlantic and Australia (Φ(ST) = 0.0744, P < 0.01; F(ST) = 0.0015, P > 0.05) further suggested that this species may be able to maintain some degree of genetic connectivity even across ocean basins. In contrast, sharks from the Mediterranean Sea exhibited marked genetic differentiation from all other localities studied (Φ(ST) = 0.3808, F(ST) = 0.1149, P < 0.001). This finding suggests that the shallow depth of the Strait of Gibraltar acts as a barrier to dispersal and that isolation and genetic drift may have had an important role shaping the Mediterranean shark population over time. Analyses of life history traits allowed the direct comparison among regions providing a complete characterization of this shark's populations. Sharks from the Mediterranean had markedly smaller adult body size and size at maturity compared to Atlantic and Pacific individuals. Together, these results suggest the existence of an isolated and unique population of C. coelolepis inhabiting the Mediterranean that most likely became separated from the Atlantic in the late Pleistocene.

To Madagascar and back: long-distance, return migration across open ocean by a pregnant female bull shark Carcharhinus leucas

A large, pregnant, female bull shark Carcharhinus leucas was tracked migrating from Seychelles across open ocean to south-east Madagascar, c. 2000 km away, and back again. In Madagascar, the shark spent a prolonged period shallower than 5 m, consistent with entering estuarine habitat to pup, and upon return to Seychelles the shark was slender and no longer gravid. This represents an unprecedented return migration across the open ocean for a C. leucas and highlights the need for international collaboration to manage the regional C. leucas population sustainably.

Microsatellite loci confirm a lack of population connectivity among globally distributed populations of the tope shark Galeorhinus galeus (Triakidae)

This study used 11 polymorphic nuclear microsatellite loci to determine the population connectivity of five geographically isolated populations of tope shark Galeorhinus galeus (Africa, Australia, North America, South America and western Europe). Genetic analyses revealed significant structure among all populations indicating a lack of population connectivity. These findings indicate that globally distributed populations of G. galeus are isolated and should be managed as distinct, independent stocks.

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.

There and back again: a review of residency and return migrations in sharks, with implications for population structure and management

The overexploitation of sharks has become a global environmental issue in need of a comprehensive and multifaceted management response. Tracking studies are beginning to elucidate how shark movements shape the internal dynamics and structure of populations, which determine the most appropriate scale of these management efforts. Tracked sharks frequently either remain in a restricted geographic area for an extended period of time (residency) or return to a previously resided-in area after making long-distance movements (site fidelity). Genetic studies have shown that some individuals of certain species preferentially return to their exact birthplaces (natal philopatry) or birth regions (regional philopatry) for either parturition or mating, even though they make long-distance movements that would allow them to breed elsewhere. More than 80 peer-reviewed articles, constituting the majority of published shark tracking and population genetic studies, provide evidence of at least one of these behaviors in a combined 31 shark species from six of the eight extant orders. Residency, site fidelity, and philopatry can alone or in combination structure many coastal shark populations on finer geographic scales than expected based on their potential for dispersal. This information should therefore be used to scale and inform assessment, management, and conservation activities intended to restore depleted shark populations.

Contemporary population structure and post-glacial genetic demography in a migratory marine species, the blacknose shark, Carcharhinus acronotus

Patterns of population structure and historical genetic demography of blacknose sharks in the western North Atlantic Ocean were assessed using variation in nuclear-encoded microsatellites and sequences of mitochondrial (mt)DNA. Significant heterogeneity and/or inferred barriers to gene flow, based on microsatellites and/or mtDNA, revealed the occurrence of five genetic populations localized to five geographic regions: the southeastern U.S Atlantic coast, the eastern Gulf of Mexico, the western Gulf of Mexico, Bay of Campeche in the southern Gulf of Mexico and the Bahamas. Pairwise estimates of genetic divergence between sharks in the Bahamas and those in all other localities were more than an order of magnitude higher than between pairwise comparisons involving the other localities. Demographic modelling indicated that sharks in all five regions diverged after the last glacial maximum and, except for the Bahamas, experienced post-glacial, population expansion. The patterns of genetic variation also suggest that the southern Gulf of Mexico may have served as a glacial refuge and source for the expansion. Results of the study demonstrate that barriers to gene flow and historical genetic demography contributed to contemporary patterns of population structure in a coastal migratory species living in an otherwise continuous marine habitat. The results also indicate that for many marine species, failure to properly characterize barriers in terms of levels of contemporary gene flow could in part be due to inferences based solely on equilibrium assumptions. This could lead to erroneous conclusions regarding levels of connectivity in species of conservation concern.

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.

Genetic structure of Barbus spp. populations in the Marches Region of central Italy and its relevance to conservation actions

A genetic survey of Barbus spp. populations in the Marches Region (Adriatic River basins), central Italy, was carried out using mitochondrial and nuclear markers (partial D-loop, cyt b sequences and microsatellite loci) in order to ascertain their systematic position and to address their genetic structure which is key to conservation action planning. Analyses were conducted on sequences obtained from 91 individuals collected from eight sampling sites in five different rivers, from two specimens provided by the Ichthyological Centre of Rome and mitochondrial sequences of Barbus spp. retrieved from GenBank. Presumptive classification based on external morphological characters was not confirmed by genetic analysis, by means of which all specimens collected in the Marches Region were ascribed to Barbus plebejus. Genetic diversity values (h and π) of sampling groups were all different from 0 except the one sample collected from the upper reaches of the River Tenna, above a hydroelectric dam. Population connectivity and colonization patterns of the studied area were inferred from an analysis of molecular variance distribution and evolutionary relationships among haplotypes. The results point to different levels of isolation among sampling groups due to ecological and anthropogenic factors and the effect of an artificial barrier on genetic variability and conservation status of the population. Finally, this study confirms the uncertainty associated with systematic classification of Barbus spp. based on morphological characters due to the phenotypic plasticity of the species.

Natural or artificial? Habitat-use by the bull shark, Carcharhinus leucas

Background

Despite accelerated global population declines due to targeted and illegal fishing pressure for many top-level shark species, the impacts of coastal habitat modification have been largely overlooked. We present the first direct comparison of the use of natural versus artificial habitats for the bull shark, Carcharhinus leucas, an IUCN ‘Near-threatened’ species - one of the few truly euryhaline sharks that utilises natural rivers and estuaries as nursery grounds before migrating offshore as adults. Understanding the value of alternate artificial coastal habitats to the lifecycle of the bull shark is crucial for determining the impact of coastal development on this threatened but potentially dangerous species.

Methodology/Findings

We used longline surveys and long-term passive acoustic tracking of neonate and juvenile bull sharks to determine the ontogenetic value of natural and artificial habitats to bull sharks associated with the Nerang River and adjoining canals on the Gold Coast, Australia. Long-term movements of tagged sharks suggested a preference for the natural river over artificial habitat (canals). Neonates and juveniles spent the majority of their time in the upper tidal reaches of the Nerang River and undertook excursions into adjoining canals. Larger bull sharks ranged further and frequented the canals closer to the river mouth.

Conclusions/Significance

Our work suggests with increased destruction of natural habitats, artificial coastal habitat may become increasingly important to large juvenile bull sharks with associated risk of attack on humans. In this system, neonate and juvenile bull sharks utilised the natural and artificial habitats, but the latter was not the preferred habitat of neonates. The upper reaches of tidal rivers, often under significant modification pressure, serve as nursery sites for neonates. Analogous studies are needed in similar systems elsewhere to assess the spatial and temporal generality of this research.

Evidence for reproductive philopatry in the bull shark Carcharhinus leucas

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

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

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

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

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

Global phylogeography with mixed-marker analysis reveals male-mediated dispersal in the endangered scalloped hammerhead shark (Sphyrna lewini)

Background

The scalloped hammerhead shark, Sphyrna lewini, is a large endangered predator with a circumglobal distribution, observed in the open ocean but linked ontogenetically to coastal embayments for parturition and juvenile development. A previous survey of maternal (mtDNA) markers demonstrated strong genetic partitioning overall (global Φ_ST = 0.749) and significant population separations across oceans and between discontinuous continental coastlines.

Methodology/Principal Findings

We surveyed the same global range with increased sample coverage (N = 403) and 13 microsatellite loci to assess the male contribution to dispersal and population structure. Biparentally inherited microsatellites reveal low or absent genetic structure across ocean basins and global genetic differentiation (F_ST = 0.035) over an order of magnitude lower than the corresponding measures for maternal mtDNA lineages (Φ_ST = 0.749). Nuclear allelic richness and heterozygosity are high throughout the Indo-Pacific, while genetic structure is low. In contrast, allelic diversity is low while population structure is higher for populations at the ends of the range in the West Atlantic and East Pacific.

Conclusions/Significance

These data are consistent with the proposed Indo-Pacific center of origin for S. lewini, and indicate that females are philopatric or adhere to coastal habitats while males facilitate gene flow across oceanic expanses. This study includes the largest sampling effort and the most molecular loci ever used to survey the complete range of a large oceanic predator, and findings emphasize the importance of incorporating mixed-marker analysis into stock assessments of threatened and endangered shark species.

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.