Genetic stock structure of New Zealand fish and the use of genomics in fisheries management: an overview and outlook

Peering into the gaps: Long-read sequencing illuminates structural variants and genomic evolution in the Australasian snapper

Even before genome sequencing, genetic resources have supported species management and breeding programs. Current technologies, such as long-read sequencing, resolve complex genomic regions, like those rich in repeats or high in GC content. Improved genome contiguity enhances accuracy in identifying structural variants (SVs) and transposable elements (TEs). We present an improved genome assembly and SV catalogue for the Australasian snapper (Chrysophrys auratus). The new assembly is more contiguous, allowing for putative identification of 14 centromeres and transfer of 26,115 gene annotations from yellowfin seabream. Compared to the previous assembly, 35,000 additional SVs, including larger and more complex rearrangements, were annotated. SVs and TEs exhibit a distribution pattern skewed towards chromosome ends, likely influenced by recombination. Some SVs overlap with growth-related genes, underscoring their significance. This upgraded genome serves as a foundation for studying natural and artificial selection, offers a reference for related species, and sheds light on genome dynamics shaped by evolution.

Genetic Structure and Diversity of Hatchery and Wild Populations of Yellow Catfish Tachysurus fulvidraco (Siluriformes: Bagridae) from Korea

Yellow catfish Tachysurus fulvidraco is an important commercial fish species in South Korea. However, due to their current declines in its distribution area and population size, it is being released from hatchery populations into wild populations. Hatchery populations also produced from wild broodstocks are used for its captive breeding. We reported 15 new microsatellite DNA markers of T. fulvidraco to identify the genetic diversity and structure of its hatchery and wild populations, providing baseline data for useful resource development strategies. The observed heterozygosity of the hatchery populations ranged from 0.816 to 0.873, and that of the wild populations ranged from 0.771 to 0.840. Their inbreeding coefficient ranged from -0.078 to 0.024. All populations experienced a bottleneck (p < 0.05), with effective population sizes ranging from 21 to infinity. Their gene structure was divided into two groups with STRUCTURE results of K = 2. It was confirmed that each hatchery population originated from a different wild population. This study provides genetic information necessary for the future development and conservation of fishery resources for T. fulvidraco.

Genomic analyses indicate resilience of a commercially and culturally important marine gastropod snail to climate change

Genomic vulnerability analyses are being increasingly used to assess the adaptability of species to climate change and provide an opportunity for proactive management of harvested marine species in changing oceans. Southeastern Australia is a climate change hotspot where many marine species are shifting poleward. The turban snail, Turbo militaris is a commercially and culturally harvested marine gastropod snail from eastern Australia. The species has exhibited a climate-driven poleward range shift over the last two decades presenting an ongoing challenge for sustainable fisheries management. We investigate the impact of future climate change on T. militaris using genotype-by-sequencing to project patterns of gene flow and local adaptation across its range under climate change scenarios. A single admixed, and potentially panmictic, demographic unit was revealed with no evidence of genetic subdivision across the species range. Significant genotype associations with heterogeneous habitat features were observed, including associations with sea surface temperature, ocean currents, and nutrients, indicating possible adaptive genetic differentiation. These findings suggest that standing genetic variation may be available for selection to counter future environmental change, assisted by widespread gene flow, high fecundity and short generation time in this species. We discuss the findings of this study in the content of future fisheries management and conservation.

Estimating microhaplotype allele frequencies from low-coverage or pooled sequencing data

BACKGROUND

Microhaplotypes have the potential to be more cost-effective than SNPs for applications that require genetic panels of highly variable loci. However, development of microhaplotype panels is hindered by a lack of methods for estimating microhaplotype allele frequency from low-coverage whole genome sequencing or pooled sequencing (pool-seq) data.

RESULTS

We developed new methods for estimating microhaplotype allele frequency from low-coverage whole genome sequence and pool-seq data. We validated these methods using datasets from three non-model organisms. These methods allowed estimation of allele frequency and expected heterozygosity at depths routinely achieved from pooled sequencing.

CONCLUSIONS

These new methods will allow microhaplotype panels to be designed using low-coverage WGS and pool-seq data to discover and evaluate candidate loci. The python script implementing the two methods and documentation are available at https://www.github.com/delomast/mhFromLowDepSeq .

A multiplexed plant-animal SNP array for selective breeding and species conservation applications

Reliable and high-throughput genotyping platforms are of immense importance for identifying and dissecting genomic regions controlling important phenotypes, supporting selection processes in breeding programs, and managing wild populations and germplasm collections. Amongst available genotyping tools, single nucleotide polymorphism arrays have been shown to be comparatively easy to use and generate highly accurate genotypic data. Single-species arrays are the most commonly used type so far; however, some multi-species arrays have been developed for closely related species that share single nucleotide polymorphism markers, exploiting inter-species cross-amplification. In this study, the suitability of a multiplexed plant-animal single nucleotide polymorphism array, including both closely and distantly related species, was explored. The performance of the single nucleotide polymorphism array across species for diverse applications, ranging from intra-species diversity assessments to parentage analysis, was assessed. Moreover, the value of genotyping pooled DNA of distantly related species on the single nucleotide polymorphism array as a technique to further reduce costs was evaluated. Single nucleotide polymorphism performance was generally high, and species-specific single nucleotide polymorphisms proved suitable for diverse applications. The multi-species single nucleotide polymorphism array approach reported here could be transferred to other species to achieve cost savings resulting from the increased throughput when several projects use the same array, and the pooling technique adds another highly promising advancement to additionally decrease genotyping costs by half.

Oceans of plenty? Challenges, advancements, and future directions for the provision of evidence-based fisheries management advice

Marine population modeling, which underpins the scientific advice to support fisheries interventions, is an active research field with recent advancements to address modern challenges (e.g., climate change) and enduring issues (e.g., data limitations). Based on discussions during the 'Land of Plenty' session at the 2021 World Fisheries Congress, we synthesize current challenges, recent advances, and interdisciplinary developments in biological fisheries models (i.e., data-limited, stock assessment, spatial, ecosystem, and climate), management strategy evaluation, and the scientific advice that bridges the science-policy interface. Our review demonstrates that proliferation of interdisciplinary research teams and enhanced data collection protocols have enabled increased integration of spatiotemporal, ecosystem, and socioeconomic dimensions in many fisheries models. However, not all management systems have the resources to implement model-based advice, while protocols for sharing confidential data are lacking and impeding research advances. We recommend that management and modeling frameworks continue to adopt participatory co-management approaches that emphasize wider inclusion of local knowledge and stakeholder input to fill knowledge gaps and promote information sharing. Moreover, fisheries management, by which we mean the end-to-end process of data collection, scientific analysis, and implementation of evidence-informed management actions, must integrate improved communication, engagement, and capacity building, while incorporating feedback loops at each stage. Increasing application of management strategy evaluation is viewed as a critical unifying component, which will bridge fisheries modeling disciplines, aid management decision-making, and better incorporate the array of stakeholders, thereby leading to a more proactive, pragmatic, transparent, and inclusive management framework-ensuring better informed decisions in an uncertain world.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s11160-022-09726-7.

Genome assembly and isoform analysis of a highly heterozygous New Zealand fisheries species, the tarakihi (Nemadactylus macropterus)

Although being some of the most valuable and heavily exploited wild organisms, few fisheries species have been studied at the whole-genome level. This is especially the case in New Zealand, where genomics resources are urgently needed to assist fisheries management. Here, we generated 55 Gb of short Illumina reads (92× coverage) and 73 Gb of long Nanopore reads (122×) to produce the first genome assembly of the marine teleost tarakihi [Nemadactylus macropterus (Forster, 1801)], a highly valuable fisheries species in New Zealand. An additional 300 Mb of Iso-Seq reads were obtained to assist in gene annotation. The final genome assembly was 568 Mb long with an N50 of 3.37 Mb. The genome completeness was high, with 97.8% of complete Actinopterygii Benchmarking Universal Single-Copy Orthologs. Heterozygosity values estimated through k-mer counting (1.00%) and bi-allelic SNPs (0.64%) were high compared with the same values reported for other fishes. Iso-Seq analysis recovered 91,313 unique transcripts from 15,515 genes (mean ratio of 5.89 transcripts per gene), and the most common alternative splicing event was intron retention. This highly contiguous genome assembly and the isoform-resolved transcriptome will provide a useful resource to assist the study of population genomics and comparative eco-evolutionary studies in teleosts and related organisms.

Salmonidae Genome: Features, Evolutionary and Phylogenetic Characteristics

The salmon family is one of the most iconic and economically important fish families, primarily possessing meat of excellent taste as well as irreplaceable nutritional and biological value. One of the most common and, therefore, highly significant members of this family, the Atlantic salmon (Salmo salar L.), was not without reason one of the first fish species for which a high-quality reference genome assembly was produced and published. Genomic advancements are becoming increasingly essential in both the genetic enhancement of farmed salmon and the conservation of wild salmon stocks. The salmon genome has also played a significant role in influencing our comprehension of the evolutionary and functional ramifications of the ancestral whole-genome duplication event shared by all Salmonidae species. Here we provide an overview of the current state of research on the genomics and phylogeny of the various most studied subfamilies, genera, and individual salmonid species, focusing on those studies that aim to advance our understanding of salmonid ecology, physiology, and evolution, particularly for the purpose of improving aquaculture production. This review should make potential researchers pay attention to the current state of research on the salmonid genome, which should potentially attract interest in this important problem, and hence the application of new technologies (such as genome editing) in uncovering the genetic and evolutionary features of salmoniforms that underlie functional variation in traits of commercial and scientific importance.

Genomic Stock Structure of the Marine Teleost Tarakihi (Nemadactylus macropterus) Provides Evidence of Potential Fine-Scale Adaptation and a Temperature-Associated Cline Amid Panmixia

Tarakihi (Nemadactylus macropterus) is an important fishery species with widespread distribution around New Zealand and off the southern coasts of Australia. However, little is known about whether the populations are locally adapted or genetically structured. To address this, we conducted whole-genome resequencing of 175 tarakihi from around New Zealand and Tasmania (Australia) to obtain a dataset of 7.5 million genome-wide and high-quality single nucleotide polymorphisms (SNPs). Variant filtering, FST-outlier analysis, and redundancy analysis (RDA) were used to evaluate population structure, adaptive structure, and locus-environment associations. A weak but significant level of neutral genetic differentiation was found between tarakihi from New Zealand and Tasmania (FST = 0.0054–0.0073, P ≤ 0.05), supporting the existence of at least two separate reproductive stocks. No clustering was detected among the New Zealand populations (ΦST &lt; 0.001, P = 0.77). Outlier-based, presumably adaptive variation suggests fine-scale adaptive structure between locations around central New Zealand off the east (Wairarapa, Cape Campbell, and Hawke’s Bay) and the west coast (Tasman Bay/Golden Bay and Upper West Coast of South Island). Allele frequencies from 55 loci were associated with at least one of six environmental variables, of which 47 correlated strongly with yearly mean water temperature. Although genes associated with these loci are linked to various functions, the most common functions were integral components of membrane and cilium assembly. Projection of the RDA indicates the existence of a latitudinal temperature cline. Our work provides the first genomic insights supporting panmixia of tarakihi in New Zealand and evidence of a genomic cline that appears to be driven by the temperature gradients, together providing crucial information to inform the stock assessment of this species, and to widen the insights of the ecological drivers of adaptive variation in a marine species.

Genomic prediction of growth in a commercially, recreationally, and culturally important marine resource, the Australian snapper (Chrysophrys auratus)

Growth is one of the most important traits of an organism. For exploited species, this trait has ecological and evolutionary consequences as well as economical and conservation significance. Rapid changes in growth rate associated with anthropogenic stressors have been reported for several marine fishes, but little is known about the genetic basis of growth traits in teleosts. We used reduced genome representation data and genome-wide association approaches to identify growth-related genetic variation in the commercially, recreationally, and culturally important Australian snapper (Chrysophrys auratus, Sparidae). Based on 17,490 high-quality single-nucleotide polymorphisms and 363 individuals representing extreme growth phenotypes from 15,000 fish of the same age and reared under identical conditions in a sea pen, we identified 100 unique candidates that were annotated to 51 proteins. We documented a complex polygenic nature of growth in the species that included several loci with small effects and a few loci with larger effects. Overall heritability was high (75.7%), reflected in the high accuracy of the genomic prediction for the phenotype (small vs large). Although the single-nucleotide polymorphisms were distributed across the genome, most candidates (60%) clustered on chromosome 16, which also explains the largest proportion of heritability (16.4%). This study demonstrates that reduced genome representation single-nucleotide polymorphisms and the right bioinformatic tools provide a cost-efficient approach to identify growth-related loci and to describe genomic architectures of complex quantitative traits. Our results help to inform captive aquaculture breeding programs and are of relevance to monitor growth-related evolutionary shifts in wild populations in response to anthropogenic pressures.

Population structure of lake whitefish (Coregonus clupeaformis) from the Mississippian lineage in North America

The lake whitefish ( Coregonus clupeaformis) is a commercially valuable freshwater species with a broad distribution in North America. Some phylogeographic work has been done on this species, but little is known about genetic population subdivision among populations of the widely dispersed Mississippian lineage. We used 3,173 single nucleotide polymorphisms in 508 lake whitefish from 22 different lakes to examine population structure across central Canada and the United States. Bayesian clustering, ordination, and fixation indices identified population subdivision that largely reflected geographic distance and hydrological connectivity, with greater differentiation between lakes that are farther apart. Population subdivision was hierarchical, with greater differentiation between Canadian provinces and less differentiation based on river basins within provincial boundaries. Interestingly, isolation by distance alone was not sufficient to account for all of the observed genetic differentiation among populations. We conclude that important components of lake whitefish genetic diversity are present at different spatial scales, and that populations within the Mississippian lineage have differentiated widely across their range.

Genome-wide analysis reveals the genetic stock structure of hoki (Macruronus novaezelandiae)

The assessment of the genetic structuring of biodiversity is crucial for management and conservation. This is particularly critical for widely distributed and highly mobile deep-water teleosts, such as hoki (Macruronus novaezelandiae). This species is significant to Māori people and supports the largest commercial fishery in New Zealand, but uncertainty about its stock structure presents a challenge for management. Here, we apply a comprehensive genomic analysis to shed light on the demographic structure of this species by (1) assembling the genome, (2) generating a catalogue of genome-wide SNPs to infer the stock structure and (3) identifying regions of the genome under selection. The final genome assembly used short and long reads and is near complete, representing 93.8% of BUSCO genes, and consisting of 566 contigs totalling 501 Mb. Whole-genome re-sequencing of 510 hoki sampled from 14 locations around New Zealand and Australia, at a read depth greater than 10×, produced 227,490 filtered SNPs. Analyses of these SNPs were able to resolve the stock structure of hoki into two genetically and geographically distinct clusters, one including the Australian and the other one all New Zealand locations, indicating genetic exchange between these regions is limited. Location differences within New Zealand samples were much more subtle (global F ST  = 0.0006), and while small and significant differences could be detected, they did not conclusively identify additional substructures. Ten putative adaptive SNPs were detected within the New Zealand samples, but these also did not geographically partition the dataset further. Contemporary and historical N e estimation suggest the current New Zealand population of hoki is large yet declining. Overall, our study provides the first genomic resources for hoki and provides detailed insights into the fine-scale population genetic structure to inform the management of this species.

The genome of New Zealand trevally (Carangidae: Pseudocaranx georgianus) uncovers a XY sex determination locus

BACKGROUND

The genetic control of sex determination in teleost species is poorly understood. This is partly because of the diversity of mechanisms that determine sex in this large group of vertebrates, including constitutive genes linked to sex chromosomes, polygenic constitutive mechanisms, environmental factors, hermaphroditism, and unisexuality. Here we use a de novo genome assembly of New Zealand silver trevally (Pseudocaranx georgianus) together with sex-specific whole genome sequencing data to detect sexually divergent genomic regions, identify candidate genes and develop molecular makers.

RESULTS

The de novo assembly of an unsexed trevally (Trevally_v1) resulted in a final assembly of 579.4 Mb in length, with a N50 of 25.2 Mb. Of the assembled scaffolds, 24 were of chromosome scale, ranging from 11 to 31 Mb in length. A total of 28,416 genes were annotated after 12.8 % of the assembly was masked with repetitive elements. Whole genome re-sequencing of 13 wild sexed trevally (seven males and six females) identified two sexually divergent regions located on two scaffolds, including a 6 kb region at the proximal end of chromosome 21. Blast analyses revealed similarity between one region and the aromatase genes cyp19 (a1a/b) (E-value < 1.00E-25, identity > 78.8 %). Males contained higher numbers of heterozygous variants in both regions, while females showed regions of very low read-depth, indicative of male-specificity of this genomic region. Molecular markers were developed and subsequently tested on 96 histologically-sexed fish (42 males and 54 females). Three markers amplified in absolute correspondence with sex (positive in males, negative in females).

CONCLUSIONS

The higher number of heterozygous variants in males combined with the absence of these regions in females support a XY sex-determination model, indicating that the trevally_v1 genome assembly was developed from a male specimen. This sex system contrasts with the ZW sex-determination model documented in closely related carangid species. Our results indicate a sex-determining function of a cyp19a1a-like gene, suggesting the molecular pathway of sex determination is somewhat conserved in this family. The genomic resources developed here will facilitate future comparative work, and enable improved insights into the varied sex determination pathways in teleosts. The sex marker developed in this study will be a valuable resource for aquaculture selective breeding programmes, and for determining sex ratios in wild populations.

Charting a course for genetic diversity in the UN Decade of Ocean Science

The health of the world's oceans is intrinsically linked to the biodiversity of the ecosystems they sustain. The importance of protecting and maintaining ocean biodiversity has been affirmed through the setting of the UN Sustainable Development Goal 14 to conserve and sustainably use the ocean for society's continuing needs. The decade beginning 2021-2030 has additionally been declared as the UN Decade of Ocean Science for Sustainable Development. This program aims to maximize the benefits of ocean science to the management, conservation, and sustainable development of the marine environment by facilitating communication and cooperation at the science-policy interface. A central principle of the program is the conservation of species and ecosystem components of biodiversity. However, a significant omission from the draft version of the Decade of Ocean Science Implementation Plan is the acknowledgment of the importance of monitoring and maintaining genetic biodiversity within species. In this paper, we emphasize the importance of genetic diversity to adaptive capacity, evolutionary potential, community function, and resilience within populations, as well as highlighting some of the major threats to genetic diversity in the marine environment from direct human impacts and the effects of global climate change. We then highlight the significance of ocean genetic diversity to a diverse range of socioeconomic factors in the marine environment, including marine industries, welfare and leisure pursuits, coastal communities, and wider society. Genetic biodiversity in the ocean, and its monitoring and maintenance, is then discussed with respect to its integral role in the successful realization of the 2030 vision for the Decade of Ocean Science. Finally, we suggest how ocean genetic diversity might be better integrated into biodiversity management practices through the continued interaction between environmental managers and scientists, as well as through key leverage points in industry requirements for Blue Capital financing and social responsibility.

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.