Genome-wide SNPs resolve spatiotemporal patterns of connectivity within striped marlin (Kajikia audax), a broadly distributed and highly migratory pelagic species

Genomics and conservation: Guidance from training to analyses and applications

Environmental change is intensifying the biodiversity crisis and threatening species across the tree of life. Conservation genomics can help inform conservation actions and slow biodiversity loss. However, more training, appropriate use of novel genomic methods and communication with managers are needed. Here, we review practical guidance to improve applied conservation genomics. We share insights aimed at ensuring effectiveness of conservation actions around three themes: (1) improving pedagogy and training in conservation genomics including for online global audiences, (2) conducting rigorous population genomic analyses properly considering theory, marker types and data interpretation and (3) facilitating communication and collaboration between managers and researchers. We aim to update students and professionals and expand their conservation toolkit with genomic principles and recent approaches for conserving and managing biodiversity. The biodiversity crisis is a global problem and, as such, requires international involvement, training, collaboration and frequent reviews of the literature and workshops as we do here.

Unidirectional trans-Atlantic gene flow and a mixed spawning area shape the genetic connectivity of Atlantic bluefin tuna

The commercially important Atlantic bluefin tuna (Thunnus thynnus), a large migratory fish, has experienced notable recovery aided by accurate resource assessment and effective fisheries management efforts. Traditionally, this species has been perceived as consisting of eastern and western populations, spawning respectively in the Mediterranean Sea and the Gulf of Mexico, with mixing occurring throughout the Atlantic. However, recent studies have challenged this assumption by revealing weak genetic differentiation and identifying a previously unknown spawning ground in the Slope Sea used by Atlantic bluefin tuna of uncertain origin. To further understand the current and past population structure and connectivity of Atlantic bluefin tuna, we have assembled a unique dataset including thousands of genome-wide single-nucleotide polymorphisms (SNPs) from 500 larvae, young of the year and spawning adult samples covering the three spawning grounds and including individuals of other Thunnus species. Our analyses support two weakly differentiated but demographically connected ancestral populations that interbreed in the Slope Sea. Moreover, we also identified signatures of introgression from albacore (Thunnus alalunga) into the Atlantic bluefin tuna genome, exhibiting varied frequencies across spawning areas, indicating strong gene flow from the Mediterranean Sea towards the Slope Sea. We hypothesize that the observed genetic differentiation may be attributed to increased gene flow caused by a recent intensification of westward migration by the eastern population, which could have implications for the genetic diversity and conservation of western populations. Future conservation efforts should consider these findings to address potential genetic homogenization in the species.

Genomic assessment reveals signal of adaptive selection in populations of the Spotted rose snapper Lutjanus guttatus from the Tropical Eastern Pacific

Background

The lack of barriers in the marine environment has promoted the idea of panmixia in marine organisms. However, oceanographic conditions and habitat characteristics have recently been linked to genetic structure in marine species. The Tropical Eastern Pacific (TEP) is characterized by dynamic current systems and heterogeneous oceanographic conditions. The Gulf of Panama (part of the equatorial segment for the TEP) is influenced by a complex current system and heterogeneous environment, which has been shown to limit the gene flow for shoreline species. Next Generation Sequencing (NGS) has contributed to detect genetic differences in previously reported panmictic species by the assessment of loci associated with selection and to understand how selection acts affects marine populations. Lutjanus guttatus is a species distributed in the TEP for which previous studies using mitochondrial data recovered a panmictic pattern along its distributional range. In this study, we used SNP data of L. guttatus individuals sampled along its range to evaluate population genetic structure and investigate whether oceanographic factors influence the species’ genetic architecture. Finally, we assessed the role of adaptive selection by evaluating the contribution of outlier and neutral loci to genetic divergence.

Methods

The RADcap method was used to obtain 24 million paired reads for 123 individuals of L. guttatus covering nearly all its distributional area. Genetic variation was assessed using both spatial and non-spatial methods by comparing three different data sets: (i) a Combined Loci (CL dataset = 2003 SNPs); a search for putative loci under selection allowed the evaluation of (ii) Neutral Loci (NL dataset = 1858 SNPs) and (iii) Outlier Loci (OL dataset = 145 SNPs). We used the estimating effective migration surface (EEMS) approach to detect possible barriers to gene flow.

Results

Genetic differences were found in the OL dataset, showing two clusters (Northern and Southern), whereas NL showed no differences. This result may be related to the Selection-Migration balance model. The limit between the Northern and Southern groups was in the Gulf of Panama, which has been previously identified as a barrier to gene flow for other species, mainly due to its heterogeneous oceanographic conditions. The results suggest that selection plays an important role in generating genetic differences in Lutjanus guttatus. A migration corridor was detected that coincides with the Costa Rica Coastal Current that flows from Central America to the Gulf of California, allowing the homogenization of the northern population. In the Southern cluster, a migration corridor was observed with the OL from Panama to Colombia, which could be associated with the currents found in the Gulf of Panama. Genetic variation found in the OL of Lutjanus guttatus highlights the usefulness of NGS data in evaluating the role of selection in population differentiation.

Stepping up to genome scan allows stock differentiation in the worldwide distributed blue shark Prionace glauca

The blue shark Prionace glauca is a top predator with one of the widest geographical distributions of any shark species. It is classified as Critically Endangered in the Mediterranean Sea, and Near Threatened globally. Previous genetic studies did not reject the null hypothesis of a single global population. The blue shark was proposed as a possible archetype of the "grey zone of population differentiation," coined to designate cases where population structure may be too recent or too faint to be detected using a limited set of markers. Here, blue shark samples collected throughout its global range were sequenced using a specific RAD method (DArTseq), which recovered 37,655 genome-wide single nucleotide polymorphisms (SNPs). Two main groups emerged, with Mediterranean Sea and northern Atlantic samples (Northern population) differentiated significantly from the Indo-west Pacific samples (Southern population). Significant pairwise F_ST values indicated further genetic differentiation within the Atlantic Ocean, and between the Atlantic Ocean and the Mediterranean Sea. Reconstruction of recent demographic history suggested divergence between Northern and Southern populations occurred about 500 generations ago and revealed a drastic reduction in effective population size from a large ancestral population. Our results illustrate the power of genome scans to detect population structure and reconstruct demographic history in highly migratory marine species. Given that the management plans of the blue shark (targeted or bycatch) fisheries currently assume panmictic regional stocks, we strongly recommend that the results presented here be considered in future stock assessments and conservation strategies.

Are genomic updates of well-studied species worth the investment for conservation? A case study of the Critically Endangered Magdalena River turtle

As genomic-scale data sets become economically feasible for most organisms, a key question for conservation biology is whether the increased resolution offered by new genomic approaches justifies repeating earlier studies based on traditional markers, rather than investing those same time and monetary resources in less-known species. Genomic studies offer clear advantages when the objective is to identify adaptive loci that may be critical to conservation policy-makers. However, the answer is far less certain for the population and landscape studies based on neutral loci that dominate the conservation genetics research agenda. We used RADseq to revisit earlier molecular studies of the IUCN Critically Endangered Magdalena River turtle (Podocnemis lewyana), documenting the conservation insights gained by increasing the number of neutral markers by several orders of magnitude. Earlier research indicated that P. lewyana has the lowest genetic diversity known for any chelonian, and little or no population differentiation among independent rivers. In contrast, the RADseq data revealed discrete population structure with isolation-by-distance within river segments and identified precise population breaks clearly delineating management units. It also confirmed that the species does not have extremely low heterozygosity and that effective population sizes are probably sufficient to maintain long-term evolutionary potential. Contrary to earlier inferences from more limited population genetic markers, our genomic data suggest that management strategies should shift from active genetic rescue to more passive protection without extreme interventions. We conclude with a list of examples of conservation studies in other vertebrates indicating that for many systems a genomic update is worth the investment.

Genomic assessment of global population structure in a highly migratory and habitat versatile apex predator, the tiger shark (Galeocerdo cuvier)

Understanding the population dynamics of highly mobile, widely distributed, oceanic sharks, many of which are overexploited, is necessary to aid their conservation management. We investigated the global population genomics of tiger sharks (Galeocerdo cuvier), a circumglobally distributed, apex predator displaying remarkable behavioral versatility in its diet, habitat use (near coastal, coral reef, pelagic), and individual movement patterns (spatially resident to long-distance migrations). We genotyped 242 tiger sharks from 10 globally distributed locations at more than 2000 single nucleotide polymorphisms. Although this species often conducts massive distance migrations, the data show strong genetic differentiation at both neutral (FST=0.125-0.144) and candidate outlier loci (FST=0.570-0.761) between western Atlantic and Indo-Pacific sharks, suggesting the potential for adaptation to the environments specific to these oceanic regions. Within these regions, there was mixed support for population differentiation between northern and southern hemispheres in the western Atlantic, and none for structure within the Indian Ocean. Notably, the results demonstrate a low level of population differentiation of tiger sharks from the remote Hawaiian archipelago compared to sharks from the Indian Ocean (FST=0.003-0.005, P<0.01). Given concerns about biodiversity loss and marine ecosystem impacts caused by overfishing of oceanic sharks in the midst of rapid environmental change, our results suggest it imperative that international fishery management prioritize conservation of the evolutionary potential of the highly genetically differentiated Atlantic and Indo-Pacific populations of this unique apex predator. Furthermore, we suggest targeted management attention to tiger sharks in the Hawaiian archipelago based on a precautionary biodiversity conservation perspective.

Adaptive markers distinguish North and South Pacific Albacore amid low population differentiation

Albacore (Thunnus alalunga) support an economically valuable global fishery, but surprisingly little is known about the population structure of this highly migratory species. Physical tagging data suggest that Albacore from the North and South Pacific Ocean are separate stocks, but results from previous genetic studies did not support this two stock hypothesis. In addition, observed biological differences among juveniles suggest that there may be population substructure in the North Pacific. We used double-digest restriction site-associated DNA sequencing to assess population structure among 308 Albacore caught in 12 sample areas across the Pacific Ocean (10 North, 2 South). Since Albacore are highly migratory and spawning areas are unknown, sample groups were not assumed to be equivalent to populations and the genetic data were analyzed iteratively. We tested for putatively adaptive differences among groups and for genetic variation associated with sex. Results indicated that Albacore in the North and South Pacific can be distinguished using 84 putatively adaptive loci, but not using the remaining 12,788 presumed neutral sites. However, two individuals likely represent F1 hybrids between the North and South Pacific populations, and 43 Albacore potentially exhibit lower degrees of mixed ancestry. In addition, four or five cross-hemisphere migrants were potentially identified. No genetic evidence was found for population substructure within the North Pacific, and no loci appeared to distinguish males from females. Potential functions for the putatively adaptive loci were identified, but an annotated Albacore genome is required for further exploration. Future research should try to locate spawning areas so that life history, demography, and genetic population structure can be linked and spatiotemporal patterns can be investigated.

Genomic signatures of divergent selection are associated with social behaviour for spinner dolphin ecotypes

Understanding the genomic basis of adaptation is critical for understanding evolutionary processes and predicting how species will respond to environmental change. Spinner dolphins in the eastern tropical Pacific (ETP) present a unique system for studying adaptation. Within this large geographical region are four spinner dolphin ecotypes with weak neutral genetic divergence and no obvious barriers to gene flow, but strong spatial variation in morphology, behaviour and habitat. These ecotypes have large population sizes, which could reduce the effects of drift and facilitate selection. To identify genomic regions putatively under divergent selective pressures between ecotypes, we used genome scans with 8994 RADseq single nucleotide polymorphisms (SNPs) to identify population differentiation outliers and genotype-environment association outliers. Gene ontology enrichment analyses indicated that outlier SNPs from both types of analyses were associated with multiple genes involved in social behaviour and hippocampus development, including 15 genes associated with the human social disorder autism. Evidence for divergent selection on social behaviour is supported by previous evidence that these spinner dolphin ecotypes differ in mating systems and associated social behaviours. In particular, three of the ETP ecotypes probably have a polygynous mating system characterized by strong premating competition among males, whereas the fourth ecotype probably has a polygynandrous mating system characterized by strong postmating competition such as sperm competition. Our results provide evidence that selection for social behaviour may be an evolutionary force driving diversification of spinner dolphins in the ETP, potentially as a result of divergent sexual selection associated with different mating systems. Future studies should further investigate the potential adaptive role of the candidate genes identified here, and could probably find further signatures of selection using whole genome sequence data.

Current status and potential contributions of fisheries statistics from artisanal fisheries for managing juvenile istiophorid billfishes in Southern Brazil

Billfishes are considered important fishery resources and the identification of aggregation sites is imperative for proper management. Here we present evidence of a seasonal aggregation site for juvenile istiophorid billfishes in southern Brazil. We discuss the results as they relate with the need for management at a local scale and participatory monitoring with artisanal fishing communities as a way to ensure access to data on the occurrence and population status in the long term.