Robert S. Gregory

Parallel adaptive evolution of Atlantic cod on both sides of the Atlantic Ocean in response to temperature

Despite the enormous economic and ecological importance of marine organisms, the spatial scales of adaptation and biocomplexity remain largely unknown. Yet, the preservation of local stocks that possess adaptive diversity is critical to the long-term maintenance of productive stable fisheries and ecosystems. Here, we document genomic evidence of range-wide adaptive differentiation in a broadcast spawning marine fish, Atlantic cod (Gadus morhua), using a genome survey of single nucleotide polymorphisms. Of 1641 gene-associated polymorphisms examined, 70 (4.2%) tested positive for signatures of selection using a Bayesian approach. We identify a subset of these loci (n = 40) for which allele frequencies show parallel temperature-associated clines (p < 0.001, r2 = 0.89) in the eastern and western north Atlantic. Temperature associations were robust to the statistical removal of geographic distance or latitude effects, and contrasted ‘neutral’ loci, which displayed no temperature association. Allele frequencies at temperature-associated loci were significantly correlated, spanned three linkage groups and several were successfully annotated supporting the involvement of multiple independent genes. Our results are consistent with the evolution and/or selective sweep of multiple genes in response to ocean temperature, and support the possibility of a new conservation paradigm for non-model marine organisms based on genomic approaches to resolving functional and adaptive diversity.

Genomic islands of divergence and their consequences for the resolution of spatial structure in an exploited marine fish

As populations diverge, genomic regions associated with adaptation display elevated differentiation. These genomic islands of adaptive divergence can inform conservation efforts in exploited species, by refining the delineation of management units, and providing genomic tools for more precise and effective population monitoring and the successful assignment of individuals and products. We explored heterogeneity in genomic divergence and its impact on the resolution of spatial population structure in exploited populations of Atlantic cod, Gadus morhua, using genome wide expressed sequence derived single nucleotide polymorphisms in 466 individuals sampled across the range. Outlier tests identified elevated divergence at 5.2% of SNPs, consistent with directional selection in one‐third of linkage groups. Genomic regions of elevated divergence ranged in size from a single position to several cM. Structuring at neutral loci was associated with geographic features, whereas outlier SNPs revealed genetic discontinuities in both the eastern and western Atlantic. This fine‐scale geographic differentiation enhanced assignment to region of origin, and through the identification of adaptive diversity, fundamentally changes how these populations should be conserved. This work demonstrates the utility of genome scans for adaptive divergence in the delineation of stock structure, the traceability of individuals and products, and ultimately a role for population genomics in fisheries conservation.

Evaluating SNP ascertainment bias and its impact on population assignment in Atlantic cod, Gadus morhua.

The increasing use of single nucleotide polymorphisms (SNPs) in studies of nonmodel organisms accentuates the need to evaluate the influence of ascertainment bias on accurate ecological or evolutionary inference. Using a panel of 1641 expressed sequence tag–derived SNPs developed for northwest Atlantic cod (Gadus morhua), we examined the influence of ascertainment bias and its potential impact on assignment of individuals to populations ranging widely in origin. We hypothesized that reductions in assignment success would be associated with lower diversity in geographical regions outside the location of ascertainment. Individuals were genotyped from 13 locations spanning much of the contemporary range of Atlantic cod. Diversity, measured as average sample heterozygosity and number of polymorphic loci, declined (c. 30%) from the western (He = 0.36) to eastern (He = 0.25) Atlantic, consistent with a signal of ascertainment bias. Assignment success was examined separately for pools of loci representing differing degrees of reductions in diversity. SNPs displaying the largest declines in diversity produced the most accurate assignment in the ascertainment region (c. 83%) and the lowest levels of correct assignment outside the ascertainment region (c. 31%). Interestingly, several isolated locations showed no effect of assignment bias and consistently displayed 100% correct assignment. Contrary to expectations, estimates of accurate assignment range‐wide using all loci displayed remarkable similarity despite reductions in diversity. Our results support the use of large SNP panels in assignment studies of high geneflow marine species. However, our evidence of significant reductions in assignment success using some pools of loci suggests that ascertainment bias may influence assignment results and should be evaluated in large‐scale assignment studies.

Dispersal Patterns, Active Behaviour, and Flow Environment during Early Life History of Coastal Cold Water Fishes

During the pelagic larval phase, fish dispersal may be influenced passively by surface currents or actively determined by swimming behaviour. In situ observations of larval swimming are few given the constraints of field sampling. Active behaviour is therefore often inferred from spatial patterns in the field, laboratory studies, or hydrodynamic theory, but rarely are these approaches considered in concert. Ichthyoplankton survey data collected during 2004 and 2006 from coastal Newfoundland show that changes in spatial heterogeneity for multiple species do not conform to predictions based on passive transport. We evaluated the interaction of individual larvae with their environment by calculating Reynolds number as a function of ontogeny. Typically, larvae hatch into a viscous environment in which swimming is inefficient, and later grow into more efficient intermediate and inertial swimming environments. Swimming is therefore closely related to length, not only because of swimming capacity but also in how larvae experience viscosity. Six of eight species sampled demonstrated consistent changes in spatial patchiness and concomitant increases in spatial heterogeneity as they transitioned into more favourable hydrodynamic swimming environments, suggesting an active behavioural element to dispersal. We propose the tandem assessment of spatial heterogeneity and hydrodynamic environment as a potential approach to understand and predict the onset of ecologically significant swimming behaviour of larval fishes in the field.

Operationalizing integrated ecosystem assessments within a multidisciplinary team: lessons learned from a worked example

Operationalizing integrated ecosystem assessments within a multidisciplinary team: lessons learned from a worked example Geret S. DePiper*, Sarah K. Gaichas, Sean M. Lucey, Patricia Pinto da Silva, M. Robin Anderson, Heather Breeze, Alida Bundy, Patricia M. Clay, Gavin Fay, Robert J. Gamble, Robert S. Gregory, Paula S. Fratantoni, Catherine L. Johnson, Mariano Koen-Alonso, Kristin M. Kleisner, Julia Olson, Charles T. Perretti, Pierre Pepin, Fred Phelan, Vincent S. Saba, Laurel A. Smith, Jamie C. Tam, Nadine D. Templeman, and Robert P. Wildermuth NOAA Northeast Fisheries Science Center, 166 Water Street, Woods Hole, MA 02543, USA Fisheries and Oceans Canada, Northwest Atlantic Fisheries Centre, 80 East White Hills, St. John’s, NL A1C 5X1, Canada Fisheries and Oceans Canada, Bedford Institute of Oceanography, 1 Challenger Drive, Dartmouth, NS B2Y 4A2, Canada School for Marine Science & Technology, University of Massachusetts Dartmouth, 200 Mill Road, Suite 30, Fairhaven, MA 02719, USA Environmental Defense Fund, Floor 28, 123 Mission Street, San Francisco, CA 94105, USA National Oceanic and Atmospheric Administration, National Marine Fisheries Service, Northeast Fisheries Science Center, Geophysical Fluid Dynamics Laboratory, Princeton University Forrestal Campus, 201 Forrestal Road, Princeton, NJ 08540, USA

Factors Regulating Early Life History Dispersal of Atlantic Cod (Gadus morhua) from Coastal Newfoundland

To understand coastal dispersal dynamics of Atlantic cod (Gadus morhua), we examined spatiotemporal egg and larval abundance patterns in coastal Newfoundland. In recent decades, Smith Sound, Trinity Bay has supported the largest known overwintering spawning aggregation of Atlantic cod in the region. We estimated spawning and dispersal characteristics for the Smith Sound-Trinity Bay system by fitting ichthyoplankton abundance data to environmentally-driven, simplified box models. Results show protracted spawning, with sharply increased egg production in early July, and limited dispersal from the Sound. The model for the entire spawning season indicates egg export from Smith Sound is 13%•day−1 with a net mortality of 27%•day–1. Eggs and larvae are consistently found in western Trinity Bay with little advection from the system. These patterns mirror particle tracking models that suggest residence times of 10–20 days, and circulation models indicating local gyres in Trinity Bay that act in concert with upwelling dynamics to retain eggs and larvae. Our results are among the first quantitative dispersal estimates from Smith Sound, linking this spawning stock to the adjacent coastal waters. These results illustrate the biophysical interplay regulating dispersal and connectivity originating from inshore spawning of coastal northwest Atlantic.

Structural and functional connectivity of marine fishes within a semi-enclosed Newfoundland fjord

The interplay between structural connectivity (i.e. habitat continuity) and functional connectivity (i.e. dispersal probability) in marine fishes was examined in a coastal fjord (Holyrood Pond, Newfoundland, Canada) that is completely isolated from the North Atlantic Ocean for most of the year. Genetic differentiation was described in three species (rainbow smelt Osmerus mordax, white hake Urophycis tenuis and Atlantic cod Gadus morhua) with contrasting life histories using seven to 10 microsatellite loci and a protein-coding locus, PanI (G. morhua). Analysis of microsatellite differentiation indicated clear genetic differences between the fjord and coastal regions; however, the magnitude of difference was no more elevated than adjacent bays and was not enhanced by the fjords isolation. Osmerus mordax was characterized by the highest structure overall with moderate differentiation between the fjord and St Marys Bay (F(ST)c.0.047). In contrast, U. tenuis and G. morhua displayed weak differentiation (F(ST) < 0.01). Nonetheless, these populations did demonstrate high rates (< 75%) of Bayesian self-assignment. Furthermore, elevated differentiation was observed at the PanI locus in G. morhua between the fjord and other coastal locations. Interestingly, locus-specific genetic differentiation and expected heterozygosity were negatively associated in O. mordax, in contrast to the positive associations observed in U. tenuis and G. morhua. Gene flow in these species is apparently unencumbered by limited structural connectivity, yet the observed differentiation suggests that population structuring exists over small scales despite high dispersal potential.