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.

Ventilation of the deep Gulf of Mexico and potential insights to the Atlantic Meridional Overturning Circulation.

Ventilation in the deep Gulf of Mexico (GoM), its connection to the North Atlantic, and its susceptibility to changes of the Atlantic Meridional Overturning Circulation are investigated by combining observations of radiocarbon and volume transport with a Coupled Model Intercomparison Project phase 6 (CMIP6) General Circulation Model (GCM) ensemble output. Radiocarbon data and multiannual volume transport through the Yucatan Channel suggest deep water residence times ~100 years for the GoM. Comparisons to previous radiocarbon observations suggest that the deep GoM has aged in the recent past, consistent with observed raising temperatures and the CMIP6 GCM simulations. The distribution of radiocarbon indicates a time frame of ~160 years between North Atlantic convection and complete ventilation of the deep GoM. This suggests that aging and warming of GoM deep waters were initiated in the North Atlantic before ~1890 consistent with reported rapid and persistent AMOC weakening since the Little Ice Age.

Distribution and densities of fish larvae species with contrasting life histories as a function of oceanographic variables in the deep-water region of the southern Gulf of Mexico.

We describe the larval occurrence and density of six fish species with contrasting life histories and examine their relationships with oceanographic variables during two seasons in the deep-water region (> 1000 m) of the southern Gulf of Mexico based on 12 cruises (2011-2018). Given that Caranx crysos adults are neritic, larval presence close to the continental shelf indicates offshore cross-shelf transport to oceanic waters, which likely leads to mortality. Generalized additive models indicated that C. crysos density was not related with oceanographic variables, while Auxis spp. (with neritic and oceanic adults) was related to wind speed, sea surface temperature, sea surface height, and surface chlorophyll a. The mesopelagic Benthosema suborbitale, Notolychnus valdiviae and Bregmaceros atlanticus were more abundant and broadly distributed, and higher density was found in conditions indicative of higher nutrient availability and productivity, suggesting greater feeding success and survival. The distribution of the epi- and mesopelagic Cubiceps pauciradiatus extended through the southern Gulf of Mexico, and was related to wind speed, sea surface temperature, stratification and chlorophyll a. Our results suggest that the density of the neritic species in oceanic waters could be mediated by regional cross-shelf transport, while for oceanic species is linked with productivity.

Environmental conditions drive zooplankton community structure in the epipelagic oceanic water of the southern Gulf of Mexico: A molecular approach.

Zooplankton plays a pivotal role in sustaining the majority of marine ecosystems. The distribution patterns and diversity of zooplankton provide key information for understanding the functioning of these ecosystems. Nevertheless, due to the numerous cryptic and sibling species and the lack of diagnostic characteristics for early developmental stages, the identification of the global-to-local patterns of zooplankton biodiversity and biogeography remains challenging in different research fields. The spatial and temporal changes in the zooplankton community in the open waters of the southern Gulf of Mexico were assessed using metabarcoding analysis of the V9 region of 18S rRNA and mitochondrial cytochrome oxidase c subunit I (COI). Additionally, a multiscale analysis was implemented to evaluate which environmental predictors may explain the variability in the structure of the zooplankton community. Our findings suggest that the synergistic effects of dissolved oxygen concentration, temperature, and longitude (intended as a proxy for still unidentified predictors) may explain both spatial and temporal zooplankton variability even with low contribution. Furthermore, the zooplankton distribution probably reflects the coexistence of three heterogeneous ecoregions and a bio-physical partitioning of the studied area. Finally, some taxa were either exclusive or predominant with either 18S or COI markers. This may suggest that comprehensive assessments of the zooplankton community may be more accurately met by the use of multilocus approaches.

Yellowfin tuna (Thunnus albacares) foraging habitat and trophic position in the Gulf of Mexico based on intrinsic isotope tracers.

Yellowfin tuna (YFT, Thunnus albacares) is a commercially important species targeted by fisheries in the Gulf of Mexico (GM). Previous studies suggest a high degree of residency in the northern GM, although part of the population performs movements to southern Mexican waters. Whether YFT caught in southern waters also exhibit residency or migrate to the northern gulf is currently uncertain, and little is known regarding their trophic ecology. The isotopic composition (bulk & amino acids) of YFT muscle and liver tissues were compared to a zooplankton-based synoptic isoscape from the entire GM to infer feeding areas and estimate Trophic Position (TP). The spatial distribution of δ15Nbulk and δ15NPhe values of zooplankton indicated two distinct isotopic baselines: one with higher values in the northern GM likely driven by denitrification over the continental shelf, and another in the central-southern gulf, where nitrogen fixation predominates. Based on the contribution of the two regional isotopic baselines to YFT tissues, broad feeding areas were inferred, with a greater contribution of the northern GM (over a one-year time scale by muscle), and to a lesser extent in the central-southern GM (over the ca. 6-month scale by liver). This was corroborated by similarities in δ15NPhe values between YFT and the northern GM. TP estimates were calculated based on stable isotope analysis of bulk (SIA) and compound-specific isotope analysis (CSIA-AA) of the canonical source and trophic amino acids. Mean TP based on SIA was 4.9 ± 1.0 and mean TP based on CSIA-A was 3.9 ± 0.2. YFT caught within the Mexican region seem to feed in northern and in central and southern GM, while feeding in the northern GM has a temporal component. Thus, management strategies need to consider that YFT caught in US and Mexican waters are a shared binational resource that exhibit feeding migrations within the GM.

Trophic ecology of Mexican Pacific harbor seal colonies using carbon and nitrogen stable isotopes.

There is limited information that provides a comprehensive understanding of the trophic ecology of Mexican Pacific harbor seal (Phoca vitulina richardii) colonies. While scat analysis has been used to determine the diet of some colonies, the integrative characterization of its feeding habits on broader temporal and spatial scales remains limited. We examined potential feeding grounds, trophic niche width, and overlap, and inferred the degree of dietary specialization using stable carbon and nitrogen isotope ratios (δ13C and δ15N) in this subspecies. We analyzed δ13C and δ15N on fur samples from pups collected at five sites along the western coast of the Baja California Peninsula, Mexico. Fur of natal coat of Pacific harbor seal pups begins to grow during the seventh month in utero until the last stage of gestation. Therefore pup fur is a good proxy for the mother's feeding habits in winter (~December to March), based on the timing of gestation for the subspecies in this region. Our results indicated that the δ13C and δ15N values differed significantly among sampling sites, with the highest mean δ15N value occurring at the southernmost site, reflecting a well-characterized north to south latitudinal 15N-enrichment in the food web. The tendency identified in δ13C values, in which the northern colonies showed the most enriched values, suggests nearshore and benthic-demersal feeding habits. A low variance in δ13C and δ15N values for each colony (<1‰) and relatively small standard ellipse areas suggest a specialized foraging behavior in adult female Pacific harbor seals in Mexican waters.

Metazoan parasite infracommunities of the dusky flounder (Syacium papillosum) as bioindicators of environmental conditions in the continental shelf of the Yucatan Peninsula, Mexico.

BACKGROUND: We assessed metrics of the metazoan parasite infracommunities of the dusky flounder (Syacium papillosum) as indicators of aquatic environmental health of the Yucatan Shelf (YS) prior to oil extraction. We sampled the dusky flounder and its parasites along the YS, mostly during the 2015 north wind season (November-April). Our aims were: (i) to determine whether the parasite infracommunity metrics of S. papillosum exhibit significant differences among YS subregions; (ii) to determine whether the probability of the occurrence of its parasite species and individuals were affected by environmental variables, nutrients, heavy metals and hydrocarbons at the seascape level; and (iii) to determine whether there were statistical differences between the parasite infracommunity metrics of S. papillosum from YS and those of Syacium gunteri from the Campeche Sound. Multivariate statistical analyses and generalised additive models (GAMs) were used to examine the potential statistical associations between the contaminants, environmental variables and parasite community metrics, and the maximum entropy algorithm (MaxEnt) was used to characterise the habitat's suitability for the parasite's probability of occurrence. RESULTS: We recovered 48 metazoan parasite species from 127 S. papillosum, with larval cestodes and digeneans being the most numerically-dominant. Multivariate analyses showed significant differences in parasite infracommunity metrics among Western YS, Mid YS and Caribbean subregions, with the latter being the richest in species but not in individuals. The GAM and MaxEnt results indicated a negative effect of top predators (e.g. sharks and rays) removal on parasite metrics. The parasite infracommunities of S. papillosum were twice as rich in the number of species and individuals as those reported for S. gunteri from the Campeche Sound. CONCLUSIONS: The significant differences among subregions in parasite metrics were apparently due to the interruption of the Yucatan current during the north wind season. The fishing of top predators in combination with an influx of nutrients and hydrocarbons in low concentrations coincides with an increase in larval cestodes and digeneans in S. papillosum. The dusky flounder inhabits a region (YS) with a larger number of metazoan parasite species compared with those available for S. gunteri in the Campeche Sound, suggesting better environmental conditions for transmission in the YS.

Elements of time and place: manganese and barium in shark vertebrae reflect age and upwelling histories.

As upper-level predators, sharks are important for maintaining marine food web structure, but populations are threatened by fishery exploitation. Sustainable management of shark populations requires improved understanding of migration patterns and population demographics, which has traditionally been sought through physical and/or electronic tagging studies. The application of natural tags such as elemental variations in mineralized band pairs of elasmobranch vertebrae cartilage could also reveal endogenous and exogenous processes experienced by sharks throughout their life histories. Here, elemental profiles were characterized in vertebrae encompassing complete life histories (birth-to-death) of shortfin mako (Isurus oxyrinchus), common thresher (Alopias vulpinus) and blue shark (Prionace glauca) of known tag and recapture locations in the eastern North Pacific Ocean. All sharks were injected with oxytetracycline at initial capture, released and subsequently recaptured, with individual liberty times ranging from 215 days to 6 years. Vertebral band pairs forming over the liberty intervals were verified by counting the number of band pairs deposited since the oxytetracycline band. Regular oscillations in vertebrae manganese (Mn) content corresponded well with the number of validated band pairs, suggesting that Mn variation could be used to age sharks. Increases in vertebrae barium concentration were correlated with times when individuals occupied areas with high coastal upwelling indices, the timing and spatial intensity of which varied from year to year. Interspecific relationships were probably influenced by behavioural differences in horizontal and vertical habitat use, feeding habits and thermoregulatory physiology. These results indicate that vertebral sclerochronology has the potential to advance our knowledge of elasmobranch life history including age and growth estimation and environmental reconstruction.

Amino acid-specific δ15N trophic enrichment factors in fish fed with formulated diets varying in protein quantity and quality.

Compound-specific isotope analysis (CSIA) of amino acids (AAs) in consumer tissues is a developing technique with wide-ranging applications for identifying nitrogen (N) sources and estimating animal trophic level. Controlled experiments are essential for determining which dietary conditions influence variability in N stable isotopes (δ15N) trophic enrichment factors in bulk tissue (TEFbulk) and AAs (TEFAA). To date, however, studies have not independently evaluated the effect of protein quantity and quality (digestibility) on TEFs, complicating the application of AA-δ15N values for estimating trophic levels. We conducted a 98-d feeding experiment using five formulated isoenergetic feeds prepared with a high-quality protein source to evaluate the effect of protein quantity and quality on TEFs of liver and muscle tissues of juvenile Pacific yellowtail (Seriola lalandi), a carnivorous fish species. We decreased protein digestibility using well-established protocols that do not change AA profiles. Growth rates were higher in diets with higher protein content, and isotopic equilibrium was reached for both fish tissues and all treatments. Protein quantity and quality influenced isotope discrimination depending on tissue type and AA. In liver tissue, bulk TEFs showed a limited but significant relationship with protein quality, but did not differ with protein quantity or quality in muscle. None of the pre-established source AAs (Lys, Met, Phe, and Gly) TEFs varied significantly with protein quantity or quality in liver tissue. However, in muscle tissue, TEFPhe increased significantly with protein content and decreased in response to reduced digestibility, indicating it may not serve as proxy for baseline isotopic values used to calculate trophic level. Among trophic AAs, TEFLeu decreased significantly with increasing protein quantity in liver tissue, while both Leu and Ile TEFs decreased with lower protein digestibility in muscle tissue. Our results indicate that CSIA-AA in liver tissue provides more robust source and trophic AA-δ15N values than in muscle.

Protein retention assessment of four levels of poultry by-product substitution of fishmeal in rainbow trout (Oncorhynchus mykiss) diets using stable isotopes of nitrogen (δ15N) as natural tracers.

This is second part from an experiment where the nitrogen retention of poultry by-product meal (PBM) compared to fishmeal (FM) was evaluated using traditional indices. Here a quantitative method using stable isotope ratios of nitrogen (δ(15)N values) as natural tracers of nitrogen incorporation into fish biomass is assessed. Juvenile rainbow trout (Oncorhynchus mykiss) were fed for 80 days on isotopically distinct diets in which 0, 33, 66 and 100% of FM as main protein source was replaced by PBM. The diets were isonitrogenous, isolipidic and similar in gross energy content. Fish in all treatments reached isotopic equilibrium by the end of the experiment. Two-source isotope mixing models that incorporated the isotopic composition of FM and PBM as well as that of formulated feeds, empirically derived trophic discrimination factors and the isotopic composition of fish that had reached isotopic equilibrium to the diets were used to obtain a quantitative estimate of the retention of each source of nitrogen. Fish fed the diets with 33 and 66% replacement of FM by PBM retained poultry by-product meal roughly in proportion to its level of inclusion in the diets, whereas no differences were detected in the protein efficiency ratio. Coupled with the similar biomass gain of fishes fed the different diets, our results support the inclusion of PBM as replacement for fishmeal in aquaculture feeds. A re-feeding experiment in which all fish were fed a diet of 100% FM for 28 days indicated isotopic turnover occurred very fast, providing further support for the potential of isotopic ratios as tracers of the retention of specific protein sources into fish tissues. Stable isotope analysis is a useful tool for studies that seek to obtain quantitative estimates of the retention of different protein sources.