The San Nicolas Island fox (Urocyon littoralis dickeyi) is genetically the most monomorphic sexually reproducing animal population yet reported and has no variation in hypervariable genetic markers. Such low levels of variation imply lower resistance to pathogens, reduced fitness, and problems in distinguishing kin from non-kin. In vertebrates, the MHC contains genes that influence disease resistance and kin recognition and may be under intense balancing selection in some populations. Hence, genetic variation at the MHC might persist despite the extreme monomorphism shown by neutral markers. We examine variation of five loci within the MHC of San Nicolas Island foxes and find remarkably high levels of variation. Further, we show by simulation that genetic monomorphism at neutral loci and high MHC variation could arise only through an extreme population bottleneck of <10 individuals, Ϸ10 -20 generations ago, accompanied by unprecedented selection coefficients of >0.5 on MHC loci. These results support the importance of balancing selection as a mechanism to maintain variation in natural populations and expose the difficulty of using neutral markers as surrogates for variation in fitness-related loci.T he island fox (Urocyon littoralis) is an endemic North American canid that inhabits six of the eight Channel Islands off the coast of southern California (Fig. 1). As suggested by the archeological record and molecular genetic data, foxes colonized the three northern Channel Islands (San Miguel, Santa Rosa, and Santa Cruz) Ϸ16,000 years ago and, subsequently, were transported by Native Americans to the three southern Channel Islands (San Nicolas, Santa Catalina, and San Clemente) 800 to 4,300 years ago (1-4) ( Fig. 1). Effective population size varies with island area and ranges from Ͻ200 to Ϸ1,000 individuals (Table 1). Levels of genetic variation reflect population size and colonization history, with the San Nicolas Island population having the second smallest effective population size and a recent colonization history (4) ( Fig. 1 and Table 1). No variation has been discovered for any of four independent genetic marker classes in the San Nicolas Island population, including supposedly neutral hypervariable microsatellite loci (5) and multilocus fingerprints (2), for which the probability of genetic identity is commonly Ͻ1 in several million (6). Recently, because of dramatic declines, populations on the three Northern Islands and Santa Catalina Island have been proposed for listing under the U.S. Endangered Species Act (7).The MHC contains the most variable set of coding genes in vertebrates, with as many as 349 alleles described for a single locus (8) and heterozygosity values that generally exceed those predicted by neutrality (9). Class I and II MHC molecules are responsible for the presentation to T cells of intracellular (endogenous) and extracellular (exogenous) peptides, respectively (10). High levels of heterozygosity at the MHC may be maintained by balancing selection through pathogen-mediated selection, ...

The major histocompatibility complex (MHC) contains genes integral to immune response in vertebrates. MHC genes have been shown to be under selection in a number of vertebrate taxa, making them intriguing for population genetic studies. We have conducted a survey of genetic variation in an MHC class II gene for steelhead trout from 24 sites in coastal California and compared this variation to that observed at 16 presumably neutral microsatellite loci. A high amount of allelic variation was observed at the MHC when compared to previously published studies on other Pacific salmonids. Elevated nonsynonymous substitutions, relative to synonymous substitutions, were detected at the MHC gene, indicating the signature of historical balancing selection. The MHC data were tested for correlations to and deviations from the patterns found with the microsatellite data. Estimates of allelic richness for the MHC gene and for the microsatellites were positively correlated, as were estimates of population differentiation (F(ST)). An analysis for F(ST) outliers indicates that the MHC locus has an elevated F(ST) relative to the neutral expectation, although a significant result was found for only one particular geographical subgroup. Relatively uniform allele frequency distributions were detected in four populations, although this finding may be partially due to recent population bottlenecks. These results indicate that, at the scale studied here, drift and migration play a major role in the observed geographical variability of MHC genes in steelhead, and that contemporary selection is relatively weak and difficult to detect.

Efforts to taxonomically delineate species are often confounded with conflicting information and subjective interpretation. Advances in genomic methods have resulted in a new approach to taxonomic identification that stands to greatly reduce much of this conflict. This approach is ideal for species complexes, where divergence times are recent (evolutionarily) and lineages less well defined. The California Roach/Hitch fish species complex is an excellent example, experiencing a convoluted geologic history, diverse habitats, conflicting species designations and potential admixture between species. Here we use this fish complex to illustrate how genomics can be used to better clarify and assign taxonomic categories. We performed restriction-site associated DNA (RAD) sequencing on 255 Roach and Hitch samples collected throughout California to discover and genotype thousands of single nucleotide polymorphism (SNPs). Data were then used in hierarchical principal component, admixture, and FST analyses to provide results that consistently resolved a number of ambiguities and provided novel insights across a range of taxonomic levels. At the highest level, our results show that the CA Roach/Hitch complex should be considered five species split into two genera (4 + 1) as opposed to two species from distinct genera (1 +1). Subsequent levels revealed multiple subspecies and distinct population segments within identified species. At the lowest level, our results indicate Roach from a large coastal river are not native but instead introduced from a nearby river. Overall, this study provides a clear demonstration of the power of genomic methods for informing taxonomy and serves as a model for future studies wishing to decipher difficult species questions. By allowing for systematic identification across multiple scales, taxonomic structure can then be tied to historical and contemporary ecological, geographic or anthropogenic factors.

Species endemic to alpine environments can evolve via steep ecological selection gradients between lowland and upland environments. Additionally, many alpine environments have faced repeated glacial episodes over the past two million years, fracturing these endemics into isolated populations. In this “glacial pulse” model of alpine diversification, cycles of allopatry and ecologically divergent glacial refugia play a role in generating biodiversity, including novel admixed (“fused”) lineages. We tested for patterns of glacial pulse lineage diversification in the Yosemite toad (Anaxyrus [Bufo] canorus), an alpine endemic tied to glacially influenced meadow environments. Using double‐digest RADseq on populations densely sampled from a portion of the species range, we identified nine distinct lineages with divergence times ranging from 18 to 724 thousand years ago (ka), coinciding with multiple Sierra Nevada glacial events. Three lineages have admixed origins, and demographic models suggest these fused lineages have persisted throughout past glacial cycles. Directionality indices supported the hypothesis that some lineages recolonized Yosemite from east of the ice sheet, whereas other lineages remained in western refugia. Finally, refugial niche reconstructions suggest that low‐ and high‐elevation lineages have convergently adapted to similar climatic niches. Our results suggest glacial cycles and refugia may be important crucibles of adaptive diversity across deep evolutionary time.

No abstract

The amount and distribution of population genetic variation is crucial information for the design of effective conservation strategies for endangered species and can also be used to provide inference about demographic processes and patterns of migration. Here, we describe variation at a large number of nuclear genes in sea otters Enhydra lutris ssp. We surveyed 14 variable microsatellite loci and two genes of the major histocompatibility complex (MHC) in up to 350 California sea otters Enhydra lutris nereis, which represents $10% of the subspecies' population, and 46 otters from two Alaskan sites. We utilized methods for detecting past reductions in effective population size to examine the effects of near extinction from the fur trade. Summary statistic tests largely failed to find a signal of a recent population size reduction (within the past 200 years), but a Bayesian method found a signal of a strong reduction over a longer time scale (up to 500 years ago). These results indicate that the reduction in size began long enough ago that much genetic variation was lost before the 19th century fur trade. A comparison of geographic distance and pairwise relatedness for individual otters found no evidence of kinbased spatial clustering for either gender. This indicates that there is no population structure, due to extended family groups, within the California population. A survey of population genetic variation found that two of the MHC genes, DQB and DRB, had two alleles present and one of the genes, DRA, was monomorphic in otters. This contrasts with other mammals, where they are often the most variable coding genes known. Genetic variation in the sea otter is among the lowest observed for a mammal and raises concerns about the long-term viability of the species, particularly in the face of future environmental changes.

We have isolated major histocompatibility complex (MHC) class II beta loci from the little greenbul (Andropadus virens), an African songbird. We utilized preexisting information about conserved regions of the avian MHC to design primers to amplify a pool of sequences representing multiple loci. From this pool, a unique locus spanning 1109 bp that we designate as Anvi-DAB1 was cloned and sequenced. We designed locus-specific primers based on this sequence information and amplified six alleles from seven individuals. Compared to other A. virens MHC sequences obtained from genomic DNA or cDNA, the variability of sequences from Anvi-DAB1 was low and the ratio of nonsynonymous to synonymous substitution was much less than one, suggesting that Anvi-DAB1 may either be a pseudogene or a nonclassical MHC locus. Phylogenetic analysis revealed that the Anvi-DAB1 locus was highly divergent when compared with other passerine or A. virens genomic or transcribed MHC sequences. The use of conserved MHC primers followed by analysis of cloned sequences allows rapid isolation of MHC loci from exotic species and avoids laborious large-scale cloning and sequencing.

We used molecular genetic techniques to identify the source population of northern pike (Esox lucius) illegally introduced into Lake Davis, California, and to study the dynamics of genetic change since the introduction. We typed 10 tetranucleotide microsatellite loci from 11 populations and used measures of population differentiation, phylogenetic analysis, and assignment tests to determine the origins of the pre- and post-rotenone-treatment populations. We found that levels of genetic differentiation were low between the pre- and both post-treatment Lake Davis samples and phylogenetic analysis revealed that they were all closely related. Assignment tests classified a high proportion of posttreatment individuals to the original Lake Davis population. Consequently, our results indicate that the fish population that reappeared in 1999 and a subsequent sample in 2002 were descendents or survivors of the initial population. We were unable to assign the pre- or post-treatment individuals consistently to any one population from a panel of potential source populations, indicating that the ultimate source of the Lake Davis population was likely not sampled. Genetic signatures of a recent population bottleneck are evident in the samples from Lake Davis and the effective population size over the sampling period was approximately 7–12 individuals.