Reference genome choice and filtering thresholds jointly influence phylogenomic analyses.

Molecular phylogenies are a cornerstone of modern comparative biology and are commonly employed to investigate a range of biological phenomena, such as diversification rates, patterns in trait evolution, biogeography, and community assembly. Recent work has demonstrated that significant biases may be introduced into downstream phylogenetic analyses from processing genomic data; however, it remains unclear whether there are interactions among bioinformatic parameters or biases introduced through the choice of reference genome for sequence alignment and variant-calling. We address these knowledge gaps by employing a combination of simulated and empirical data sets to investigate to what extent the choice of reference genome in upstream bioinformatic processing of genomic data influences phylogenetic inference, as well as the way that reference genome choice interacts with bioinformatic filtering choices and phylogenetic inference method. We demonstrate that more stringent minor allele filters bias inferred trees away from the true species tree topology, and that these biased trees tend to be more imbalanced and have a higher center of gravity than the true trees. We find greatest topological accuracy when filtering sites for minor allele count >3-4 in our 51-taxa data sets, while tree center of gravity was closest to the true value when filtering for sites with minor allele count >1-2. In contrast, filtering for missing data increased accuracy in the inferred topologies; however, this effect was small in comparison to the effect of minor allele filters and may be undesirable due to a subsequent mutation spectrum distortion. The bias introduced by these filters differs based on the reference genome used in short read alignment, providing further support that choosing a reference genome for alignment is an important bioinformatic decision with implications for downstream analyses. These results demonstrate that attributes of the study system and dataset (and their interaction) add important nuance for how best to assemble and filter short read genomic data for phylogenetic inference.

Variable hybridization between two Lake Tanganyikan cichlid species in recent secondary contact.

Closely related taxa frequently exist in sympatry before the evolution of robust reproductive barriers, which can lead to substantial gene flow. Post-divergence gene flow can promote several disparate trajectories of divergence ranging from the erosion of distinctiveness and eventual collapse of the taxa to the strengthening of reproductive isolation. Among many relevant factors, understanding the demographic history of divergence (e.g. divergence time and extent of historical gene flow) can be particularly informative when examining contemporary gene flow between closely related taxa because this history can influence gene flow's prevalence and consequences. Here, we used genotyping-by-sequencing data to investigate speciation and contemporary hybridization in two closely related and sympatrically distributed Lake Tanganyikan cichlid species in the genus Petrochromis. Demographic modelling supported a speciation scenario involving divergence in isolation followed by secondary contact with bidirectional gene flow. Further investigation of this recent gene flow found evidence of ongoing hybridization between the species that varied in extent between different co-occurring populations. Relationships between abundance and the degree of admixture across populations suggest that the availability of conspecific mates may influence patterns of hybridization. These results, together with the observation that sets of recently diverged cichlid taxa are generally geographically separated in the lake, suggest that ongoing speciation in Lake Tanganyikan cichlids relies on initial spatial isolation. Additionally, the spatial heterogeneity of admixture between the Petrochromis species illustrates the complexities of hybridization when species are in recent secondary contact.

Hybridization decreases native cutthroat trout reproductive fitness.

Examining natural selection in wild populations is challenging, but crucial to understanding many ecological and evolutionary processes. Additionally, in hybridizing populations, natural selection may be an important determinant of the eventual outcome of hybridization. We characterized several components of relative fitness in hybridizing populations of Yellowstone cutthroat trout and rainbow trout in an effort to better understand the prolonged persistence of both parental species despite predictions of extirpation. Thousands of genomic loci enabled precise quantification of hybrid status in adult and subsequent juvenile generations; a subset of those data also identified parent-offspring relationships. We used linear models and simulations to assess the effects of ancestry on reproductive output and mate choice decisions. We found a relatively low number of late-stage (F3+) hybrids and an excess of F2 juveniles relative to the adult generation in one location, which suggests the presence of hybrid breakdown decreasing the fitness of F2+ hybrids later in life. Assessments of reproductive output showed that Yellowstone cutthroat trout are more likely to successfully reproduce and produce slightly more offspring than their rainbow trout and hybrid counterparts. Mate choice appeared to be largely random, though we did find statistical support for slight female preference for males of similar ancestry. Together, these results show that native Yellowstone cutthroat trout are able to outperform rainbow trout in terms of reproduction and suggest that management action to exclude rainbow trout from spawning locations may bolster the now-rare Yellowstone cutthroat trout.

The Genetic Population Structure of Lake Tanganyika's Lates Species Flock, an Endemic Radiation of Pelagic Top Predators.

Understanding genetic connectivity plays a crucial role in species conservation decisions, and genetic connectivity is an important component of modern fisheries management. In this study, we investigated the population genetics of four endemic Lates species of Lake Tanganyika (Lates stappersii, L. microlepis, L. mariae, and L. angustifrons) using reduced-representation genomic sequencing methods. We find the four species to be strongly differentiated from one another (mean interspecific FST = 0.665), with no evidence for contemporary admixture. We also find evidence for strong genetic structure within L. mariae, with the majority of individuals from the most southern sampling site forming a genetic group that is distinct from the individuals at other sampling sites. We find evidence for much weaker structure within the other three species (L. stappersii, L. microlepis, and L. angustifrons). Our ability to detect this weak structure despite small and unbalanced sample sizes and imprecise geographic sampling locations suggests the possibility for further structure undetected in our study. We call for further research into the origins of the genetic differentiation in these four species-particularly that of L. mariae-which may be important for conservation and management of this culturally and economically important clade of fishes.

Structural genomic variation leads to genetic differentiation in Lake Tanganyika's sardines.

Identifying patterns in genetic structure and the genetic basis of ecological adaptation is a core goal of evolutionary biology and can inform the management and conservation of species that are vulnerable to population declines exacerbated by climate change. We used reduced-representation genomic sequencing methods to gain a better understanding of genetic structure among and within populations of Lake Tanganyika's two sardine species, Limnothrissa miodon and Stolothrissa tanganicae. Samples of these ecologically and economically important species were collected across the length of Lake Tanganyika, as well as from nearby Lake Kivu, where L. miodon was introduced in 1959. Our results reveal differentiation within both S. tanganicae and L. miodon that is not explained by geography. Instead, this genetic differentiation is due to the presence of large sex-specific regions in the genomes of both species, but involving different polymorphic sites in each species. Our results therefore indicate rapidly evolving XY sex determination in the two species. Additionally, we found evidence of a large chromosomal rearrangement in L. miodon, creating two homokaryotypes and one heterokaryotype. We found all karyotypes throughout Lake Tanganyika, but the frequencies vary along a north-south gradient and differ substantially in the introduced Lake Kivu population. We do not find evidence for significant isolation by distance, even over the hundreds of kilometres covered by our sampling, but we do find shallow population structure.

Genomics and the origin of species.

Speciation is a fundamental evolutionary process, the knowledge of which is crucial for understanding the origins of biodiversity. Genomic approaches are an increasingly important aspect of this research field. We review current understanding of genome-wide effects of accumulating reproductive isolation and of genomic properties that influence the process of speciation. Building on this work, we identify emergent trends and gaps in our understanding, propose new approaches to more fully integrate genomics into speciation research, translate speciation theory into hypotheses that are testable using genomic tools and provide an integrative definition of the field of speciation genomics.

The genomic substrate for adaptive radiation in African cichlid fish.

Cichlid fishes are famous for large, diverse and replicated adaptive radiations in the Great Lakes of East Africa. To understand the molecular mechanisms underlying cichlid phenotypic diversity, we sequenced the genomes and transcriptomes of five lineages of African cichlids: the Nile tilapia (Oreochromis niloticus), an ancestral lineage with low diversity; and four members of the East African lineage: Neolamprologus brichardi/pulcher (older radiation, Lake Tanganyika), Metriaclima zebra (recent radiation, Lake Malawi), Pundamilia nyererei (very recent radiation, Lake Victoria), and Astatotilapia burtoni (riverine species around Lake Tanganyika). We found an excess of gene duplications in the East African lineage compared to tilapia and other teleosts, an abundance of non-coding element divergence, accelerated coding sequence evolution, expression divergence associated with transposable element insertions, and regulation by novel microRNAs. In addition, we analysed sequence data from sixty individuals representing six closely related species from Lake Victoria, and show genome-wide diversifying selection on coding and regulatory variants, some of which were recruited from ancient polymorphisms. We conclude that a number of molecular mechanisms shaped East African cichlid genomes, and that amassing of standing variation during periods of relaxed purifying selection may have been important in facilitating subsequent evolutionary diversification.

Playing out Liem's Paradox: Opportunistic Piscivory across Lake Tanganyikan Cichlids.

Trophic specialization is a key feature of the diversity of cichlid fish adaptive radiations. However, K. F. Liem observed that even species with highly specialized trophic morphologies have dietary flexibility, enabling them to exploit episodic food resources opportunistically. Evidence for dietary flexibility comes largely from laboratory studies, and it is unclear whether cichlid fishes undergo diet shifts in the wild. We report observations of diet switching by multiple cichlid species in Lake Tanganyika as a consequence of unusual concentrations of schooling juvenile clupeid fishes. Fish species with varying degrees of trophic specialization converged on a single prey: juvenile sardines that are also endemic to Lake Tanganyika (Stolothrissa tanganicae and Limnothrissa miodon). We provide evidence for cichlid species acting as jacks-of-all-trades and discuss this evidence in the framework of Liem's classic paradox: that trophic specialization does not preclude dietary flexibility.

Variable hybridization outcomes in trout are predicted by historical fish stocking and environmental context.

Hybridization can profoundly affect the genomic composition and phenotypes of closely related species, and provides an opportunity to identify mechanisms that maintain reproductive isolation between species. Recent evidence suggests that hybridization outcomes within a species pair can vary across locations. However, we still do not know how variable outcomes of hybridization are across geographic replicates, and what mechanisms drive that variation. In this study, we described hybridization outcomes across 27 locations in the North Fork Shoshone River basin (Wyoming, USA) where native Yellowstone cutthroat trout and introduced rainbow trout co-occur. We used genomic data and hierarchical Bayesian models to precisely identify ancestry of hybrid individuals. Hybridization outcomes varied across locations. In some locations, only rainbow trout and advanced backcrossed hybrids towards rainbow trout were present, while trout in other locations had a broader range of ancestry, including both parental species and first-generation hybrids. Later-generation intermediate hybrids were rare relative to backcrossed hybrids and rainbow trout individuals. Using an individual-based simulation, we found that outcomes of hybridization in the North Fork Shoshone River basin deviate substantially from what we would expect under null expectations of random mating and no selection against hybrids. Since this deviation implies that some mechanisms of reproductive isolation function to maintain parental taxa and a diversity of hybrid types, we then modelled hybridization outcomes as a function of environmental variables and stocking history that are likely to affect prezygotic barriers to hybridization. Variables associated with history of fish stocking were the strongest predictors of hybridization outcomes, followed by environmental variables that might affect overlap in spawning time and location.

Genomics of Rapid Incipient Speciation in Sympatric Threespine Stickleback.

Ecological speciation is the process by which reproductively isolated populations emerge as a consequence of divergent natural or ecologically-mediated sexual selection. Most genomic studies of ecological speciation have investigated allopatric populations, making it difficult to infer reproductive isolation. The few studies on sympatric ecotypes have focused on advanced stages of the speciation process after thousands of generations of divergence. As a consequence, we still do not know what genomic signatures of the early onset of ecological speciation look like. Here, we examined genomic differentiation among migratory lake and resident stream ecotypes of threespine stickleback reproducing in sympatry in one stream, and in parapatry in another stream. Importantly, these ecotypes started diverging less than 150 years ago. We obtained 34,756 SNPs with restriction-site associated DNA sequencing and identified genomic islands of differentiation using a Hidden Markov Model approach. Consistent with incipient ecological speciation, we found significant genomic differentiation between ecotypes both in sympatry and parapatry. Of 19 islands of differentiation resisting gene flow in sympatry, all were also differentiated in parapatry and were thus likely driven by divergent selection among habitats. These islands clustered in quantitative trait loci controlling divergent traits among the ecotypes, many of them concentrated in one region with low to intermediate recombination. Our findings suggest that adaptive genomic differentiation at many genetic loci can arise and persist in sympatry at the very early stage of ecotype divergence, and that the genomic architecture of adaptation may facilitate this.

The smelly path to sympatric speciation?

In this issue of Molecular Ecology, Poelstra et al. () use genomic data to show that cichlid species in the adaptive radiation in Lake Ejagham, Cameroon, experienced gene flow with a riverine relative, primarily prior to their diversification. Intriguingly, this introgression brought with it olfactory alleles that the authors suggest may play a sensory role in speciation. As a classic example of sympatric speciation due to the highly restricted geography of this small (0.49 km2 ) crater lake (Schliewen et al., ), this result sheds new light on the history of this fascinating radiation. As genomic data for cichlid radiations accumulate, finding evidence of introgression increasingly appears to be the rule rather than the exception (Kautt et al., ; Meier et al., , ). However, two points are pressing as evidence for introgression mounts: (a) It is crucial to understand the assumptions of analytical approaches used and (b) it is important to think clearly about the divergent contexts in which evidence for introgression has been invoked as an important feature of adaptive radiation.

Divergent parasite infections in sympatric cichlid species in Lake Victoria.

Parasitism has been proposed as a factor in host speciation, as an agent affecting coexistence of host species in species-rich communities and as a driver of post-speciation diversification. Young adaptive radiations of closely related host species of varying ecological and genomic differentiation provide interesting opportunities to explore interactions between patterns of parasitism, divergence and coexistence of sympatric host species. Here, we explored patterns in ectoparasitism in a community of 16 fully sympatric cichlid species at Makobe Island in Lake Victoria, a model system of vertebrate adaptive radiation. We asked whether host niche, host abundance or host genetic differentiation explains variation in infection patterns. We found significant differences in infections, the magnitude of which was weakly correlated with the extent of genomic divergence between the host species, but more strongly with the main ecological gradient, water depth. These effects were most evident with infections of Cichlidogyrus monogeneans, whereas the only host species with a strictly crevice-dwelling niche, Pundamilia pundamilia, deviated from the general negative relationship between depth and parasitism. In accordance with the Janzen-Connell hypothesis, we also found that host abundance tended to be positively associated with infections in some parasite taxa. Data on the Pundamilia sister species pairs from three other islands with variable degrees of habitat (crevice) specialization suggested that the lower parasite abundance of P. pundamilia at Makobe could result from both habitat specialization and the evolution of specific resistance. Our results support influences of host genetic differentiation and host ecology in determining infections in this diverse community of sympatric cichlid species.

Ecological opportunity and sexual selection together predict adaptive radiation.

A fundamental challenge to our understanding of biodiversity is to explain why some groups of species undergo adaptive radiations, diversifying extensively into many and varied species, whereas others do not. Both extrinsic environmental factors (for example, resource availability, climate) and intrinsic lineage-specific traits (for example, behavioural or morphological traits, genetic architecture) influence diversification, but few studies have addressed how such factors interact. Radiations of cichlid fishes in the African Great Lakes provide some of the most dramatic cases of species diversification. However, most cichlid lineages in African lakes have not undergone adaptive radiations. Here we compile data on cichlid colonization and diversification in 46 African lakes, along with lake environmental features and information about the traits of colonizing cichlid lineages, to investigate why adaptive radiation does and does not occur. We find that extrinsic environmental factors related to ecological opportunity and intrinsic lineage-specific traits related to sexual selection both strongly influence whether cichlids radiate. Cichlids are more likely to radiate in deep lakes, in regions with more incident solar radiation and in lakes where there has been more time for diversification. Weak or negative associations between diversification and lake surface area indicate that cichlid speciation is not constrained by area, in contrast to diversification in many terrestrial taxa. Among the suite of intrinsic traits that we investigate, sexual dichromatism, a surrogate for the intensity of sexual selection, is consistently positively associated with diversification. Thus, for cichlids, it is the coincidence between ecological opportunity and sexual selection that best predicts whether adaptive radiation will occur. These findings suggest that adaptive radiation is predictable, but only when species traits and environmental factors are jointly considered.

Pattern and Process in the Comparative Study of Convergent Evolution.

Understanding processes that have shaped broad-scale biodiversity patterns is a fundamental goal in evolutionary biology. The development of phylogenetic comparative methods has yielded a tool kit for analyzing contemporary patterns by explicitly modeling processes of change in the past, providing neontologists tools for asking questions previously accessible only for select taxa via the fossil record or laboratory experimentation. The comparative approach, however, differs operationally from alternative approaches to studying convergence in that, for studies of only extant species, convergence must be inferred using evolutionary process models rather than being directly measured. As a result, investigation of evolutionary pattern and process cannot be decoupled in comparative studies of convergence, even though such a decoupling could in theory guard against adaptationist bias. Assumptions about evolutionary process underlying comparative tools can shape the inference of convergent pattern in sometimes profound ways and can color interpretation of such patterns. We discuss these issues and other limitations common to most phylogenetic comparative approaches and suggest ways that they can be avoided in practice. We conclude by promoting a multipronged approach to studying convergence that integrates comparative methods with complementary tests of evolutionary mechanisms and includes ecological and biogeographical perspectives. Carefully employed, the comparative method remains a powerful tool for enriching our understanding of convergence in macroevolution, especially for investigation of why convergence occurs in some settings but not others.

Demographic modelling with whole-genome data reveals parallel origin of similar Pundamilia cichlid species after hybridization.

Modes and mechanisms of speciation are best studied in young species pairs. In older taxa, it is increasingly difficult to distinguish what happened during speciation from what happened after speciation. Lake Victoria cichlids in the genus Pundamilia encompass a complex of young species and polymorphic populations. One Pundamilia species pair, P. pundamilia and P. nyererei, is particularly well suited to study speciation because sympatric population pairs occur with different levels of phenotypic differentiation and reproductive isolation at different rocky islands within the lake. Genetic distances between allopatric island populations of the same nominal species often exceed those between the sympatric species. It thus remained unresolved whether speciation into P. nyererei and P. pundamilia occurred once, followed by geographical range expansion and interspecific gene flow in local sympatry, or if the species pair arose repeatedly by parallel speciation. Here, we use genomic data and demographic modelling to test these alternative evolutionary scenarios. We demonstrate that gene flow plays a strong role in shaping the observed patterns of genetic similarity, including both gene flow between sympatric species and gene flow between allopatric populations, as well as recent and early gene flow. The best supported model for the origin of P. pundamilia and P. nyererei population pairs at two different islands is one where speciation happened twice, whereby the second speciation event follows shortly after introgression from an allopatric P. nyererei population that arose earlier. Our findings support the hypothesis that very similar species may arise repeatedly, potentially facilitated by introgressed genetic variation.

Genomic landscape of early ecological speciation initiated by selection on nuptial colour.

Ecological speciation is the evolution of reproductive isolation as a consequence of direct divergent natural selection or ecologically mediated divergent sexual selection. While the genomic signature of the former has been extensively studied in recent years, only few examples exist for genomic differentiation where environment-dependent sexual selection has played an important role. Here, we describe a very young (~90 years old) population of threespine sticklebacks exhibiting phenotypic and genomic differentiation between two habitats within the same pond. We show that differentiation among habitats is limited to male throat colour and nest type, traits known to be subject to sexual selection. Divergence in these traits mirrors divergence in much older benthic and limnetic stickleback species pairs from North American west coast lakes, which also occur in sympatry but are strongly reproductively isolated from each other. We demonstrate that in our population, differences in throat colour and breeding have been stable over a decade, but in contrast to North American benthic and limnetic stickleback species, these mating trait differences are not accompanied by divergence in morphology related to feeding, predator defence or swimming performance. Using genomewide SNP data, we find multiple genomic islands with moderate differentiation spread across several chromosomes, whereas the rest of the genome is undifferentiated. The islands contain potential candidate genes involved in visual perception of colour. Our results suggest that phenotypic and multichromosome genomic divergence of these morphs was driven by environment-dependent sexual selection, demonstrating incipient speciation after only a few decades of divergence in sympatry.

Cichlid species-area relationships are shaped by adaptive radiations that scale with area.

A positive relationship between species richness and island size is thought to emerge from an equilibrium between immigration and extinction rates, but the influence of species diversification on the form of this relationship is poorly understood. Here, we show that within-lake adaptive radiation strongly modifies the species-area relationship for African cichlid fishes. The total number of species derived from in situ speciation increases with lake size, resulting in faunas orders of magnitude higher in species richness than faunas assembled by immigration alone. Multivariate models provide evidence for added influence of lake depth on the species-area relationship. Diversity of clades representing within-lake radiations show responses to lake area, depth and energy consistent with limitation by these factors, suggesting that ecological factors influence the species richness of radiating clades within these ecosystems. Together, these processes produce lake fish faunas with highly variable composition, but with diversities that are well predicted by environmental variables.

Genome-wide RAD sequence data provide unprecedented resolution of species boundaries and relationships in the Lake Victoria cichlid adaptive radiation.

Although population genomic studies using next generation sequencing (NGS) data are becoming increasingly common, studies focusing on phylogenetic inference using these data are in their infancy. Here, we use NGS data generated from reduced representation genomic libraries of restriction-site-associated DNA (RAD) markers to infer phylogenetic relationships among 16 species of cichlid fishes from a single rocky island community within Lake Victoria's cichlid adaptive radiation. Previous attempts at sequence-based phylogenetic analyses in Victoria cichlids have shown extensive sharing of genetic variation among species and no resolution of species or higher-level relationships. These patterns have generally been attributed to the very recent origin (<15,000 years) of the radiation, and ongoing hybridization between species. We show that as we increase the amount of sequence data used in phylogenetic analyses, we produce phylogenetic trees with unprecedented resolution for this group. In trees derived from our largest data supermatrices (3 to >5.8 million base pairs in width), species are reciprocally monophyletic with high bootstrap support, and the majority of internal branches on the tree have high support. Given the difficulty of the phylogenetic problem that the Lake Victoria cichlid adaptive radiation represents, these results are striking. The strict interpretation of the topologies we present here warrants caution because many questions remain about phylogenetic inference with very large genomic data set and because we can with the current analysis not distinguish between effects of shared ancestry and post-speciation gene flow. However, these results provide the first conclusive evidence for the monophyly of species in the Lake Victoria cichlid radiation and demonstrate the power that NGS data sets hold to resolve even the most difficult of phylogenetic challenges.

Parallel shifts in trout feeding morphology suggest rapid adaptation to alpine lake environments.

Eco-evolutionary interactions following ecosystem change provide critical insight into the ability of organisms to adapt to shifting resource landscapes. Here we explore evidence for the rapid parallel evolution of trout feeding morphology following eco-evolutionary interactions with zooplankton in alpine lakes stocked at different points in time in the Wind River Range (Wyoming, USA). In this system, trout predation has altered the zooplankton species community and driven a decrease in average zooplankton size. In some lakes that were stocked decades ago, we find shifts in gill raker traits consistent with the hypothesis that trout have rapidly adapted to exploit available smaller-bodied zooplankton more effectively. We explore this morphological response in multiple lake populations across two species of trout (cutthroat trout, Oncorhynchus clarkii, and golden trout Oncorhynchus aguabonita) and examine the impact of resource availability on morphological variation in gill raker number among lakes. Furthermore, we present genetic data to provide evidence that historically stocked cutthroat trout populations likely derive from multiple population sources, and incorporate variation from genomic relatedness in our exploration of environmental predictors of feeding morphology. These findings describe rapid adaptation and eco-evolutionary interactions in trout and document an evolutionary response to novel, contemporary ecosystem change.

Isotopic signatures induced by upwelling reveal regional fish stocks in Lake Tanganyika.

Lake Tanganyika's pelagic fish sustain the second largest inland fishery in Africa and are under pressure from heavy fishing and global warming related increases in stratification. The strength of water column stratification varies regionally, with a more stratified north and an upwelling-driven, biologically more productive south. Only little is known about whether such regional hydrodynamic regimes induce ecological or genetic differences among populations of highly mobile, pelagic fish inhabiting these different areas. Here, we examine whether the regional contrasts leave distinct isotopic imprints in the pelagic fish of Lake Tanganyika, which may reveal differences in diet or lipid content. We conducted two lake-wide campaigns during different seasons and collected physical, nutrient, chlorophyll, phytoplankton and zooplankton data. Additionally, we analyzed the pelagic fish-the clupeids Stolothrissa tanganicae, Limnothrissa miodon and four Lates species-for their isotopic and elemental carbon (C) and nitrogen (N) compositions. The δ13C values were significantly higher in the productive south after the upwelling/mixing period across all trophic levels, implying that the fish have regional foraging grounds, and thus record these latitudinal isotope gradients. By combining our isotope data with previous genetic results showing little geographic structure, we demonstrate that the fish reside in a region for a season or longer. Between specimens from the north and south we found no strong evidence for varying trophic levels or lipid contents, based on their bulk δ15N and C:N ratios. We suggest that the development of regional trophic or physiological differences may be inhibited by the lake-wide gene flow on the long term. Overall, our findings show that the pelagic fish species, despite not showing evidence for genetic structure at the basin scale, form regional stocks at the seasonal timescales. This implies that sustainable management strategies may consider adopting regional fishing quotas.