Complex patterns of genetic population structure in the mouthbrooding marine catfish, Bagre marinus, in the Gulf of Mexico and U.S. Atlantic.

Patterns of genetic variation reflect interactions among microevolutionary forces that vary in strength with changing demography. Here, patterns of variation within and among samples of the mouthbrooding gafftopsail catfish (Bagre marinus, Family Ariidae) captured in the U.S. Atlantic and throughout the Gulf of Mexico were analyzed using genomics to generate neutral and non-neutral SNP data sets. Because genomic resources are lacking for ariids, linkage disequilibrium network analysis was used to examine patterns of putatively adaptive variation. Finally, historical demographic parameters were estimated from site frequency spectra. The results show four differentiated groups, corresponding to the (1) U.S. Atlantic, and the (2) northeastern, (3) northwestern, and (4) southern Gulf of Mexico. The non-neutral data presented two contrasting signals of structure, one due to increases in diversity moving west to east and north to south, and another to increased heterozygosity in the Atlantic. Demographic analysis suggested that recently reduced long-term effective population size in the Atlantic is likely an important driver of patterns of genetic variation and is consistent with a known reduction in population size potentially due to an epizootic. Overall, patterns of genetic variation resemble that of other fishes that use the same estuarine habitats as nurseries, regardless of the presence/absence of a larval phase, supporting the idea that adult/juvenile behavior and habitat are important predictors of contemporary patterns of genetic structure.

Mesoscale activity drives the habitat suitability of yellowfin tuna in the Gulf of Mexico.

Yellowfin tuna, Thunnus albacares, represents an important component of commercial and recreational fisheries in the Gulf of Mexico (GoM). We investigated the influence of environmental conditions on the spatiotemporal distribution of yellowfin tuna using fisheries' catch data spanning 2012-2019 within Mexican waters. We implemented hierarchical Bayesian regression models with spatial and temporal random effects and fixed effects of several environmental covariates to predict habitat suitability (HS) for the species. The best model included spatial and interannual anomalies of the absolute dynamic topography of the ocean surface (ADTSA and ADTIA, respectively), bottom depth, and a seasonal cyclical random effect. High catches occurred mainly towards anticyclonic features at bottom depths > 1000 m. The spatial extent of HS was higher in years with positive ADTIA, which implies more anticyclonic activity. The highest values of HS (> 0.7) generally occurred at positive ADTSA in oceanic waters of the central and northern GoM. However, high HS values (> 0.6) were observed in the southern GoM, in waters with cyclonic activity during summer. Our results highlight the importance of mesoscale features for the spatiotemporal distribution of yellowfin tunas and could help to develop dynamic fisheries management strategies in Mexico and the U.S. for this valuable resource.

Nursery origin of yellowfin tuna in the western Atlantic Ocean: significance of Caribbean Sea and trans-Atlantic migrants.

Natural geochemical markers in the otolith of yellowfin tuna (Thunnus albacares) were used to establish nursery-specific signatures for investigating the origin of fish captured in the western Atlantic Ocean (WAO). Two classes of chemical markers (trace elements, stable isotopes) were used to first establish nursery-specific signatures of age-0 yellowfin tuna from four primary production zones in the Atlantic Ocean: Gulf of Mexico, Caribbean Sea, Cape Verde, and Gulf of Guinea. Next, mixture and individual assignment methods were applied to predict the origin of sub-adult and adult yellowfin tuna from two regions in the WAO (Gulf of Mexico, Mid Atlantic Bight) by relating otolith core signatures (corresponding to age-0 period) to baseline signatures of age-0 fish from each nursery. Significant numbers of migrants from Caribbean Sea and eastern Atlantic Ocean (EAO) production zones (Gulf of Guinea, Cape Verde) were detected in the WAO, suggesting that fisheries in this region were subsidized by outside spawning/nursery areas. Contributions from local production (Gulf of Mexico) were also evident in samples from both WAO fisheries, but highly variable from year to year. High levels of mixing by yellowfin tuna from the different production zones and pronounced interannual trends in nursery-specific contribution rates in the WAO emphasize the complex and dynamic nature of this species' stock structure and population connectivity. Given that geographic shifts in distribution across national or political boundaries leads to governance and management challenges, this study highlights the need for temporally resolved estimates of nursery origin to refine assessment models and promote the sustainable harvest of this species.

Movement patterns and habitat use of tiger sharks (Galeocerdo cuvier) across ontogeny in the Gulf of Mexico.

The tiger shark (Galeocerdo cuvier) is globally distributed with established coastal and open-ocean movement patterns in many portions of its range. While all life stages of tiger sharks are known to occur in the Gulf of Mexico (GoM), variability in habitat use and movement patterns over ontogeny have never been quantified in this large marine ecosystem. To address this data gap we fitted 56 tiger sharks with Smart Position and Temperature transmitting tags between 2010 and 2018 and examined seasonal and spatial distribution patterns across the GoM. Additionally, we analyzed overlap of core habitats (i.e., 50% kernel density estimates) among individuals relative to large benthic features (oil and gas platforms, natural banks, bathymetric breaks). Our analyses revealed significant ontogenetic and seasonal differences in distribution patterns as well as across-shelf (i.e., regional) and sex-linked variability in movement rates. Presumably sub-adult and adult sharks achieved significantly higher movement rates and used off-shelf deeper habitats at greater proportions than juvenile sharks, particularly during the fall and winter seasons. Further, female maximum rate of movement was higher than males when accounting for size. Additionally, we found evidence of core regions encompassing the National Oceanographic and Atmospheric Administration designated Habitat Areas of Particular Concern (i.e., shelf-edge banks) during cooler months, particularly by females, as well as 2,504 oil and gas platforms. These data provide a baseline for future assessments of environmental impacts, such as climate variability or oil spills, on tiger shark movements and distribution in the region. Future research may benefit from combining alternative tracking tools, such as acoustic telemetry and genetic approaches, which can facilitate long-term assessment of the species' movement dynamics and better elucidate the ecological significance of the core habitats identified here.

Capture stress and post-release mortality of blacktip sharks in recreational charter fisheries of the Gulf of Mexico.

Understanding the stress responses of sharks to recreational catch and release fishing has important management and conservation implications. The blacktip shark Carcharhinus limbatus is a popular recreational species targeted throughout the western, central and eastern Gulf of Mexico (Gulf) yet it is unclear what levels of physiological stress result from catch-release fishing practices with hook and line gear and if the stress levels result in post-release mortality. This study correlates physiological response to stress through blood chemistry analysis and examines post-release behaviour of adult blacktip sharks caught to determine post-release mortality rates. Release behaviour was determined by pop-up satellite archival transmitting (PSAT) tags that record temperature, depth and light level data. To quantify physiological stress levels, blood samples were collected from 52 blacktip sharks and a suite of metabolic and osmotic markers were measured. Thirty-six of those blacktip sharks were also outfitted with a PSAT tag yielding time-at-large from 3 to 180 days. Of the 36 tags, 22 (61%) provided sufficient data to confirm post-release fate and 11 (31%) were recovered providing high-resolution data. Tag data suggests a post-release morality rate of 22.7% (95% confidence interval 7.8-45.4%), with mortality occurring within minutes (immediate mortality) to over 12 h post-release (delayed mortality). Compared to survivors, immediate mortalities exhibited significantly higher lactate (median 2.8 mmol/Lsurvivor vs 5.9 mmol/Limmediate mortality) and significantly lower hematocrit (median 24.4% survivor vs 14% immediate mortality) levels, but no difference was detected between survivors and delayed mortalities. Higher mortality in the western (30%) compared to the central (20%) Gulf may be due to shark handling. All PSATs from mortalities (N = 5) were recovered, and archived data revealed evidence of tag ingestion by predators. Results suggest reduced fight time, decreased handling time and limited air exposure provide blacktip sharks the best survival chances after release by recreational anglers.

Population connectivity of pelagic megafauna in the Cuba-Mexico-United States triangle.

The timing and extent of international crossings by billfishes, tunas, and sharks in the Cuba-Mexico-United States (U.S.) triangle was investigated using electronic tagging data from eight species that resulted in >22,000 tracking days. Transnational movements of these highly mobile marine predators were pronounced with varying levels of bi- or tri-national population connectivity displayed by each species. Billfishes and tunas moved throughout the Gulf of Mexico and all species investigated (blue marlin, white marlin, Atlantic bluefin tuna, yellowfin tuna) frequently crossed international boundaries and entered the territorial waters of Cuba and/or Mexico. Certain sharks (tiger shark, scalloped hammerhead) displayed prolonged periods of residency in U.S. waters with more limited displacements, while whale sharks and to a lesser degree shortfin mako moved through multiple jurisdictions. The spatial extent of associated movements was generally associated with their differential use of coastal and open ocean pelagic ecosystems. Species with the majority of daily positions in oceanic waters off the continental shelf showed the greatest tendency for transnational movements and typically traveled farther from initial tagging locations. Several species converged on a common seasonal movement pattern between territorial waters of the U.S. (summer) and Mexico (winter).

Spatial, temporal, and habitat-related variation in abundance of pelagic fishes in the Gulf of Mexico: potential implications of the deepwater horizon oil spill.

Time-series data collected over a four-year period were used to characterize patterns of abundance for pelagic fishes in the northern Gulf of Mexico (GoM) before (2007-2009) and after (2010) the Deepwater Horizon oil spill. Four numerically dominant pelagic species (blackfin tuna, blue marlin, dolphinfish, and sailfish) were included in our assessment, and larval density of each species was lower in 2010 than any of the three years prior to the oil spill, although larval abundance in 2010 was often statistically similar to other years surveyed. To assess potential overlap between suitable habitat of pelagic fish larvae and surface oil, generalized additive models (GAMs) were developed to evaluate the influence of ocean conditions on the abundance of larvae from 2007-2009. Explanatory variables from GAMs were then linked to environmental data from 2010 to predict the probability of occurrence for each species. The spatial extent of surface oil overlapped with early life habitat of each species, possibly indicating that the availability of high quality habitat was affected by the DH oil spill. Shifts in the distribution of spawning adults is another factor known to influence the abundance of larvae, and the spatial occurrence of a model pelagic predator (blue marlin) was characterized over the same four-year period using electronic tags. The spatial extent of oil coincided with areas used by adult blue marlin from 2007-2009, and the occurrence of blue marlin in areas impacted by the DH oil spill was lower in 2010 relative to pre-spill years.