Sex-specific diet differences in harbor seals (Phoca vitulina) via spatial assortment.

The lack of recovery of Chinook salmon (Oncorhynchus tshawytscha) in the Pacific Northwest has been blamed in part on predation by pinnipeds, particularly the harbor seal (Phoca vitulina). Previous work at a limited number of locations has shown that male seal diet contains more salmon than that of female seals and that sex ratios at haul-out sites differ spatiotemporally. This intrapopulation variation in predation may result in greater effects on salmon than suggested by models assuming equal spatial distribution and diet proportion. To address the generality of these patterns, we examined the sex ratios and diet of male and female harbor seals from 13 haul-out sites in the inland waters of Washington State and the province of British Columbia during 2012-2018. DNA metabarcoding was conducted to determine prey species proportions of individual scat samples. The sex of harbor seals was then determined from each scat matrix sample with the use of quantitative polymerase chain reaction (qPCR). We analyzed 2405 harbor seal scat samples using generalized linear mixed models (GLMMs) to examine the factors influencing harbor seal sex ratio at haul-out sites and permutational multivariate analysis of variance (PERMANOVA) to examine the influence of sex and haul-out site on harbor seal diet composition. We found that the overall sex ratio was 1:1.02 (female:male) with notable spatiotemporal variation. Salmoniformes were about 2.6 times more abundant in the diet of males than in the diet of females, and Chinook salmon comprised ca. three times more of the average male harbor seal's diet than the average female's diet. Based on site-specific sex ratios and diet data, we identified three haul-out sites where Chinook salmon appear to be under high predation pressure by male harbor seals: Cowichan Bay, Cutts Area, and Fraser River. Our study indicates that combining sex-specific pinniped diet data with the sex ratio of haul-out sites can help identify priority sites of conservation concern.

Time-resolved ion energy measurements using a retarding potential analyzer for electric propulsion applications.

To completely characterize the evolving state of a plasma, diagnostic tools that enable measurements of the time-resolved behavior are required. In this study, a gridded ion source with superimposed oscillations was utilized to verify the functionality of a high-speed retarding potential analyzer (HSRPA), at frequencies equivalent to the low frequency oscillations occurring in Hall effect thrusters (HETs). The verification of this device provides an effective alternative to existing diagnostics for measuring time-resolved ion energies. Retarding potential analyzers (RPAs) have established themselves as a fundamental diagnostic in the field of electric propulsion (EP), enabling the measurement of ion energy distributions within the plumes of EP thrusters. The work presented here has demonstrated the capability of a standard RPA in conjunction with high-speed circuitry and data fusion techniques to produce time-resolved ion energy distribution functions (IEDFs) at higher frequencies and beam potentials than have previously been investigated. Tested frequencies ranged between 20 and 80 kHz with 10 V peak-to peak oscillations at a mean beam potential of 570 V. In addition, measurements were conducted with several waveforms, functioning as the superimposed oscillation, including a sine wave, triangle wave, and noisy sine wave. Data from the HSRPA were successfully reconstructed into time series utilizing two data fusion techniques: the empirical transfer function method and shadow manifold interpolation. Time-resolved IEDFs were produced at all frequency set points and waveforms. This investigation has demonstrated the HSRPA effectiveness at producing time-resolved measurements under conditions similar to those occurring in HETs.

Net overboard: Comparing marine eDNA sampling methodologies at sea to unravel marine biodiversity.

Environmental DNA (eDNA) analyses are powerful for describing marine biodiversity but must be optimized for their effective use in routine monitoring. To maximize eDNA detection probabilities of sparsely distributed populations, water samples are usually concentrated from larger volumes and filtered using fine-pore membranes, often a significant cost-time bottleneck in the workflow. This study aimed to streamline eDNA sampling by investigating plankton net versus bucket sampling, direct versus sequential filtration including self-preserving filters. Biodiversity was assessed using metabarcoding of the small ribosomal subunit (18S rRNA) and mitochondrial cytochrome c oxidase I (COI) genes. Multispecies detection probabilities were estimated for each workflow using a probabilistic occupancy modelling approach. Significant workflow-related differences in biodiversity metrics were reported. Highest amplicon sequence variant (ASV) richness was attained by the bucket sampling combined with self-preserving filters, comprising a large portion of microplankton. Less diversity but more metazoan taxa were captured in the net samples combined with 5 µm pore size filters. Prefiltered 1.2 µm samples yielded few or no unique ASVs. The highest average (~32%) metazoan detection probabilities in the 5 µm pore size net samples confirmed the effectiveness of preconcentration plankton for biodiversity screening. These results contribute to streamlining eDNA sampling protocols for uptake and implementation in marine biodiversity research and surveillance.

Data on the diets of Salish Sea harbour seals from DNA metabarcoding.

Marine trophic ecology data are in high demand as natural resource agencies increasingly adopt ecosystem-based management strategies that account for complex species interactions. Harbour seal (Phoca vitulina) diet data are of particular interest because the species is an abundant predator in the northeast Pacific Ocean and Salish Sea ecosystem that consumes Pacific salmon (Oncorhynchus spp.). A multi-agency effort was therefore undertaken to produce harbour seal diet data on an ecosystem scale using, 1) a standardized set of scat collection and analysis methods, and 2) a newly developed DNA metabarcoding diet analysis technique designed to identify prey species and quantify their relative proportions in seal diets. The DNA-based dataset described herein contains records from 4,625 harbour seal scats representing 52 haulout sites, 7 years, 12 calendar months, and a total of 11,641 prey identifications. Prey morphological hard parts analyses were conducted alongside, resulting in corresponding hard parts data for 92% of the scat DNA samples. A custom-built prey DNA sequence database containing 201 species (192 fishes, 9 cephalopods) is also provided.

Using in-situ environmental DNA sampling to detect the invasive New Zealand Mud Snail (Potamopyrgus antipodarum) in freshwaters.

Environmental DNA (eDNA) detection of aquatic invasive species is currently at the forefront of aquatic conservation efforts because the methodology provides a cost effective and sensitive means to detect animals at low densities. Developments in eDNA technologies have improved detection probabilities for rare, indicator, and invasive species over the past decade. However, standard lab analysis can take days or weeks before results are available and is prohibitive when rapid management decisions are required for mitigation. Here, we investigated the performance of a real-time quantitative PCR system for on-site eDNA detection of New Zealand mud snails (Potamopyrgus antipodarum). Six sites in western Washington, USA were sampled using the rapid eDNA technique and traditional methods, with five samples per site. On-site eDNA detection of mud snails resulted in a 10% increase in positive sites (16/30 = 53% positive) relative to visual surveys (13/30 = 43% positive). In addition, positive associations were observed between mud snail eDNA concentration (eDNA copies per reaction) and the number of mud snail individuals at each site (R 2 = 0.78). We show that the rapid on-site eDNA technology can be effective for detection and quantification of New Zealand mud snails in freshwaters. This on-site eDNA detection approach could possibly be used to initiate management protocols that allow for more rapid responses during the onset of biological invasions.

Large-scale molecular barcoding of prey DNA reveals predictors of intrapopulation feeding diversity in a marine predator.

Predator-prey interactions are critical in understanding how communities function. However, we need to describe intraspecific variation in diet to accurately depict those interactions. Harbor seals (Phoca vitulina) are an abundant marine predator that prey on species of conservation concern. We estimated intrapopulation feeding diversity (variation in feeding habits between individuals of the same species) of harbor seals in the Salish Sea. Estimates of feeding diversity were examined relative to sex, month, and location using a novel approach that combined molecular techniques, repeated cross-sectional sampling of scat, and a specialization metric (within-individual consistency in diet measured by the Proportional Similarity Index ( P S i )). Based on 1,083 scat samples collected from five haul-out sites during four nonsequential years, we quantified diet using metabarcoding techniques and determined the sex of the scat depositor using a molecular assay. Results suggest that intrapopulation feeding diversity was present. Specialization was high over short periods (24-48 hr, P S i  = 0.392, 95% CI = 0.013, R = 100,000) and variable in time and space. Females showed more specialization than males, particularly during summer and fall. Additionally, demersal and benthic prey species were correlated with more specialized diets. The latter finding suggests that this type of prey likely requires specific foraging strategies and that there are trade-offs between pelagic and benthic foraging styles for harbor seals. This differential feeding on prey species, as well as between sexes of harbor seals, indicates that predator-prey interactions in harbor seals are complex and that each sex may have a different impact on species of conservation concern. As such, describing intrapopulation feeding diversity may unravel hitherto unknown complex predator-prey interactions in the community.

Counting with DNA in metabarcoding studies: How should we convert sequence reads to dietary data?

Advances in DNA sequencing technology have revolutionized the field of molecular analysis of trophic interactions, and it is now possible to recover counts of food DNA sequences from a wide range of dietary samples. But what do these counts mean? To obtain an accurate estimate of a consumer's diet should we work strictly with data sets summarizing frequency of occurrence of different food taxa, or is it possible to use relative number of sequences? Both approaches are applied to obtain semi-quantitative diet summaries, but occurrence data are often promoted as a more conservative and reliable option due to taxa-specific biases in recovery of sequences. We explore representative dietary metabarcoding data sets and point out that diet summaries based on occurrence data often overestimate the importance of food consumed in small quantities (potentially including low-level contaminants) and are sensitive to the count threshold used to define an occurrence. Our simulations indicate that using relative read abundance (RRA) information often provides a more accurate view of population-level diet even with moderate recovery biases incorporated; however, RRA summaries are sensitive to recovery biases impacting common diet taxa. Both approaches are more accurate when the mean number of food taxa in samples is small. The ideas presented here highlight the need to consider all sources of bias and to justify the methods used to interpret count data in dietary metabarcoding studies. We encourage researchers to continue addressing methodological challenges and acknowledge unanswered questions to help spur future investigations in this rapidly developing area of research.

Large-scale molecular diet analysis in a generalist marine mammal reveals male preference for prey of conservation concern.

Sex-specific diet information is important in the determination of predator impacts on prey populations. Unfortunately, the diet of males and females can be difficult to describe, particularly when they are marine predators. We combined two molecular techniques to describe haul-out use and prey preferences of male and female harbor seals (Phoca vitulina) from Comox and Cowichan Bay (Canada) during 2012-2013. DNA metabarcoding quantified the diet proportions comprised of prey species in harbor seal scat, and qPCR determined the sex of the individual that deposited each scat. Using 287 female and 260 male samples, we compared the monthly sex ratio with GLMs and analyzed prey consumption relative to sex, season, site, and year with PERMANOVA. The sex ratio between monthly samples differed widely in both years (range = 12%-79% males) and showed different patterns at each haul-out site. Male and female diet differed across both years and sites: Females consumed a high proportion of demersal fish species while males consumed more salmonid species. Diet composition was related to both sex and season (PERMANOVA: R 2 = 27%, p < 0.001; R 2 = 24%, p < 0.001, respectively) and their interaction (PERMANOVA: R 2 = 11%, p < 0.001). Diet differences between males and females were consistent across site and year, suggesting fundamental foraging differences, including that males may have a larger impact on salmonids than females. Our novel combination of techniques allowed for both prey taxonomic and spatiotemporal resolution unprecedented in marine predators.

Rapid Detection and Monitoring of Flavobacterium psychrophilum in Water by Using a Handheld, Field-Portable Quantitative PCR System.

Advances in technology are making it easier for rapid field detection of microbes in aquaculture. Specifically, real-time quantitative PCR (qPCR) analysis, which has traditionally been confined to laboratory-based protocols, is now available in a handheld, field-portable system. The feasibility of using the Biomeme handheld qPCR system for rapid (<50 min) on-site detection and monitoring of Flavobacterium psychrophilum from filtered water samples was evaluated. Paired water samples were collected over a 23-d period from microcosm tanks that housed fish injected with known levels of F. psychrophilum. Water samples were filtered through 0.45-µm nitrocellulose filters and were analyzed with both the Biomeme qPCR platform and a traditional bench qPCR protocol. The two methods identified similar fluctuations in F. psychrophilum DNA throughout the study. Standard curves relating quantification cycles to the number of F. psychrophilum colony-forming units (CFU) were constructed and analyzed; results indicated that CFU increased rapidly between days 6 and 8 of the trial and then progressively decreased during the remaining 15 d. Average calculated log10 (CFU/mL) values were significantly correlated between the two platforms. Rapid, field-based qPCR can be incorporated into daily water quality monitoring protocols to help detect and monitor microbes in aquaculture systems.

Application of real-time quantitative PCR assays for detecting marine Brucella spp. in fish.

Brucella ceti and Brucella pinnipedialis have been documented as occurring in marine mammals, and B. ceti has been identified in 3 naturally acquired human cases. Seroconversion and infection patterns in Pacific Northwest harbor seals ( Phoca vitulina richardii) and North Atlantic hooded seals ( Cystophora cristata) indicate post-weaning exposure through prey consumption or lungworm infection, suggesting fish and possibly invertebrates play an epizootiologic role in marine Brucella transmission and possible foodborne risk to humans. We determined if real-time quantitative PCR (qPCR) assays can detect marine Brucella DNA in fish DNA. Insertion sequence (IS) 711 gene and sequence type (ST)27 primer-probe sets were used to detect Brucella associated with marine mammals and human zoonotic infections, respectively. First, DNA extracts from paired-species fish (containing 2 species) samples were tested and determined to be Brucella DNA negative using both IS 711 and ST27 primer-probe sets. A representative paired-species fish DNA sample was spiked with decreasing concentrations of B. pinnipedialis DNA to verify Brucella detection by the IS 711 primer-probe within fish DNA. A standard curve, developed using isolated DNA from B. pinnipedialis, determined the limit of detection. Finally, the IS 711 primer-probe was used to test Atlantic cod ( Gadus morhua) DNA extracts experimentally infected with the B. pinnipedialis hooded seal strain. In culture-positive cod tissue, the IS 711 limit of detection was ~1 genome copy of Brucella. Agreement between culture and PCR results for the 9 positive and 9 negative cod tissues was 100%. Although a larger sample set is required for validation, our study shows that qPCR can detect marine Brucella in fish.

Competing tradeoffs between increasing marine mammal predation and fisheries harvest of Chinook salmon.

Many marine mammal predators, particularly pinnipeds, have increased in abundance in recent decades, generating new challenges for balancing human uses with recovery goals via ecosystem-based management. We used a spatio-temporal bioenergetics model of the Northeast Pacific Ocean to quantify how predation by three species of pinnipeds and killer whales (Orcinus orca) on Chinook salmon (Oncorhynchus tshawytscha) has changed since the 1970s along the west coast of North America, and compare these estimates to salmon fisheries. We find that from 1975 to 2015, biomass of Chinook salmon consumed by pinnipeds and killer whales increased from 6,100 to 15,200 metric tons (from 5 to 31.5 million individual salmon). Though there is variation across the regions in our model, overall, killer whales consume the largest biomass of Chinook salmon, but harbor seals (Phoca vitulina) consume the largest number of individuals. The decrease in adult Chinook salmon harvest from 1975-2015 was 16,400 to 9,600 metric tons. Thus, Chinook salmon removals (harvest + consumption) increased in the past 40 years despite catch reductions by fisheries, due to consumption by recovering pinnipeds and endangered killer whales. Long-term management strategies for Chinook salmon will need to consider potential conflicts between rebounding predators or endangered predators and prey.

Quantitative DNA metabarcoding: improved estimates of species proportional biomass using correction factors derived from control material.

DNA metabarcoding is a powerful new tool allowing characterization of species assemblages using high-throughput amplicon sequencing. The utility of DNA metabarcoding for quantifying relative species abundances is currently limited by both biological and technical biases which influence sequence read counts. We tested the idea of sequencing 50/50 mixtures of target species and a control species in order to generate relative correction factors (RCFs) that account for multiple sources of bias and are applicable to field studies. RCFs will be most effective if they are not affected by input mass ratio or co-occurring species. In a model experiment involving three target fish species and a fixed control, we found RCFs did vary with input ratio but in a consistent fashion, and that 50/50 RCFs applied to DNA sequence counts from various mixtures of the target species still greatly improved relative abundance estimates (e.g. average per species error of 19 ± 8% for uncorrected vs. 3 ± 1% for corrected estimates). To demonstrate the use of correction factors in a field setting, we calculated 50/50 RCFs for 18 harbour seal (Phoca vitulina) prey species (RCFs ranging from 0.68 to 3.68). Applying these corrections to field-collected seal scats affected species percentages from individual samples (Δ 6.7 ± 6.6%) more than population-level species estimates (Δ 1.7 ± 1.2%). Our results indicate that the 50/50 RCF approach is an effective tool for evaluating and correcting biases in DNA metabarcoding studies. The decision to apply correction factors will be influenced by the feasibility of creating tissue mixtures for the target species, and the level of accuracy needed to meet research objectives.

Improving accuracy of DNA diet estimates using food tissue control materials and an evaluation of proxies for digestion bias.

Ecologists are increasingly interested in quantifying consumer diets based on food DNA in dietary samples and high-throughput sequencing of marker genes. It is tempting to assume that food DNA sequence proportions recovered from diet samples are representative of consumer's diet proportions, despite the fact that captive feeding studies do not support that assumption. Here, we examine the idea of sequencing control materials of known composition along with dietary samples in order to correct for technical biases introduced during amplicon sequencing and biological biases such as variable gene copy number. Using the Ion Torrent PGM(©) , we sequenced prey DNA amplified from scats of captive harbour seals (Phoca vitulina) fed a constant diet including three fish species in known proportions. Alongside, we sequenced a prey tissue mix matching the seals' diet to generate tissue correction factors (TCFs). TCFs improved the diet estimates (based on sequence proportions) for all species and reduced the average estimate error from 28 ± 15% (uncorrected) to 14 ± 9% (TCF-corrected). The experimental design also allowed us to infer the magnitude of prey-specific digestion biases and calculate digestion correction factors (DCFs). The DCFs were compared with possible proxies for differential digestion (e.g. fish protein%, fish lipid%) revealing a strong relationship between the DCFs and percent lipid of the fish prey, suggesting prey-specific corrections based on lipid content would produce accurate diet estimates in this study system. These findings demonstrate the value of parallel sequencing of food tissue mixtures in diet studies and offer new directions for future research in quantitative DNA diet analysis.

Quantifying sequence proportions in a DNA-based diet study using Ion Torrent amplicon sequencing: which counts count?

A goal of many environmental DNA barcoding studies is to infer quantitative information about relative abundances of different taxa based on sequence read proportions generated by high-throughput sequencing. However, potential biases associated with this approach are only beginning to be examined. We sequenced DNA amplified from faeces (scats) of captive harbour seals (Phoca vitulina) to investigate whether sequence counts could be used to quantify the seals' diet. Seals were fed fish in fixed proportions, a chordate-specific mitochondrial 16S marker was amplified from scat DNA and amplicons sequenced using an Ion Torrent PGM™. For a given set of bioinformatic parameters, there was generally low variability between scat samples in proportions of prey species sequences recovered. However, proportions varied substantially depending on sequencing direction, level of quality filtering (due to differences in sequence quality between species) and minimum read length considered. Short primer tags used to identify individual samples also influenced species proportions. In addition, there were complex interactions between factors; for example, the effect of quality filtering was influenced by the primer tag and sequencing direction. Resequencing of a subset of samples revealed some, but not all, biases were consistent between runs. Less stringent data filtering (based on quality scores or read length) generally produced more consistent proportional data, but overall proportions of sequences were very different than dietary mass proportions, indicating additional technical or biological biases are present. Our findings highlight that quantitative interpretations of sequence proportions generated via high-throughput sequencing will require careful experimental design and thoughtful data analysis.

Great shearwater (Puffinus gravis) mortality events along the eastern coast of the United States.

The Great Shearwater (Puffinus gravis) is an abundant pelagic seabird that undertakes transequatorial migrations between the North and South Atlantic Ocean. This species is a useful indicator of large-scale alterations in marine dynamics due to its wide geographic range, long-distance migrations, and relative abundance. From 1993 to 2011, 12 separate mortality events, with 4,961 Great Shearwaters recovered, were documented along the eastern coast of the United States. Of these, seven events (n=4,885) occurred in the Southeast (SE) and five (n=76) in the Northeast (NE) United States. The cause of death was determined either by necropsy (n=60) or external examination (n=4,901). All Great Shearwaters stranded along the SE United States were emaciated while 58% were emaciated in the NE United States. No plastic was observed in Great Shearwaters in the SE US (n=27), but the gastrointestinal tract of 82% (n=27) of all stranded birds along the NE United States had at least one plastic bead. There was no evidence of infectious disease or heavy metals in stranded Great Shearwaters examined (n=14, from the 2005 SE event). Stable isotope analysis of feathers (n=9, from a 2007 SE event) suggests dietary differences between emaciated stranded birds and live-caught healthy birds. The temporal distribution of stranding detections suggests a general increase in the number of observed Great Shearwater strandings over the past two decades. From 1993 to 2000 there were a total of three mortality events with 296 individual Great Shearwaters. However, there was a threefold increase in the number of mortality events from 2001 to 2011 (nine events involving 4,665 individuals). The causes of this apparent increase in strandings are unknown but may be due to an increase in reporting effort over the past two decades combined with changing oceanographic conditions in the South Atlantic Ocean, leading to large-scale mortality of emaciated Great Shearwaters along the east coast of the United States.

Baseline health parameters and species comparisons among free-ranging Atlantic sharpnose (Rhizoprionodon terraenovae), bonnethead (Sphyrna tiburo), and spiny dogfish (Squalus acanthias) sharks in Georgia, Florida, and Washington, USA.

Sharks are of commercial, research, conservation, and exhibition importance but we know little regarding health parameters and population status for many species. Here we present health indicators and species comparisons for adults of three common wild-caught species: 30 Atlantic sharpnose sharks (Rhizoprionodon terraenovae) and 31 bonnethead sharks (Sphyrna tiburo) from the western Atlantic, and 30 spiny dogfish sharks (Squalus acanthias) from the eastern Pacific. All animals were captured during June-July 2009 and 2010. Median values and preliminary reference intervals were calculated for hematology, plasma biochemistry, trace nutrients, and vitamin A, E, and D concentrations. Significant differences, attributable to physiologic differences among the species, were found in the basic hematologic and plasma biochemistry variables. Significant species differences in arsenic and selenium plasma concentrations were found and appear to coincide with diet and habitat variability among these three species. Vitamin E was significantly higher in the bonnethead shark, again related to the foraging ecology and ingestion of plant material by this species. The Atlantic sharpnose had significantly higher vitamin A concentrations, supported by the higher proportion of teleosts in the diet. Vitamin D was below the limit of quantification in all three species. These preliminary reference intervals for health variables can be used to assess and monitor the population health and serve as indicators of nutritional status in these populations of wild elasmobranchs.