Genetic diversity loss in a biodiversity hotspot: ancient DNA quantifies genetic decline and former connectivity in a critically endangered marsupial

The extent of genetic diversity loss and former connectivity between fragmented populations are often unknown factors when studying endangered species. While genetic techniques are commonly applied in extant populations to assess temporal and spatial demographic changes, it is no substitute for directly measuring past diversity using ancient DNA (aDNA). We analysed both mitochondrial DNA (mtDNA) and nuclear microsatellite loci from 64 historical fossil and skin samples of the critically endangered Western Australian woylie (Bettongia penicillata ogilbyi), and compared them with 231 (n = 152 for mtDNA) modern samples. In modern woylie populations 15 mitochondrial control region (CR) haplotypes were identified. Interestingly, mtDNA CR data from only 29 historical samples demonstrated 15 previously unknown haplotypes and detected an extinct divergent clade. Through modelling, we estimated the loss of CR mtDNA diversity to be between 46% and 91% and estimated this to have occurred in the past 2000–4000 years in association with a dramatic population decline. In addition, we obtained near‐complete 11‐loci microsatellite profiles from 21 historical samples. In agreement with the mtDNA data, a number of ‘new’ microsatellite alleles was only detected in the historical populations despite extensive modern sampling, indicating a nuclear genetic diversity loss >20%. Calculations of genetic diversity (heterozygosity and allelic rarefaction) showed that these were significantly higher in the past and that there was a high degree of gene flow across the woylie's historical range. These findings have an immediate impact on how the extant populations are managed and we recommend the implementation of an assisted migration programme to prevent further loss of genetic diversity. Our study demonstrates the value of integrating aDNA data into current‐day conservation strategies.     

Mitochondrial DNA and nuclear microsatellites reveal high diversity and genetic structure in an avian top predator,the white‐tailed sea eagle,in central Europe

We analysed 123 white‐tailed sea eagles (Haliaeetus albicilla) from (primarily central) Europe with respect to variability and differentiation based on 499 bp of the mitochondrial control region and genotypes at seven unlinked nuclear microsatellites. Variability was high (overall expected heterozygosity, haplotype and nucleotide diversity being 0.70, 0.764 and 0.00698, respectively) and both marker systems showed a subdivision into two main genetic clusters (microsatellites) or haplogroups (mtDNA). In line with earlier analyses focusing on populations from northern and eastern Europe, as well as from Asia, we found a high level of admixture in Europe and no signs of a bottleneck – despite a severe decline of white‐tailed sea eagle populations during the 20^th century. Europe is thus a global stronghold for this species not only with respect to the number of breeding pairs but also regarding the proportion of species‐wide genetic diversity. Our dense sampling revealed a possibly clinal variation within central Europe from north‐west to south‐east that was reflected by the distribution of mtDNA haplotypes as well as the two microsatellite‐based clusters. This population differentiation in central Europe probably originated from a geographically structured postglacial colonization and was later enhanced by recent demographic fluctuations. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 99 , 727–737.     

Positive relationships between genetic diversity and abundance in fishes

Molecular markers, such as mitochondrial DNA and microsatellite loci, are widely studied to assess population genetics and phylogeography; however, the selective neutrality of these markers is increasingly being questioned. Given the importance of molecular markers in fisheries science and conservation, we evaluated the neutrality of both mtDNA and microsatellite loci through their associations with population size. We surveyed mtDNA and microsatellite data from the primary literature and determined whether genetic diversity increased with abundance across a total of 105 marine and freshwater fishes, with both global fisheries catch data and body size as proxies for abundance (with an additional 57 species for which only body size data were assessed). We found that microsatellite data generally yielded higher associations with abundance than mtDNA data, and within mtDNA analyses, number of haplotypes and haplotype diversity were more strongly associated with abundance than nucleotide diversity, particularly for freshwater fishes. We compared genetic diversity between freshwater and marine fishes and found that marine fishes had higher values of all measures of genetic diversity than freshwater fishes. Results for both mtDNA and microsatellites generally conformed to neutral expectations, although weaker relationships were often found between mtDNA nucleotide diversity and ‘abundance’ compared to any other genetic statistic. We speculate that this is because of historical events unrelated to natural selection, although a role for selection cannot be ruled out.     

Genetic structure and demographic history should inform conservation: Chinese cobras currently treated as homogenous show population divergence

An understanding of population structure and genetic diversity is crucial for wildlife conservation and for determining the integrity of wildlife populations. The vulnerable Chinese cobra (Naja atra) has a distribution from the mouth of the Yangtze River down to northern Vietnam and Laos, within which several large mountain ranges and water bodies may influence population structure. We combined 12 microsatellite loci and 1117 bp of the mitochondrial cytochrome b gene to explore genetic structure and demographic history in this species, using 269 individuals from various localities in Mainland China and Vietnam. High levels of genetic variation were identified for both mtDNA and microsatellites. mtDNA data revealed two main (Vietnam + southern China + southwestern China; eastern + southeastern China) and one minor (comprising only two individuals from the westernmost site) clades. Microsatellite data divided the eastern + southeastern China clade further into two genetic clusters, which include individuals from the eastern and southeastern regions, respectively. The Luoxiao and Nanling Mountains may be important barriers affecting the diversification of lineages. In the haplotype network of cytchrome b, many haplotypes were represented within a "star" cluster and this and other tests suggest recent expansion. However, microsatellite analyses did not yield strong evidence for a recent bottleneck for any population or genetic cluster. The three main clusters identified here should be considered as independent management units for conservation purposes. The release of Chinese cobras into the wild should cease unless their origin can be determined, and this will avoid problems arising from unnatural homogenization.     

The impacts of inbreeding,drift and selection on genetic diversity in captive breeding populations

The goal of captive breeding programmes is often to maintain genetic diversity until re‐introductions can occur. However, due in part to changes that occur in captive populations, approximately one‐third of re‐introductions fail. We evaluated genetic changes in captive populations using microsatellites and mtDNA. We analysed six populations of white‐footed mice that were propagated for 20 generations using two replicates of three protocols: random mating (RAN), minimizing mean kinship (MK) and selection for docility (DOC). We found that MK resulted in the slowest loss of microsatellite genetic diversity compared to RAN and DOC. However, the loss of mtDNA haplotypes was not consistent among replicate lines. We compared our empirical data to simulated data and found no evidence of selection. Our results suggest that although the effects of drift may not be fully mitigated, MK reduces the loss of alleles due to inbreeding more effectively than random mating or docility selection. Therefore, MK should be preferred for captive breeding. Furthermore, our simulations show that incorporating microsatellite data into the MK framework reduced the magnitude of drift, which may have applications in long‐term or extremely genetically depauperate captive populations.     

Genetic structure and phylogeography of a European flagship species,the white‐tailed sea eagle Haliaeetus albicilla

We analysed 120 white‐tailed sea eagles Haliaeetus albicilla from eastern (Poland and Estonia) and southeastern (Serbian Danube population) Europe for genetic variability and structuring at the mitochondrial control region and seven nuclear microsatellite loci. We combined this new dataset with sequence and genotype data from previous analyses covering Greenland and Eurasia (total sample sizes of 420 and 186 individuals for mtDNA and microsatellites, respectively) to address the following questions: 1) does the large eastern population in Europe add significantly to the species‘ overall genetic diversity? 2) Do the new sequence data match the clinal distribution pattern (west to east) of the two major mtDNA lineages? 3) Does the preliminary hypothesis of two nuclear genetic clusters recently found in this species hold for the whole of Europe, and do these clusters show a geographic pattern? Our results confirmed Europe as a stronghold of genetic diversity in white‐tailed sea eagles, and the east of the continent contributed disproportionately to this, the reason being the admixture of eagles with different genetic background. As hypothesised, both mitochondrial lineages were recovered also in eastern Europe, but the globally more eastern lineage was dominant. The presence of two microsatellite clusters was also confirmed, and these groups, too, show a non‐random geographic distribution, with, except for Poland, a high proportion of ‘eastern‐type’ eagles in the populations of east–central and eastern Europe.     

Genetic diversity of the Chinese goat in the littoral zone of the Yangtze River as assessed by microsatellite and mtDNA

The objective of this study was to assess the genetic diversity and population structure of goats in the Yangtze River region using microsatellite and mtDNA to better understand the current status of those goat genetic diversity and the effects of natural landscape in fashion of domestic animal genetic diversity. The genetic variability of 16 goat populations in the littoral zone of the Yangtze River was estimated using 21 autosomal microsatellites, which revealed high diversity and genetic population clustering with a dispersed geographical distribution. A phylogenetic analysis of the mitochondrial D‐loop region (482 bp) was conducted in 494 goats from the Yangtze River region. In total, 117 SNPs were reconstructed, and 173 haplotypes were identified, 94.5% of which belonged to lineages A and B. Lineages C, D, and G had lower frequencies (5.2%), and lineage F haplotypes were undetected. Several high‐frequency haplotypes were shared by different ecogeographically distributed populations, and the close phylogenetic relationships among certain low‐frequency haplotypes indicated the historical exchange of genetic material among these populations. In particular, the lineage G haplotype suggests that some west Asian goat genetic material may have been transferred to China via Muslim migration.     

Locals,resettlers, and pilgrims: A genetic portrait of three pre‐Columbian Andean populations

The common practice of resettlement and the development of administrative and ceremonial systems shaped the population landscape of the Andean region under the Inca rule. The area surrounding Coropuna and Solimana volcanoes, in the Arequipa region (Peru), carried a high‐density, multiethnic population. We studied the genetic variation among three pre‐Columbian populations from three functionally diverse archaeological sites excavated in this region. By analyzing the genetic composition of a large ceremonial center (Acchaymarca), an isolated pastoral settlement (Tompullo 2), and an agricultural settlement characterized by architectural features rare in the region (Puca), we investigated the patterns of population movements and the distribution of genetic diversity. We obtained mitochondrial DNA sequences for 25 individuals and autosomal microsatellite profiles for 20 individuals from Acchaymarca and Puca sites. These were compared with previously published genetic data for Tompullo 2 and other pre‐Columbian populations. We found differences among the genetic portraits of the three populations, congruent with the archaeologically described functions and characteristics of the sites. The Acchaymarca population had the highest genetic diversity and possessed the lowest number of unique mtDNA haplotypes. The Tompullo 2 population exhibited the lowest level of genetic diversity. The Puca population was distinct from the other two populations owing to a high frequency of haplogroup A haplotypes, what potentially explains the non‐local character of the burial architecture. Our analyses of microsatellite data suggest that gene flow between sites was mostly mediated by females, which is consistent with ethnohistorical knowledge of the social organization of the pre‐Columbian communities. Am J Phys Anthropol 154:402–412, 2014. © 2014 Wiley Periodicals, Inc.     

Three divergent lineages within an Australian marsupial (Petrogale penicillata) suggest multiple major refugia for mesic taxa in southeast Australia

Mesic southeastern Australia represents the continent's ancestral biome and is highly biodiverse, yet its phylogeographic history remains poorly understood. Here, we examine mitochondrial DNA (mtDNA) control region and microsatellite diversity in the brush‐tailed rock‐wallaby (Petrogale penicillata; n = 279 from 31 sites), to assess historic evolutionary and biogeographic processes in southeastern Australia. Our results (mtDNA, microsatellites) confirmed three geographically discrete and genetically divergent lineages within brush‐tailed rock‐wallabies, whose divergence appears to date to the mid‐Pleistocene. These three lineages had been hypothesized previously but data were limited. While the Northern and Central lineages were separated by a known biogeographic barrier (Hunter Valley), the boundary between the Central and Southern lineages was not. We propose that during particularly cool glacial cycles, the high peaks of the Great Dividing Range and the narrow adjacent coastal plain resulted in a more significant north–south barrier for mesic taxa in southeastern Australia than has been previously appreciated. Similarly, located phylogeographic breaks in codistributed species highlight the importance of these regions in shaping the distribution of biodiversity in southeastern Australia and suggest the existence of three major refuge areas during the Pleistocene. Substructuring within the northern lineage also suggests the occurrence of multiple local refugia during some glacial cycles. Within the three major lineages, most brush‐tailed rock‐wallaby populations were locally highly structured, indicating limited dispersal by both sexes. The three identified lineages represent evolutionarily significant units and should be managed to maximize the retention of genetic diversity within this threatened species.     

Y‐chromosome evidence supports widespread signatures of three‐species Canis hybridization in eastern North America

There has been considerable discussion on the origin of the red wolf and eastern wolf and their evolution independent of the gray wolf. We analyzed mitochondrial DNA (mtDNA) and a Y‐chromosome intron sequence in combination with Y‐chromosome microsatellites from wolves and coyotes within the range of extensive wolf–coyote hybridization, that is, eastern North America. The detection of divergent Y‐chromosome haplotypes in the historic range of the eastern wolf is concordant with earlier mtDNA findings, and the absence of these haplotypes in western coyotes supports the existence of the North American evolved eastern wolf (Canis lycaon). Having haplotypes observed exclusively in eastern North America as a result of insufficient sampling in the historic range of the coyote or that these lineages subsequently went extinct in western geographies is unlikely given that eastern‐specific mtDNA and Y‐chromosome haplotypes represent lineages divergent from those observed in extant western coyotes. By combining Y‐chromosome and mtDNA distributional patterns, we identified hybrid genomes of eastern wolf, coyote, gray wolf, and potentially dog origin in Canis populations of central and eastern North America. The natural contemporary eastern Canis populations represent an important example of widespread introgression resulting in hybrid genomes across the original C. lycaon range that appears to be facilitated by the eastern wolf acting as a conduit for hybridization. Applying conventional taxonomic nomenclature and species‐based conservation initiatives, particularly in human‐modified landscapes, may be counterproductive to the effective management of these hybrids and fails to consider their evolutionary potential.     

Conservation implications of the evolutionary history and genetic diversity hotspots of the snowshoe hare

With climate warming, the ranges of many boreal species are expected to shift northward and to fragment in southern peripheral ranges. To understand the conservation implications of losing southern populations, we examined range‐wide genetic diversity of the snowshoe hare (Lepus americanus), an important prey species that drives boreal ecosystem dynamics. We analysed microsatellite (8 loci) and mitochondrial DNA sequence (cytochrome b and control region) variation in almost 1000 snowshoe hares. A hierarchical structure analysis of the microsatellite data suggests initial subdivision in two groups, Boreal and southwestern. The southwestern group further splits into Greater Pacific Northwest and U.S. Rockies. The genealogical information retrieved from mtDNA is congruent with the three highly differentiated and divergent groups of snowshoe hares. These groups can correspond with evolutionarily significant units that might have evolved in separate refugia south and east of the Pleistocene ice sheets. Genetic diversity was highest at mid‐latitudes of the species' range, and genetic uniqueness was greatest in southern populations, consistent with substructuring inferred from both mtDNA and microsatellite analyses at finer levels of analysis. Surprisingly, snowshoe hares in the Greater Pacific Northwest mtDNA lineage were more closely related to black‐tailed jackrabbits (Lepus californicus) than to other snowshoe hares, which may result from secondary introgression or shared ancestral polymorphism. Given the genetic distinctiveness of southern populations and minimal gene flow with their northern neighbours, fragmentation and loss of southern boreal habitats could mean loss of many unique alleles and reduced evolutionary potential.     

High genetic diversity in a small population: the case of Chilean blue whales

It is generally assumed that species with low population sizes have lower genetic diversities than larger populations and vice versa. However, this would not be the case for long‐lived species with long generation times, and which populations have declined due to anthropogenic effects, such as the blue whale (Balaenoptera musculus). This species was intensively decimated globally to near extinction during the 20th century. Along the Chilean coast, it is estimated that at least 4288 blue whales were hunted from an apparently pre‐exploitation population size (k) of a maximum of 6200 individuals (Southeastern Pacific). Thus, here, we describe the mtDNA (control region) and nDNA (microsatellites) diversities of the Chilean blue whale aggregation site in order to verify the expectation of low genetic diversity in small populations. We then compare our findings with other blue whale aggregations in the Southern Hemisphere. Interestingly, although the estimated population size is small compared with the pre‐whaling era, there is still considerable genetic diversity, even after the population crash, both in mitochondrial (N = 46) and nuclear (N = 52) markers (H_d = 0.890 and H_o = 0.692, respectively). Our results suggest that this diversity could be a consequence of the long generation times and the relatively short period of time elapsed since the end of whaling, which has been observed in other heavily‐exploited whale populations. The genetic variability of blue whales on their southern Chile feeding grounds was similar to that found in other Southern Hemisphere blue whale feeding grounds. Our phylogenetic analysis of mtDNA haplotypes does not show extensive differentiation of populations among Southern Hemisphere blue whale feeding grounds. The present study suggests that although levels of genetic diversity are frequently used as estimators of population health, these parameters depend on the biology of the species and should be taken into account in a monitoring framework study to obtain a more complete picture of the conservation status of a population.     

Global genetic analysis reveals the putative native source of the invasive termite,Reticulitermes flavipes,in France

Biological invasions are recognized as a major threat to both natural and managed ecosystems. Phylogeographic and population genetic analyses can provide information about the geographical origins and patterns of introduction and explain the causes and mechanisms by which introduced species have become successful invaders. Reticulitermes flavipes is a North American subterranean termite that has been introduced into several areas, including France where introduced populations have become invasive. To identify likely source populations in the USA and to compare the genetic diversity of both native and introduced populations, an extensive molecular genetic study was undertaken using the COII region of mtDNA and 15 microsatellite loci. Our results showed that native northern US populations appeared well differentiated from those of the southern part of the US range. Phylogenetic analysis of both mitochondrial and nuclear markers showed that French populations probably originated from southeastern US populations, and more specifically from Louisiana. All of the mtDNA haplotypes shared between the United States and France were found in Louisiana. Compared to native populations in Louisiana, French populations show lower genetic diversity at both mtDNA and microsatellite markers. These findings are discussed along with the invasion routes of R. flavipes as well as the possible mechanisms by which French populations have evolved after their introduction.     

Reconstruction of caribou evolutionary history in Western North America and its implications for conservation

The role of Beringia as a refugium and route for trans-continental exchange of fauna during glacial cycles of the past 2million years are well documented; less apparent is its contribution as a significant reservoir of genetic diversity. Using mitochondrial DNA sequences and 14 microsatellite loci, we investigate the phylogeographic history of caribou (Rangifer tarandus) in western North America. Patterns of genetic diversity reveal two distinct groups of caribou. Caribou classified as a Northern group, of Beringian origin, exhibited greater number and variability in mtDNA haplotypes compared to a Southern group originating from refugia south of glacial ice. Results indicate that subspecies R. t. granti of Alaska and R. t. groenlandicus of northern Canada do not constitute distinguishable units at mtDNA or microsatellites, belying their current status as separate subspecies. Additionally, the Northern Mountain ecotype of woodland caribou (presently R. t. caribou) has closer kinship to caribou classified as granti or groenlandicus. Comparisons of mtDNA and microsatellite data suggest that behavioural and ecological specialization is a more recently derived life history characteristic. Notably, microsatellite differentiation among Southern herds is significantly greater, most likely as a result of human-induced landscape fragmentation and genetic drift due to smaller population sizes. These results not only provide important insight into the evolutionary history of northern species such as caribou, but also are important indicators for managers evaluating conservation measures for this threatened species.     

Y‐chromosome evidence supports asymmetric dog introgression into eastern coyotes

Hybridization has played an important role in the evolutionary history of Canis species in eastern North America. Genetic evidence of coyote–dog hybridization based on mitochondrial DNA (mtDNA) is lacking compared to that based on autosomal markers. This discordance suggests dog introgression into coyotes has potentially been male biased, but this hypothesis has not been formally tested. Therefore, we investigated biparentally, maternally, and paternally inherited genetic markers in a sample of coyotes and dogs from southeastern Ontario to assess potential asymmetric dog introgression into coyotes. Analysis of autosomal microsatellite genotypes revealed minimal historical and contemporary admixture between coyotes and dogs. We observed only mutually exclusive mtDNA haplotypes in coyotes and dogs, but we observed Y‐chromosome haplotypes (Y‐haplotypes) in both historical and contemporary coyotes that were also common in dogs. Species‐specific Zfy intron sequences of Y‐haplotypes shared between coyotes and dogs confirmed their homology and indicated a putative origin from dogs. We compared Y‐haplotypes observed in coyotes, wolves, and dogs profiled in multiple studies, and observed that the Y‐haplotypes shared between coyotes and dogs were either absent or rare in North American wolves, present in eastern coyotes, but absent in western coyotes. We suggest the eastern coyote has experienced asymmetric genetic introgression from dogs, resulting from predominantly historical hybridization with male dogs and subsequent backcrossing of hybrid offspring with coyotes. We discuss the temporal and spatial dynamics of coyote–dog hybridization and the conditions that may have facilitated the introgression of dog Y‐chromosomes into coyotes. Our findings clarify the evolutionary history of the eastern coyote.     

Evaluating DNA degradation rates in faecal pellets of the endangered pygmy rabbit

Noninvasive genetic sampling of faecal pellets can be a valuable method for monitoring rare and cryptic wildlife populations, like the pygmy rabbit (Brachylagus idahoensis). To investigate this method's efficiency for pygmy rabbit monitoring, we evaluated the effect of sample age on DNA degradation in faecal pellets under summer field conditions. We placed 275 samples from known individuals in natural field conditions for 1–60 days and assessed DNA quality by amplifying a 294‐base‐pair (bp) mitochondrial DNA (mtDNA) locus and five nuclear DNA (nDNA) microsatellite loci (111–221 bp). DNA degradation was influenced by sample age, DNA type, locus length and rabbit sex. Both mtDNA and nDNA exhibited high PCR success rates (94.4%) in samples <1 day old. Success rates for microsatellite loci declined rapidly from 80.0% to 42.7% between days 5 and 7, likely due to increased environmental temperature. Success rates for mtDNA amplification remained higher than nDNA over time, with moderate success (66.7%) at 21 days. Allelic dropout rates were relatively high (17.6% at <1 day) and increased to 100% at 60 days. False allele rates ranged from 0 to 30.0% and increased gradually over time. We recommend collecting samples as fresh as possible for individual identification during summer field conditions. Our study suggests that this method can be useful for future monitoring efforts, including occupancy surveys, individual identification, population estimation, parentage analysis and monitoring of genetic diversity both of a re‐introduced population in central Washington and across their range.     

Mitochondrial DNA D‐loop sequence variation in maternal lineages of Iranian native horses

To understand the origin and genetic diversity of Iranian native horses, mitochondrial DNA (mtDNA) D‐loop sequences were generated for 95 horses from five breeds sampled in eight geographical locations in Iran. Sequence analysis of a 247‐bp segment revealed a total of 27 haplotypes with 38 polymorphic sites. Twelve of 19 mtDNA haplogroups were identified in the samples. The most common haplotypes were found within haplogroup X2. Within‐population haplotype and nucleotide diversities of the five breeds ranged from 0.838 ± 0.056 to 0.974 ± 0.022 and 0.011 ± 0.002 to 0.021 ± 0.001 respectively, indicating a relatively high genetic diversity in Iranian horses. The identification of several ancient sequences common between the breeds suggests that the lineage of the majority of Iranian horse breeds is old and obviously originated from a vast number of mares. We found in all native Iranian horse breeds lineages of the haplogroups D and K, which is concordant with the previous findings of Asian origins of these haplogroups. The presence of haplotypes E and K in our study also is consistent with a geographical west–east direction of increasing frequency of these haplotypes and a genetic fusion in Iranian horse breeds.     

Population structure in the endangered cyprinid fish Pararasbora moltrechti in Taiwan, based on mitochondrial and microsatellite DNAs

The genetic structure of four populations of Pararasbora moltrechti, an endemic species of the Cyprinidae in Taiwan, was investigated based on the genetic variation of mtDNA Cyt b gene and five microsatellite loci. High haplotype diversity (h = 0.92) but low nucleotide diversity (0.004) in mtDNA was detected in this endangered species. In total, 33 haplotypes and four clusters were identified in its mtDNA. Nevertheless, low correspondence was found between geographical division and mtDNA clusters. In contrast, Bayesian cluster analysis of the microsatellite data identified four genetic groups and revealed highly structured populations. Significantly negative Tajima's D statistics and mismatch distribution analyses suggest that P. moltrechti populations may have experienced a demographic expansion. In light of the results of a nested clade analysis of mtDNA haplotypes, we conclude that recent population fluctuations and restricted gene flow played major roles in shaping the spatial genetic structure of P. moltrechti populations.     

Genetic characterization and population bottleneck in the Hucul horse based on microsatellite and mitochondrial data

The aim of this work was to gather information about the origin and genetic characterization of the Central European Hucul horse based on 71 horses using 17 microsatellites and the D‐loop region of mtDNA. Their genetic relationship to the Polish Konik (N = 7), German (N = 4) and Hungarian wild Przewalski horses (N = 4) and 200 horse sequences from GenBank was also analysed. Both microsatellite and mtDNA analysis showed a high genetic variation in the Hucul. A total of 130 alleles were detected, the mean number of observed alleles per microsatellite was 7.647, and the number of effective alleles was 4.401. The average observed and expected heterozygosity were 0.706 and 0.747, respectively. The high heterozygosity values and Wright's fixation index (F_IS) (?0.128) indicated a low level of inbreeding, low or no selection pressure, and large number of alleles. mtDNA analysis revealed 18 haplotypes for the Hucul population with a total of 23 variable sites. Haplotype and nucleotide diversities were 0.935 ± 0.011 and 0.022 ± 0.012, respectively. Neutrality tests (Tajima's D and Fu's Fs) were non‐significant, and mismatch distribution was ragged, indicating that the Hucul population is in genetic equilibrium. The most frequent mtDNA D‐loop region belonged to haplogroup A (48%), which was also present in Przewalski Wild horse samples, while Polish Konik samples belonged to three haplotypes and C, F, and G haplogroups. Large and significant pairwise Φ_ST values along with a small number of common haplotypes indicated a low level of gene flow and lack of genetic structure among the three studied breeds (Hucul, Konik, and Przewalski Wild horse). The present work contributes to our knowledge of the genetic diversity of the Hucul horse and helps to define its genetic conservation. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2013, 109 , 54–65.     

Contrasting patterns of genetic diversity at three different genetic markers in a marine mammal metapopulation

Many studies use genetic markers to explore population structure and variability within species. However, only a minority use more than one type of marker and, despite increasing evidence of a link between heterozygosity and individual fitness, few ask whether diversity correlates with population trajectory. To address these issues, we analysed data from the Steller's sea lion, Eumetiopias jubatus , where three stocks are distributed over a vast geographical range and where both genetic samples and detailed demographic data have been collected from many diverse breeding colonies. To previously published mitochondrial DNA (mtDNA) and microsatellite data sets, we have added new data for amplified fragment length polymorphism (AFLP) markers, comprising 238 loci scored in 285 sea lions sampled from 23 natal rookeries. Genotypic diversity was low relative to most vertebrates, with only 37 loci (15.5%) being polymorphic. Moreover, contrasting geographical patterns of genetic diversity were found at the three markers, with Nei's gene diversity tending to be higher for AFLPs and microsatellites in rookeries of the western and Asian stocks, while the highest mtDNA values were found in the eastern stock. Overall, and despite strongly contrasting demographic histories, after applying phylogenetic correction we found little correlation between genetic diversity and either colony size or demography. In contrast, we were able to show a highly significant positive relationship between AFLP diversity and current population size across a range of pinniped species, even though equivalent analyses did not reveal significant trends for either microsatellites or mtDNA.