Coming out of the starting blocks: extended lag time rearranges genetic diversity in introduced marine fishes of Hawai'i

Biological invasions with known histories are rare, especially in the sea, and empirical studies of the genetic consequences are even rarer. Fifty-five years ago, the state of Hawai'i began a remarkable, if unintentional, 'experiment' with the introduction of three reef fishes, Lutjanus fulvus, Cephalopholis argus and Lutjanus kasmira. All have since expanded from the initial introduction of 2204 to 3163 individuals; however, historical records show that initially L. fulvus remained scarce, C. argus had modest population expansion and L. kasmira experienced rapid population growth. The consequences of differential population growth rates are apparent in F-statistics: Hawaiian L. fulvus demonstrate strong and significant haplotype frequency shifts from the founder location (F(ST) = 0.449), C. argus shows low but significant differentiation (F(ST) = 0.066) and L. kasmira is nearly identical to the founder location (F(ST) = 0.008). All three species had higher mtDNA diversity in the introduced range, which can be explained by multiple sources for L. fulvus and L. kasmira, but not for C. argus. We conclude that lag time before population expansion, in conjunction with genetic drift, has defined the genetic architecture of these three species in the introduced range.     

Effective population size associated with self-fertilization: lessons from temporal changes in allele frequencies in the selfing annual Medicago truncatula

Despite its significance in evolutionary and conservation biology, few estimates of effective population size (N(e)) are available in plant species. Self-fertilization is expected to affect N(e), through both its effect on homozygosity and population dynamics. Here, we estimated N(e) using temporal variation in allele frequencies for two contrasted populations of the selfing annual Medicago truncatula: a large and continuous population and a subdivided population. Estimated N(e) values were around 5-10% of the population census size suggesting that other factors than selfing must contribute to variation in allele frequencies. Further comparisons between monolocus allelic variation and changes in the multilocus genotypic composition of the populations show that the local dynamics of inbred lines can play an important role in the fluctuations of allele frequencies. Finally, comparing N(e) estimates and levels of genetic variation suggest that H(e) is a poor estimator of the contemporaneous variance effective population size.     

On the Potential for Estimating the Effective Number of Breeders from Heterozygote-Excess in Progeny

The important parameter of effective population size is rarely estimable directly from demographic data. Indirect estimates of effective population size may be made from genetic data such as temporal variation of allelic frequencies or linkage disequilibrium in cohorts. We suggest here that an indirect estimate of the effective number of breeders might be based on the excess of heterozygosity expected in a cohort of progeny produced by a limited number of males and females. In computer simulations, heterozygote excesses for 30 unlinked loci having various numbers of alleles and allele-frequency profiles were obtained for cohorts produced by samples of breeders drawn from an age-structured population and having known variance in reproductive success and effective number. The 95% confidence limits around the estimate contained the true effective population size in 70 of 72 trials and the Spearman rank correlation of estimated and actual values was 0.991. An estimate based on heterozygote excess might have certain advantages over the previous estimates, requiring only single-locus and single-cohort data, but the sampling error among individuals and the effect of departures from random union of gametes still need to be explored.     

Molecular evidence for historical and recent population size reductions of tiger salamanders (Ambystoma tigrinum) in Yellowstone National Park

Population declines caused by natural and anthropogenic factors can quickly erode genetic diversity in natural populations. In this study, we examined genetic variation within 10 tiger salamander populations across northern Yellowstone National Park in Wyoming and Montana, USA using eight microsatellite loci. We tested for the genetic signature of population decline using heterozygosity excess, shifts in allele frequencies, and low ratios of allelic number to allelic size range (M-ratios). We found different results among the three tests. All 10 populations had low M-ratios, five had shifts in allele frequencies and only two had significant heterozygosity excesses. These results support theoretical expectations of different temporal signatures among bottleneck tests and suggest that both historical fish stocking, recent, sustained drought, and possibly an emerging amphibian disease have contributed to declines in effective population size.     

Temporal estimates of effective population size in species with overlapping generations

The standard temporal method for estimating effective population size (N(e)) assumes that generations are discrete, but it is routinely applied to species with overlapping generations. We evaluated bias in the estimates N(e) caused by violation of this assumption, using simulated data for three model species: humans (type I survival), sparrow (type II), and barnacle (type III). We verify a previous proposal by Felsenstein that weighting individuals by reproductive value is the correct way to calculate parametric population allele frequencies, in which case the rate of change in age-structured populations conforms to that predicted by discrete-generation models. When the standard temporal method is applied to age-structured species, typical sampling regimes (sampling only newborns or adults; randomly sampling the entire population) do not yield properly weighted allele frequencies and result in biased N(e). The direction and magnitude of the bias are shown to depend on the sampling method and the species' life history. Results for populations that grow (or decline) at a constant rate paralleled those for populations of constant size. If sufficient demographic data are available and certain sampling restrictions are met, the Jorde-Ryman modification of the temporal method can be applied to any species with overlapping generations. Alternatively, spacing the temporal samples many generations apart maximizes the drift signal compared to sampling biases associated with age structure.     

High genetic diversity despite drastic bottleneck in a critically endangered,long‐lived seabird,the Mascarene Petrel Pseudobulweria aterrima

The Mascarene Petrel Pseudobulweria aterrima is a critically endangered seabird endemic to Reunion Island, with an extremely small population suffering several threats. Fifteen polymorphic microsatellite loci were isolated from this species to analyse genetic diversity, estimate contemporary effective population size, search for evidence of a population bottleneck and see whether results support the hypothesis that life history traits could preserve allelic diversity in small populations. Results from 22 individuals found grounded as a consequence of light pollution highlight a surprisingly high genetic diversity, an absence of inbreeding, a contemporary effective population size estimated at approximately 1211 individuals and a probable bottleneck around 10 000 generations ago. Additional studies on genetic diversity and structure from a larger number of samples are thus required to evaluate the evolutionary potential of this critically endangered species.     

How Informative Is Wright's Estimator of the Number of Genes Affecting a Quantitative Character?

S. Wright suggested an estimator, m, of the number of loci, m, contributing to the difference in a quantitative character between two differentiated populations, which is calculated from the phenotypic means and variances in the two parental populations and their F1 and F2 hybrids. The same method can also be used to estimate m contributing to the genetic variance within a single population, by using divergent selection to create differentiated lines from the base population. In this paper we systematically examine the utility and problems of this technique under the influences of unequal allelic effects and initial allele frequencies, and linkage, which are known to lead m to underestimate m. In addition, we examine the effects of population size and selection intensity during the generations of selection. During selection, the estimator m rapidly approaches its expected value at the selection limit. With reasonable assumptions about unequal allelic effects and initial allele frequencies, the expected value of m without linkage is likely to be on the order of one-third of the number of genes. The estimates suffer most seriously from linkage. The practical maximum expectation of m is just about the number of chromosomes, considerably less than the "recombination index" which has been assumed to be the upper limit. The estimates are also associated with large sampling variances. An estimator of the variance of m derived by R. Lande substantially underestimates the actual variance. Modifications to the method can ameliorate some of the problems. These include using F3 or later generation variances or the genetic variance in the base population, and replicating the experiments and estimation procedure. However, even in the best of circumstances, information from m is very limited and can be misleading.     

Where are they now? The fates of two genetic lineages in an introduced Hawaiian reef fish

Humans, both wittingly and unwittingly, have been transporting marine organisms beyond their native ranges for centuries ( Ruiz et al. 1997 ). A central challenge of invasion biology is to identify the factors that determine whether introduced species fail to become established, become benign members of a community, or spread so far and reach such densities as to be considered invasive. Organismal features such as physiological tolerance, niche breadth and fecundity are critical, but by themselves are inaccurate predictors of the fates of introduced species ( Sakai et al. 2001 ). The size, age distribution, and genetic makeup of founder populations are also important, but because they are usually unknown they are most often viewed as sources of uncertainty. For marine species with planktonic larvae, the challenge is even greater because the consequences of a planktonic phase for dispersal and population viability are not well understood. In this issue, Gaither et al. (2010a) present a remarkable account of the introduction of a reef fish for which the number and genetic makeup of the founders are known. Between 1956 and 1961, the Division of Fish and Game for the Territory of Hawaii introduced 12 non‐indigenous fish species into Hawaiian waters to establish commercial and sport fisheries. The introduction of Lutjanus kasmira, the bluestriped snapper, was the most successful ( Fig. 1 ). There were two releases of fish from French Polynesia. In 1958, 2431 fish from the Marquesas Islands were released on Oahu, followed in 1961 with an additional 728 fish from the Society Islands. The blue striped snapper rapidly spread to the other Hawaiian Islands, reaching the northwestern end of the archipelago by 1992. The choice of the Marquesas as one of two sources for the introduction was fortuitous. Gaither et al. (2010b) found that the Marquesas population is genetically distinct from all other Indo‐Pacific populations of L. kasmira. Mitochondrial cytochrome b sequences of fish from the Marquesas belong to a separate lineage that diverged from others in the species roughly half a Ma. Allele frequencies for several nuclear loci are also distinct. This provided Gaither et al. (2010a) with an extraordinary opportunity to examine what became of the mixed genetic heritage of Hawaiian blue striped snappers after 50 years.

Figure 1 Open in figure viewer PowerPoint The bluestriped snapper, Lutjanus kasmira, introduced to Hawaii 50 years ago and now an abundant reef fish expanded from a small founder population with minimal changes in the diversity or frequencies of mitochondrial and nuclear genetic markers.     

Ecology of the introduced snapper Lutjanus kasmiva (Forsskal) in the reef fish assemblage of a Hawaiian bay

Since its first introduction in 1955, the snapper Lutjanus kasmira (taape) has developed large populations in shallow coastal waters of Hawaii. Visual abundance estimates and trap catches of taape were different between the main habitat types of the reef. Overall, taape was the second most abundant species by numbers and biomass over a hard substratum and was patchily distributed, even within habitat type. The largest individuals occurred singly or in small groups on shallow reef slopes; on the deep slopes and in spur-and-groove habitat, taape occurred in larger groups of smaller sized fish. Catches in fish traps showed a negative relationship with the complexity of the adjacent natural substratum, suggesting that traps are more productive where they offer additional cover in areas with less natural cover. In all areas in which traps were used, from fringing reef to open sand, the catch was dominated by taape. Tagging confirmed the theory that taape tend to maintain a limited range over long periods. Tagging results also indicated occasional long-range movements among major habitats and provided data confirming patterns of movement between reef habitat by day and sand habitat by night. Taape collected over open sand substratum by trap and by hook-and-line were smaller than those collected by the same means at the reef. Observations and collections suggested an ontogenetic trend in habitat use by taape from initial settlement in fringe areas toward ultimate residence on the main reef. Other common fish species in the reef assemblages were tested for statistical associations with taape. Although several positive species associations were found, multivariate analysis did not reveal patterns that indicated strong ecological relationships.     

Age and growth of Lutjanus kasmira (Forskål) in Hawaiian waters

The growth of Hawaiian taape, Lutjanus kasmira , was studied by examining otoliths and by analysing length-frequency distribution. Annual hyaline and opaque markings were visible in whole mounts of sagittae, which were verified by enumeration of daily increments with a scanning electron microscope (SEM) and through marginal increment analysis. The von Bertalanffy growth curve was fitted to the data, resulting in: where t. l . is total length (cm) and t is age (years). SEM observations revealed that the slowgrowth hyaline zones were composed of daily increments too small (0.4–0.8 μm) to be resolved optically. Thus, age estimates derived by numerically integrating otolith growth rate data obtained with a light microscope showed a negative bias, resulting in overestimation of growth rates. Parameter estimates obtained from three different types of length-frequency analysis were also unstable. This was due, at least in part, to differences in the size composition of fish sampled with different fishing gears and from different depths.
The growth rate registered in Hawaii falls within the reported growth coefficients of lutjanids, whereas it is one of the highest in the Pacific and clearly higher than a deep-water lutjanid species growth in Hawaii. Probably, this high growth rate may have been enhanced by the relative lack of competitors in the depauperate Hawaiian marine fish community.     

Temporal Variation in Genetic Structure of the Chinese Rare Minnow (Gobiocypris rarus) in Its Type Locality Revealed by Microsatellite Markers

Gobiocypris rarus, an endemic cyprinid fish with high fecundity, lives mainly in small water systems easily influenced by changes in natural surroundings. This study used 11 polymorphic microsatellite primers to identify the temporal variation of its topotype population. Moderate genetic diversity, inbreeding phenomena, and limited temporal variation between 1997 and 2006 were revealed in the topotype population. The main temporal fluctuations involved only the change of allelic frequencies over two loci and allelic richness. The effective population size was estimated to be 645. The authors argue that inbreeding did not induce dramatic depression effects on the topotype population, and the forces to maintain genetic diversity were mainly from environmental fluctuations and life history traits. Considering that the topotype population is facing increased habitat loss, destruction, and disturbance due to human activities, the authors suggest that a habitat and species management area be established in the type locality.     

Using maximum likelihood to estimate population size from temporal changes in allele frequencies.

We develop a maximum-likelihood framework for using temporal changes in allele frequencies to estimate the number of breeding individuals in a population. We use simulations to compare the performance of this estimator to an F-statistic estimator of variance effective population size. The maximum-likelihood estimator had a lower variance and smaller bias. Taking advantage of the likelihood framework, we extend the model to include exponential growth and show that temporal allele frequency data from three or more sampling events can be used to test for population growth.     

Experimental evaluation of the usefulness of microsatellite DNA for detecting demographic bottlenecks

Evolutionary and conservation biologists often use molecular markers to evaluate whether populations have experienced demographic bottlenecks that resulted in a loss of genetic variation. We evaluated the utility of microsatellites for detection of recent, severe bottlenecks and compared the amounts of genetic diversity lost in bottlenecks of different sizes. In experimental mesocosms, we established replicate populations by releasing 1, 2, 4 or 8 pairs of the western mosquitofish, Gambusia affinis (Poeciliidae). Using eight polymorphic microsatellite loci, we quantified seven indices of genetic diversity or change that have been used to assess the effects of demographic bottlenecks on populations. We compared indices for the experimentally bottlenecked populations to those for the source population and examined differences between populations established with different numbers of founders. Direct count heterozygosity and the proportion of polymorphic loci were not very sensitive to genetic changes that resulted from the experimental bottlenecks. Heterozygosity excess and expected heterozygosity were useful to varying degrees in the detection of bottlenecks. Allelic diversity and temporal variance in allele frequencies were most sensitive to genetic changes that resulted from the bottlenecks, and the temporal variance method was slightly more correlated with bottleneck size than was allelic diversity. Based on comparisons to a previous study with allozymes, heterozygosity, temporal variance in allele frequencies and allelic diversity, but not proportion of polymorphic loci, appear to be more sensitive to demographic bottlenecks when quantified using microsatellites. We found that analysis of eight highly polymorphic loci was sufficient to detect a recent demographic bottleneck and to obtain an estimate of the magnitude of bottleneck severity.     

Samples from subdivided populations yield biased estimates of effective size that overestimate the rate of loss of genetic variation

Many empirical studies estimating effective population size apply the temporal method that provides an estimate of the variance effective size through the amount of temporal allele frequency change under the assumption that the study population is completely isolated. This assumption is frequently violated, and the magnitude of the resulting bias is generally unknown. We studied how gene flow affects estimates of effective size obtained by the temporal method when sampling from a population system and provide analytical expressions for the expected estimate under an island model of migration. We show that the temporal method tends to systematically underestimate both local and global effective size when populations are connected by gene flow, and the bias is sometimes dramatic. The problem is particularly likely to occur when sampling from a subdivided population where high levels of gene flow obscure identification of subpopulation boundaries. In such situations, sampling in a manner that prevents biased estimates can be difficult. This phenomenon might partially explain the frequently reported unexpectedly low effective population sizes of marine populations that have raised concern regarding the genetic vulnerability of even exceptionally large populations.     

Local and Temporal Variation at the MHC I-<Emphasis Type="Italic">A1</Emphasis> Locus in Pink Salmon <Emphasis Type="Italic">Oncorhynchus gorbuscha</Emphasis> of Southern Sakhalin

In order to test the MHC I-A1 locus of the main histocompatiblity complex as an effective molecular marker for identification of populations of Far Eastern pink salmon, Oncorhynchus gorbuscha, the comparison of the allelic and genotypic frequencies between two neighboring populations of the Lyutoga and Taranay rivers (Aniva Bay, Sakhalin Island), as well as within the populations: within generation, between adjacent and distant generations, have been performed. It was revealed that temporal variation at MHC I-A1 is statistically insignificant in comparison between the samples of the juveniles sampled at different times during the downstream migration. It is also insignificant between the samples of juvenile and adult fish of the same generation, and between samples divided by a single and seven generations. At the same time, the samples of neighboring populations differ significantly. The possibility of differentiation of the populations on the local scale and temporal stability of the allele frequencies indicate the MHC I-A1 as valuable marker for the development of the population data bases. Inclusion of the MHC I-A1 into the panel of other markers (microsatellites, in particular) allows for considerable increase in the resolution power for identification as to which population an individual in the mixed sea catch belongs.     

Historical Analysis of Genetic Variation Reveals Low Effective Population Size in a Northern Pike (Esox Lucius) Population

Effective population size (N(e)) of a natural fish population was estimated from temporal changes in allele frequencies at seven microsatellite loci. Use of a historical collection of fish scales made it possible to increase the precision of estimates by increasing the time interval between samples and to use an equation developed for discrete generations without correcting for demographic parameters. Estimates of N(e) for the time intervals 1961-1977 and 1977-1993 were 35 and 72, respectively. For the entire interval, 1961-1993, the estimate of N(e) was 48 when based on a weighted mean derived from the above two estimates or 125 when calculated from 1961 and 1993 samples only. Corresponding ratios of effective size to adult census size ranged from 0.03 to 0.14. An N(e) of 48 over a 32-year period would imply that this population lost as much as 8% of its heterozygosity in that time. Results suggest the potential for using genetic methods based on microsatellite loci data to compare historical trends in N(e) with population dynamic parameters. Such comparisons will help to evaluate the relationship between genetic diversity and long-term persistence of natural populations.     

A note on effective population size with overlapping generations

A simple derivation is given for a formula obtained previously for the effective size of random-mating populations with overlapping generations. The effective population size is the same as that for a population with discrete generations having the same variance of lifetime family size and the same number of individuals entering the population per generation.     

Gene Genealogy and Variance of Interpopulational Nucleotide Differences

A mathematical theory is developed for computing the probability that m genes sampled from one population (species) and n genes sampled from another are derived from l genes that existed at the time of population splitting. The expected time of divergence between the two most closely related genes sampled from two different populations and the time of divergence (coalescence) of all genes sampled are studied by using this theory. It is shown that the time of divergence between the two most closely related genes can be used as an approximate estimate of the time of population splitting (T) only when T identical to t/(2N) is small, where t and N are the number of generations and the effective population size, respectively. The variance of Nei and Li's estimate (d) of the number of net nucleotide differences between two populations is also studied. It is shown that the standard error (Sd) of d is larger than the mean when T is small (T much less than 1). In this case, Sd is reduced considerably by increasing sample size. When T is large (T greater than 1), however, a large proportion of the variance of d is caused by stochastic factors, and increase in the sample size does not help to reduce Sd. To reduce the stochastic variance of d, one must use data from many independent unlinked gene loci.     

Effective size of a wild salmonid population is greatly reduced by hatchery supplementation

Many declining and commercially important populations are supplemented with captive-born individuals that are intentionally released into the wild. These supplementation programs often create large numbers of offspring from relatively few breeding adults, which can have substantial population-level effects. We examined the genetic effects of supplementation on a wild population of steelhead (Oncorhynchus mykiss) from the Hood River, Oregon, by matching 12 run-years of hatchery steelhead back to their broodstock parents. We show that the effective number of breeders producing the hatchery fish (broodstock parents; N(b)) was quite small (harmonic mean N(b)=25 fish per brood-year vs 373 for wild fish), and was exacerbated by a high variance in broodstock reproductive success among individuals within years. The low N(b) caused hatchery fish to have decreased allelic richness, increased average relatedness, more loci in linkage disequilibrium and substantial levels of genetic drift in comparison with their wild-born counterparts. We also documented a substantial Ryman-Laikre effect whereby the additional hatchery fish doubled the total number of adult fish on the spawning grounds each year, but cut the effective population size of the total population (wild and hatchery fish combined) by nearly two-thirds. We further demonstrate that the Ryman-Laikre effect is most severe in this population when (1) >10% of fish allowed onto spawning grounds are from hatcheries and (2) the hatchery fish have high reproductive success in the wild. These results emphasize the trade-offs that arise when supplementation programs attempt to balance disparate goals (increasing production while maintaining genetic diversity and fitness).     

The genetic effective and adult census size of an Australian population of tiger prawns (Penaeus esculentus)

This study compares estimates of the census size of the spawning population with genetic estimates of effective current and long-term population size for an abundant and commercially important marine invertebrate, the brown tiger prawn (Penaeus esculentus). Our aim was to focus on the relationship between genetic effective and census size that may provide a source of information for viability analyses of naturally occurring populations. Samples were taken in 2001, 2002 and 2003 from a population on the east coast of Australia and temporal allelic variation was measured at eight polymorphic microsatellite loci. Moments-based and maximum-likelihood estimates of current genetic effective population size ranged from 797 to 1304. The mean long-term genetic effective population size was 9968. Although small for a large population, the effective population size estimates were above the threshold where genetic diversity is lost at neutral alleles through drift or inbreeding. Simulation studies correctly predicted that under these experimental conditions the genetic estimates would have non-infinite upper confidence limits and revealed they might be overestimates of the true size. We also show that estimates of mortality and variance in family size may be derived from data on average fecundity, current genetic effective and census spawning population size, assuming effective population size is equivalent to the number of breeders. This work confirms that it is feasible to obtain accurate estimates of current genetic effective population size for abundant Type III species using existing genetic marker technology.