Improved noninvasive genotyping method: application to brown bear (Ursus arctos) faeces

We redesigned new microsatellite primers and one sex‐specific primer for amplification of faecal DNA from brown bears (Ursus arctos). We also combined a semi‐nested polymerase chain reaction (PCR) with a newly developed multiplex preamplification method in order to increase the quality of the amplified DNA fragments. In comparison with a conventional PCR approach, the genotyping error rate was substantially reduced and the amplification rate was increased. This new approach could be transposed to other species where conventional PCR methods experience low success due to limited DNA concentration and/or quality.     

An evaluation of long-term preservation methods for brown bear (Ursus arctos) faecal DNA samples

Relatively few large-scale faecal DNA studieshave been initiated due to difficulties inamplifying low quality and quantity DNAtemplate. To improve brown bear faecal DNA PCRamplification success rates and to determinepost collection sample longevity, fivepreservation methods were evaluated: 90%ethanol, DETs buffer, silica-dried, oven-driedstored at room temperature, and oven-driedstored at –20 °C. Preservationeffectiveness was evaluated for 50 faecalsamples by PCR amplification of a mitochondrialDNA (mtDNA) locus (146 bp) and a nuclear DNA(nDNA) locus (200 bp) at time points of oneweek, one month, three months and six months. Preservation method and storage timesignificantly impacted mtDNA and nDNAamplification success rates. For mtDNA, allpreservation methods had 75% success atone week, but storage time had a significantimpact on the effectiveness of the silicapreservation method. Ethanol preserved sampleshad the highest success rates for both mtDNA(86.5%) and nDNA (84%). Nuclear DNAamplification success rates ranged from 26–88%, and storage time had a significant impacton all methods but ethanol. Preservationmethod and storage time should be importantconsiderations for researchers planningprojects utilizing faecal DNA. We recommendpreservation of faecal samples in 90% ethanolwhen feasible, although when collecting inremote field conditions or for both DNA andhormone assays a dry collection method may beadvantageous.     

Monitoring the effective population size of a brown bear (Ursus arctos) population using new single-sample approaches

The effective population size (N(e) ) could be the ideal parameter for monitoring populations of conservation concern as it conveniently summarizes both the evolutionary potential of the population and its sensitivity to genetic stochasticity. However, tracing its change through time is difficult in natural populations. We applied four new methods for estimating N(e) from a single sample of genotypes to trace temporal change in N(e) for bears in the Northern Dinaric Mountains. We genotyped 510 bears using 20 microsatellite loci and determined their age. The samples were organized into cohorts with regard to the year when the animals were born and yearly samples with age categories for every year when they were alive. We used the Estimator by Parentage Assignment (EPA) to directly estimate both N(e) and generation interval for each yearly sample. For cohorts, we estimated the effective number of breeders (N(b) ) using linkage disequilibrium, sibship assignment and approximate Bayesian computation methods and extrapolated these estimates to N(e) using the generation interval. The N(e) estimate by EPA is 276 (183-350 95% CI), meeting the inbreeding-avoidance criterion of N(e) > 50 but short of the long-term minimum viable population goal of N(e) > 500. The results obtained by the other methods are highly consistent with this result, and all indicate a rapid increase in N(e) probably in the late 1990s and early 2000s. The new single-sample approaches to the estimation of N(e) provide efficient means for including N(e) in monitoring frameworks and will be of great importance for future management and conservation.     

Non-invasive genetic study of the endangered Cantabrian brown bear (Ursus arctos)

The Brown Bear (Ursus arctos) population present in the Cantabrian Mountains has suffered a dramatic decline in recent centuries and is now threatened with extinction. This situation has led to the development and implementation of a species recovery plan. To accomplish this plan, we need to improve our knowledge about the ecology, demography and genetics of this population. This paper presents the genetic analysis of the Cantabrian brown bear population using non-invasive samples (faeces and hairs) collected between 2004 and 2006. It was necessary to optimize a set of 18 microsatellite loci and a sex marker (several new multiplex reactions were developed) to obtain a suitable probability of identity among genotypes to work with this small, deeply structured population. Genotyping of 48 individuals was carried out using a two-step PCR protocol to increase the quality of the multilocus genotypes. Validation of genotypes was performed using a multi-tube approach combined with different software programmes to measure their error rate and reliability. Diversity in the Cantabrian population was low (H _e = 0.51) and the population was markedly subdivided into two subpopulations (western and eastern) without current gene flow between them. The level of divergence between the two subpopulations (F _st = 0.41) and the extremely low diversity in the eastern group (H _e = 0.25) indicate that this has had an extremely low effective population size and had been isolated from the main group during the last century. Connectivity between the two subpopulations will be of prime importance for the long-term survival of this species in the Cantabrian Mountains.     

Geographic and genetic boundaries of brown bear (Ursus arctos) population in the Caucasus

The taxonomic status of brown bears in the Caucasus remains unclear. Several morphs or subspecies have been identified from the morphological (craniological) data, but the status of each of these subspecies has never been verified by molecular genetic methods. We analysed mitochondrial DNA sequences (control region) to reveal phylogenetic relationships and infer divergence time between brown bear subpopulations in the Caucasus. We estimated migration and gene flow from both mitochondrial DNA and microsatellite allele frequencies, and identified possible barriers to gene flow among the subpopulations. Our suggestion is that all Caucasian bears belong to the nominal subspecies of Ursus arctos. Our results revealed two genetically and geographically distinct maternal haplogroups: one from the Lesser Caucasus and the other one from the Greater Caucasus. The genetic divergence between these haplogroups dates as far back as the beginning of human colonization of the Caucasus. Our analysis of the least‐cost distances between the subpopulations suggests humans as a major barrier to gene flow. The low genetic differentiation inferred from microsatellite allele frequencies indicates that gene flow between the two populations in the Caucasus is maintained through the movements of male brown bears. The Likhi Ridge that connects the Greater and Lesser Caucasus mountains is the most likely corridor for this migration.     

Phylogeography and mitochondrial diversity of extirpated brown bear (Ursus arctos) populations in the contiguous United States and Mexico

The fossil record indicates that the brown bear (Ursus arctos) colonized North America from Asia over 50 000 years ago. The species historically occupied the western United States and northern Mexico but has been extirpated from over 99% of this range in the last two centuries. To evaluate colonization hypotheses, subspecific classifications, and historical patterns and levels of genetic diversity in this region, we sequenced 229 nucleotides of the mitochondrial DNA control region in 108 museum specimens. The work was set in a global context by synthesizing all previous brown bear control region sequences from around the world. In mid-latitude North America a single moderately diverse clade is observed, represented by 23 haplotypes with up to 3.5% divergence. Only eight of 23 haplotypes (35%) are observed in the extensively sampled extant populations suggesting a substantial loss of genetic variability. The restriction of all haplotypes from mid-latitude North America to a single clade suggests that this region was founded by bears with a similar maternal ancestry. However, the levels and distributions of diversity also suggest that the colonizing population was not a small founder event, and that expansion occurred long enough ago for local mutations to accrue. Our data are consistent with recent genetic evidence that brown bears were south of the ice prior to the last glacial maximum. There is no support for previous subspecies designations, although bears of the southwestern United States may have had a distinctive, but recent, pattern of ancestry.     

Seasonal patterns in the physiology of the European brown bear (Ursus arctos arctos) in Finland

The physiological indicators such as body temperature, blood chemistry and hematology of seven European brown bears (Ursus arctos arctos) were used in the present study. They were kept in either the Zoological Garden of University of Oulu (65°N, 25°24'E) or the Ranua Zoological Garden approx. 150 km NE of Oulu. Transmitters with a temperature-dependent pulse rate were implanted subcutaneously or into the abdominal cavity under anesthesia. Our data indicate that the body temperature of the bear decreases during the winter sleep to 4–5°C below the normal level (37.0–37.5°C). The lowest values, 33.1–33.3°C, were measured several times in midwinter. Hematocrit, hemoglobin and erythrocyte counts seem to be higher, and the leucocyte count lower during the denning period than in the awake bear. Plasma N-wastes were lower during the winter sleep than before or after it. The analysed blood parameters showed that plasma catecholamines and thyroid hormones decreased in the fall.     

Nuclear DNA microsatellite analysis of genetic diversity and gene flow in the Scandinavian brown bear (Ursus arctos)

In the 1930s, the Scandinavian brown bear was close to extinction due to vigorous extermination programmes in Norway and Sweden. Increased protection of the brown bear in Scandinavia has resulted in the recovery of four subpopulations, which currently contain close to 1000 individuals. Effective conservation and management of the Scandinavian brown bear requires knowledge of the current levels of genetic diversity and gene flow among the four subpopulations. Earlier studies of mitochondrial DNA (mtDNA) diversity revealed extremely low levels of genetic variation, and population structure that grouped the three northern subpopulations in one genetic clade and the southernmost subpopulation in a second highly divergent clade. In this study, we extended the analysis of genetic diversity and gene flow in the Scandinavian brown bear using data from 19 nuclear DNA microsatellite loci. Results from the nuclear loci were strikingly different than the mtDNA results. Genetic diversity levels in the four subpopulations were equivalent to diversity levels in nonbottlenecked populations from North America, and significantly higher than levels in other bottlenecked and isolated brown bear populations. Gene flow levels between subpopulations ranged from low to moderate and were correlated with geographical distance. The substantial difference in results obtained using mtDNA and nuclear DNA markers stresses the importance of collecting data from both types of genetic markers before interpreting data and making recommendations for the conservation and management of natural populations. Based on the results from the mtDNA and nuclear DNA data sets, we propose one evolutionarily significant unit and four management units for the brown bear in Scandinavia.     

Split parturition observed in a captive North American brown bear (Ursus arctos)

Reproductive physiology in North American ursids is characterized by mating from spring to early summer, delayed implantation, and birth during hibernation. During spring 2008, a captive adult female brown bear was mated with two adult males. Pregnancy was determined by elevated progesterone concentrations during late fall before hibernation. Two male cubs were born on December 31, 2008, and a third female cub was born 17 days later on January 16. All were successfully raised and all were confirmed to have identical paternity. When normalized to age, cub growth rates did not differ. To our knowledge, this is the first documented case of markedly different birth dates in a single litter of brown bear cubs.     

Integrating sign surveys and telemetry data for estimating brown bear (Ursus arctos) density in the Romanian Carpathians

Accurate population size estimates are important information for sustainable wildlife management. The Romanian Carpathians harbor the largest brown bear (Ursus arctos) population in Europe, yet current management relies on estimates of density that lack statistical oversight and ignore uncertainty deriving from track surveys. In this study, we investigate an alternative approach to estimate brown bear density using sign surveys along transects within a novel integration of occupancy models and home range methods. We performed repeated surveys along 2‐km segments of forest roads during three distinct seasons: spring 2011, fall‐winter 2011, and spring 2012, within three game management units and a Natura 2000 site. We estimated bears abundances along transects using the number of unique tracks observed per survey occasion via N‐mixture hierarchical models, which account for imperfect detection. To obtain brown bear densities, we combined these abundances with the effective sampling area of the transects, that is, estimated as a function of the median (± bootstrapped SE) of the core home range (5.58 ± 1.08 km²) based on telemetry data from 17 bears tracked for 1‐month periods overlapping our surveys windows. Our analyses yielded average brown bear densities (and 95% confidence intervals) for the three seasons of: 11.5 (7.8–15.3), 11.3 (7.4–15.2), and 12.4 (8.6–16.3) individuals/100 km². Across game management units, mean densities ranged between 7.5 and 14.8 individuals/100 km². Our method incorporates multiple sources of uncertainty (e.g., effective sampling area, imperfect detection) to estimate brown bear density, but the inference fundamentally relies on unmarked individuals only. While useful as a temporary approach to monitor brown bears, we urge implementing DNA capture–recapture methods regionally to inform brown bear management and recommend increasing resources for GPS collars to improve estimates of effective sampling area.     

Genetic diversity of endangered brown bear (Ursus arctos) populations at the crossroads of Europe, Asia and Africa

Aim  Middle East brown bears ( Ursus arctos syriacus Hemprich and Ehrenberg, 1828) are presently on the edge of extinction. However, little is known of their genetic diversity. This study investigates that question as well as that of Middle East brown bear relationships to surrounding populations of the species.
Location  Middle East region of south-western Asia.
Methods  We performed DNA analyses on 27 brown bear individuals. Twenty ancient bone samples (Late Pleistocene to 20th century) from natural populations and seven present-day samples obtained from captive individuals were analysed.
Results  Phylogenetic analyses of the mitochondrial sequences obtained from seven ancient specimens identify three distinct maternal clades, all unrelated to one recently described from North Africa. Brown bears from Iran exhibit striking diversity (three individuals, three haplotypes) and form a unique clade that cannot be linked to any extant one. Individuals from Syria belong to the Holarctic clade now observed in Eastern Europe, Turkey, Japan and North America. Specimens from Lebanon surprisingly appear as tightly linked to the clade of brown bears now in Western Europe. Moreover, we show that U. a. syriacus in captivity still harbour haplotypes closely linked to those found in ancient individuals.
Main conclusion  This study brings important new information on the genetic diversity of brown bear populations at the crossroads of Europe, Asia and Africa. It reveals a high level of diversity in Middle East brown bears and extends the historical distribution of the Western European clade to the East. Our analyses also suggest the value of a specific breeding programme for captive populations.     

Microevolution of the mitochondrial DNA control region in the Japanese brown bear (Ursus arctos) population.

We investigated nucleotide sequences of the mitochondrial DNA control region to describe natural genetic variations and to assess the relationships between subpopulations of the brown bear Ursus arctos on Hokkaido Island, Japan. Using the polymerase chain reaction product-direct sequencing technique, partial sequences (about 930 bases) of the control region were determined for 56 brown bears sampled throughout Hokkaido Island. A sequence alignment revealed that the brown bear control region included a variable sequence on the 5' side and a repetitive region on the 3' side. Phylogenetic trees reconstructed from the 5' variable region (696-702 bases) exhibited 17 haplotypes, which were clustered into three groups (Clusters A, B, and C). The distribution of each group did not overlap with those of the others, and the three different areas were located in separate mountainous forests of Hokkaido Island. Furthermore, most of the phylogenetically close haplotypes within each group were distributed geographically close to each other. In addition, the 3' repetitive region (arrays of 10 bases) exhibited a much faster mutation rate than the 5' variable region, resulting in heteroplasmy. Such mitochondrial DNA divergence in each group could have occurred after the brown bears migrated from the continent to Hokkaido and became fixed in the different areas.     

Mitogenetic structure of brown bears (Ursus arctos L.) in northeastern Europe and a new time frame for the formation of European brown bear lineages

We estimated the phylogenetic relationships of brown bear maternal haplotypes from countries of northeastern Europe (Estonia, Finland and European Russia), using sequences of mitochondrial DNA (mtDNA) control region of 231 bears. Twenty-five mtDNA haplotypes were identified. The brown bear population in northeastern Europe can be divided into three haplogroups: one with bears from all three countries, one with bears from Finland and Russia, and the third composed almost exclusively of bears from European Russia. Four haplotypes from Finland and European Russia matched exactly with haplotypes from Slovakia, suggesting the significance of the current territory of Slovakia in ancient demographic processes of brown bears. Based on the results of this study and those from the recent literature, we hypothesize that the West Carpathian Mountains have served either as one of the northernmost refuge areas or as an important movement corridor for brown bears of the Eastern lineage towards northern Europe during or after the last ice age. Bayesian analyses were performed to investigate the temporal framework of brown bear lineages in Europe. The molecular clock was calibrated using Beringian brown bear sequences derived from radiocarbon-dated ancient samples, and the estimated mutation rate was 29.8% (13.3%-47.6%) per million years. The whole European population and Western and Eastern lineages formed about 175,000, 70,000 and 25,000 years before present, respectively. Our approach to estimating the time frame of brown bear evolution demonstrates the importance of using an appropriate mutation rate, and this has implications for other studies of Pleistocene populations.     

Multisource noninvasive genetics of brown bears (Ursus arctos) in Greece reveals a highly structured population and a new matrilineal contact zone in southern Europe

In human‐dominated landscapes, connectivity is crucial for maintaining demographically stable mammalian populations. Here, we provide a comprehensive noninvasive genetic study for the brown bear population in the Hellenic Peninsula. We analyze its population structuring and connectivity, estimate its population size throughout its distribution, and describe its phylogeography in detail for the first time. Our results, based on 150 multilocus genotypes and on 244‐bp sequences of the mtDNA control region, show the population is comprised by three highly differentiated genetic clusters, consistent with geographical populations of Pindos, Peristeri, and Rhodope. By detecting two male bears with Rhodopean ancestry in the western demes, we provide strong evidence for the ongoing genetic connectivity of the geographically fragmented eastern and western distributions, which suggests connectivity of the larger East Balkan and Pindos‐Dinara populations. Total effective population size (N _e) was estimated to be 199 individuals, and total combined population size (N _C) was 499, with each cluster showing a relatively high level of genetic variability, suggesting that migration has been sufficient to counteract genetic erosion. The mtNDA results were congruent with the microsatellite data, and the three genetic clusters were matched predominantly with an equal number of mtDNA haplotypes that belong to the brown bear Western mitochondrial lineage (Clade 1), with two haplotypes being globally new and endemic. The detection of a fourth haplotype that belongs to the Eastern lineage (Clade 3a1) in three bears from the western distribution places the southernmost secondary contact zone between the Eastern and Western lineages in Greece and generates new hypotheses about postglacial maxima migration routes. This work indicates that the genetic composition and diversity of Europe''s low‐latitude fringe population are the outcome of ancient and historical events and highlight its importance for the connectivity and long‐term persistence of the species in the Balkans.     

Mitochondrial cytochrome b gene variation in brown bear (Ursus arctos Linnaeus, 1758) from southern part of Russian Far East

The genetic variability of brown bear Ursus arctos from the southern part of the Russian Far East was first examined based on the variations in the mitochondrial DNA cytochrome b sequence. The presence of two phylogenetic groups of haplotypes described previously for other parts of the species range was demonstrated. Part of the samples belonged to the haplotype group distributed across the whole range, while another part belonged to the rare group previously only reported for Japan and Alaska. These findings partially clarify the pattern of brown-bear colonization on the territory of the Russian Far East and Japan.     

The influence of diet on faecal DNA amplification and sex identification in brown bears (Ursus arctos)

To evaluate the influence of diet on faecal DNA amplification, 11 captive brown bears (Ursus arctos) were placed on six restricted diets: grass (Trifolium spp., Haplopappus hirtus and Poa pratensis), alfalfa (Lupinus spp.), carrots (Daucus spp.), white-tailed deer (Odocoileus virginianus), blueberries (Vaccinium spp.) and salmon (Salmo spp.). DNA was extracted from 50 faecal samples of each restricted diet, and amplification of brown bear DNA was attempted for a mitochondrial DNA (mtDNA) locus and nuclear DNA (nDNA) locus. For mtDNA, no significant differences were observed in amplification success rates across diets. For nDNA, amplification success rates for salmon diet extracts were significantly lower than all other diet extracts (P < 0.001). To evaluate the accuracy of faecal DNA sex identification when female carnivores consume male mammalian prey, female bears were fed male white-tailed deer. Four of 10 extracts amplified, and all extracts were incorrectly scored as male due to amplification of X and Y-chromosome fragments. The potential biases highlighted in this study have broad implications for researchers using faecal DNA for individual and sex identification, and should be evaluated in other species.