Genetic structuring and recent demographic history of red pandas (Ailurus fulgens) inferred from microsatellite and mitochondrial DNA

Clarification of the genetic structure and population history of a species can shed light on the impacts of landscapes, historical climate change and contemporary human activities and thus enables evidence‐based conservation decisions for endangered organisms. The red panda (Ailurus fulgens) is an endangered species distributing at the edge of the Qinghai‐Tibetan Plateau and is currently subject to habitat loss, fragmentation and population decline, thus representing a good model to test the influences of the above‐mentioned factors on a plateau edge species. We combined nine microsatellite loci and 551 bp of mitochondrial control region (mtDNA CR) to explore the genetic structure and demographic history of this species. A total of 123 individuals were sampled from 23 locations across five populations. High levels of genetic variation were identified for both mtDNA and microsatellites. Phylogeographic analyses indicated little geographic structure, suggesting historically wide gene flow. However, microsatellite‐based Bayesian clustering clearly identified three groups (Qionglai‐Liangshan, Xiaoxiangling and Gaoligong‐Tibet). A significant isolation‐by‐distance pattern was detected only after removing Xiaoxiangling. For mtDNA data, there was no statistical support for a historical population expansion or contraction for the whole sample or any population except Xiaoxiangling where a signal of contraction was detected. However, Bayesian simulations of population history using microsatellite data did pinpoint population declines for Qionglai, Xiaoxiangling and Gaoligong, demonstrating significant influences of human activity on demography. The unique history of the Xiaoxiangling population plays a critical role in shaping the genetic structure of this species, and large‐scale habitat loss and fragmentation is hampering gene flow among populations. The implications of our findings for the biogeography of the Qinghai‐Tibetan Plateau, subspecies classification and conservation of red pandas are discussed.     

Short-term impact of disturbance on genetic diversity and structure of Indonesian populations of the butterfly Drupadia theda in East Kalimantan

We investigated the short‐term impact of disturbance on genetic diversity and structure of the tropical butterfly Drupadia theda Felder (Lepidoptera: Lycaenidae). Populations were sampled from five landscapes in East Kalimantan (Borneo, Indonesia) which were differentially disturbed by selective logging and the 1997/1998 El Niño Southern Oscillation (ENSO)‐induced drought and fires. Sampling occurred before (in 1997) and after the forest fires (in 1998, 1999, 2000, and 2004). Drupadia theda populations underwent serious population size reductions following the 1997/1998 ENSO event. For a total of 208 individuals, we sequenced a 509‐bp segment of mtDNA containing the control region plus the 5’ end of the 12S rDNA gene. Haplotype diversity in D. theda populations ranged from 0.468 to 0.953. Just after the 1997/1998 ENSO event, number of recorded individuals and genetic diversity were very low in D. theda populations sampled in the two severely burned areas and in a small pristine forest fragment that was surrounded by burned forest and thereby affected by drought. Interestingly, higher levels of genetic diversity were observed in logged forest compared to proximate pristine forest. After 1998, the genetic composition within the three ENSO‐disturbed areas diverged. In the twice‐burned forest, the genetic diversity in 1999 already approached pre‐fire levels, while it remained nearly unchanged in proximate once‐burned forest. Our data suggest that the 1997/1998 ENSO‐induced drought and fires caused massive reductions in the genetic diversity of D. theda and that population recoveries were linked to their geographical position relative to patches of unburned forest (and thus to source populations).     

Discrepancies in population differentiation at microsatellites, mitochondrial DNA and plumage colour in the pied flycatcher--inferring evolutionary processes

Genetic differentiation between three populations of the pied flycatcher Ficedula hypoleuca (Norway, Czech Republic and Spain, respectively) was investigated at microsatellite loci and mitochondrial DNA (mtDNA) sequences and compared with the pattern of differentiation of male plumage colour. The Czech population lives sympatrically with the closely related collared flycatcher (F. albicollis) whereas the other two are allopatric. Allopatric populations are on average more conspicuously coloured than sympatric ones, a pattern that has been explained by sexual selection for conspicuous colour in allopatry and a character displacement on breeding plumage colour in sympatry that reduces the rate of hybridization with the collared flycatcher. The Czech population was genetically indistinguishable from the Norwegian population at microsatellite loci and mtDNA sequences. Recent isolation and/or gene flow may explain the lack of genetic differentiation. Accordingly, different selection on plumage colour in the two populations is either sufficiently strong so that gene flow has little impact on the pattern of colour variation, or differentiation of plumage colour occurred so recently that the (presumably) neutral, fast evolving markers employed here are unable to reflect the differentiation. Genetically, the Spanish population was significantly differentiated from the other populations, but the divergence was much more pronounced at mtDNA compared to microsatellites. This may reflect increased rate of differentiation by genetic drift at the mitochondrial, compared with the nuclear genome, caused by the smaller effective population size of the former genome. In accordance with this interpretation, a genetic pattern consistent with effects of small population size in the Spanish population (genetic drift and inbreeding) were also apparent at the microsatellites, namely reduced allelic diversity and heterozygous deficiency.     

Relative information content of polymorphic microsatellites and mitochondrial DNA for inferring dispersal and population genetic structure in the olive sea snake, Aipysurus laevis

Polymorphic microsatellites are widely considered more powerful for resolving population structure than mitochondrial DNA (mtDNA) markers, particularly for recently diverged lineages or geographically proximate populations. Weaker population subdivision for biparentally inherited nuclear markers than maternally inherited mtDNA may signal male-biased dispersal but can also be attributed to marker-specific evolutionary characteristics and sampling properties. We discriminated between these competing explanations with a population genetic study on olive sea snakes, Aipysurus laevis. A previous mtDNA study revealed strong regional population structure for A. laevis around northern Australia, where Pleistocene sea-level fluctuations have influenced the genetic signatures of shallow-water marine species. Divergences among phylogroups dated to the Late Pleistocene, suggesting recent range expansions by previously isolated matrilines. Fine-scale population structure within regions was, however, poorly resolved for mtDNA. In order to improve estimates of fine-scale genetic divergence and to compare population structure between nuclear and mtDNA, 354 olive sea snakes (previously sequenced for mtDNA) were genotyped for five microsatellite loci. F statistics and Bayesian multilocus genotype clustering analyses found similar regional population structure as mtDNA and, after standardizing microsatellite F statistics for high heterozygosities, regional divergence estimates were quantitatively congruent between marker classes. Over small spatial scales, however, microsatellites recovered almost no genetic structure and standardized F statistics were orders of magnitude smaller than for mtDNA. Three tests for male-biased dispersal were not significant, suggesting that recent demographic expansions to the typically large population sizes of A. laevis have prevented microsatellites from reaching mutation-drift equilibrium and local populations may still be diverging.     

Patterns of population subdivision, gene flow and genetic variability in the African wild dog (Lycaon pictus)

African wild dogs are large, highly mobile carnivores that are known to disperse over considerable distances and are rare throughout much of their geographical range. Consequently, genetic variation within and differentiation between geographically separated populations is predicted to be minimal. We determined the genetic diversity of mitochondrial DNA (mtDNA) control region sequences and microsatellite loci in seven populations of African wild dogs. Analysis of mtDNA nucleotide diversity suggests that, historically, wild dog populations have been small relative to other large carnivores. However, population declines due to recent habitat loss have not caused a dramatic reduction in genetic diversity. We found one historical and eight recent mtDNA genotypes in 280 individuals that defined two highly divergent clades. In contrast to a previous, more limited, mtDNA analysis, sequences from these clades are not geographically restricted to eastern or southern African populations. Rather, we found a large admixture zone spanning populations from Botswana, Zimbabwe and south-eastern Tanzania. Mitochondrial and microsatellite differentiation between populations was significant and unique mtDNA genotypes and alleles characterized the populations. However, gene flow estimates (Nm) based on microsatellite data were generally greater than one migrant per generation. In contrast, gene flow estimates based on the mtDNA control region were lower than expected given differences in the mode of inheritance of mitochondrial and nuclear markers which suggests a male bias in long-distance dispersal.     

Range-wide population genetic analysis of the endangered northern riffleshell mussel, <Emphasis Type="Italic">Epioblasma torulosa rangiana</Emphasis> (Bivalvia: Unionoida)

The northern riffleshell (Epioblasma torulosa rangiana) is a critically endangered unionoid species in need of conservation throughout its range. It is the only unionoid to be federally protected in both Canada and the U.S. We use sequences from two mtDNA genes and 15 microsatellite loci to assess genetic variation among 86 individuals from the four populations in the three remaining drainages in which E. t. rangiana is known to be reproducing. All of these populations are in formerly glaciated landscapes that emerged <10 kya. The mtDNA sequence data do not indicate significant geographic structure among the populations. However, allelic data from the microsatellite loci reveal highly significant population structuring. Individuals of E. t. rangiana can be assigned to their own river of origin with 98.8% accuracy. Significant isolation-by-distance occurs. This analysis will be useful to conservation managers in documenting the genetic structure, patterns of isolation, and genetic variability within and among populations of E. t. rangiana.     

Microsatellite analysis supports mitochondrial phylogeography of the hellbender (Cryptobranchus alleganiensis)

We investigated genetic diversity of the hellbender (Cryptobranchus alleganiensis) throughout its range in the eastern US using nuclear markers and compared our results to a previously published mitochondrial analysis. A variety of nuclear markers, including protein-coding gene introns and microsatellites were tested but only microsatellites were variable enough for population level analysis. Microsatellite loci showed moderate among population sharing of alleles, in contrast to the reciprocal monophyly exhibited by mitochondrial DNA. However, analyses using F-statistics and Bayesian clustering algorithms showed considerable population subdivision and clustered hellbender populations into the same major groups as the mtDNA. The microsatellites combined with the mtDNA data suggest that gene flow is severely restricted or non-existent among eight major groups, and potentially among populations (rivers) within groups. The combined mtDNA and microsatellite data suggest that the currently recognized hellbender subspecies are paraphyletic. We suggest that the eight independent groups identified in our study should be managed as such, rather than basing conservation decisions on the two named subspecies of hellbender.     

Population genetic characterization of the endangered silver skipper,Leptalina unicolor (Lepidoptera: Hesperiidae), using novel microsatellite markers and mitochondrial DNA gene sequences

The silver stripped skipper, Leptalina unicolor, is listed as an endangered species with a paucity of leftover populations in Korea. Prior population genetic analysis is essential for the establishment of effective conservation strategies. In the present study, we investigated the genetic diversity and the relationships of L. unicolor populations in Korea using 12 newly developed L. unicolor-specific microsatellite markers and two mitochondrial DNA (mtDNA) sequences (a total of 1283 bp from COI and CytB). Due to the endangered status and rarity, a total of 33 individuals from two localities in Gangwon-do Province (Seohwa and Girin) and one in Gyeongsangnam-do Province (Danjang) were sampled. The concatenated mtDNA sequences revealed very low genetic diversity in each population, with only one to three independent haplotypes per population. F_ST, principal coordinate analysis, and independent structural analysis of the concatenated mtDNA sequences and the 12 microsatellite loci developed in this study showed that each L. unicolor population in Korea is nearly completely isolated, although inbreeding has not yet occurred. Long-term conservation habitat recovery that increases the available population appears to be important in the prevention of genetic drift and inbreeding, which occur in small isolated populations and reduce the viability of populations under fluctuating environmental conditions.     

Population stock structure of leatherback turtles (Dermochelys coriacea) in the Atlantic revealed using mtDNA and microsatellite markers

This study presents a comprehensive genetic analysis of stock structure for leatherback turtles (Dermochelys coriacea), combining 17 microsatellite loci and 763 bp of the mtDNA control region. Recently discovered eastern Atlantic nesting populations of this critically endangered species were absent in a previous survey that found little ocean-wide mtDNA variation. We added rookeries in West Africa and Brazil and generated longer sequences for previously analyzed samples. A total of 1,417 individuals were sampled from nine nesting sites in the Atlantic and SW Indian Ocean. We detected additional mtDNA variation with the longer sequences, identifying ten polymorphic sites that resolved a total of ten haplotypes, including three new variants of haplotypes previously described by shorter sequences. Population differentiation was substantial between all but two adjacent rookery pairs, and F _ST values ranged from 0.034 to 0.676 and 0.004 to 0.205 for mtDNA and microsatellite data respectively, suggesting that male-mediated gene flow is not as widespread as previously assumed. We detected weak (F _ST = 0.008 and 0.006) but significant differentiation with microsatellites between the two population pairs that were indistinguishable with mtDNA data. POWSIM analysis showed that our mtDNA marker had very low statistical power to detect weak structure (F _ST < 0.005), while our microsatellite marker array had high power. We conclude that the weak differentiation detected with microsatellites reflects a fine scale level of demographic independence that warrants recognition, and that all nine of the nesting colonies should be considered as demographically independent populations for conservation. Our findings illustrate the importance of evaluating the power of specific genetic markers to detect structure in order to correctly identify the appropriate population units to conserve.     

Population structure of the yellow-footed rock-wallaby Petrogale xanthopus (Gray, 1854) inferred from mtDNA sequences and microsatellite loci

The yellow-footed rock-wallaby Petrogale xanthopus is considered to be potentially vulnerable to extinction. This wallaby inhabits naturally disjunct rocky outcrops which could restrict dispersal between populations, but the extent to which that occurs is unknown. Genetic differences between populations were assessed using mitochondrial DNA (control region) sequencing and analysis of variation at four microsatellite loci among three geographically close sites in south-west Queensland (P. x. celeris) and, for mtDNA only, samples from South Australia (P. x. xanthopus) as well. Populations from South Australia and Queensland had phylogenetically distinct mtDNA, supporting the present classification of these two groups as evolutionarily distinct entities. Within Queensland, populations separated by 70 km of unsuitable habitat differed significantly for mtDNA and at microsatellite loci. Populations separated by 10 km of apparently suitable habitat had statistically homogeneous mtDNA, but a significant difference in allele frequency at one microsatellite locus. Tests for Hardy-Weinberg equilibrium and micro-geographical variation at microsatellite loci did not detect any substructuring between two wallaby aggregations within a colony encircling a single rock outcrop. Although the present study was limited by small sample sizes at two of the three Queensland locations examined, the genetic results suggest that dispersal between colonies is limited, consistent with an ecological study of dispersal at one of the sites. Considering both the genetic and ecological data, we suggest that management of yellow-footed rock-wallabies should treat each colony as an independent unit and that conservation of the Queensland and South Australian populations as separate entities is warranted.     

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.     

Changes in genetic diversity and differentiation in Red‐cockaded woodpeckers (Dryobates borealis) over the past century

Red‐cockaded woodpeckers (RCW; Dryobates borealis) declined after human activities reduced their fire‐maintained pine ecosystem to <3% of its historical range in the southeastern United States and degraded remaining habitat. An estimated 1.6 million RCW cooperative breeding groups declined to about 3,500 groups with no more than 10,000 birds by 1978. Management has increased RCW population abundances since they were at their lowest in the 1990s. However, no range‐wide study has been undertaken since then to investigate the impacts of this massive bottleneck or infer the effects of conservation management and recent demographic recoveries. We used mitochondrial DNA sequences (mtDNA) and nine nuclear microsatellite loci to determine if range‐wide demographic declines resulted in changes to genetic structure and diversity in RCW by comparing samples collected before 1970 (mtDNA data only), between 1992 and 1995 (mtDNA and microsatellites), and between 2010 and 2014 (mtDNA and microsatellites). We show that genetic diversity has been lost as detected by a reduction in the number of mitochondrial haplotypes. This reduction was apparent in comparisons of pre‐1970 mtDNA data with data from the 1992–1995 and 2010–2014 time points, with no change between the latter two time points in mtDNA and microsatellite analyses. The mtDNA data also revealed increases in range‐wide genetic differentiation, with a genetically panmictic population present throughout the southeastern United States in the pre‐1970s data and subsequent development of genetic structure that has remained unchanged since the 1990s. Genetic structure was also uncovered with the microsatellite data, which like the mtDNA data showed little change between the 1992–1995 and 2010–2014 data sets. Temporal haplotype networks revealed a consistent, star‐like phylogeny, suggesting that despite the overall loss of haplotypes, no phylogenetically distinct mtDNA lineages were lost when the population declined. Our results may suggest that management during the last two decades has prevented additional losses of genetic diversity.     

Population structure and impact of supportive breeding inferred from mitochondrial and microsatellite DNA analyses in land-locked Atlantic salmon Salmo salar L.

Four tributaries of Lake St-Jean (Québec, Canada) are used for spawning and juvenile habitat by land-locked Atlantic salmon. Spawning runs have drastically declined since the mid-1980s, and consequently, a supportive-breeding programme was undertaken in 1990. In this study, we analysed seven microsatellite loci and mtDNA, and empirically estimated effective population sizes to test the hypotheses that (i) fish spawning in different tributaries form genetically distinct populations and (ii) the supportive breeding programme causes genetic perturbations on wild populations. Allele frequency distribution, molecular variance and genetic distance estimates all supported the hypothesis of genetic differentiation among salmon from different tributaries. Gene flow among some populations was much more restricted than previously reported for anadromous populations despite the small geographical scale (40 km) involved. Both mtDNA and microsatellites revealed a more pronounced differentiation between populations from two tributaries of a single river compared with their differentiation with a population from a neighbouring river. The comparison of wild and F₁-hatchery fish (produced from breeders originating from the same river) indicated significant changes in allele frequencies and losses of low-frequency alleles but no reduction in heterozygosity. Estimates of variance and inbreeding population size indicated that susceptibility to genetic drift and inbreeding in one population increased by twofold after only one generation of supplementation.     

Using 2 Genetic Markers to Discriminate Among Canada Goose Populations in Ohio

Abstract: Canada goose (Branta canadensis) harvest management depends on reliable estimates of harvest composition, and established genetic methods provide an alternative to traditional methods. We expanded upon previous genetic studies by comparing the utility of 6 nuclear microsatellite loci and mitochondrial (mtDNA) control region sequences for discriminating among giant (B. c. maxima) and interior (B. c. interior) populations in Ohio (USA) Canada goose harvests at both individual and population levels. Subspecies and populations exhibited greater differentiation in mtDNA (F_ST = 0.202) than microsatellites (F_ST = 0.021), as would be expected based on differences in effective population size. Neither microsatellites nor mtDNA alone were sufficient for estimating harvest composition at the subspecies or population level in simulations and empirical blind tests using individuals of known origin; however, a combined microsatellite + mtDNA dataset yielded accurate and precise harvest derivations at the subspecies level. Both population-level mixed stock analysis and individual-level assignment tests provided accurate results, but a large proportion of birds could not be assigned with confidence at the individual level. We applied mixed stock analysis and the combined microsatellite + mtDNA dataset to Ohio's 2003–2004 harvest and found that interior populations accounted for 4.9% (95% CI = 1.7–8.0%) of the statewide early season and 9.3% (95% CI = 6.9–11.6%) of the regular and late-season harvested sample. These results suggest that maximum likelihood harvest derivations are highly dependent on the choice of genetic markers. Studies should only employ markers that exhibit sufficient variation and have been shown through simulations and empirical testing to accurately discriminate among the subspecies or management populations of interest.     

What Maintains the Central North Pacific Genetic Discontinuity in Pacific Herring?

Pacific herring show an abrupt genetic discontinuity in the central North Pacific that represents secondary contact between refuge populations previously isolated during Pleistocene glaciations. Paradoxically, high levels of gene flow produce genetic homogeneity among ocean-type populations within each group. Here, we surveyed variability in mtDNA control-region sequences (463 bp) and nine microsatellite loci in Pacific herring from sites across the North Pacific to further explore the nature of the genetic discontinuity around the Alaska Peninsula. Consistent with previous studies, little divergence (Φ_ST  = 0.011) was detected between ocean-type populations of Pacific herring in the North West Pacific, except for a population in the Yellow Sea (Φ_ST  = 0.065). A moderate reduction in genetic diversity for both mtDNA and microsatellites in the Yellow Sea likely reflects founder effects during the last colonization of this sea. Reciprocal monophyly between divergent mtDNA lineages (Φ_ST  = 0.391) across the Alaska Peninsula defines the discontinuity across the North Pacific. However, microsatellites did not show a strong break, as eastern Bering Sea (EBS) herring were more closely related to NE Pacific than to NW Pacific herring. This discordance between mtDNA and microsatellites may be due to microsatellite allelic convergence or to sex-biased dispersal across the secondary contact zone. The sharp discontinuity between Pacific herring populations may be maintained by high-density blocking, competitive exclusion or hybrid inferiority.     

Genetic population structure and introgression in Anopheles dirus mosquitoes in South-east Asia

Genetic structure and species relationships were studied in three closely related mosquito species, Anopheles dirus A, C and D in Thailand using 11 microsatellite loci and compared with previous mitochondrial DNA (mtDNA) data on the same populations. All three species were well differentiated from each other at the microsatellite loci. Given the almost complete absence of mtDNA differentiation between An. dirus A and D, this endorses the previous suggestion of mtDNA introgression between these species. The high degree of differentiation between the northern and southern population of An. dirus C (RST = 0.401), in agreement with mtDNA data, is suggestive of incipient species. The lack of genetic structure indicated by microsatellites in four populations of An. dirus A across northern Thailand also concurs with mtDNA data. However, in An. dirus D a limited but significant level of structure was detected by microsatellites over ~400 km in northern Thailand, whereas the mtDNA detected no population differentiation over a much larger area (>1200 km). There is prior evidence for population expansion in the mtDNA. If this is due to a selective sweep originating in An. dirus D, the microsatellite data may indicate greater barriers to gene flow within An. dirus D than in species A. Alternatively, there may have been historical introgression of mtDNA and subsequent demographic expansion which occurred first in An. dirus D so enabling it to accumulate some population differentiation. In the latter case the lack of migration-drift equilibrium precludes the inference of absolute or relative values of gene flow in An. dirus A and D.     

Population structure of the African savannah elephant inferred from mitochondrial control region sequences and nuclear microsatellite loci

Two hundred and thirty-six mitochondrial DNA nucleotide sequences were used in combination with polymorphism at four nuclear microsatellite loci to assess the amount and distribution of genetic variation within and between African savannah elephants. They were sampled from 11 localities in eastern, western and southern Africa. In the total sample, 43 haplotypes were identified and an overall nucleotide diversity of 2.0% was observed. High levels of polymorphism were also observed at the microsatellite loci both at the level of number of alleles and gene diversity. Nine to 14 alleles per locus across populations and 44 alleles in the total sample were found. The gene diversity ranged from 0.51 to 0.72 in the localities studied. An analysis of molecular variance showed significant genetic differentiation between populations within regions and also between regions. The extent of subdivision between populations at the mtDNA control region was approximately twice as high as shown by the microsatellite loci (mtDNA F(ST) = 0.59; microsatellite R(ST) = 0.31). We discuss our results in the light of Pleistocene refugia and attribute the observed pattern to population divergence in allopatry accompanied by a recent population admixture following a recent population expansion.     

Nuclear and mitochondrial markers reveal distinctiveness of a small population of bottlenose whales (Hyperoodon ampullatus) in the western North Atlantic

Small populations at the edge of a species' distribution can represent evolutionary relics left behind after range contractions due to climate change or human exploitation. The distinctiveness and genetic diversity of a small population of bottlenose whales in the Gully, a submarine canyon off Nova Scotia, was quantified by comparison to other North Atlantic populations using 10 microsatellites and mitrochondrial DNA (mtDNA) control region sequences (434 bp). Both markers confirmed the distinctiveness of the Gully (n = 34) from the next nearest population, off Labrador (n = 127; microsatellites -F(ST)= 0.0243, P < 0.0001; mtDNA -Phi(ST) = 0.0456, P < 0.05). Maximum likelihood microsatellite estimates suggest that less than two individuals per generation move between these areas, refuting the hypothesis of population links through seasonal migration. Both males and females appear to be philopatric, based on significant differentiation at both genomes and similar levels of structuring among the sexes for microsatellites. mtDNA diversity was very low in all populations (h = 0.51, pi = 0.14%), a pattern which may be due to selective sweeps associated with this species' extreme deep-diving ecology. Whaling had a substantial impact on bottlenose whale abundance, with over 65 000 animals killed before the hunt ceased in the early 1970s. Genetic diversity was similar among all populations, however, and no signal for bottlenecks was detected, suggesting that the Gully is not a relic of a historically wider distribution. Instead, this unique ecosystem appears to have long provided a stable year-round habitat for a distinct population of bottlenose whales.     

Inferring population history and demography using microsatellites, mitochondrial DNA, and major histocompatibility complex (MHC) genes

Microsatellites and mitochondrial DNA (mtDNA) have traditionally been used in population genetics because of their variability and presumed neutrality, whereas genes of the major histocompatibility complex (MHC) are increasingly of interest because strong selective pressures shape their standing variation. Despite the potential for MHC genes, microsatellites, and mtDNA sequences to complement one another in deciphering population history and demography, the three are rarely used in tandem. Here we report on MHC, microsatellite, and mtDNA variability in a single large population of the eastern tiger salamander (Ambystoma tigrinum tigrinum). We use the mtDNA mismatch distribution and, on microsatellite data, the imbalance index and bottleneck tests to infer aspects of population history and demography. Haplotype and allelic variation was high at all loci surveyed, and heterozygosity was high at the nuclear loci. We find concordance among neutral molecular markers that suggests our study population originated from post-Pleistocene expansions of multiple, fragmented sources that shared few migrants. Differences in N(e) estimates derived from haploid and diploid genetic markers are potentially attributable to secondary contact among source populations that experienced rapid mtDNA divergence and comparatively low levels of nuclear DNA divergence. We find strong evidence of natural selection acting on MHC genes and estimate long-term effective population sizes (N(e)) that are very large, making small selection intensities significant evolutionary forces in this population.