Intense Habitat-Specific Fisheries-Induced Selection at the Molecular Pan I Locus Predicts Imminent Collapse of a Major Cod Fishery

Predation is a powerful agent in the ecology and evolution of predator and prey. Prey may select multiple habitats whereby different genotypes prefer different habitats. If the predator is also habitat-specific the prey may evolve different habitat occupancy. Drastic changes can occur in the relation of the predator to the evolved prey. Fisheries exert powerful predation and can be a potent evolutionary force. Fisheries-induced selection can lead to phenotypic changes that influence the collapse and recovery of the fishery. However, heritability of the phenotypic traits involved and selection intensities are low suggesting that fisheries-induced evolution occurs at moderate rates at decadal time scales. The Pantophysin I (Pan I) locus in Atlantic cod (Gadus morhua), representing an ancient balanced polymorphism predating the split of cod and its sister species, is under an unusual mix of balancing and directional selection including current selective sweeps. Here we show that Pan I alleles are highly correlated with depth with a gradient of 0.44% allele frequency change per meter. AA fish are shallow-water and BB deep-water adapted in accordance with behavioral studies using data storage tags showing habitat selection by Pan I genotype. AB fish are somewhat intermediate although closer to AA. Furthermore, using a sampling design covering space and time we detect intense habitat-specific fisheries-induced selection against the shallow-water adapted fish with an average 8% allele frequency change per year within year class. Genotypic fitness estimates (0.08, 0.27, 1.00 of AA, AB, and BB respectively) predict rapid disappearance of shallow-water adapted fish. Ecological and evolutionary time scales, therefore, are congruent. We hypothesize a potential collapse of the fishery. We find that probabilistic maturation reaction norms for Atlantic cod at Iceland show declining length and age at maturing comparable to changes that preceded the collapse of northern cod at Newfoundland, further supporting the hypothesis. We speculate that immediate establishment of large no-take reserves may help avert collapse.


Introduction
Predation is a powerful agent in the ecology and evolution of predator and prey. Prey may select multiple habitats whereby different genotypes prefer different habitats. If the predator is also habitat-specific the prey may evolve different habitat occupancy. Drastic changes can occur in the relation of the predator to the evolved prey. For example the predator may exterminate the prey in the habitat available to the predator and thus lose its prey. The prey selecting the alternative habitat may be released from predation and might even evolve to become a new species free of predation. The stronger the predation and more efficient the predator, the more drastic would be the evolutionary changes. It is, therefore, important to study the way predation acts on genes in species that select multiple habitats by genotype.
Man the hunter has become a mechanized techno-beast, a highly efficient predator. In particular, commercial fisheries searching for fish with computerized fish-finders and airplanes and scooping up fish with several thousand-ton capacities with ships powered by several thousand horsepower engines are a case in point [1]. Modern fisheries are uncontrolled experiments in evolution [2,3]. Fisheries target certain phenotypes and, therefore, can be a powerful agent of natural selection. They also frequently are habitat-specific by concentrating on the most accessible locations.
Fisheries-induced selection [4][5][6][7][8][9][10] can lead to phenotypic changes that influence the collapse and recovery of fisheries [4,5,11]. Fisheries-induced selection is primarily discussed as selective mortality directly targeting certain phenotypes such as length-or weight-at-age [4,6,7] by size-selective fishing or indirectly targeting phenotypes such as age-at-maturation due to age-at-entry into the fisheries or by location of fishery [4]. However, the genetic determination and heritability of the quantitative traits involved is largely unknown [4,8,12]. General estimates of the heritability of life-history traits are low as are estimated selection intensities although fishing mortality can be very high [13,14]. From these observations emerges the view that fisheries-induced evolution occurs at moderate rates with significant changes only observable at decadal time scales [13]. Nevertheless, the recommendation is that the management adopt short-term conservation/management measures that also meet concerns about long-term evolutionary impact. Thus conservation of old, big fish is promoted as a combined short-and long-term conservation strategy [13]. Evolutionary changes were implicated in the collapse and non-recovery of the northern cod of Newfoundland and in the near collapse of North Sea cod [11,13,14]. Following the collapse of the northern cod fishery in Newfoundland [11], the cod fishery at Iceland has remained as one of the worlds major and ostensibly sustainable cod fishery.
The single molecular Pantophysin I (Pan I) locus represents an ancient balanced polymorphism predating the split of Atlantic cod (Gadus morhua) and its sister species [15]. The A and B alleles at the locus differ by multiple nucleotide and six amino-acid substitutions. The locus exhibits local adaptation and is under an unusual mix of balancing and directional selection including current selective sweeps [16,17]. Pan I variation also has been taken as evidence for population sub-structuring caused by breeding structure [18,19]. However, the precise roles that fisheries-induced selection and natural selection play in generating that structure are unknown.
Here we test the hypothesis that fisheries can exert powerful habitat-specific selection on a fish species with consequent evolutionary changes. We report exceptionally strong selective changes at the Pan I locus which is involved in habitat selection of depth by genotype in populations of Atlantic cod at Iceland. With a heritability of 100% and indirect but intense selection due to fishing in preferred habitat we observe significant changes on a yearly basis. Thus, fisheries-induced evolution is short term [20]. We identify a threat of collapse due to the selection imposed by the fishery. We consider that if management acted immediately it may be possible to avert collapse. Our study demonstrates the importance of molecular population genetics of targeted loci for studies of fisheries-induced selection and highlights the importance of evolutionary thinking for both short-and long-term management of exploited fish populations.

Allele frequency and depth of sampling
The relationship of allele frequency and depth was highly regular and significant ( Figure 1, Table 1 and Figure S1; and see Table S1 for an overview of supplementary data) with frequency of the A allele decreasing rapidly to a depth of less than 200 m but staying relatively level in deeper waters. A linear regression equation of allele frequency of A on depth for depth less than 200 m was P A = 0.80620.0044D (t = 214.6, P%0.001) or a 0.44% drop in allele frequency per meter. A generalized linear model (glm) fit gave almost identical results ( Figure 1).
In spring, significant heterozygote deficiency characterized the three top 25 m depth layers, an apparent Wahlund effect [21] due to the convergence in shallow waters of groups of fish that differed in allele frequencies. AA fish were rare in deep water in spring but were found at all depths in fall although they preferred shallow water (Table 1). In fall significant deviations from Hardy-Weinberg were found at various depths representing heterozygote excess in all instances. For fall and spring combined the sign of F IS was most often negative (sign test, P,0.01) indicating a general tendency for heterozygote excess.
Population differentiation in a two and a three dimensional habitat Spatial population differentiation exhibited by Pan I [19] is more fully explained by differences in depth among localities although depth and locality are confounded. We found large differences in allele frequencies among localities defined by pooling sampling stations in squares within areas ( Figure 2; we pooled to increase sample size). Close inspection of the figure, however, shows that there can be large allele frequency differences among neighboring localities within divisions. These may be described as an inshore/offshore difference, however, depth is a more important explanatory variable. We found apparent spatial differentiation significant in all instances (P = 0.001 based on 1000 permutations in all instances. Table 2). Interestingly, the overall F ST = 0.074 was considerably lower than found for Northeast-Southwest comparison [19]. Also the differentiation among sampling stations within divisions was higher than among divisions ( Table 2). Considering depth of sampling there is clearly a confounding of depth and geographic location. Shallow stations are located in the Northeast and North whereas very deep-water stations are exclusively in the South and Southwest ( Figure S2 and Figure  S3). Thus, the greatest contrast in depth was between the Northeast and South/Southwest. Given the gradient of 0.44% change in allele frequency with 1 m change in depth (Figure 1 and Figure S1), depth differences between the Northeast and South/ Southwest very likely contributed to the apparent spatial differentiation of Icelandic cod [19]. There were in some instances large differences in allele frequencies among neighboring stations within divisions ( Figure 2) particularly in divisions that included sampling stations of different depths ( Figure S3). Because depth and location are essentially a factorial or crossed design [22] we cannot use hierarchical or nested F-statistics to test their effects. However, a crossed factor can be tested independently within a level of the other factor [23] and probabilities combined in meta analysis [24]. Differentiation by depth within divisions (Table 3) was relatively high and significant in all instances except division 3. The test combining probabilities was also significant (X 2~{ 2 X log P i ð Þ~97:37, df = 16, P%0.001). In contrast, considering divisions within depth classes differentiation was considerably lower (Table 4) but nevertheless significant in several instances and overall (X 2~{ 2 X log P i ð Þ~103:20, df = 22, P%0.001).
Overall, therefore, there was greater differentiation by depth than by geographic locality but the two factors remain statistically confounded. Depth of course is a proxy for some biological and environmental factors [25]. For the purposes of this paper the relationship with depth is of great importance. Even if significant spatial variation existed that was not confounded by depth it would not alter implications of the selective effects discussed here.

Fishing pressure
Data on Atlantic cod catches show that, in general, the brunt of the fishing is carried by 5-7 year old fish. By eight years of age a year class is severely reduced and by nine years it is all but fished out ( Figure S4; and see Table S1 for key points of supplementary data). Furthermore, analysis of catch by gear shows that the heavy fishing occurs in shallow water ( Figure S5). Catch by long line is to some extent conducted in deeper water (200-300 m) with a shift towards shallower waters in recent years as shown by the higher density of catch at depths less than 100 m and lower density at 200-300 m depths ( Figure S5). Although also taking considerable catch in shallow waters, fishing by bottom trawl is distributed over the greatest range of depth ( Figure S5). Bottom trawl which targets greatest depth range brings in more than 40% of the catch ( Figure S6). From 1997 to 2007 the catch taken by gear targeting shallow water has increased. Catch with hand line, Danish seine, and gill nets increased up until 2003. However, their catch have decreased since then. Catch with bottom trawl has decreased but catch with long line has increased substantially and after 2005 nearly matches the catch taken by bottom trawl. Thus, fishing mortality is heaviest in shallow water [26] and at least after 2000 fishing targeting shallow water has increased.
Total catch is regulated by government issued quotas of total allowable catch (TAC) which are based on advice from the Marine Research Institute (MRI) of recommended total allowable catch (recommended TAC). Current management strategies of Icelandic cod are based on a catch limitation system where each vessel is allowed a certain share of the TAC. Annual recommended TACs  are based on scientific assessment of state of fish stocks and ecosystem condition, but have been reduced lately. For example, the TAC for the 2007/2008 fishing year was 130,000 tons, a 63,000 ton reduction from the previous year. The recommended TAC for the following fishing year was a further 7000 ton reduction but 160,000 was issued. There is a tight correlation of total catch, issued TAC, and recommended TAC. However, total catch almost invariably exceeds issued TACs which in turn are always greater than recommended TAC except for 1996 to 2007 when a catch rule has been in effect ( Figure S7). Annual catch, effort, and catch per unit effort for different gear show a complex interaction with each other and with quotas and TACs ( Figure 3 and Figure S7). The stock reached an all time low population numbers in 1993-1995 [27] ( Figure S7). In the following years the MRI estimated an increase in numbers and increased quotas were issued. Catches of most gear increased and increased effort followed immediately or with a lag (of up to one year) for some gear. Both catch and effort peaked around 2000 and 2001. Catch per unit effort peaked earlier. Catch and effort of both long line and hand line started to increase later than other gear and peaked later or in 2004 and decreased somewhat after that ( Figure 3). This, at least partly, accounts for the relative increase of catch by long line in the total catch ( Figure S6). As catch diminished after reaching a peak, effort was reduced more rapidly and thus the relative measure of catch per unit effort increased. Using the relative measure of catch per unit effort as an indicator of stock abundance must take total catch into account and whether it is increasing or decreasing [28].
Following the 1995 crash government issued a new catch rule limiting annual quotas to 25% of average fishable biomass. In 2000 and again 2001, with the benefit of hindsight, the MRI re-estimated the population numbers for the previous years. A considerable stock size overestimation and underestimation of fishing mortality was apparent [29,30], amounting to 25-50%. Issued TACs were based on the overestimated stocks and catches have been 27-40% of the fishable stock, far exceeding the target of 25%. Thus, for example, for 2000 fishing mortality was estimated at 0.86 compared to approximately 0.4 if the catch had been at the 25% target [30]. Consequently the population has experienced increased predation pressure through increased fishing mortality [31].
Part of the overestimation is explained by a lower than predicted mean weight at age [30]. Changes in catchability  Table 4. F statistics among divisions within depth classes.  [28,31] are, however, the main explanation [30]. There were also changes towards fishing older and larger fish. As an example the gill net fleet changed most of its nets from 7 inch to 9 inch mesh size from 1994-1998 [30].

Probabilistic maturation reaction norms
We applied the principles of estimating probabilistic maturation reaction norms [32][33][34] to evaluate the potential effects of fisheries on changes in life-history traits. We used data on mean maturity, length and age [35] and estimated maturity ogives. From these we estimated the probability of maturing m at length or age ( Figure 4). Quantiles of length (and age) at both 50% and 95% probability of maturing show a significant downward trend (slope of linear regression: 20.9 and 22.2 cm per year respectively, P%0.001). Length (and age) at 5% probability of maturing increased (slope of linear regression: 0.4 cm per year, P%0.001).

Genic and genotypic frequencies at age
There were large and highly regular allele and genotypic frequency changes on age within year classes ( Figure 5 and Table 5). Year classes are independent realizations of birth and death processes yet changes were in the same direction. Allele and genotypic frequencies in older year classes appeared as continuation of frequencies among more recent year classes ( Figure 5). Comparing the 2002 to 1996 year classes, frequencies of A allele on age decreased both within and among year classes (Dp A~{ 8:4% + 4:8% per year). Similarly, frequencies of AA decreased (Df AA~{ 7:5% + 5:0%), AB increased at young ages but decreased slightly on average (Df AB~{ 1:7% + 3:0%), and BB increased (Df BB~9 :3% + 5:0%). There was a reversal with BB becoming the most common genotype by about eight years of age. The 1995 and 1994 year classes had higher A allele frequencies compared to the more recent year classes. However, comparing these two year classes the same pattern of decrease of AA and increase of BB with age held. Evolutionary changes are observable on a yearly basis [13] and thus ecological and evolutionary time scales are congruent.
We sampled three consecutive years and thus, except for the most recent and the oldest year classes, a year class entered our samples at three ages. Adjacent year classes have two overlapping age classes. The large observed evolutionary changes predict that the A allele would be decreasing in frequency. This is exactly what we observed as can be seen by comparing allele frequencies at age ( Figure 5 and Figure S8). In general allele frequencies at age are lower among the more recent year classes with a negative slope of a regression of allele frequency at age on year class ( Figure S8). This shows that A allele frequencies have decreased with time.

Fitness estimation and prediction of changes
Genotypic frequencies changed significantly between years within year class (Table 5). Frequency changes can be used to estimate relative fitness of Pan I genotypes. Overall, the genotypic frequencies changed rapidly to about eight years of age ( Figure 6) and stayed relatively level after that, an age at which the brunt of fishing of a year class is over. Catch-at-age data show that by eight years of age, a year class is severely reduced and by nine years is almost fished out ( Figure S4). Therefore, the selective pressure of the fishery is mostly over by these ages. Taking notice of this fact we took the ratio of the gam predicted frequencies ( Figure 6) among 8 year old (post-selection) to 4 year old (entering the fishery) as weights to estimate fitness (Table 6). Relative AA fitness is only 8% and AB 27% showing partial dominance. We used the upper and lower confidence limits to predict best-case and worstcase scenarios. Similar low fitness was obtained using slightly different methods for estimation (Table S2 and Table S3).
Plugging the fitness estimates (Table 6) into an equation for allele frequency change under a constant-viability selection model showed that the A allele would be eliminated in 4-5 generations, assuming continued selection of this magnitude. Given a generation time in Atlantic cod of 4.8 years [36], this gives two generations or about 10 years until predicted near disappearance of shallow-water adapted AA fish and four generations or 20 years until disappearance of heterozygous AB fish (Figure 7 panel 1). Fitness estimates from Table S2 gave similar results (Figure 7 panel 2). Using the most optimistic fitness values (Table 6) doubled the time in a best-case scenario (Figure 7 panel 3) with AB heterozygotes still making up some 15% of the population after seven generations. Under a worst-case scenario the most pessimistic values ( Table 6) predicted disappearance of AA in one generation and of AB in two generations or 10 years (Figure 7 panel 4). Similar results (not shown) were obtained with the fitness estimates in Table S2 and Table S3 obtained in slightly different ways.

Depth and population differentiation
There is an apparent spatial differentiation of Atlantic cod in Icelandic waters with differentiation between the Northeast and Southwest at F ST = 0.261 and F ST = 0.003 for the Pan I and nine microsatellite loci respectively [19]. It is important to examine whether spatial differentiation contributes to the observed selective patterns. The differentiation is clearly driven by the Pan I locus. The very low (but significant) differentiation observed at the microsatellite loci is of questionable biological significance [37] especially considering the fact that two of the nine microsatellite loci used (Gmo8 and Gmo34) are outliers in genetic differentiation [38,39] being influenced either by direct or hitch-hiking selection. In fact, in a study in northern Norway Gmo34 shows strong linkage disequilibrium with Pan I [39]. The Pan I locus is generally acknowledged to be under selection [15][16][17][18]40,41]. A differentiation of Pan I similar to the Icelandic case [19] also is observed for Arctic and Coastal cod in Norway [18]. Both studies [18,19] concluded that breeding structure with reduced gene flow was the most likely explanation for the observed spatial differentiation, yet both noticed a relationship of Pan I allele frequency and depth of sampling similar to the above (Figure 1 and Figure S1).
The apparent Wahlund effects (deficiency of heterozygotes) in shallow water in spring could be interpreted as signs of population structure by depth similar to that observed between coastal and Arctic cod in Norway [18]. To the extent that such segregation by depth actually occurs, the fisheries-induced habitat-specific selection discussed further below would be more efficient in removing shallow-water adapted fish and the effect would be harder to reverse. Conversely, if, as we suggest, there is near Figure 4. Probabilistic maturation reaction norms: length and age at 50% probability of maturing on cohort. Solid dots and solid lines represent length or age at 50% probability of maturing. Upper and lower open dots and dashed lines represent length or age at 95 and 5% probability of maturing respectively. Lines are linear regression of length or age on cohort. Based on data on mean length, age and maturity ratio from table 3.1.4 in Anonymous [27] (and see [74]). doi:10.1371/journal.pone.0005529.g004 Fisheries and Selection in Cod panmixia on breeding grounds, selection would be operating in the face of free gene flow between niches. This would reduce the efficiency of selection, it would take longer to lose the shallowwater adapted fish, and the effect would not be as difficult to reverse. A counter argument is based on the observed heterozygote excess at some depths. Overdominance in fitness is one possible explanation for heterozygote excess. However, it is very difficult to detect selection as a deviation from Hardy-Weinberg equilibrium [42]. If fitnesses form a geometric series genotypic frequencies will be in Hardy-Weinberg before and after selection [43] and directional selection may give rise to spurious overdominance if sampling is done before selection is complete [44]. If overdominance was strong enough to produce a significant heterozygote excess we would expect to see it in our fitness estimation. Such strong overdominance would also be expected to counter the heterozygote deficiency observed in spring. Unequal allele frequencies among male and female parents is another textbook explanation for heterozygote excess [45]. But neither overdominance nor male/female allele frequency differences are satisfactory explanations for the observed patterns for it is hard to see why significant excess would be restricted to certain depths at certain times of year. Negative assortative mating is another textbook explanation for heterozygote excess. The intense fisheries-induced selection observed might lead to secondary selection for apostatic mating that would reduce the effects of the selection among their offspring. Thus there could be an advantage for shallow-water adapted phenotypes to be attracted to phenotypes from the opposite deep-water habitat. Instead we consider the generation of heterozygote deficiency in spring (apparent Wahlund effect) and its disappearance in fall as well as the generation of heterozygote excess in fall and its disappearance in spring as a sign of a dynamical system. We consider the patterns most likely to be due to behavioral responses [46]. Fish fitted with data-storage tags can be classified into deep-water and shallowwater behavioral types that either forage and stay mainly in shallow water or migrate to deeper and colder waters foraging at thermal fronts and going to shallow water only for breeding [47].  These types are correlated with Pan I: AA are shallow-water type, BB are deep-water type that move to shallow waters during breeding and AB are somewhat intermediate although closer to AA than BB [46]. Pan I genotypes thus are functionally related to behavioral types that select their habitat by depth. Under natural conditions the polymorphism may be balanced by shallow-water vs. deep-water niche-variation specialization. If heterozygous AB fish move between the shallow-water and deep-water habitats to a greater extent than homozygous AA and BB fish there could be at any time an excess of heterozygotes at any depth. However, the role of various mating and feeding behaviors in generating heterozygote excess remains to be studied.

Selection at Pan I
The mechanism of balancing selection at Pan I under natural conditions is not known [16] but we suggest specialization to shallow and deep-water niches with free gene flow between the niches as the fish converge on the shallow water breeding grounds. The fitness estimates and relationships observed here, however, very likely tip the balance [16,17] leading to directional selection in favor of BB deep-water adapted fish. It is possible that some unknown fitness component would overcome the mortality leading to overdominance and balanced polymorphism but there is no evidence for that in our results. The selective pressure is most certainly due to habitat-specific fishing mortality that is heavily directed against shallow water  fish [26] ( Figure S5). If it continues unabated shallow-water fish will disappear rapidly with consequent collapse of the population and the fishery. We have not considered here the effects of selection and agestructure [48,49]. In the absence of selection a population attains a stable age-distribution. Selection will alter l x and m x schedules and will rip a population out of a stable age-distribution. Demographic processes will tend to restore a stable age-distribution and with continued selection a tug of forces will ensue. Models of intense selection in age-structured populations [49] show that agestructure and age-dependency of selection can either increase or decrease the intensity of selection. In particular, if selective mortality hits reproductive ages (as is true here) age-structure may intensify selection. Therefore, there is hardly reason to think that age-structure will alleviate the threat of collapse. However, this remains to be investigated.
Fisheries-induced selection is most often considered sizeselective mortality [12] directly targeting specific phenotypes. Here, however, selection is indirect. Fisheries are to a very large extent conducted in shallow rather than deep water [26] ( Figure  S5) and, therefore, fishing mortality will generate selection against genotypes adapted to shallow water although there is no direct targeting of specific phenotypes. There are general lessons in this for population and conservation genetics that changes in habitat can lead to intense selection even if the mortality is non-selective in the habitat in which it occurs. Sequence variation studies show that the Pan I A and B alleles are ancient predating the split of Atlantic cod from its sister species, the Pacific walleye pollock Theragra chalcogramma [15,16]. Such studies provide a deep window into the species past evolutionary and selective history. However, here we observe a steep allele frequency gradient with depth and intense fisheries-induced selection. These factors clearly would influence measures of population differentiation such as F ST . We, therefore, question the use of Pan I and other loci under such intense selection as markers for analysis of population breeding structure due to reduced gene flow [18,19], particularly if depth, the confounding of depth and geographic location, and habitatspecific fishing pressures are not controlled for. In particular, we question the practice of combining results from strongly selected loci such as Pan I with variation of supposedly neutral microsatellite loci to study breeding structure [18,19] due to population isolation. The strongly selected loci will drive the overall measure of differentiation even in the face of considerable gene flow. Fisheries-induced selection can cause differentiation of  Table 6; top right panel based on fitness estimates from Supplementary Table S2; lower left panel based on highest fitness estimates from Table 6 (a best case scenario); lower right panel based on lowest fitness estimates from Table 6 (a worst case scenario). Starting frequency of 0.738 assumed based on intercept of gam fit in Figure 6. Years based on a generation time of 4.8 years [36]. Color codes are black for the A allele and red, magenta, and blue for the AA, AB, and BB genotypes respectively. doi:10.1371/journal.pone.0005529.g007 Pan I among local groups and using that differentiation to argue for local population structure and special management of local populations may be circular.
There is a great need to understand selection and local molecular adaptation and how fisheries can indirectly and inadvertently generate intense selection as in this case. Conservation measures such as conserving large fish [13] would not be enough if they failed to protect genetic variants adapted to local niches.

Fisheries-induced changes and evolution of life-history traits
Data on catch per unit effort are difficult to interpret as indicators of population abundance [28]. For the Icelandic cod fishery total catch is highly correlated with TAC and increases and decreases with it. Total catch always exceeds the TAC and thus the fishing mortality is greater than assumed with the recommended TAC. As expected fishing effort increases with catch and TAC and so does catch per unit effort. As catch diminishes for some gear, such as gill net, effort is also reduced and at a low total catch the catch per unit effort may increase again. Other gear behave differently. For long line in particular, which lately is taking a larger share of the total catch, both catch and effort have increased and stayed high. There are signs of changes in catchability [31,50] and increased fishing pressure and effort as has been observed in the decline and collapse of other cod populations [11,31].
Rapid changes in maturation preceded the collapse of the northern cod at Newfoundland [6]. Similarly we have indications of maturation changes occurring in Icelandic cod. The caveat is that we have not studied differential maturation of the sexes or potential effects of geographic location (e.g. north vs. south [19]) or environment and we are using averages as data. We interpret the trends in estimated probabilistic maturation reaction norms to mean that the sigmoid (logistic) maturation curves on length (and/ or age) are changing. As their inflection points are pushed towards shorter lengths or lower ages they are also changing shape, becoming a stepped function. Probabilistic maturation reaction norms are useful for assessing genetic changes in the presence of environmental variation and phenotypic plasticity [32]. We, therefore, hypothesize that these are selective changes. Small and young fish may be evolving to delay reproduction while larger and older fish evolve to mature earlier and the fish become mature in a narrower window of length (or age). Overall, therefore, we have signatures for Icelandic cod of changes in effort and in lifehistory traits that are comparable to changes observed in the collapse of other cod stocks.
Olsen et al. [6] cautioned that ''Although eroding maturation reaction norms can thus signal extreme exploitation pressures, they are not to be misinterpreted as signs of imminent stock collapse. But exploitation pressures so strong that they overturn a species' natural pattern of life-history adaptation certainly ought to be cause for concern.'' Our study certainly appears to meet the criteria for concern. The strength of selection imposed by the fishery in Iceland is extremely high with selection coefficients of 92% and 73% against AA and AB genotypes respectively. This is in the high end in the distribution of known selection coefficients [51,52]. The life-history changes coupled with the Pan I changes are perhaps even more dramatic than that documented for the northern cod [6].

Future of population and fishery
We hypothesize that a collapse of the fishery is imminent if Pan I genotypic frequencies change as predicted. This hypothesis is supported also by changes in life-history. Considering fate of the fishery, deep-water fish are harder to catch and as they increase in frequency the fishery may become commercially in-viable. Fishing mortality in the preferred habitat would then cease before exterminating the A allele and fitness would revert back to natural values. Under that scenario selection pressures will diminish and the A allele may not go extinct completely which presumably would help subsequent recovery [6]. An alternative, and more likely, scenario is contraction of habitat use from tertiary and secondary to the most suitable primary habitat as a densitydependent response similar to the collapse of the northern cod stock in Newfoundland [53] and of the North Sea cod [54]. Fishers will go after smaller and smaller but equally dense clusters of fish in the primary habitat that still allow a profitable commercial fishery until the shallow water fish disappear from all habitats. Fisheriesinduced selection at the Pan I locus may have contributed to the collapse of the northern cod and other threatened cod stocks. If fishing mortality causes large decreases in frequencies of the A allele it is unlikely to revert quickly back to previous values after fishing ceases. The fact that we see rapid changes in frequency means that back selection for A from natural causes clearly is much lower than the intense selection against A caused by the fishing mortality. Also, the lack of changes among fish greater than eight years of age ( Figure 6) shows no evidence for back selection. Therefore, upon collapse the fishery would take a longer time to recover than it takes to collapse [11]. With AA and AB fish decimated by fishery we can inquire whether BB fish could invade the shallow-water niche and support a commercial fishery. This is unlikely as BB fish are deep-water adapted types [46]. We have found no evidence for historical separation or population structure by depth. However, to the extent that such a structure exists with limitation on interbreeding between shallow water and deep water fish a collapse would be more rapid and its effects would be harder to reverse.
Current management strategies of Icelandic cod are based on a catch limitation system where each vessel is allowed a certain share of the total allowable catch (TAC). Annual TACs are based on scientific assessment of state of fish stocks and ecosystem condition, but have been reduced lately. In addition, special measures for protecting small fish and the ecosystem are implemented. Thus relevant areas may be closed for short periods, if the percentage of small fish or by-catch exceed set limits. Additionally, cod spawning grounds are closed annually at the height of the spawning season in March and April to protect spawning fish. Apparently, however, these measures do not protect shallow-water fish as selection is similarly affecting all year classes that entered our samples ( Figure 5). Furthermore, some current management measures may actually intensify selection. For example, during the March/ April stop, deep-water BB fish move to shallow waters to spawn [46] (Table 1) whereupon they return to deep water and relative safety from fishing mortality. On re-opening, fishing mortality will hit the AA and AB fish that stay in shallow-water. Thus without knowledge of local adaptation good conservation intentions may exacerbate the problem.

Averting collapse with no-take reserves
We consider that our study meets criteria for concern that the Icelandic cod stock is imperiled. Can anything be done to avert collapse? Upon collapse of the northern cod of Newfoundland the Canadian government imposed a moratorium on fishing [11]. Such a drastic measure if imposed in Iceland doubtless would avert collapse. Alternatively management measures that shifted fishing from shallow-water to deep-water or measures that distributed fishing effort evenly over all depth ranges by controlling fishing by different gear also could possibly help avert collapse. However, we consider that such strategies would be difficult to implement. Alternatively we speculate and suggest that it may be possible to avert collapse by adopting a different strategy of removing selection pressures against shallow-water adapted AA and AB fish. This highlights the use of evolutionary thinking for management and conservation issues. Given that current practices are ineffective in protecting shallow-water adapted fish, we suggest that immediate action is required. We suggest that establishment of large no-take marine reserves that range from the shoreline down to the very deep waters of at least 500 meters or more would protect all genotypes. In the case of Icelandic cod an obvious area is Selvogsbanki and Faxafloi, the main spawning grounds in the Southwest [26,55]. Additional areas would be the shallow-waters in the Northeast which were closed for some years with good results but subsequently re-opened [56]. The advantage of no-take reserves would be to relieve selection pressures against the shallowwater adapted AA and AB fish. Although there are gaps in our knowledge of no-take reserves [57,58] we predict standard benefits of spillover of adults from prime into secondary and tertiary habitat and export of pelagic eggs and larvae that will ultimately benefit the fishery [59,60].

Sampling and measurement
To assess temporal and spatial variation in Pan I frequencies we sampled cod measured and aged at all predetermined sampling stations during the Marine Research Institute spring spawning surveys in 2005, 2006, and 2007 in eight of nine divisions revised from definitions in the METACOD project [19,61] (Figure S2). At each station a set with 12 gill nets of alternating six to nine inch mesh size made with mono-filament or multi-filament yarn was laid out. The sets stayed in place for at least 12 hours. Each net was 50 meters long and a set of 12 thus was 600 m long. The height of nets was 50-60 meshes or about 12-15 m. In the steepslope deep waters off the South coast ( Figure S2) a double set of 24 gill nets 1200 m long was laid out. A set of gill nets thus could cover a range of space and depths. We used mean location and mean depth in the analysis. We similarly took stratified random samples of stations taken during the MRI fall ground fish surveys in 2004, 2005, and 2006.
From each net up to 25 fish were taken for measurement of various individual traits. Otoliths were taken for age determination from a single fish from each net and a sample of gill tissue was taken for genetic analysis from these and preserved in 96% ethanol. The year class (cohort) of a fish was determined from the sample year and age read from otolith. Based on our sampling design most year classes entered our sampling for three consecutive ages except the very recent and old year classes which entered for one or two ages.
Commercial catches and effort for 1997-2007 by different gear were obtained from logbook data and from official statistics. Effort of bottom trawl is trawling time in minutes, effort of Danish seine is number of throws, effort of gill net is number of sets of nets, effort of hand line is number of hours at sea, and effort of long line is number of lines times number of hooks per line. Log book results on effort and catch per unit effort for most of the gear have also been presented in figures 9.3.1 and 9.3.2 in [62].

Molecular analysis
The Pan I locus of Atlantic cod has two alleles, A and B, defined by the presence or absence of a DraI restriction site [63]. We used a proteinase K digestion/chelex 100 method [64] for DNA isolation from tissue. We used primers 3 and 20 [15] to amplify a 489 base pair fragment of the Pan I gene and digested that with DraI to reveal diagnostic bands of the three genotypes on an agarose gel [63]. Altogether we genotyped over 8100 individuals.

Estimating probabilistic maturation reaction norms
To estimate probabilistic maturation reaction norms [32] we used methods for estimating probabilities of maturing from maturity ogives [33,34]. We used data on sexual maturity at age in Marine Research Institute spring surveys (Table 3.1.4 in [35]) and mean length and age. We estimated maturity ogives o, probability of being mature at mean length and age, with a generalized linear model (glm) logistic regression: logit½o(a : l l)~a z b(a : l l) with sample sizes as weights. The data are mean maturity and mean length l¯at age a and thus fitting a full model [34] is not possible. Instead the interaction term a:lt akes into account potential non-linearity of age and length. Following Barot et al. [33,34]  To parameterize the reaction norms we fitted a logistic model of logit(m) on mean length or age and calculated the quantiles of length and age at 5%, 50% and 95% m using the dose.p function of the MASS library under R [65].

Statistics and fitness estimation
We mostly used R [66] and various in house functions and packages under R for statistical and genetic analysis. In particular we used the LATTICE package [67,68] that implements Trellis graphics in R, the MGCV package [69,70] for fitting generalized additive models (gam) and the HIERFSTAT package [71] for R that implements an algorithm [72] for estimating F-statistics at any level of a nested or hierarchical structure.
Fitness is considered a weight, U, that transforms a genotypic frequency at one age into a frequency at another and higher age. To estimate the weight we took ratios of genotypes frequencies at two ages. We used the generalized additive model (gam) fit ( Figure 6), with a quasibinomial link function to model overdispersion, to predict genotypic frequencies and approximate 95% confidence intervals at 4 and 8 years of age (Table 6). We took the ratio of these predictions to estimate the weights. Furthermore, we took the ratios of upper to lower and lower to upper confidence intervals for predictions of best-case and worst-case scenarios. Relative fitnesses, W i , are the weights scaled to the most fit Pan I BB genotype.
In supplementary material (Table S2) we also took the pooled observed frequencies among 3 and 4 year old (entering the fishery; ''pre-selection'') and among 8-13 year old (ages at which frequencies do not change much and the brunt of the fishery is over; ''post-selection'') to estimate the weights and standard errors based on variance of ratios. In Table S3 we used median genotypic frequencies at age within the various year classes to estimate weights as yearly transitions. Assuming independent action of weights in time we multiplied the yearly weights to get an overall weight.
We used our estimated fitness values to plug into a constantfitness viability model [73] and assuming non-overlapping generations [49] Dp A~p q p W AA { W AB ð Þ z q W AB { W BB ð Þ ½ W for predicting allele and genotypic frequency changes at Pan I. The effect of selection and age-structure [48,49] was not considered here.

Supporting Information
Figure S1 Generalized additive model smooth fit of frequency of A allele on depth of sampling. Shaded region represents two standard errors above and below fit; points are partial residuals; the gam smooth was fitted using a quasibinomial link function to model overdispersion; estimated degrees of freedom edf = 3.2; approximate significance of smooth term