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

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

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

Genomic evidence indicates small island-resident populations and sex-biased behaviors of Hawaiian reef Manta Rays

BACKGROUND

Reef manta rays (Mobula alfredi) are globally distributed in tropical and subtropical seas. Their life history traits (slow growth, late maturity, low reproductive output) make them vulnerable to perturbations and therefore require informed management strategies. Previous studies have reported wide-spread genetic connectivity along continental shelves suggesting high gene flow along continuous habitats spanning hundreds of kilometers. However, in the Hawaiian Islands, tagging and photo-identification evidence suggest island populations are isolated despite proximity, a hypothesis that has not yet been evaluated with genetic data.

RESULTS

This island-resident hypothesis was tested by analyzing whole mitogenome haplotypes and 2048 nuclear single nucleotide polymorphisms (SNPs) between M. alfredi (n = 38) on Hawai'i Island and Maui Nui (the 4-island complex of Maui, Moloka'i, Lāna'i and Kaho'olawe). Strong divergence in the mitogenome (ΦST = 0.488) relative to nuclear genome-wide SNPs (neutral FST = 0.003; outlier FST = 0.186), and clustering of mitochondrial haplotypes among islands provides robust evidence that female reef manta rays are strongly philopatric and do not migrate between these two island groups. Combined with restricted male-mediated migration, equivalent to a single male moving between islands every 2.2 generations (~ 64 years), we provide evidence these populations are significantly demographically isolated. Estimates of contemporary effective population size (Ne) are 104 (95% CI: 99-110) in Hawai'i Island and 129 (95% CI: 122-136) in Maui Nui.

CONCLUSIONS

Concordant with evidence from photo identification and tagging studies, these genetic results indicate reef manta rays in Hawai'i have small, genetically-isolated resident island populations. We hypothesize that due to the Island Mass Effect, large islands provide sufficient resources to support resident populations, thereby making crossing deep channels separating island groups unnecessary. Small effective population size, low genetic diversity, and k-selected life history traits make these isolated populations vulnerable to region-specific anthropogenic threats, which include entanglement, boat strikes, and habitat degradation. The long-term persistence of reef manta rays in the Hawaiian Islands will require island-specific management strategies.

Exploring fishing threat at fleet segment and subregional scale: Least expert knowledge and a resilience versus disturbance-based approach as conservation's tools for cartilaginous fish

Based on an explorative but rigorous elicitation framework, we obtained the bycatch fishing probability at the fishing fleet segment level using expert estimates. Based on the knowledge of three scientific experts, we developed a new and creative structured method for smart and fast fishery-related risk assessments for species of high conservation concern. In order to test the method here propose, we applied it to 76 cartilaginous fish species (included in the IUCN Red Lists) and on five different fishing segments at both Italian and Mediterranean scale. The method produced qualitative results specific to the threat posed by fishing for each species and each segment with information between and within the segments. Based on the interpretation of resilience-disturbance interactions developed for ecological systems, the quantitative results provided reliable cumulative metrics, measuring the extinction risk due to fishing and the response to overfishing for the species considered. Additionally, the results highlight that the method perform best on a small geographic scale. Therefore, the application of this new method on other subregional or local scales where very few data are available (e.g., fishing effort) could be a valuable tool for the preliminary assessment for species of conservation concern. In fact, despite the absence of detailed catch data at local geographic scales, the flexibility of this method could help to highlight potential fishery-related conservation problems and thus redirect conservation strategies for threatened marine species such as many sharks and rays species.

Characterization of 35 new microsatellite markers for the blacktip reef shark (Carcharhinus melanopterus) and cross-species amplification in eight other shark species

BACKGROUND

Shark species are overfished at a global scale, as they are poached for the finning industry or are caught as bycatch. Efficient conservation measures require fine-scale spatial and temporal studies to characterize shark habitat use, infer migratory habits, analyze relatedness, and detect population genetic differentiation. Gathering these types of data is costly and time-consuming, especially when it requires collection of shark tissue samples.

METHODS AND RESULTS

Genetic tools, such as microsatellite markers, are the most economical sampling method for collecting genetic data, as they enable the estimation of genetic diversity, population structure and parentage relationships and are thus an efficient way to inform conservation strategies. Here, a set of 45 microsatellite loci was tested on three blacktip reef shark (Carcharhinus melanopterus) populations from three Polynesian islands: Moorea, Morane and Tenararo. The set was composed of 10 previously published microsatellite markers and 35 microsatellite markers that were developed specifically for C. melanopterus as part of the present study. The 35 novel and 10 existing loci were cross-amplified on eight additional shark species (Carcharhinus amblyrhynchos, C. longimanus, C. sorrah, Galeocerdo cuvier, Negaprion acutidens, Prionacea glauca, Rhincodon typus and Sphyrna lewini). These species had an average of 69% of successful amplification, considered if at least 50% of the individual samples being successfully amplified per species and per locus.

CONCLUSIONS

This novel microsatellite marker set will help address numerous knowledge gaps that remain, concerning genetic stock identification, shark behavior and reproduction via parentage analysis.

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.

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.

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.

Lack of multiple paternity in the oceanodromous tiger shark (Galeocerdo cuvier)

Multiple paternity has been documented as a reproductive strategy in both viviparous and ovoviviparous elasmobranchs, leading to the assumption that multiple mating may be ubiquitous in these fishes. However, with the majority of studies conducted on coastal and nearshore elasmobranchs that often form mating aggregations, parallel studies on pelagic, semi-solitary species are lacking. The tiger shark (Galeocerdo cuvier) is a large pelagic shark that has an aplacental viviparous reproductive mode which is unique among the carcharhinids. A total of 112 pups from four pregnant sharks were genotyped at nine microsatellite loci to assess the possibility of multiple paternity or polyandrous behaviour by female tiger sharks. Only a single pup provided evidence of possible multiple paternity, but with only seven of the nine loci amplifying for this individual, results were inconclusive. In summary, it appears that the tiger sharks sampled in this study were genetically monogamous. These findings may have implications for the genetic diversity and future sustainability of this population.

Variability in multiple paternity rates for grey reef sharks (Carcharhinus amblyrhynchos) and scalloped hammerheads (Sphyrna lewini)

This study assessed the presence and prevalence of multiple paternity (MP) in litters of grey reef sharks (Carcharhinus amblyrhynchos) and scalloped hammerheads (Sphyrna lewini) opportunistically caught in Papua New Guinea (PNG). Litter size between species were significantly different with an average of 3.3 pups for grey reef sharks and 17.2 pups for scalloped hammerhead. Using 14 and 10 microsatellite loci respectively, we identified MP in 66% of grey reef sharks (4 out of 6 litters) and 100% MP in scalloped hammerheads (5 litters). We found high paternal skew (the uneven contribution of sires per litter) and a positive correlation between female adult size and litter size in scalloped hammerheads but not in grey reef sharks. Differences in the frequency of MP between species and the identification of paternal skew may be linked with mating strategies and post-copulatory mechanisms. Multiple paternity is thought to benefit populations by enhancing genetic diversity therefore increasing the population’s genetic resilience to extrinsic pressures. The identification of MP in two shark species reported here, further elucidates the complex breeding strategies elasmobranchs undertake.

Switch from sexual to parthenogenetic reproduction in a zebra shark

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

Assessing multiple paternity in three commercially exploited shark species: Mustelus mustelus, Carcharhinus obscurus and Sphyrna lewini

In this study, multiple paternity (MP) was investigated in three commercially important shark species, common smoothhound Mustelus mustelus, dusky shark Carcharhinus obscurus and scalloped hammerhead Sphyrna lewini occurring in southern Africa. Reduced marker panels of between five and six microsatellite loci were constructed for each species and used to genotype and assess the presence of MP in a total of 60 M. mustelus individuals from six litters, 90 C. obscurus individuals from 14 litters and 54 S. lewini individuals from 13 litters. Analysis in GERUD and COLONY revealed the presence of MP in all three species. Multiple paternities were observed in 67, 35 and 46% of the litters of M. mustelus, C. obscurus and S. lewini, with corresponding average sire size of 1·6, 1·4 and 2·0, respectively. The variation in the rate of MP among the three species is in accordance with previous studies whilst the comparatively high frequency of MP observed for M. mustelus, matches what has previously been reported for shark species demonstrating aggregation behaviour.

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

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

A case of genetic polyandry in the shortfin mako Isurus oxyrinchus

This article documents a case of genetic polyandry in the oceanic and pelagic shortfin mako Isurus oxyrinchus and briefly comments on the implications of this finding.

Multiple paternity and hybridization in two smooth-hound sharks

Multiple paternity appears to be a common trait of elasmobranch mating systems, with its occurrence likely driven by convenience, due to females seeking to minimize the stress of male harassment. Here we use molecular markers to analyse the frequency of multiple paternity in two related viviparous sharks, Mustelus mustelus and Mustelus punctulatus. We first applied molecular methods to assign pregnant females, embryos and additional reference adults (N = 792) to one of the two species. Paternity analysis was performed using a total of 9 polymorphic microsatellites on 19 females and 204 embryos of M. mustelus, and on 13 females and 303 embryos of M. punctulatus. Multiple paternity occurs in both species, with 47% of M. mustelus and 54% of M. punctulatus litters sired by at least two fathers. Female fecundity is not influenced by multiple mating and in 56% of polyandrous litters paternity is skewed, with one male siring most of the pups. Genetic analyses also revealed hybridization between the two species, with a M. punctulatus female bearing pups sired by a M. mustelus male. The frequency of polyandrous litters in these species is consistent with aspects of their reproductive biology, such as synchronous ovulation and possible occurrence of breeding aggregations.

Insights into the antiviral immunity against grass carp (Ctenopharyngodon idella) reovirus (GCRV) in grass carp

Global fish production from aquaculture has rapidly grown over the past decades, and grass carp shares the largest portion. However, hemorrhagic disease caused by grass carp reovirus (GCRV) results in tremendous loss of grass carp (Ctenopharyngodon idella) industry. During the past years, development of molecular biology and cellular biology technologies has promoted significant advances in the understanding of the pathogen and the immune system. Immunoprophylaxis based on stimulation of the immune system of fish has also got some achievements. In this review, authors summarize the recent progresses in basic researches on GCRV; viral nucleic acid sensors, high-mobility group box proteins (HMGBs); pattern recognition receptors (PRRs), Toll-like receptors (TLRs) and retinoic acid inducible gene I- (RIG-I-) like receptors (RLRs); antiviral immune responses induced by PRRs-mediated signaling cascades of type I interferon (IFN-I) and IFN-stimulated genes (ISGs) activation. The present review also notices the potential applications of molecule genetic markers. Additionally, authors discuss the current preventive and therapeutic strategies (vaccines, RNAi, and prevention medicine) and highlight the importance of innate immunity in long term control for grass carp hemorrhagic disease.

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.

Frequency of multiple paternity in gummy shark, Mustelus antarcticus, and rig, Mustelus lenticulatus, and the implications of mate encounter rate, postcopulatory influences, and reproductive mode

Elasmobranchs (sharks and rays) show an amazing diversity of reproductive modes and behaviors. Multiple paternity (MP) has been identified in all species where more than 1 litter has been investigated; yet neither direct nor indirect benefits from MP have been determined in elasmobranchs. This has led to the suggestion that MP in this group may simply be a product of convenience polyandry with variation in the frequency of MP driven by differences in mate encounter rates. Here, we use molecular markers to investigate polyandry and MP in 2 closely related and commercially important species of shark, Mustelus antarcticus and Mustelus lenticulatus. In total, 328 M. antarcticus embryos originating from 29 different mothers and 75 M. lenticulatus embryos originating from 19 different mothers were genotyped using 8 microsatellite loci. We find that MP occurs in both species. However, in both species, the majority of litters were sired by a single father. Our results do not support increased fecundity per se as a driver of MP. Further, our results do not suggest that high population densities with resulting high mate encounter rates generated by breeding aggregations necessarily lead to high frequencies of MP. Importantly, we note evidence of reproductive skew within polyandrous litters, which is a predicted outcome of postcopulatory mechanisms.

Direct genetic evidence for reproductive philopatry and associated fine-scale migrations in female blacktip reef sharks (Carcharhinus melanopterus) in French Polynesia

Conservation of top predators has been emphasized as essential in an ecosystem due to their role in trophic chain regulation. Optimizing conservation strategies for these endangered marine top predators requires direct estimates of breeding patterns and connectivity as these are essential to understanding the population dynamics. There have been some attempts to investigate breeding patterns of reef sharks from litter reconstruction using molecular analyses. However, direct fine-scale migrations of female sharks for parturition as well as connectivity at a medium scale like between islands remain mostly unknown. We used microsatellite DNA markers and a likelihood-based parentage analysis to determine breeding patterns of female blacktip reef sharks in Moorea (Society Islands, French Polynesia). Most females gave birth at their home island but some migrated to specific nursery areas outside the area they are attached to, sometimes going to another island 50 km away across deep ocean. Our analysis also revealed that females migrated to the same nursery for every birthing event. Many offspring showed a high level of inbreeding indicating an overall reduced population size, restricted movements and dispersal, or specific mating behaviour. Females represent the vectors that transport the genes at nursery grounds, and their fidelity should thus define reproductive units. As females seem to be philopatric, males could be the ones dispersing genes between populations. These results highlight the need to conserve coastal zones where female reef sharks seem to exhibit philopatry during the breeding season.