Contrasting effect of hybridization on genetic differentiation in three rockfish species with similar life history.

Hybridization can provide evolutionary benefits (e.g., population resilience to climate change) through the introduction of adaptive alleles and increase of genetic diversity. Nevertheless, management strategies may be designed based only on the parental species within a hybrid zone, without considering the hybrids. This can lead to ineffective spatial management of species, which can directly harm population diversity and negatively impact food webs. Three species of rockfish (Brown Rockfish (Sebastes caurinus), Copper Rockfish (S. auriculatus), and Quillback Rockfish (S. maliger)) are known to hybridize within Puget Sound, Washington, but genetic data from these species are used to infer population structure in the entire genus, including in species that do not hybridize. The goal of this project was to estimate the hybridization rates within the region and determine the effect of hybridization on geographic patterns of genetic structure. We sequenced 290 Brown, Copper, and Quillback rockfish using restriction-site associated DNA sequencing (RADseq) from four regions within and outside Puget Sound, Washington. We show that (i) hybridization within Puget Sound was asymmetrical, not recent, widespread among individuals, and relatively low level within the genome, (ii) hybridization affected population structure in Copper and Brown rockfish, but not in Quillback Rockfish and (iii) after taking hybridization into account we found limited directional dispersal in Brown and Copper rockfish, and evidence for two isolated populations in Quillback Rockfish. Our results suggest that rockfish population structure is species-specific, dependent on the extent of hybridization, and cannot be inferred from one species to another despite similar life history.

Five decades of change in somatic growth of Pacific hake from Puget Sound and Strait of Georgia.

Declines in fish body size have been reported in many populations and these changes likely have important ramifications for the sustainability of harvested species and ecosystem function. Pacific hake, Merluccius productus, have shown declines in size over the last several decades for populations located in Puget Sound (PS), Washington, USA, and Strait of Georgia (SoG), British Columbia, Canada. To examine this decrease in size, we used archived otoliths from both populations to assess when the decrease in somatic growth occurred and explored what factors and processes might explain the decline, including otolith microchemistry to infer the environment experienced by fish at different ages. Results indicated that substantial changes in juvenile somatic growth have occurred across decades. The divergence in body size occurred in the second summer, whereby SoG fish grew, on average, 18% more than PS fish. Within the PS population, somatic growth differed significantly among fish that hatched in the 1980s, 1990s, and 2010s, such that the more recently hatched fish grew 26% more in their first summer and 71% less in their second summer relative to those that hatched in the 1980s. In comparison, growth of SoG fish did not differ between those that hatched in 1970s and 1990s. For both populations growth in the first and third summer was positively and negatively related, respectively, to the abundance of harbor seals, while growth in the first and second summer was negatively related to salinity. Overall, this study highlights the complicated nature of Pacific hake population recovery under dynamic, and typically uncontrollable, variation in biotic and abiotic conditions.

Population assessment using multivariate time-series analysis: A case study of rockfishes in Puget Sound.

Estimating a population's growth rate and year-to-year variance is a key component of population viability analysis (PVA). However, standard PVA methods require time series of counts obtained using consistent survey methods over many years. In addition, it can be difficult to separate observation and process variance, which is critical for PVA. Time-series analysis performed with multivariate autoregressive state-space (MARSS) models is a flexible statistical framework that allows one to address many of these limitations. MARSS models allow one to combine surveys with different gears and across different sites for estimation of PVA parameters, and to implement replication, which reduces the variance-separation problem and maximizes informational input for mean trend estimation. Even data that are fragmented with unknown error levels can be accommodated. We present a practical case study that illustrates MARSS analysis steps: data choice, model set-up, model selection, and parameter estimation. Our case study is an analysis of the long-term trends of rockfish in Puget Sound, Washington, based on citizen science scuba surveys, a fishery-independent trawl survey, and recreational fishery surveys affected by bag-limit reductions. The best-supported models indicated that the recreational and trawl surveys tracked different, temporally independent assemblages that declined at similar rates (an average of -3.8% to -3.9% per year). The scuba survey tracked a separate increasing and temporally independent assemblage (an average of 4.1% per year). Three rockfishes (bocaccio, canary, and yelloweye) are listed in Puget Sound under the US Endangered Species Act (ESA). These species are associated with deep water, which the recreational and trawl surveys sample better than the scuba survey. All three ESA-listed rockfishes declined as a proportion of recreational catch between the 1970s and 2010s, suggesting that they experienced similar or more severe reductions in abundance than the 3.8-3.9% per year declines that were estimated for rockfish populations sampled by the recreational and trawl surveys.