Temporal dynamics of linkage disequilibrium in two populations of bighorn sheep

Linkage disequilibrium (LD) is the nonrandom association of alleles at two markers. Patterns of LD have biological implications as well as practical ones when designing association studies or conservation programs aimed at identifying the genetic basis of fitness differences within and among populations. However, the temporal dynamics of LD in wild populations has received little empirical attention. In this study, we examined the overall extent of LD, the effect of sample size on the accuracy and precision of LD estimates, and the temporal dynamics of LD in two populations of bighorn sheep (Ovis canadensis) with different demographic histories. Using over 200 microsatellite loci, we assessed two metrics of multi‐allelic LD, D′, and χ^′2. We found that both populations exhibited high levels of LD, although the extent was much shorter in a native population than one that was founded via translocation, experienced a prolonged bottleneck post founding, followed by recent admixture. In addition, we observed significant variation in LD in relation to the sample size used, with small sample sizes leading to depressed estimates of the extent of LD but inflated estimates of background levels of LD. In contrast, there was not much variation in LD among yearly cross‐sections within either population once sample size was accounted for. Lack of pronounced interannual variability suggests that researchers may not have to worry about interannual variation when estimating LD in a population and can instead focus on obtaining the largest sample size possible.     

Testing marker-based estimates of heritability in the wild

Marker-based estimates of heritability are an attractive alternative to pedigree-based methods for estimating quantitative genetic parameters in field studies where it is difficult or impossible to determine relationships and pedigrees. Here I test the ability of the marker-based method to estimate heritability of a suite of traits in a wild population of bighorn sheep (Ovis canadensis) using marker data from 32 microsatellite loci. I compared marker-based estimates with estimates obtained using a pedigree and the animal model. Marker-based estimates of heritability were imprecise and downwardly biased. The high degree of uncertainty in marker-based estimates suggests that the method may be sufficient to detect the presence of genetic variance for highly heritable traits, but not sufficiently reliable to estimate genetic parameters.     

A quantitative trait locus analysis of personality in wild bighorn sheep

Personality, the presence of persistent behav105 ioral differences among individuals over time or contexts, potentially has important ecological and evolutionary consequences. However, a lack of knowledge about its genetic architecture limits our ability to understand its origin, evolution, and maintenance. Here, we report on a genome‐wide quantitative trait locus (QTL) analysis for two personality traits, docility and boldness, in free‐living female bighorn sheep from Ram Mountain, Alberta, Canada. Our variance component linkage analysis based on 238 microsatellite loci genotyped in 310 pedigreed individuals identified suggestive docility and boldness QTL on sheep chromosome 2 and 6, respectively. A lack of QTL overlap indicated that genetic covariance between traits was not modulated by pleiotropic effects at a major locus and may instead result from linkage disequilibrium or pleiotropic effects at QTL of small effects. To our knowledge, this study represents the first attempt to dissect the genetic architecture of personality in a free‐living wildlife population, an important step toward understanding the link between molecular genetic variation in personality and fitness and the evolutionary processes maintaining this variation.     

Sexing the Sciuridae: a simple and accurate set of molecular methods to determine sex in tree squirrels, ground squirrels and marmots

We determined the sequence of the male-specific minor histocompatibility complex antigen (Smcy) from the Y chromosome of seven squirrel species (Sciuridae, Rodentia). Based on conserved regions inside the Smcy intron sequence, we designed PCR primers for sex determination in these species that can be co-amplified with nuclear loci as controls. PCR co-amplification yields two products for males and one for females that are easily visualized as bands by agarose gel electrophoresis. Our method provides simple and reliable sex determination across a wide range of squirrel species.     

Comparing measures of breeding inequality and opportunity for selection with sexual selection on a quantitative character in bighorn rams

The reliability and consistency of the many measures proposed to quantify sexual selection have been questioned for decades. Realized selection on quantitative characters measured by the selection differential i was approximated by metrics based on variance in breeding success, using either the opportunity for sexual selection I_s or indices of inequality. There is no consensus about which metric best approximates realized selection on sexual characters. Recently, the opportunity for selection on character mean OSM was proposed to quantify the maximum potential selection on characters. Using 21 years of data on bighorn sheep (Ovis canadensis), we investigated the correlations between seven indices of inequality, I_s, OSM and i on horn length of males. Bighorn sheep are ideal for this comparison because they are highly polygynous and sexually dimorphic, ram horn length is under strong sexual selection, and we have detailed knowledge of individual breeding success. Different metrics provided conflicting information, potentially leading to spurious conclusions about selection patterns. I_δ, an index of breeding inequality, and, to a lesser extent, I_s showed the highest correlation with i on horn length, suggesting that these indices document breeding inequality in a selection context. OSM on horn length was strongly correlated with i, I_s and indices of inequality. By integrating information on both realized sexual selection and breeding inequality, OSM appeared to be the best proxy of sexual selection and may be best suited to explore its ecological bases.     

Assessment of identity disequilibrium and its relation to empirical heterozygosity fitness correlations: a meta‐analysis

Heterozygosity fitness correlations (HFCs) have frequently been used to detect inbreeding depression, under the assumption that genome‐wide heterozygosity is a good proxy for inbreeding. However, meta‐analyses of the association between fitness measures and individual heterozygosity have shown that often either no correlations are observed or the effect sizes are small. One of the reasons for this may be the absence of variance in inbreeding, a requisite for generating general‐effect HFCs. Recent work has highlighted identity disequilibrium (ID) as a measure that may capture variance in the level of inbreeding within a population; however, no thorough assessment of ID in natural populations has been conducted. In this meta‐analysis, we assess the magnitude of ID (as measured by the g₂ statistic) from 50 previously published HFC studies and its relationship to the observed effect sizes of those studies. We then assess how much power the studies had to detect general‐effect HFCs, and the number of markers that would have been needed to generate a high expected correlation (r² = 0.9) between observed heterozygosity and inbreeding. Across the majority of studies, g₂ values were not significantly different than zero. Despite this, we found that the magnitude of g₂ was associated with the average effect sizes observed in a population, even when point estimates were nonsignificant. These low values of g₂ translated into low expected correlations between heterozygosity and inbreeding and suggest that many more markers than typically used are needed to robustly detect HFCs.     

HOT SPOTS OF GENETIC DIVERSITY DESCENDED FROM MULTIPLE PLEISTOCENE REFUGIA IN AN ALPINE UNGULATE

Species that inhabit naturally fragmented environments are expected to be spatially structured and exhibit reduced genetic diversity at the periphery of their range. Patterns of differentiation may also reflect historical processes such as recolonization from glacial refugia. We examined the relative importance of these factors in shaping the spatial patterns of genetic differentiation across the range of an alpine specialist, the North American mountain goat (Oreamnos americanus). Contrary to fossil evidence that suggests a single southern refugium, we detected evidence for additional refugia in northern British Columbia and the Alaskan coast using both mitochondrial and microsatellite DNA. A core area of elevated genetic diversity characterized both regions, and molecular dating suggested a recent Pleistocene split was followed by demographic expansion. Across their range, mountain goats were highly genetically structured and displayed the expected pattern of declining diversity toward the periphery. Gene flow was high within contiguous mountain ranges, but cross‐assignments paradoxically suggest that long‐distance contemporary dispersal movements are not uncommon. These results improve our understanding of how historical vicariance and contemporary fragmentation influence population differentiation, and have implications for conserving the adaptive potential of alpine populations and habitat.     

Molecular phylogeny of North American Branchiobdellida (Annelida: Clitellata)

Branchiobdellidans, or crayfish worms, are ectosymbiotic clitellate annelids associated primarily with freshwater crayfishes. The main objectives of our study were to infer a molecular phylogeny for the North American Branchiobdellida, examine its congruence with morphology-based hypotheses of relationships at the subfamily and genus level, and use our dataset to assess consistency of GenBank-archived branchiobdellidan sequences. We used nucleotide sequence data from two mtDNA genes (COI and 16S rDNA) and three nuclear genes (28S rDNA, 18S rDNA, and ITS1) to estimate phylogenetic relationships among 47 described and one undescribed species of Branchiobdellida. We recovered a monophyletic branchiobdellidan clade with generally short branch lengths, suggesting that a large portion of the taxon has likely undergone a recent and rapid radiation in North America. Results from our phylogenetic analyses indicate that current taxonomic groupings are largely unsupported by the molecular data. All four subfamilies are either paraphyletic or polyphyletic, and only three of seven sampled non-monotypic genera were monophyletic. We found a high rate (49%) of inconsistency in GenBank-archived sequences, over 70% of which can be attributed to field- or laboratory-based error.