Evidence of discrete yellowfin tuna (Thunnus albacares) populations demands rethink of management for this globally important resource

Estimating Demographic Parameters for Bearded Seals, Erignathus barbatus, in Alaska Using Close-Kin Mark-Recapture Methods

Reliable estimates of population abundance and demographics are essential for managing harvested species. Ice-associated phocids, "ice seals," are a vital resource for subsistence-dependent coastal Native communities in western and northern Alaska, USA. In 2012, the Beringia distinct population segment of the bearded seal, Erignathus barbatus nauticus, was listed as "threatened" under the US Endangered Species Act requiring greater scrutiny for management assessments. We sought to estimate requisite population parameters from harvested seals by using close-kin mark-recapture (CKMR) methods, the first such application for marine mammals. Samples from 1758 bearded seals harvested by Bering, Chukchi, and Beaufort Sea communities during 1998-2020 were genotyped, genetically sexed, and aged by tooth annuli. After rigorous quality control, kin relationships were established for 1484 seals including two parent-offspring pairs (POPs) and 25 potential second-order kin pairs. Most of the second-order kin were half-sibling pairs (HSPs), but four were potential grandparent-grandchild pairs (GGPs). There were no full sibling pairs, suggesting a lack of mate fidelity. Mitochondrial DNA analysis identified 17 potential HSPs as paternally related, providing substantial evidence of persistent heterogeneity in reproductive success among adult males. The statistical CKMR model incorporates probabilities associated with POPs, HSPs, and GGPs and assumes known ages and a stable population. Our top model accommodates heterogeneity in adult male breeding success and yields an abundance estimate of ~409,000 with a coefficient of variation (CV) = 0.35, which is substantially greater than the "non-heterogeneity" model estimate of ~232,000 (CV = 0.21), an important difference for managing a harvested species. Using CKMR methods with harvested species provides estimates of abundance with the added opportunity to acquire information about adult survival, fecundity, and breeding success that could be applied to other species of concern, marine and terrestrial.

Implications of past and present genetic connectivity for management of the saltwater crocodile (Crocodylus porosus)

Effective management of protected species requires information on appropriate evolutionary and geographic population boundaries and knowledge of how the physical environment and life-history traits combine to shape the population structure and connectivity. Saltwater crocodiles (Crocodylus porosus) are the largest and most widely distributed of living crocodilians, extending from Sri Lanka to Southeast Asia and down to northern Australia. Given the long-distance movement capabilities reported for C. porosus, management units are hypothesised to be highly connected by migration. However, the magnitude, scale, and consistency of connection across managed populations are not fully understood. Here we used an efficient genotyping method that combines DArTseq and sequence capture to survey≈ 3000 high-quality genome-wide single nucleotide polymorphisms from 1176 C. porosus sampled across nearly the entire range of the species in Queensland, Australia. We investigated historical and present-day connectivity patterns using fixation and diversity indices coupled with clustering methods and the spatial distribution of kin pairs. We inferred kinship using forward simulation coupled with a kinship estimation method that is robust to unspecified population structure. The results demonstrated that the C. porosus population has substantial genetic structure with six broad populations correlated with geographical location. The rate of gene flow was highly correlated with spatial distance, with greater differentiation along the east coast compared to the west. Kinship analyses revealed evidence of reproductive philopatry and limited dispersal, with approximately 90% of reported first and second-degree relatives showing a pairwise distance of <50 km between sampling locations. Given the limited dispersal, lack of suitable habitat, low densities of crocodiles and the high proportion of immature animals in the population, future management and conservation interventions should be considered at regional and state-wide scales.

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.

Lack of genetic differentiation in yellowfin tuna has conservation implications in the Eastern Pacific Ocean

Yellowfin tuna, Thunnus albacares, is an important global fishery and of particular importance in the Eastern Pacific Ocean (EPO). According to the 2019 Inter-American Tropical Tuna Commission (IATTC) assessment, yellowfin tuna within the EPO is a single stock, and is being managed as one stock. However, previous studies indicate site fidelity, or limited home ranges, of yellowfin tuna which suggests the potential for multiple yellowfin tuna stocks within the EPO, which was supported by a population genetic study using microsatellites. If numerous stocks are present, management at the wrong spatial scales could cause the loss of minor yellowfin tuna populations in the EPO. In this study we used double digestion RADseq to assess the genetic structure of yellowfin tuna in the EPO. A total of 164 yellowfin tuna from Cabo San Lucas, México, and the Galápagos Islands and Santa Elena, Ecuador, were analysed using 18,011 single nucleotide polymorphisms. Limited genetic differentiation (F_ST = 0.00058–0.00328) observed among the sampling locations (México, Ecuador, Peru, and within Ecuador) is consistent with presence of a single yellowfin tuna population within the EPO. Our findings are consistent with the IATTC assessment and provide further evidence of the need for transboundary cooperation for the successful management of this important fishery throughout the EPO.

Sample size requirements for genetic studies on yellowfin tuna

In population genetics, the amount of information for an analytical task is governed by the number of individuals sampled and the amount of genetic information measured on each of those individuals. In this work, we assessed the numbers of individual yellowfin tuna (Thunnus albacares) and genetic markers required for ocean-basin scale inferences. We assessed this for three distinct data analysis tasks that are often employed: testing for differences between genetic profiles; stock delineation, and; assignment of individuals to stocks. For all analytical tasks, we used real (not simulated) data from four sampling locations that span the tropical Pacific Ocean. Whilst spatially separated, the genetic differences between the sampling sites were not substantial, a maximum of approximately F_st = 0.02, which is quite typical of large pelagic fish. We repeatedly sub-sampled the data, mimicking a new survey, and performed the analyses. False positive rates were also assessed by re-sampling and randomly assigning fish to groups. Varying the sample sizes indicated that some analytical tasks, namely profile testing, required relatively few individuals per sampling location (n ≳ 10) and single nucleotide polymorphisms (SNPs, m ≳ 256). Stock delineation required more individuals per sampling location (n ≳ 25). Assignment of fish to sampling locations required substantially more individuals, more in fact than we had available (n > 50), although this sample size could be reduced to n ≳ 30 when individual fish were assumed to belong to one of the groups sampled. With these results, designers of molecular ecological surveys for yellowfin tuna, and users of information from them, can assess whether the information content is adequate for the required inferential task.

Climate-assisted persistence of tropical fish vagrants in temperate marine ecosystems

Rising temperatures and extreme climate events are propelling tropical species into temperate marine ecosystems, but not all species can persist. Here, we used the heatwave-driven expatriation of tropical Black Rabbitfish (Siganus fuscescens) to the temperate environments of Western Australia to assess the ecological and evolutionary mechanisms that may entail their persistence. Population genomic assays for this rabbitfish indicated little genetic differentiation between tropical residents and vagrants to temperate environments due to high migration rates, which were likely enhanced by the marine heatwave. DNA metabarcoding revealed a diverse diet for this species based on phytoplankton and algae, as well as an ability to feed on regional resources, including kelp. Irrespective of future climate scenarios, these macroalgae-consuming vagrants may self-recruit in temperate environments and further expand their geographic range by the year 2100. This expansion may compromise the health of the kelp forests that form Australia's Great Southern Reef. Overall, our study demonstrates that projected favourable climate conditions, continued large-scale genetic connectivity between populations, and diet versatility are key for tropical range-shifting fish to establish in temperate ecosystems.

Genomic methods reveal independent demographic histories despite strong morphological conservatism in fish species

Human overexploitation of natural resources has placed conservation and management as one of the most pressing challenges in modern societies, especially in regards to highly vulnerable marine ecosystems. In this context, cryptic species are particularly challenging to conserve because they are hard to distinguish based on morphology alone, and thus it is often unclear how many species coexist in sympatry, what are their phylogenetic relationships and their demographic history. We answer these questions using morphologically similar species of the genus Mugil that are sympatric in the largest coastal Marine Protected Area in the Tropical Southwestern Atlantic marine province. Using a sub-representation of the genome, we show that individuals are assigned to five highly differentiated genetic clusters that are coincident with five mitochondrial lineages, but discordant with morphological information, supporting the existence of five species with conserved morphology in this region. A lack of admixed individuals is consistent with strong genetic isolation between sympatric species, but the most likely species tree suggests that in one case speciation has occurred in the presence of interspecific gene flow. Patterns of genetic diversity within species suggest that effective population sizes differ up to two-fold, probably reflecting differences in the magnitude of population expansions since species formation. Together, our results show that strong morphologic conservatism in marine environments can lead to species that are difficult to distinguish morphologically but that are characterized by an independent evolutionary history, and thus that deserve species-specific management strategies.

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.

The effects of weather variability on patterns of genetic diversity in Tasmanian bettongs

While the effects of climate (long-term, prevailing weather) on species abundance, range and genetic diversity have been widely studied, short-term, localized variations in atmospheric conditions (i.e., weather) can also rapidly alter species' geographical ranges and population sizes, but little is known about how they affect genetic diversity. We investigated the relationship between weather and range-wide genetic diversity in a marsupial, Bettongia gaimardi, using dynamic species distribution models (SDMs). Genetic diversity was lower in parts of the range where the weather-based SDM predicted high variability in probability of B. gaimardi occurrence during 1950-2009. This is probably an effect of lower population sizes and extinction-recolonization cycles in places with highly variable weather. Spatial variation in genetic diversity was also better predicted by mean probabilities of B. gaimardi occurrence from weather- than climate-based SDMs. Our results illustrate the importance of weather in driving population dynamics and species distributions on decadal timescales and thereby in affecting genetic diversity. Modelling the links between changing weather patterns, species distributions and genetic diversity will allow researchers to better forecast biological impacts of climate change.

Distinguishing Between Nile Tilapia Strains Using a Low-Density Single-Nucleotide Polymorphism Panel

Nile tilapia (Oreochromis niloticus) is among the most important finfish in aquaculture, particularly in Asia. Numerous genetically improved strains of Nile tilapia have been developed and disseminated through formal and informal channels to hatcheries, many of which operate at a relatively small scale in developing countries. The primary objective of this study was to assess the extent to which molecular genetic tools can identify different and interrelated strains of Nile tilapia in Bangladesh and the Philippines, two globally significant producers. A tool was developed using a low-density panel of single-nucleotide polymorphisms (SNPs), genotyping-by-sequencing and discriminant analysis of principal components (DAPC). When applied to 2,057 samples from 205 hatcheries in Bangladesh and the Philippines, for hatcheries where the hatchery-identified strain was one of the sampled core populations used to develop the tool, hatchery-identified and DAPC-assigned hatchery-level strains were in agreement in 74.1% of cases in Bangladesh and 80.6% of cases in the Philippines. The dominant hatchery-identified and DAPC-assigned strains were GIFT, in Bangladesh, and GET-ExCEL-a composite strain partially derived from GIFT-in the Philippines.

Ethiopia's transforming wheat landscape: tracking variety use through DNA fingerprinting

Ethiopia is the largest wheat producer in sub-Saharan Africa yet remains a net importer. Increasing domestic wheat production is a national priority. Improved varieties provide an important pathway to enhancing productivity and stability of production. Reliably tracking varietal use and dynamics is a challenge, and the value of conventional recall surveys is increasingly questioned. We report the first nationally representative, large-scale wheat DNA fingerprinting study undertaken in Ethiopia. Plot level comparison of DNA fingerprinting with farmer recall from nearly 4000 plots in the 2016/17 season indicates that only 28% of farmers correctly named wheat varieties grown. The DNA study reveals that new, rust resistant bread wheat varieties are now widely adopted. Germplasm originating from CGIAR centres has made a significant contribution. Corresponding productivity gains and economic benefits have been substantial, indicating high returns to investments in wheat improvement. The study provides an accurate assessment of wheat varietal status and sets a benchmark for national policy-makers and donors. In recent decades, the Ethiopian wheat landscape has transformed from local tetraploid varieties to widespread adoption of high yielding, rust resistant bread wheat. We demonstrate that DNA fingerprinting can be applied at scale and is likely to transform future crop varietal adoption studies.

Population genomics reveals a mismatch between management and biological units in green abalone (Haliotis fulgens)

Effective fishery management strategies should be based on stock delimitation and knowledge of the spatial scale at which species are distributed. However, a mismatch often occurs between biological and management units of fishery resources. The green abalone (Haliotis fulgens) supports an important artisanal fishery in the west coast of the Baja California Peninsula (BCP), Mexico, which has shown a declining tendency despite the several management measures. Thus, the aim of this study was to characterize the spatial patterns of neutral genomic variation of green abalone along the BCP to test whether the genomic structure patterns support the current green abalone management areas. To test this hypothesis, a set of 2,170 putative neutral single nucleotide polymorphisms discovered by a double digest restriction-site associated DNA approach was used on 10 locations along the BCP. The results revealed a population structure with three putative groups: Guadalupe Island and northern and southern BCP locations. The contemporary gene flow might be explained by local oceanographic features, where it is bidirectional within the southern region but with a predominant southward flow from the northern region. These findings indicated that the administrative areas did not match the biological units of H. fulgens fishery; hence, the stock assessment and management areas should be revised.

Genetic diversity and structure of circumtropical almaco jack, Seriola rivoliana: tool for conservation and management

The almaco jack, Seriola rivoliana, is a circumtropical pelagic fish of importance both in commercial fisheries and in aquaculture. To understand levels of genetic diversity within and among populations in the wild, population genetic structure and the relative magnitude of migration were assessed using mtDNA sequence data and single nucleotide polymorphisms (SNPs) from individuals sampled from locations in the Pacific and Atlantic Oceans. A total of 25 variable sites of cytochrome c oxidase subunit 1 and 3678 neutral SNPs were recovered. Three genetic groups were identified, with both marker types distributed in different oceanic regions: Pacific-1 in central Pacific, Pacific-2 in eastern Pacific and Atlantic in western Atlantic. Nonetheless, the analysis of SNP identified a fourth population in the Pacific coast of Baja California Sur, Mexico (Pacific-3), whereas that of mtDNA did not. This mito-nuclear discordance is likely explained by a recently diverged Pacific-3 population. In addition, two mtDNA haplogroups were found within the western Atlantic, likely indicating that the species came into the Atlantic from the Indian Ocean with historical gene flow from the eastern Pacific. Relative gene flow among ocean basins was low with _r m < 0.2, whereas in the eastern Pacific it was asymmetric and higher from south to north (_r m > 0.79). The results reflect the importance of assessing genetic structure and gene flow of natural populations for the purposes of sustainable management.

One panel to rule them all: DArTcap genotyping for population structure, historical demography, and kinship analyses, and its application to a threatened shark

With recent advances in sequencing technology, genomic data are changing how important conservation management decisions are made. Applications such as Close-Kin Mark-Recapture demand large amounts of data to estimate population size and structure, and their full potential can only be realised through ongoing improvements in genotyping strategies. Here we introduce DArTcap, a cost-efficient method that combines DArTseq and sequence capture, and illustrate its use in a high resolution population analysis of Glyphis garricki, a rare, poorly known and threatened euryhaline shark. Clustering analyses and spatial distribution of kin pairs from four different regions across northern Australia and one in Papua New Guinea, representing its entire known range, revealed that each region hosts at least one distinct population. Further structuring is likely within Van Diemen Gulf, the region that included the most rivers sampled, suggesting additional population structuring would be found if other rivers were sampled. Coalescent analyses and spatially explicit modelling suggest that G. garricki experienced a recent range expansion during the opening of the Gulf of Carpentaria following the conclusion of the Last Glacial Maximum. The low migration rates between neighbouring populations of a species that is found only in restricted coastal and riverine habitats show the importance of managing each population separately, including careful monitoring of local and remote anthropogenic activities that may affect their environments. Overall we demonstrated how a carefully chosen SNP panel combined with DArTcap can provide highly accurate kinship inference and also support population structure and historical demography analyses, therefore maximising cost-effectiveness.

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

Genomic methodologies offer unprecedented opportunities for statistically robust studies of species broadly distributed in environments conducive to high gene flow, providing valuable information for wildlife conservation and management. Here, we sequence restriction site-associated DNA to characterize genome-wide single nucleotide polymorphisms (SNPs) in a broadly distributed and highly migratory large pelagic fish, striped marlin (Kajikia audax). Assessment of over 4,000 SNPs resolved spatiotemporal patterns of genetic connectivity throughout the species range in the Pacific and, for the first time, Indian oceans. Individual-based cluster analyses identified six genetically distinct populations corresponding with the western Indian, eastern Indian, western South Pacific, and eastern central Pacific oceans, as well as two populations in the North Pacific Ocean (F ST = 0.0137-0.0819). F ST outlier analyses identified a subset of SNPs (n = 59) putatively under the influence of natural selection and capable of resolving populations separated by comparatively high degrees of genetic differentiation. Temporal collections available for some regions demonstrated the stability of allele frequencies over three to five generations of striped marlin. Relative migration rates reflected lower levels of genetic connectivity between Indian Ocean populations (m R ≤ 0.37) compared with most populations in the Pacific Ocean (m R ≥ 0.57) and highlight the importance of the western South Pacific in facilitating gene flow between ocean basins. Collectively, our results provide novel insights into rangewide population structure for striped marlin and highlight substantial inconsistencies between genetically distinct populations and stocks currently recognized for fisheries management. More broadly, we demonstrate that species capable of long-distance dispersal in environments lacking obvious physical barriers to movement can display substantial population subdivision that persists over multiple generations and that may be facilitated by both neutral and adaptive processes. Importantly, surveys of genome-wide markers enable inference of population-level relationships using sample sizes practical for large pelagic fishes of conservation concern.

Impact of two of the world's largest protected areas on longline fishery catch rates

Two of the largest protected areas on earth are U.S. National Monuments in the Pacific Ocean. Numerous claims have been made about the impacts of these protected areas on the fishing industry, but there has been no ex post empirical evaluation of their effects. We use administrative data documenting individual fishing events to evaluate the economic impact of the expansion of these two monuments on the Hawaii longline fishing fleet. Surprisingly, catch and catch-per-unit-effort are higher since the expansions began. To disentangle the causal effect of the expansions from confounding factors, we use unaffected control fisheries to perform a difference-in-differences analysis. We find that the monument expansions had little, if any, negative impacts on the fishing industry, corroborating ecological models that have predicted minimal impacts from closing large parts of the Pacific Ocean to fishing.

Bioregions in Marine Environments: Combining Biological and Environmental Data for Management and Scientific Understanding

Bioregions are important tools for understanding and managing natural resources. Bioregions should describe locations of relatively homogenous assemblages of species occur, enabling managers to better regulate activities that might affect these assemblages. Many existing bioregionalization approaches, which rely on expert-derived, Delphic comparisons or environmental surrogates, do not explicitly include observed biological data in such analyses. We highlight that, for bioregionalizations to be useful and reliable for systems scientists and managers, the bioregionalizations need to be based on biological data; to include an easily understood assessment of uncertainty, preferably in a spatial format matching the bioregions; and to be scientifically transparent and reproducible. Statistical models provide a scientifically robust, transparent, and interpretable approach for ensuring that bioregions are formed on the basis of observed biological and physical data. Using statistically derived bioregions provides a repeatable framework for the spatial representation of biodiversity at multiple spatial scales. This results in better-informed management decisions and biodiversity conservation outcomes.

Indications of strong adaptive population genetic structure in albacore tuna (Thunnus alalunga) in the southwest and central Pacific Ocean

Albacore tuna (Thunnus alalunga) has a distinctly complex life history in which juveniles and adults separate geographically but at times inhabit the same spaces sequentially. The species also migrates long distances and presumably experiences varied regimes of physical stress over a lifetime. There are, therefore, many opportunities for population structure to arise based on stochastic differences or environmental factors that promote local adaptation. However, with the extent of mobility consistently demonstrated by tagged individuals, there is also a strong argument for panmixia within an ocean basin. It is important to confirm such assumptions from a population genetics standpoint for this species in particular because albacore is one of the principal market tuna species that sustains massive global fisheries and yet is also a slow-growing temperate tuna. Consequently, we used 1,837 neutral SNP loci and 89 loci under potential selection to analyze population genetic structure among five sample groups collected from the western and central South Pacific. We found no evidence to challenge panmixia at neutral loci, but strong indications of structuring at adaptive loci. One population sample, from French Polynesia in 2004, was particularly differentiated. Unfortunately, the current study cannot infer whether the divergence is geographic or temporal, or possibly caused by sample distribution. We encourage future studies to include potentially adaptive loci and to continue fine-scale observations within an ocean basin, and not to assume genome-wide panmixia.

Presence of kin-biased social associations in a lizard with no parental care: the eastern water dragon (Intellagama lesueurii)

Numerous studies have observed kin-biased social associations in a variety of species. Many of these studies have focused on species exhibiting parental care, which may facilitate the transmission of the social environment from parents to offspring. This becomes problematic when disentangling whether kin-biased associations are driven by kin recognition, or are a product of transmission of the social environment during ontogeny, or a combination of both. Studying kin-biased associations in systems that lack parental care may aid in addressing this issue. Furthermore, when studying kin-biased social associations, it is important to differentiate whether these originate from preferential choice or occur randomly as a result of habitat use or limited dispersal. Here, we combined high-resolution single-nucleotide polymorphism data with a long-term behavioral data set of a reptile with no parental care to demonstrate that eastern water dragons (Intellagama lesueurii) bias their nonrandom social associations toward their kin. In particular, we found that although the overall social network was not linked to genetic relatedness, individuals associated with kin more than expected given availability in space and also biased social preferences toward kin. This result opens important opportunities for the study of kinship-driven associations without the confounding effect of vertical transmission of social environments. Furthermore, we present a robust multiple-step approach for determining whether kin-biased social associations are a result of active social decisions or random encounters resulting from habitat use and dispersal patterns.

Standing genomic variation within coding and regulatory regions contributes to the adaptive capacity to climate in a foundation tree species

Global climate is rapidly changing, and the ability for tree species to adapt is dependent on standing genomic variation; however, the distribution and abundance of functional and adaptive variants are poorly understood in natural systems. We test key hypotheses regarding the genetics of adaptive variation in a foundation tree: genomic variation is associated with climate, and genomic variation is more likely to be associated with temperature than precipitation or aridity. To test these hypotheses, we used 9,593 independent, genomic single-nucleotide polymorphisms (SNPs) from 270 individuals sampled from Corymbia calophylla's entire distribution in south-western Western Australia, spanning orthogonal temperature and precipitation gradients. Environmental association analyses returned 537 unique SNPs putatively adaptive to climate. We identified SNPs associated with climatic variation (i.e., temperature [458], precipitation [75] and aridity [78]) across the landscape. Of these, 78 SNPs were nonsynonymous (NS), while 26 SNPs were found within gene regulatory regions. The NS and regulatory candidate SNPs associated with temperature explained more deviance (27.35%) than precipitation (5.93%) and aridity (4.77%), suggesting that temperature provides stronger adaptive signals than precipitation. Genes associated with adaptive variants include functions important in stress responses to temperature and precipitation. Patterns of allelic turnover of NS and regulatory SNPs show small patterns of change through climate space with the exception of an aldehyde dehydrogenase gene variant with 80% allelic turnover with temperature. Together, these findings provide evidence for the presence of adaptive variation to climate in a foundation species and provide critical information to guide adaptive management practices.

Developing noninvasive methodologies to assess koala population health through detecting Chlamydia from scats

Wildlife diseases are a recognized driver of global biodiversity loss, have substantial economic impacts, and are increasingly becoming a threat to human health. Disease surveillance is critical but remains difficult in the wild due to the substantial costs and potential biases associated with most disease detection methods. Noninvasive scat surveys have been proposed as a health monitoring methodology to overcome some of these limitations. Here, we use the known threat of Chlamydia disease to the iconic, yet vulnerable, koala Phascolarctos cinereus to compare three methods for Chlamydia detection in scats: multiplex quantitative PCR, next generation sequencing, and a detection dog specifically trained on scats from Chlamydia-infected koalas. All three methods demonstrated 100% specificity, while sensitivity was variable. Of particular interest is the variable sensitivity of these diagnostic tests to detect sick individuals (i.e., not only infection as confirmed by Chlamydia-positive swabs, but with observable clinical signs of the disease); for koalas with urogenital tract disease signs, sensitivity was 78% with quantitative PCR, 50% with next generation genotyping and 100% with the detection dog method. This may be due to molecular methods having to rely on high-quality DNA whereas the dog most likely detects volatile organic compounds. The most appropriate diagnostic test will vary with disease prevalence and the specific aims of disease surveillance. Acknowledging that detection dogs might not be easily accessible to all, the future development of affordable and portable "artificial noses" to detect diseases from scats in the field might enable cost-effective, rapid and large-scale disease surveillance.

Accounting for kin sampling reveals genetic connectivity in Tasmanian and New Zealand school sharks, Galeorhinus galeus

Fishing represents a major problem for conservation of chondrichthyans, with a quarter of all species being overexploited. School sharks, Galeorhinus galeus, are targeted by commercial fisheries in Australia and New Zealand. The Australian stock has been depleted to below 20% of its virgin biomass, and the species is recorded as Conservation Dependent within Australia. Individuals are known to move between both countries, but it is disputed whether the stocks are reproductively linked. Accurate and unbiased determination of stock and population connectivity is crucial to inform effective management. In this study, we assess the genetic composition and population connectivity between Australian and New Zealand school sharks using genome-wide SNPs, while accounting for non-random kin sampling. Between 2009 and 2013, 88 neonate and juvenile individuals from Tasmanian and New Zealand nurseries were collected and genotyped. Neutral loci were analyzed to detect fine-scale signals of reproductive connectivity. Seven full-sibling groups were identified and removed for unbiased analysis. Based on 6,587 neutral SNPs, pairwise genetic differentiation from Tasmanian and New Zealand neonates was non-significant (F ST = 0.0003, CI95 = [-0.0002, 0.0009], p = 0.1163; D est = 0.0006 ± 0.0002). This pattern was supported by clustering results. In conclusion, we show a significant effect of non-random sampling of kin and identify fine-scale reproductive connectivity between Australian and New Zealand school sharks.

OPEN RESEARCH BADGES

This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://doi.org/10.5061/dryad.pd8612j.

Genetic analyses reveal complex dynamics within a marine fish management area

Genetic data have great potential for improving fisheries management by identifying the fundamental management units-that is, the biological populations-and their mixing. However, so far, the number of practical cases of marine fisheries management using genetics has been limited. Here, we used Atlantic cod in the Baltic Sea to demonstrate the applicability of genetics to a complex management scenario involving mixing of two genetically divergent populations. Specifically, we addressed several assumptions used in the current assessment of the two populations. Through analysis of 483 single nucleotide polymorphisms (SNPs) distributed across the Atlantic cod genome, we confirmed that a model of mechanical mixing, rather than hybridization and introgression, best explained the pattern of genetic differentiation. Thus, the fishery is best monitored as a mixed-stock fishery. Next, we developed a targeted panel of 39 SNPs with high statistical power for identifying population of origin and analyzed more than 2,000 tissue samples collected between 2011 and 2015 as well as 260 otoliths collected in 2003/2004. These data provided high spatial resolution and allowed us to investigate geographical trends in mixing, to compare patterns for different life stages and to investigate temporal trends in mixing. We found similar geographical trends for the two time points represented by tissue and otolith samples and that a recently implemented geographical management separation of the two populations provided a relatively close match to their distributions. In contrast to the current assumption, we found that patterns of mixing differed between juveniles and adults, a signal likely linked to the different reproductive dynamics of the two populations. Collectively, our data confirm that genetics is an operational tool for complex fisheries management applications. We recommend focussing on developing population assessment models and fisheries management frameworks to capitalize fully on the additional information offered by genetically assisted fisheries monitoring.

The population genomics of yellowfin tuna (Thunnus albacares) at global geographic scale challenges current stock delineation

Yellowfin tuna, Thunnus albacares, is one of the most important seafood commodities in the world. Despite its great biological and economic importance, conflicting evidence arises from classical genetic and tagging studies concerning the yellowfin tuna population structure at local and global oceanic scales. Access to more powerful and cost effective genetic tools would represent the first step towards resolving the population structure of yellowfin tuna across its distribution range. Using a panel of 939 neutral Single Nucleotide Polymorphisms (SNPs), and the most comprehensive data set of yellowfin samples available so far, we found genetic differentiation among the Atlantic, Indian and Pacific oceans. The genetic stock structure analysis carried out with 33 outlier SNPs, putatively under selection, identified discrete populations within the Pacific Ocean and, for the first time, also within the Atlantic Ocean. Stock assessment approaches that consider genetic differences at neutral and adaptive genomic loci should be routinely implemented to check the status of the yellowfin tuna, prevent illegal trade, and develop more sustainable management measures.

High seas fisheries play a negligible role in addressing global food security

Recent international negotiations have highlighted the need to protect marine diversity on the high seas-the ocean area beyond national jurisdiction. However, restricting fishing access on the high seas raises many concerns, including how such restrictions would affect food security. We analyze high seas catches and trade data to determine the contribution of the high seas catch to global seafood production, the main species caught on the high seas, and the primary markets where these species are sold. By volume, the total catch from the high seas accounts for 4.2% of annual marine capture fisheries production and 2.4% of total seafood production, including freshwater fisheries and aquaculture. Thirty-nine fish and invertebrate species account for 99.5% of the high seas targeted catch, but only one species, Antarctic toothfish, is caught exclusively on the high seas. The remaining catch, which is caught both on the high seas and in national jurisdictions, is made up primarily of tunas, billfishes, small pelagic fishes, pelagic squids, toothfish, and krill. Most high seas species are destined for upscale food and supplement markets in developed, food-secure countries, such as Japan, the European Union, and the United States, suggesting that, in aggregate, high seas fisheries play a negligible role in ensuring global food security.

Reliably discriminating stock structure with genetic markers: Mixture models with robust and fast computation

Delineating naturally occurring and self-sustaining subpopulations (stocks) of a species is an important task, especially for species harvested from the wild. Despite its central importance to natural resource management, analytical methods used to delineate stocks are often, and increasingly, borrowed from superficially similar analytical tasks in human genetics even though models specifically for stock identification have been previously developed. Unfortunately, the analytical tasks in resource management and human genetics are not identical-questions about humans are typically aimed at inferring ancestry (often referred to as "admixture") rather than breeding stocks. In this article, we argue, and show through simulation experiments and an analysis of yellowfin tuna data, that ancestral analysis methods are not always appropriate for stock delineation. In this work, we advocate a variant of a previously introduced and simpler model that identifies stocks directly. We also highlight that the computational aspects of the analysis, irrespective of the model, are difficult. We introduce some alternative computational methods and quantitatively compare these methods to each other and to established methods. We also present a method for quantifying uncertainty in model parameters and in assignment probabilities. In doing so, we demonstrate that point estimates can be misleading. One of the computational strategies presented here, based on an expectation-maximization algorithm with judiciously chosen starting values, is robust and has a modest computational cost.

Reconciling differences in natural tags to infer demographic and genetic connectivity in marine fish populations

Processes regulating population connectivity are complex, ranging from extrinsic environmental factors to intrinsic individual based features, and are a major force shaping the persistence of fish species and population responses to harvesting and environmental change. Here we developed an integrated assessment of demographic and genetic connectivity of European flounder Platichthys flesus in the northeast Atlantic (from the Norwegian to the Portuguese coast) and Baltic Sea. Specifically, we used a Bayesian infinite mixture model to infer the most likely number of natal sources of individuals based on otolith near core chemical composition. Simultaneously, we characterised genetic connectivity via microsatellite DNA markers, and evaluated how the combined use of natural tags informed individual movement and long-term population exchange rates. Individual markers provided different insights on movement, with otolith chemistry delineating Norwegian and Baltic Sea sources, whilst genetic markers showed a latitudinal pattern which distinguished southern peripheral populations along the Iberian coast. Overall, the integrated use of natural tags resulted in outcomes that were not readily anticipated by individual movement or gene flow markers alone. Our ecological and evolutionary approach provided a synergistic view on connectivity, which will be paramount to align biological and management units and safeguard species' biocomplexity.

From Chromosomes to Genome: Insights into the Evolutionary Relationships and Biogeography of Old World Knifefishes (Notopteridae; Osteoglossiformes)

In addition to its wide geographical distribution, osteoglossiform fishes represent one of the most ancient freshwater teleost lineages; making it an important group for systematic and evolutionary studies. These fishes had a Gondwanan origin and their past distribution may have contributed to the diversity present in this group. However, cytogenetic and genomic data are still scarce, making it difficult to track evolutionary trajectories within this order. In addition, their wide distribution, with groups endemic to different continents, hinders an integrative study that allows a globalized view of its evolutionary process. Here, we performed a detailed chromosomal analysis in Notopteridae fishes, using conventional and advanced molecular cytogenetic methods. Moreover, the genetic distances of examined species were assessed by genotyping using diversity arrays technology sequencing (DArTseq). These data provided a clear picture of the genetic diversity between African and Asian Notopteridae species, and were highly consistent with the chromosomal, geographical, and historical data, enlightening their evolutionary diversification. Here, we discuss the impact of continental drift and split of Pangea on their recent diversity, as well as the contribution to biogeographical models that explain their distribution, highlighting the role of the Indian subcontinent in the evolutionary process within the family.

Fresh is best: Accurate SNP genotyping from koala scats

Maintaining genetic diversity is a crucial component in conserving threatened species. For the iconic Australian koala, there is little genetic information on wild populations that is not either skewed by biased sampling methods (e.g., sampling effort skewed toward urban areas) or of limited usefulness due to low numbers of microsatellites used. The ability to genotype DNA extracted from koala scats using next‐generation sequencing technology will not only help resolve location sample bias but also improve the accuracy and scope of genetic analyses (e.g., neutral vs. adaptive genetic diversity, inbreeding, and effective population size). Here, we present the successful SNP genotyping (1272 SNP loci) of koala DNA extracted from scat, using a proprietary DArTseq^™ protocol. We compare genotype results from two‐day‐old scat DNA and 14‐day‐old scat DNA to a blood DNA template, to test accuracy of scat genotyping. We find that DNA from fresher scat results in fewer loci with missing information than DNA from older scat; however, 14‐day‐old scat can still provide useful genetic information, depending on the research question. We also find that a subset of 209 conserved loci can accurately identify individual koalas, even from older scat samples. In addition, we find that DNA sequences identified from scat samples through the DArTseq^™ process can provide genetic identification of koala diet species, bacterial and viral pathogens, and parasitic organisms.

Widespread gene flow between oceans in a pelagic seabird species complex

Global-scale gene flow is an important concern in conservation biology as it has the potential to either increase or decrease genetic diversity in species and populations. Although many studies focus on the gene flow between different populations of a single species, the potential for gene flow and introgression between species is understudied, particularly in seabirds. The only well-studied example of a mixed-species, hybridizing population of petrels exists on Round Island, in the Indian Ocean. Previous research assumed that Round Island represents a point of secondary contact between Atlantic (Pterodroma arminjoniana) and Pacific species (Pterodroma neglecta and Pterodroma heraldica). This study uses microsatellite genotyping and tracking data to address the possibility of between-species hybridization occurring outside the Indian Ocean. Dispersal and gene flow spanning three oceans were demonstrated between the species in this complex. Analysis of migration rates estimated using bayesass revealed unidirectional movement of petrels from the Atlantic and Pacific into the Indian Ocean. Conversely, structure analysis revealed gene flow between species of the Atlantic and Pacific oceans, with potential three-way hybrids occurring outside the Indian Ocean. Additionally, geolocation tracking of Round Island petrels revealed two individuals travelling to the Atlantic and Pacific. These results suggest that interspecific hybrids in Pterodroma petrels are more common than was previously assumed. This study is the first of its kind to investigate gene flow between populations of closely related Procellariiform species on a global scale, demonstrating the need for consideration of widespread migration and hybridization in the conservation of threatened seabirds.

Identifying hybridization and admixture using SNPs: application of the DArTseq platform in phylogeographic research on vertebrates

Next-generation sequencing (NGS) approaches are increasingly being used to generate multi-locus data for phylogeographic and evolutionary genetics research. We detail the applicability of a restriction enzyme-mediated genome complexity reduction approach with subsequent NGS (DArTseq) in vertebrate study systems at different evolutionary and geographical scales. We present two case studies using SNP data from the DArTseq molecular marker platform. First, we used DArTseq in a large phylogeographic study of the agamid lizard Ctenophorus caudicinctus, including 91 individuals and spanning the geographical range of this species across arid Australia. A low-density DArTseq assay resulted in 28 960 SNPs, with low density referring to a comparably reduced set of identified and sequenced markers as a cost-effective approach. Second, we applied this approach to an evolutionary genetics study of a classic frog hybrid zone (Litoria ewingii-Litoria paraewingi) across 93 individuals, which resulted in 48 117 and 67 060 SNPs for a low- and high-density assay, respectively. We provide a docker-based workflow to facilitate data preparation and analysis, then analyse SNP data using multiple methods including Bayesian model-based clustering and conditional likelihood approaches. Based on comparison of results from the DArTseq platform and traditional molecular approaches, we conclude that DArTseq can be used successfully in vertebrates and will be of particular interest to researchers working at the interface between population genetics and phylogenetics, exploring species boundaries, gene exchange and hybridization.

Null alleles are ubiquitous at microsatellite loci in the Wedge Clam (Donax trunculus)

Recent studies have reported an unusually high frequency of nonamplifying alleles at microsatellite loci in bivalves. Null alleles have been associated with heterozygous deficits in many studies. While several studies have tested for its presence using different analytical tools, few have empirically tested for its consequences in estimating population structure and differentiation. We characterised 16 newly developed microsatellite loci and show that null alleles are ubiquitous in the wedge clam, Donax trunculus. We carried out several tests to demonstrate that the large heterozygous deficits observed in the newly characterised loci were most likely due to null alleles. We tested the robustness of microsatellite genotyping for population assignment by showing that well-recognised biogeographic regions of the south Atlantic and south Mediterranean coast of Spain harbour genetically different populations.

Genomic Differentiation and Demographic Histories of Atlantic and Indo-Pacific Yellowfin Tuna (Thunnus albacares) Populations

Recent developments in the field of genomics have provided new and powerful insights into population structure and dynamics that are essential for the conservation of biological diversity. As a commercially highly valuable species, the yellowfin tuna (Thunnus albacares) is intensely exploited throughout its distribution in tropical oceans around the world, and is currently classified as near threatened. However, conservation efforts for this species have so far been hampered by limited knowledge of its population structure, due to incongruent results of previous investigations. Here, we use whole-genome sequencing in concert with a draft genome assembly to decipher the global population structure of the yellowfin tuna, and to investigate its demographic history. We detect significant differentiation of Atlantic and Indo-Pacific yellowfin tuna populations as well as the possibility of a third diverged yellowfin tuna group in the Arabian Sea. We further observe evidence for past population expansion as well as asymmetric gene flow from the Indo-Pacific to the Atlantic.

Inferring contemporary and historical genetic connectivity from juveniles

Measuring population connectivity is a critical task in conservation biology. While genetic markers can provide reliable long-term historical estimates of population connectivity, scientists are still limited in their ability to determine contemporary patterns of gene flow, the most practical time frame for management. Here, we tackled this issue by developing a new approach that only requires juvenile sampling at a single time period. To demonstrate the usefulness of our method, we used the Speartooth shark (Glyphis glyphis), a critically endangered species of river shark found only in tropical northern Australia and southern Papua New Guinea. Contemporary adult and juvenile shark movements, estimated with the spatial distribution of kin pairs across and within three river systems, was contrasted with historical long-term connectivity patterns, estimated from mitogenomes and genome-wide SNP data. We found strong support for river fidelity in juveniles with the within-cohort relationship analysis. Male breeding movements were highlighted with the cross-cohort relationship analysis, and female reproductive philopatry to the river systems was revealed by the mitogenomic analysis. We show that accounting for juvenile river fidelity and female philopatry is important in population structure analysis and that targeted sampling in nurseries and juvenile aggregations should be included in the genomic toolbox of threatened species management.