Footprints of local adaptation span hundreds of linked genes in the Atlantic silverside genome.

The study of local adaptation in the presence of ongoing gene flow is the study of natural selection in action, revealing the functional genetic diversity most relevant to contemporary pressures. In addition to individual genes, genome-wide architecture can itself evolve to enable adaptation. Distributed across a steep thermal gradient along the east coast of North America, Atlantic silversides (Menidia menidia) exhibit an extraordinary degree of local adaptation in a suite of traits, and the capacity for rapid adaptation from standing genetic variation, but we know little about the patterns of genomic variation across the species range that enable this remarkable adaptability. Here, we use low-coverage, whole-transcriptome sequencing of Atlantic silversides sampled along an environmental cline to show marked signatures of divergent selection across a gradient of neutral differentiation. Atlantic silversides sampled across 1371 km of the southern section of its distribution have very low genome-wide differentiation (median F ST = 0.006 across 1.9 million variants), consistent with historical connectivity and observations of recent migrants. Yet almost 14,000 single nucleotide polymorphisms (SNPs) are nearly fixed (F ST > 0.95) for alternate alleles. Highly differentiated SNPs cluster into four tight linkage disequilibrium (LD) blocks that span hundreds of genes and several megabases. Variants in these LD blocks are disproportionately nonsynonymous and concentrated in genes enriched for multiple functions related to known adaptations in silversides, including variation in lipid storage, metabolic rate, and spawning behavior. Elevated levels of absolute divergence and demographic modeling suggest selection maintaining divergence across these blocks under gene flow. These findings represent an extreme case of heterogeneity in levels of differentiation across the genome, and highlight how gene flow shapes genomic architecture in continuous populations. Locally adapted alleles may be common features of populations distributed along environmental gradients, and will likely be key to conserving variation to enable future responses to environmental change.

Contrasting genomic shifts underlie parallel phenotypic evolution in response to fishing.

Humans cause widespread evolutionary change in nature, but we still know little about the genomic basis of rapid adaptation in the Anthropocene. We tracked genomic changes across all protein-coding genes in experimental fish populations that evolved pronounced shifts in growth rates due to size-selective harvest over only four generations. Comparisons of replicate lines show parallel allele frequency shifts that recapitulate responses to size-selection gradients in the wild across hundreds of unlinked variants concentrated in growth-related genes. However, a supercluster of genes also rose rapidly in frequency and dominated the evolutionary dynamic in one replicate line but not in others. Parallel phenotypic changes thus masked highly divergent genomic responses to selection, illustrating how contingent rapid adaptation can be in the face of strong human-induced selection.

Mercury bioaccumulation increases with latitude in a coastal marine fish (Atlantic silverside, Menidia menidia).

Human exposure to the neurotoxic methylmercury (MeHg) occurs primarily via the consumption of marine fish, but the processes underlying large-scale spatial variations in fish MeHg concentrations [MeHg], which influence human exposure, are not sufficiently understood. We used the Atlantic silverside (Menidia menidia), an extensively studied model species and important forage fish, to examine latitudinal patterns in total Hg [Hg] and [MeHg]. Both [Hg] and [MeHg] significantly increased with latitude (0.014 and 0.048 µg MeHg g-1 dw per degree of latitude in juveniles and adults, respectively). Four known latitudinal trends in silverside traits help explain these patterns: latitudinal increase in MeHg assimilation efficiency, latitudinal decrease in MeHg efflux, latitudinal increase in weight loss due to longer and more severe winters, and latitudinal increase in food consumption as an adaptation to decreasing length of the growing season. Given the absence of a latitudinal pattern in particulate MeHg, a diet proxy for zooplanktivorous fish, we conclude that large-scale spatial variation in growth is the primary control of Hg bioaccumulation in this and potentially other fish species.

Pattern and scale of geographic variation in environmental sex determination in the Atlantic silverside, Menidia menidia.

The Atlantic silverside, Menidia menidia (Pisces: Atherinidae), exhibits an exceptionally high level of clinal variation in sex determination across its geographic range. Previous work suggested linear changes in the level of temperature-dependent sex determination (TSD) with increasing latitude. Based on comparisons at 31 sites encompassing the entire species' range, we find that the change in level of TSD with latitude is instead highly nonlinear. The level of TSD is uniformly high in the south (Florida to New Jersey), then declines rapidly into the northern Gulf of Maine where genotypic sex determination (GSD) predominates and then rebounds to moderate levels of TSD in the northern-most populations of the Gulf of St. Lawrence. Major latitudinal breakpoints occur in central New Jersey (40(o)N) and the northern Gulf of Maine (44(o)N). No populations display pure TSD or GSD. Length of the growing season is the likely agent of selection driving variation in TSD with a threshold at 210 days. Because gene flow among populations is high, such distinct patterns of geographic variation in TSD/GSD are likely maintained by contemporary selection thereby demonstrating the adaptive fine tuning of sex determining mechanisms.

Circulating levels of plasma IGF-I during recovery from size-selective harvesting in Menidia menidia.

Selection for growth-related traits in domesticated fishes often results in predictable changes within the growth hormone-insulin-like growth factor (GH-IGF-1) axis. Little is known about the mechanisms controlling changes in growth capacity resulting from fishery-induced evolution. We took advantage of a long-term study where Menidia menidia were selected for size at age over multiple generations to mimic fisheries-induced selection. This selection regime produced three populations with significant differences in intrinsic growth rate. These growth differences partially rebounded, but persisted even after selection was relaxed, resulting in fast, intermediate, and slow-growing lines. Plasma IGF-1 was measured in these populations as a potential target of selection on growth. IGF-1 was significantly correlated with current length and mass, and was positively correlated with growth rate (g d(-1)) in two lines, indicating it may be an appropriate indicator of growth capacity. The slow-growing line exhibited higher overall IGF-1 levels relative to the depressed IGF-1 seen in the fast-growing line, contrary to our prediction. We offer possible explanations for this unusual pattern and argue that somatic growth is likely to be under control of mechanism(s) downstream to IGF-1. IGF-1 provides an interesting basis for understanding endocrine control of growth in response to artificial selection and recovery.

The response of correlated traits following cessation of fishery-induced selection.

The application of evolutionary principles to the management of fisheries has gained considerable attention recently. Harvesting of fish may apply directional or disruptive selection to key life-history traits, and evidence for fishery-induced evolution is growing. The traits that are directly selected upon are often correlated (genetically or phenotypically) with a suite of interrelated physiological, behavioral, and morphological characters. A question that has received comparatively little attention is whether or not, after cessation of fishery-induced selection, these correlated traits revert back to previous states. Here, we empirically examine this question. In experiments with the Atlantic silverside, Menidia menidia, we applied size-selective culling for five generations and then maintained the lines a further five generations under random harvesting. We found that some traits do return to preharvesting levels (e.g., larval viability), some partially recover (e.g., egg volume, size-at-hatch), and others show no sign of change (e.g., food consumption rate, vertebral number). Such correlations among characters could, in theory, greatly accelerate or decelerate the recovery of fish populations. These results may explain why some fish stocks fail to recover after fishing pressure is relaxed.

Spatial scale and divergent patterns of variation in adapted traits in the ocean.

The geography of adaptive genetic variation is crucial to species conservation yet poorly understood in marine systems. We analyse the spatial scale of genetic variation in traits that broadly display adaptation throughout the range of a highly dispersive marine species. We conducted common garden experiments on the Atlantic silverside, Menidia menidia, from 39 locations along its 3000 km range thereby mapping genetic variation for growth rate, vertebral number and sex determination. Each trait displayed unique clinal patterns, with significant differences (adaptive or not) occurring over very small distances. Breakpoints in the cline differed among traits, corresponding only partially with presumed eco-geographical boundaries. Because clinal patterns are unique to each selected character, neutral genes or those coding for a single character cannot serve as proxies for the genetic structure as a whole. Conservation plans designed to protect essential genetic subunits of a species will need to account for such complex spatial structures.

Adaptation to climate change: contrasting patterns of thermal-reaction-norm evolution in Pacific versus Atlantic silversides.

How organisms may adapt to rising global temperatures is uncertain, but concepts can emerge from studying adaptive physiological trait variations across existing spatial climate gradients. Many ectotherms, particularly fish, have evolved increasing genetic growth capacities with latitude (i.e. countergradient variation (CnGV) in growth), which are thought to be an adaptation primarily to strong gradients in seasonality. In contrast, evolutionary responses to gradients in mean temperature are often assumed to involve an alternative mode, 'thermal adaptation'. We measured thermal growth reaction norms in Pacific silverside populations (Atherinops affinis) occurring across a weak latitudinal temperature gradient with invariant seasonality along the North American Pacific coast. Instead of thermal adaptation, we found novel evidence for CnGV in growth, suggesting that CnGV is a ubiquitous mode of reaction-norm evolution in ectotherms even in response to weak spatial and, by inference, temporal climate gradients. A novel, large-scale comparison between ecologically equivalent Pacific versus Atlantic silversides (Menidia menidia) revealed how closely growth CnGV patterns reflect their respective climate gradients. While steep growth reaction norms and increasing growth plasticity with latitude in M. menidia mimicked the strong, highly seasonal Atlantic coastal gradient, shallow reaction norms and much smaller, latitude-independent growth plasticity in A. affinis resembled the weak Pacific latitudinal temperature gradient.

Regional differentiation and post-glacial expansion of the Atlantic silverside, Menidia menidia, an annual fish with high dispersal potential.

The coastal marine environment of the Northwest Atlantic contains strong environmental gradients that create distinct marine biogeographic provinces by limiting dispersal, recruitment, and survival. This region has also been subjected to numerous Pleistocene glacial cycles, resulting in repeated extirpations and recolonizations in northern populations of marine organisms. In this study, we examined patterns of genetic structure and historical demography in the Atlantic silverside, Menidia menidia, an annual marine fish with high dispersal potential but with well-documented patterns of clinal phenotypic adaptation along the environmental gradients of the Northwest Atlantic. Contrary to previous studies indicating genetic homogeneity that should preclude regional adaptation, results demonstrate subtle but significant (FST = 0.07; P < 0.0001) genetic structure among three phylogeographic regions that partially correspond with biogeographic provinces, suggesting regional limits to gene flow. Tests for non-equilibrium population dynamics and latitudinal patterns in genetic diversity indicate northward population expansion from a single southern refugium following the last glacial maximum, suggesting that phylogeographic and phenotypic patterns have relatively recent origins. The recovery of phylogeographic structure and the partial correspondence of these regions to recognized biogeographic provinces suggest that the environmental gradients that shape biogeographic patterns in the Northwest Atlantic may also limit gene flow in M. menidia, creating phylogeographic structure and contributing to the creation of latitudinal phenotypic clines in this species.

Ontogenesis of gonadal aromatase gene expression in atlantic silverside (Menidia menidia) populations with genetic and temperature-dependent sex determination.

Cytochrome P450 aromatase (P450arom), an enzyme that converts testosterone to 17beta-estradiol, is an important mediator of sex determination in teleosts with genetic sex determination (GSD) and temperature-dependent sex determination (TSD). We compared the ontogenetic expression of P450arom in two populations of Atlantic silversides, Menidia menidia, which exhibit TSD (South Carolina) or GSD (Nova Scotia, Canada) using quantitative, real-time polymerase chain reaction (qRT-PCR). Embryos and newly hatched larvae were reared at an intermediate sex ratio-producing temperature (21 degrees C), and older larvae and juveniles were reared at temperatures that feminize (15 degrees C) and masculinize (28 degrees C) to assess the temperature response of P450arom during development. Before sex determination, embryos and newly-hatched larvae displayed negligible P450arom expression, indicating minimal upregulation of this gene before sex determination. Gene expression increased in both populations during sex differentiation. Nova Scotia fish with GSD exhibited presumptive male- and female-like expression levels during early sex differentiation that were not influenced by temperature. South Carolina fish displayed low levels of expression at 28 degrees C with significantly heightened expression in some individuals at 15 degrees C, indicating that P450arom is temperature sensitive in the population with TSD. Populations also differed in the timing and maximal levels of P450arom expression, with fish from Nova Scotia exhibiting both the highest and earliest increase in expression in presumptive females. Our results support the hypothesis that P450arom is involved in female sex differentiation in this species, but is only responsive to temperature in M. menidia populations that exhibit TSD.

High connectivity among locally adapted populations of a marine fish (Menidia menidia).

Patterns of connectivity are important in understanding the geographic scale of local adaptation in marine populations. While natural selection can lead to local adaptation, high connectivity can diminish the potential for such adaptation to occur. Connectivity, defined as the exchange of individuals among subpopulations, is presumed to be significant in most marine species due to life histories that include widely dispersive stages. However, evidence of local adaptation in marine species, such the Atlantic silverside, Menidia menidia, raises questions concerning the degree of connectivity. We examined geochemical signatures in the otoliths, or ear bones, of adult Atlantic silversides collected in 11 locations along the northeastern coast of the United States from New Jersey to Maine in 2004 and eight locations in 2005 using laser ablation inductively coupled plasma mass spectrometry (ICP-MS) and isotope ratio monitoring mass spectrometry (irm-MS). These signatures were then compared to baseline signatures of juvenile fish of known origin to determine natal origin of these adult fish. We then estimated migration distances and the degree of mixing from these data. In both years, fish generally had the highest probability of originating from the same location in which they were captured (0.01-0.80), but evidence of mixing throughout the sample area was present. Furthermore, adult M. menidia exhibit highly dispersive behavior with some fish migrating over 700 km. The probability of adult fish returning to natal areas differed between years, with the probability being, on average, 0.2 higher in the second year. These findings demonstrate that marine species with largely open populations are capable of local adaptation despite apparently high gene flow.

The covariance between genetic and environmental influences across ecological gradients: reassessing the evolutionary significance of countergradient and cogradient variation.

Patterns of phenotypic change across environmental gradients (e.g., latitude, altitude) have long captivated the interest of evolutionary ecologists. The pattern and magnitude of phenotypic change is determined by the covariance between genetic and environmental influences across a gradient. Cogradient variation (CoGV) occurs when covariance is positive: that is, genetic and environmental influences on phenotypic expression are aligned and their joint influence accentuates the change in mean trait value across the gradient. Conversely, countergradient variation (CnGV) occurs when covariance is negative: that is, genetic and environmental influences on phenotypes oppose one another, thereby diminishing the change in mean trait expression across the gradient. CnGV has so far been found in at least 60 species, with most examples coming from fishes, amphibians, and insects across latitudinal or altitudinal gradients. Traits that display CnGV most often involve metabolic compensation, that is, the elevation of various physiological rates processes (development, growth, feeding, metabolism, activity) to counteract the dampening effect of reduced temperature, growing season length, or food supply. Far fewer examples of CoGV have been identified (11 species), and these most often involve morphological characters. Increased knowledge of spatial covariance patterns has furthered our understanding of Bergmann size clines, phenotypic plasticity, species range limits, tradeoffs in juvenile growth rate, and the design of conservation strategies for wild species. Moreover, temporal CnGV explains some cases of an apparent lack of phenotypic response to directional selection and provides a framework for predicting evolutionary responses to climate change.

Reversal of evolutionary downsizing caused by selective harvest of large fish.

Evolutionary responses to the long-term exploitation of individuals from a population may include reduced growth rate, age at maturation, body size and productivity. Theoretical models suggest that these genetic changes may be slow or impossible to reverse but rigorous empirical evidence is lacking. Here, we provide the first empirical demonstration of a genetically based reversal of fishing-induced evolution. We subjected six populations of silverside fish (Menidia menidia) to three forms of size-selective fishing for five generations, thereby generating twofold differences among populations in mean weight and yield (biomass) at harvest. This was followed by an additional five generations during which size-selective harvest was halted. We found that evolutionary changes were reversible. Populations evolving smaller body size when subjected to size-selective fishing displayed a slow but significant increase in size when fishing ceased. Neither phenotypic variance in size nor juvenile survival was reduced by the initial period of selective fishing, suggesting that sufficient genetic variation remained to allow recovery. By linear extrapolation, we predict full recovery in about 12 generations, although the rate of recovery may taper off near convergence. The recovery rate in any given wild population will also depend on other agents of selection determined by the specifics of life history and environment. By contrast, populations that in the first five generations evolved larger size and yield showed little evidence of reversal. These results show that populations have an intrinsic capacity to recover genetically from harmful evolutionary changes caused by fishing, even without extrinsic factors that reverse the selection gradient. However, harvested species typically have generation times of 3-7 years, so recovery may take decades. Hence, the need to account for evolution in managing fisheries remains.

The role of experiments in understanding fishery-induced evolution.

Evidence of fishery-induced evolution has been accumulating rapidly from various avenues of investigation. Here we review the knowledge gained from experimental approaches. The strength of experiments is in their ability to disentangle genetic from environmental differences. Common garden experiments have provided direct evidence of adaptive divergence in the wild and therefore the evolvability of various traits that influence production in numerous species. Most of these cases involve countergradient variation in physiological, life history, and behavioral traits. Selection experiments have provided examples of rapid life history evolution and, more importantly, that fishery-induced selection pressures cause simultaneous divergence of not one but a cluster of genetically and phenotypically correlated traits that include physiology, behavior, reproduction, and other life history characters. The drawbacks of experiments are uncertainties in the scale-up from small, simple environments to larger and more complex systems; the concern that taxons with short life cycles used for experimental research are atypical of those of harvested species; and the difficulty of adequately simulating selection due to fishing. Despite these limitations, experiments have contributed greatly to our understanding of fishery-induced evolution on both empirical and theoretical levels. Future advances will depend on integrating knowledge from experiments with those from modeling, field studies, and molecular genetic approaches.

Coevolution of foraging behavior with intrinsic growth rate: risk-taking in naturally and artificially selected growth genotypes of Menidia menidia.

Although there is accumulating evidence of growth-rate optimization by natural selection, the coevolution of growth rate and risk-taking behavior has not been sufficiently documented. The Atlantic silverside fish, Menidia menidia, displays countergradient variation in growth across a latitudinal gradient: genotypes from Nova Scotia (NS), for example, grow in length twofold faster than those from South Carolina (SC). Past work has established that fast growth is adaptive in northern climates, but the trade-off is poorer swimming performance and higher susceptibility to predators. We compared escape behavior and willingness to forage under threat of predation among growth genotypes reared and tested under common-garden conditions. When chased with a predator model, NS fish occupied shelter more quickly than SC fish. When food was supplied after a chase, NS fish reemerged from the shelter much more quickly than SC fish and immediately commenced feeding, whereas many SC fish displayed timid behavior and did not feed. When food was absent following a chase, however, NS fish remained in the shelter longer than did SC fish and both displayed timid behavior. Hence, the fast-growing NS genotype was bolder than SC fish in the presence of food, but shyer in the absence of food. These behaviors are adaptive given the physiological constraints intrinsic to each genotype. Experiments on captive populations of silversides that had been artificially selected for fast or slow growth confirmed that foraging behavior is genetically correlated with intrinsic growth rate, although in these trials the fast-growth genotype was always more bold, regardless of food availability, as would be expected in the absence of predators. We conclude that risk-taking foraging behavior coevolves adaptively with intrinsic growth rate in M. menidia.

Maladaptive changes in multiple traits caused by fishing: impediments to population recovery.

Some overharvested fish populations fail to recover even after considerable reductions in fishing pressure. The reasons are unclear but may involve genetic changes in life history traits that are detrimental to population growth when natural environmental factors prevail. We empirically modelled this process by subjecting populations of a harvested marine fish, the Atlantic silverside, to experimental size-biased fishing regimes over five generations and then measured correlated responses across multiple traits. Populations where large fish were selectively harvested (as in most fisheries) displayed substantial declines in fecundity, egg volume, larval size at hatch, larval viability, larval growth rates, food consumption rate and conversion efficiency, vertebral number, and willingness to forage. These genetically based changes in numerous traits generally reduce the capacity for population recovery.

Tests of bioaccumulation models for polychlorinated biphenyl compounds: a study of young-of-the-year bluefish in the Hudson River estuary, USA.

A field-based study regarding uptake of polychlorinated biphenyl compounds (PCBs) by young-of-the-year (YOY) bluefish (Pomatomus saltatrix) was initiated to test a steady-state model of bioaccumulation and trophic transfer in a rapidly growing fish. Determination of prey composition as well as size-dependent growth and specific consumption rates for YOY bluefish from separate field and laboratory studies enabled the input of these species-specific parameters into the model. Furthermore, the time and duration of the exposure of YOY bluefish to dissolved PCBs from a well-characterized system (Hudson River, USA) was well known. Patterns of accumulation of individual PCB congeners differed relative to the accumulation of total PCBs, with the greatest net accumulation occurring for the higher-molecular-weight congeners. Comparison of lipid-normalized bioaccumulation factors (BAFs) with the octanol-water partition coefficients of individual PCB congeners revealed bluefish to be above the BAFs predicted by lipid-based equilibrium partitioning, suggesting that uptake from food is an important source of PCBs in YOY bluefish. Comparison of measured BAFs with values predicted by a steady-state, food-chain model showed good first-order agreement.

Evolution of intrinsic growth rate: metabolic costs drive trade-offs between growth and swimming performance in Menidia menidia.

There is strong evidence that genetic capacity for growth evolves toward an optimum rather than an absolute maximum. This implies that fast growth has a cost and that trade-offs occur between growth and other life-history traits, but the fundamental mechanisms are poorly understood. Previous work on the Atlantic silverside fish Menidia menidia has demonstrated a trade-off between growth and swimming performance. We hypothesize that the trade-off derives from the competing metabolic demands associated with growth and swimming activity. We tested this by measuring standard metabolic rate (M(STD)), maximum sustainable metabolic rate (M(ACT)) and metabolic scope of laboratory-reared silversides originating from two geographically distinct populations with well-documented differences in genetic capacity for growth. The fast-growth genotype had a significantly greater M(STD) than the slow-growth genotype, but a similar MACT when swum to near exhaustion. The scope for activity of the fast-growth genotype was lower than that of the slow-growth genotype. Furthermore, the fast-growth genotype eats larger meals, thereby incurring a greater postprandial oxygen demand. We conclude that a metabolic trade-off occurs between growth and other metabolic demands and that this trade-off provides a general mechanism underlying the evolution of growth rate.