Effective number of white shark (Carcharodon carcharias, Linnaeus) breeders is stable over four successive years in the population adjacent to eastern Australia and New Zealand

Population size is a central parameter for conservation; however, monitoring abundance is often problematic for threatened marine species. Despite substantial investment in research, many marine species remain data-poor presenting barriers to the evaluation of conservation management outcomes and the modeling of future solutions. Such is the case for the white shark (Carcharodon carcharias), a highly mobile apex predator for whom recent and substantial population declines have been recorded in many globally distributed populations. Here, we estimate the effective number of breeders that successfully contribute offspring in one reproductive cycle (Nb) to provide a snapshot of recent reproductive effort in an east Australian-New Zealand population of white shark. Nb was estimated over four consecutive age cohorts (2010, 2011, 2012, and 2013) using two genetic estimators (linkage disequilibrium; LD and sibship assignment; SA) based on genetic data derived from two types of genetic markers (single nucleotide polymorphisms; SNPs and microsatellite loci). While estimates of Nb using different marker types produced comparable estimates, microsatellite loci were the least precise. The LD and SA estimates of Nb within cohorts using SNPs were comparable; for example, the 2013 age cohort Nb(SA) was 289 (95% CI 200-461) and Nb(LD) was 208.5 (95% CI 116.4-712.7). We show that over the time period studied, Nb was stable and ranged between 206.1 (SD ± 45.9) and 252.0 (SD ± 46.7) per year using a combined estimate of Nb(LD+SA) from SNP loci. In addition, a simulation approach showed that in this population the effective population size (Ne) per generation can be expected to be larger than Nb per reproductive cycle. This study demonstrates how breeding population size can be monitored over time to provide insight into the effectiveness of recovery and conservation measures for the white shark, where the methods described here may be applicable to other data-poor species of conservation concern.

Sweepstakes reproductive success is absent in a New Zealand snapper (Chrysophrus auratus) population protected from fishing despite "tiny" Ne /N ratios elsewhere

A landmark study published in 2002 estimated a very small N_e /N ratio (around 10^-5 ) in a population of pink snapper (Chrysophrys auratus, Forster, 1801) in the Hauraki Gulf in New Zealand. It epitomized the tiny N_e /N ratios (<10^-3 ) reported in marine species due to the hypothesized operation of sweepstakes reproductive success (SRS). Here we re-evaluate the occurrence of SRS in marine species and the potential effect of fishing on the N_e /N ratio by studying the same species in the same region, but in a population that has been protected from fishing since 1975. We combine empirical, simulation and model-based approaches to estimate N_e (and N_b ) from genotypes of 1,044 adult fish and estimate N using recapture-probabilities. The estimated N_e /N ratio was much larger (0.33, SE: 0.14) than expected. The magnitude of estimates of population-wide variance in individual lifetime reproductive success (10-18) suggested that the sweepstakes effect was negligible in the study population. After evaluating factors that could explain the contrast between studies - experimental design, life history differences, environmental effects and the influence of exploitation on the N_e /N ratio - we conclude that the low N_e of the Hauraki Gulf population is associated with demographic instability in the harvested compared to the protected population despite circumstantial evidence that the 2002 study may have underestimated N_e . This study has broad implications for the prevailing view that reproductive success in the sea is largely driven by chance, and for genetic monitoring of populations using the N_e /N ratio and N_b .

Population genetics of Southern Hemisphere tope shark (Galeorhinus galeus): Intercontinental divergence and constrained gene flow at different geographical scales

The tope shark (Galeorhinus galeus Linnaeus, 1758) is a temperate, coastal hound shark found in the Atlantic and Indo-Pacific oceans. In this study, the population structure of Galeorhinus galeus was determined across the entire Southern Hemisphere, where the species is heavily targeted by commercial fisheries, as well as locally, along the South African coastline. Analysis was conducted on a total of 185 samples using 19 microsatellite markers and a 671 bp fragment of the NADH dehydrogenase subunit 2 (ND2) gene. Across the Southern Hemisphere, three geographically distinct clades were recovered, including one from South America (Argentina, Chile), one from Africa (all the South African collections) and an Australia-New Zealand clade. Nuclear data revealed significant population subdivisions (F_ST = 0.192 to 0.376, p<0.05) indicating limited gene flow for tope sharks across ocean basins. Marked population connectivity was however evident across the Indian Ocean based on Bayesian clustering analysis. More locally in South Africa, F-statistics and multivariate analysis supported moderate to high gene flow across the Atlantic/Indian Ocean boundary (F_ST = 0.035 to 0.044, p<0.05), with exception of samples from Struisbaai and Port Elizabeth which differed significantly from the rest. Discriminant and Bayesian clustering analysis indicated admixture in all sampling populations, decreasing from west to east, corroborating possible restriction to gene flow across regional oceanographic barriers. Mitochondrial sequence data recovered seven haplotypes (h = 0.216, π = 0.001) for South Africa, with one major haplotype shared by 87% of the individuals and at least one private haplotype for each sampling location except Port Elizabeth. As with many other coastal shark species with cosmopolitan distribution, this study confirms the lack of both historical dispersal and inter-oceanic gene flow while also implicating contemporary factors such as oceanic currents and thermal fronts to drive local genetic structure of G. galeus on a smaller spatial scale.

Application of environmental DNA to detect an endangered marine skate species in the wild

Environmental DNA (eDNA) techniques have only recently been applied in the marine environment to detect the presence of marine species. Species-specific primers and probes were designed to detect the eDNA of the endangered Maugean skate (Zearaja maugeana) from as little as 1 L of water collected at depth (10–15 m) in Macquarie Harbour (MH), Tasmania. The identity of the eDNA was confirmed as Z. maugeana by sequencing the qPCR products and aligning these with the target sequence for a 100% match. This result has validated the use of this eDNA technique for detecting a rare species, Z. maugeana, in the wild. Being able to investigate the presence, and possibly the abundance, of Z. maugeana in MH and Bathurst harbour (BH), would be addressing a conservation imperative for the endangered Z. maugeana. For future application of this technique in the field, the rate of decay was determined for Z. maugeana eDNA under ambient dissolved oxygen (DO) levels (55% saturation) and lower DO (20% saturation) levels, revealing that the eDNA can be detected for 4 and 16 hours respectively, after which eDNA concentration drops below the detection threshold of the assay. With the rate of decay being influenced by starting eDNA concentrations, it is recommended that samples be filtered as soon as possible after collection to minimize further loss of eDNA prior to and during sample processing.

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.

Characterization, development and multiplexing of microsatellite markers in three commercially exploited reef fish and their application for stock identification

Thirty-four microsatellite loci were isolated from three reef fish species; golden snapper Lutjanus johnii, blackspotted croaker Protonibea diacanthus and grass emperor Lethrinus laticaudis using a next generation sequencing approach. Both IonTorrent single reads and Illumina MiSeq paired-end reads were used, with the latter demonstrating a higher quality of reads than the IonTorrent. From the 1–1.5 million raw reads per species, we successfully obtained 10–13 polymorphic loci for each species, which satisfied stringent design criteria. We developed multiplex panels for the amplification of the golden snapper and the blackspotted croaker loci, as well as post-amplification pooling panels for the grass emperor loci. The microsatellites characterized in this work were tested across three locations of northern Australia. The microsatellites we developed can detect population differentiation across northern Australia and may be used for genetic structure studies and stock identification.

Population Genetic Diversity in the Australian 'Seascape': A Bioregion Approach

Genetic diversity within species may promote resilience to environmental change, yet little is known about how such variation is distributed at broad geographic scales. Here we develop a novel Bayesian methodology to analyse multi-species genetic diversity data in order to identify regions of high or low genetic diversity. We apply this method to co-distributed taxa from Australian marine waters. We extracted published summary statistics of population genetic diversity from 118 studies of 101 species and > 1000 populations from the Australian marine economic zone. We analysed these data using two approaches: a linear mixed model for standardised data, and a mixed beta-regression for unstandardised data, within a Bayesian framework. Our beta-regression approach performed better than models using standardised data, based on posterior predictive tests. The best model included region (Integrated Marine and Coastal Regionalisation of Australia (IMCRA) bioregions), latitude and latitude squared. Removing region as an explanatory variable greatly reduced model performance (delta DIC 23.4). Several bioregions were identified as possessing notably high genetic diversity. Genetic diversity increased towards the equator with a 'hump' in diversity across the range studied (-9.4 to -43.7°S). Our results suggest that factors correlated with both region and latitude play a role in shaping intra-specific genetic diversity, and that bioregion can be a useful management unit for intra-specific as well as species biodiversity. Our novel statistical model should prove useful for future analyses of within species genetic diversity at broad taxonomic and geographic scales.

The complete mitochondrial genome of the black jewfish Protonibea diacanthus (Perciformes: Sciaenidae)

We describe the complete mitochondrial genome of the black Jewfish Protonibea diacanthus. It was assembled from approximately 1.6 million reads produced by Ion Torrent next generation sequencing. The complete genome was 16,521 bp in length consisting of 13 protein-coding regions, 22 tRNA, 12S and 16S rRNA as well as two non-coding regions. The A+T base content (52.8%) is similar to other teleosts.

The complete mitochondrial genome of the golden snapper Lutjanus johnii (Perciformes: Lutjanidae)

We describe the complete mitochondrial genome of the golden snapper Lutjanus johnii. It was assembled from approximately 1.4 million reads produced by Ion Torrent next generation sequencing. The complete genome was 16,596 bp in length consisting of 13 protein-coding regions, 22 tRNA, 12S and 16S rRNA as well as two non-coding regions. The A+T base content (52.8%) is similar to other teleosts.

Differential use of Alcian blue and toluidine blue dyes for the quantification and isolation of anionic glycoconjugates from cell cultures: application to proteoglycans and a high-molecular-weight glycoprotein synthesized by articular chondrocytes

Alcian blue and toluidine blue dyes form complexes with anionic glycoconjugates (AG) such as proteoglycans (PG) and glycosaminoglycans (GAG). However, the Alcian blue-AG complexes do not readily dissociate, while the toluidine blue-AG complexes do so in salt solutions. This differential dissociation of the dye-AG complexes has been utilized in the analysis and isolation of radiolabeled AG elaborated by articular chondrocyte cultures incubated with the radiolabeled precursors of AG. For the rapid quantification of newly synthesized (35)S-labeled PG, small replicate aliquots of the radiolabeled culture media were applied directly to cellulose acetate strips, stained with Alcian blue and the stained immobilized radiolabeled PG was quantified by liquid scintillation counting. Comparison of anionic glycoconjugates quantified in the culture media employing toluidine blue and Alcian blue staining on cellulose acetate trips gave similar results. Staining on cellulose acetate strips using these two dyes is particularly suited for the simultaneous processing of large numbers of samples, as illustrated by the screening of the effects of biological materials and drugs on AG synthesis, in cultures labeled with [(35)S]-sulfate and [(3)H]-glucosamine. The Alcian blue and toluidine blue precipitation methods yielded similar results for the total AG recovered from the media of TGF-beta-stimulated chondrocytes. Electrophoretic analysis of toluidine blue- and Alcian blue-precipitated AG followed by autoradiography and Alcian blue staining in combination with silver nitrate demonstrated that both dyes yielded similar pattern of bands on gels. However, some AG from Alcian blue precipitate did not enter the gel, suggesting incomplete dissociation of Alcian blue-AG complex. The application of the toluidine blue precipitation method, in combination with enzymatic digestion of the GAG chains of the PGs, is illustrated by the isolation of a non-PG high-molecular-weight AG, as well as the PGs from the media of chondrocyte cultures stimulated by TGF-beta.