Mark R. Christie

Genetic adaptation to captivity can occur in a single generation

Captive breeding programs are widely used for the conservation and restoration of threatened and endangered species. Nevertheless, captive-born individuals frequently have reduced fitness when reintroduced into the wild. The mechanism for these fitness declines has remained elusive, but hypotheses include environmental effects of captive rearing, inbreeding among close relatives, relaxed natural selection, and unintentional domestication selection (adaptation to captivity). We used a multigenerational pedigree analysis to demonstrate that domestication selection can explain the precipitous decline in fitness observed in hatchery steelhead released into the Hood River in Oregon. After returning from the ocean, wild-born and first-generation hatchery fish were used as broodstock in the hatchery, and their offspring were released into the wild as smolts. First-generation hatchery fish had nearly double the lifetime reproductive success (measured as the number of returning adult offspring) when spawned in captivity compared with wild fish spawned under identical conditions, which is a clear demonstration of adaptation to captivity. We also documented a tradeoff among the wild-born broodstock: Those with the greatest fitness in a captive environment produced offspring that performed the worst in the wild. Specifically, captive-born individuals with five (the median) or more returning siblings (i.e., offspring of successful broodstock) averaged 0.62 returning offspring in the wild, whereas captive-born individuals with less than five siblings averaged 2.05 returning offspring in the wild. These results demonstrate that a single generation in captivity can result in a substantial response to selection on traits that are beneficial in captivity but severely maladaptive in the wild.

Self‐recruitment and sweepstakes reproduction amid extensive gene flow in a coral‐reef fish

Identifying patterns of larval dispersal within marine metapopulations is vital for effective fisheries management, appropriate marine reserve design, and conservation efforts. We employed genetic markers (microsatellites) to determine dispersal patterns in bicolour damselfish (Pomacentridae: Stegastes partitus). Tissue samples of 751 fish were collected in 2004 and 2005 from 11 sites encompassing the Exuma Sound, Bahamas. Bayesian parentage analysis identified two parent–offspring pairs, which is remarkable given the large population sizes and 28 day pelagic larval duration of bicolour damselfish. The two parent–offspring pairs directly documented self‐recruitment at the two northern‐most sites, one of which is a long‐established marine reserve. Principal coordinates analyses of pair‐wise relatedness values further indicated that self‐recruitment was common in all sampled populations. Nevertheless, measures of genetic differentiation (FST) and results from assignment methods suggested high levels of gene flow among populations. Comparisons of heterozygosity and relatedness among samples of adults and recruits indicated spatially and temporally independent sweepstakes events, whereby only a subset of adults successfully contribute to subsequent generations. These results indicate that self‐recruitment and sweepstakes reproduction are the predominant, ecologically‐relevant processes that shape patterns of larval dispersal in this system.

Who are the missing parents? Grandparentage analysis identifies multiple sources of gene flow into a wild population

In order to increase the size of declining salmonid populations, supplementation programmes intentionally release fish raised in hatcheries into the wild. Because hatchery‐born fish often have lower fitness than wild‐born fish, estimating rates of gene flow from hatcheries into wild populations is essential for predicting the fitness cost to wild populations. Steelhead trout (Oncorhynchus mykiss) have both freshwater resident and anadromous (ocean‐going) life history forms, known as rainbow trout and steelhead, respectively. Juvenile hatchery steelhead that ‘residualize’ (become residents rather than go to sea as intended) provide a previously unmeasured route for gene flow from hatchery into wild populations. We apply a combination of parentage and grandparentage methods to a three‐generation pedigree of steelhead from the Hood River, Oregon, to identify the missing parents of anadromous fish. For fish with only one anadromous parent, 83% were identified as having a resident father while 17% were identified as having a resident mother. Additionally, we documented that resident hatchery males produced more offspring with wild anadromous females than with hatchery anadromous females. One explanation is the high fitness cost associated with matings between two hatchery fish. After accounting for all of the possible matings involving steelhead, we find that only 1% of steelhead genes come from residualized hatchery fish, while 20% of steelhead genes come from wild residents. A further 23% of anadromous steelhead genes come from matings between two resident parents. If these matings mirror the proportion of matings between residualized hatchery fish and anadromous partners, then closer to 40% of all steelhead genes come from wild trout each generation. These results suggest that wild resident fish contribute substantially to endangered steelhead ‘populations’ and highlight the need for conservation and management efforts to fully account for interconnected Oncorhynchus mykiss life histories.

Parentage in natural populations: novel methods to detect parent-offspring pairs in large data sets

Parentage analysis in natural populations presents a valuable yet unique challenge because of large numbers of pairwise comparisons, marker set limitations and few sampled true parent–offspring pairs. These limitations can result in the incorrect assignment of false parent–offspring pairs that share alleles across multi‐locus genotypes by chance alone. I first define a probability, Pr(δ), to estimate the expected number of false parent–offspring pairs within a data set. This probability can be used to determine whether one can accept all putative parent–offspring pairs with strict exclusion. I next define the probability Pr(φ|λ), which employs Bayes’ theorem to determine the probability of a putative parent–offspring pair being false given the frequencies of shared alleles. This probability can be used to separate true parent–offspring pairs from false pairs that occur by chance when a data set lacks sufficient numbers of loci to accept all putative parent–offspring pairs. Finally, I propose a method to quantitatively determine how many loci to let mismatch for study‐specific error rates and demonstrate that few data sets should need to allow more than two loci to mismatch. I test all theoretical predictions with simulated data and find that, first, Pr(δ) and Pr(φ|λ) have very low bias, and second, that power increases with lower sample sizes, uniform allele frequency distributions, and higher numbers of loci and alleles per locus. Comparisons of Pr(φ|λ) to strict exclusion and CERVUS demonstrate that this method may be most appropriate for large natural populations when supplemental data (e.g. genealogies, candidate parents) are absent.

On the reproductive success of early-generation hatchery fish in the wild

Large numbers of hatchery salmon spawn in wild populations each year. Hatchery fish with multiple generations of hatchery ancestry often have heritably lower reproductive success than wild fish and may reduce the fitness of an entire population. Whether this reduced fitness also occurs for hatchery fish created with local‐ and predominantly wild‐origin parents remains controversial. Here, we review recent studies on the reproductive success of such ‘early‐generation’ hatchery fish that spawn in the wild. Combining 51 estimates from six studies on four salmon species, we found that (i) early‐generation hatchery fish averaged only half the reproductive success of their wild‐origin counterparts when spawning in the wild, (ii) the reduction in reproductive success was more severe for males than for females, and (iii) all species showed reduced fitness due to hatchery rearing. We review commonalities among studies that point to possible mechanisms (e.g., environmental versus genetic effects). Furthermore, we illustrate that sample sizes typical of these studies result in low statistical power to detect fitness differences unless the differences are substantial. This review demonstrates that reduced fitness of early‐generation hatchery fish may be a general phenomenon. Future research should focus on determining the causes of those fitness reductions and whether they lead to long‐term reductions in the fitness of wild populations.

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; Nb) was quite small (harmonic mean Nb=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 Nb 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).

A single generation of domestication heritably alters the expression of hundreds of genes

The genetic underpinnings associated with the earliest stages of plant and animal domestication have remained elusive. Because a genome-wide response to selection can take many generations, the earliest detectable changes associated with domestication may first manifest as heritable changes to global patterns of gene expression. Here, to test this hypothesis, we measured differential gene expression in the offspring of wild and first-generation hatchery steelhead trout (Oncorhynchus mykiss) reared in a common environment. Remarkably, we find that there were 723 genes differentially expressed between the two groups of offspring. Reciprocal crosses reveal that the differentially expressed genes could not be explained by maternal effects or by chance differences in the background levels of gene expression among unrelated families. Gene-enrichment analyses reveal that adaptation to the novel hatchery environment involved responses in wound healing, immunity and metabolism. These findings suggest that the earliest stages of domestication may involve adaptation to highly crowded conditions.

Bayesian parentage analysis with systematic accountability of genotyping error, missing data and false matching

MOTIVATION The goal of any parentage analysis is to identify as many parent-offspring relationships as possible, while minimizing incorrect assignments. Existing methods can achieve these ends, but they require additional information in the form of demographic data, thousands of markers and/or estimates of genotyping error rates. For many non-model systems, it is simply not practical, cost-effective or logistically feasible to obtain this information. Here, we develop a Bayesian parentage method that only requires the sampled genotypes to account for genotyping error, missing data and false matches. RESULTS Extensive testing with microsatellite and SNP datasets reveals that our Bayesian parentage method reliably controls for the number of false assignments, irrespective of the genotyping error rate. When the number of loci is limiting, our approach maximizes the number of correct assignments by accounting for the frequencies of shared alleles. Comparisons with exclusion and likelihood-based methods on an empirical salmon dataset revealed that our Bayesian method had the highest ratio of correct to incorrect assignments.

QUANTIFYING EVOLUTIONARY POTENTIAL OF MARINE FISH LARVAE: HERITABILITY, SELECTION, AND EVOLUTIONARY CONSTRAINTS

For many marine fish, intense larval mortality may provide considerable opportunity for selection, yet much less is known about the evolutionary potential of larval traits. We combined field demographic studies and manipulative experiments to estimate quantitative genetic parameters for both larval size and swimming performance for a natural population of a common coral‐reef fish, the bicolor damselfish (Stegastes partitus). We also examined selection on larval size by synthesizing information from published estimates of selective mortality. We introduce a method that uses the Lande–Arnold framework for examining selection on quantitative traits to empirically reconstruct adaptive landscapes. This method allows the relationship between phenotypic value and fitness components to be described across a broad range of trait values. Our results suggested that despite strong viability selection for large larvae and moderate heritability (h2= 0.29), evolutionary responses of larvae would likely be balanced by reproductive selection favoring mothers that produce more, smaller offspring. Although long‐term evolutionary responses of larval traits may be constrained by size‐number trade‐offs, our results suggest that phenotypic variation in larval size may be an ecologically important source of variability in population dynamics through effects on larval survival and recruitment to benthic populations.

Spatial and temporal patterns of larval dispersal in a coral-reef fish metapopulation: evidence of variable reproductive success

Many marine organisms can be transported hundreds of kilometres during their pelagic larval stage, yet little is known about spatial and temporal patterns of larval dispersal. Although traditional population‐genetic tools can be applied to infer movement of larvae on an evolutionary timescale, large effective population sizes and high rates of gene flow present serious challenges to documenting dispersal patterns over shorter, ecologically relevant, timescales. Here, we address these challenges by combining direct parentage analysis and indirect genetic analyses over a 4‐year period to document spatial and temporal patterns of larval dispersal in a common coral‐reef fish: the bicolour damselfish (Stegastes partitus). At four island locations surrounding Exuma Sound, Bahamas, including a long‐established marine reserve, we collected 3278 individuals and genotyped them at 10 microsatellite loci. Using Bayesian parentage analysis, we identified eight parent‐offspring pairs, thereby directly documenting dispersal distances ranging from 0 km (i.e., self‐recruitment) to 129 km (i.e., larval connectivity). Despite documenting substantial dispersal and gene flow between islands, we observed more self‐recruitment events than expected if the larvae were drawn from a common, well‐mixed pool (i.e., a completely open population). Additionally, we detected both spatial and temporal variation in signatures of sweepstakes and Wahlund effects. The high variance in reproductive success (i.e., ‘sweepstakes’) we observed may be influenced by seasonal mesoscale gyres present in the Exuma Sound, which play a prominent role in shaping local oceanographic patterns. This study documents the complex nature of larval dispersal in a coral‐reef fish, and highlights the importance of sampling multiple cohorts and coupling both direct and indirect genetic methods in order disentangle patterns of dispersal, gene flow and variable reproductive success.