Misidentification of istiophorid billfishes by fisheries observers raises uncertainty over stock status.

We investigated whether commercially landed black marlin Istiompax indica were being misidentified by fisheries observers operating throughout the Pacific Ocean. Of 83 samples reported by observers as I. indica, 77% were genetically identified to be blue marlin Makaira nigricans and 2% to be striped marlin Kajikia audax. The high rate of misidentification by observers places considerable uncertainty over historic catch ratios of Indo-Pacific marlin and stock assessments relying on the validity of these data.

Evolution of diadromy in fish: insights from a tropical genus (Kuhlia species).

Diadromous species undergo regular migration between fresh and marine waters. This behavior is found in many species, including fish, mollusks, and crustaceans, some of which are commercially valuable species. Several attempts to trace the evolution of this behavior have been made in Salmonidae and Galaxiidae, but ambiguous phylogenies and multiple character state changes prevented unequivocal conclusions. The Kuhliidae family consists of 12 fish species that inhabit tropical islands in the Indo-Pacific region. The species have marine, partially catadromous, or fully catadromous life histories (i.e., they migrate from rivers to the sea to reproduce). The evolution of migratory behavior was traced on a well-resolved phylogeny. Catadromous Kuhlia species were basal, and partially catadromous and marine species formed derived monophyletic groups. This is, to our knowledge, the first time that a clear origin and polarity for the diadromous character has been demonstrated. We propose that the relative lack of resources in tropical, inshore, marine habitats and the ephemeral and isolated nature of freshwater environments of tropical islands, combined with phenotypic plasticity of migratory traits, play key roles in driving the evolution of diadromy in the Kuhliidae and probably in other groups. This work is an important starting point to understand the role of diadromy in speciation and adaptation in unstable habitats.

The genetic signature of recent speciation in manta rays (Manta alfredi and M. birostris).

Manta rays have been taxonomically revised as two species, Manta alfredi and M. birostris, on the basis of morphological and meristic data, yet the two species occur in extensive mosaic sympatry. We analysed the genetic signatures of the species boundary using a portion of the nuclear RAG1 (681 base pairs), mitochondrial CO1 (574 bp) and ND5 genes (1188 bp). The assay with CO1 sequences, widely used in DNA barcoding, failed to distinguish the two species. The two species were clearly distinguishable, however, with no shared RAG1 or ND5 haplotypes. The species were reciprocally monophyletic for RAG1, but paraphyletic for ND5 sequences. Qualitative evidence and statistical inferences using the 'Isolation-with-Migration models' indicated that these results were better explained with post-divergence gene flow in the recent past rather than incomplete lineage sorting with zero gene flow since speciation. An estimate of divergence time was less than 0.5 Ma with an upper confidence limit of within 1 Ma. Recent speciation of highly mobile species in the marine environment is of great interest, as it suggests that speciation may have occurred in the absence of long-term physical barriers to gene flow. We propose that the ecologically driven forces such as habitat choice played a significant role in speciation in manta rays.

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.

Fishing for DNA? Designing baits for population genetics in target enrichment experiments: Guidelines, considerations and the new tool supeRbaits.

Targeted sequencing is an increasingly popular next-generation sequencing (NGS) approach for studying populations that involves focusing sequencing efforts on specific parts of the genome of a species of interest. Methodologies and tools for designing targeted baits are scarce but in high demand. Here, we present specific guidelines and considerations for designing capture sequencing experiments for population genetics for both neutral genomic regions and regions subject to selection. We describe the bait design process for three diverse fish species: Atlantic salmon, Atlantic cod and tiger shark, which was carried out in our research group, and provide an evaluation of the performance of our approach across both historical and modern samples. The workflow used for designing these three bait sets has been implemented in the R-package supeRbaits, which encompasses our considerations and guidelines for bait design for the benefit of researchers and practitioners. The supeRbaits R-package is user-friendly and versatile. It is written in C++ and implemented in R. supeRbaits and its manual are available from Github: https://github.com/BelenJM/supeRbaits.

Likelihood-based genetic mark-recapture estimates when genotype samples are incomplete and contain typing errors.

Genotypes produced from samples collected non-invasively in harsh field conditions often lack the full complement of data from the selected microsatellite loci. The application to genetic mark-recapture methodology in wildlife species can therefore be prone to misidentifications leading to both 'true non-recaptures' being falsely accepted as recaptures (Type I errors) and 'true recaptures' being undetected (Type II errors). Here we present a new likelihood method that allows every pairwise genotype comparison to be evaluated independently. We apply this method to determine the total number of recaptures by estimating and optimising the balance between Type I errors and Type II errors. We show through simulation that the standard error of recapture estimates can be minimised through our algorithms. Interestingly, the precision of our recapture estimates actually improved when we included individuals with missing genotypes, as this increased the number of pairwise comparisons potentially uncovering more recaptures. Simulations suggest that the method is tolerant to per locus error rates of up to 5% per locus and can theoretically work in datasets with as little as 60% of loci genotyped. Our methods can be implemented in datasets where standard mismatch analyses fail to distinguish recaptures. Finally, we show that by assigning a low Type I error rate to our matching algorithms we can generate a dataset of individuals of known capture histories that is suitable for the downstream analysis with traditional mark-recapture methods.

Spatial and temporal genetic variation in an exploited reef fish: The effects of exploitation on cohort genetic structure.

Many coral reef fishes are fished, often resulting in detrimental genetic effects; however, reef fishes often show unpredictable patterns of genetic variation, which potentially mask the effects of fishing. Our goals were to characterize spatial and temporal genetic variation and determine the effects of fishing on an exploited reef fish, Plectropomus leopardus, Lacepède (the common coral trout). To determine population structure, we genotyped 417 Great Barrier Reef coral trout from four populations sampled in 2 years (1996 and 2004) at nine microsatellite loci. To test for exploitation effects, we additionally genotyped 869 individuals from a single cohort (ages 3-5) across eight different reefs, including fished and control populations. Genetic structure differed substantially in the two sampled years, with only 1 year exhibiting isolation by distance. Thus, genetic drift likely plays a role in shaping population genetic structure in this species. Although we found no loss of genetic diversity associated with exploitation, our relatedness patterns show that pulse fishing likely affects population genetics. Additionally, genetic structure in the cohort samples likely reflected spatial variation in recruitment contributing to genetic structure at the population level. Overall, we show that fishing does impact coral reef fishes, highlighting the importance of repeated widespread sampling to accurately characterize the genetic structure of reef fishes, as well as the power of analysing cohorts to avoid the impacts of recruitment-related genetic swamping. The high temporal and spatial variability in genetic structure, combined with possible selection effects, will make conservation/management of reef fish species complex.

Ocean currents and the population genetic signature of fish migrations.

Animal migrations are a fascinating and global phenomenon, yet they are often difficult to study and sometimes poorly understood. Here, we build on classic ecological theory by hypothesizing that some enigmatic spawning migrations across coastal marine habitats can be inferred from the population genetic signature of larval dispersal by ocean currents. We test this assumption by integrating spatially realistic simulations of alternative spawning migration routes, associated patterns of larval dispersal, and associated variation in the population genetic structure of eastern Australian sea mullet (Mugil cephalus). We then use simulation results to assess the implications of alternative spawning destinations for larval replenishment, and we contrast simulated against measured population genetic variation. Both analyses suggest that the spawning migrations of M. cephalus in eastern Australia are likely to be localized (approximately 100 km along the shore), and that spawning is likely to occur in inshore waters. Our conclusions are supported by multiple lines of evidence available through independent studies, but they challenge the more traditional assumption of a single, long-distance migration event with subsequent offshore spawning in the East Australian Current. More generally, our study operationalizes classic theory on the relationship between fish migrations, ocean currents, and reproductive success. However, rather than confirming the traditionally assumed adaptation of migratory behavior to dominant ocean current flow, our findings support the concept of a genetically measurable link between fish migrations and local oceanographic conditions, specifically water temperature and coastal retention of larvae. We believe that future studies using similar approaches for high resolution and spatially realistic ecological-genetic scenario testing can help rapidly advance our understanding of key ecological processes in many other marine species.

neogen: A tool to predict genetic effective population size (Ne ) for species with generational overlap and to assist empirical Ne study design.

Molecular genetic estimates of population effective size (Ne ) lose accuracy and precision when insufficient numbers of samples or loci are used. Ideally, researchers would like to forecast the necessary power when planning their project. neogen (genetic Ne for Overlapping Generations) enables estimates of precision and accuracy in advance of empirical investigation and allows exploration of the power available in different user-specified age-structured sampling schemes. neogen provides a population simulation and genetic power analysis framework that simulates the demographics, genetic composition, and Ne , from species-specific life history, mortality, population size, and genetic priors. neogen guides the user to establish a tractable sampling regime and to determine the numbers of samples and microsatellite or SNP loci required for accurate and precise genetic Ne estimates when sampling a natural population. neogen is useful at multiple stages of a study's life cycle: when budgeting, as sampling and locus development progresses, and for corroboration when empirical Ne estimates are available. The underlying model is applicable to a wide variety of iteroparous species with overlapping generations (e.g., mammals, birds, reptiles, long-lived fishes). In this paper, we describe the neogen model, detail the workflow for the point-and-click software, and explain the graphical results. We demonstrate the use of neogen with empirical Australian east coast zebra shark (Stegostoma fasciatum) data. For researchers wishing to make accurate and precise genetic Ne estimates for overlapping generations species, neogen facilitates planning for sample and locus acquisition, and with existing empirical genetic Ne estimates neogen can corroborate population demographic and life history properties.

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.

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.

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.

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.

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.

Strong population structure deduced from genetics, otolith chemistry and parasite abundances explains vulnerability to localized fishery collapse in a large Sciaenid fish, Protonibea diacanthus.

As pressure on coastal marine resources is increasing globally, the need to quantitatively assess vulnerable fish stocks is crucial in order to avoid the ecological consequences of stock depletions. Species of Sciaenidae (croakers, drums) are important components of tropical and temperate fisheries and are especially vulnerable to exploitation. The black-spotted croaker, Protonibea diacanthus, is the only large sciaenid in coastal waters of northern Australia where it is targeted by commercial, recreational and indigenous fishers due to its food value and predictable aggregating behaviour. Localized declines in the abundance of this species have been observed, highlighting the urgent requirement by managers for information on fine- and broad-scale population connectivity. This study examined the population structure of P. diacanthus across north-western Australia using three complementary methods: genetic variation in microsatellite markers, otolith elemental composition and parasite assemblage composition. The genetic analyses demonstrated that there were at least five genetically distinct populations across the study region, with gene flow most likely restricted by inshore biogeographic barriers such as the Dampier Peninsula. The otolith chemistry and parasite analyses also revealed strong spatial variation among locations within broad-scale regions, suggesting fine-scale location fidelity within the lifetimes of individual fish. The complementarity of the three techniques elucidated patterns of connectivity over a range of spatial and temporal scales. We conclude that fisheries stock assessments and management are required at fine scales (100 s of km) to account for the restricted exchange among populations (stocks) and to prevent localized extirpations of this species. Realistic management arrangements may involve the successive closure and opening of fishing areas to reduce fishing pressure.

Harnessing the Power of Genomics to Secure the Future of Seafood.

Best use of scientific knowledge is required to maintain the fundamental role of seafood in human nutrition. While it is acknowledged that genomic-based methods allow the collection of powerful data, their value to inform fisheries management, aquaculture, and biosecurity applications remains underestimated. We review genomic applications of relevance to the sustainable management of seafood resources, illustrate the benefits of, and identify barriers to their integration. We conclude that the value of genomic information towards securing the future of seafood does not need to be further demonstrated. Instead, we need immediate efforts to remove structural roadblocks and focus on ways that support integration of genomic-informed methods into management and production practices. We propose solutions to pave the way forward.

The complete validated mitochondrial genome of the silver gemfish Rexea solandri (Cuvier, 1832) (Perciformes, Gempylidae).

The silver gemfish Rexea solandri is an important economic resource but Vulnerable to overfishing in Australian waters. The complete mitochondrial genome sequence is described from 1.6 million reads obtained via next generation sequencing. The total length of the mitogenome is 16,350 bp comprising 2 rRNA, 13 protein-coding genes, 22 tRNA and 2 non-coding regions. The mitogenome sequence was validated against sequences of PCR fragments and BLAST queries of Genbank. Gene order was equivalent to that found in marine fishes.

The complete mitochondrial genome of the sandbar shark Carcharhinus plumbeus.

The sandbar shark, Carcharhinus plumbeus, a major representative species in shark fisheries worldwide is now considered vulnerable to overfishing. A pool of 774,234 Roche 454 shotgun sequences from one individual were assembled into a 16,706 bp mitogenome with 33× average coverage depth. It comprised 13 protein coding genes, 22 transfer RNA's, 2 ribosomal genes and 2 non-coding regions, typical of a vertebrate mitogenome. As expected for sharks, an A-T nucleotide bias was evident. This adds to rapidly growing number of mitogenome assemblies for the economically important Carcharhinidae family. The C. plumbeus mitogenome will assist researchers, fisheries and conservation managers interested in shark molecular systematics, phylogeography, conservation genetics, population and stock structure.

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 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.