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

Whole-Genome Resequencing to Evaluate Life History Variation in Anadromous Migration of Oncorhynchus mykiss

Anadromous fish experience physiological modifications necessary to migrate between vastly different freshwater and marine environments, but some species such as Oncorhynchus mykiss demonstrate variation in life history strategies with some individuals remaining exclusively resident in freshwater, whereas others undergo anadromous migration. Because there is limited understanding of genes involved in this life history variation across populations of this species, we evaluated the genomic difference between known anadromous (n = 39) and resident (n = 78) Oncorhynchus mykiss collected from the Klickitat River, WA, USA, with whole-genome resequencing methods. Sequencing of these collections yielded 5.64 million single-nucleotide polymorphisms that were tested for significant differences between resident and anadromous groups along with previously identified candidate gene regions. Although a few regions of the genome were marginally significant, there was one region on chromosome Omy12 that provided the most consistent signal of association with anadromy near two annotated genes in the reference assembly: COP9 signalosome complex subunit 6 (CSN6) and NACHT, LRR, and PYD domain–containing protein 3 (NLRP3). Previously identified candidate genes for anadromy within the inversion region of chromosome Omy05 in coastal steelhead and rainbow trout were not informative for this population as shown in previous studies. Results indicate that the significant region on chromosome Omy12 may represent a minor effect gene for male anadromy and suggests that this life history variation in Oncorhynchus mykiss is more strongly driven by other mechanisms related to environmental rearing such as epigenetic modification, gene expression, and phenotypic plasticity. Further studies into regulatory mechanisms of this trait are needed to understand drivers of anadromy in populations of this protected species.

Endoplasmic reticulum stress as a key mechanism in stunted growth of seawater rainbow trout (Oncorhynchus mykiss)

BACKGROUND

Rainbow trout (Oncorhynchus mykiss) is a salmonid species with a complex life-history. Wild populations are naturally divided into freshwater residents and sea-run migrants. Migrants undergo an energy-demanding adaptation for life in seawater, known as smoltification, while freshwater residents display these changes in an attenuated magnitude and rate. Despite this, in seawater rainbow trout farming all fish are transferred to seawater. Under these circumstances, weeks after seawater transfer, a significant portion of the fish die (around 10%) or experience growth stunting (GS; around 10%), which represents an important profitability and welfare issue. The underlying causes leading to GS in seawater-transferred rainbow trout remain unknown. In this study, we aimed at characterising the GS phenotype in seawater-transferred rainbow trout using untargeted and targeted approaches. To this end, the liver proteome (LC-MS/MS) and lipidome (LC-MS) of GS and fast-growing phenotypes were profiled to identify molecules and processes that are characteristic of the GS phenotype. Moreover, the transcription, abundance or activity of key proteins and hormones related to osmoregulation (Gill Na+, K + -ATPase activity), growth (plasma IGF-I, and liver igf1, igfbp1b, ghr1 and ctsl) and stress (plasma cortisol) were measured using targeted approaches.

RESULTS

No differences in Gill Na+, K + -ATPase activity and plasma cortisol were detected between the two groups. However, a significant downregulation in plasma IGF-I and liver igf1 transcription pointed at this growth factor as an important pathomechanism for GS. Changes in the liver proteome revealed reactive-oxygen-species-mediated endoplasmic reticulum stress as a key mechanism underlying the GS phenotype. From the lipidomic analysis, key observations include a reduction in triacylglycerols and elevated amounts of cardiolipins, a characteristic lipid class associated with oxidative stress, in GS phenotype.

CONCLUSION

While the triggers to the activation of endoplasmic reticulum stress are still unknown, data from this study point towards a nutritional deficiency as an underlying driver of this phenotype.

An evaluation of the potential factors affecting lifetime reproductive success in salmonids

Lifetime reproductive success (LRS), the number of offspring produced over an organism's lifetime, is a fundamental component of Darwinian fitness. For taxa such as salmonids with multiple species of conservation concern, understanding the factors affecting LRS is critical for the development and implementation of successful conservation management practices. Here, we reviewed the published literature to synthesize factors affecting LRS in salmonids including significant effects of hatchery rearing, life history, and phenotypic variation, and behavioral and spawning interactions. Additionally, we found that LRS is affected by competitive behavior on the spawning grounds, genetic compatibility, local adaptation, and hybridization. Our review of existing literature revealed limitations of LRS studies, and we emphasize the following areas that warrant further attention in future research: (1) expanding the range of studies assessing LRS across different life-history strategies, specifically accounting for distinct reproductive and migratory phenotypes; (2) broadening the variety of species represented in salmonid fitness studies; (3) constructing multigenerational pedigrees to track long-term fitness effects; (4) conducting LRS studies that investigate the effects of aquatic stressors, such as anthropogenic effects, pathogens, environmental factors in both freshwater and marine environments, and assessing overall body condition, and (5) utilizing appropriate statistical approaches to determine the factors that explain the greatest variation in fitness and providing information regarding biological significance, power limitations, and potential sources of error in salmonid parentage studies. Overall, this review emphasizes that studies of LRS have profoundly advanced scientific understanding of salmonid fitness, but substantial challenges need to be overcome to assist with long-term recovery of these keystone species in aquatic ecosystems.

The future of parentage analysis: From microsatellites to SNPs and beyond

Parentage analysis is a cornerstone of molecular ecology that has delivered fundamental insights into behaviour, ecology and evolution. Microsatellite markers have long been the king of parentage, their hypervariable nature conferring sufficient power to correctly assign offspring to parents. However, microsatellite markers have seen a sharp decline in use with the rise of next-generation sequencing technologies, especially in the study of population genetics and local adaptation. The time is ripe to review the current state of parentage analysis and see how it stands to be affected by the emergence of next-generation sequencing approaches. We find that single nucleotide polymorphisms (SNPs), the typical next-generation sequencing marker, remain underutilized in parentage analysis but are gaining momentum, with 58 SNP-based parentage analyses published thus far. Many of these papers, particularly the earlier ones, compare the power of SNPs and microsatellites in a parentage context. In virtually every case, SNPs are at least as powerful as microsatellite markers. As few as 100-500 SNPs are sufficient to resolve parentage completely in most situations. We also provide an overview of the analytical programs that are commonly used and compatible with SNP data. As the next-generation parentage enterprise grows, a reliance on likelihood and Bayesian approaches, as opposed to strict exclusion, will become increasingly important. We discuss some of the caveats surrounding the use of next-generation sequencing data for parentage analysis and conclude that the future is bright for this important realm of molecular ecology.

Dam trout: Genetic variability in Oncorhynchus mykiss above and below barriers in three Columbia River systems prior to restoring migrational access

Restoration of access to lost habitat for threatened and endangered fishes above currently impassable dams represents a major undertaking. Biological monitoring is critical to understand the dynamics and success of anadromous recolonization as, in the case of Oncorhynchus mykiss, anadromous steelhead populations are reconnected with their conspecific resident rainbow trout counterparts. We evaluate three river systems in the Lower Columbia River basin: the White Salmon, Sandy, and Lewis rivers that are in the process of removing and/or providing passage around existing human-made barriers in O. mykiss riverine habitat. In these instances, now isolated resident rainbow trout populations will be exposed to competition and/or genetic introgression with steelhead and vice versa. Our genetic analyses of 2,158 fish using 13 DNA microsatellite (mSAT) loci indicated that within each basin anadromous O. mykiss were genetically distinct from and significantly more diverse than their resident above-dam trout counterparts. Above long-standing natural impassable barriers, each of these watersheds also harbors unique rainbow trout gene pools with reduced levels of genetic diversity. Despite frequent releases of non-native steelhead and rainbow trout in each river, hatchery releases do not appear to have had a significant genetic effect on the population structure of O. mykiss in any of these watersheds. Simulation results suggest there is a high likelihood of identifying anadromous x resident individuals in the Lewis and White Salmon rivers, and slightly less so in the Sandy River. These genetic data are a prerequisite for informed monitoring, managing, and conserving the different life history forms during upstream recolonization when sympatry of life history forms of O. mykiss is restored.

Life history variation is maintained by fitness trade-offs and negative frequency-dependent selection

The maintenance of diverse life history strategies within and among species remains a fundamental question in ecology and evolutionary biology. By using a near-complete 16-year pedigree of 12,579 winter-run steelhead (Oncorhynchus mykiss) from the Hood River, Oregon, we examined the continued maintenance of two life history traits: the number of lifetime spawning events (semelparous vs. iteroparous) and age at first spawning (2-5 years). We found that repeat-spawning fish had more than 2.5 times the lifetime reproductive success of single-spawning fish. However, first-time repeat-spawning fish had significantly lower reproductive success than single-spawning fish of the same age, suggesting that repeat-spawning fish forego early reproduction to devote additional energy to continued survival. For single-spawning fish, we also found evidence for a fitness trade-off for age at spawning: older, larger males had higher reproductive success than younger, smaller males. For females, in contrast, we found that 3-year-old fish had the highest mean lifetime reproductive success despite the observation that 4- and 5-year-old fish were both longer and heavier. This phenomenon was explained by negative frequency-dependent selection: as 4- and 5-year-old fish decreased in frequency on the spawning grounds, their lifetime reproductive success became greater than that of the 3-year-old fish. Using a combination of mathematical and individual-based models parameterized with our empirical estimates, we demonstrate that both fitness trade-offs and negative frequency-dependent selection observed in the empirical data can theoretically maintain the diverse life history strategies found in this population.

Grandparentage assignments identify unexpected adfluvial life history tactic contributing offspring to a reintroduced population

Diversity in life history tactics contributes to the persistence of a population because it helps to protect against stochastic environments by varying individuals in space and time. However, some life history tactics may not be accounted for when assessing the demographic viability of a population. One important factor in demographic viability assessments is cohort replacement rate (CRR), which is defined as the number of future adults produced by an adult. We assessed if precocial resident males (<age‐3) and adfluvial Chinook salmon (Oncorhynchus tshawytscha), adults that reside in freshwater their entire lives, contributed offspring to a reintroduced population from 2008 to 2013. We found that 9 ± 5% of offspring with an unassigned parent remained unexplained after accounting for sources of human error. Using grandparentage assignments, we identified 31 precocial resident males and 48 probable adfluvial Chinook salmon produced by anadromous mate pairs from 2007 to 2012. Previously published CRR estimates for the 2007 and 2008 reintroduced adults, based on only anadromous returning adult offspring, were 0.40 and 0.31, respectively. By incorporating adfluvial females, we found CRR estimates increased by 17% (CRR: 0.46) and 13% (CRR: 0.35) for the 2007 and 2008 cohorts, respectively.

Bayesian pedigree inference with small numbers of single nucleotide polymorphisms via a factor-graph representation

We develop a computational framework for addressing pedigree inference problems using small numbers (80-400) of single nucleotide polymorphisms (SNPs). Our approach relaxes the assumptions, which are commonly made, that sampling is complete with respect to the pedigree and that there is no genotyping error. It relies on representing the inferred pedigree as a factor graph and invoking the Sum-Product algorithm to compute and store quantities that allow the joint probability of the data to be rapidly computed under a large class of rearrangements of the pedigree structure. This allows efficient MCMC sampling over the space of pedigrees, and, hence, Bayesian inference of pedigree structure. In this paper we restrict ourselves to inference of pedigrees without loops using SNPs assumed to be unlinked. We present the methodology in general for multigenerational inference, and we illustrate the method by applying it to the inference of full sibling groups in a large sample (n=1157) of Chinook salmon typed at 95 SNPs. The results show that our method provides a better point estimate and estimate of uncertainty than the currently best-available maximum-likelihood sibling reconstruction method. Extensions of this work to more complex scenarios are briefly discussed.

Patterns of variation and covariation in the shapes of mandibular bones of juvenile salmonids in the genus Oncorhynchus

What is the nature of evolutionary divergence of the jaw skeleton within the genus Oncorhynchus? How can two associated bones evolve new shapes and still maintain functional integration? Here, we introduce and test a "concordance" hypothesis, in which an extraordinary matching of the evolutionary shape changes of the dentary and angular articular serves to preserve their fitting together. To test this hypothesis, we examined morphologies of the dentary and angular articular at parr (juvenile) stage, and at three levels of biological organization—between salmon and trout, between sister species within both salmon and trout, and among three types differing in life histories within one species, Oncorhynchus mykiss. The comparisons show bone shape divergences among the groups at each level; morphological divergence between salmon and trout is marked even at this relatively early life history stage. We observed substantial matching between the two mandibular bones in both pattern and amount of shape variation, and in shape covariation across species. These findings strongly support the concordance hypothesis, and reflect functional and/or developmental constraint on morphological evolution. We present evidence for developmental modularity within both bones. The locations of module boundaries were predicted from the patterns of evolutionary divergences, and for the dentary, at least, would appear to facilitate its functional association with the angular articular. The modularity results suggest that development has biased the course of evolution.

The effects of riverine physical complexity on anadromy and genetic diversity in steelhead or rainbow trout Oncorhynchus mykiss around the Pacific Rim

This study explored the relationship between riverine physical complexity, as determined from remotely sensed metrics, and anadromy and genetic diversity in steelhead or rainbow trout Oncorhynchus mykiss. The proportion of anadromy (estimated fraction of individuals within a drainage that are anadromous) was correlated with riverine complexity, but this correlation appeared to be driven largely by a confounding negative relationship between drainage area and the proportion of anadromy. Genetic diversity decreased with latitude, was lower in rivers with only non-anadromous individuals and also decreased with an increasing ratio of floodplain area to total drainage area. Anadromy may be less frequent in larger drainages due to the higher cost of migration associated with reaches farther from the ocean, and the negative relationship between genetic diversity and floodplain area may be due to lower effective population size resulting from greater population fluctuations associated with higher rates of habitat turnover. Ultimately, the relationships between riverine physical complexity and migratory life history or genetic diversity probably depend on the spatial scale of analysis.

Molecular pedigree reconstruction and estimation of evolutionary parameters in a wild Atlantic salmon river system with incomplete sampling: a power analysis

Background

Pedigree reconstruction using genetic analysis provides a useful means to estimate fundamental population biology parameters relating to population demography, trait heritability and individual fitness when combined with other sources of data. However, there remain limitations to pedigree reconstruction in wild populations, particularly in systems where parent-offspring relationships cannot be directly observed, there is incomplete sampling of individuals, or molecular parentage inference relies on low quality DNA from archived material. While much can still be inferred from incomplete or sparse pedigrees, it is crucial to evaluate the quality and power of available genetic information a priori to testing specific biological hypotheses. Here, we used microsatellite markers to reconstruct a multi-generation pedigree of wild Atlantic salmon (Salmo salar L.) using archived scale samples collected with a total trapping system within a river over a 10 year period. Using a simulation-based approach, we determined the optimal microsatellite marker number for accurate parentage assignment, and evaluated the power of the resulting partial pedigree to investigate important evolutionary and quantitative genetic characteristics of salmon in the system.

Results

We show that at least 20 microsatellites (ave. 12 alleles/locus) are required to maximise parentage assignment and to improve the power to estimate reproductive success and heritability in this study system. We also show that 1.5 fold differences can be detected between groups simulated to have differing reproductive success, and that it is possible to detect moderate heritability values for continuous traits (h² ~ 0.40) with more than 80% power when using 28 moderately to highly polymorphic markers.

Conclusion

The methodologies and work flow described provide a robust approach for evaluating archived samples for pedigree-based research, even where only a proportion of the total population is sampled. The results demonstrate the feasibility of pedigree-based studies to address challenging ecological and evolutionary questions in free-living populations, where genealogies can be traced only using molecular tools, and that significant increases in pedigree assignment power can be achieved by using higher numbers of markers.

Gene flow between sympatric life history forms of Oncorhynchus mykiss located above and below migratory barriers

Oncorhynchus mykiss have a diverse array of life history types, and understanding the relationship among types is important for management of the species. Patterns of gene flow between sympatric freshwater resident O. mykiss, commonly known as rainbow trout, and anadromous O. mykiss, commonly known as steelhead, populations are complex and poorly understood. In this study, we attempt to determine the occurrence and pathways of gene flow and the degree of genetic similarity between sympatric resident and anadromous O. mykiss in three river systems, and investigate whether resident O. mykiss are producing anadromous offspring in these rivers, two of which have complete barriers to upstream migration. We found that the population structure of the O. mykiss in these rivers appears to be influenced more by the presence of a barrier to upstream migration than by life history type. The sex ratio of resident O. mykiss located above a barrier, and smolts captured in screw traps was significantly skewed in favor of females, whereas the reverse was true below the barriers, suggesting that male resident O. mykiss readily migrate downstream over the barrier, and that precocious male maturation may be occurring in the anadromous populations. Through paternity analyses, we also provide direct confirmation that resident O. mykiss can produce offspring that become anadromous. Most (89%) of the resident O. mykiss that produced anadromous offspring were males. Our results add to the growing body of evidence that shows that gene flow does readily occur between sympatric resident and anadromous O. mykiss life history types, and indicates that resident O. mykiss populations may be a potential repository of genes for the anadromous life history type.

How much does inbreeding contribute to the reduced fitness of hatchery-born steelhead (Oncorhynchus mykiss) in the wild?

Many declining populations are supplemented with captive-born individuals that are released directly into the wild. Because captive-born individuals can have lower fitness in the wild than their wild-born counterparts, a comprehensive understanding of the mechanisms responsible for the reduced fitness of these individuals is required for appropriate conservation and management decisions. Inbreeding among captive-born individuals is one plausible mechanism because captive breeding programs frequently use small numbers of breeders to create large numbers of siblings that are subsequently released together into the wild. We tested this hypothesis in a supplementation program for steelhead (Oncorhynchus mykiss) from the Hood River, Oregon, for which first-generation hatchery fish were demonstrated to have lower fitness in the wild than their wild-born counterparts. To determine the contribution of inbreeding to this fitness decline, we first assigned 11 run-years of hatchery steelhead (3005 fish) back to their broodstock parents (462 fish) using 8 polymorphic microsatellite loci. By combining pedigree analyses with species-specific estimates of genetic load, we found that inbreeding could at most account for a 1-4% reduction in the fitness of hatchery fish relative to wild fish. Thus, inbreeding alone cannot adequately explain the 15% average fitness decline observed in first-generation hatchery fish from this population.

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.

Can interbreeding of wild and artificially propagated animals be prevented by using broodstock selected for a divergent life history?

Two strategies have been proposed to avoid negative genetic effects of artificially propagated individuals on wild populations: (i) integration of wild and captive populations to minimize domestication selection and (ii) segregation of released individuals from the wild population to minimize interbreeding. We tested the efficacy of the strategy of segregation by divergent life history in a steelhead trout, Oncorhynchus mykiss, system, where hatchery fish were selected to spawn months earlier than the indigenous wild population. The proportion of wild ancestry smolts and adults declined by 10–20% over the three generations since the hatchery program began. Up to 80% of the naturally produced steelhead in any given year were hatchery/wild hybrids. Regression model selection analysis showed that the proportion of hatchery ancestry smolts was lower in years when stream discharge was high, suggesting a negative effect of flow on reproductive success of early-spawning hatchery fish. Furthermore, proportions of hybrid smolts and adults were higher in years when the number of naturally spawning hatchery-produced adults was higher. Divergent life history failed to prevent interbreeding when physical isolation was ineffective, an inadequacy that is likely to prevail in many other situations.

Signatures of natural selection among lineages and habitats in Oncorhynchus mykiss

Recent advances in molecular interrogation techniques now allow unprecedented genomic inference about the role of adaptive genetic divergence in wild populations. We used high-throughput genotyping to screen a genome-wide panel of 276 single nucleotide polymorphisms (SNPs) for the economically and culturally important salmonid Oncorhynchus mykiss. Samples included 805 individuals from 11 anadromous and resident populations from the northwestern United States and British Columbia, and represented two major lineages including paired populations of each life history within single drainages of each lineage. Overall patterns of variation affirmed clear distinctions between lineages and in most instances, isolation by distance within them. Evidence for divergent selection at eight candidate loci included significant landscape correlations, particularly with temperature. High diversity of two nonsynonymous mutations within the peptide-binding region of the major histocompatibility complex (MHC) class II (DAB) gene provided signatures of balancing selection. Weak signals for potential selection between sympatric resident and anadromous populations were revealed from genome scans and allele frequency comparisons. Our results suggest an important adaptive role for immune-related functions and present a large genomic resource for future studies

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).

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.

A new version of PRT software for sibling groups reconstruction with comments regarding several issues in the sibling reconstruction problem

Pedigree reconstruction using genotypic markers has become an important tool for the study of natural populations. The nonstandard nature of the underlying statistical problems has led to the necessity of developing specialized statistical and computational methods. In this article, a new version of pedigree reconstruction tools (PRT 2.0) is presented. The software implements algorithms proposed in Almudevar & Field (Journal of Agricultural Biological and Environmental Statistics, 4, 1999, 136) and Almudevar (Biometrics, 57, 2001a, 757) for the reconstruction of single generation sibling groups (SG). A wider range of enumeration algorithms is included, permitting improved computational performance. In particular, an iterative version of the algorithm designed for larger samples is included in a fully automated form. The new version also includes expanded simulation utilities, as well as extensive reporting, including half-sibling compatibility, parental genotype estimates and flagging of potential genotype errors. A number of alternative algorithms are described and demonstrated. A comparative discussion of the underlying methodologies is presented. Although important aspects of this problem remain open, we argue that a number of methodologies including maximum likelihood estimation (COLONY 1.2 and 2.0) and the set cover formulation (KINALYZER) exhibit undesirable properties in the sibling reconstruction problem. There is considerable evidence that large sets of individuals not genetically excluded as siblings can be inferred to be a true sibling group, but it is also true that unrelated individuals may be genetically compatible with a true sibling group by chance. Such individuals may be identified on a statistical basis. PRT 2.0, based on these sound statistical principles, is able to efficiently match or exceed the highest reported accuracy rates, particularly for larger SG. The new version is available at http://www.urmc.rochester.edu/biostat/people/faculty/almudevar.cfm.

Little impact of hatchery supplementation that uses native broodstock on the genetic structure and diversity of steelhead trout revealed by a large-scale spatio-temporal microsatellite survey

Artificial breeding programs initiated to enhance the size of animal populations are often motivated by the desire to increase harvest opportunities. The introduction of non-native genotypes, however, can have negative evolutionary impacts. These may be direct, such as introgressive hybridization, or indirect via competition. Less is known about the effects of stocking with native genotypes. We assayed variation at nine microsatellite loci in 902 steelhead trout (Oncorhynchus mykiss) from five rivers in British Columbia, Canada. These samples were collected over 58 years, a time period that spanned the initiation of native steelhead trout broodstock hatchery supplementation in these rivers. We detected no changes in estimates of effective population size, genetic variation or temporal genetic structure within any population, nor of altered genetic structure among them. Genetic interactions with nonmigratory O. mykiss, the use of substantial numbers of primarily native broodstock with an approximate 1:1 male-to-female ratio, and/or poor survival and reproductive success of hatchery fish may have minimized potential genetic changes. Although no genetic changes were detected, ecological effects of hatchery programs still may influence wild population productivity and abundance. Their effects await the design and implementation of a more comprehensive evaluation program.

Grandfathering in a new era of parentage analysis

The advent of DNA fingerprinting and microsatellite techniques has revolutionized the way in which we investigate genetic pedigrees in the wild (Pemberton 2008). With large and often incomplete data sets consisting of hundreds to thousands of individuals over multiple generations becoming commonplace, new methods in parentage analysis are being developed to rise to the next generation of questions and challenges. In this issue, Christie et al. (2011) provide a simple yet elegant solution to the problem of identifying missing parents and assessing hybrid fitness in a mixed population of wild and hatchery steelhead trout (Oncorhynchus mykiss) where not all individuals can be sampled effectively. They develop a new method of grandparent analysis where parental genotypes can be reconstructed using data from candidate grandparent crosses and F2 offspring genotypes, allowing for new explorations of hybridization, migration and gene flow in wild populations.