Effects of fishing on growth traits: a simulation analysis.AbstractFisheries often target individuals based on size. Sizeselective fishing can create selection differentials on lifehistory traits and, when those traits have a genetic basis, may cause evolution. The evolution of lifehistory traits affects potential yield and sustainability of fishing, and it is therefore an issue for fishery management. Yet fishery managers usually disregard the possibility of evolution, because little guidance is available to predict evolutionary consequences of management strategies. We attempt to provide some generic guidance. We develop an individualbased model of a population with overlapping generations
selectivity 1. curve, age at 50% selectivity, and duration of fishing season. We found that each characteristic affected the magnitude of selection differentials. The most vulnerable stocks were those with a short spawning or fishing season. Under almost all lifehistory and fishery characteristics examined, selection differentials created by realistic fishing mortality rates are considerable. ********** Fishing is typically size selective. It almost always targets the larger individuals of a population and can thus shift the spawning stock towards smaller, slowergrowing individuals. If somatic somatic /so·mat·ic/ (somat´ik) 1. pertaining to or characteristic of the soma or body. 2. pertaining to the body wall in contrast to the viscera. so·mat·ic adj. growth has some genetic basis, sizeselective fishing may cause evolution toward a smaller sizeatage. Changes in somatic growth are well documented in field data, and several studies implicate im·pli·cate tr.v. im·pli·cat·ed, im·pli·cat·ing, im·pli·cates 1. To involve or connect intimately or incriminatingly: evidence that implicates others in the plot. 2. fishing (Ricker, 1981; Harris and McGovern, 1997; Haugen and Vollestad, 2001; Sinclair et al., 2002). However, with typical field data, it is difficult to rule out other explanations. Changes in growth could result from fluctuations in population density or the environment. Furthermore, they may not be evolutionary, but instead expressions of phenotypic phe·no·type n. 1. a. The observable physical or biochemical characteristics of an organism, as determined by both genetic makeup and environmental influences. b. variability. Because of such possibilities, the idea that fishing can cause evolution has often been accepted because of compelling theoretical arguments, rather than on empirical support. However, the laboratory experiments of Conover and Munch (2002) demonstrated that size selection can cause evolution of growth traits. More and more, fishinginduced evolution is considered not just possible, but prevalent (Law, 2000; Stockwell et al., 2003). The evolution of growth traits, despite wide acknowledgement of the potential for evolution of these traits, is usually a low priority in fishery management. However, it raises at least four management concerns. First, any reduction in growth rate or maximum size can decrease recreational and economic value (Miller and Kapuscinski, 1994). Second, size selection could reduce genetic variability Introduction Genetic Variability
prison where former professor Farragut, who had killed his brother, witnesses the torments and chaos of the penal system. [Am. Lit.: Cheever Falconer in Weiss, 151] See : Imprisonment and Mackay, 1996), unpredictably altering correlated cor·re·late v. cor·re·lat·ed, cor·re·lat·ing, cor·re·lates v.tr. 1. To put or bring into causal, complementary, parallel, or reciprocal relation. 2. traits and population fitness. Third, evolution may not easily be reversed, even with afterthefact management. Fourth, the evolution of growth and other lifehistory traits can modify population dynamics Population dynamics is the study of marginal and longterm changes in the numbers, individual weights and age composition of individuals in one or several populations, and biological and environmental processes influencing those changes. (Bronikowski et al., 2002; Shertzer and Ellner, 2002) and therefore potential yield (Edley and Law, 1988; Heino 1998). Evolution in fishes can be rapid (Reznick et al., 1997; Hendry et al., 2000; Quinn et al., 2001), so that evolutionary, population, and fishery dynamics occur on similar timescales (Sinervo et al., 2000; Shertzer et al., 2002; Yoshida et al., 2003). These dynamics imply that evolution matters for fishery management on the timescale of years or decades. For fishing to cause evolution, two conditions must be met. There must be a selection differential on a phenotypic trait trait (trat) 1. any genetically determined characteristic; also, the condition prevailing in the heterozygous state of a recessive disorder, as the sickle cell trait. 2. a distinctive behavior pattern. and a genetic basis must exist for the trait's expression (i.e., the trait must be heritable her·i·ta·ble adj. 1. Capable of being passed from one generation to the next; hereditary. 2. Capable of inheriting or taking by inheritance. ). Selection differential is defined as the difference in the mean phenotypic trait value of parents before and after selection (e.g., sizeselective fishing). Stokes Stokes , William 18041878. British physician. Known especially for his studies of diseases of the chest and heart, he expanded on the observations of John Cheyne in describing the breathing irregularity now known as CheyneStokes respiration. and Law (2000) argued that, under exploitation levels in many of today's fisheries, "selection differentials on body size should be substantial and measurable." Even so, attempts to estimate selection differentials of actual fish stocks have been rare (but see Law and Rowell, 1993; Miller and Kapuscinski, 1994). This lack of estimates is surprising, given that the data needed are often available, as noted by Law (2001). The second necessary condition, heritability heritability /her·i·ta·bil·i·ty/ (her?itahbil´ite) the quality of being heritable; a measure of the extent to which a phenotype is influenced by the genotype. her·i·ta·bil·i·ty n. 1. , is defined as the proportion of phenotypic variability in offspring that is due to the genotypes of the parents. It can range from zero to one, with a higher value potentially speeding the evolutionary response to selection. Field estimates of heritability in fish size are uncommon because in nature it is difficult (although not impossible; McAllister et al., 1992) to separate genetic and environmental effects on phenotypes. Almost all estimates come from laboratory experiments (e.g., Hadley et al., 1991; Conover and Munch, 2002; Vandeputte et al., 2002), mostly on populations from aquaculture aquaculture, the raising and harvesting of fresh and saltwater plants and animals. The most economically important form of aquaculture is fish farming, an industry that accounts for an ever increasing share of world fisheries production. breeding programs (e.g., Gjedrem, 1983; Jarayabhand and Thavornyutikarn, 1995; Henryon et al., 2002). One might expect laboratory experiments to overestimate natural heritabilities, because experiments tend to reduce environmental effects on total phenotypic variance, but estimates from the laboratory have been similar to those from the field (Weigensberg and Roff, 1996). The laboratory experiments indicate that heritabilities in fish growth traits may vary widely among populations but are high enough to allow rapid evolution, given a large enough selection differential. Models of evolutionary response to selective harvest have usually taken one of two approaches: quantitative genetics quantitative genetics The scientific study of the statistical analysis of the effects that heredity and environment have on phenotypic variation. (e.g., Law, 1991; Ratner and Lande, 2001) or lifehistory optimization optimization Field of applied mathematics whose principles and methods are used to solve quantitative problems in disciplines including physics, biology, engineering, and economics. (e.g., Blythe and Stokes, 1999). In the present study, we take a different approach. Rather than attempt to predict evolution explicitly, we focus on selection differentials, a necessary (but not sufficient) condition for an evolutionary response. We use simulation analyses to compute To perform mathematical operations or general computer processing. For an explanation of "The 3 C's," or how the computer processes data, see computer. selection differentials caused by fishing. The simulation model is one common in fisheries. It consists of an agestructured population following von Bertalanffy growth, with fishing and reproduction modeled as continuous processes. Our goal is to compare selection differentials across a variety of lifehistory and fishery characteristics. We quantify selection differentials on growth parameters and body size. If growth traits are heritable, those lifehistory and fishery characteristics with the largest selection differentials are most likely to generate an evolutionary response. Armed with such knowledge, fishery managers can weigh potential evolutionary effects when choosing a fishing strategy. Materials and methods To compute selection differentials caused by sizeselective fishing we used an individualbased model (Fig. 1). To initialize To start anew, which typically involves clearing all or some part of memory or disk. the model, 250,000 individual phenotypes were generated. Each was assigned a set of lifehistory parameters and then duplicated. One copy entered an unfished population that experienced only natural mortality; the other copy entered a fished population that experienced both natural and fishing mortality. Growth, survival, and reproductive success Reproductive success is defined as the passing of genes onto the next generation in a way that they too can pass those genes on. In practice, this is often a tally of the number of offspring produced by an individual. of individuals were simulated with monthly time steps for a single year. At the end of the simulation, selection differentials on growth parameters were computed as the percent change between the mean values of spawners spawners see broodfish. in the two populations. [FIGURE 1 OMITTED] Model structure The model comprised three basic lifehistory functions: growth, survival, and reproduction. For each individual, size was assumed a function of age (a) and followed the von Bertalanffy model, (1) [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII ASCII or American Standard Code for Information Interchange, a set of codes used to represent letters, numbers, a few symbols, and control characters. Originally designed for teletype operations, it has found wide application in computers. .], where l(a)= the lengthatage of an individual; [L.sub.[infinity]] = the theoretical maximum length; K = the growth rate, and [t.sub.0] = the theoretical age when size would have been zero. In our study, each individual's age and size were updated at each monthly time step. Survival was computed differently for the two populations. In the unfished population, individuals survived with a probability depending only on the natural mortality rate (M/yr). In the fished population, individuals survived with a probability depending on both the natural mortality rate and the sizespecific fishing mortality rate. Size selectivity [s(l)] by the fishery increased with length according to according to prep. 1. As stated or indicated by; on the authority of: according to historians. 2. In keeping with: according to instructions. 3. the logistic lo·gis·tic also lo·gis·ti·cal adj. 1. Of or relating to symbolic logic. 2. Of or relating to logistics. [Medieval Latin logisticus, of calculation equation (2) [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII.], where [[beta].sub.s] = the slope of the selectivity curve; and [L.sub.s] = the length at 50% selectivity. The function s(l) describes the proportion of the fullyselected fishing mortality rate (F) experienced by individuals of length l. The sizespecific fishing mortality rate, therefore, is s(l)F per year. Fishing was applied over a fishing season of duration [D.sub.F]. The probability of reproduction was assumed equal to the probability of maturity [m(a)]. In the model, maturity increases with age and is independent of length. Although maturity likely relates to length through bioenergetics bioenergetics, n 1. system in which natural healing is enhanced by creating harmony between the patient's body and the natural environment. 2. , the relationship was not modeled here because it is, in general, poorly understood. Like selectivity (Eq. 2), re(a) was modeled by a logistic equation, but with a slope parameter, [[beta].sub.m], and age at 50% maturity, [A.sub.m]. In nature, values of lifehistory parameters K and [A.sub.m] are related to a stock's natural mortality rate. A higher natural mortality rate reduces the expected lifespan and consequently tends to be associated with a higher growth rate (K) and a younger age at maturity ([A.sub.m]). In the simulation, K and [A.sub.m] were related to natural mortality by lifehistory invariants (detailed later). Lifehistory invariants have a strong theoretical and empirical basis (Roff, 1984; Beverton, 1992; Charnov, 1993) and have been valuable in other fishery applications (Mangel, 1996; Charnov and Skuladottir, 2000; Frisk A term used in Criminal Law to refer to the superficial running of the hands over the body of an individual by a law enforcement agent or official in order to determine whether such individual is holding an illegal object, such as a weapon or narcotics. et al., 2001; Williams and Shertzer, 2003). Simulation To initialize the simulation, individuals were assigned at random to a cohort cohort /co·hort/ (ko´hort) 1. in epidemiology, a group of individuals sharing a common characteristic and observed over time in the group. 2. . The number of cohorts was determined as the age at which approximately 1% of the population would be expected to remain under natural mortality [ln(0.01)/M, rounded to the nearest integer integer: see number; number theory ]. Probabilities of cohort membership decayed exponentially ex·po·nen·tial adj. 1. Of or relating to an exponent. 2. Mathematics a. Containing, involving, or expressed as an exponent. b. with age according to M; the probability of the oldest cohort was adjusted to include the remaining fraction of fish (i.e., a plus group). The probabilities were scaled to sum to one, and a uniform random number was drawn to determine an individual's cohort. Next, individuals were assigned parameter values for von Bertalanffy growth. The value of [t.sub.0] was fixed at 0.5 yr. Values of [L.sub.[infinity]] and K were chosen uniquely for each individual. Following Xiao (1994), [L.sub.[infinity]] and K were assumed to follow a bivariate bi·var·i·ate adj. Mathematics Having two variables: bivariate binomial distribution. Adj. 1. normal distribution with standard deviations [[sigma].sub.L] and [[sigma].sub.K], respectively, and correlation [rho]. Finally, individuals were assigned a time step (month) within the year to attempt spawning. The time step was chosen from months distributed uniformly over a spawning season of duration, [D.sub.s]. Once assigned parameter values, each individual was duplicated. One copy entered the unfished population, the other the fished population. The populations were simulated in parallel over a single model year. The simulation iterated each individual through monthly time steps. At each step, the simulation computed growth and checked for survival and reproduction. In the unfished population, the monthly probability of survival was exp exp abbr. 1. exponent 2. exponential (M/12). In the fished population, the monthly probability of survival during the fishing season depended on natural mortality and on the sizespecific fishing mortality. For simplicity, we assumed size within a month was fixed so that that the probability of survival was exp[(M/12s([l.sub.0])F)/[D.sub.F]], where [l.sub0] was an individual's size at the beginning of the month. Outside the fishing season, only natural mortality applied. To check for survival, a uniform random number was drawn and compared to the survival probability appropriate for the population. Each individual surviving to its assigned spawning time had the opportunity to reproduce. In that case, a uniform random number was drawn and compared to the probability of reproduction. If reproduction was successful, the individual's growth parameters went into a pool of parents used to compute selection differentials. Growth parameters [L.sub.[infinity]] and K jointly determine sizeatage, and it is on these parameters that we describe selection differentials. At the end of the simulation year, we computed a selection differential on each growth parameter as the percent difference between mean trait values ([L.sub.[infinity]] or K) of the unfished and fished parents. Based on the differences in [L.sub.[infinity]] and K, we also computed upper and lower bounds This article is about order theory and lattice theory. For analysis of algorithms in computational complexity, see Big O notation. In mathematics, especially in order theory, an upper bound of a subset S of some partially ordered set (P of selection differentials on sizeatage. The bounds occur where age approaches [t.sub.0] or [infinity]. Because each population consisted of the same set of individuals at the beginning of the year, any difference in growth traits between parents at the end of the year could be attributed solely to fishing. Base model and variations We began with a base model built on parameter values chosen or computed to represent common lifehistory and fishery characteristics (Table 1). We then conducted a variety of sensitivity analyses. In the base model, the natural mortality rate (M) was set at 0.2/yr, a value common for many fish species. Sensitivity analyses used M = 0.1, 0.4, or 0.8. The value of M affects the values of K, [A.sub.m], and [L.sub.s], according to the lifehistory invariant (programming) invariant  A rule, such as the ordering of an ordered list or heap, that applies throughout the life of a data structure or procedure. Each change to the data structure must maintain the correctness of the invariant. relationships (Table 1). The relationship between M and K is often referred to as the M/K ratio. Charnov (1993) suggested a central value for fishes of M/K=l.65, which we used in the base model. Beverton (1992) examined the M/K ratio for fishes and found a range of 0.5 to 2.5. We used this range in our sensitivity analyses to examine the effect of the M/K ratio on selection differentials (Table 2). The base model treated [L.sub.[infinity]] and K as independent variables ([rho]=0, Table 1). Often these parameters are correlated. A metaanalysis by He and Stewart (2001) of 235 fish populations indicated a correlation value of 0.28. The negative correlation Noun 1. negative correlation  a correlation in which large values of one variable are associated with small values of the other; the correlation coefficient is between 0 and 1 indirect correlation could be expected from a tradeoff between growth rate (represented by K) and maximum size (represented by [L.sub.[infinity]]), as has been suggested in studies of bioenergetics (Stearns, 1992; Hutchings, 1993; Mangel, 1996). Our sensitivity analyses considered negative values of correlation that range from 0.25 to 1. With the base model, selectivity and maturity were assumed to be "knifeedge," a functional form often used in fisheries for convenience. Also, in the base model the size at 50% selectivity ([L.sub.s]) was assumed to occur at an age equal to the age at 50% maturity ([A.sub.m]). Although these fishery characteristics are common, selectivity and maturity may not be knifeedge or coincide. In sensitivity analyses, we examined different shapes of selectivity and maturity curves (Fig. 2). We also examined the affect of shifting the age at 50% selectivity from 2 to 2, in relation to the base case. This shift corresponds to a range in [L.sub.s] values from 574 to 738. For simplicity, we held F constant for these sensitivity analyses, implying constant effort but resulting in different amounts of removals. [FIGURE 2 OMITTED] Under logistic selectivity, the oldest, largest fish receive the highest rate of exploitation The rate of exploitation is a concept in Marxian political economy. It usually refers to the ratio of the total amount of unpaid labor done (surplusvalue) to the total amount of wages paid (the value of labour power). . Yet often the largest fish are unavailable to a fishery because of migration patterns or regulations (e.g., a maximum size limit). Thus our sensitivity analyses included a cap on susceptible sizes. The cap was set at 70, 80, or 90% of [L.sub.[infinity]]. Using the base model, we examined the effects of annual fishing mortality rate over values that range from F=0 to F=10/yr, which is 0 to 50 times the natural mortality rate. Fishing mortality was applied continuously throughout the year (i.e., [D.sub.F]=1). In sensitivity analyses, we examined shorter fishing seasons ranging from one to six months. The F was still an annual rate but was applied over fewer months and adjusted so that the number of fish removed was the same as when [D.sub.F]=1. For seasons shorter than a full year, fishing was assumed to occur at the beginning of the year. Like the fishing season, the duration of the spawning season was a full year in the base model ([D.su.s]=1). In sensitivity analyses, the spawning season ranged from one to six months and was assumed to occur at the end of the year. A selection differential cannot exist without phenotypic variation. The base model assumed a coefficient of variation Coefficient of Variation A measure of investment risk that defines risk as the standard deviation per unit of expected return. (CV) of 20c2 in both [L.sub.[infinity]] and K. For sensitivity analyses, combinations of 0%, 10%, and 20% CV in [L.sub.[infinity]] and K were examined for the influence of growth variability on selection differentials of [L.sub.[infinity]] and K. Results Changes in growth parameters [L.sub.[infinity]] and K affect sizeatage jointly, resulting in nonuniform selection differentials across ages (Fig. 3). The selection differentials on size are bounded by the differentials at the extreme ages, [t.sub.0] and [infinity]. At the youngest age, the selection differential on size is limited by the sum of the selection differentials on [L.sub.[infinity]] and K plus their product. (At age [t.sub.0], the selection differential on size is undefined.) As age increases, the selection differential on size increases or decreases monotonically toward an asymptote asymptote In mathematics, a line or curve that acts as the limit of another line or curve. For example, a descending curve that approaches but does not reach the horizontal axis is said to be asymptotic to that axis, which is the asymptote of the curve. , the selection differential on [L.sub.[infinity]]. Thus selection differentials on size across all ages are bounded by those at [L.sub.[infinity]] + K + [L.sub.[infinity]] K and [L.sub.[infinity]]. The selection differential on the smallest fish (age approaching [t.sub.0]) is an upper bound when the selection differential on K is positive, and a lower bound when negative. These properties are important for interpreting how selection differentials on sizeatage correspond to differentials on [L.sub.[infinity]] and K. [FIGURE 3 OMITTED] Using the base model, we computed selection differentials on [L.sub.[infinity]] and K as functions of fishing mortality, over the range F=0 to F=10/yr. The selection differentials increased with F nonlinearly, resulting in a concave Concave Property that a curve is below a straight line connecting two end points. If the curve falls above the straight line, it is called convex. relationship (Fig. 4). However for F<2.0, the relationship is nearly linear. [FIGURE 4 OMITTED] The alternative models also revealed linear relationships between selection differentials and F, for F<2.0 (figures not shown). In addition, those relationships have a zero intercept (by definition, no fishing, no selection differential). Because the relationships are (nearly) linear and have a common intercept, the rank of selection differentials among models does not change across values of F. A model that bears the highest selection differential at F=0.2 does so at F=2.0. We therefore present results of sensitivity analyses for a single value of F (F=0.8/yr), with the understanding that for other values of F (up to 2.0), magnitudes of selection differentials can be inferred and ranks among models are maintained. Increased variation in [L.sub.[infinity]] and K tended to increase the selection differentials, and interaction between the two growth parameters (Tables 2 and 3). Selection differentials on [L.sub.[infinity]] were generally larger than those on K. In the base model, the largest selection differential on each growth parameter occurred when variation in the focal parameter In mathematics, the focal parameter of a conic section is the distance from the focus to the directrix. It is denoted p or k. conic section equation focal parameter circle was highest and variation in the other parameter was zero. The selection differentials on sizeatage were largest when variation in both parameters was highest (20% CV for both [L.sub.[infinity]] and K). Lifehistory parameters The correlation ([rho]) between [L.sub.[infinity]] and K was assumed to be zero in the base model and negative in sensitivity analyses. The effect of correlation depended on variation in the growth parameters. When the CV was zero for either parameter, correlation had no effect on selection differentials (Tables 2 and 3). When the CV was positive for both, a negative correlation decreased selection differentials in relation to those from the base model (Tables 2 and 3). For decreased values of the correlation coefficient Correlation Coefficient A measure that determines the degree to which two variable's movements are associated. The correlation coefficient is calculated as: (i.e., stronger negative correlation), the percent selection differentials on K decreased, whereas the percent selection differentials on [L.sub.[infinity]] either decreased or remained constant. The percent selection differentials on the size near age [t.sub.0] ranged from 3.7% to 0.1% for values of (jargon) for values of  A common rhetorical maneuver at MIT is to use any of the canonical random numbers as placeholders for variables. "The max function takes 42 arguments, for arbitrary values of 42". "There are 69 ways to leave your lover, for 69 = 50". [rho] from 0 to 1. The percent selection differentials on [L.sub.[infinity]] remained relatively constant, ranging from 2.1% to 2.5%, with the highest at [rho]=0 (Fig. 5). [FIGURE 5 OMITTED] Knifeedge maturity ([[beta].sub.m] = [infinity]) resulted in larger selection differentials than did other maturity curves (Tables 2 and 3). As the slope of the maturity curve became more gradual, the selection differentials decreased. For [[beta].sub.m] values greater than 1, the selection differentials on size were similar to those of the knifeedge case (Fig. 5). The effect of M on selection differentials was relatively small (Tables 2 and 3). Changes in M from 0.1 to 0.8 led to small changes in selection differentials (Fig. 5). The largest selection differentials tended to occur near intermediate values of M (Tables 2 and 3, Fig. 5). This nonlinear A system in which the output is not a uniform relationship to the input. nonlinear  (Scientific computation) A property of a system whose output is not proportional to its input. response in the selection differentials is not surprising because changes in M affected the values of K, [A.sub.m], and maximum age nonlinearly (Table 1). Changes in the M/K ratio did not reveal a clear trend (Tables 2 and 3, Fig. 5). As with M, the M/K ratio affects other parameters; therefore changes in M/K could be expected to produce a nonlinear response in the selection differentials. The percent selection differential on [L.sub.[infinity]] was lowest at an intermediate value of M/ K=2 (Table 3). The percent selection differentials on K showed no consistent trend (Table 2). For M/K values from 0.5 to 2.5, the selection differentials on size across ages ranged from 2.3% to 4.0% (Fig. 5). Decreases in the spawning season duration ([D.sub.s]) caused a near linear increase in the selection differentials (Tables 2 and 3, Fig. 5). A compressed spawning duration of one month resulted in a range of 5.0% to 7.4% selection differential on size across ages (Fig. 5). Of all the lifehistory parameters examined in this analysis, spawning duration had the greatest effect. Fishery parameters A limit ([L.sub.u]) on sizes susceptible to the fishery decreased the selection differentials (Tables 2 and 3, Fig. 5). The percent selection differential at all ages was zero for [L.sub.u] = 800 and 0.1% for [L.sub.u] = 700 (Fig. 5). In these analyses, F was held constant. Consequently, smaller values of [L.sub.u] correspond to fewer fish removed. An alternative approach would have been to maintain constant catch by increasing F, which would have led to selection differentials larger than those in Tables 2 and 3. Knifeedge selectivity ([[beta].sub.s] = [infinity]) caused larger selection differentials than did selectivity curves with more gradual slopes (Tables 2 and 3). For [[beta].sub.s] greater than 0.1, the selection differential rapidly converged to that of the knifeedge case (Fig. 5). As with [L.sub.u], F was held constant across [[beta].sub.s], sensitivity analyses. A change in the ages of fishery selectivity had little effect on selection differentials (Tables 2 and 3, Fig. 5). When selectivity was set to a larger age or size, the selection differential decreased slightly. In this case, selectivity was occurring after maturity, allowing more fish to reproduce before reaching sizes selected by the fishery. However if harvest had been held constant instead of F, the selection differentials would have been larger. When selectivity was set to a smaller age or size, the selection differential decreased slightly or remained constant. This result is due to a reduction in the time exposed to differential fishing mortality. Differential fishing mortality occurs only on the sizes where selectivity is less than one; otherwise fishing mortality is constant for all individuals. Under von Bertalanffy growth, younger fish grow more quickly. A decrease in the age or size of selectivity shifts the fishing pressure to ages with quicker growth, reducing the time individuals experience differential fishing pressure and therefore the potential for selection differentials. If harvest had been held constant instead of F, the selection differentials would have been larger. The fishing season duration ([D.sub.F]) affected selection differentials in ways similar to the spawning season duration (Tables 2 and 3, Fig. 5). A fishing season of one month resulted in an upper bound of selection differentials that ranged from 4.8% to 7.3% over all ages (Fig. 5). Of all the fishery parameters examined in this analysis, a concentrated fishing season resulted in the largest selection differentials. Discussion The individualbased simulation approach used here simplifies computation of selection differentials and isolates the causefishing. Yet with any simulation analysis (language, simulation) SIMulation ANalysis  (SIMAN) A simulation language, especially for manufacturing systems, developed by C. Dennis Pegden in 1983. ["Introduction to Simulation using SIMAN", C.D. Pegden et al, McGrawHill 1990]. , one must interpret results in light of model assumptions. With our model maturity was assumed to be a function of age, and the computation of selection differentials were consequently focused to those on growth traits and size. If maturity were considered a function of size, it too would have been subject to a selection differential. Changes in size or age at maturity have been considered in other studies (Stokes and Blythe, 1993; Haugen and Vollestad, 2001; Olsen et al., 2004) and are likely connected to growth parameters through bioenergetic constraints. A central assumption is that somatic growth follows the von Bertalanffy model. That model was chosen because of its successful track record (Chen et al., 1992; Quinn and Deriso, 1999). Lifehistory characteristics other than growth are assumed to follow lifehistory invariant relationships. The invariants constrain con·strain tr.v. con·strained, con·strain·ing, con·strains 1. To compel by physical, moral, or circumstantial force; oblige: felt constrained to object. See Synonyms at force. 2. biological parameters to values that represent an "average stock." Of course, no stock is truly average, and therefore our sensitivity analyses incorporate considerable deviation from lifehistory invariants. In our simulation, the largest selection differentials occurred when the spawning or fishing seasons were compressed. We modeled fishing seasons at the beginning of the year and spawning seasons at the end of the year, and in a singleyear simulation, the annual timing of the fishing and spawning seasons will affect selection differentials. For example, if the onemonth fishing season had been modeled at the end of the year, the selection differential would be smaller because of the 11 months of spawning prior to fishing mortality. Over multiple years, however, the annual timing of the fishing and spawning seasons is less important than their duration and overlap. Our model simulated selection differentials at the onset of a fishery. As a fishery progresses, selection differentials should decrease as lifehistory parameters shift in the direction of selection. A multiyear simulation of evolution would require knowledge or assumptions about heritability and trait distributions, both of which are likely to be dynamic. Even so, a shortterm simulation, where selection differentials and heritability are assumed to be static, may be an informative approximation approximation /ap·prox·i·ma·tion/ (ahprok?sima´shun) 1. the act or process of bringing into proximity or apposition. 2. a numerical value of limited accuracy. . We simulated evolution of the basemodel population, assuming a static heritability of 0.2 and selection differentials of 2.5% for [L.sub.[infinity]] and 1.2% for K (values from Tables 2 and 3 with 20% CV's in both parameters). Two simulations were conducted with different values for fishing mortality. With F = 4M, five years of evolution led to a 9.0% decrease in the capacity of spawning biomass. With F = M, five years led to a 2.3% decrease. With real fishery data it is often impossible to document conclusively con·clu·sive adj. Serving to put an end to doubt, question, or uncertainty; decisive. See Synonyms at decisive. con·clusive·ly adv. that fishing causes a genetic change in growth. Any such change may be hard to measure, fall within the range of statistical variability due to sampling, or be masked by strong year classes. Selection for reduced growth may be compensated by densitydependent effects (for example, lower abundance leaving more resources for survivors to allocate towards growth). Even when a change can be demonstrated, fishing is just one potential explanation. Alternative explanations include environmentally driven evolution and reaction norms (i.e., phenotypic expressions of a genotypeenvironment interaction). Nonetheless, sizeselective fishing is widespread and often accompanies changes in somatic growth rates Growth Rates The compounded annualized rate of growth of a company's revenues, earnings, dividends, or other figures. Notes: Remember, historically high growth rates don't always mean a high rate of growth looking into the future. (Ricker, 1981; Harris and McGovern, 1997; Haugen and Vollestad, 2001; Sinclair et al., 2002). Until recently, the question was whether fishing can cause changes in growth that are evolutionary, and the answer was "yes ... probably." The laboratory experiments of Conover and Munch (2002) removed any doubt. However, those experiments represented an extreme fishery in terms of its potential to inflict a selection differential: high F compressed in time (90% of population removed in one day), knifeedge selectivity, nonoverlapping generations, and a population where all individuals are susceptible. The goal of our study was to shed light on selection differentials created by fishing under realistic ranges of lifehistory and fishery characteristics. Understanding how lifehistory characteristics affect selection differentials is important for identifying which stocks are most susceptible to evolution of growth traits. For example, susceptibility susceptibility the state of being susceptible. Refers usually to infectious disease but may be to physical factors such as wetting or to psychological factors such as harassment. increases with compression of the spawning season. Fish species with compressed spawning seasons, such as many anadromous anadromous said of fish; those living most of their lives in the sea but entering rivers to spawn. species, may be at higher risk of evolution from sizeselective fisheries. Understanding how fishery patterns affect selection differentials has direct management implications because it is the fishery parameters that can be controlled. For example, our results indicate that sizeselective fisheries compressed in time are apt to cause high selection differentials. Managers should avoid "derby" style harvests, such as the annual Pacific herring Noun 1. Pacific herring  important food fish of the northern Pacific Clupea harengus pallasii Clupea, genus Clupea  type genus of the Clupeidae: typical herrings sacroe fisheries, which are completed in only a few days. Other management strategies could reduce selection differentials, such as slot limits, reduction in the slope of selectivity curves, and partial selectivity after the age at maturity. However, because no sizeselective fishing pattern can preclude some directional selection Immunology < Genetics In population genetics, directional selection on growth, management by area closures may be the best option for avoiding fisheryinduced evolution of growth traits. As fishing technology improves, so does the ability to fully and rapidly exploit fish populations, and thus increase the potential for evolutionary responses. Still, when overfishing Overfishing occurs when fishing activities reduce fish stocks below an acceptable level. This can occur in any body of water from a pond to the oceans. More precise biological and bioeconomic terms define 'acceptable level'. depletes a stock, low abundance is usually the paramount concern. With appropriate management, stock abundance may recover, but prefishing growth capacity may recover more slowly or not at all if genetic variation is lost. Given plausible heritabilities of growth traits, this analysis shows that under a wide variety of lifehistory and fishery characteristics, selection differentials are large enough to allow for rapid evolution. Table 1 Parameter values used in the base model. Formulas for the growth rate (K) and the age at 50% maturity ([A.sub.m]) are lifehistory invariant relationships from Charnov (1993) and Beverton (1992), respectively. The formula for [L.sub.s] is the length at age [A.sub.m] according to von Bertalanffy growth. A value of [infinity] for slope parameters corresponds to a knifeedge curve. Parameter Description M Natural mortality rate (per year) F Fishing mortality rate (per year) [[bar.L].sub.[infinity]] Mean asymptotic size in growth function [bar.K] Mean growth rate in growth function [t.sub.0] Location parameter in growth function C[V.sub.L] Coefficient of variation in [L.sub.[infinity]] C[V.sub.K] Coefficient of variation in K [rho] Correlation between [L.sub.[infinity]] and K [[beta].sub.s] Slope of the size selectivity curve [[beta].sub.m] Slope of the maturity curve [A.sub.m] Age at 50% maturity [L.sub.s] Length at 50% selectivity [D.sub.S] Duration of spawning season (yr) [D.sub.F] Duration of fishing season (yr) Parameter Formula Value M Fixed 0.2 F Fixed 0 to 10 [[bar.L].sub.[infinity]] Fixed 1000 [bar.K] M/1.65 0.12 [t.sub.0] Fixed 0.5 C[V.sub.L] Fixed 20% C[V.sub.K] Fixed 20% [rho] Fixed 0 [[beta].sub.s] Fixed [infinity] [[beta].sub.m] Fixed [infinity] [A.sub.m] log[(3[bar.K] + M)/M]/ [bar.K] 8.55 [L.sub.s] [[bar.L].sub.[infinity]][1  666 exp([bar.K][[A.sub.m]  [t.sub.o]])] [D.sub.S] Fixed 1 [D.sub.F] Fixed 1 Table 2 Percent selection differential on the von Bertalanffy growth coefficient (K) at fishing mortality = 0.8/yr. Columns correspond to the levels of the coefficient of variation (CV=0%, 10%, 20%) in K and in the asymptotic length ([L.sub.[infinity]]). Any combination with 0% CV in K is not presented because it results in zero selection differential. The first row corresponds to the base model and subsequent rows correspond to changes in the base model: correlation between [L.sub.[infinity]] and K ([rho]), slope of maturity curve ([[beta].sub.m]), natural mortality (M), M/K ratio, duration of annual spawning season ([D.sub.S]), maximum size limit ([L.sub.u]), slope of selectivity curve ([[beta].sub.s]), change in age at 50% selectivity ([A.sub.s]) in relation to the base case, and duration of annual fishing season ([D.sub.F]). Parameter values [L.sub.[infinity]]: [L.sub.[infinity]]: 0%CV 10%CV K: 10%CV K: 10%CV Base 0.7 0.5 [rho] = 1 0.7 0.7 [rho] = 0.75 0.7 0.3 [rho] = 0.5 0.7 0.0 [rho] = 0.25 0.7 0.3 [[beta].sub.m] = 0.25 0.2 0.2 [[beta].sub.m] = 0.5 0.3 0.3 [[beta].sub.m] = 1 0.5 0.4 M = 0.1 0.4 0.4 M = 0.4 0.7 0.4 M = 0.8 0.6 0.3 M/K = 0.5 0.6 0.3 M/K = 1 0.4 0.3 M/K = 2 0.5 0.4 M/K = 2.5 0.8 0.6 [D.sub.S] = 1/12 1.6 1.0 [D.sub.S] = 3/12 1.4 0.9 [D.sub.S] = 6/12 1.1 0.8 [L.sub.u] = 700 0.0 0.0 [L.sub.u] = 800 0.2 0.0 [L.sub.u] = 900 0.4 0.2 [[beta].sub.s] = 0.01 0.3 0.2 [[beta].sub.s] = 0.05 0.6 0.4 [[beta].sub.s] = 0.1 0.6 0.5 [A.sub.s] = 2 0.1 0.2 [A.sub.s] = 1 0.4 0.4 [A.sub.s] = 1 0.6 0.5 [A.sub.s] = 2 0.5 0.4 [D.sub.F] = 1/12 1.5 1.0 [D.sub.F] = 3/12 1.3 0.9 [D.sub.F] = 6/12 1.2 0.7 Parameter values [L.sub.[infinity]]: [L.sub.[infinity]]: 20%CV 0%CV K: 10%CV K: 10%CV Base 0.3 2.1 [rho] = 1 1.3 2.1 [rho] = 0.75 0.8 2.1 [rho] = 0.5 0.4 2.1 [rho] = 0.25 0.0 2.1 [[beta].sub.m] = 0.25 0.1 0.7 [[beta].sub.m] = 0.5 0.2 1.1 [[beta].sub.m] = 1 0.3 1.7 M = 0.1 0.3 1.6 M = 0.4 0.3 2.0 M = 0.8 0.2 1.6 M/K = 0.5 0.1 1.9 M/K = 1 0.2 1.5 M/K = 2 0.3 1.9 M/K = 2.5 0.4 2.4 [D.sub.S] = 1/12 0.6 4.5 [D.sub.S] = 3/12 0.5 4.1 [D.sub.S] = 6/12 0.5 3.3 [L.sub.u] = 700 0.0 0.1 [L.sub.u] = 800 0.0 0.3 [L.sub.u] = 900 0.1 1.1 [[beta].sub.s] = 0.01 0.2 1.1 [[beta].sub.s] = 0.05 0.3 1.9 [[beta].sub.s] = 0.1 0.3 2.0 [A.sub.s] = 2 0.2 1.1 [A.sub.s] = 1 0.3 1.7 [A.sub.s] = 1 0.3 2.1 [A.sub.s] = 2 0.3 1.9 [D.sub.F] = 1/12 0.6 4.3 [D.sub.F] = 3/12 0.6 3.9 [D.sub.F] = 6/12 0.5 3.3 Parameter values [L.sub.[infinity]]: [L.sub.[infinity]]: 0%CV 20%CV K:20%CV K:20%CV Base 1.7 1.2 [rho] = 1 0.2 2.3 [rho] = 0.75 0.8 0.9 [rho] = 0.5 1.1 0.1 [rho] = 0.25 1.4 0.6 [[beta].sub.m] = 0.25 0.7 0.6 [[beta].sub.m] = 0.5 1.1 0.9 [[beta].sub.m] = 1 1.4 1.1 M = 0.1 1.5 1.2 M = 0.4 1.5 1.0 M = 0.8 1.1 0.7 M/K = 0.5 1.0 0.6 M/K = 1 1.3 0.8 M/K = 2 1.6 1.2 M/K = 2.5 2.1 1.5 [D.sub.S] = 1/12 3.6 2.3 [D.sub.S] = 3/12 3.3 2.2 [D.sub.S] = 6/12 2.7 1.8 [L.sub.u] = 700 0.1 0.0 [L.sub.u] = 800 0.2 0.0 [L.sub.u] = 900 0.8 0.3 [[beta].sub.s] = 0.01 1.0 0.8 [[beta].sub.s] = 0.05 1.6 1.2 [[beta].sub.s] = 0.1 1.7 1.2 [A.sub.s] = 2 1.2 1.0 [A.sub.s] = 1 1.6 1.1 [A.sub.s] = 1 1.7 1.2 [A.sub.s] = 2 1.6 1.1 [D.sub.F] = 1/12 3.5 2.3 [D.sub.F] = 3/12 3.1 2.1 [D.sub.F] = 6/12 2.6 1.8 Table 3 Percent selection differential on the von Bertalanffy asymptotic length ([L.sub.[infinity]]) at fishing mortality = 0.8/yr. Columns correspond to the levels of the coefficient of variation (CV=O%, 10%, 20%) in [L.sub.[infinity]] and in the growth coefficient (K). Any combination with 0% CV in [L.sub.[infinity]] is not presented because it results in zero selection differential. The first row corresponds to the base model and subsequent rows correspond to changes in the base model: correlation between [L.sub.[infinity]] and K ([rho]), slope of maturity curve ([[beta].sub.m]), natural mortality (M), M/K ratio, duration of annual spawning season ([D.sub.S]), maximum size limit ([L.sub.u]), slope of selectivity curve ([[beta].sub.s]), change in age at 50% selectivity ([A.sub.s]) in relation to the base case, and duration of annual fishing season ([D.sub.F]). Parameter values [L.sub.[infinity]]: [L.sub.[infinity]]: 0%CV 10%CV K: 10%CV K: 10%CV Base 1.0 2.7 [rho] = 1 1.0 2.8 [rho] = 0.75 1.0 2.7 [rho] = 0.5 1.0 2.8 [rho] = 0.25 1.0 2.8 [[beta].sub.m] = 0.25 0.3 1.2 [[beta].sub.m] = 0.5 0.5 1.8 [[beta].sub.m] = 1 0.8 2.4 M = 0.1 0.8 2.6 M = 0.4 1.0 2.8 M = 0.8 1.0 2.6 M/K = 0.5 1.4 3.2 M/K = 1 0.9 2.7 M/K = 2 0.9 2.5 M/K = 2.5 1.1 2.8 [D.sub.S] = 1/12 2.2 5.6 [D.sub.S] = 3/12 2.0 5.1 [D.sub.S] = 6/12 1.6 4.4 [L.sub.u] = 700 0.1 0.1 [L.sub.u] = 800 0.0 0.1 [L.sub.u] = 900 0.3 0.5 [[beta].sub.s] = 0.01 0.5 1.9 [[beta].sub.s] = 0.05 0.9 2.7 [[beta].sub.s] = 0.1 1.0 2.7 [A.sub.s] = 2 0.4 2.1 [A.sub.s] = 1 0.7 2.5 [A.sub.s] = 1 1.0 2.8 [A.sub.s] = 2 1.0 2.7 [D.sub.F] = 1/12 2.2 5.5 [D.sub.F] = 3/12 1.9 4.9 [D.sub.F] = 6/12 1.6 4.2 Parameter values [L.sub.[infinity]]: [L.sub.[infinity]]: 20%CV 0%CV K: 10%CV K: 10%CV Base 0.9 2.7 [rho] = 1 0.7 2.7 [rho] = 0.75 0.7 2.6 [rho] = 0.5 0.8 2.7 [rho] = 0.25 0.9 2.7 [[beta].sub.m] = 0.25 0.3 1.1 [[beta].sub.m] = 0.5 0.5 1.8 [[beta].sub.m] = 1 0.7 2.4 M = 0.1 0.7 2.5 M = 0.4 1.0 2.7 M = 0.8 1.0 2.6 M/K = 0.5 1.4 3.2 M/K = 1 0.9 2.7 M/K = 2 0.8 2.5 M/K = 2.5 1.0 2.7 [D.sub.S] = 1/12 2.0 5.5 [D.sub.S] = 3/12 1.8 5.0 [D.sub.S] = 6/12 1.5 4.3 [L.sub.u] = 700 0.1 0.1 [L.sub.u] = 800 0.0 0.1 [L.sub.u] = 900 0.2 0.5 [[beta].sub.s] = 0.01 0.5 1.9 [[beta].sub.s] = 0.05 0.9 2.6 [[beta].sub.s] = 0.1 0.9 2.7 [A.sub.s] = 2 0.4 2.1 [A.sub.s] = 1 0.7 2.5 [A.sub.s] = 1 1.0 2.7 [A.sub.s] = 2 0.9 2.6 [D.sub.F] = 1/12 2.0 5.3 [D.sub.F] = 3/12 1.7 4.8 [D.sub.F] = 6/12 1.5 4.1 Parameter values [L.sub.[infinity]]: [L.sub.[infinity]]: 10%CV 20%CV K: 10%CV K: 10%CV Base 0.8 2.5 [rho] = 1 0.1 2.3 [rho] = 0.75 0.2 2.2 [rho] = 0.5 0.4 2.3 [rho] = 0.25 0.6 2.4 [[beta].sub.m] = 0.25 0.3 1.1 [[beta].sub.m] = 0.5 0.5 1.7 [[beta].sub.m] = 1 0.7 2.2 M = 0.1 0.7 2.4 M = 0.4 0.8 2.6 M = 0.8 0.8 2.4 M/K = 0.5 1.3 3.1 M/K = 1 0.8 2.5 M/K = 2 0.7 2.3 M/K = 2.5 0.8 2.5 [D.sub.S] = 1/12 1.7 5.0 [D.sub.S] = 3/12 1.5 4.6 [D.sub.S] = 6/12 1.3 3.9 [L.sub.u] = 700 0.1 0.1 [L.sub.u] = 800 0.1 0.1 [L.sub.u] = 900 0.1 0.4 [[beta].sub.s] = 0.01 0.5 1.8 [[beta].sub.s] = 0.05 0.8 2.5 [[beta].sub.s] = 0.1 0.8 2.5 [A.sub.s] = 2 0.4 2.1 [A.sub.s] = 1 0.7 2.4 [A.sub.s] = 1 0.9 2.5 [A.sub.s] = 2 0.8 2.4 [D.sub.F] = 1/12 1.6 4.8 [D.sub.F] = 3/12 1.5 4.4 [D.sub.F] = 6/12 1.2 3.7 Acknowledgments We thank R. Munoz, M. Prager, and D. Vaughan for comments on the manuscript. This work was supported by the National Marine Fisheries Service through its Southeast Fisheries Science Center. Literature cited Beverton, R. J. H. 1992. Patterns of reproductive strategy parameters in some marine teleost teleost fish of the class Osteichthyes, having the skeleton completely ossified. fishes. J. Fish. Biol. 41(suppl. B):137160. Blythe, S. P., and T. K. Stokes. 1993. Sizeselective haresting and atatmaturity. I: Some theoretical implications for management of evolving resources. In The exploitation of evolving resources (T. M. Stokes, J. M. McGlade, and R. Law, eds.), p. 232247. Lecture Notes in Biomathematics bi·o·math·e·mat·ics n. (used with a sing. verb) The application of mathematical principles to biological processes. bi 99. SpringerVerlag, Berlin. Bronikowski, A. M., M. E. Clark, F. H. Rodd, and D. N. Reznick. 2002. Populationdynamic consequences of predatorinduced life history variation in the guppy (Poecilia reticulata). Ecology 83:21942204. Charnov, E. L. 1993. Life history invariants: some explorations of symmetry symmetry, generally speaking, a balance or correspondence between various parts of an object; the term symmetry is used both in the arts and in the sciences. in evolutionary ecology Evolutionary ecology lies at the intersection of ecology and evolutionary biology. It approaches the study of ecology in a way that explicitly considers the evolutionary histories of species and the interactions between them. , 159 p. Oxford Univ. Press, Oxford, England. Charnov, E. L., and U. Skuladottir. 2000. Dimensionless invariants for the optimal size (age) of sex change. Evol. Ecol. Res. 2:10671071. Chen, Y., D. A. Jackson, and H. H. Harvey. 1992. A comparison of von Bertalanffy and polynomial polynomial, mathematical expression which is a finite sum, each term being a constant times a product of one or more variables raised to powers. With only one variable the general form of a polynomial is a_{0}x^{n}+a_{} functions in modeling fish growth data. Can. J. Fish. Aquat. Sci. 49:12281235. Conover, D. O., and S. B. Munch. 2002. Sustaining fisheries yields over evolutionary time scales. Science 297:946. Edley, M. T., and R. Law. 1988. Evolution of life histories and yields in experimental populations of Daphnia magna. Biol. J. Linn linn n. Scots 1. A waterfall. 2. A steep ravine. [Scottish Gaelic linne, pool, waterfall.] . Soc. 34:309326. Falconer, D. S., and T. F. C. Mackay. 1996. An introduction to quantitative genetics, 4th ed., 245 p. Longman Group Ltd., Harlow, Essex, England. Frisk, M. G., T. J. Miller, and M. J. Fogarty. 2001. Estimation and analysis of biological parameters in elasmobranch elasmobranch (ĭlăs`məbrăngk), cartilaginous fish, member of the subclass Elasmobranchii of the vertebrate class Chondrichthyes (see Chordata). This group includes sharks, skates, and rays. fishes: a comparative life history study. Can. J. Fish. Aquat. Sci. 58:969981. Gjedrem, T. 1983. Genetic variation in quantitative traits and selective breeding
Selective breeding in domesticated animals is the process of developing a cultivated breed over time. in fish and shellfish shellfish, popular name for certain edible mollusks (see Mollusca), e.g., oysters, clams, and scallops, and for certain edible crustaceans, e.g., crabs, lobsters, and shrimps. All are aquatic invertebrates with shells; they are not fish. . Aquaculture 33:5172. Hadley, N. H., R. T. Dillon, and J. J. Manzi. 1991. Realized heritability of growth rate in the hard clam Mercenaria mercenaria. Aquaculture 93:109119. Harris, P. J., and J. C. McGovern. 1997. Changes in the life history of red porgy porgy (pôr`gē), common name for members of the Sparidae, a family of smallmouthed fishes with strong teeth adapted for crushing their food of shellfish and crustaceans. , Pagrus pagrus, from the southeastern United States United States, officially United States of America, republic (2005 est. pop. 295,734,000), 3,539,227 sq mi (9,166,598 sq km), North America. The United States is the world's third largest country in population and the fourth largest country in area. , 19721994. Fish. Bull. 95:732747. Haugen, T. O., and L. A. Vollestad. 2001. A century of lifehistory evolution in grayling grayling, common name for a brilliantly colored fish belonging to the genus Thymallus, of the family Salmonidae (salmon family), and closely allied to the smelt. Graylings are found chiefly in clear, cold, fresh waters of the Northern Hemisphere. . Genetica 112113 and 475491. He, J. X., and D. J. Stewart. 2001. Age and size at first reproduction of fishes: predictive models based only on growth trajectories. Ecology 82:784791. Heino, M. 1998. Management of evolving fish stocks. Can. J. Fish. Aquat. Sci. 55:19711982. Henryon, M., A. Jokumsen, P. Berg, I. Lund, P. B. Pederson, N. J. Olesen, W. J. Slierendrecht. 2002. Genetic variation for growth rate, feed conversion efficiency, and disease resistance exists within a farmed population of rainbow trout rainbow trout Species (Oncorhynchus mykiss) of fish in the salmon family (Salmonidae) noted for spectacular leaps and hard fighting when hooked. It has been introduced from western North America to many other countries. . Aquaculture 209: 5976. Hendry, A. P., J. K Wenburg, P. Bentzen, E. C. Volk, T. P. Quinn. 2000. Rapid evolution of reproductive isolation An important concept in evolutionary biology, reproductive isolation is a category of mechanisms that prevent two or more populations from exchanging genes. The separation of the gene pools of populations, under some conditions, can lead to the genesis of distinct species. in the wild: evidence from introduced salmon. Science 290: 516518. Hutchings, J. A. 1993. Adaptive life histories effected by agespecific survival and growth rate. Ecology 74:673684. Jarayabhand, P., and M. Thavornyutikarn. 1995. Realized heritability estimation on growth rate of oyster oyster, edible bivalve mollusk found in beds in shallow, warm waters of all oceans. The shell is made up of two valves, the upper one flat and the lower convex, with variable outlines and a rough outer surface. , Saccostrea cucullata Born, 1778. Aquaculture 138:111118. Law, R. 1991. On the quantitative genetics of correlated characters under directional selection in agestructured populations. Phil. Trans. R. Soc. Lond. B 331:213223. 2000. Fishing, selection, and phenotypic evolution. ICES J. Mar. Sci. 57:659668. 2001. Phenotypic and genetic changes due to selective exploitation. In Conservation of exploited species (J. D. Reynolds, G. M. Mace, K. H. Redford, and J. G. Robinson John George Robinson was chief mechanical engineer of the Great Central Railway from 1900 to 1922. Prior to this, he designed locomotives for the Waterford and Limerick Railway (between 1884 and 1900). , eds.), p. 323342. Cambridge Univ. Press, Cambridge, England. Law, R., and C. A. Rowell. 1993. Cohortstructured populations, selection responses, and exploitation of the North Sea cod. In The exploitation of evolving resources (T. K. Stokes, J. M. McGlade, and R. Law, eds.), p. 155173. Lecture Notes in Biomathematics 99. SpringerVerlag, Berlin, Germany. Mangel, M. 1996. Life history invariants, age at maturity and the ferox trout. Evol. Ecol. 10:249263. McAllister, M. K., R. M. Peterman Pe´ter`man n. 1. A fisherman;  so called after the apostle Peter. , and D. M. Gillis. 1992. Statistical evaluation of a largescale fishing experiment designed to test for a genetic effect of sizeselective fishing on the British Columbia British Columbia, province (2001 pop. 3,907,738), 366,255 sq mi (948,600 sq km), including 6,976 sq mi (18,068 sq km) of water surface, W Canada. Geography pink salmon pink salmon Food fish (Oncorhynchus gorbuscha, family Salmonidae) of the North Pacific that constitutes half of the commercial fishery of Pacific salmon. It weighs about 4.5 lbs (2 kg) and is marked with large, irregular spots. Pink salmon often spawn on tidal flats. (Oncorhynchus gorbuscha). Can. J. Fish. Sci. 49:12941304. Miller, L. M., and A. R. Kapuscinski. 1994. Estimation of selection differentials from fish scales: a step towards evaluating genetic alteration of fish size in exploited populations. Can. J. Fish. Sci. 51:774783. Olsen, E. M., M. Heino, G. R. Lilly, M. J. Morgan, J. Brattey, B. Ernande, and U. Dieckmann. 2004. Maturation maturation /mat·u·ra·tion/ (machura´shun) 1. the process of becoming mature. 2. attainment of emotional and intellectual maturity. 3. trends indicative of rapid evolution preceded the collapse of northern cod. Nature 428: 932 935. Quinn, T. J., II, and R. B. Deriso. 1999. Quantitative fish dynamics, 542 p. Oxford Univ. Press, New York New York, state, United States New York, Middle Atlantic state of the United States. It is bordered by Vermont, Massachusetts, Connecticut, and the Atlantic Ocean (E), New Jersey and Pennsylvania (S), Lakes Erie and Ontario and the Canadian province of , NY. Quinn, T. P., M. T. Kinnison, and M. J. Unwin. 2001. Evolution of Chinook salmon chinook salmon or king salmon Prized North Pacific food and sport fish (Oncorhynchus tshawytscha) of the salmon family. The average weight is about 22 lbs (10 kg), but individuals of 50–80 lbs (22–36 kg) are not unusual. (Oncorhynchus tshawytscha) populations in New Zealand New Zealand (zē`lənd), island country (2005 est. pop. 4,035,000), 104,454 sq mi (270,534 sq km), in the S Pacific Ocean, over 1,000 mi (1,600 km) SE of Australia. The capital is Wellington; the largest city and leading port is Auckland. : pattern, rate, and process. Genetica 112113:493513. Ratner, S., and R. Lande. 2001. Demographic and evolutionary responses to selective harvesting in populations with discrete generations. Ecology 82:30933104. Reznick, D. N., F. H. Shaw, F. H. Rodd, and R. G. Shaw. 1997. Evaluation of the rate of evolution in natural populations of guppies ''This article is about an American popculture term. For the fish, see Guppy Guppies is an acronym which stands for Generation X Yuppies. The combination of the two nelogistic generational terms is used to loosely identify anyone who was in their twenties during the 1990s, (Poecilia reticulata). Science 275:19341937. Ricker, W. E. 1981. Changes in the average size and average age of Pacific salmon. Can. J. Fish. Aquat. Sci 38: 16361656. Roff, D. A. 1984. The evolution of lifehistory parameters in teleosts. Can. J. Fish. Aquat. Sci. 41:9841000. Shertzer, K. W., and S. P. Ellner. 2002. Energy storage and the evolution of population dynamics. J. Theor. Biol. 215:183200. Shertzer, K. W., S. P. Ellner, G. F. Fussmann, and N. G. Hairston Jr. 2002. Predatorprey cycles in an aquatic microcosm mi·cro·cosm n. A small, representative system having analogies to a larger system in constitution, configuration, or development: "He sees the auto industry as a microcosm of the U.S. : testing hypotheses of mechanism. J. Anim. Ecol. 71:802815. Sinclair, A. F., D. P. Swain, and J. M. Hanson. 2002. Disentangling the effects of sizeselective mortality, density, and temperature on lengthatage. Can. J. Fish. Aquat. Sci. 59:372382. Sinervo, B., E. Svensson, and T. Comendant. 2000. Density cycles and an offspring quantity and quality game driven by natural selection. Nature 406:985988. Stearns, S. C. 1992. The evolution of life histories, 379 p. Oxford Univ. Press, Oxford, England. Stockwell, C. A., A. P. Hendry, and M. T. Kinnison. 2003. Contemporary evolution meets conservation biology conservation biology n. The branch of biology that deals with the effects of humans on the environment and with the conservation of biological diversity. . Trends Ecol. Evol. 18:94101. Stokes, T. K., and S. P. Blythe. 1993. Sizeselective haresting and atatmaturity. II: Real populations and management options. In The exploitation of evolving resources (T. K. Stokes, J. M. McGlade, and R. Law (eds.), p. 232237. Lecture Notes in Biomathematics 99, SpringerVerlag, Berlin. Stokes, K., and R. Law. 2000. Fishing as an evolutionary force. Mar. Ecol. Prog. Ser. 208:307309. Vandeputte, M., E. Quillet, and B. Chevassus. 2002. Early development and survival in brown trout brown trout Prized and wary European game fish (Salmo trutta, family Salmonidae) that is favoured for food. The species includes several varieties (e.g., the Loch Leven trout of Britain). The brown trout is recognized by the lightringed black spots on its brown body. (Salmo trutta fario L.): indirect effects of selection for growth rate and estimation of genetic parameters. Aquaculture 204:435445. Weigensberg, I., and D. A. Roff. 1996. Natural heritabilities: can they be reliably estimated in the laboratory. Evolution 50:21492157. Williams, E. H., and K. W. Shertzer. 2003. Implications of lifehistory invariants for biological reference points used in fishery management. Can. J. Fish. Aquat. Sci. 60:710720. Yoshida, T., L. E. Jones, S. P. Ellner, G. F. Fussmann, and N. G. Hairson Jr. 2003. Rapid evolution drives ecological dynamics in a predatorprey system. Nature 424:303306. Xiao, Y. 1994. Von Bertalanffy growth models with variability in, and correlation between, K and [L.sub.[infinity]]. Can. J. Fish. Aquat. Sci. 51:15851590. Manuscript submitted 16 April 2004 to the Scientific Editor's Office Manuscript approved for publication 20 December 2004 by the Scientific Editor. Erik H. Williams Kyle W. Shertzer Center for Coastal Fisheries and Habitat Research 101 Pivers Island Road Beaufort, North Carolina Beaufort (pronounced "BOfurt" / IPA: ˈbo.fɚt) is a town in Carteret County, North Carolina, United States. 28516 Email address: Erik.Williams@noaa.gov 

Reader Opinion