Extracting DNA from 'jaws': high yield and quality from archived tiger shark (Galeocerdo cuvier) skeletal material.

Archived specimens are highly valuable sources of DNA for retrospective genetic/genomic analysis. However, often limited effort has been made to evaluate and optimize extraction methods, which may be crucial for downstream applications. Here, we assessed and optimized the usefulness of abundant archived skeletal material from sharks as a source of DNA for temporal genomic studies. Six different methods for DNA extraction, encompassing two different commercial kits and three different protocols, were applied to material, so-called bio-swarf, from contemporary and archived jaws and vertebrae of tiger sharks (Galeocerdo cuvier). Protocols were compared for DNA yield and quality using a qPCR approach. For jaw swarf, all methods provided relatively high DNA yield and quality, while large differences in yield between protocols were observed for vertebrae. Similar results were obtained from samples of white shark (Carcharodon carcharias). Application of the optimized methods to 38 museum and private angler trophy specimens dating back to 1912 yielded sufficient DNA for downstream genomic analysis for 68% of the samples. No clear relationships between age of samples, DNA quality and quantity were observed, likely reflecting different preparation and storage methods for the trophies. Trial sequencing of DNA capture genomic libraries using 20 000 baits revealed that a significant proportion of captured sequences were derived from tiger sharks. This study demonstrates that archived shark jaws and vertebrae are potential high-yield sources of DNA for genomic-scale analysis. It also highlights that even for similar tissue types, a careful evaluation of extraction protocols can vastly improve DNA yield.

A novel field method to distinguish between cryptic carcharhinid sharks, Australian blacktip shark Carcharhinus tilstoni and common blacktip shark C. limbatus, despite the presence of hybrids.

Multivariate and machine-learning methods were used to develop field identification techniques for two species of cryptic blacktip shark. From 112 specimens, precaudal vertebrae (PCV) counts and molecular analysis identified 95 Australian blacktip sharks Carcharhinus tilstoni and 17 common blacktip sharks Carcharhinus limbatus. Molecular analysis also revealed 27 of the 112 were C. tilstoni × C. limbatus hybrids, of which 23 had C. tilstoni PCV counts and four had C. limbatus PCV counts. In the absence of further information about hybrid phenotypes, hybrids were assigned as either C. limbatus or C. tilstoni based on PCV counts. Discriminant analysis achieved 80% successful identification, but machine-learning models were better, achieving 100% successful identification, using six key measurements (fork length, caudal-fin peduncle height, interdorsal space, second dorsal-fin height, pelvic-fin length and pelvic-fin midpoint to first dorsal-fin insertion). Furthermore, pelvic-fin markings could be used for identification: C. limbatus has a distinct black mark >3% of the total pelvic-fin area, while C. tilstoni has markings with diffuse edges, or has smaller or no markings. Machine learning and pelvic-fin marking identification methods were field tested achieving 87 and 90% successful identification, respectively. With further refinement, the techniques developed here will form an important part of a multi-faceted approach to identification of C. tilstoni and C. limbatus and have a clear management and conservation application to these commercially important sharks. The methods developed here are broadly applicable and can be used to resolve species identities in many fisheries where cryptic species exist.

Can estimates of genetic effective population size contribute to fisheries stock assessments?

Sustainable exploitation of fisheries populations is challenging to achieve when the size of the population prior to exploitation and the actual numbers removed over time and across fishing zones are not clearly known. Quantitative fisheries' modeling is able to address this problem, but accurate and reliable model outcomes depend on high quality input data. Much of this information is obtained through the operation of the fishery under consideration, but while this seems appropriate, biases may occur. For example, poorly quantified changes in fishing methods that increase catch rates can erroneously suggest that the overall population size is increasing. Hence, the incorporation of estimates of abundance derived from independent data sources is preferable. We review and evaluate a fisheries-independent method of indexing population size; inferring adult abundance from estimates of the genetic effective size of a population (Ne ). Recent studies of elasmobranch species have shown correspondence between Ne and ecologically determined estimates of the population size (N). Simulation studies have flagged the possibility that the range of Ne /N ratios across species may be more restricted than previously thought, and also show that declines in Ne track declines in the abundance of model fisheries species. These key developments bring this new technology closer to implementation in fisheries science, particularly for data-poor fisheries or species of conservation interest.

Improved confidence intervals for the linkage disequilibrium method for estimating effective population size.

The linkage disequilibrium method is currently the most widely used single sample estimator of genetic effective population size. The commonly used software packages come with two options, referred to as the parametric and jackknife methods, for computing the associated confidence intervals. However, little is known on the coverage performance of these methods, and the published data suggest there may be some room for improvement. Here, we propose two new methods for generating confidence intervals and compare them with the two in current use through a simulation study. The new confidence interval methods tend to be conservative but outperform the existing methods for generating confidence intervals under certain circumstances, such as those that may be encountered when making estimates using large numbers of single-nucleotide polymorphisms.

NeEstimator v2: re-implementation of software for the estimation of contemporary effective population size (Ne ) from genetic data.

NeEstimator v2 is a completely revised and updated implementation of software that produces estimates of contemporary effective population size, using several different methods and a single input file. NeEstimator v2 includes three single-sample estimators (updated versions of the linkage disequilibrium and heterozygote-excess methods, and a new method based on molecular coancestry), as well as the two-sample (moment-based temporal) method. New features include the following: (i) an improved method for accounting for missing data; (ii) options for screening out rare alleles; (iii) confidence intervals for all methods; (iv) the ability to analyse data sets with large numbers of genetic markers (10 000 or more); (v) options for batch processing large numbers of different data sets, which will facilitate cross-method comparisons using simulated data; and (vi) correction for temporal estimates when individuals sampled are not removed from the population (Plan I sampling). The user is given considerable control over input data and composition, and format of output files. The freely available software has a new JAVA interface and runs under MacOS, Linux and Windows.

Stretched to the limit; can a short pelagic larval duration connect adult populations of an Indo-Pacific diadromous fish (Kuhlia rupestris)?

Freshwater species on tropical islands face localized extinction and the loss of genetic diversity. Their habitats can be ephemeral due to variability in freshwater run-off and erosion. Even worse, anthropogenic effects on these ecosystems are intense. Most of these species are amphidromous or catadromous (i.e. their life cycle includes a marine larval phase), which buffers them against many of these effects. A long pelagic larval duration (PLD) was thought to be critical to ensure the colonization and persistence in tropical islands, but recent findings indicated that several species with short PLDs are successful in those ecosystems. To test the potential of a short PLD in maintaining genetic connectivity and forestalling extirpation, we studied Kuhlia rupestris, a catadromous fish species with an extensive distribution in the western Pacific and Indian Oceans. Using a combination of molecular genetic markers (13 microsatellite loci and two gene regions from mtDNA) and modelling of larval dispersal, we show that a short PLD constrains genetic connectivity over a wide geographical range. Molecular markers showed that the short PLD did not prevent genetic divergence through evolutionary time and speciation has occurred or is occurring. Modelling of larvae dispersal suggested limited recent connectivity between genetically homogeneous populations across the Coral Sea. However, a short PLD can maintain connectivity on a subocean basin scale. Conservation and management of tropical diadromous species needs to take into account that population connectivity may be more limited than previously suspected in those species.

Accounting for missing data in the estimation of contemporary genetic effective population size (N(e) ).

Theoretical models are often applied to population genetic data sets without fully considering the effect of missing data. Researchers can deal with missing data by removing individuals that have failed to yield genotypes and/or by removing loci that have failed to yield allelic determinations, but despite their best efforts, most data sets still contain some missing data. As a consequence, realized sample size differs among loci, and this poses a problem for unbiased methods that must explicitly account for random sampling error. One commonly used solution for the calculation of contemporary effective population size (N(e) ) is to calculate the effective sample size as an unweighted mean or harmonic mean across loci. This is not ideal because it fails to account for the fact that loci with different numbers of alleles have different information content. Here we consider this problem for genetic estimators of contemporary effective population size (N(e) ). To evaluate bias and precision of several statistical approaches for dealing with missing data, we simulated populations with known N(e) and various degrees of missing data. Across all scenarios, one method of correcting for missing data (fixed-inverse variance-weighted harmonic mean) consistently performed the best for both single-sample and two-sample (temporal) methods of estimating N(e) and outperformed some methods currently in widespread use. The approach adopted here may be a starting point to adjust other population genetics methods that include per-locus sample size components.

Comparison of the reproductive ecology of two sympatric blacktip sharks (Carcharhinus limbatus and Carcharhinus tilstoni) off north-eastern Australia with species identification inferred from vertebral counts.

Precaudal vertebral counts were used to distinguish between 237 morphologically similar Carcharhinus limbatus and Carcharhinus tilstoni and were congruent with differences in reproductive ecology between the species. In addition to differing lengths at maturity and adult body size, the two species had asynchronous parturition, were born at different sizes and the relative frequencies of neonates differed in two coastal nursery areas. Despite evidence that hybridization can occur, these differences suggest the species are largely reproductively isolated.

Evidence for reproductive philopatry in the bull shark Carcharhinus leucas.

Reproductive philopatry in bull sharks Carcharhinus leucas was investigated by comparing mitochondrial (NADH dehydrogenase subunit 4, 797 base pairs and control region genes 837 base pairs) and nuclear (three microsatellite loci) DNA of juveniles sampled from 13 river systems across northern Australia. High mitochondrial and low microsatellite genetic diversity among juveniles sampled from different rivers (mitochondrial φ(ST) = 0·0767, P < 0·05; microsatellite F(ST) = -0·0022, P > 0·05) supported female reproductive philopatry. Genetic structure was not further influenced by geographic distance (P > 0·05) or long-shore barriers to movement (P > 0·05). Additionally, results suggest that C. leucas in northern Australia has a long-term effective population size of 11 000-13 000 females and has undergone population bottlenecks and expansions that coincide with the timing of the last ice-ages.

A review of the application of molecular genetics for fisheries management and conservation of sharks and rays.

Since the first investigation 25 years ago, the application of genetic tools to address ecological and evolutionary questions in elasmobranch studies has greatly expanded. Major developments in genetic theory as well as in the availability, cost effectiveness and resolution of genetic markers were instrumental for particularly rapid progress over the last 10 years. Genetic studies of elasmobranchs are of direct importance and have application to fisheries management and conservation issues such as the definition of management units and identification of species from fins. In the future, increased application of the most recent and emerging technologies will enable accelerated genetic data production and the development of new markers at reduced costs, paving the way for a paradigm shift from gene to genome-scale research, and more focus on adaptive rather than just neutral variation. Current literature is reviewed in six fields of elasmobranch molecular genetics relevant to fisheries and conservation management (species identification, phylogeography, philopatry, genetic effective population size, molecular evolutionary rate and emerging methods). Where possible, examples from the Indo-Pacific region, which has been underrepresented in previous reviews, are emphasized within a global perspective.

Genetic population structure of grey mackerel Scomberomorus semifasciatus in northern Australia.

This study used mtDNA sequence and microsatellite markers to elucidate the population structure of Scomberomorus semifasciatus collected from 12 widespread sampling locations in Australia. Samples (n = 544) were genotyped with nine microsatellite loci, and 353 were sequenced for the control (384 bp) and ATPase (800 bp) mtDNA gene regions. Combined interpretation of microsatellite and mtDNA data identified four genetic stocks of S. semifasciatus: Western Australia, north-west coast of the Northern Territory, Gulf of Carpentaria and the eastern coast of Queensland. Connectivity among stocks across northern Australia from the Northern Territory to the eastern coast of Queensland was high (mean F(ST) = 0·003 for the microsatellite data and Φ(ST) = 0·033 and 0·009 for control region and ATPase, respectively) leading to some uncertainty about stock boundaries. In contrast, there was a clear genetic break between the stock in Western Australia compared to the rest of northern Australia (mean F(ST) = 0·132 for the microsatellite data and Φ(ST) = 0·135 and 0·188 for control region and ATPase, respectively). This indicates a restriction to gene flow possibly associated with suboptimal habitat along the Kimberley coast (north Western Australia). The appropriate scale of management for this species corresponds to the jurisdictions of the three Australian states, except that authorities in Queensland and Northern Territory should co-ordinate the management of the Gulf of Carpentaria stock.

Population structure, effective population size and adverse effects of stocking in the endangered Australian eastern freshwater cod Maccullochella ikei.

Microsatellite markers were used to examine spatio-temporal genetic variation in the endangered eastern freshwater cod Maccullochella ikei in the Clarence River system, eastern Australia. High levels of population structure were detected. A model-based clustering analysis of multilocus genotypes identified four populations that were highly differentiated by F-statistics (F(ST) = 0·09 - 0·49; P < 0·05), suggesting fragmentation and restricted dispersal particularly among upstream sites. Hatchery breeding programmes were used to re-establish locally extirpated populations and to supplement remnant populations. Bayesian and frequency-based analyses of hatchery fingerling samples provided evidence for population admixture in the hatchery, with the majority of parental stock sourced from distinct upstream sites. Comparison between historical and contemporary wild-caught samples showed a significant loss of heterozygosity (21%) and allelic richness (24%) in the Mann and Nymboida Rivers since the commencement of stocking. Fragmentation may have been a causative factor; however, temporal shifts in allele frequencies suggest swamping with hatchery-produced M. ikei has contributed to the genetic decline in the largest wild population. This study demonstrates the importance of using information on genetic variation and population structure in the management of breeding and stocking programmes, particularly for threatened species.

IUCN classification zones concord with, but underestimate, the population genetic structure of the zebra shark Stegostoma fasciatum in the Indo-West Pacific.

The Indo-West Pacific (IWP), from South Africa in the western Indian Ocean to the western Pacific Ocean, contains some of the most biologically diverse marine habitats on earth, including the greatest biodiversity of chondrichthyan fishes. The region encompasses various densities of human habitation leading to contrasts in the levels of exploitation experienced by chondrichthyans, which are targeted for local consumption and export. The demersal chondrichthyan, the zebra shark, Stegostoma fasciatum, is endemic to the IWP and has two current regional International Union for the Conservation of Nature (IUCN) Red List classifications that reflect differing levels of exploitation: 'Least Concern' and 'Vulnerable'. In this study, we employed mitochondrial ND4 sequence data and 13 microsatellite loci to investigate the population genetic structure of 180 zebra sharks from 13 locations throughout the IWP to test the concordance of IUCN zones with demographic units that have conservation value. Mitochondrial and microsatellite data sets from samples collected throughout northern Australia and Southeast Asia concord with the regional IUCN classifications. However, we found evidence of genetic subdivision within these regions, including subdivision between locations connected by habitat suitable for migration. Furthermore, parametric F(ST) analyses and Bayesian clustering analyses indicated that the primary genetic break within the IWP is not represented by the IUCN classifications but rather is congruent with the Indonesian throughflow current. Our findings indicate that recruitment to areas of high exploitation from nearby healthy populations in zebra sharks is likely to be minimal, and that severe localized depletions are predicted to occur in zebra shark populations throughout the IWP region.

A PCR assay for gender assignment in dugong (Dugong dugon) and West Indian manatee (Trichechus manatus).

Gender assignment for some aquatic mammals in the field is difficult. Molecular sexing from tissue biopsies is possible as males are heterogametic. Here we describe a multiplex PCR assay that amplifies the male specific SRY gene and differentiates ZFX and ZFY gametologues in two sirenian species, dugong (Dugong dugon) and West Indian manatee (Trichechus manatus). The assay was validated with animals of known gender and proved accurate and robust to experimental failure.

Mitochondrial and allozyme genetics of incipient speciation in a landlocked population of Galaxias truttaceus (Pisces: Galaxiidae).

Galaxias truttaceus is found in coastal rivers and streams in south-eastern Australia. It spawns at the head of estuaries in autumn and the larvae spend 3 months of winter at sea before returning to fresh water. In Tasmania there are landlocked populations of G. truttaceus in a cluster of geologically young lakes on the recently glaciated Central Plateau. These populations have no marine larval stage and spawn in the lakes in spring. Speciation due to land locking is thought to be a frequent occurrence within Galaxias. To investigate the nature of the speciation event which may be occurring within lake populations of G. truttaceus we studied the mitochondrial DNA (mtDNA) and allozyme diversity of both lake and stream populations. Using the presence or absence of restriction sites recognized by 13 six-base restriction endonucleases, we found 58 mtDNA haplotypes among 150 fish collected from 13 Tasmanian and one south-east Australian mainland stream populations. The most parsimonious network relating the haplotypes by site loss or gain was starlike in shape. We argue that this arrangement is best explained by selection upon slightly beneficial mutations within the mitochondrial genome. Gene diversity analysis under Wright's island model showed that the populations in each drainage were not genetically subdivided. Only two of these stream haplotypes were found among the 66 fish analyzed from four lake populations. Despite the extreme lack of mtDNA diversity in lake populations, the observed nuclear DNA heterozygosity of 40 lake fish (0.10355) was only slightly less than that of 82 stream fish (0.11635). In the short time (3000-7000 years) that the lake fish have been landlocked, random genetic drift in a finite, stable-sized population was probably not responsible for the lack of mtDNA diversity in the lake populations. We infer the lake populations have probably experienced at least one, severe, but transitory bottleneck possibly induced by natural selection for life-history characters essential for survival in the lacustrine habitat. If speciation is occurring in the landlocked populations of G. truttaceus, then it may be driven by genetic transilience.