Connectivity and population subdivision at the fringe of a large brown bear (Ursus arctos) population in North Western Europe

Identifying priority core habitats and corridors for effective conservation of brown bears in Iran

Iran lies at the southernmost range limit of brown bears globally. Therefore, understanding the habitat associations and patterns of population connectivity for brown bears in Iran is relevant for the species’ conservation. We applied species distribution modeling to predict habitat suitability and connectivity modeling to identify population core areas and corridors. Our results showed that forest density, topographical roughness, NDVI and human footprint were the most influential variables in predicting brown bear distribution. The most crucial core areas and corridor networks for brown bear are concentrated in the Alborz and Zagros Mountains. These two core areas were predicted to be fragmented into a total of fifteen isolated patches if dispersal of brown bear across the landscape is limited to 50,000 cost units, and aggregates into two isolated habitat patches if the species is capable of dispersing 400,000 cost units. We found low overlap between corridors, and core habitats with protected areas, suggesting that the existing protected area network may not be adequate for the conservation of brown bear in Iran. Our results suggest that effective conservation of brown bears in Iran requires protection of both core habitats and the corridors between them, especially outside Iran’s network of protected areas.

Phylogeography, genetic diversity, and connectivity of brown bear populations in Central Asia

Knowledge of genetic diversity and population structure is critical for conservation and management planning at the population level within a species' range. Many brown bear populations in Central Asia are small and geographically isolated, yet their phylogeographic relationships, genetic diversity, and contemporary connectivity are poorly understood. To address this knowledge gap, we collected brown bear samples from the Gobi Desert (n = 2360), Altai, Sayan, Khentii, and Ikh Khyangan mountains of Mongolia (n = 79), and Deosai National Park in the Himalayan Mountain Range of Pakistan (n = 5) and generated 927 base pairs of mitochondrial DNA (mtDNA) sequence data and genotypes at 13 nuclear DNA microsatellite loci. We documented high levels of mtDNA and nDNA diversity in the brown bear populations of northern Mongolia (Altai, Sayan, Buteeliin nuruu and Khentii), but substantially lower diversity in brown bear populations in the Gobi Desert and Himalayas of Pakistan. We detected 3 brown bear mtDNA phylogeographic groups among bears of the region, with clade 3a1 in Sayan, Khentii, and Buteeliin nuruu mountains, clade 3b in Altai, Sayan, Buteeliin nuruu, Khentii, and Ikh Khyangan, and clade 6 in Gobi and Pakistan. Our results also clarified the phylogenetic relationships and divergence times with other brown bear mtDNA clades around the world. The nDNA genetic structure analyses revealed distinctiveness of Gobi bears and different population subdivisions compared to mtDNA results. For example, genetic distance for nDNA microsatellite loci between the bears in Gobi and Altai (FST = 0.147) was less than that of the Gobi and Pakistan (FST = 0.308) suggesting more recent male-mediated nuclear gene flow between Gobi and Altai than between Gobi and the Pakistan bears. Our results provide valuable information for conservation and management of bears in this understudied region of Central Asia and highlight the need for special protection and additional research on Gobi brown bears.

Genetic substructure and admixture as important factors in linkage disequilibrium-based estimation of effective number of breeders in recovering wildlife populations

The number of effective breeders (Nb ) and effective population size (Ne ) are population parameters reflective of evolutionary potential, susceptibility to stochasticity, and viability. We have estimated these parameters using the linkage disequilibrium-based approach with LDNE through the latest phase of population recovery of the brown bears (Ursus arctos) in Finland (1993-2010; N = 621). This phase of the recovery was recently documented to be associated with major changes in genetic composition. In particular, differentiation between the northern and the southern genetic cluster declined rapidly within 1.5 generations. Based on this, we have studied effects of the changing genetic structure on Nb and Ne , by comparing estimates for whole Finland with the estimates for the two genetic clusters. We expected a potentially strong relationship between estimate sizes and genetic differentiation, which should disappear as the population recovers and clusters merge. Consistent with this, our estimates for whole Finland were lower than the sum of the estimates of the two genetic clusters and both approaches produced similar estimates in the end. Notably, we also found that admixed genotypes strongly increased the estimates. In all analyses, our estimates for Ne were larger than Nb and likely reflective for brown bears of the larger region of Finland and northwestern Russia. Conclusively, we find that neglecting genetic substructure may lead to a massive underestimation of Nb and Ne . Our results also suggest the need for further empirical analysis focusing on individuals with admixed genotypes and their potential high influence on Nb and Ne .

Wolf population genetics in Europe: a systematic review, meta-analysis and suggestions for conservation and management

The grey wolf (Canis lupus) is an iconic large carnivore that has increasingly been recognized as an apex predator with intrinsic value and a keystone species. However, wolves have also long represented a primary source of human-carnivore conflict, which has led to long-term persecution of wolves, resulting in a significant decrease in their numbers, genetic diversity and gene flow between populations. For more effective protection and management of wolf populations in Europe, robust scientific evidence is crucial. This review serves as an analytical summary of the main findings from wolf population genetic studies in Europe, covering major studies from the 'pre-genomic era' and the first insights of the 'genomics era'. We analyse, summarize and discuss findings derived from analyses of three compartments of the mammalian genome with different inheritance modes: maternal (mitochondrial DNA), paternal (Y chromosome) and biparental [autosomal microsatellites and single nucleotide polymorphisms (SNPs)]. To describe large-scale trends and patterns of genetic variation in European wolf populations, we conducted a meta-analysis based on the results of previous microsatellite studies and also included new data, covering all 19 European countries for which wolf genetic information is available: Norway, Sweden, Finland, Estonia, Latvia, Lithuania, Poland, Czech Republic, Slovakia, Germany, Belarus, Russia, Italy, Croatia, Bulgaria, Bosnia and Herzegovina, Greece, Spain and Portugal. We compared different indices of genetic diversity in wolf populations and found a significant spatial trend in heterozygosity across Europe from south-west (lowest genetic diversity) to north-east (highest). The range of spatial autocorrelation calculated on the basis of three characteristics of genetic diversity was 650-850 km, suggesting that the genetic diversity of a given wolf population can be influenced by populations up to 850 km away. As an important outcome of this synthesis, we discuss the most pressing issues threatening wolf populations in Europe, highlight important gaps in current knowledge, suggest solutions to overcome these limitations, and provide recommendations for science-based wolf conservation and management at regional and Europe-wide scales.

Genetic and demographic recovery of an isolated population of brown bear Ursus arctos L., 1758

The brown bear Ursus arctos L., 1758 population of the Cantabrian Mountains (northwestern Spain) became isolated from other bear populations in Europe about 500 years ago and has declined due to hunting and habitat degradation. At the beginning of the 20th century, the Cantabrian population split into eastern and western subpopulations, and genetic exchange between them ceased. In the early 1990s, total population size was estimated to be < 100 bears. Subsequently, reduction in human-caused mortality has brought about an increase in numbers, mainly in the western subpopulation, likely promoting male-mediated migration and gene flow from the western nucleus to the eastern. To evaluate the possible genetic recovery of the small and genetically depauperate eastern subpopulation, in 2013 and 2014 we genotyped hair and faeces samples (116 from the eastern subpopulation and 36 from the western) for 18 microsatellite markers. Data from the annual count of females with cubs of the year (COY) during the past twenty-six years was used to analyze demographic changes. The number of females with COY fell to a minimum of seven in the western and three in eastern subpopulations in the biennium 1993-1994 and reached a respective maximum of 54 and 10 individuals in 2013-2014. We also observed increased bear dispersal and gene flow, mainly from the western to the eastern subpopulation. Of the 26 unique genotypes detected in the eastern subpopulation, 14 (54%) presented an admixture composition, and seven (27%) were determined to be migrants from the western subpopulation. Hence, the two separated and clearly structured subpopulations identified in the past currently show some degree of genetic admixture. This research shows the partial demographic recovery and a change in genetic composition due to migration process in a population of bears that has been isolated for several centuries.

Evidence of rapid change in genetic structure and diversity during range expansion in a recovering large terrestrial carnivore

Recovery of natural populations occurs often with simultaneous or subsequent range expansions. According to population genetic theory, genetic structuring emerges at the expansion front together with decreasing genetic diversity, owing to multiple founder events. Thereupon, as the expansion proceeds and connectivity among populations is established, homogenization and a resurgence of genetic diversity are to be expected. Few studies have used a fine temporal scale combined with genetic sampling to track range expansions as they proceed in wild animal populations. As a natural experiment, the historical eradication of large terrestrial carnivores followed by their recovery and recolonization may facilitate empirical tests of these ideas. Here, using brown bear (Ursus arctos) as model species, we tested predictions from genetic theory of range expansion. Individuals from all over Finland were genotyped for every year between 1996 and 2010 using 12 validated autosomal microsatellite markers. A latitudinal shift of about 110 km was observed in the distribution and delineation of genetic clusters during this period. As the range expansion proceeded, we found, as theory predicts, that the degree of genetic structure decreased, and that both genetic variation and admixture increased. The genetic consequences of range expansions may first be detected after multiple generations, but we found major changes in genetic composition after just 1.5 generations, accompanied by population growth and increased migration. These rapid genetic changes suggest an ongoing concerted action of geographical and demographic expansion combined with substantial immigration of bears from Russia during the recovery of brown bears within the large ecosystem of northern Europe.

A European Concern? Genetic Structure and Expansion of Golden Jackals (Canis aureus) in Europe and the Caucasus

In the first continent-wide study of the golden jackal (Canis aureus), we characterised its population genetic structure and attempted to identify the origin of European populations. This provided a unique insight into genetic characteristics of a native carnivore population with rapid large-scale expansion. We analysed 15 microsatellite markers and a 406 base-pair fragment of the mitochondrial control region. Bayesian-based and principal components methods were applied to evaluate whether the geographical grouping of samples corresponded with genetic groups. Our analysis revealed low levels of genetic diversity, reflecting the unique history of the golden jackal among Europe’s native carnivores. The results suggest ongoing gene flow between south-eastern Europe and the Caucasus, with both contributing to the Baltic population, which appeared only recently. The population from the Peloponnese Peninsula in southern Greece forms a common genetic cluster with samples from south-eastern Europe (ΔK approach in STRUCTURE, Principal Components Analysis [PCA]), although the results based on BAPS and the estimated likelihood in STRUCTURE indicate that Peloponnesian jackals may represent a distinct population. Moreover, analyses of population structure also suggest either genetic distinctiveness of the island population from Samos near the coast of Asia Minor (BAPS, most STRUCTURE, PCA), or possibly its connection with the Caucasus population (one analysis in STRUCTURE). We speculate from our results that ancient Mediterranean jackal populations have persisted to the present day, and have merged with jackals colonising from Asia. These data also suggest that new populations of the golden jackal may be founded by long-distance dispersal, and thus should not be treated as an invasive alien species, i.e. an organism that is “non-native to an ecosystem, and which may cause economic or environmental harm or adversely affect human health”. These insights into the genetic structure and ancestry of Baltic jackals have important implications for management and conservation of jackals in Europe. The golden jackal is listed as an Annex V species in the EU Habitats Directive and as such, considering also the results presented here, should be legally protected in all EU member states.

Y-chromosomal testing of brown bears (Ursus arctos): Validation of a multiplex PCR-approach for nine STRs suitable for fecal and hair samples

High-resolution Y-chromosomal markers have been applied to humans and other primates to study population genetics, migration, social structures and reproduction. Y-linked markers allow the direct assessment of the genetic structure and gene flow of uniquely male inherited lineages and may also be useful for wildlife conservation and forensics, but have so far been available only for few wild species. Thus, we have developed two multiplex PCR reactions encompassing nine Y-STR markers identified from the brown bear (Ursus arctos) and tested them on hair, fecal and tissue samples. The multiplex PCR approach was optimized and analyzed for species specificity, sensitivity and stutter-peak ratios. The nine Y-STRs also showed specific STR-fragments for male black bears and male polar bears, while none of the nine markers produced any PCR products when using DNA from female bears or males from 12 other mammals. The multiplex PCR approach in two PCR reactions could be amplified with as low as 0.2 ng template input. Precision was high in DNA templates from hairs, fecal scats and tissues, with standard deviations less than 0.14 and median stutter ratios from 0.04 to 0.63. Among the eight di- and one tetra-nucleotide repeat markers, we detected simple repeat structures in seven of the nine markers with 9-25 repeat units. Allelic variation was found for eight of the nine Y-STRs, with 2-9 alleles for each marker and a total of 36 alleles among 453 male brown bears sampled mainly from Northern Europe. We conclude that the multiplex PCR approach with these nine Y-STRs would provide male bear Y-chromosomal specificity and evidence suited for samples from conservation and wildlife forensics.

Recovery of large carnivores in Europe's modern human-dominated landscapes

The conservation of large carnivores is a formidable challenge for biodiversity conservation. Using a data set on the past and current status of brown bears (Ursus arctos), Eurasian lynx (Lynx lynx), gray wolves (Canis lupus), and wolverines (Gulo gulo) in European countries, we show that roughly one-third of mainland Europe hosts at least one large carnivore species, with stable or increasing abundance in most cases in 21st-century records. The reasons for this overall conservation success include protective legislation, supportive public opinion, and a variety of practices making coexistence between large carnivores and people possible. The European situation reveals that large carnivores and people can share the same landscape.

Admixture and gene flow from Russia in the recovering Northern European brown bear (Ursus arctos)

Large carnivores were persecuted to near extinction during the last centuries, but have now recovered in some countries. It has been proposed earlier that the recovery of the Northern European brown bear is supported by migration from Russia. We tested this hypothesis by obtaining for the first time continuous sampling of the whole Finnish bear population, which is located centrally between the Russian and Scandinavian bear populations. The Finnish population is assumed to experience high gene flow from Russian Karelia. If so, no or a low degree of genetic differentiation between Finnish and Russian bears could be expected. We have genotyped bears extensively from all over Finland using 12 validated microsatellite markers and compared their genetic composition to bears from Russian Karelia, Sweden, and Norway. Our fine masked investigation identified two overlapping genetic clusters structured by isolation-by-distance in Finland (pairwise FST = 0.025). One cluster included Russian bears, and migration analyses showed a high number of migrants from Russia into Finland, providing evidence of eastern gene flow as an important driver during recovery. In comparison, both clusters excluded bears from Sweden and Norway, and we found no migrants from Finland in either country, indicating that eastern gene flow was probably not important for the population recovery in Scandinavia. Our analyses on different spatial scales suggest a continuous bear population in Finland and Russian Karelia, separated from Scandinavia.

Brown and polar bear Y chromosomes reveal extensive male-biased gene flow within brother lineages

Brown and polar bears have become prominent examples in phylogeography, but previous phylogeographic studies relied largely on maternally inherited mitochondrial DNA (mtDNA) or were geographically restricted. The male-specific Y chromosome, a natural counterpart to mtDNA, has remained underexplored. Although this paternally inherited chromosome is indispensable for comprehensive analyses of phylogeographic patterns, technical difficulties and low variability have hampered its application in most mammals. We developed 13 novel Y-chromosomal sequence and microsatellite markers from the polar bear genome and screened these in a broad geographic sample of 130 brown and polar bears. We also analyzed a 390-kb-long Y-chromosomal scaffold using sequencing data from published male ursine genomes. Y chromosome evidence support the emerging understanding that brown and polar bears started to diverge no later than the Middle Pleistocene. Contrary to mtDNA patterns, we found 1) brown and polar bears to be reciprocally monophyletic sister (or rather brother) lineages, without signals of introgression, 2) male-biased gene flow across continents and on phylogeographic time scales, and 3) male dispersal that links the Alaskan ABC islands population to mainland brown bears. Due to female philopatry, mtDNA provides a highly structured estimate of population differentiation, while male-biased gene flow is a homogenizing force for nuclear genetic variation. Our findings highlight the importance of analyzing both maternally and paternally inherited loci for a comprehensive view of phylogeographic history, and that mtDNA-based phylogeographic studies of many mammals should be reevaluated. Recent advances in sequencing technology render the analysis of Y-chromosomal variation feasible, even in nonmodel organisms.

A sensitive and specific multiplex PCR approach for sex identification of ursine and tremarctine bears suitable for non-invasive samples

We report a new approach for molecular sex identification of extant Ursinae and Tremarctinae bears. Two Y-specific fragments (SMCY and 318.2) and one X-specific fragment (ZFX) are amplified in a multiplex PCR, yielding a double test for male-specific amplification and an internal positive control. The primers were designed and tested to be bear-specific, thereby minimizing the risk of cross-amplification in other species including humans. The high sensitivity and small amplicon sizes (100, 124, 160 base pairs) facilitate analysis of non-invasively obtained DNA material. DNA from tissue and blood as well as from 30 non-invasively collected hair and faeces yielded clear and easily interpretable results. The fragments were detected both by standard gel electrophoresis and automated capillary electrophoresis.

Lessons learned from microsatellite development for nonmodel organisms using 454 pyrosequencing

Microsatellites, also known as simple sequence repeats (SSRs), are among the most commonly used marker types in evolutionary and ecological studies. Next Generation Sequencing techniques such as 454 pyrosequencing allow the rapid development of microsatellite markers in nonmodel organisms. 454 pyrosequencing is a straightforward approach to develop a high number of microsatellite markers. Therefore, developing microsatellites using 454 pyrosequencing has become the method of choice for marker development. Here, we describe a user friendly way of microsatellite development from 454 pyrosequencing data and analyse data sets of 17 nonmodel species (plants, fungi, invertebrates, birds and a mammal) for microsatellite repeats and flanking regions suitable for primer development. We then compare the numbers of successfully lab-tested microsatellite markers for the various species and furthermore describe diverse challenges that might arise in different study species, for example, large genome size or nonpure extraction of genomic DNA. Successful primer identification was feasible for all species. We found that in species for which large repeat numbers are uncommon, such as fungi, polymorphic markers can nevertheless be developed from 454 pyrosequencing reads containing small repeat numbers (five to six repeats). Furthermore, the development of microsatellite markers for species with large genomes was also with Next Generation Sequencing techniques more cost and time-consuming than for species with smaller genomes. In this study, we showed that depending on the species, a different amount of 454 pyrosequencing data might be required for successful identification of a sufficient number of microsatellite markers for ecological genetic studies.

Limited gene flow among brown bear populations in far Northern Europe? Genetic analysis of the east-west border population in the Pasvik Valley

Noninvasively collected genetic data can be used to analyse large-scale connectivity patterns among populations of large predators without disturbing them, which may contribute to unravel the species' roles in natural ecosystems and their requirements for long-term survival. The demographic history of brown bears (Ursus arctos) in Northern Europe indicates several extinction and recolonization events, but little is known about present gene flow between populations of the east and west. We used 12 validated microsatellite markers to analyse 1580 hair and faecal samples collected during six consecutive years (2005-2010) in the Pasvik Valley at 70°N on the border of Norway, Finland and Russia. Our results showed an overall high correlation between the annual estimates of population size (N(c) ), density (D), effective size (N(e) ) and N(e) /N(c) ratio. Furthermore, we observed a genetic heterogeneity of ∼0.8 and high N(e) /N(c) ratios of ∼0.6, which suggests gene flow from the east. Thus, we expanded the population genetic study to include Karelia (Russia, Finland), Västerbotten (Sweden) and Troms (Norway) (477 individuals in total) and detected four distinct genetic clusters with low migration rates among the regions. More specifically, we found that differentiation was relatively low from the Pasvik Valley towards the south and east, whereas, in contrast, moderately high pairwise F(ST) values (0.91-0.12) were detected between the east and the west. Our results indicate ongoing limits to gene flow towards the west, and the existence of barriers to migration between eastern and western brown bear populations in Northern Europe.