Conservation genetics of the fisher (Martes pennanti) based on mitochondrial DNA sequencing

Translocation of animals to re-establish extirpated populations or to maintain declining ones has often been carried out without genetic information on source or target populations, or adequate consideration of the potential effects of mixing genetic stocks. We consider the conservation status of the fisher (Martes pennanti) and evaluate the potential genetic consequences of past and future translocations on this medium-sized carnivore by examining population variation in mitochondrial control-region sequences. We sampled populations throughout the fisher's range in North America including five populations unaffected by translocations and two western populations that had received long-distance translocations. Twelve haplotypes showed little sequence divergence. Haplotype frequencies differed significantly among subspecies and between populations within subspecies. Analysis of molecular variance (amova) and neighbour-joining analyses of haplotype relationships revealed population subdivision similar to current subspecies designations, but which may reflect an isolation-by-distance pattern. Populations in Oregon and in Montana and Idaho received several translocations and each showed greater similarity to the populations where translocations originated than to adjacent populations. Additional sequences obtained from museum specimens collected prior to any translocations suggest historical gene flow among populations in British Columbia, Washington, Oregon, and California. Anthropogenic impacts in that region have greatly reduced and isolated extant populations in Oregon and California. Future translocations may be necessary to recover populations in Washington and portions of Oregon and California; our results indicate that British Columbia would be the most appropriate source population.     

Genetic diversity, structure, and demographic change in tanoak, Lithocarpus densiflorus (Fagaceae), the most susceptible species to sudden oak death in California

Knowledge of population genetic structure of tanoak (Lithocarpus densiflorus) is of interest to pathologists seeking natural variation in resistance to sudden oak death disease, to resource managers who need indications of conservation priorities in this species now threatened by the introduced pathogen (Phytophthora ramorum), and to biologists with interests in demographic processes that have shaped plant populations. We investigated population genetic structure using nuclear and chloroplast DNA (cpDNA) and inferred the effects of past population demographic processes and contemporary gene flow. Our cpDNA results revealed a strong pattern of differentiation of four regional groups (coastal California, southern Oregon, Klamath mountains, and Sierra Nevada). The chloroplast haplotype phylogeny suggests relatively deep divergence of Sierra Nevada and Klamath populations from those of coastal California and southern Oregon. A widespread coastal California haplotype may have resulted from multiple refugial sites during the Last Glacial Maximum or from rapid recolonization from few refugia. Analysis of nuclear microsatellites suggests two major groups: (1) central coastal California and (2) Sierra Nevada/Klamath/southern Oregon and an area of admixture in north coastal California. The low level of nuclear differentiation is likely to be due to pollen gene flow among populations during postglacial range expansion.     

Phylogeography and spatial genetic structure of the Southern torrent salamander: implications for conservation and management

The Southern torrent salamander (Rhyacotriton variegatus) was recently found not warranted for listing under the US Endangered Species Act due to lack of information regarding population fragmentation and gene flow. Found in small-order streams associated with late-successional coniferous forests of the US Pacific Northwest, threats to their persistence include disturbance related to timber harvest activities. We conducted a study of genetic diversity throughout this species' range to 1) identify major phylogenetic lineages and phylogeographic barriers and 2) elucidate regional patterns of population genetic and spatial phylogeographic structure. Cytochrome b sequence variation was examined for 189 individuals from 72 localities. We identified 3 major lineages corresponding to nonoverlapping geographic regions: a northern California clade, a central Oregon clade, and a northern Oregon clade. The Yaquina River may be a phylogeographic barrier between the northern Oregon and central Oregon clades, whereas the Smith River in northern California appears to correspond to the discontinuity between the central Oregon and northern California clades. Spatial analyses of genetic variation within regions encompassing major clades indicated that the extent of genetic structure is comparable among regions. We discuss our results in the context of conservation efforts for Southern torrent salamanders.     

Characterizing the molecular variation among American marten (Martes americana) subspecies from Oregon and California

We investigated the subspecific identity of a rediscovered population of American martens within the range of a presumed extinct subspecies (Martes americana humboldtensis) by comparing mitochondrial DNA sequence diversity from contemporary individuals within the described ranges of M. a. humboldtensis, nearby ranges of M. a. caurina and M. a. sierrae, and a museum specimen of M. a. humboldtensis. Martens from the rediscovered population shared a haplotype (#2) with the museum specimen. This haplotype was found only in the coastal regions of Oregon and California, suggesting that the rediscovered population represents descendants of a relictual population that previously existed in coastal California. The subspecific boundary between M. a. humboldtensis and M. a. caurina may not be valid, because haplotype #2 was shared between coastal Oregon and coastal California populations and no known contemporary or historical biogeographic barriers prevent north–south movement. Thus, marten populations currently located in coastal forests of California and Oregon should be managed collectively to preserve the connectivity that our data suggest occurred historically. M. a. sierrae differed substantially from both M. a. humboldtensis and M. a. caurina, suggesting marten populations were not a historically genetically homogeneous population and divergence may have occurred in separate glacial refugia.     

Seascape drivers of Macrocystis pyrifera population genetic structure in the northeast Pacific

At small spatial and temporal scales, genetic differentiation is largely controlled by constraints on gene flow, while genetic diversity across a species' distribution is shaped on longer temporal and spatial scales. We assess the hypothesis that oceanographic transport and other seascape features explain different scales of genetic structure of giant kelp, Macrocystis pyrifera. We followed a hierarchical approach to perform a microsatellite‐based analysis of genetic differentiation in Macrocystis across its distribution in the northeast Pacific. We used seascape genetic approaches to identify large‐scale biogeographic population clusters and investigate whether they could be explained by oceanographic transport and other environmental drivers. We then modelled population genetic differentiation within clusters as a function of oceanographic transport and other environmental factors. Five geographic clusters were identified: Alaska/Canada, central California, continental Santa Barbara, California Channel Islands and mainland southern California/Baja California peninsula. The strongest break occurred between central and southern California, with mainland Santa Barbara sites forming a transition zone between the two. Breaks between clusters corresponded approximately to previously identified biogeographic breaks, but were not solely explained by oceanographic transport. An isolation‐by‐environment (IBE) pattern was observed where the northern and southern Channel Islands clustered together, but not with closer mainland sites, despite the greater distance between them. The strongest environmental association with this IBE pattern was observed with light extinction coefficient, which extends suitable habitat to deeper areas. Within clusters, we found support for previous results showing that oceanographic connectivity plays an important role in the population genetic structure of Macrocystis in the Northern hemisphere.     

Evidence of spatial genetic structure in a California bunchgrass population

We investigated the scale of genetic variation of purple needlegrass (Nassella pulchra), a species commonly used in California for grassland restoration. Common garden and field data revealed evidence of genetic differentiation between two intermixed microhabitats characterized by differences in soil depth and community composition. We assessed the genetic variation within a single population using randomly amplified polymorphic DNA (RAPD) data collected from clusters of five individuals in 40 locations. We found no evidence for genetic structure at the whole population level. At smaller spatial scales, however, we found strong evidence that genetic subdivision of the population occurs at the level of the maternal neighborhood. We suggest that the interaction between widespread pollen dispersal and restricted seed dispersal may be the primary factor generating these results; panmictic pollen dispersal will make detection of genetic patterning difficult at larger spatial scales while limited seed dispersal will generate local genetic structure. As a result, the detection of population genetic structure will depend on the spatial scale of analysis. Local selection gradients related to topography and soil depth are also likely to play a role in structuring local genetic variation. Since N. pulchra is widely used in California in grassland and woodland habitat restoration, we suggest that, as a general rule, care should be exercised in transferring germplasm for the purposes of conservation when little is known about the within-population genetic subdivision of a plant species. Received: 23 December 1996 / Accepted: 20 May 1997     

Major biogeographic barriers in eastern Australia have shaped the population structure of widely distributed Eucalyptus moluccana and its putative subspecies

We have investigated the impact of recognized biogeographic barriers on genetic differentiation of grey box (Eucalyptus moluccana), a common and widespread tree species of the family Myrtaceae in eastern Australian woodlands, and its previously proposed four subspecies moluccana, pedicellata, queenslandica, and crassifolia. A range of phylogeographic analyses were conducted to examine the population genetic differentiation and subspecies genetic structure in E. moluccana in relation to biogeographic barriers. Slow evolving markers uncovering long term processes (chloroplast DNA) were used to generate a haplotype network and infer phylogeographic barriers. Additionally, fast evolving, hypervariable markers (microsatellites) were used to estimate demographic processes and genetic structure among five geographic regions (29 populations) across the entire distribution of E. moluccana. Morphological features of seedlings, such as leaf and stem traits, were assessed to evaluate population clusters and test differentiation of the putative subspecies. Haplotype network analysis revealed twenty chloroplast haplotypes with a main haplotype in a central position shared by individuals belonging to the regions containing the four putative subspecies. Microsatellite analysis detected the genetic structure between Queensland (QLD) and New South Wales (NSW) populations, consistent with the McPherson Range barrier, an east‐west spur of the Great Dividing Range. The substructure was detected within QLD and NSW in line with other barriers in eastern Australia. The morphological analyses supported differentiation between QLD and NSW populations, with no difference within QLD, yet some differentiation within NSW populations. Our molecular and morphological analyses provide evidence that several geographic barriers in eastern Australia, including the Burdekin Gap and the McPherson Range have contributed to the genetic structure of E. moluccana. Genetic differentiation among E. moluccana populations supports the recognition of some but not all the four previously proposed subspecies, with crassifolia being the most differentiated.     

DNA from bird-dispersed seed and wind-disseminated pollen provides insights into postglacial colonization and population genetic structure of whitebark pine (Pinus albicaulis)

Uniparentally inherited mitochondrial (mt)DNA and chloroplast (cp)DNA microsatellites (cpSSRs) were used to examine population genetic structure and biogeographic patterns of bird-dispersed seed and wind-disseminated pollen of whitebark pine (Pinus albicaulis Engelm.). Sampling was conducted from 41 populations throughout the range of the species. Analyses provide evidence for an ancestral haplotype and two derived mtDNA haplotypes with distinct regional distributions. An abrupt contact zone between mtDNA haplotypes in the Cascade Range suggests postglacial biogeographic movements. Among three cpSSR loci, 42 haplotypes were detected within 28 cpSSR sample populations that were aggregated into six regions. Analysis of molecular variance (amova) was used to determine the hierarchical genetic structure of cpSSRs. amova and population pairwise comparisons (FST ) of cpSSR, and geographical distribution of mtDNA haplotypes provide insights into historical changes in biogeography. The genetic data suggest that whitebark pine has been intimately tied to climatic change and associated glaciation, which has led to range movements facilitated by seed dispersal by Clark's nutcracker (Nucifraga columbiana Wilson). The two hypotheses proposed to explain the genetic structure are: (i) a northward expansion into Canada and the northern Cascades in the early Holocene; and (ii) historical gene flow between Idaho and the Oregon Cascades when more continuous habitat existed in Central Oregon during the late Pleistocene. Genetic structure and insights gained from historical seed movements provide a basis on which to develop recovery plans for a species that is at risk from multiple threats.     

Microsatellite variation and fine-scale population structure in the wood frog (Rana sylvatica)

We investigated genetic population structure in wood frogs (Rana sylvatica) from a series of Prairie Pothole wetlands in the northern Great Plains. Amphibians are often thought to exist in demographic metapopulations, which require some movement between populations, yet genetic studies have revealed strong subdivision among populations, even at relatively fine scales (several km). Wood frogs are highly philopatric and studies of dispersal suggest that they may exhibit subdivision on a scale of approximately 1-2 km. We used microsatellites to examine population structure among 11 breeding assemblages separated by as little as 50 m up to approximately 5.5 km, plus one population separated from the others by 20 km. We found evidence for differentiation at the largest distances we examined and among a few neighbouring ponds, but most populations were strikingly similar in allele frequencies, suggesting high gene flow among all but the most distant populations. We hypothesize that the few significant differences among neighbouring populations at the finest scale may be a transient effect of extinction-recolonization founder events, driven by periodic drying of wetlands in this hydrologically dynamic landscape.     

海南岛海菖蒲种群克隆多样性和遗传结构

海草是生长在潮间带和潮下带的单子叶植物,由海草植物组成的海草床是生态系统服务价值最高的生态系统之一.然而,近几十年人类活动干扰、全球气候变化等因素导致海草床衰退严重.海菖蒲是分布于热带、体型最大的雌雄异株海草,我国位于该物种的分布北缘,本文对其克隆多样性和遗传结构进行研究,以期为该海草的保护提供参考.采用4对多态微卫星标记对采自海南岛4个地点的现存海菖蒲种群的样品进行基因型分型.结果表明:海菖蒲种群克隆多样性和遗传多样性较低,这与所研究种群处于分布区北缘有关;种群间遗传分化值范围较大(0.073~0.309),这可能是由于分布于不同港湾的种群间距离范围较大以及局域绝灭/再拓殖的遗传漂变效应所致;各种群未发现近期经历种群瓶颈的信号,很可能是由于种群内遗传多样性已经很低,种群减小未能导致遗传多样性明显降低.根据种群遗传特征,提出了重点保护种群的建议,鉴于目前我国海菖蒲等海草快速衰退的局面,应强化海草保护并实施海草床生态恢复.     

Genetic diversity and biogeography of red turpentine beetle Dendroctonus valens in its native and invasive regions

Sequences of 479 bp region of the mitochondrial COI gene were applied to detect population genetic diversity and structure of Dendroctonus valens populations. By comparing the genetic diversity between native and invasive populations, it was shown that the genetic diversity of Chinese populations was obviously lower than that of native populations with both indices of haplotype diversity and Nei's genetic diversity, suggesting genetic bottleneck occurred in the invasive process of D. valens , and was then followed by a relatively quick population buildup. According to phylogenetic analyses of haplotypes, we suggested that the origin of the Chinese population was from California, USA. Phylogenetic and network analysis of native populations of D. valens revealed strong genetic structure at two distinct spatial and temporal scales in North America. The main cause resulting in current biogeographic pattern was supposedly due to recycled glacial events. Meanwhile, a cryptic species might exist in the Mexican and Guatemalan populations.     

Chloroplast DNA variation and population structure in the widespread forest tree, Eucalyptus grandis

Recognition of genetic structure of populations and the ability to identify vulnerable populations is useful for the formation of conservation management strategies for plants. Eucalyptus grandis is a tall forest tree that has a major area of occurrence in subtropical eastern Australia, with smaller populations located in the east coast tropics. Many widespread forest species exhibit population differentiation that corresponds to geographic regions. However, Eucalyptus grandis appears to be an exception based on isozyme and morphological data. This is intriguing given a large discontinuity between northern populations and those in the southern part of the species range. In this study, the distribution of a maternally inherited chloroplast locus was examined because it was more likely to reveal genetic structure due to the slower evolution of the chloroplast genome and limited dispersal of seed in eucalypts. As expected, the G _ST for chloroplast DNA was higher than that for nuclear DNA but indicated low population differentiation for a forest tree species. Phylogeographic analysis indicated that the 15 populations grouped into three broad geographical regions; however, overall population structure was weak suggesting that the large geographical disjunction in the distribution of E. grandis may be relatively recent. A paradigm for conservation management of E. grandis based on chloroplast DNA haplotype distribution would take into account the low differentiation among populations.     

Environmental heterogeneity and not vicariant biogeographic barriers generate community‐wide population structure in desert‐adapted snakes

Genetic structure can be influenced by local adaptation to environmental heterogeneity and biogeographic barriers, resulting in discrete population clusters. Geographic distance among populations, however, can result in continuous clines of genetic divergence that appear as structured populations. Here, we evaluate the relevant importance of these three factors over a landscape characterized by environmental heterogeneity and the presence of a hypothesized biogeographic barrier in producing population genetic structure within 13 codistributed snake species using a genomic data set. We demonstrate that geographic distance and environmental heterogeneity across western North America contribute to population genomic divergence. Surprisingly, landscape features long thought to contribute to biogeographic barriers play little role in divergence community wide. Our results suggest that isolation by environment is the most important contributor to genomic divergence. Furthermore, we show that models of population clustering that incorporate spatial information consistently outperform nonspatial models, demonstrating the importance of considering geographic distances in population clustering. We argue that environmental and geographic distances as drivers of community‐wide divergence should be explored before assuming the role of biogeographic barriers.     

Population structure and genetic diversity in the endangered bluemask darter (Etheostoma akatulo)

Bluemask darters (Etheostoma akatulo) were sampled from the four drainages where extant populations of this narrowly endemic freshwater fish are known to exist. Population genetic diversity and structure were assessed at 10 microsatellite loci. All populations exhibited low levels of genetic variation, with expected heterozygosity ranging from 0.2 to 0.35. Significant population subdivision was found among most tributaries, and genetic divergence was strongly correlated with geographic distance. Bayesian population assignment and pairwise population differentiation measures both identified a lack of differentiation between E. akatulo populations inhabiting Cane Creek and the Caney Fork. This observation reduced the number of distinct breeding populations of this species to three. We also used approximate Bayesian computation to compare three models of demographic history in this species. A constant population size model was favored over models that included historic or recent population reductions. Our results suggest that impoundment of the Caney Fork and its tributaries, by completion of Great Falls Dam in 1916, was not responsible for the reduced genetic diversity in the sampled populations. Given the low levels of genetic diversity within populations and the limited geographic distribution, future conservation efforts should seek to maximize available habitat while simultaneously limiting the influences of anthropogenic stressors in the system.     

Genetic structure of California's elk: a legacy of extirpations,reintroductions, population expansions,and admixture

California, USA, is home to 3 subspecies of North American elk (Cervus canadensis): Roosevelt (C. c. roosevelti), Rocky Mountain (C. c. nelsoni), and tule (C. c. nannodes). Effective management requires a baseline understanding of each subspecies' range, admixture zones, and geographic patterns of genetic diversity. To address these questions, we genotyped 1,271 individual elk from California (n = 1,204) and reference populations of Rocky Mountain and Roosevelt elk from Nevada (n = 32) and Oregon (n = 35), USA. Using 19 polymorphic microsatellite loci, we detected admixture between Roosevelt and Rocky Mountain elk at a contact zone in northern California, and between Roosevelt and tule elk in north-coastal California and central-coastal California. We identified a genetically distinct population of Roosevelt elk in northwestern California, likely reflecting the remnant population that survived a large demographic decline from overhunting during the 1800s. Tule elk exhibited lower levels of heterozygosity (0.44 ± 0.03 [SD]) and allelic richness (2.9 ± 0.2) than Rocky Mountain (0.58 ± 0.05, 4.9 ± 0.4, respectively) and Roosevelt (0.50 ± 0.06, 4.4 ± 0.6, respectively) elk. Among tule elk populations, heterozygosity varied, with the lowest heterozygosity (0.23 ± 0.05) corresponding to the oldest enclosed herd used over the past century as a source of translocations. Among tule elk populations, genetic structure revealed several cases of successful and unsuccessful reintroduction or augmentation attempts. Results provide an essential baseline for future monitoring and decisions about harvest management, translocations to preserve genetic diversity, and landscape-level conservation planning to maintain, enhance, or obstruct connectivity of elk populations. Genome-wide sequencing and analyses are needed to quantify inbreeding absolutely and assess genetic load and the age of admixture where subspecies currently exchange genes.     

Historical and contemporary population genetic connectivity of the European short-snouted seahorse Hippocampus hippocampus and implications for management

This first genetic study of Hippocampus hippocampus covers the species' entire geographic range and employs two mtDNA markers (control region and cytochrome b) to establish patterns of population structuring. A total of 255 specimens from 21 locations were used to obtain 89 concatenated haplotypes. The common haplotype was present in all but one population, however, most haplotypes were unique. The haplotype network had a star-like construction, suggesting expansion from a bottleneck event. F(ST) and AMOVA revealed population subdivision into three geographic regions (English Channel + Bay of Biscay, Mediterranean Sea + Atlantic Ocean Iberian coast + Macaronesian Islands, and West Africa) with barriers to gene flow indentified at Cape Finisterre and the Cape Verde frontal zone. Neutrality tests and nested clade analysis suggest a complex demographic history, with both historic events and contemporary processes shaping patterns of genetic differentiation. The genetic population subdivision detected in this study indicates that H. hippocampus should be managed as three separate units. This is especially pertinent as H. hippocampus populations within the West African region are the only ones known to be specifically targeted for exploitation.     

Genetic structure of desert ground squirrels over a 20-degree-latitude transect from Oregon through the Baja California peninsula

The genetic structure of populations over a wide geographical area should reflect the demographic and evolutionary processes that have shaped a species across its range. We examined the population genetic structure of antelope ground squirrels (Ammospermophilus leucurus) across the complex of North American deserts from the Great Basin of Oregon to the cape region of the Baja California peninsula. We sampled 73 individuals from 13 major localities over this 2500-km transect, from 43 to 22 degrees north. Our molecular phylogeographical analysis of 555 bp of the mitochondrial cytochrome b gene and 510 bp of the control region revealed great genetic uniformity in a single clade that extends from Oregon to central Baja California. A second distinct clade occupies the southern half of the peninsula. The minimal geographical structure of the northern clade, its low haplotype diversity and the distribution of pairwise differences between haplotypes suggest a rapid northward expansion of the population that must have followed a northward desert habitat shift associated with the most recent Quaternary climate warming and glacial retreat. The higher haplotype diversity within the southern clade and distribution of pairwise differences between haplotypes suggest that the southern clade has a longer, more stable history associated with a southern peninsular refugium. This system, as observed, reflects both historical and contemporary ecological and evolutionary responses to physical environmental gradients within genetically homogeneous populations.     

Genotyping‐in‐Thousands by sequencing reveals marked population structure in Western Rattlesnakes to inform conservation status

Delineation of units below the species level is critical for prioritizing conservation actions for species at‐risk. Genetic studies play an important role in characterizing patterns of population connectivity and diversity to inform the designation of conservation units, especially for populations that are geographically isolated. The northernmost range margin of Western Rattlesnakes (Crotalus oreganus) occurs in British Columbia, Canada, where it is federally classified as threatened and restricted to five geographic regions. In these areas, Western Rattlesnakes hibernate (den) communally, raising questions about connectivity within and between den complexes. At present, Western Rattlesnake conservation efforts are hindered by a complete lack of information on genetic structure and degree of isolation at multiple scales, from the den to the regional level. To fill this knowledge gap, we used Genotyping‐in‐Thousands by sequencing (GT‐seq) to genotype an optimized panel of 362 single nucleotide polymorphisms (SNPs) from individual samples (n = 461) collected across the snake's distribution in western Canada and neighboring Washington (USA). Hierarchical STRUCTURE analyses found evidence for population structure within and among the five geographic regions in BC, as well as in Washington. Within these regions, 11 genetically distinct complexes of dens were identified, with some regions having multiple complexes. No significant pattern of isolation‐by‐distance and generally low levels of migration were detected among den complexes across regions. Additionally, snakes within dens generally were more related than those among den complexes within a region, indicating limited movement. Overall, our results suggest that the single, recognized designatable unit for Western Rattlesnakes in Canada should be re‐assessed to proactively focus conservation efforts on preserving total genetic variation detected range‐wide. More broadly, our study demonstrates a novel application of GT‐seq for investigating patterns of diversity in wild populations at multiple scales to better inform conservation management.     

Fine scale population genetic structure of pumas in the Intermountain West

In this study, I examined the population genetic structure of subpopulations of pumas (Puma concolor) in Idaho and surrounding states. Patterns of genetic diversity, population structure, levels of inbreeding, and the relationship between genetic differentiation and dispersal distance within and between 15 subpopulations of pumas were compared. Spatial analyses revealed that the Snake River plain was an important barrier to movement between northern and southern regions of Idaho. In addition, subpopulations south of the Snake River plain exhibited lower levels of genetic diversity, higher levels of inbreeding, and a stronger pattern of isolation by distance relative to subpopulations north of the Snake River plain. Lower levels of diversity and restricted gene flow are likely the result of historically lower population sizes in conjunction with more recent changes in habitat use and available dispersal corridors for movement. The subdivision of puma populations north and south of the Snake River plain, along with the patterns of genetic diversity within regions, indicate that landscape features are affecting the population genetic structure of pumas in Idaho. These results indicate that information about the effects of landscape features on the distribution of genetic diversity should be considered when designing plans for the management and conservation of pumas.     

Genetic diversity and population parameters of sea otters, Enhydra lutris, before fur trade extirpation from 1741-1911

All existing sea otter, Enhydra lutris, populations have suffered at least one historic population bottleneck stemming from the fur trade extirpations of the eighteenth and nineteenth centuries. We examined genetic variation, gene flow, and population structure at five microsatellite loci in samples from five pre-fur trade populations throughout the sea otter's historical range: California, Oregon, Washington, Alaska, and Russia. We then compared those values to genetic diversity and population structure found within five modern sea otter populations throughout their current range: California, Prince William Sound, Amchitka Island, Southeast Alaska and Washington. We found twice the genetic diversity in the pre-fur trade populations when compared to modern sea otters, a level of diversity that was similar to levels that are found in other mammal populations that have not experienced population bottlenecks. Even with the significant loss in genetic diversity modern sea otters have retained historical structure. There was greater gene flow before extirpation than that found among modern sea otter populations but the difference was not statistically significant. The most dramatic effect of pre fur trade population extirpation was the loss of genetic diversity. For long term conservation of these populations increasing gene flow and the maintenance of remnant genetic diversity should be encouraged.