Retrospective genomics highlights changes in genetic composition of tiger sharks (Galeocerdo cuvier) and potential loss of a south-eastern Australia population.

Over the last century, many shark populations have declined, primarily due to overexploitation in commercial, artisanal and recreational fisheries. In addition, in some locations the use of shark control programs also has had an impact on shark numbers. Still, there is a general perception that populations of large ocean predators cover wide areas and therefore their diversity is less susceptible to local anthropogenic disturbance. Here we report on temporal genomic analyses of tiger shark (Galeocerdo cuvier) DNA samples that were collected from eastern Australia over the past century. Using Single Nucleotide Polymorphism (SNP) loci, we documented a significant change in genetic composition of tiger sharks born between ~1939 and 2015. The change was most likely due to a shift over time in the relative contribution of two well-differentiated, but hitherto cryptic populations. Our data strongly indicate a dramatic shift in the relative contribution of these two populations to the overall tiger shark abundance on the east coast of Australia, possibly associated with differences in direct or indirect exploitation rates.

Diving into the vertical dimension of elasmobranch movement ecology.

Knowledge of the three-dimensional movement patterns of elasmobranchs is vital to understand their ecological roles and exposure to anthropogenic pressures. To date, comparative studies among species at global scales have mostly focused on horizontal movements. Our study addresses the knowledge gap of vertical movements by compiling the first global synthesis of vertical habitat use by elasmobranchs from data obtained by deployment of 989 biotelemetry tags on 38 elasmobranch species. Elasmobranchs displayed high intra- and interspecific variability in vertical movement patterns. Substantial vertical overlap was observed for many epipelagic elasmobranchs, indicating an increased likelihood to display spatial overlap, biologically interact, and share similar risk to anthropogenic threats that vary on a vertical gradient. We highlight the critical next steps toward incorporating vertical movement into global management and monitoring strategies for elasmobranchs, emphasizing the need to address geographic and taxonomic biases in deployments and to concurrently consider both horizontal and vertical movements.

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.

Novel signature fatty acid profile of the giant manta ray suggests reliance on an uncharacterised mesopelagic food source low in polyunsaturated fatty acids.

Traditionally, large planktivorous elasmobranchs have been thought to predominantly feed on surface zooplankton during daytime hours. However, the recent application of molecular methods to examine long-term assimilated diets, has revealed that these species likely gain the majority from deeper or demersal sources. Signature fatty acid analysis (FA) of muscle tissue was used to examine the assimilated diet of the giant manta ray Mobula birostris, and then compared with surface zooplankton that was collected during feeding and non-feeding events at two aggregation sites off mainland Ecuador. The FA profiles of M. birostris and surface zooplankton were markedly different apart from similar proportions of arachidonic acid, which suggests daytime surface zooplankton may comprise a small amount of dietary intake for M. birostris. The FA profile of M. birostris muscle was found to be depleted in polyunsaturated fatty acids, and instead comprised high proportions of 18:1ω9 isomers. While 18:1ω9 isomers are not explicitly considered dietary FAs, they are commonly found in high proportions in deep-sea organisms, including elasmobranch species. Overall, the FA profile of M. birostris suggests a diet that is mesopelagic in origin, but many mesopelagic zooplankton species also vertically migrate, staying deep during the day and moving to shallower waters at night. Here, signature FA analysis is unable to resolve the depth at which these putative dietary items were consumed and how availability of this prey may drive distribution and movements of this large filter-feeder.

Population structure and connectivity of tiger sharks (Galeocerdo cuvier) across the Indo-Pacific Ocean basin.

Population genetic structure using nine polymorphic nuclear microsatellite loci was assessed for the tiger shark (Galeocerdo cuvier) at seven locations across the Indo-Pacific, and one location in the southern Atlantic. Genetic analyses revealed considerable genetic structuring (FST > 0.14, p < 0.001) between all Indo-Pacific locations and Brazil. By contrast, no significant genetic differences were observed between locations from within the Pacific or Indian Oceans, identifying an apparent large, single Indo-Pacific population. A lack of differentiation between tiger sharks sampled in Hawaii and other Indo-Pacific locations identified herein is in contrast to an earlier global tiger shark nDNA study. The results of our power analysis provide evidence to suggest that the larger sample sizes used here negated any weak population subdivision observed previously. These results further highlight the need for cross-jurisdictional efforts to manage the sustainable exploitation of large migratory sharks like G. cuvier.

Prey Density Threshold and Tidal Influence on Reef Manta Ray Foraging at an Aggregation Site on the Great Barrier Reef.

Large tropical and sub-tropical marine animals must meet their energetic requirements in a largely oligotrophic environment. Many planktivorous elasmobranchs, whose thermal ecologies prevent foraging in nutrient-rich polar waters, aggregate seasonally at predictable locations throughout tropical oceans where they are observed feeding. Here we investigate the foraging and oceanographic environment around Lady Elliot Island, a known aggregation site for reef manta rays Manta alfredi in the southern Great Barrier Reef. The foraging behaviour of reef manta rays was analysed in relation to zooplankton populations and local oceanography, and compared to long-term sighting records of reef manta rays from the dive operator on the island. Reef manta rays fed at Lady Elliot Island when zooplankton biomass and abundance were significantly higher than other times. The critical prey density threshold that triggered feeding was 11.2 mg m-3 while zooplankton size had no significant effect on feeding. The community composition and size structure of the zooplankton was similar when reef manta rays were feeding or not, with only the density of zooplankton changing. Higher zooplankton biomass was observed prior to low tide, and long-term (~5 years) sighting data confirmed that more reef manta rays are also observed feeding during this tidal phase than other times. This is the first study to examine prey availability at an aggregation site for reef manta rays and it indicates that they feed in locations and at times of higher zooplankton biomass.

Morphology of the ampullae of Lorenzini in juvenile freshwater Carcharhinus leucas.

Ampullae of Lorenzini were examined from juvenile Carcharhinus leucas (831-1,045 mm total length) captured from freshwater regions of the Brisbane River. The ampullary organ structure differs from all other previously described ampullae in the canal wall structure, the general shape of the ampullary canal, and the apically nucleated supportive cells. Ampullary pores of 140-205 µm in diameter are distributed over the surface of the head region with 2,681 and 2,913 pores present in two sharks that were studied in detail. The primary variation of the ampullary organs appears in the canal epithelial cells which occur as either flattened squamous epithelial cells or a second form of pseudostratified contour-ridged epithelial cells; both cell types appear to release material into the ampullary lumen. Secondarily, this ampullary canal varies due to involuted walls that form a clover-like canal wall structure. At the proximal end of the canal, contour-ridged cells abut a narrow region of cuboidal epithelial cells that verge on the constant, six alveolar sacs of the ampulla. The alveolar sacs contain numerous receptor and supportive cells bound by tight junctions and desmosomes. Pear-shaped receptor cells that possess a single apical kinocilium are connected basally by unmyelinated neural boutons. Opposed to previously described ampullae of Lorenzini, the supportive cells have an apical nucleus, possess a low number of microvilli, and form a unique, jagged alveolar wall. A centrally positioned centrum cap of cuboidal epithelial cells overlies a primary afferent lateral line nerve.

A comparative analysis of feeding and trophic level ecology in stingrays (Rajiformes; Myliobatoidei) and electric rays (Rajiformes: Torpedinoidei).

Standardised diets and trophic level (T L) estimates were calculated for 75 ray species from the suborders Myliobatoidei (67 spp.) and Torpedinoidei (8 spp.). Decapod crustaceans (31.71 ± 3.92%) and teleost fishes (16.45 ± 3.43%) made the largest contribution to the standardised diet of the Myliobatoidei. Teleost fishes (37.40 ± 16.09%) and polychaete worms (31.96 ± 14.22%) were the most prominent prey categories in the standardised diet of the suborder Torpedinoidei. Cluster analysis identified nine major trophic guilds the largest of which were decapod crustaceans (24 species), teleost fishes (11 species) and molluscs (11 species). Trophic level estimates for rays ranged from 3.10 for Potamotrygon falkneri to 4.24 for Gymnura australis, Torpedo marmorata and T. nobiliana. Secondary consumers with a T L <4.00 represented 84% of the species examined, with the remaining 12 species (16%) classified as tertiary consumers (T L ≥ 4.00). Tertiary consumers included electric rays (Torpedo, 3 spp. and Hypnos, 1 sp.), butterfly rays (Gymnura, 4 spp.), stingrays (2 spp.) and Potamotrygonid stingrays (2 spp.). Feeding strategies were identified as the primary factor of influence with respect to Myliobatoidei and Torpedinoidei T L estimates with inter-family comparisons providing the greatest insight into Myliobatoidei and Torpedinoidei relationships.