Individual variation in life-history timing: synchronous presence, asynchronous events and phenological compensation in a wild mammal

Many animals and plants have species-typical annual cycles, but individuals vary in their timing of life-history events. Individual variation in fur replacement (moult) timing is poorly understood in mammals due to the challenge of repeated observations and longitudinal sampling. We examined factors that influence variation in moult duration and timing among elephant seals (Mirounga angustirostris). We quantified the onset and progression of fur loss in 1178 individuals. We found that an exceptionally rapid visible moult (7 days, the shortest of any mammals or birds), and a wide range of moult start dates (spanning 6-10× the event duration) facilitated high asynchrony across individuals (only 20% of individuals in the population moulting at the same time). Some of the variation was due to reproductive state, as reproductively mature females that skipped a breeding season moulted a week earlier than reproductive females. Moreover, individual variation in timing and duration within age-sex categories far outweighed (76-80%) variation among age-sex categories. Individuals arriving at the end of the moult season spent 50% less time on the beach, which allowed them to catch up in their annual cycles and reduce population-level variance during breeding. These findings underscore the importance of individual variation in annual cycles.

Cross-shore transport and eddies promote large scale response to urban eutrophication

A key control on the magnitude of coastal eutrophication is the degree to which currents quickly transport nitrogen derived from human sources away from the coast to the open ocean before eutrophication develops. In the Southern California Bight (SCB), an upwelling-dominated eastern boundary current ecosystem, anthropogenic nitrogen inputs increase algal productivity and cause subsurface acidification and oxygen (O2 ) loss along the coast. However, the extent of anthropogenic influence on eutrophication beyond the coastal band, and the physical transport mechanisms and biogeochemical processes responsible for these effects are still poorly understood. Here, we use a submesoscale-resolving numerical model to document the detailed biogeochemical mass balance of nitrogen, carbon and oxygen, their physical transport, and effects on offshore habitats. Despite management of terrestrial nutrients that has occurred in the region over the last 20 years, coastal eutrophication continues to persist. The input of anthropogenic nutrients promote an increase in productivity, remineralization and respiration offshore, with recurrent O2 loss and pH decline in a region located 30-90 km from the mainland. During 2013 to 2017, the spatially averaged 5-year loss rate across the Bight was 1.3 mmol m- 3 O2 , with some locations losing on average up to 14.2 mmol m- 3 O2 . The magnitude of loss is greater than model uncertainty assessed from data-model comparisons and from quantification of intrinsic variability. This phenomenon persists for 4 to 6 months of the year over an area of 278,40 km2 ( ∼ 30% of SCB area). These recurrent features of acidification and oxygen loss are associated with cross-shore transport of nutrients by eddies and plankton biomass and their accumulation and retention within persistent eddies offshore within the SCB.

Genomic resources for the Yellowfin tuna Thunnus albacares

BACKGROUND

The Yellowfin tuna (Thunnus albacares) is a large tuna exploited by major fisheries in tropical and subtropical waters of all oceans except the Mediterranean Sea. Genomic studies of population structure, adaptive variation or of the genetic basis of phenotypic traits are needed to inform fisheries management but are currently limited by the lack of a reference genome for this species. Here we report a draft genome assembly and a linkage map for use in genomic studies of T. albacares.

METHODS AND RESULTS

Illumina and PacBio SMRT sequencing were used in combination to generate a hybrid assembly that comprises 743,073,847 base pairs contained in 2,661 scaffolds. The assembly has a N50 of 351,587 and complete and partial BUSCO scores of 86.47% and 3.63%, respectively. Double-digest restriction associated DNA (ddRAD) was used to genotype the 2 parents and 164 of their F1 offspring resulting from a controlled breeding cross, retaining 19,469 biallelic single nucleotide polymorphism (SNP) loci. The SNP loci were used to construct a linkage map that features 24 linkage groups that represent the 24 chromosomes of yellowfin tuna. The male and female maps span 1,243.8 cM and 1,222.9 cM, respectively. The map was used to anchor the assembly in 24 super-scaffolds that contain 79% of the yellowfin tuna genome. Gene prediction identified 46,992 putative genes 20,203 of which could be annotated via gene ontology.

CONCLUSIONS

The draft reference will be valuable to interpret studies of genome wide variation in T. albacares and other Scombroid species.

Caribbean deepwater snappers: Application of the bomb radiocarbon age estimation validation in understanding aspects of ecology and life history

Reef fishes have been utilized as food fish throughout the U.S. Caribbean and Gulf of Mexico waters for centuries, with increasing fishing effort in recent decades. As a result, many species have experienced declines in landings, including deepwater snappers such as queen snapper Etelis oculatus and blackfin snapper Lutjanus buccanella. However, little to no peer-reviewed published research exists on basic life history parameters for either species. Confirming the accuracy of an age estimation method for a fish species is essential for ensuring sustainable fisheries management. This is because in the assessment of fisheries species population age-based parameters, including longevity, age at sexual maturity, growth rate, mortality, age-specific reproductive output, and lifetime reproductive output, are important in understanding overall life history strategies of managed stocks. The past stock assessment on U.S. Caribbean queen snapper utilized an estimated longevity of 8 y, derived from length frequencies for fish from St. Lucia. Blackfin snapper has an estimated longevity of 27 y based on a relatively small study from offshore waters of the southeastern U.S. The focus of our investigation was to estimate maximum longevity of two data-poor species in the U.S. Caribbean. The accuracy of ageing methods was tested via bomb radiocarbon age estimation validation and effects of depth on Δ14C in otolith cores and eye lens core values were examined. Results from our work indicate a maximum validated age of 45 y for queen snapper, and 43 y for blackfin snapper. Our findings indicate queen snapper and blackfin snapper are long-lived (> 40 y). The resulting Δ14C comparison between eye lens cores and otolith cores has important implications for the study of age validation, specifically when deepwater species are involved.

Genetic diversity and demographic history of the leopard seal: A Southern Ocean top predator

Leopard seals (Hydrurga leptonyx) are top predators that can exert substantial top-down control of their Antarctic prey species. However, population trends and genetic diversity of leopard seals remain understudied, limiting our understanding of their ecological role. We investigated the genetic diversity, effective population size and demographic history of leopard seals to provide fundamental data that contextualizes their predatory influence on Antarctic ecosystems. Ninety leopard seals were sampled from the northern Antarctic Peninsula during the austral summers of 2008-2019 and a 405bp segment of the mitochondrial control region was sequenced for each individual. We uncovered moderate levels of nucleotide (π = 0.013) and haplotype (Hd = 0.96) diversity, and the effective population size was estimated at around 24,000 individuals (NE = 24,376; 95% CI: 16,876-33,126). Consistent with findings from other ice-breeding pinnipeds, Bayesian skyline analysis also revealed evidence for population expansion during the last glacial maximum, suggesting that historical population growth may have been boosted by an increase in the abundance of sea ice. Although leopard seals can be found in warmer, sub-Antarctic locations, the species' core habitat is centered on the Antarctic, making it inherently vulnerable to the loss of sea ice habitat due to climate change. Therefore, detailed assessments of past and present leopard seal population trends are needed to inform policies for Antarctic ecosystems.

Eye lens-derived Δ14C signatures validate extreme longevity in the deepwater scorpaenid blackbelly rosefish (Helicolenus dactylopterus)

Many members of the scorpaenid subfamily: Sebastinae (rockfishes and their relatives) exhibit slow growth and extreme longevity (> 100 y), thus are estimated to be vulnerable to overfishing. Blackbelly rosefish (Helicolenus dactylopterus) is a deepwater sebastine whose longevity estimates range widely, possibly owing to different regional levels of fisheries exploitation across its Atlantic Ocean range. However, age estimation has not been validated for this species and ageing for sebastines in general is uncertain. We performed age validation of northern Gulf of Mexico blackbelly rosefish via an application of the bomb radiocarbon chronometer which utilized eye lens cores instead of more traditional otolith cores as the source of birth year Δ¹⁴C signatures. The correspondence of eye lens core Δ¹⁴C with a regional reference series was tested with a novel Bayesian spline analysis, which revealed otolith opaque zone counts provide accurate age estimates. Maximum observed longevity was 90 y, with 17.5% of individuals aged to be > 50 y. Bayesian growth analysis, with estimated length-at-birth included as a prior, revealed blackbelly rosefish exhibit extremely slow growth (k = 0.08 y^-1). Study results have important implications for the management of blackbelly rosefish stocks, as extreme longevity and slow growth imply low resilience to fishing pressure.

Trophic ontogeny of a generalist predator is conserved across space

Consumers can influence ecological patterns and processes through their trophic roles and contributions to the flow of energy through ecosystems. However, the diet and associated trophic roles of consumers commonly change during ontogeny. Despite the prevalence of ontogenetic variation in trophic roles of most animals, we lack an understanding of whether they change consistently across local populations and broad geographic gradients. We examined how the diet and trophic position of a generalist marine predator varied with ontogeny across seven broadly separated locations (~ 750 km). We observed a high degree of heterogeneity in prey consumed without evidence of spatial structuring in this variability. However, compound-specific isotope analysis of amino acids revealed remarkably consistent patterns of increasing trophic position through ontogeny across local populations, suggesting that the roles of this generalist predator scaled with its body size across space. Given the high degree of diet heterogeneity we observed, this finding suggests that even though the dietary patterns differed, the underlying food web architecture transcended variation in prey species across locations for this generalist consumer. Our research addresses a gap in empirical field work regarding the interplay between stage-structured populations and food webs, and suggests ontogenetic changes in trophic position can be consistent in generalist consumers.

Population structure and adaptive differentiation in the sea cucumber Apostichopus californicus and implications for spatial resource management

A growing body of evidence suggests that spatial population structure can develop in marine species despite large population sizes and high gene flow. Characterizing population structure is important for the effective management of exploited species, as it can be used to identify appropriate scales of management in fishery and aquaculture contexts. The California sea cucumber, Apostichopus californicus, is one such exploited species whose management could benefit from further characterization of population structure. Using restriction site-associated DNA (RAD) sequencing, we developed 2075 single nucleotide polymorphisms (SNPs) to quantify genetic structure over a broad section of the species' range along the North American west coast and within the Salish Sea, a region supporting the Washington State A. californicus fishery and developing aquaculture production of the species. We found evidence for population structure (global fixation index (FST) = 0.0068) with limited dispersal driving two patterns of differentiation: isolation-by-distance and a latitudinal gradient of differentiation. Notably, we found detectable population differences among collection sites within the Salish Sea (pairwise FST = 0.001-0.006). Using FST outlier detection and gene-environment association, we identified 10.2% of total SNPs as putatively adaptive. Environmental variables (e.g., temperature, salinity) from the sea surface were more correlated with genetic variation than those same variables measured near the benthos, suggesting that selection on pelagic larvae may drive adaptive differentiation to a greater degree than selection on adults. Our results were consistent with previous estimates of and patterns in population structure for this species in other extents of the range. Additionally, we found that patterns of neutral and adaptive differentiation co-varied, suggesting that adaptive barriers may limit dispersal. Our study provides guidance to decision-makers regarding the designation of management units for A. californicus and adds to the growing body of literature identifying genetic population differentiation in marine species despite large, nominally connected populations.

Bacterial microbiomes from mucus and breath of southern resident killer whales (Orcinus orca)

Opportunities to assess odontocete health are restricted due to their limited time at the surface, relatively quick movements and large geographic ranges. For endangered populations such as the southern resident killer whales (SKRWs) of the northeast Pacific Ocean, taking advantage of non-invasive samples such as expelled mucus and exhaled breath is appealing. Over the past 12 years, such samples were collected, providing a chance to analyse and assess their bacterial microbiomes using amplicon sequencing. Based on operational taxonomic units, microbiome communities from SRKW and transient killer whales showed little overlap between mucus, breath and seawater from SRKW habitats and six bacterial phyla were prominent in expelled mucus but not in seawater. Mollicutes and Fusobacteria were common and abundant in mucus, but not in breath or seawater, suggesting these bacterial classes may be normal constituents of the SRKW microbiome. Out of 134 bacterial families detected, 24 were unique to breath and mucus, including higher abundances of Burkholderiaceae, Moraxellaceae and Chitinophagaceae. Although there were multiple bacterial genera in breath or mucus that include pathogenic species (e.g. Campylobacter, Hemophilus, Treponema), the presence of these bacteria is not necessarily evidence of disease or infection. Future emphasis on genotyping mucus samples to the individual animal will allow further assessment in the context of that animal's history, including body condition index and prior contaminants burden. This study is the first to examine expelled mucus from cetaceans for microbiomes and demonstrates the value of analysing these types of non-invasive samples.

Estimating population size when individuals are asynchronous: A model illustrated with northern elephant seal breeding colonies

Our aim was to develop a method for estimating the number of animals using a single site in an asynchronous species, meaning that not all animals are present at once so that no one count captures the entire population. This is a common problem in seasonal breeders, and in northern elephant seals, we have a model for quantifying asynchrony at the Año Nuevo colony. Here we test the model at several additional colonies having many years of observations and demonstrate how it can account for animals not present on any one day. This leads to correction factors that yield total population from any single count throughout a season. At seven colonies in California for which we had many years of counts of northern elephant seals, we found that female arrival date varied < 2 days between years within sites and by < 5 days between sites. As a result, the correction factor for any one day was consistent, and at each colony, multiplying a female count between 26 and 30 Jan by 1.15 yielded an estimate of total population size that minimized error. This provides a method for estimating the female population size at colonies not yet studied. Our method can produce population estimates with minimal expenditure of time and resources and will be applicable to many seasonal species with asynchronous breeding phenology, particularly colonial birds and other pinnipeds. In elephant seals, it will facilitate monitoring the population over its entire range.

Trophic ecology of yellowtail rockfish (Sebastes flavidus) during a marine heat wave off central California, USA

The yellowtail rockfish, Sebastes flavidus, is a widespread and abundant mesopredator in the California Current Large Marine Ecosystem. We utilized stomach content and stable isotope analyses to investigate the trophic ecology of this species at three sites off central California just before (August–October 2013) and during (August and October 2014) a marine heat wave. Sebastes flavidus largely consumed pelagic prey (zooplankton and micronekton). Diets were dominated by tunicates (salps and pyrosomes), pelagic crustaceans (euphausiids, hyperid amphipods, larval decapods), and fishes, with the relative contribution of these prey taxa varying spatially (sample location, longitude, depth) and temporally (year, month), based on complementary multivariate analyses. Prey-specific indices demonstrated that individual S. flavidus diet composition typically was dominated by one of these prey groups, and that prey switching occurred based on the relative availability of prey and their energetic importance. Stable isotope analysis of δ¹⁵N indicated that the S. flavidus populations sampled in 2014 had been feeding at an elevated trophic position and more variable prey spectrum relative to 2013, probably as a consequence of greater piscivory and the incorporation of temporal changes in diet composition. Because its opportunistic feeding behavior reflects the dynamism and heterogeneity of the pelagic forage preyscape, S. flavidus may be an important ecosystem indicator species. For example, the novel incorporation of pyrosomes as a large portion of the diet of S. flavidus during 2013–2014 directly related to the massive increase in pyrosome abundance in the California Current during the 2014 marine heat wave.

Using claws to compare reproduction, stress and diet of female bearded and ringed seals in the Bering and Chukchi seas, Alaska, between 1953-1968 and 1998-2014

Rapid climate warming is decreasing sea ice thickness, extent and duration. Marine mammals such as bearded (Erignathus barbatus) and ringed (Pusa hispida) seals, which use sea ice for pupping, molting and resting, may be negatively affected. Claws from bearded and ringed seals store up to 14 and 12 years of sequential analyte data, respectively. These data can be used to compare reproduction, stress and diet across decades. In this study, we compare progesterone, cortisol and carbon and nitrogen stable isotopes in female bearded and ringed seals during 1953-1968 (pre-1968, a period prior to sea ice decline) to 1998-2014 (post-1998, a period during sea ice decline). When comparing these periods, bearded seals had statistically higher cortisol concentrations post-1998, and for both species δ13C was more negative post-1998, while progesterone and δ15N did not change. There was a positive relationship between progesterone and cortisol Z-scores for both species, except for ringed seals post-1998. There was a negative relationship between cortisol Z-scores and δ13C for bearded seals evident in post-1998 indicating that higher cortisol Z-scores are associated with more negative δ13C in bearded seals in recent years. This negative relationship between cortisol and δ13C in bearded seals suggests a shift to higher prey diversity, possibly due to changes in sea ice in the Pacific Arctic evident post 1998. Progesterone Z-scores corresponded to expected differences among non-pregnant, unimplanted, implanted and post-partum individuals. Using these data, pregnancy history was determined for reproductive years for each individual female sampled, which could allow for yearly pregnancy rates to be calculated given a large enough representative sample of the population. These results combine decades of observational studies with hormones and stable isotopes to infer changes in reproduction, stress and diet, as well as the connection between these life history parameters.

Predator type influences the frequency of functional responses to prey in marine habitats

The functional response of a consumer to a gradient of resource density is a widespread and consistent framework used to quantify the importance of consumption to population dynamics and stability. Within benthic marine ecosystems, both crustaceans and fishes can provide strong top-down pressure on prey populations. Taxon-specific differences in biomechanics or habitat use, among other factors, may lead to variable functional response forms or parameter values (attack rate, handling time). Based on a review of 189 individual functional response fits, we find that these predator guilds differ in their frequency distribution of functional response types, with crustaceans exhibiting nearly double the proportion of sigmoidal, density-dependent functional responses (Holling type III) as predatory fishes. The implications of this finding for prey population stability are significant because type III responses allow prey persistence while type II responses are de-stabilizing and can lead to extinction. Comparing per capita predation rates across diverse taxa can provide integrative insights into predatory effects and the ability of predation to drive community structure.

Validation of a novel method to create temporal records of hormone concentrations from the claws of ringed and bearded seals

Ringed (Pusa hispida) and bearded seals (Erignathus barbatus) inhabit vast and often remote areas in the Arctic, making it difficult to obtain long-term physiological information concerning health and reproduction. These seals are experiencing climate-driven changes in their habitat that could result in physiological stress. Chronic physiological stress can lead to immunosuppression, decreased reproduction and decreased growth. Recently, keratin has become a popular matrix to measure steroid hormones, such as stress-related cortisol and reproduction-related progesterone. We developed and validated methods to extract cortisol and progesterone from the claws of adult female ringed (n = 20) and bearded (n = 3) seals using enzyme immunosorbent assays. As ringed and bearded seal claws grow, a pair of dark- and light-colored bands of keratin is deposited annually providing a guide for sampling. Two processing methods were evaluated, removal of claw material with a grinding bit or grinding followed by mechanical pulverization (102 paired samples from six claws, two each from three seals). Adding the mechanical pulverization step resulted in a 1.5-fold increase in hormone extraction. Progesterone from the proximal claw band was evaluated to biologically validate claw material as a measure of pregnancy in ringed seals (n = 14). Claws from pregnant seals had significantly higher claw progesterone concentrations than from non-pregnant seals. This suggests that the elevated progesterone associated with gestation was reflected in the claws, and that the most proximal claw band was indicative of pregnancy status at time of death. Thus, although the sample size was low and the collection dates unbalanced, this study demonstrates the potential to use claws to examine an extended time series (up to 12 yrs) of cortisol and progesterone concentrations in ringed and bearded seal claws.

Migratory culture, population structure and stock identity in North Pacific beluga whales (Delphinapterus leucas)

The annual return of beluga whales, Delphinapterus leucas, to traditional seasonal locations across the Arctic may involve migratory culture, while the convergence of discrete summering aggregations on common wintering grounds may facilitate outbreeding. Natal philopatry and cultural inheritance, however, has been difficult to assess as earlier studies were of too short a duration, while genetic analyses of breeding patterns, especially across the beluga's Pacific range, have been hampered by inadequate sampling and sparse information on wintering areas. Using a much expanded sample and genetic marker set comprising 1,647 whales, spanning more than two decades and encompassing all major coastal summering aggregations in the Pacific Ocean, we found evolutionary-level divergence among three geographic regions: the Gulf of Alaska, the Bering-Chukchi-Beaufort Seas, and the Sea of Okhotsk (Φst = 0.11-0.32, Rst = 0.09-0.13), and likely demographic independence of (Fst-mtDNA = 0.02-0.66), and in many cases limited gene flow (Fst-nDNA = 0.0-0.02; K = 5-6) among, summering groups within regions. Assignment tests identified few immigrants within summering aggregations, linked migrating groups to specific summering areas, and found that some migratory corridors comprise whales from multiple subpopulations (PBAYES = 0.31:0.69). Further, dispersal is male-biased and substantial numbers of closely related whales congregate together at coastal summering areas. Stable patterns of heterogeneity between areas and consistently high proportions (~20%) of close kin (including parent-offspring) sampled up to 20 years apart within areas (G = 0.2-2.9, p>0.5) is the first direct evidence of natal philopatry to migration destinations in belugas. Using recent satellite telemetry findings on belugas we found that the spatial proximity of winter ranges has a greater influence on the degree of both individual and genetic exchange than summer ranges (rwinter-Fst-mtDNA = 0.9, rsummer-Fst-nDNA = 0.1). These findings indicate widespread natal philopatry to summering aggregation and entire migratory circuits, and provide compelling evidence that migratory culture and kinship helps maintain demographically discrete beluga stocks that can overlap in time and space.

Spatial Distribution of Reef Fish Species along the Southeast US Atlantic Coast Inferred from Underwater Video Survey Data

Marine fish abundance and distribution often varies across spatial scales for a variety of reasons, and this variability has significant ecological and management consequences. We quantified the distribution of reef-associated fish species along the southeast United States Atlantic coast using underwater video survey samples (N = 4,855 in 2011-2014) to elucidate variability within species across space, depths, and habitats, as well as describe broad-scale patterns in species richness. Thirty-two species were seen at least 10 times on video, and the most commonly observed species were red porgy (Pagrus pagrus; 41.4% of videos), gray triggerfish (Balistes capriscus; 31.0%), black sea bass (Centropristis striata; 29.1%), vermilion snapper (Rhomboplites aurorubens; 27.7%), and red snapper (Lutjanus campechanus; 22.6%). Using generalized additive models, we found that most species were non-randomly distributed across space, depths, and habitats. Most rare species were observed along the continental shelf break, except for goliath grouper (Epinephelus itajara), which was found on the continental shelf in Florida and Georgia. We also observed higher numbers of species in shelf-break habitats from southern North Carolina to Georgia, and fewer in shallower water and at the northern and southern ends of the southeast United States Atlantic coast. Our study provides the first broad-scale description of the spatial distribution of reef fish in the region to be based on fishery-independent data, reinforces the utility of underwater video to survey reef fish, and can help improve the management of reef fish in the SEUS, for example, by improving indices of abundance.