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

Coupled social and ecological change drove the historical extinction of the California grizzly bear (Ursus arctos californicus)

Large carnivores (order Carnivora) are among the world's most threatened mammals due to a confluence of ecological and social forces that have unfolded over centuries. Combining specimens from natural history collections with documents from archival records, we reconstructed the factors surrounding the extinction of the California grizzly bear (Ursus arctos californicus), a once-abundant brown bear subspecies last seen in 1924. Historical documents portrayed California grizzlies as massive hypercarnivores that endangered public safety. Yet, morphological measurements on skulls and teeth generate smaller body size estimates in alignment with extant North American grizzly populations (approx. 200 kg). Stable isotope analysis (δ13C, δ15N) of pelts and bones (n = 57) revealed that grizzlies derived less than 10% of their nutrition from terrestrial animal sources and were therefore largely herbivorous for millennia prior to the first European arrival in this region in 1542. Later colonial land uses, beginning in 1769 with the Mission era, led grizzlies to moderately increase animal protein consumption (up to 26% of diet), but grizzlies still consumed far less livestock than otherwise claimed by contemporary accounts. We show how human activities can provoke short-term behavioural shifts, such as heightened levels of carnivory, that in turn can lead to exaggerated predation narratives and incentivize persecution, triggering rapid loss of an otherwise widespread and ecologically flexible animal.

The taxonomic basis of subspecies listed as threatened and endangered under the endangered species act

More than 170 subspecies are listed as threatened or endangered under the US Endangered Species Act. Most of these subspecies were described decades ago on the basis of geographical variation in morphology using relatively primitive taxonomic methods. The US Fish and Wildlife Service defaults to subspecies descriptions by taxonomists working with specific groups of organisms, but there is no single definition of subspecies across plants and animals. Valid tests today usually entail molecular analyses of variation within and among populations, although there is no reason that behavioral, ecological or molecular characters could not be used, and include tests for significant differences between samples of the putative endangered subspecies and its nearest geographic relatives. We evaluated data gathered since subspecies listed under the ESA were described finding about one-third are valid (distinct evolutionary taxa), one-third are not, and one-third have not been tested. Therefore, it should not be assumed that because a subspecies occurs in a checklist, it is taxonomically valid. If the US Fish and Wildlife Service intends to continue listing subspecies, we suggest that they convene taxonomic experts representing various groups of organisms to provide a minimal set of criteria for a subspecies to be listed under the ESA.

Phylogeography of ancient and modern brown bears from eastern Eurasia

The brown bear (Ursus arctos) is an iconic carnivoran species of the Northern Hemisphere. Its population history has been studied extensively using mitochondrial markers, which demonstrated signatures of multiple waves of migration, arguably connected with glaciation periods. Among Eurasian brown bears, Siberian populations remain understudied. We have sequenced complete mitochondrial genomes of four ancient (~4.5-40 kya) bears from South Siberia and 19 modern bears from South Siberia and the Russian Far East. Reconstruction of phylogenetic relationships between haplotypes and evaluation of modern population structure have demonstrated that all the studied samples belong to the most widespread Eurasian clade 3. One of the ancient haplotypes takes a basal position relative to the whole of clade 3; the second is basal to the haplogroup 3a (the most common subclade), and two others belong to clades 3a1 and 3b. Modern Siberian bears retain at least some of this diversity; apart from the most common haplogroup 3a, we demonstrate the presence of clade 3b, which was previously found mainly in mainland Eurasia and Northern Japan. Our findings highlight the importance of South Siberia as a refugium for northern Eurasian brown bears and further corroborate the hypothesis of several waves of migration in the Pleistocene.

Ancient DNA reveals multiple origins and migration waves of extinct Japanese brown bear lineages

Little is known about how mammalian biogeography on islands was affected by sea-level fluctuations. In the Japanese Archipelago, brown bears (Ursus arctos) currently inhabit only Hokkaido, the northern island, but Pleistocene fossils indicate a past distribution throughout Honshu, Japan's largest island. However, the difficulty of recovering ancient DNA from fossils in temperate East Asia has limited our understanding of their evolutionary history. Here, we analysed mitochondrial DNA from a 32 500-year-old brown bear fossil from Honshu. Our results show that this individual belonged to a previously unknown lineage that split approximately 160 Ka from its sister lineage, the southern Hokkaido clade. This divergence time and fossil record suggest that brown bears migrated from the Eurasian continent to Honshu at least twice; the first population was an early-diverging lineage (greater than 340 Ka), and the second migrated via Hokkaido after approximately 160 Ka, during the ice age. Thus, glacial-age sea-level falls might have facilitated migrations of large mammals more frequently than previously thought, which may have had a substantial impact on ecosystem dynamics in these isolated islands.

Genetic assessment of captive red pandas (Ailurus fulgens) in American zoos to address management separation by putative subspecies

Red pandas (Ailurus fulgens) are small charismatic mammals native across montane southern Asia, now endangered by human impacts. They are considered "living fossils" as the sole member of a distinct family, warranting higher conservation priority. Therefore, ex situ breeding programs were initiated to prevent extinction and act as genetic reservoirs for reintroduction, although complicated by apparent taxonomic subdivision. This study investigated whether the separation of captive red pandas in the North American Species Survival Plan® by putative subspecies was justified. A 383-bp segment of mitochondrial DNA control region was therefore sequenced from 67 members representing matriarchal lineages of both groups. A network analysis placed the 11 haplotypes found into separate but closely connected clusters, with one group more strongly related than the other. Statistical analyses and diversity indices corroborated differentiation between the two management units. Phylogenetic analyses employing multiple outgroups confirmed, although not robustly, reciprocal monophyly of the four- and seven-haplotype clades representing putative subspecies Ailurus fulgens fulgens and Ailurus fulgens styani, respectively. These empirical results are adequate to justify continued independent management of these zoo subpopulations, but cannot be definitive for taxonomic classification due to limited sampling from their native range. They will, however, be useful in evaluating long-term genetic diversity changes, focusing management efforts on newly revealed evolutionary limitations, and comparing with an assessment of wild red pandas to determine how representative zoo populations are for reintroduction purposes. Maintaining genetic diversity and population structure of endangered species is essential to protect evolutionary potential and adaptations for long-term sustainability.

The Functional Impact of Mitochondrial Structure Across Subcellular Scales

Mitochondria are key determinants of cellular health. However, the functional role of mitochondria varies from cell to cell depending on the relative demands for energy distribution, metabolite biosynthesis, and/or signaling. In order to support the specific needs of different cell types, mitochondrial functional capacity can be optimized in part by modulating mitochondrial structure across several different spatial scales. Here we discuss the functional implications of altering mitochondrial structure with an emphasis on the physiological trade-offs associated with different mitochondrial configurations. Within a mitochondrion, increasing the amount of cristae in the inner membrane improves capacity for energy conversion and free radical-mediated signaling but may come at the expense of matrix space where enzymes critical for metabolite biosynthesis and signaling reside. Electrically isolating individual cristae could provide a protective mechanism to limit the spread of dysfunction within a mitochondrion but may also slow the response time to an increase in cellular energy demand. For individual mitochondria, those with relatively greater surface areas can facilitate interactions with the cytosol or other organelles but may be more costly to remove through mitophagy due to the need for larger phagophore membranes. At the network scale, a large, stable mitochondrial reticulum can provide a structural pathway for energy distribution and communication across long distances yet also enable rapid spreading of localized dysfunction. Highly dynamic mitochondrial networks allow for frequent content mixing and communication but require constant cellular remodeling to accommodate the movement of mitochondria. The formation of contact sites between mitochondria and several other organelles provides a mechanism for specialized communication and direct content transfer between organelles. However, increasing the number of contact sites between mitochondria and any given organelle reduces the mitochondrial surface area available for contact sites with other organelles as well as for metabolite exchange with cytosol. Though the precise mechanisms guiding the coordinated multi-scale mitochondrial configurations observed in different cell types have yet to be elucidated, it is clear that mitochondrial structure is tailored at every level to optimize mitochondrial function to meet specific cellular demands.

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 differentiation and asymmetric gene flow among Carpathian brown bear (Ursus arctos) populations-Implications for conservation of transboundary populations

The abundance and distribution of large carnivores in Europe have been historically reduced. Their recovery requires multilevel coordination, especially regarding transboundary populations. Here, we apply nuclear and mitochondrial genetic markers to test for admixture level and its impact on population genetic structure of contemporary brown bears (Ursus arctos) from the Eastern, Southern, and Western Carpathians. Carpathian Mountains (Europe). Nearly 400 noninvasive brown bear DNA samples from the Western (Poland) and Eastern Carpathians (Bieszczady Mountains in Poland, Slovakia, Ukraine) were collected. Together with DNA isolates from Slovakia and Romania, they were analyzed using the set of eight microsatellite loci and two mtDNA regions (control region and cytochrome b). A set of 113 individuals with complete genotypes was used to investigate genetic differentiation across national boundaries, genetic structuring within and between populations, and movement between populations. Transboundary brown bear subpopulations (Slovakia and Poland) did not show significant internal genetic structure, and thus were treated as cohesive units. All brown bears from the Western Carpathians carried mitochondrial haplotypes from the Eastern lineage, while the Western lineage prevailed in the brown bears from the Bieszczady Mountains. Despite similar levels of microsatellite variability, we documented significant differentiation among the studied populations for nuclear markers and mtDNA. We also detected male-biased and asymmetrical movement into the Bieszczady Mountains population from the Western Carpathians. Our findings suggest initial colonization of the Western Carpathians by brown bears possessing mtDNA from the Eastern lineage. Genetic structuring among populations at microsatellite loci could be a result of human-mediated alterations. Detected asymmetric gene flow suggests ongoing expansion from more abundant populations into the Bieszczady Mountains and thus supports a metapopulation model. The knowledge concerning this complex pattern can be implemented in a joint Carpathian brown bear management plan that should allow population mixing by dispersing males.

The Genome of the North American Brown Bear or Grizzly: Ursus arctos ssp. horribilis

The grizzly bear (Ursus arctos ssp. horribilis) represents the largest population of brown bears in North America. Its genome was sequenced using a microfluidic partitioning library construction technique, and these data were supplemented with sequencing from a nanopore-based long read platform. The final assembly was 2.33 Gb with a scaffold N50 of 36.7 Mb, and the genome is of comparable size to that of its close relative the polar bear (2.30 Gb). An analysis using 4104 highly conserved mammalian genes indicated that 96.1% were found to be complete within the assembly. An automated annotation of the genome identified 19,848 protein coding genes. Our study shows that the combination of the two sequencing modalities that we used is sufficient for the construction of highly contiguous reference quality mammalian genomes. The assembled genome sequence and the supporting raw sequence reads are available from the NCBI (National Center for Biotechnology Information) under the bioproject identifier PRJNA493656, and the assembly described in this paper is version QXTK01000000.

Howling from the past: historical phylogeography and diversity losses in European grey wolves

Genetic bottlenecks resulting from human-induced population declines make alarming symbols for the irreversible loss of our natural legacy worldwide. The grey wolf (Canis lupus) is an iconic example of extreme declines driven by anthropogenic factors. Here, we assessed the genetic signatures of 150 years of wolf persecution throughout the Western Palaearctic by high-throughput mitochondrial DNA sequencing of historical specimens in an unprecedented spatio-temporal framework. Despite Late Pleistocene bottlenecks, we show that historical genetic variation had remained high throughout Europe until the last several hundred years. In Western Europe, where wolves nearly got fully exterminated, diversity dramatically collapsed at the turn of the twentieth century and recolonization from few homogeneous relict populations induced drastic shifts of genetic composition. By contrast, little genetic displacement and steady levels of diversity were maintained in Eastern European regions, where human persecution had lesser effects on wolf demography. By comparing prehistoric, historic and modern patterns of genetic diversity, our study hence traces the timeframe and the active human role in the decline of the grey wolf, an emblematic yet controversial animal which symbolizes the complex relationship between human societies and nature conservation.

Consequences of a demographic bottleneck on genetic structure and variation in the Scandinavian brown bear

The Scandinavian brown bear went through a major decline in population size approximately 100 years ago, due to intense hunting. After being protected, the population subsequently recovered and today numbers in the thousands. The genetic diversity in the contemporary population has been investigated in considerable detail, and it has been shown that the population consists of several subpopulations that display relatively high levels of genetic variation. However, previous studies have been unable to resolve the degree to which the demographic bottleneck impacted the contemporary genetic structure and diversity. In this study, we used mitochondrial and microsatellite DNA markers from pre- and postbottleneck Scandinavian brown bear samples to investigate the effect of the bottleneck. Simulation and multivariate analysis suggested the same genetic structure for the historical and modern samples, which are clustered into three subpopulations in southern, central and northern Scandinavia. However, the southern subpopulation appears to have gone through a marked change in allele frequencies. When comparing the mitochondrial DNA diversity in the whole population, we found a major decline in haplotype numbers across the bottleneck. However, the loss of autosomal genetic diversity was less pronounced, although a significant decline in allelic richness was observed in the southern subpopulation. Approximate Bayesian computations provided clear support for a decline in effective population size during the bottleneck, in both the southern and northern subpopulations. These results have implications for the future management of the Scandinavian brown bear because they indicate a recent loss in genetic diversity and also that the current genetic structure may have been caused by historical ecological processes rather than recent anthropogenic persecution.

The use of museum specimens with high-throughput DNA sequencers

Natural history collections have long been used by morphologists, anatomists, and taxonomists to probe the evolutionary process and describe biological diversity. These biological archives also offer great opportunities for genetic research in taxonomy, conservation, systematics, and population biology. They allow assays of past populations, including those of extinct species, giving context to present patterns of genetic variation and direct measures of evolutionary processes. Despite this potential, museum specimens are difficult to work with because natural postmortem processes and preservation methods fragment and damage DNA. These problems have restricted geneticists' ability to use natural history collections primarily by limiting how much of the genome can be surveyed. Recent advances in DNA sequencing technology, however, have radically changed this, making truly genomic studies from museum specimens possible. We review the opportunities and drawbacks of the use of museum specimens, and suggest how to best execute projects when incorporating such samples. Several high-throughput (HT) sequencing methodologies, including whole genome shotgun sequencing, sequence capture, and restriction digests (demonstrated here), can be used with archived biomaterials.

Revealing the maternal demographic history of Panthera leo using ancient DNA and a spatially explicit genealogical analysis

Background

Understanding the demographic history of a population is critical to conservation and to our broader understanding of evolutionary processes. For many tropical large mammals, however, this aim is confounded by the absence of fossil material and by the misleading signal obtained from genetic data of recently fragmented and isolated populations. This is particularly true for the lion which as a consequence of millennia of human persecution, has large gaps in its natural distribution and several recently extinct populations.

Results

We sequenced mitochondrial DNA from museum-preserved individuals, including the extinct Barbary lion (Panthera leo leo) and Iranian lion (P. l. persica), as well as lions from West and Central Africa. We added these to a broader sample of lion sequences, resulting in a data set spanning the historical range of lions. Our Bayesian phylogeographical analyses provide evidence for highly supported, reciprocally monophyletic lion clades. Using a molecular clock, we estimated that recent lion lineages began to diverge in the Late Pleistocene. Expanding equatorial rainforest probably separated lions in South and East Africa from other populations. West African lions then expanded into Central Africa during periods of rainforest contraction. Lastly, we found evidence of two separate incursions into Asia from North Africa, first into India and later into the Middle East.

Conclusions

We have identified deep, well-supported splits within the mitochondrial phylogeny of African lions, arguing for recognition of some regional populations as worthy of independent conservation. More morphological and nuclear DNA data are now needed to test these subdivisions.

Molecular phylogeography of the brown bear (Ursus arctos) in Northeastern Asia based on analyses of complete mitochondrial DNA sequences

To further elucidate the migration history of the brown bears (Ursus arctos) on Hokkaido Island, Japan, we analyzed the complete mitochondrial DNA (mtDNA) sequences of 35 brown bears from Hokkaido, the southern Kuril Islands (Etorofu and Kunashiri), Sakhalin Island, and the Eurasian Continent (continental Russia, Bulgaria, and Tibet), and those of four polar bears. Based on these sequences, we reconstructed the maternal phylogeny of the brown bear and estimated divergence times to investigate the timing of brown bear migrations, especially in northeastern Eurasia. Our gene tree showed the mtDNA haplotypes of all 73 brown and polar bears to be divided into eight divergent lineages. The brown bear on Hokkaido was divided into three lineages (central, eastern, and southern). The Sakhalin brown bear grouped with eastern European and western Alaskan brown bears. Etorofu and Kunashiri brown bears were closely related to eastern Hokkaido brown bears and could have diverged from the eastern Hokkaido lineage after formation of the channel between Hokkaido and the southern Kuril Islands. Tibetan brown bears diverged early in the eastern lineage. Southern Hokkaido brown bears were closely related to North American brown bears.

Impact of Quaternary climatic changes and interspecific competition on the demographic history of a highly mobile generalist carnivore, the coyote

Recurrent cycles of climatic change during the Quaternary period have dramatically affected the population genetic structure of many species. We reconstruct the recent demographic history of the coyote (Canis latrans) through the use of Bayesian techniques to examine the effects of Late Quaternary climatic perturbations on the genetic structure of a highly mobile generalist species. Our analysis reveals a lack of phylogeographic structure throughout the range but past population size changes correlated with climatic changes. We conclude that even generalist carnivorous species are very susceptible to environmental changes associated with climatic perturbations. This effect may be enhanced in coyotes by interspecific competition with larger carnivores.

Nondestructive DNA extraction from museum specimens

Natural history museums around the world hold millions of animal and plant specimens that are potentially amenable to genetic analyses. With more and more populations and species becoming extinct, the importance of these specimens for phylogenetic and phylogeographic analyses is rapidly increasing. However, as most DNA extraction methods damage the specimens, nondestructive extraction methods are useful to balance the demands of molecular biologists, morphologists, and museum curators. Here, I describe a method for nondestructive DNA extraction from bony specimens (i.e., bones and teeth). In this method, the specimens are soaked in extraction buffer, and DNA is then purified from the soaking solution using adsorption to silica. The method reliably yields mitochondrial and often also nuclear DNA. The method has been adapted to DNA extraction from other types of specimens such as arthropods.

A rapid column-based ancient DNA extraction method for increased sample throughput

Genetic analyses using museum specimens and ancient DNA from fossil samples are becoming increasingly important in phylogenetic and especially population genetic studies. Recent progress in ancient DNA sequencing technologies has substantially increased DNA sequence yields and, in combination with barcoding methods, has enabled large-scale studies using any type of DNA. Moreover, more and more studies now use nuclear DNA sequences in addition to mitochondrial ones. Unfortunately, nuclear DNA is, due to its much lower copy number in living cells compared to mitochondrial DNA, much more difficult to obtain from low-quality samples. Therefore, a DNA extraction method that optimizes DNA yields from low-quality samples and at the same time allows processing many samples within a short time frame is immediately required. In fact, the major bottleneck in the analysis process using samples containing low amounts of degraded DNA now lies in the extraction of samples, as column-based methods using commercial kits are fast but have proven to give very low yields, while more efficient methods are generally very time-consuming. Here, we present a method that combines the high DNA yield of batch-based silica extraction with the time-efficiency of column-based methods. Our results on Pleistocene cave bear samples show that DNA yields are quantitatively comparable, and in fact even slightly better than with silica batch extraction, while at the same time the number of samples that can conveniently be processed in parallel increases and both bench time and costs decrease using this method. Thus, this method is suited for harvesting the power of high-throughput sequencing using the DNA preserved in the millions of paleontological and museums specimens.

The influence of climatic oscillations during the Quaternary Era on the genetic structure of Asian black bears in Japan

The Asian black bear (Ursus thibetanus) inhabits two of the main islands, Honshu and Shikoku, in Japan. To determine how climatic oscillations during the Quaternary Era affected the genetic structure of the black bear populations in Japan, we examined their phylogeographic relationships and compared their genetic structure. We analysed an approximately 700-bp sequence in the D-loop region of mitochondrial DNA collected from 589 bears in this study with 108 bears from a previous study. We observed a total of 57 haplotypes and categorized them into three clusters (Eastern, Western and Southern) based on the spatial distribution of the haplotypes. All but 2 of the 41 haplotypes in the Eastern cluster were distributed locally. Genetic diversity was generally low in northern Japan and high in central Japan. Demographic tests rejected the expansion model in northern populations. Haplotypes of the Western and Southern clusters were unique to local populations. We conclude that the extant genetic structure of the Asian black bear populations arose as follows: first, populations became small and genetic drift decreased genetic diversity in the northern area during the last glacial period, whereas large continuous populations existed in the southern part of central Japan. These patterns were essentially maintained until the present time. In western and southern Japan, the effects of climatic oscillations were smaller, and thus, local structure was maintained.

Unequal rates of Y chromosome gene divergence during speciation of the family Ursidae

Evolution of the bear family Ursidae is well investigated in terms of morphological, paleontological, and genetic features. However, several phylogenetic ambiguities occur within the subfamily Ursinae (the family Ursidae excluding the giant panda and spectacled bear), which may correlate with behavioral traits of female philopatry and male-biased dispersal which form the basis of the observed matriarchal population structure in these species. In the process of bear evolution, we investigate the premise that such behavioral traits may be reflected in patterns of variation among genes with different modes of inheritance: matrilineal mitochondrial DNA (mtDNA), patrilineal Y chromosome, biparentally inherited autosomes, and the X chromosome. In the present study, we sequenced 3 Y-linked genes (3,453 bp) and 4 X-linked genes (4,960 bp) and reanalyzed previously published sequences from autosome genes (2,347 bp) in ursid species to investigate differences in evolutionary rates associated with patterns of inheritance. The results describe topological incongruence between sex-linked genes and autosome genes and between nuclear DNA and mtDNA. In more ancestral branches within the bear phylogeny, Y-linked genes evolved faster than autosome and X-linked genes, consistent with expectations based on male-driven evolution. However, this pattern changes among branches leading to each species within the lineage of Ursinae whereby the evolutionary rates of Y-linked genes have fewer than expected substitutions. This inconsistency between more recent nodes of the bear phylogeny with more ancestral nodes may reflect the influences of sex-biased dispersal as well as molecular evolutionary characteristics of the Y chromosome, and stochastic events in species natural history, and phylogeography unique to ursine bears.