Measuring individual inbreeding in the age of genomics: marker-based measures are better than pedigrees.

Inbreeding (mating between relatives) can dramatically reduce the fitness of offspring by causing parts of the genome to be identical by descent. Thus, measuring individual inbreeding is crucial for ecology, evolution and conservation biology. We used computer simulations to test whether the realized proportion of the genome that is identical by descent (IBDG) is predicted better by the pedigree inbreeding coefficient (FP) or by genomic (marker-based) measures of inbreeding. Genomic estimators of IBDG included the increase in individual homozygosity relative to mean Hardy-Weinberg expected homozygosity (FH), and two measures (FROH and FE) that use mapped genetic markers to estimate IBDG. IBDG was more strongly correlated with FH, FE and FROH than with FP across a broad range of simulated scenarios when thousands of SNPs were used. For example, IBDG was more strongly correlated with FROH, FH and FE (estimated with ⩾10 000 SNPs) than with FP (estimated with 20 generations of complete pedigree) in populations with a recent reduction in the effective populations size (from Ne=500 to Ne=75). FROH, FH and FE generally explained >90% of the variance in IBDG (among individuals) when 35 K or more SNPs were used. FP explained <80% of the variation in IBDG on average in all simulated scenarios, even when pedigrees included 20 generations. Our results demonstrate that IBDG can be more precisely estimated with large numbers of genetic markers than with pedigrees. We encourage researchers to adopt genomic marker-based measures of IBDG as thousands of loci can now be genotyped in any species.

Estimation of effective population size in continuously distributed populations: there goes the neighborhood.

Use of genetic methods to estimate effective population size (Ne) is rapidly increasing, but all approaches make simplifying assumptions unlikely to be met in real populations. In particular, all assume a single, unstructured population, and none has been evaluated for use with continuously distributed species. We simulated continuous populations with local mating structure, as envisioned by Wright's concept of neighborhood size (NS), and evaluated performance of a single-sample estimator based on linkage disequilibrium (LD), which provides an estimate of the effective number of parents that produced the sample (Nb). Results illustrate the interacting effects of two phenomena, drift and mixture, that contribute to LD. Samples from areas equal to or smaller than a breeding window produced estimates close to the NS. As the sampling window increased in size to encompass multiple genetic neighborhoods, mixture LD from a two-locus Wahlund effect overwhelmed the reduction in drift LD from incorporating offspring from more parents. As a consequence, never approached the global Ne, even when the geographic scale of sampling was large. Results indicate that caution is needed in applying standard methods for estimating effective size to continuously distributed populations.

Species and cultural conservation in New Zealand: maori traditional ecological knowledge of tuatara.

Traditional ecological knowledge can be highly informative and integrated with complementary scientific knowledge to improve species management. This is especially true for abundant species with which indigenous peoples have frequent interactions (e.g., through harvest), but has been studied less frequently in isolated or declining species. We examined Maori traditional ecological knowledge of tuatara (Sphenodon spp., reptiles that resemble lizards but are the last living representatives of the order Sphenodontia) through semidirected interviews of elders of Te Atiawa, Ngati Koata, and Ngati Wai Iwi (similar to tribes), the guardians of several islands currently inhabited by tuatara. Maori are indigenous to New Zealand, having settled 800 to 1000 years ago. Tuatara are endemic to New Zealand, have declined in numbers since human settlement, and are now restricted to 37 offshore islands. The detail and volume of tuatara traditional ecological knowledge were less than that recorded in studies of more abundant or accessible species. In addition, traditional knowledge of the cultural significance of tuatara was more common and detailed among the elders than traditional knowledge of tuatara biology or ecology. The traditional knowledge collected, however, provided the first evidence of seven former sites of tuatara occupation, suggested five additional sites tuatara may currently occupy, contained novel hypotheses for scientific testing, and described tuatara cultural roles that have not been reported previously. We conclude that, in at least some cases, traditional ecological knowledge may persist as species decline and may serve as a valuable source of ecological information for conservation.

Genetic basis of variation in morphological and life-history traits of a wild population of pink salmon.

Understanding the genetic basis of phenotypic variation is essential for predicting the direction and rate of phenotypic evolution. We estimated heritabilities and genetic correlations of morphological (fork length, pectoral and pelvic fin ray counts, and gill arch raker counts) and life-history (egg number and individual egg weight) traits of pink salmon (Oncorhynchus gorbuscha) from Likes Creek, Alaska, in order to characterize the genetic basis of phenotypic variation in this species. Families were created from wild-caught adults, raised to the fry stage in the lab, released into the wild, and caught as returning adults and assigned to families using microsatellite loci and a growth hormone locus. Morphological traits were all moderately to highly heritable, but egg number and egg weight were not heritable, suggesting that past selection has eliminated additive genetic variation in egg number and egg weight or that there is high environmental variance in these traits. Genetic correlations were similar for nonadjacent morphological traits and adjacent traits. Genetic correlations predicted phenotypic correlations fairly accurately, but some pairs of traits with low genetic correlations had high phenotypic correlations, and vice versa, emphasizing the need to use caution when using phenotypic correlations as indices of genetic correlations. This is one of only a handful of studies to estimate heritabilities and genetic correlations for a wild population.

Linkage disequilibrium between the pseudoautosomal PEPB-1 locus and the sex-determining region of chinook salmon.

Allele frequency differences between sexes and an excess of heterozygotes in males had suggested that the PEPB-1 locus is sex-linked in chinook salmon (Oncorhynchus tshawytscha). We here estimate less than 1% recombination between PEPB-1 and a growth hormone pseudogene known to be in the sex-determining region (SEX) in 374 progeny from eight experimental matings. We present modified maximum likelihood methods for estimating haplotype frequencies from population samples at a sex-linked locus in which functional alleles occur on both the X and Y chromosomes (pseudoautosomal loci). We find nearly complete linkage disequilibrium between PEPB-1 and SEX in 20 population samples from the Puget Sound region of Washington and southern British Columbia. However, allele frequencies at PEPB-1 were similar in males and females in 35 population samples from the coast of Washington and the Columbia River basin. Pseudoautosomal regions have been described in a broad taxonomic array of vertebrates and invertebrates, and they are likely candidate regions to find genes associated with differences in life history, morphology, or behavior between males and females.

Founding events influence genetic population structure of sockeye salmon (Oncorhynchus nerka) in Lake Clark, Alaska.

Bottlenecks can have lasting effects on genetic population structure that obscure patterns of contemporary gene flow and drift. Sockeye salmon are vulnerable to bottleneck effects because they are a highly structured species with excellent colonizing abilities and often occupy geologically young habitats. We describe genetic divergence among and genetic variation within spawning populations of sockeye salmon throughout the Lake Clark area of Alaska. Fin tissue was collected from sockeye salmon representing 15 spawning populations of Lake Clark, Six-mile Lake, and Lake Iliamna. Allele frequencies differed significantly at 11 microsatellite loci in 96 of 105 pairwise population comparisons. Pairwise estimates of FST ranged from zero to 0.089. Six-mile Lake and Lake Clark populations have historically been grouped together for management purposes and are geographically proximate. However, Six-mile Lake populations are genetically similar to Lake Iliamna populations and are divergent from Lake Clark populations. The reduced allelic diversity and strong divergence of Lake Clark populations relative to Six-mile Lake and Lake Iliamna populations suggest a bottleneck associated with the colonization of Lake Clark by sockeye salmon. Geographic distance and spawning habitat differences apparently do not contribute to isolation and divergence among populations. However, temporal isolation based on spawning time and founder effects associated with ongoing glacial retreat and colonization of new spawning habitats contribute to the genetic population structure of Lake Clark sockeye salmon. Nonequilibrium conditions and the strong influence of genetic drift caution against using estimates of divergence to estimate gene flow among populations of Lake Clark sockeye salmon.

Landscape location affects genetic variation of Canada lynx (Lynx canadensis).

The effect of a population's location on the landscape on genetic variation has been of interest to population genetics for more than half a century. However, most studies do not consider broadscale biogeography when interpreting genetic data. In this study, we propose an operational definition of a peripheral population, and then explore whether peripheral populations of Canada lynx (Lynx canadensis) have less genetic variation than core populations at nine microsatellite loci. We show that peripheral populations of lynx have fewer mean numbers of alleles per population and lower expected heterozygosity. This is surprising, given the lynx's capacity to move long distances, but can be explained by the fact that peripheral populations often have smaller population sizes, limited opportunities for genetic exchange and may be disproportionately affected by ebbs and flows of species' geographical range.

A consolidated linkage map for rainbow trout (Oncorhynchus mykiss).

Androgenetic doubled haploid progeny produced from a cross between the Oregon State University and Arlee clonal rainbow trout (Oncorhynchus mykiss) lines, used for a previous published rainbow trout map, were used to update the map with the addition of more amplified fragment length polymorphic (AFLP) markers, microsatellites, type I and allozyme markers. We have added more than 900 markers, bringing the total number to 1359 genetic markers and the sex phenotype including 799 EcoRI AFLPs, 174 PstI AFLPs, 226 microsatellites, 72 VNTR, 38 SINE markers, 29 known genes, 12 minisatellites, five RAPDs, and four allozymes. Thirty major linkage groups were identified. Synteny of linkage groups in our map with the outcrossed microsatellite map has been established for all except one linkage group in this doubled haploid cross. Putative homeologous relationships among linkage groups, resulting from the autotetraploid nature of the salmonid genome, have been revealed based on the placement of duplicated microsatellites and type I loci.

Insights into recently fragmented vole populations from combined genetic and demographic data.

We combined demographic and genetic data to evaluate the effects of habitat fragmentation on the population structure of the California red-backed vole (Clethrionomys californicus). We analysed variation in the mitochondrial DNA (mtDNA) control region and five nuclear microsatellite loci in small samples collected from two forest fragments and an unfragmented control site in 1990-91. We intensively sampled the same forest fragments and two different control sites in 1998 and 1999. Vole abundances fluctuated greatly at sizes below 50 individuals per fragment. Fragment populations had significantly lower mtDNA allelic diversity than controls, but not nuclear heterozygosity or numbers of alleles. The use of only trapping and/or mtDNA marker data would imply that fragment populations are at least partially isolated and vulnerable to inbreeding depression. In contrast, the abundance estimates combined with microsatellite data show that small fragment populations must be linked to nearby forests by high rates of migration. These results provide evidence for the usefulness of combining genetic and demographic data to understand nonequilibrium population structure in recently fragmented habitats.

Heterozygosity-fitness correlations in rainbow trout: effects of allozyme loci or associative overdominance?

Previous studies with rainbow trout (Oncorhynchus mykiss) have shown that increased heterozygosity at allozyme loci is correlated with several phenotypic traits associated with fitness. We expected to find a similar effect of heterozygosity at other nuclear loci if these associations are due to loci in linkage disequilibrium with the allozyme loci (i.e., associative overdominance), rather than the allozymes themselves. We examined the association between multiple locus heterozygosity and condition factor at 10 allozyme and 10 microsatellite loci. Individuals that were more heterozygous at allozyme loci had significantly greater condition factor in two hatchery cohorts of rainbow trout (1996 P = 0.006; 1997 P < 0.001). In contrast, there was no evidence at microsatellite loci that increased heterozygosity was associated with greater condition factor. Our results suggest that the observed relationship between heterozygosity at allozyme loci and condition factor in rainbow trout appears to be due to the allozyme loci themselves, rather than associative overdominance. We cannot, however, rule out that differences in the mutation process between allozymes and microsatellites may be responsible for these observations. Regardless of the underlying mechanism, these results support the view that allozymes and microsatellites are differentially affected by natural selection.

Mitochondrial mutations may decrease population viability.

It is ironic that an organelle that is pivotal for the function of male gametes is inherited by sons only from their mother. A recent study of human fertility has provided the first definitive evidence that mitochondrial DNA mutations can reduce male fertility but with little or no effect on females. These results present the possibility that the viability of small populations might be reduced by increases in the frequency of mtDNA genotypes that lower the fitness of males.

Concordance of genetic divergence among sockeye salmon populations at allozyme, nuclear DNA, and mitochondrial DNA markers.

We examined genetic variation at 21 polymorphic allozyme loci, 15 nuclear DNA loci, and mitochondrial DNA in four spawning populations of sockeye salmon (Oncorhynchus nerka) from Cook Inlet, Alaska, to test for differences in the patterns of divergence among different types of markers. We were specifically interested in testing the suggestion that natural selection at allozyme loci compromises the effectiveness of these markers for describing the amount and patterns of gene flow among populations. We found concordance among markers in the amount of genetic variation within and among populations, with the striking exception of one allozyme locus (sAH), which exhibited more than three times the amount of among-population differentiation as other loci. A consideration of reports of discordance between allozymes and other loci indicates that these differences usually result from one or two exceptional loci. We conclude that it is important to examine many loci when estimating genetic differentiation to infer historical amounts of gene flow and patterns of genetic exchange among populations. It is less important whether those loci are allozymes or nuclear DNA markers.

Gene-centromere mapping of 312 loci in pink salmon by half-tetrad analysis.

We estimated recombination rates between 312 loci and their centromeres in gynogenetic diploid pink salmon (Oncorhynchus gorbuscha) that we produced by initiating development with irradiated sperm and blocking the maternal second meiotic division. Amplified fragment length polymorphisms (AFLPs) were significantly more centromeric than loci identified by three other techniques (allozymes, microsatellites, and PCR using primer sequences from interspersed nuclear elements). The near absence of AFLPs in distal regions could limit their utility in constructing linkage maps. A large proportion of loci had frequency of second division segregation (y) values approaching 1.0, indicating near complete crossover interference on many chromosome arms. As predicted from models of chromosomal evolution in salmonids based upon results with allozyme loci, all duplicated microsatellite loci that shared alleles (isoloci) had y values of nearly 1.0.

Small effective population size in the long-toed salamander.

The effective population sizes (Ne) of six populations of the long-toed salamander (Ambystoma macrodactylum) from Montana and Idaho, USA were estimated from allozyme data from samples collected in 1978, 1996 and 1997 using the temporal allele frequency method. Five of the six estimates ranged from 23 to 207 (mean = 123 +/- 79); one estimate was indistinguishable from infinity. In order to infer the actual Ne of salamander populations, we compared the frequency distribution of our observed Ne estimates with distributions obtained from simulated populations of known Ne. Our observed Ne estimate distribution was consistent with distributions from simulated populations with Ne values of 10, 25, and 50, suggesting an actual Ne for each of the six salamander populations of less than 100. This Ne estimate agrees with most other Ne estimates for amphibians. We conclude by discussing the conservation implications of small Ne values in amphibians in the context of increasing isolation of populations due to habitat fragmentation.

Inheritance of nuclear DNA markers in gynogenetic haploid pink salmon.

We describe the inheritance of 460 PCR-based loci in the polyploid-derived pink salmon (Oncorhynchus gorbuscha) genome using gynogenetic haploid embryos. We detected a length polymorphism in a growth hormone gene (GH-2) intron that is caused by an 81 bp insertion homologous to the 3' end of the salmonid short interspersed repetitive element (SINE) SmaI. Such insertion polymorphisms within species bring into question the use of SINEs as phylogenetic markers. We confirmed that a microsatellite locus encodes a PCR-null allele that is responsible for an apparent deficit of heterozygotes in a population sample from Prince William Sound. Another set of microsatellite primers amplified alleles of the same molecular weight from both loci of a duplicated pair. In our analysis of several PCR-based multilocus techniques, we failed to detect evidence of comigrating fragments produced by duplicated loci. Segregation analysis of PCR-based markers using gynogenetic haploid embryos ensures that the interpretation of molecular variation is not complicated by heterozygosity, diploidy, or gene duplication. We urge investigators to test the inheritance of polymorphisms in salmonids prior to using them to measure genetic variation.

Usefulness of molecular markers for detecting population bottlenecks via monitoring genetic change.

It is important to detect population bottlenecks in threatened and managed species because bottlenecks can increase the risk of population extinction. Early detection is critical and can be facilitated by statistically powerful monitoring programs for detecting bottleneck-induced genetic change. We used Monte Carlo computer simulations to evaluate the power of the following tests for detecting genetic changes caused by a severe reduction in a population's effective size (Ne): a test for loss of heterozygosity, two tests for loss of alleles, two tests for change in the distribution of allele frequencies, and a test for small Ne based on variance in allele frequencies (the 'variance test'). The variance test was most powerful; it provided an 85% probability of detecting a bottleneck of size Ne = 10 when monitoring five microsatellite loci and sampling 30 individuals both before and one generation after the bottleneck. The variance test was almost 10-times more powerful than a commonly used test for loss of heterozygosity, and it allowed for detection of bottlenecks before 5% of a population's heterozygosity had been lost. The second most powerful tests were generally the tests for loss of alleles. However, these tests had reduced power for detecting genetic bottlenecks caused by skewed sex ratios. We provide guidelines for the number of loci and individuals needed to achieve high-power tests when monitoring via the variance test. We also illustrate how the variance test performs when monitoring loci that have widely different allele frequency distributions as observed in five wild populations of mountain sheep (Ovis canadensis).

Distortion of allele frequency distributions provides a test for recent population bottlenecks.

We use population genetics theory and computer simulations to demonstrate that population bottlenecks cause a characteristic mode-shift distortion in the distribution of allele frequencies at selectively neutral loci. Bottlenecks cause alleles at low frequency (< 0.1) to become less abundant than alleles in one or more intermediate allele frequency class (e.g., 0.1-0.2). This distortion is transient and likely to be detectable for only a few dozen generations. Consequently only recent bottlenecks are likely to be detected by tests for distortions in distributions of allele frequencies. We illustrate and evaluate a qualitative graphical method for detecting a bottleneck-induced distortion of allele frequency distributions. The simple novel method requires no information on historical population sizes or levels of genetic variation; it requires only samples of 5 to 20 polymorphic loci and approximately 30 individuals. The graphical method often differentiates between empirical datasets from bottlenecked and nonbottlenecked natural populations. Computer simulations show that the graphical method is likely (P > .80) to detect an allele frequency distortion after a bottleneck of < or = 20 breeding individuals when 8 to 10 polymorphic microsatellite loci are analyzed.

Secondary tetrasomic segregation of MDH-B and preferential pairing of homeologues in rainbow trout.

We examined the inheritance of allelic variation at an isozyme locus, MDH-B, duplicated by ancestral polyploidy in salmonid fishes. We detected only disomic segregation in females. Segregation ratios in males were best explained by a mixture of disomic and tetrasomic inheritance. We propose a two-stage model of pairing in male meiosis in which, first, homologous chromosomes pair and recombine in the proximal region of the chromosome. Next, homeologous chromosomes pair and recombine distally. We suggest that this type of tetrasomic inheritance in which centromeres segregate disomically should be referred to as "secondary tetrasomy" to distinguish it from tetrasomy involving entire chromosomes (i.e., "primary tetrasomy"). Differences in segregation ratios between males indicate differences between individuals in the amount of recombination between homeologous chromosomes. We also consider the implication of these results for estimation of allele frequencies at duplicated loci in salmonid populations.

Microsatellite evolution in congeneric mammals: domestic and bighorn sheep.

We compared genotypes at eight (AC)n microsatellite loci in domestic sheep (Ovis aries) and wild Rocky Mountain bighorn sheep (O. canadensis). The domestic sheep had greater genetic variation, higher allele-size variances, and larger allele sizes than the wild sheep. Accumulating evidence from higher taxonomic comparisons shows that these parameters are biased if microsatellite loci are selected in one taxon and used in another. Our results demonstrate similar biases between congeneric species. We compared standard measures of genetic variation, differentiation, and distance within and between species (H, D, FST) to newer measures based on allele-size variance (SW, SB, RST). The size-based distances better detected species-level divergence, but standard measures better distinguished allopatric populations. Empirical calibration of these measures at the subspecies level is needed to establish their useful ranges.