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

Comparative genomics of dusky kob (Argyrosomus japonicus, Sciaenidae) conspecifics: Evidence for speciation and the genetic mechanisms underlying traits

Dusky kob (Argyrosomus japonicus) is a commercially important finfish, indigenous to South Africa, Australia, and China. Previous studies highlighted differences in genetic composition, life history, and morphology of the species across geographic regions. A draft genome sequence of 0.742 Gb (N50 = 5.49 Mb; BUSCO completeness = 97.8%) and 22,438 predicted protein-coding genes was generated for the South African (SA) conspecific. A comparison with the Chinese (CN) conspecific revealed a core set of 32,068 orthologous protein clusters across both genomes. The SA genome exhibited 440 unique clusters compared to 1928 unique clusters in the CN genome. Transportation and immune response processes were overrepresented among the SA accessory genome, whereas the CN accessory genome was enriched for immune response, DNA transposition, and sensory detection (FDR-adjusted p < 0.01). These unique clusters may represent an adaptive component of the species' pangenome that could explain population divergence due to differential environmental specialisation. Furthermore, 700 single-copy orthologues (SCOs) displayed evidence of positive selection between the SA and CN genomes, and globally these genomes shared only 92% similarity, suggesting they might be distinct species. These genes primarily play roles in metabolism and digestion, illustrating the evolutionary pathways that differentiate the species. Understanding these genomic mechanisms underlying adaptation and evolution within and between species provides valuable insights into growth and maturation of kob, traits that are particularly relevant to commercial aquaculture.

Vertebral microchemistry as an indicator of habitat use of the oceanic whitetip shark Carcharhinus longimanus in the central and eastern Pacific Ocean

The oceanic whitetip shark, Carcharhinus longimanus, is a highly migratory, epipelagic top predator that is classified as critically endangered. Although this species is widely distributed throughout the world's tropical oceans, its assumed mobility and pelagic behavior limit studies to derive required lifetime data for management. To address this data deficiency, we assessed variation in the habitat use of C. longimanus by oceanic region and over ontogeny through time series trace element and stable isotope values conserved along the vertebral centra (within translucent annulus bands) of 13 individuals sampled from the central and eastern Pacific Ocean. Elemental ratios of Mg:Ca, Mn:Ca, Fe:Ca, Zn:Ca, and Ba:Ca varied significantly among individuals from both sampling regions while principal component analysis of combined standardized elements revealed minimal overlap between the two areas. The limited overlap was also in agreement with stable isotope niches. These findings indicate that C. longimanus exhibit a degree of fidelity to sampling regions but also connectivity in a proportion of the population. The relatively stable Sr:Ca ratio supports its occurrence in oceanic environments. The decreasing trends in Ba:Ca, Mn:Ca, and Zn:Ca ratios, as well as in carbon and nitrogen isotope values along vertebral transects, indicate that C. longimanus undergo a directional habitat shift with age. Combined elemental and stable isotope values in vertebral centra provide a promising tool for elucidating lifetime data for complex pelagic species. For C. longimanus, management will need to consider subpopulation movement behavior in the Pacific to minimize the potential for localized depletions. Further work is now required to sample individuals across the entire Pacific and to link these findings with genetic and movement data to define population structure.

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.

Population structure and genetic connectivity of the scalloped hammerhead shark (Sphyrna lewini) across nursery grounds from the Eastern Tropical Pacific: Implications for management and conservation

Defining demographically independent units and understanding patterns of gene flow between them is essential for managing and conserving exploited populations. The critically endangered scalloped hammerhead shark, Sphyrna lewini, is a coastal semi-oceanic species found worldwide in tropical and subtropical waters. Pregnant females give birth in shallow coastal estuarine habitats that serve as nursery grounds for neonates and small juveniles, whereas adults move offshore and become highly migratory. We evaluated the population structure and connectivity of S. lewini in coastal areas and one oceanic island (Cocos Island) across the Eastern Tropical Pacific (ETP) using both sequences of the mitochondrial DNA control region (mtCR) and 9 nuclear-encoded microsatellite loci. The mtCR defined two genetically discrete groups: one in the Mexican Pacific and another one in the central-southern Eastern Tropical Pacific (Guatemala, Costa Rica, Panama, and Colombia). Overall, the mtCR data showed low levels of haplotype diversity ranging from 0.000 to 0.608, while nucleotide diversity ranged from 0.000 to 0.0015. More fine-grade population structure was detected using microsatellite loci where Guatemala, Costa Rica, and Panama differed significantly. Relatedness analysis revealed that individuals within nursery areas were more closely related than expected by chance, suggesting that S. lewini may exhibit reproductive philopatric behaviour within the ETP. Findings of at least two different management units, and evidence of philopatric behaviour call for intensive conservation actions for this highly threatened species in the ETP.

A priority nursery area for the conservation of the scalloped hammerhead shark Sphyrna lewini in Mexico

This study describes the genetic diversity and level of genetic differentiation of the scalloped hammerhead shark Sphyrna lewini from eight putative nursery areas in the Mexican Pacific Ocean, using the mtDNA control region. Genetic analyses revealed a particular spatial divergence between La Reforma and all the remaining sites, with five exclusive haplotypes and the highest genetic diversity. This pattern may be interpreted as the signature of regional female philopatry, relatively to a particular female-mediated gene flow for La Reforma, which shows a strong subdivision in the Gulf of California.

Molecular identification and first demographic insights of sharks based on artisanal fisheries bycatch in the Pacific Coast of Colombia: implications for conservation

The Pacific coast of Colombia is characterized by mangrove ecosystems which play a crucial role as possible nurseries for juvenile sharks. However, trophic food webs from coastal ecosystems are heavily disturbed by increased fishing pressure, which affects numerous shark species. In this region of the Eastern Tropical Pacific (ETP), fisheries' data from coastal areas are scarce and unspecific, as most sharks from artisanal fisheries are landed decapitated and finless, making their morphological identification difficult. For the establishment and implementation of effective regional conservation and management policies, information on the diversity and population dynamics of shark species is crucial. We therefore sequenced the mitochondrial NADH2 gene of 696 samples taken from fishermen's landings of shark's bycatch along the Colombian north Pacific coast. We were able to identify 14 species of sharks, two of the most abundant species were Sphyrna lewini and Carcharhinus falciformis, both evaluated on IUCN the Red List of Threatened species (Critically Endangered and Vulnerable) and CITES regulated. We found low genetic diversity in the sampled area increasing the concern for both species in the region, even more considering that the majority of individuals were juveniles. Our results showed the importance of genetic markers for first population genetic insights as a complementary tool during the decision-making process in management plans. For this specific region, strategies such as the delimitation of conservation priority areas or the regulation of fishing gears could help improve the sustainability of shark populations in the Colombian Pacific.

Using stable isotopes analysis to understand ontogenetic trophic variations of the scalloped hammerhead shark at the Galapagos Marine Reserve

Changes in life-history requirements drive trophic variations, particularly in large marine predators. The life history of many shark species is still poorly known and understanding their dietary ontogeny is a challenging task, especially for highly migratory species. Stable isotope analysis has proven as a useful method for examining the foraging strategies of sharks and other marine predators. We assessed the foraging strategies and ontogenetic changes of scalloped hammerhead sharks, Sphyrna lewini, at Galapagos Marine Reserve (GMR), by analysing δ¹³C and δ¹⁵N signatures in different maturity stages. Our isotopic results suggest ontogenetic shifts in resource use between sub-adult and adult stages, but not between adult and juvenile stages. Carbon isotopic signatures found in the juvenile stage were enriched in contrast to sub-adults (~0.73‰) suggesting a combination of the maternal input and the use of coastal resources around the Galapagos Islands. Adult female sharks also showed enrichment in δ¹³C (~0.53‰) in comparison to sub-adult stages that suggest feeding in high primary productivity areas, such as the GMR. This study improves the understanding of the trophic ecology and ontogenetic changes of a highly migratory shark that moves across the protected and unprotected waters of the Eastern Tropical Pacific.

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

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

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

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

Arsenic in Tissues and Prey Species of the Scalloped Hammerhead (Sphyrna lewini) from the SE Gulf of California

The bioaccumulation of arsenic (As) in the muscle, liver, kidneys, and brain of the shark Sphyrna lewini was measured in 40 juvenile specimens from southeast Gulf of California. Additionally, the biomagnification factor was calculated through prey items from stomach contents of the analyzed specimens. The concentrations of As (mg kg^-1, wet weight) were higher in the muscle (10.1 ± 0.3) and liver (9.4 ± 0.5) than in the brain (4.5 ± 0.3) and kidneys (4.2 ± 0.2), which may be attributed to the biological functions of each tissue. Positive correlations were found between the levels of As in muscle and liver with the biological parameters of S. lewini. Hammerhead sharks feed mainly of teleost fishes with low As values (Clupeidae fishes, 1.1 ± 0.5; Sciaenidae fishes, 1.0 ± 0.6; Scomber japonicus, 1.2 ± 0.6; and Etropus crossotus 2.1 ± 0.4) compared with the predator, indicating biomagnification. Inorganic arsenic (As_i) in muscle was estimated as 3% of the total As, although muscle consumption is unlikely to represent a risk (HQ < 1) in humans. Moreover, the probabilities of developing cancer were estimated as low (3.99 × 10^-5 to 3.32 × 10^-6). To avoid health risks related to As, a weekly ration must not exceed 69.3 and 484.8 g in children and adults, respectively.

A globally threatened shark, Carcharias taurus, shows no population decline in South Africa

Knowledge about the demographic histories of natural populations helps to evaluate their conservation status, and potential impacts of natural and anthropogenic pressures. In particular, estimates of effective population size obtained through molecular data can provide useful information to guide management decisions for vulnerable populations. The spotted ragged-tooth shark, Carcharias taurus (also known as the sandtiger or grey nurse shark), is widely distributed in warm-temperate and subtropical waters, but has suffered severe population declines across much of its range as a result of overexploitation. Here, we used multilocus genotype data to investigate the demographic history of the South African C. taurus population. Using approximate Bayesian computation and likelihood-based importance sampling, we found that the population underwent a historical range expansion that may have been linked to climatic changes during the late Pleistocene. There was no evidence for a recent anthropogenic decline. Together with census data suggesting a stable population, these results support the idea that fishing pressure and other threats have so far not been detrimental to the local C. taurus population. The results reported here indicate that South Africa could possibly harbour the last remaining, relatively pristine population of this widespread but vulnerable top predator.

Conservation genetics of elasmobranchs of the Mexican Pacific Coast, trends and perspectives

One of the most critical threats to biodiversity is the high extinction rate driven by human activities. Reducing extinction rates requires the implementation of conservation programmes based on robust scientific data. Elasmobranchs are important ecological components of the ocean, and several species sustain substantial economic activities. Unfortunately, elasmobranchs are one of the most threatened and understudied animal taxa. The Mexican Pacific Coast (MPC) is a region with high elasmobranch diversity and is the seat of major elasmobranch fisheries. But it is also a developing region with several conservation and management challenges which require national and international attention. Here, we review the conservation genetics literature of elasmobranchs from the MPC. We present a synthesis of the works using samples from the region and emphasize the main gaps and biases in these data. In addition, we discuss the benefits and challenges of generating genomic information to improve the management and conservation of an elasmobranch biodiversity hotspot in a developing country. We found 47 elasmobranch genetic articles that cover <30% of the elasmobranch diversity in the region. These studies mainly used mitochondrial DNA sequences to analyse the genetic structure of commercially important and abundant species of the order Carcharhiniformes. Some of these papers also assessed mating systems, demographic parameters, and taxonomic uncertainties, all of which are important topics for efficient management decisions. In terms of conservation genetics, elasmobranchs from the MPC remain understudied. However, high-throughput sequencing technologies have increased the power and accessibility of genomic tools, even in developing countries such as Mexico. The tools described here provide information relevant for biodiversity conservation. Therefore, we strongly suggest that investment in genomic research will assist implementation of efficient management strategies. In time, this will reduce the extinction risk of the unique elasmobranch biodiversity from the MPC.

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.

Possible female philopatry of the smooth hammerhead shark Sphyrna zygaena revealed by genetic structure patterns

We assessed the spatial pattern of genetic structure of smooth hammerhead shark Sphyrna zygaena in 10 localities from the Northern Mexican Pacific. A total of 35 haplotypes were identified in 129 sequences of the mtDNA control region. The results showed slight but significant genetic structure among localities (Φ_ST = 0.044, P < 0.001). In addition, the localities with highest number of juveniles were genetically different (Φ_ST = 0.058, P < 0.024), which may be representative of nursery areas. The genetic differentiation pattern can be associated to female philopatry and preference for particular birthing sites. Finally, historical demography shows that S. zygaena populations present a recent demographic expansion that occurred during glacial events in the late Pleistocene to early Holocene.

Evolutionary insights from comparative transcriptome and transcriptome-wide coalescence analyses in Tetrastigma hemsleyanum

BACKGROUND

Tetrastigma hemsleyanum is of great medicinal importance and used as a model system to address the evolutionary history of warm-temperate evergreen (WTE) forest biomes in East Asia over Neogene time scales. However, further studies on the neutral and adaptive divergence processes of T. hemsleyanum are currently impeded by a lack of genomic resources. In this study, we de novo assembled and annotated a reference transcriptome for two cpDNA lineages (Central-South-East vs. Southwest) of T. hemsleyanum. We further used comparative genomic and multilocus coalescent approaches to investigate the tempo and mode of lineage diversification in T. hemsleyanum.

RESULTS

A total of 52,838 and 65,197 unigenes with an N50 of 1,667 and 1,841 bp for Central-South-East (CSE) and Southwest (SW) lineages, respectively, were recovered, and 6,692 putative orthologs were identified between the two lineages. Estimation of K_a/K_s ratios for these orthologs revealed that ten genes had K_a/K_s values significantly greater than 0.5 (P < 0.05), whereas 2,099 (K_a/K_s < 0.5, P < 0.05) were inferred to be under purifying selection. Based on three bioinformatic strategies, we identified a total of 1,018 single-copy nuclear genes (SCNGs) from the orthologs. We successfully designed eight nuclear gene primer pairs with high intraspecific variation (e.g. h_T = 0.923, π_T = 1.68×10^-3), when surveyed across a subset of T. hemsleyanum individuals. Concordant with the previous cpDNA data, the haplotype networks constructed for most nuclear gene loci clearly identified the two lineages. A multilocus coalescence analysis suggested that the separation between the two lineages appears to have occurred during the mid-Pliocene. Despite their ancient divergence, both lineages experienced expansion at rather localized scales and have continued to exchange genes at a low rate.

CONCLUSIONS

This study demonstrated the utility of transcriptome sequencing as a basis for SCNG development in non-model species and the advantages of integrating multiple nuclear loci for phylogeographic and phylogenetic studies.

Weak population structure of the Spot-tail shark Carcharhinus sorrah and the Blacktip shark C. limbatus along the coasts of the Arabian Peninsula, Pakistan, and South Africa

The increase in demand for shark meat and fins has placed shark populations worldwide under high fishing pressure. In the Arabian region, the spot-tail shark Carcharhinus sorrah and the Blacktip shark Carcharhinus limbatus are among the most exploited species. In this study, we investigated the population genetic structure of C. sorrah (n = 327) along the coasts of the Arabian Peninsula and of C. limbatus (n = 525) along the Arabian coasts, Pakistan, and KwaZulu-Natal, South Africa, using microsatellite markers (15 and 11 loci, respectively). Our findings support weak population structure in both species. Carcharhinus sorrah exhibited a fine structure, subdividing the area into three groups. The first group comprises all samples from Bahrain, the second from the UAE and Yemen, and the third from Oman. Similarly, C. limbatus exhibited population subdivision into three groups. The first group, comprising samples from Bahrain and Kuwait, was highly differentiated from the second and third groups, comprising samples from Oman, Pakistan, the UAE, and Yemen; and South Africa and the Saudi Arabian Red Sea, respectively. Population divisions were supported by pairwise F ST values and discriminant analysis of principal components (DAPC), but not by STRUCTURE. We suggest that the mostly low but significant pairwise F ST values in our study are suggestive of fine population structure, which is possibly attributable to behavioral traits such as residency in C. sorrah and site fidelity and philopatry in C. limbatus. However, for all samples obtained from the northern parts of the Gulf (Bahrain and/or Kuwait) in both species, the higher but significant pairwise F ST values could possibly be a result of founder effects during the Tethys Sea closure. Based on DAPC and F ST results, we suggest each population to be treated as independent management unit, as conservation concerns emerge.

Population genetics of the deep-sea bluntnose sixgill shark, Hexanchus griseus, revealing spatial genetic heterogeneity

Hexanchus griseus is a globally distributed deep-water shark species. It inhabits tropical and temperate waters throughout the world, including the Mediterranean Sea where it is by-caught by small-scale fisheries in the region. In this study, we analysed the genetic variation of H. griseus specimens collected from different areas within and outside the Mediterranean region, to assess its genetic connectivity. The mitochondrial DNA (mtDNA) sequence analysed in this study ranged from cytochrome b to 16S rRNA genes including the control region, the 12S rRNA gene and the interspersed tRNA genes in the region, covering a total of 3731 to 3914 nucleotides. Results have shown that this species exhibits geographically distinct maternal lineages, indicating population structure along geographical ranges. These findings reveal population subdivisions not only between the Pacific Ocean and the Atlantic Ocean, but also within the oceans and on a smaller scale within the Mediterranean Sea. This highlights the need to consider each population subdivision separately when designing management plans for the conservation of this species.

Crossing lines: a multidisciplinary framework for assessing connectivity of hammerhead sharks across jurisdictional boundaries

Conservation and management of migratory species can be complex and challenging. International agreements such as the Convention on Migratory Species (CMS) provide policy frameworks, but assessments and management can be hampered by lack of data and tractable mechanisms to integrate disparate datasets. An assessment of scalloped (Sphyrna lewini) and great (Sphyrna mokarran) hammerhead population structure and connectivity across northern Australia, Indonesia and Papua New Guinea (PNG) was conducted to inform management responses to CMS and Convention on International Trade in Endangered Species listings of these species. An Integrated Assessment Framework (IAF) was devised to systematically incorporate data across jurisdictions and create a regional synopsis, and amalgamated a suite of data from the Australasian region. Scalloped hammerhead populations are segregated by sex and size, with Australian populations dominated by juveniles and small adult males, while Indonesian and PNG populations included large adult females. The IAF process introduced genetic and tagging data to produce conceptual models of stock structure and movement. Several hypotheses were produced to explain stock structure and movement patterns, but more data are needed to identify the most likely hypothesis. This study demonstrates a process for assessing migratory species connectivity and highlights priority areas for hammerhead management and research.

Genetic structure and signatures of selection in grey reef sharks (Carcharhinus amblyrhynchos)

With overfishing reducing the abundance of marine predators in multiple marine ecosystems, knowledge of genetic structure and local adaptation may provide valuable information to assist sustainable management. Despite recent technological advances, most studies on sharks have used small sets of neutral markers to describe their genetic structure. We used 5517 nuclear single-nucleotide polymorphisms (SNPs) and a mitochondrial DNA (mtDNA) gene to characterize patterns of genetic structure and detect signatures of selection in grey reef sharks (Carcharhinus amblyrhynchos). Using samples from Australia, Indonesia and oceanic reefs in the Indian Ocean, we established that large oceanic distances represent barriers to gene flow, whereas genetic differentiation on continental shelves follows an isolation by distance model. In Australia and Indonesia differentiation at nuclear SNPs was weak, with coral reefs acting as stepping stones maintaining connectivity across large distances. Differentiation of mtDNA was stronger, and more pronounced in females, suggesting sex-biased dispersal. Four independent tests identified a set of loci putatively under selection, indicating that grey reef sharks in eastern Australia are likely under different selective pressures to those in western Australia and Indonesia. Genetic distances averaged across all loci were uncorrelated with genetic distances calculated from outlier loci, supporting the conclusion that different processes underpin genetic divergence in these two data sets. This pattern of heterogeneous genomic differentiation, suggestive of local adaptation, has implications for the conservation of grey reef sharks; furthermore, it highlights that marine species showing little genetic differentiation at neutral loci may exhibit patterns of cryptic genetic structure driven by local selection.

Reconstructing the demographic history of divergence between European river and brook lampreys using approximate Bayesian computations

Inferring the history of isolation and gene flow during species divergence is a central question in evolutionary biology. The European river lamprey (Lampetra fluviatilis) and brook lamprey (L. planeri) show a low reproductive isolation but have highly distinct life histories, the former being parasitic-anadromous and the latter non-parasitic and freshwater resident. Here we used microsatellite data from six replicated population pairs to reconstruct their history of divergence using an approximate Bayesian computation framework combined with a random forest model. In most population pairs, scenarios of divergence with recent isolation were outcompeted by scenarios proposing ongoing gene flow, namely the Secondary Contact (SC) and Isolation with Migration (IM) models. The estimation of demographic parameters under the SC model indicated a time of secondary contact close to the time of speciation, explaining why SC and IM models could not be discriminated. In case of an ancient secondary contact, the historical signal of divergence is lost and neutral markers converge to the same equilibrium as under the less parameterized model allowing ongoing gene flow. Our results imply that models of secondary contacts should be systematically compared to models of divergence with gene flow; given the difficulty to discriminate among these models, we suggest that genome-wide data are needed to adequately reconstruct divergence history.

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.

Blacktip reef sharks, Carcharhinus melanopterus, have high genetic structure and varying demographic histories in their Indo-Pacific range

For free-swimming marine species like sharks, only population genetics and demographic history analyses can be used to assess population health/status as baseline population numbers are usually unknown. We investigated the population genetics of blacktip reef sharks, Carcharhinus melanopterus; one of the most abundant reef-associated sharks and the apex predator of many shallow water reefs of the Indian and Pacific Oceans. Our sampling includes 4 widely separated locations in the Indo-Pacific and 11 islands in French Polynesia with different levels of coastal development. Four-teen microsatellite loci were analysed for samples from all locations and two mitochondrial DNA fragments, the control region and cytochrome b, were examined for 10 locations. For microsatellites, genetic diversity is higher for the locations in the large open systems of the Red Sea and Australia than for the fragmented habitat of the smaller islands of French Polynesia. Strong significant structure was found for distant locations with FST values as high as ~0.3, and a smaller but still significant structure is found within French Polynesia. Both mitochondrial genes show only a few mutations across the sequences with a dominant shared haplotype in French Polynesia and New Caledonia suggesting a common lineage different to that of East Australia. Demographic history analyses indicate population expansions in the Red Sea and Australia that may coincide with sea level changes after climatic events. Expansions and flat signals are indicated for French Polynesia as well as a significant recent bottleneck for Moorea, the most human-impacted lagoon of the locations in French Polynesia.

Can novel genetic analyses help to identify low-dispersal marine invasive species?

Genetic methods can be a powerful tool to resolve the native versus introduced status of populations whose taxonomy and biogeography are poorly understood. The genetic study of introduced species is presently dominated by analyses that identify signatures of recent colonization by means of summary statistics. Unfortunately, such approaches cannot be used in low-dispersal species, in which recently established populations originating from elsewhere in the species' native range also experience periods of low population size because they are founded by few individuals. We tested whether coalescent-based molecular analyses that provide detailed information about demographic history supported the hypothesis that a sea squirt whose distribution is centered on Tasmania was recently introduced to mainland Australia and New Zealand through human activities. Methods comparing trends in population size (Bayesian Skyline Plots and Approximate Bayesian Computation) were no more informative than summary statistics, likely because of recent intra-Tasmanian dispersal. However, IMa2 estimates of divergence between putatively native and introduced populations provided information at a temporal scale suitable to differentiate between recent (potentially anthropogenic) introductions and ancient divergence, and indicated that all three non-Tasmanian populations were founded during the period of European settlement. While this approach can be affected by inaccurate molecular dating, it has considerable (albeit largely unexplored) potential to corroborate nongenetic information in species with limited dispersal capabilities.

Recent demographic history and present fine-scale structure in the Northwest Atlantic leatherback (Dermochelys coriacea) turtle population

The leatherback turtle Dermochelys coriacea is the most widely distributed sea turtle species in the world. It exhibits complex life traits: female homing and migration, migrations of juveniles and males that remain poorly known, and a strong climatic influence on resources, breeding success and sex-ratio. It is consequently challenging to understand population dynamics. Leatherbacks are critically endangered, yet the group from the Northwest Atlantic is currently considered to be under lower risk than other populations while hosting some of the largest rookeries. Here, we investigated the genetic diversity and the demographic history of contrasted rookeries from this group, namely two large nesting populations in French Guiana, and a smaller one in the French West Indies. We used 10 microsatellite loci, of which four are newly isolated, and mitochondrial DNA sequences of the control region and cytochrome b. Both mitochondrial and nuclear markers revealed that the Northwest Atlantic stock of leatherbacks derives from a single ancestral origin, but show current genetic structuration at the scale of nesting sites, with the maintenance of migrants amongst rookeries. Low nuclear genetic diversities are related to founder effects that followed consequent bottlenecks during the late Pleistocene/Holocene. Most probably in response to climatic oscillations, with a possible influence of early human hunting, female effective population sizes collapsed from 2 million to 200. Evidence of founder effects and high numbers of migrants make it possible to reconsider the population dynamics of the species, formerly considered as a metapopulation model: we propose a more relaxed island model, which we expect to be a key element in the currently observed recovering of populations. Although these Northwest Atlantic rookeries should be considered as a single evolutionary unit, we stress that local conservation efforts remain necessary since each nesting site hosts part of the genetic diversity and species history.

Recovery potential and conservation options for elasmobranchs

Many elasmobranchs have experienced strong population declines, which have been largely attributed to the direct and indirect effects of exploitation. Recently, however, live elasmobranchs are being increasingly valued for their role in marine ecosystems, dive tourism and intrinsic worth. Thus, management plans have been implemented to slow and ultimately reverse negative trends, including shark-specific (e.g. anti-finning laws) to ecosystem-based (e.g. no-take marine reserves) strategies. Yet it is unclear how successful these measures are, or will be, given the degree of depletion and slow recovery potential of most elasmobranchs. Here, current understanding of elasmobranch population recoveries is reviewed. The potential and realized extent of population increases, including rates of increase, timelines and drivers are evaluated. Across 40 increasing populations, only 25% were attributed to decreased anthropogenic mortality, while the majority was attributed to predation release. It is also shown that even low exploitation rates (2-6% per year) can halt or reverse positive population trends in six populations currently managed under recovery plans. Management measures that help restore elasmobranch populations include enforcement or near-zero fishing mortality, protection of critical habitats, monitoring and education. These measures are highlighted in a case study from the south-eastern U.S.A., where some evidence of recovery is seen in Pristis pectinata, Galeocerdo cuvier and Sphyrna lewini populations. It is concluded that recovery of elasmobranchs is certainly possible but requires time and a combination of strong and dedicated management actions to be successful.