Journal of Heredity Advance Access originally published online on December 7, 2006
Journal of Heredity 2006 97(6):603-606; doi:10.1093/jhered/esl042
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Brief Communications |
Estimation of Effective Population Size for the Long-Lived Darkblotched Rockfish Sebastes crameri
From the Department of Fisheries and Wildlife and Coastal Oregon Marine Experiment Station, Hatfield Marine Science Center, Oregon State University, 2030 SE Marine Science Drive, Newport, OR 97365 (Gomez-Uchida and Banks). Daniel Gomez-Uchida is now at Department of Biology, Dalhousie University, 1355 Oxford Street, Halifax, NS B3H 4JI, Canada
Address correspondence to Dr. D. Gomez-Uchida at the address above, or e-mail: daniel.gomez-uchida{at}lifetime.oregonstate.edu.
We report the variance effective population size (Ne) in darkblotched rockfish (Sebastes crameri) utilizing the temporal method for overlapping generations, which requires a combination of age-specific demography and genetic information from cohorts. Following calculations of age-specific survival and reproductive success from fishery data, we genotyped a sample (n = 1087) comprised by 6 cohorts (from 1995 to 2000) across 7 microsatellite loci. Our Ne estimate (
) plus 95% confidence interval was = 9157 [6495–12 215], showing that the breeding population number could be 3–4 orders of magnitude smaller than the census population size ( = 24 376 210). Our estimates resemble closely those found for fishes with similar life history, suggesting that the small 
ratio for S. crameri is most likely explained by a combination of high variance in reproductive success among individuals, genetic structure, and demographic perturbations such as historical fishing. Because small ratios have been commonly associated with potential loss of genetic variation, our estimates need careful consideration in rockfish management and conservation.
Although effective population size (Ne) is a central tenet in theoretical evolutionary biology, empirical estimates of Ne have been difficult to obtain directly from natural populations because key demographic parameters are not available throughout the entire life span of individuals (Waples 2002). The so-called temporal method overcomes this limitation by indirectly assessing the variance Ne from changes in allele frequencies between discrete generations (Waples 1989). Given that Ne represents the breeding population number and determines the intensity of genetic drift in a population (Hedrick 2000), genetic estimates have become crucial information in evaluating conservation risks of natural populations (Frankham 1995). For example, loss of genetic variation, inbreeding depression, and/or greater extinction threat have been commonly associated with small Ne (Lande 1988) and small Ne/N ratios (N = census population size) (Mace and Lande 1991).
One methodological refinement of the temporal method in conservation genetics has been the combination of demographic and genetic information when generations overlap. Essentially, Ne is estimated between adjacent year classes (cohorts) and corrected for overlapping generations using age-specific survivorship and reproductive success estimates (Jorde and Ryman 1995). This approach is suitable for ecologically and economically important long-lived marine species with multiple reproductive events (iteroparity). In addition, a long generation time complicates temporal sampling of allele frequencies among discrete generations, unless historically archived samples spanning several decades are analyzed (e.g., Hauser et al. 2002).
Here we combine age-specific demography and microsatellite DNA data from cohorts of the long-lived marine darkblotched rockfish (Sebastes crameri) to estimate variance effective population size. Darkblotched rockfish is one deep-water representative of the speciose genus Sebastes in the northeastern Pacific. Sebastes crameri is characterized by long life span (maximum age recorded 105 years; Love et al. 2002) and late maturity (males: 5 years, females: 8 years; Nichol and Pikitch 1994). Low levels of biomass motivated the Pacific Fishery Management Council to declare darkblotched rockfish overfished (Rogers et al. 2000). Our goal was to provide estimates of Ne and N and interpret their relevance to rebuilding plans (Rogers 2003) and conservation initiatives for the recovery of heavily exploited rockfish (Berkeley, Hixon, et al. 2004).
 |
Materials and Methods |
|---|
 Top Materials and Methods Results and DiscussionReferences |
|---|
Samples
Rockfish fin clips (n = 1206) came from an earlier study that analyzed the genetic structure of S. crameri ranging from the coast of Washington to northern California during 2001 and 2002 (Gomez-Uchida and Banks 2005). Ages of sampled fish were determined from otoliths using the break-and-burn technique (Chilton and Beamish 1982).
Molecular Analyses
Freshly extracted DNA from each individual was screened for 7 polymorphic microsatellite loci using polymerase chain reaction (PCR) as presented in Gomez-Uchida and Banks (2005). Extra amplification cycles (e.g., 5–10 cycles) and an additive containing 20 mg/ml of bovine serum albumin were necessary when PCR fragments were weak. Details on the electrophoresis procedures and sizing of PCR fragments can be found in Gomez-Uchida and Banks (2005).
Genetic Data Analyses
Multilocus genotypes were apportioned into cohorts that matched calendar years. In order to provide a global estimate of Ne (e.g., Turner et al. 2002), we pooled genotypes across sampling locations within year classes. Only genetic information from the 6 largest cohorts (1995–2000) (n = 1087) was used during the rest of the analyses. These year classes were homogeneously distributed across the geographic area (Gomez-Uchida D and Banks MA, unpublished data), suggesting no confounding effects of geographic and intercohort differentiation.
Allele frequencies and general statistics (e.g., allelic richness, and expected vs. observed heterozygosity per locus) for each cohort were estimated using GENETIX version 3.3 (Belkhir et al. 2002). The null hypotheses of 1) Hardy–Weinberg equilibrium (HWE) within cohorts and 2) gametic (linkage) equilibrium among loci within cohorts were tested through exact tests implemented in GENEPOP version 3.3 (Raymond and Rousset 1995). The probability of the joint hypothesis over multiple loci was calculated through Fisher's method implemented in GENEPOP. Departures from HWE were further investigated with MICROCHECKER (van Oosterhout et al. 2004) to exclude potential genotyping errors. Finally, selective neutrality of alleles was tested through the X2-statistic between adjacent year classes of Lewontin and Krakauer (1973) (Jorde and Ryman 1996).
Demographic Analyses and Estimation of Ne and N
Age-specific survival rates (li) were estimated from annual survival (Si) calculated for 38 age classes found in stock assessments of S. crameri (Rogers et al. 2000; Rogers 2003). We used the expression
 | (1) |
To obtain age-specific reproductive success (pi), we used the mean number of progeny per male (b
) and female (b
) for each particular age class as "proxy" estimators. These are the proportions of sexually mature rockfish at age from Nichol and Pikitch (1994) times the mean weight at age from (Rogers et al. 2000; Rogers 2003). These variables were adjusted to result in a net reproduction rate 
for each sex. Finally, we derived a value of reproductive success weighing both sexes equally: 
(Jorde and Ryman 1996).
An estimate of Ne () was obtained using the temporal method for overlapping generations as described in Jorde and Ryman (1995). Briefly, we used the equation:
 | (2) |
is the grand mean of temporal allele frequency change (among all adjacent cohorts and over all loci) derived from Pollak (1983) and corrected for sampling variance. Equations to compute C, and G are found in Jorde and Ryman (1995). In addition, the 95% confidence interval (CI) for and was calculated according to Waples (1989) assuming a chi-square distribution for F.
The estimated census population size (
) corresponded to the arithmetic mean of total numbers of fish in the population between 2001 and 2002 (Rogers et al. 2000; Rogers 2003).
|
Results and Discussion |
|---|
|
TopMaterials and MethodsResults and DiscussionReferences |
|---|
Cohort genetic statistics did not vary significantly among cohorts; in addition, no significant deviations from HWE were found and linkage equilibrium was met in all pairwise loci comparisons within each year class (Table 1). MICROCHECKER did not detect genotyping errors or null alleles. Lastly, the X2-statistic of Lewontin and Krakauer (1973) suggests that changes in microsatellite allele frequencies between cohorts are the sole effect of genetic drift (Table 2).
|
|
We observed a nearly exponential decline in survival (li) as a function of age; on the contrary, reproductive success (pi) peaked at 8–12 years (Figure 1). Both age-specific variables behaved according to theoretical expectations (e.g., Jorde and Ryman 1995). With these results, we then calculated the following values: C = 110.89 and G = 9.7. An estimate of
[95% CI] = 0.0006 [0.0005–0.0009] was obtained by averaging F' values across cohort pairs shown in Table 2. This provides an overall point estimate plus 95% CI of = 9157 [6495–12 215]. Based on = 24 376 500, the ratio could therefore vary between 0.001 and 0.0001 for darkblotched rockfish.
|
Wright (1931) predicted that under ideal conditions—even sex ratio, nonoverlapping generations, constant population size, panmixia, and Poisson variance in reproductive success among individuals—the ratio Ne/N should approach unity. However,
in natural populations average 0.1 using a variety of analytical methods (Frankham 1995), clearly indicating that these ideal conditions are commonly violated in real populations. What biological phenomena can explain the strikingly large discrepancy we observe between and values in darkblotched rockfish?
First, uneven sex ratios cannot greatly explain the small ratio found in S. crameri because extensive life-history data suggest no significant departures from 1:1 sex ratios (Nichol and Pikitch 1994; Rogers et al. 2000; Rogers 2003). Second, overlapping generations certainly violate ideal conditions. However, we have employed an ad hoc approach, with stable age structure and constant population size being the basic assumptions of the model. Ages are assumed to be errorless in the theoretical model; however, Rogers (2005) provides evidence that aging errors are commonplace in stock assessments of S. crameri, particularly for older individuals (>20 years old). Imprecise ages are nevertheless more likely to impact the accuracy of N than because 1) we used fish ranging from ages 1 to 7 to obtain the genetic information and 2) the correction factor (C) and mean generation length (G) depend primarily on survival and reproductive success, with both variables decaying rapidly for fish older than 12 years old (Figure 1).
Historical fishing has been demonstrated to truncate age structure and to diminish the population number of numerous groundfish stocks, including our species (Berkeley, Hixon, et al. 2004). From Jorde and Ryman (1996), it is clear that a reduction in age classes will decrease C and will therefore be underestimated; they further argue that monitoring the population over time can circumvent demographic disturbances affecting the correction factor C. Long-term genetic data collections during the current stage of stock recovery are thus necessary to validate C and in darkblotched rockfish.
Third, both significant spatial structure and genetic isolation by distance suggest that darkblotched rockfish is not best explained by a panmictic model (Gomez-Uchida and Banks 2005). In principle, a metapopulation with stepping-stone migration can be considered a large single panmictic population for purposes of estimating effective size if the overall genetic differentiation is small (Wang and Caballero 1999) as found in Gomez-Uchida and Banks (2005) for S. crameri. Although differences in size among subpopulations might potentially bias our downwardly (Wang and Caballero 1999), it is unlikely that this alone could account for the large difference in magnitude observed between and
Fourth, our estimates seem to agree closely with published and ratios for marine fishes with similar life history, including red drum (Turner et al. 2002), New Zealand snapper (Hauser et al. 2002), and Atlantic cod (Hutchinson et al. 2003). These fishes possess both high fecundity and mortality during early life stages (type III survivorship); this life strategy, when coupled with highly variable oceanographic conditions, could lead to a larger variance in reproductive success among individuals (Vk) than expected under a Poisson distribution, with the small minority of individuals () replacing the entire population () in every generation (Hedgecock 1994). Recent findings in black rockfish (Sebastes melanops) indicate that older females produce larvae of better quality than younger ones, hence more capable of producing competent offspring (Berkeley, Chapman, and Sogard 2004), supporting the hypothesis that Vk could vary greatly among individual female rockfish.
In conclusion, small for S. crameri is most likely explained by a combination of factors, such as high variance in reproductive success among individuals (Vk), genetic differences among subpopulations, and potential demographic disturbances such as intense fishing. Our results suggest that only several thousand darkblotched rockfish rather than millions in the whole population represent the breeding population number every generation. Given that small ratios are often found within populations of threatened species (Mace and Lande 1991), our estimates need careful consideration for management of vulnerable long-lived stocks (Musick 1999). However, further research is needed to 1) disentangle the individual effects of the aforementioned factors on Ne and 2) obtain independent estimates of key demographic parameters (e.g., Waples 2002) to validate ratios across marine species. Only then will we have confidence in extending findings made here for S. crameri to assess long-term conservation risks of further harvest of other rockfish and marine stocks.
|
Acknowledgments |
|---|
We thank J. B. Rogers at the Northwest Fisheries Science Center in Newport for her comprehensive assessments to understand the population dynamics of darkblotched rockfish and for sharing invaluable demographic information. We are indebted to L. Hauser, P. E. Jorde, R. S. Waples, D.E. Ruzzante, L.W. Seeb, and 3 anonymous referees for providing criticisms on several earlier versions of this paper. Funding for this investigation came from the Cooperative Institute of Marine Resources Studies. The primary author is grateful for a research assistantship from the Coastal Oregon Marine Experiment Station and a Mamie Markham Award from Hatfield Marine Science Center.
|
Footnotes |
|---|
Corresponding Editor: Lisa Seeb
Received December 21, 2005
Accepted October 3, 2006
|
References |
|---|
|
TopMaterials and MethodsResults and DiscussionReferences |
|---|
-
BelkhirK,BorsaP,ChikhiL,RaufasteN,BonhommeF. 2002. GENETIX 4.04, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions: CNRS UMR 5000, Université de Montpellier II, Montpellier, France [Internet]. [cited 2006 November 23] Available from: http://www.genetix.univ-montp2.fr/genetix/genetix.htm.
Berkeley SA, Chapman C, Sogard SM. (2004) Maternal age as a determinant of larval growth and survival in a marine fish, Sebastes melanops. Ecology 85:1258–1264.
Berkeley SA, Hixon MA, Larson RJ, Love MS. (2004) Fisheries sustainability via protection of age structure and spatial distribution of fish populations. Fisheries 29:23–32.
Chilton DE and Beamish RJ. (1982) Age determination methods for fishes studied by the Groundfish program at the Pacific Biological Station. Can Spec Publ Fish Aquat Sci 60:1–102.
Frankham R. (1995) Conservation genetics. Annu Rev Genet 29:305–327.[CrossRef][ISI][Medline]
Gomez-Uchida D and Banks MA. (2005) Microsatellite analyses of spatial genetic structure in darkblotched rockfish Sebastes crameri: is pooling samples safe? Can J Fish Aquat Sci 62:1874–1886.[CrossRef]
Hauser L, Adcock GJ, Smith PJ, Ramirez JHB, Carvalho GR. (2002) Loss of microsatellite diversity and low effective population size in an overexploited population of New Zealand snapper (Pagrus auratus). Proc Natl Acad Sci USA 99:11742–11747.
Hedgecock D. (1994) Does variance in reproductive success limit effective population sizes of marine organisms? In Beaumont AR (Ed.). Genetics and evolution of aquatic organisms(Chapman and Hall, London) pp. 122–134.
Hedrick PW. (2000) Genetics of populations. (Jones and Bartlett Publishers, Sudbury (MA)).
Hutchinson WF, van Oosterhaut C, Rogers SI, Carvalho GR. (2003) Temporal analysis of archived samples indicates marked genetic changes in declining North Sea cod (Gadus morhua). Proc R Soc Lond Ser B Biol Sci 270:2125–2132.[Medline]
Jorde PE and Ryman N. (1995) Temporal allele frequency change and estimation of effective size in populations with overlapping generations. Genetics 139:1077–1090.[Abstract]
Jorde PE and Ryman N. (1996) Demographic genetics of brown trout (Salmo trutta) and estimation of effective population size from temporal change of allele frequencies. Genetics 143:1369–1381.[Abstract]
Lande R. (1988) Genetics and demography in biological conservation. Science 241:1455–1460.
Lewontin RC and Krakauer J. (1973) Distribution of gene frequency as a test of the theory of the selective neutrality of polymorphisms. Genetics 74:175–195.
Love MS, Yoklavich MM, Thorsteinson L. (2002) The rockfishes of the northeast Pacific. (University of California Press, Los Angeles (CA)).
Mace GM and Lande R. (1991) Assessing extinction threats: toward a reevaluation of IUCN threatened species categories. Conserv Biol 5:148–157.[CrossRef]
Musick JA. (1999) Ecology and conservation of long-lived marine animals. Am Fish Soc Symp 23:1–10.
Nichol DG and Pikitch EK. (1994) Reproduction of darkblotched rockfish off the Oregon coast. Trans Am Fish Soc 123:469–481.[CrossRef]
Pollak E. (1983) A new method for estimating the effective population size from allele frequency changes. Genetics 104:531–548.
Raymond M and Rousset F. (1995) GENEPOP (version 1.2): population genetics software for exact tests and ecumenicism. J Hered 86:248–249.
RogersJB. 2003. Darkblotched rockfish (Sebastes crameri) 2003 stock status and rebuilding update [Internet]. [cited 2006 November 23]. Pacific Fishery Management Council, Portland, Oregon. Available from: http://www.pcouncil.org/groundfish/gfsafe0803/darkblotched.pdf.
RogersJB. 2005. Status of the darkblotched rockfish (Sebastes crameri) resource in 2005 [Internet]. [cited 2006 November 23]. Pacific Fishery Management Council 7700 NE Ambassador Place, Suite 101 Portland, Oregon 97220-1384. Available from: http://www.pcouncil.org/groundfish/gfsafe0406/darkblotched2005b_7-6-05_4SAFE.pdf.
RogersJB,MethotR,BuilderTL,PinerK,WilkinsM. 2000. Status of the darkblotched rockfish (Sebastes crameri) resource in 2000. Appendix: status of the Pacific coast groundfish fishery through 2000 and recommended acceptable biological catches for 2001 [Internet]. [cited 2006 November 23]. Pacific Fishery Management Council, Portland, Oregon. Available from: http://www.pcouncil.org/groundfish/gfsafe1000/darkblotched.pdf.
Turner TF, Wares JP, Gold JR. (2002) Genetic effective size is three orders of magnitude smaller than adult census size in an abundant, estuarine-dependent marine fish (Sciaenops ocellatus). Genetics 162:1329–1339.
van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. (2004) MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes 4:535–538.[CrossRef]
Wang JL and Caballero A. (1999) Developments in predicting the effective size of subdivided populations. Heredity 82:212–226.[CrossRef]
Waples RS. (1989) A generalized approach for estimating effective population size from temporal changes in allele frequency. Genetics 121:379–391.
Waples RS. (2002) Evaluating the effect of stage-specific survivorship on the Ne/N ratio. Mol Ecol 11:1029–1037.[CrossRef][Medline]
Wright S. (1931) Evolution in Mendelian populations. Genetics 16:114–138.
CiteULike 
Connotea 
Del.icio.us What's this?
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||