Flexibility for fuelling reproduction in a pelagic ray (Mobula eregoodoo) suggested by bioenergetic modelling.

This study investigated the measurements of energy density and bioenergetic modelling for a pelagic ray, Mobula eregoodoo, to estimate its relative allocation to various bodily processes and especially reproduction. The data revealed M. eregoodoo uses up to 21.0% and 2.5% of its annual energy budget on growth and reproduction, respectively. During pregnancy, females depleted energy reserves in the liver, which, along with their biennial reproductive cycle, aligns with general theory that ectotherms are capital breeders and thus build energy reserves before reproduction. Nonetheless, the reduction in energy reserves did not account for all reproductive costs, and therefore, gravid females supplement reproductive costs through energy derived from the diet, according to an income-breeding strategy. These characteristics imply that M. eregoodoo exhibits some flexibility in fuelling reproduction depending on energy availability throughout the reproductive cycle, which may be prevalent in other elasmobranchs. The data represent the first estimates of both the metabolic costs of gestation in elasmobranchs and the relative cost of reproduction in rays. Energy costs and plasticity associated with highly variable reproductive strategies in elasmobranchs may influence long-term population viability under a rapidly changing environment.

Reef manta rays forage on tidally driven, high density zooplankton patches in Hanifaru Bay, Maldives.

Manta rays forage for zooplankton in tropical and subtropical marine environments, which are generally nutrient-poor. Feeding often occurs at predictable locations where these large, mobile cartilaginous fishes congregate to exploit ephemeral productivity hotspots. Investigating the zooplankton dynamics that lead to such feeding aggregations remains a key question for understanding their movement ecology. The aim of this study is to investigate the feeding environment at the largest known aggregation for reef manta rays Mobula alfredi in the world. We sampled zooplankton throughout the tidal cycle, and recorded M. alfredi activity and behaviour, alongside environmental variables at Hanifaru Bay, Maldives. We constructed generalised linear models to investigate possible relationships between zooplankton dynamics, environmental parameters, and how they influenced M. alfredi abundance, behaviour, and foraging strategies. Zooplankton biomass changed rapidly throughout the tidal cycle, and M. alfredi feeding events were significantly related to high zooplankton biomass. Mobula alfredi switched from non-feeding to feeding behaviour at a prey density threshold of 53.7 mg dry mass m-3; more than double the calculated density estimates needed to theoretically meet their metabolic requirements. The highest numbers of M. alfredi observed in Hanifaru Bay corresponded to when they were engaged in feeding behaviour. The community composition of zooplankton was different when M. alfredi was feeding (dominated by copepods and crustaceans) compared to when present but not feeding (more gelatinous species present than in feeding samples). The dominant zooplankton species recorded was Undinula vulgaris. This is a large-bodied calanoid copepod species that blooms in oceanic waters, suggesting offshore influences at the site. Here, we have characterised aspects of the feeding environment for M. alfredi in Hanifaru Bay and identified some of the conditions that may result in large aggregations of this threatened planktivore, and this information can help inform management of this economically important marine protected area.

Delimiting cryptic species within the brown-banded bamboo shark, Chiloscyllium punctatum in the Indo-Australian region with mitochondrial DNA and genome-wide SNP approaches.

BACKGROUND: Delimiting cryptic species in elasmobranchs is a major challenge in modern taxonomy due the lack of available phenotypic features. Employing stand-alone genetics in splitting a cryptic species may prove problematic for further studies and for implementing conservation management. In this study, we examined mitochondrial DNA and genome-wide nuclear single nucleotide polymorphisms (SNPs) in the brown-banded bambooshark, Chiloscyllium punctatum to evaluate potential cryptic species and the species-population boundary in the group. RESULTS: Both mtDNA and SNP analyses showed potential delimitation within C. punctatum from the Indo-Australian region and consisted of four operational taxonomic units (OTUs), i.e. those from Indo-Malay region, the west coast of Sumatra, Lesser Sunda region, and the Australian region. Each OTU can be interpreted differently depending on available supporting information, either based on biological, ecological or geographical data. We found that SNP data provided more robust results than mtDNA data in determining the boundary between population and cryptic species. CONCLUSION: To split a cryptic species complex and erect new species based purely on the results of genetic analyses is not recommended. The designation of new species needs supportive diagnostic morphological characters that allow for species recognition, as an inability to recognise individuals in the field creates difficulties for future research, management for conservation and fisheries purposes. Moreover, we recommend that future studies use a comprehensive sampling regime that encompasses the full range of a species complex. This approach would increase the likelihood of identification of operational taxonomic units rather than resulting in an incorrect designation of new species.

Mutualism promotes site selection in a large marine planktivore.

Mutualism is a form of symbiosis whereby both parties benefit from the relationship. An example is cleaning symbiosis, which has been observed in terrestrial and marine environments. The most recognized form of marine cleaning symbiosis is that of cleaner fishes and their clients.Cleaner species set up cleaning stations on the reef, and other species seek out their services. However, it is not well understood how the presence of cleaning stations influence movements of large highly mobile species. We examined the role of cleaning stations as a driver of movement and habitat use in a mobile client species.Here, we used a combination of passive acoustic telemetry and in-water surveys to investigate cleaning station attendance by the reef manta ray Mobula alfredi. We employed a novel approach in the form of a fine-scale acoustic receiver array set up around a known cleaning area and tagged 42 rays. Within the array, we mapped structural features, surveyed the distribution of cleaner wrasse, and observed the habitat use of the rays.We found manta ray space use was significantly associated with blue-streak cleaner wrasse Labroides dimidiatus distribution and hard coral substrate. Cleaning interactions dominated their habitat use at this site, taking precedence over other life history traits such as feeding and courtship.This study has demonstrated that cleaning symbiosis is a driver for highly mobile, and otherwise pelagic, species to visit inshore reef environments. We suggest that targeted and long-term use of specific cleaning stations reflects manta rays having a long-term memory and cognitive map of some shallow reef environments where quality cleaning is provided. We hypothesize that animals prefer cleaning sites in proximity to productive foraging regions.

Age and growth of the tropical oviparous shark, Chiloscyllium punctatum in Indonesian waters.

The brown-banded bamboo shark, Chiloscyllium punctatum, is the most common shark caught in coastal commercial fisheries throughout Southeast Asia, yet there is a lack of the life-history information necessary for reliable stock assessments. The authors estimated growth rates and age at maturity using analysis of growth bands in vertebral centra. They trialled four different techniques to enhance the visibility and improve identification of the putative annual growth bands necessary for age estimation. The authors found that the burn method on whole vertebral centra provided the most readable and consistent results for age analysis. The logistic model was chosen as the best-fit growth model for age estimation of 330 individual C. punctatum from Indonesia. Several age verification methods, including marginal increment ratio and length-frequency analysis, were performed with the support of age validation through the use of calcein-labelled vertebrae from two sharks maintained in captivity. This study found that C. punctatum from Indonesian waters is a fast-growing species that can grow up to 18 cm year-1 , reach an estimated maximum total length of 1 m, mature at c. 6.5 years and live for up to 14 years.

Short- and long-term diets of the threatened longhorned pygmy devil ray, Mobula eregoodoo determined using stable isotopes.

Most mobulids are listed as near threatened to endangered. Nonetheless, effective conservation measures are hindered by knowledge gaps in their ecology and behaviour. In particular, few studies have assessed diets and trophic ecologies that could inform methods to avoid fishing mortality. Here, a shortfall in data for the longhorned pygmy devil ray, Mobula eregoodoo was addressed by describing temporal variability in dietary preferences using stable isotope analysis. During summer and autumn in 2017, five bather-protection gillnets were deployed off eastern Australia (29° S, 153.5° E). From the catches of these gillnets, 35 adult M. eregoodoo had liver, muscle and stomach contents sampled to determine δ13 C and δ15 N profiles. Analyses revealed that surface zooplankton and zooplanktivorous teleosts were important dietary components across short- and long-term temporal scales. Large quantities of undigested sandy sprat, Hyperlophus vittatus, in the stomachs of some specimens unequivocally confirm feeding on teleosts. A narrow isotopic niche and minimal isotopic overlap with reef manta rays, Mobula alfredi from the same geographic region in eastern Australia implies M. eregoodoo has unique and highly specialised resource use relative to other mobulids in the area. The species is clearly vulnerable to capture during inshore migrations, presumably where they feed on shallow-water shoaling teleosts. Female M. eregoodoo likely have a low annual reproductive output, so population recoveries from fishing-induced declines are likely to be slow. Measures to reduce the by catch of M. eregoodoo in local bather-protection gillnets, and artisanal fisheries more broadly, should be given priority.

The geographic distribution of reef and oceanic manta rays (Mobula alfredi and Mobula birostris) in Australian coastal waters.

The known distribution of manta rays in Australian waters is patchy, with records primarily centred around tourism hotspots. We collated 11,614 records of Mobula alfredi from photo-ID databases (n = 10,715), aerial surveys (n = 378) and online reports (n = 521). The study confirms an uninterrupted coastal distribution from north of 26°S and 31°S on the west and east coasts, respectively. More southerly M. alfredi records relate to warm-water events with a southernmost extent at 34°S. Coastal sightings of Mobula birostris were rare (n = 32), likely reflecting a preference for offshore waters, but encompass a wider latitudinal extent than M. alfredi of 10-40°S.

Powering Ocean Giants: The Energetics of Shark and Ray Megafauna.

Shark and ray megafauna have crucial roles as top predators in many marine ecosystems, but are currently among the most threatened vertebrates and, based on historical extinctions, may be highly susceptible to future environmental perturbations. However, our understanding of their energetics lags behind that of other taxa. Such knowledge is required to answer important ecological questions and predict their responses to ocean warming, which may be limited by expanding ocean deoxygenation and declining prey availability. To develop bioenergetics models for shark and ray megafauna, incremental improvements in respirometry systems are useful but unlikely to accommodate the largest species. Advances in biologging tools and modelling could help answer the most pressing ecological questions about these iconic species.

Development and characterization of 17 polymorphic microsatellite markers for the reef manta ray (Mobula alfredi).

OBJECTIVE: Limited sample sizes are often a problem for species of conservation concern when using genetic tools to make population assessments. Lack of analytical power from small sample sizes can be compensated for by use of a large marker set. Here we report on development and characterization of 17 novel microsatellite markers for the reef manta ray (Mobula alfredi). RESULTS: Loci were screened on 60 reef manta rays (M. alfredi) sampled from the east coast of Australia. The number of alleles per locus varied from 2 to 13 with observed heterozygosities ranging between 0.300 and 0.917. The development of these 17 additional markers increases the total number of microsatellite markers available for this species to 27.

Comparative morphology of the electrosensory system of the epaulette shark Hemiscyllium ocellatum and brown-banded bamboo shark Chiloscyllium punctatum.

We compared the electrosensory system of two benthic elasmobranchs Hemiscyllium ocellatum and Chiloscyllium punctatum. The distribution of the ampullary pores on the head was similar for both species, with a higher density of pores anteriorly and a lower density posteriorly, although C. punctatum generally possessed larger pores. Ampullary canals of the mandibular cluster were quasi-sinusoidal in H. ocellatum, a shape previously found in benthic rays only, whereas ampullary canals in C. punctatum were of a linear morphology as reported for many shark and ray species previously. The ampullae proper were of the lobular type, as occurs in most galean sharks. Chiloscyllium punctatum had six sensory chambers compared with the five per ampulla in H. ocellatum, which were generally smaller than those of C. punctatum. The sensory epithelium comprised flattened receptor cells, compared with the usual pear-shaped receptor cells encountered in other elasmobranchs and their apically nucleated supportive cells did not protrude markedly into the ampullary lumen, unlike those in benthic rays.

Tetraphyllidean and onchoproteocephalidean cestodes of elasmobranchs from Moreton Bay, Australia: description of two new species and new records for seven described species.

Parasitological examination of elasmobranchs of Moreton Bay, Queensland, Australia, resulted in the discovery of cestodes belonging to several armed genera of the Tetraphyllidea and Onchoproteocephalidea. Two new tetraphyllideans, Yorkeria moretonensis n. sp. and Yorkeria williamsi n. sp., are described from Chiloscyllium cf. punctatum (Hemiscylliidae). Yorkeria moretonensis n. sp. differs from its congeners in the possession of vitelline follicles that are discontinuous in the region of the ovary and in the length of its pedicels. Yorkeria williamsi n. sp. is most similar to Y. parva Southwell, 1927, but has larger, oval bothridia, longer pedicels and differences in the sizes of the scolex hooks. Yorkeria longstaffae Caira, Jensen & Rajan, 2007 is reported from Moreton Bay for the first time, and Spiniloculus mavensis Southwell, 1925 is re-reported from the type-locality and likely type-host (Moreton Bay and Chiloscyllium cf. punctatum, respectively), over 90 years after its original description. Six known onchoproteocephalideans, Acanthobothrium cannoni Campbell & Beveridge, 2002, A. chisholmae Campbell & Beveridge, 2002, A. ocallaghani Campbell & Beveridge, 2002, A. margieae Fyler, 2011, Megalonchos shawae Caira, Reyda & Mega, 2007 and M. sumansinghai Caira, Reyda & Mega, 2007, are reported from Moreton Bay for the first time, representing significant range extensions for all species.

Satellite tagging highlights the importance of productive Mozambican coastal waters to the ecology and conservation of whale sharks.

The whale shark Rhincodon typus is an endangered, highly migratory species with a wide, albeit patchy, distribution through tropical oceans. Ten aerial survey flights along the southern Mozambican coast, conducted between 2004-2008, documented a relatively high density of whale sharks along a 200 km stretch of the Inhambane Province, with a pronounced hotspot adjacent to Praia do Tofo. To examine the residency and movement of whale sharks in coastal areas around Praia do Tofo, where they may be more susceptible to gill net entanglement, we tagged 15 juveniles with SPOT5 satellite tags and tracked them for 2-88 days (mean = 27 days) as they dispersed from this area. Sharks travelled between 10 and 2,737 km (mean = 738 km) at a mean horizontal speed of 28 ± 17.1 SD km day-1. While several individuals left shelf waters and travelled across international boundaries, most sharks stayed in Mozambican coastal waters over the tracking period. We tested for whale shark habitat preferences, using sea surface temperature, chlorophyll-a concentration and water depth as variables, by computing 100 random model tracks for each real shark based on their empirical movement characteristics. Whale sharks spent significantly more time in cooler, shallower water with higher chlorophyll-a concentrations than model sharks, suggesting that feeding in productive coastal waters is an important driver of their movements. To investigate what this coastal habitat choice means for their conservation in Mozambique, we mapped gill nets during two dedicated aerial surveys along the Inhambane coast and counted gill nets in 1,323 boat-based surveys near Praia do Tofo. Our results show that, while whale sharks are capable of long-distance oceanic movements, they can spend a disproportionate amount of time in specific areas, such as along the southern Mozambique coast. The increasing use of drifting gill nets in this coastal hotspot for whale sharks is likely to be a threat to regional populations of this iconic species.

Patterns of specificity and diversity in species of Paraorygmatobothrium Ruhnke, 1994 (Cestoda: Phyllobothriidae) in Moreton Bay, Queensland, Australia, with the description of four new species.

A survey of tapeworms of galeomorph sharks from Moreton Bay (Queensland, Australia) identified a complex of species of Paraorygmatobothrium Ruhnke, 1994 infecting 11 carcharhiniform and two orectolobiform species. Combined morphological and multi-locus molecular analyses (based on the 28S nuclear ribosomal RNA and partial mitochondrial NADH dehydrogenase subunit 1 genes) revealed the presence of 12 species of Paraorygmatobothrium; four species (Paraorygmatobothrium christopheri n. sp., P. harti n. sp., P. sinclairtaylori n. sp. and P. ullmanni n. sp.) are considered to be new to science and are formally described, four represent known species, and four lack sufficient morphological data to allow definitive identification. In contrast to previous records for the genus, four of the species found in this study exhibited low host specificity [P. orectolobi (Butler, 1987) Ruhnke, 2011, P. sinclairtaylori, P. ullmanni and Paraorygmatobothrium sp. 3], three stenoxenic species were each found in two closely-related sharks (P. orectolobi, P. ullmanni and Paraorygmatobothrium sp. 3) and one euryxenic species was found in five species from two shark families (P. sinclairtaylori). One species was found to exhibit mild morphologically plasticity (P. orectolobi), with size range being associated with different shark species. Conversely, collections of almost morphologically indistinguishable specimens from single shark species were found to represent multiple species of Paraorygmatobothrium. The findings of this study indicate that the description of species of this genus on the basis of morphological data alone is problematic and that the inclusion of multi-locus molecular data is essential for future work on Paraorygmatobothrium. Host specificity, morphology and phylogenetic relatedness of species of Paraorygmatobothrium are explored.

Towards sustainable fishery management for skates in South America: The genetic population structure of Zearaja chilensis and Dipturus trachyderma (Chondrichthyes, Rajiformes) in the south-east Pacific Ocean.

The longnose skates (Zearaja chilensis and Dipturus trachyderma) are the main component of the elasmobranch fisheries in the south-east Pacific Ocean. Both species are considered to be a single stock by the fishery management in Chile however, little is known about the level of demographic connectivity within the fishery. In this study, we used a genetic variation (560 bp of the control region of the mitochondrial genome and ten microsatellite loci) to explore population connectivity at five locations along the Chilean coast. Analysis of Z. chilensis populations revealed significant genetic structure among off-shore locations (San Antonio, Valdivia), two locations in the Chiloé Interior Sea (Puerto Montt and Aysén) and Punta Arenas in southern Chile. For example, mtDNA haplotype diversity was similar across off-shore locations and Punta Arenas (h = 0.46-0.50), it was significantly different to those in the Chiloé Interior Sea (h = 0.08). These results raise concerns about the long-term survival of the species within the interior sea, as population resilience will rely almost exclusively on self-recruitment. In contrast, little evidence of genetic structure was found for D. trachyderma. Our results provide evidence for three management units for Z. chilensis, and we recommend that separate management arrangements are required for each of these units. However, there is no evidence to discriminate the extant population of Dipturus trachyderma as separate management units. The lack of genetic population subdivision for D. trachyderma appears to correspond with their higher dispersal ability and more offshore habitat preference.

Environmental influences and ontogenetic differences in vertical habitat use of black marlin (Istiompax indica) in the southwestern Pacific.

The black marlin (Istiompax indica) is a highly migratory billfish that occupies waters throughout the tropical and subtropical Indo-Pacific. To characterize the vertical habitat use of I. indica, we examined the temperature-depth profiles collected using 102 pop-up satellite archival tags deployed off the east coast of Australia. Modelling of environmental variables revealed location, sea-surface height deviation, mixed layer depth and dissolved oxygen to all be significant predictors of vertical habitat use. Distinct differences in diel movements were observed between the size classes, with larger size classes of marlin (greater than 50 kg) undertaking predictable bounce-diving activity during daylight hours, while diving behaviour of the smallest size class occurred randomly during both day and night. Overall, larger size classes of I. indica were found to use an increased thermal range and spend more time in waters below 150 m than fish of smaller size classes. The differences in the diving behaviour among size classes were suggested to reflect ontogenetic differences in foraging behaviour or physiology. The findings of this study demonstrate, for the first time to our knowledge, ontogenetic differences in vertical habitat in a species of billfish, and further the understanding of pelagic fish ecophysiology in the presence of global environmental change.

Manta birostris, predator of the deep? Insight into the diet of the giant manta ray through stable isotope analysis.

The characterization of diet for the giant manta ray Manta birostris has been problematic given their large-scale movement patterns and the difficulty in obtaining stomach contents from this species. The large majority of existing information is based on observational data limited to feeding events at the sea surface during daylight. Recently discovered aggregation sites for the giant manta ray off mainland Ecuador are some of the most accessible to date and provide a unique opportunity for researchers to gather much needed information on this elusive species. To assess how important surface zooplankton is to giant manta ray diet, we conducted stable isotope analysis (15N and 13C) on M. birostris muscle and surface zooplankton. Trophic position estimates placed M. birostris overall at a secondary consumer level of approximately 3.4 but there was large variation in δ15N and δ13C values among individuals. Manta birostris muscle tissue δ13C values were also not consistent with this species feeding predominantly on surface zooplankton and suggest that the majority of dietary intake is of mesopelagic origin. Given the conservative life history and fisheries pressure on large planktivores, knowledge of their trophic role and foraging strategies is essential to better understand their ecology and develop effective conservation measures.

The complete validated mitochondrial genome of the yellownose skate Zearaja chilensis (Guichenot 1848) (Rajiformes, Rajidae).

The yellownose skate Zearaja chilensis is endemic to South America. The species is the target of a valuable commercial fishery in Chile, but is highly susceptible to over-exploitation. The complete mitochondrial genome was described from 694,593 sequences obtained using Ion Torrent Next Generation Sequencing. The total length of the mitogenome was 16,909 bp, comprising 2 rRNAs, 13 protein-coding genes, 22 tRNAs and 2 non-coding regions. Comparison between the proposed mitogenome and one previously described from "raw fish fillets from a skate speciality restaurant in Seoul, Korea" resulted in 97.4% similarity, rather than approaching 100% similarity as might be expected. The 2.6% dissimilarity may indicate the presence of two separate stocks or two different species of, ostensibly, Z. chilensis in South America and highlights the need for caution when using genetic resources without a taxonomic reference or a voucher specimen.

Genetype and phylogenomic position of the frilled shark Chlamydoselachus anguineus inferred from the mitochondrial genome.

The genetype of Chlamydoselachus anguineus from Australian waters is described using the whole mitochondrial genome obtained from Illumina NGS technology. Total length of the mitogenome is 17 313 bp, consisting of 2 rRNAs, 13 protein-coding genes, 22 tRNA genes and 2 non-coding regions thus updating the previously available mitogenome for this species. The phylogenomic reconstruction comprising all available species of Superorder Squalomorphi supports the inclusion of C. anguineus in a divergent clade inside Order Hexanchiformes. Phyletic relationships inferred from the whole mitochondrial genomes are in agreement with traditional taxonomy. The low divergence between C. anguineus genomes (>99.9% genetic identity) is consistent with a widespread population in the west Pacific Ocean.

The phylogenetic position of the giant devil ray Mobula mobular (Bonnaterre, 1788) (Myliobatiformes, Myliobatidae) inferred from the mitochondrial genome.

The giant devil ray, Mobula mobular, is a member of one of the most distinct groups of cartilaginous fishes, the Mobulidae (manta and devil rays), and is the only mobulid assessed as Endangered due its restricted distribution, high bycatch mortality and suspected population decline. The complete mitochondrial genome is 18 913 base pairs in length and comprises 2 rRNAs, 13 protein-coding genes, 22 tRNAs and 2 non-coding regions. Comparison with the partial mitogenome of M. japanica suggests a sister-cryptic species complex and two different taxonomic units. However, the limited divergence within the species (>99.9% genetic identity) may be the result of a geographically and numerically restricted population of M. mobular within the Mediterranean Sea.

The phylogenomic position of the grey nurse shark Carcharias taurus Rafinesque, 1810 (Lamniformes, Odontaspididae) inferred from the mitochondrial genome.

The complete mitochondrial genome of the grey nurse shark Carcharias taurus is described from 25 963 828 sequences obtained using Illumina NGS technology. Total length of the mitogenome is 16 715 bp, consisting of 2 rRNAs, 13 protein-coding regions, 22 tRNA and 2 non-coding regions thus updating the previously published mitogenome for this species. The phylogenomic reconstruction inferred from the mitogenome of 15 species of Lamniform and Carcharhiniform sharks supports the inclusion of C. taurus in a clade with the Lamnidae and Cetorhinidae. This complete mitogenome contributes to ongoing investigation into the monophyly of the Family Odontaspididae.