Joseph I. Hoffman

Microsatellite genotyping errors: detection approaches, common sources and consequences for paternal exclusion

Microsatellite genotyping errors will be present in all but the smallest data sets and have the potential to undermine the conclusions of most downstream analyses. Despite this, little rigorous effort has been made to quantify the size of the problem and to identify the commonest sources of error. Here, we use a large data set comprising almost 2000 Antarctic fur seals Arctocephalus gazella genotyped at nine hypervariable microsatellite loci to explore error detection methods, common sources of error and the consequences of errors on paternal exclusion. We found good concordance among a range of contrasting approaches to error‐rate estimation, our range being 0.0013 to 0.0074 per single locus PCR (polymerase chain reaction). The best approach probably involves blind repeat‐genotyping, but this is also the most labour‐intensive. We show that several other approaches are also effective at detecting errors, although the most convenient alternative, namely mother–offspring comparisons, yielded the lowest estimate of the error rate. In total, we found 75 errors, emphasizing their ubiquitous presence. The most common errors involved the misinterpretation of allele banding patterns (n = 60, 80%) and of these, over a third (n = 22, 36.7%) were due to confusion between homozygote and adjacent allele heterozygote genotypes. A specific test for whether a data set contains the expected number of adjacent allele heterozygotes could provide a useful tool with which workers can assess the likely size of the problem. Error rates are also positively correlated with both locus polymorphism and product size, again indicating aspects where extra effort at error reduction should be directed. Finally, we conducted simulations to explore the potential impact of genotyping errors on paternity exclusion. Error rates as low as 0.01 per allele resulted in a rate of false paternity exclusion exceeding 20%. Errors also led to reduced estimates of male reproductive skew and increases in the numbers of pups that matched more than one candidate male. Because even modest error rates can be strongly influential, we recommend that error rates should be routinely published and that researchers make an attempt to calculate how robust their analyses are to errors.

Female fur seals show active choice for males that are heterozygous and unrelated.

Much debate surrounds the exact rules that influence mating behaviour, and in particular the selective forces that explain the evolution of female preferences. A key example is the lek paradox, in which female choice is expected rapidly to become ineffective owing to loss of additive genetic variability for the preferred traits. Here we exploit a remarkable system in which female fur seals exert choice by moving across a crowded breeding colony to visit largely static males. We show that females move further to maximize the balance between male high multilocus heterozygosity and low relatedness. Such a system shows that female choice can be important even in a strongly polygynous species, and at the same time may help to resolve the lek paradox because heterozygosity has low heritability and inbreeding avoidance means there is no single ‘best’ male for all females.

High-throughput sequencing reveals inbreeding depression in a natural population

Significance Many studies of wild populations reveal links between heterozygosity and fitness, with relatively heterozygous individuals carrying fewer parasites, living longer and being more attractive to mates. These patterns appear ubiquitous and are often highly significant, but heterozygosity usually accounts for very little of the total variation in fitness. However, most studies analyze only around 10 loci, representing a tiny fraction of the genome. We therefore used high-throughput DNA sequencing to estimate genome-wide heterozygosity based on over 10,000 loci and found it to accurately reflect inbreeding. Applied to wild harbor seals, genome-wide heterozygosity explained almost half of the variation in parasite infection. By implication, a greater proportion of fitness variation could be linked to genotype than previously thought. Proxy measures of genome-wide heterozygosity based on approximately 10 microsatellites have been used to uncover heterozygosity fitness correlations (HFCs) for a wealth of important fitness traits in natural populations. However, effect sizes are typically very small and the underlying mechanisms remain contentious, as a handful of markers usually provides little power to detect inbreeding. We therefore used restriction site associated DNA (RAD) sequencing to accurately estimate genome-wide heterozygosity, an approach transferrable to any organism. As a proof of concept, we first RAD sequenced oldfield mice (Peromyscus polionotus) from a known pedigree, finding strong concordance between the inbreeding coefficient and heterozygosity measured at 13,198 single-nucleotide polymorphisms (SNPs). When applied to a natural population of harbor seals (Phoca vitulina), a weak HFC for parasite infection based on 27 microsatellites strengthened considerably with 14,585 SNPs, the deviance explained by heterozygosity increasing almost fivefold to a remarkable 49%. These findings arguably provide the strongest evidence to date of an HFC being due to inbreeding depression in a natural population lacking a pedigree. They also suggest that under some circumstances heterozygosity may explain far more variation in fitness than previously envisaged.

Mitochondrial DNA barcoding detects some species that are real, and some that are not

Mimicry and extensive geographical subspecies polymorphism combine to make species in the ithomiine butterfly genus Mechanitis (Lepidoptera; Nymphalidae) difficult to determine. We use mitochondrial DNA (mtDNA) barcoding, nuclear sequences and amplified fragment length polymorphism (AFLP) genotyping to investigate species limits in this genus. Although earlier biosystematic studies based on morphology described only four species, mtDNA barcoding revealed eight well‐differentiated haplogroups, suggesting the presence of four new putative ‘cryptic species’. However, AFLP markers supported only one of these four new ‘cryptic species’ as biologically meaningful. We demonstrate that in this genus, deep genetic divisions expected on the basis of mtDNA barcoding are not always reflected in the nuclear genome, and advocate the use of AFLP markers as a check when mtDNA barcoding gives unexpected results.

Deep genetic subdivision within a continuously distributed and highly vagile marine mammal, the Steller's sea lion ( Eumetopias jubatus )

The Stellers sea lion Eumetopias jubatus is an endangered marine mammal that has experienced dramatic population declines over much of its range during the past five decades. Studies using mitochondrial DNA (mtDNA) have shown that an apparently continuous population includes a strong division, yielding two discrete stocks, western and eastern. Based on a weaker split within the western stock, a third Asian stock has also been defined. While these findings indicate strong female philopatry, a recent study using nuclear microsatellite markers found little evidence of any genetic structure, implying extensive paternal gene flow. However, this result was at odds with mark–recapture data, and both sample sizes and genetic resolution were limited. To address these concerns, we increased analytical power by genotyping over 700 individuals from across the species’ range at 13 highly polymorphic microsatellite loci. We found a clear phylogenetic break between populations of the eastern stock and those of the western and Asian stocks. However, our data provide little support for the classification of a separate Asian stock. Our findings show that mtDNA structuring is not due simply to female philopatry, but instead reflects a genuine discontinuity within the range, with implications for both the phylogeography and conservation of this important marine mammal.

Dispersal, philopatry and intergroup relatedness: fine-scale genetic structure in the white-breasted thrasher, Ramphocinclus brachyurus.

Dispersal is a fundamental process influencing evolution, social behaviour, and the long‐term persistence of populations. We use both observational and genetic data to investigate dispersal, kin‐clustering and intergroup relatedness in the white‐breasted thrasher, Ramphocinclus brachyurus, a cooperatively breeding bird that is globally endangered. Mark‐resighting data suggested sex‐biased dispersal, with females dispersing over greater distances while males remained philopatric. Accordingly, spatial autocorrelation analysis showed highly significant fine‐scale genetic structure among males, but not among females. This fine‐scale genetic structuring of the male population resulted in very high levels of relatedness between dominant males at neighbouring nests, similar to that seen within cooperative groups in many species where kin selection is cited as a cause of cooperation. By implication, between‐group as well as within‐group cooperation may be important, potentially creating a feedback loop in which short‐distance dispersal by males leads to the formation of male kin clusters that in turn facilitate nepotistic interactions and favour further local recruitment. The strength of spatial autocorrelation, as measured by the autocorrelation coefficient, r, was approximately two to three times greater than that reported in previous studies of animals. Relatively short dispersal distances by both males and females may have a negative impact on the white‐breasted thrashers ability to colonize new areas, and may influence the long‐term persistence of isolated populations. This should be taken into account when designating protected areas or selecting sites for habitat restoration.

EXPLORING THE RELATIONSHIP BETWEEN PARENTAL RELATEDNESS AND MALE REPRODUCTIVE SUCCESS IN THE ANTARCTIC FUR SEAL ARCTOCEPHALUS GAZELLA

Abstract Recent genetic studies of natural populations have shown that heterozygosity and other genetic estimates of parental relatedness correlate with a wide variety of fitness traits, from juvenile survival and parasite resistance to male reproductive success. Many of these traits involve health and survival, where the underlying mechanism may involve changes in the effectiveness of the immune system. However, for traits such as reproductive success, the likely mechanisms remain less obvious. In this paper, we examine the relationship between heterozygosity and a range of traits that contribute to male reproductive success, including time spent on territories and competitiveness. Our analysis is based on observational and genetic data from eight consecutive breeding seasons at a colony of the Antarctic fur seal, Arctocephalus gazella. Overall, male reproductive success was found to correlate strongly with internal relatedness (IR, a form of heterozygosity). When different components of success were analyzed, we found that IR correlates independently with reproductive longevity, time spent ashore, and competitive ability per unit mating opportunity on the study beach, with more heterozygous males being more successful. Behavioral observations were sufficiently detailed to allow examination of how daily mean IR values for males present on the beach varied within seasons and from year to year. Again, significant variation was found both among and within seasons, with more homozygous males appearing less able to hold territories in poor seasons when pup production is low and, within a season, at both the start of the season and to some extent around the peak of female estrus. Finally, we tested whether the benefits of high heterozygosity are due mainly to a genomewide effect (e.g. inbreeding depression) or to single locus heterosis by asking whether the relationship between IR and male success was robust to the removal of any single locus or to any pair of loci. Since the relationship remained significant in all cases, we favor a multilocus explanation for the effects we report.

MALE REPRODUCTIVE STRATEGY AND THE IMPORTANCE OF MATERNAL STATUS IN THE ANTARCTIC FUR SEAL ARCTOCEPHALUS GAZELLA

Abstract. Although mammalian mating systems are classically characterized in terms of male competition and polygyny, it is becoming increasingly apparent that alternative male strategies and female choice may play important roles. For example, females who mate with males from a dominant dynasty risk producing inbred offspring. Many pinnipeds are highly polygynous, but in some species alternative male strategies such as aquatic mating appear to be important, even when behavioral observations suggest strong polygyny. Here, we analyze male reproductive success in the Antarctic fur seal Arctocephalus gazella, an otariid described behaviorally as being highly polygynous, by combining a microsatellite paternity analysis spanning seven consecutive breeding seasons with detailed behavioral data on both sexes. Territorial males fathered 59% of 660 pups analyzed from our study colony. Male reproductive skew was considerable, with a quarter of all paternities assigned to just 12 top individuals on a beach where mean annual pup production was 635. Most males were successful for only a single season, but those able to return over successive years enjoyed rapidly increasing success with each additional season of tenure. We found no evidence of alternative male reproductive tactics such as aquatic or sneaky terrestrial mating. However, paternity was strongly influenced by maternal status. Females observed on the beach without a pup were significantly less likely to conceive to a sampled territorial male than equivalent females that did pup. In addition, their pups carried combinations of paternal alleles that were less likely to be found on the study beach and exhibited lower levels of shared paternity. Thus, from a territorial males perspective, not all females offer equal opportunities for fertilization.

Heterozygosity and lungworm burden in harbour seals (Phoca vitulina)

In several studies, heterozygosity measured at around 10 microsatellite markers correlates with parasite load. Usually the effect size is small, but while this may reflect reality, it may also be possible that too few markers are used or the measure of fitness contains too much error to reveal what is actually a much stronger underlying effect. Here, we analysed over 200 stranded harbour seals (Phoca vitulina) for an association between lungworm burden and heterozygosity, conducting thorough necropsies on the seals and genotyping the samples obtained for 27 microsatellites. We found that homozygosity predicts higher worm burdens, but only in young animals, where the worms have the greatest impact on fitness. Testing each locus separately, we found that a significant majority reveal a weak but similar trend for heterozygosity to be protective against high lungworm burden, suggesting a genome-wide effect, that is, inbreeding. This conclusion is supported by the fact that heterozygosity is correlated among markers in young animals but not in otherwise equivalent older ones. Taken as a whole, our results support the notion that homozygosity increases susceptibility to parasitic infection and suggest that parasites can be effective in removing inbred individuals from the population.

Bayesian inference of a historical bottleneck in a heavily exploited marine mammal

Emerging Bayesian analytical approaches offer increasingly sophisticated means of reconstructing historical population dynamics from genetic data, but have been little applied to scenarios involving demographic bottlenecks. Consequently, we analysed a large mitochondrial and microsatellite dataset from the Antarctic fur seal Arctocephalus gazella, a species subjected to one of the most extreme examples of uncontrolled exploitation in history when it was reduced to the brink of extinction by the sealing industry during the late eighteenth and nineteenth centuries. Classical bottleneck tests, which exploit the fact that rare alleles are rapidly lost during demographic reduction, yielded ambiguous results. In contrast, a strong signal of recent demographic decline was detected using both Bayesian skyline plots and Approximate Bayesian Computation, the latter also allowing derivation of posterior parameter estimates that were remarkably consistent with historical observations. This was achieved using only contemporary samples, further emphasizing the potential of Bayesian approaches to address important problems in conservation and evolutionary biology.