Evolutionary history of enigmatic bears in the Tibetan Plateau-Himalaya region and the identity of the yeti

Range-wide evolutionary relationships and historical demography of brown bears (Ursus arctos) revealed by whole-genome sequencing of isolated central Asian populations

Phylogeographic studies uncover hidden pathways of divergence and inform conservation. Brown bears (Ursus arctos) have one of the broadest distributions of all land mammals, ranging from Eurasia to North America, and are an important model for evolutionary studies. Although several whole genomes were available for individuals from North America, Europe and Asia, limited whole-genome data were available from Central Asia, including the highly imperilled brown bears in the Gobi Desert. To fill this knowledge gap, we sequenced whole genomes from nine Asian brown bears from the Gobi Desert of Mongolia, Northern Mongolia and the Himalayas of Pakistan. We combined these data with published brown bear sequences from Europe, Asia and North America, as well as other bear species. Our goals were to determine the evolutionary relationships among brown bear populations worldwide, their genetic diversity and their historical demography. Our analyses revealed five major lineages of brown bears based on a filtered set of 684,081 single nucleotide polymorphisms. We found distinct evolutionary lineages of brown bears in the Gobi, Himalayas, northern Mongolia, Europe and North America. The lowest level of genetic diversity and the highest level of inbreeding were found in Pakistan, the Gobi Desert and Central Italy. Furthermore, the effective population size (Ne ) for all brown bears decreased over the last 70,000 years. Our results confirm the genetic distinctiveness and ancient lineage of brown bear subspecies in the Gobi Desert of Mongolia and the Himalayas of Pakistan and highlight their importance for conservation.

Ancient bears provide insights into Pleistocene ice age refugia in Southeast Alaska

During the Late Pleistocene, major parts of North America were periodically covered by ice sheets. However, there are still questions about whether ice-free refugia were present in the Alexander Archipelago along the Southeast (SE) Alaska coast during the last glacial maximum (LGM). Numerous subfossils have been recovered from caves in SE Alaska, including American black (Ursus americanus) and brown (U. arctos) bears, which today are found in the Alexander Archipelago but are genetically distinct from mainland bear populations. Hence, these bear species offer an ideal system to investigate long-term occupation, potential refugial survival and lineage turnover. Here, we present genetic analyses based on 99 new complete mitochondrial genomes from ancient and modern brown and black bears spanning the last ~45,000 years. Black bears form two SE Alaskan subclades, one preglacial and another postglacial, that diverged >100,000 years ago. All postglacial ancient brown bears are closely related to modern brown bears in the archipelago, while a single preglacial brown bear is found in a distantly related clade. A hiatus in the bear subfossil record around the LGM and the deep split of their pre- and postglacial subclades fail to support a hypothesis of continuous occupancy in SE Alaska throughout the LGM for either species. Our results are consistent with an absence of refugia along the SE Alaska coast, but indicate that vegetation quickly expanded after deglaciation, allowing bears to recolonize the area after a short-lived LGM peak.

Arctic introgression and chromatin regulation facilitated rapid Qinghai-Tibet Plateau colonization by an avian predator

The Qinghai-Tibet Plateau (QTP), possesses a climate as cold as that of the Arctic, and also presents uniquely low oxygen concentrations and intense ultraviolet (UV) radiation. QTP animals have adapted to these extreme conditions, but whether they obtained genetic variations from the Arctic during cold adaptation, and how genomic mutations in non-coding regions regulate gene expression under hypoxia and intense UV environment, remain largely unknown. Here, we assemble a high-quality saker falcon genome and resequence populations across Eurasia. We identify female-biased hybridization with Arctic gyrfalcons in the last glacial maximum, that endowed eastern sakers with alleles conveying larger body size and changes in fat metabolism, predisposing their QTP cold adaptation. We discover that QTP hypoxia and UV adaptations mainly involve independent changes in non-coding genomic variants. Our study highlights key roles of gene flow from Arctic relatives during QTP hypothermia adaptation, and cis-regulatory elements during hypoxic response and UV protection.

First camera-trap confirmation of Tibetan Brown Bear Ursus arctos pruinosus Blyth, 1854 (Mammalia: Carnivora: Ursidae) with a review of its distribution and status in Nepal

The Tibetan Brown Bear Ursus arctos pruinosus is a large mammalian carnivore of high-altitude environments that is closely associated with the pastoral landscape. Limited information is available on this species, probably due to its rarity in the Himalaya. To date, scientific evidence of the presence of Tibetan Brown Bears has not been reported officially. The information presented here is based on data collected in the central Himalayan region of Nepal in 2003–2014 during biodiversity surveys and other research. Methods included random walks along livestock trails, transect surveys, opportunistic camera trapping, and herders’ reports & interviews. This is the first camera-trap confirmation of the Tibetan Brown Bear in the central Himalaya. The distribution map was updated based on direct observation, signs and field reports gathered from reliable sources. The presence of signs (diggings, footprints, and feces) and direct observation in the Annapurna-Manaslu landscape reveal that bears are closely associated with Himalayan marmots and other small rodents. Local folklore, legends, and cultural beliefs have played important roles in Brown Bear conservation in the central Himalaya. 

Ancient divergence of Indian and Tibetan wolves revealed by recombination-aware phylogenomics

The grey wolf (Canis lupus) expanded its range across Holarctic regions during the late Pleistocene. Consequently, most grey wolves share recent (<100,000 years ago) maternal origins corresponding to a widespread Holarctic clade. However, two deeply divergent (200,000-700,000 years ago) mitochondrial clades are restricted, respectively, to the Indian subcontinent and the Tibetan Plateau, where remaining wolves are endangered. No genome-wide analysis had previously included wolves corresponding to the mitochondrial Indian clade or attempted to parse gene flow and phylogeny. We sequenced four Indian and two Tibetan wolves and included 31 additional canid genomes to resolve the phylogenomic history of grey wolves. Genomic analyses revealed Indian and Tibetan wolves to be distinct from each other and from broadly distributed wolf populations corresponding to the mitochondrial Holarctic clade. Despite gene flow, which was reflected disproportionately in high-recombination regions of the genome, analyses revealed Indian and Tibetan wolves to be basal to Holarctic grey wolves, in agreement with the mitochondrial phylogeny. In contrast to mitochondrial DNA, however, genomic findings suggest the possibility that the Indian wolf could be basal to the Tibetan wolf, a discordance potentially reflecting selection on the mitochondrial genome. Together, these findings imply that southern regions of Asia have been important centers for grey wolf evolution and that Indian and Tibetan wolves represent evolutionary significant units (ESUs). Further study is needed to assess whether these ESUs warrant recognition as distinct species. This question is especially urgent regarding the Indian wolf, which represents one of the world's most endangered wolf populations.

Mitogenome Phylogeny Including Data from Additional Subspecies Provides New Insights into the Historical Biogeography of the Eurasian lynx Lynx lynx

Previous molecular studies of the wide-ranging Eurasian lynx Lynx lynx focused mainly on its northern Palearctic populations, with the consequence that the reconstruction of this species’ evolutionary history did not include genetic variation present in its southern Palearctic distribution. We sampled a previously not considered Asian subspecies (L. l. dinniki), added published data from another Asian subspecies (L. l. isabellinus), and reassessed the Eurasian lynx mtDNA phylogeny along with previously published data from northern Palearctic populations. Our mitogenome-based analyses revealed the existence of three major clades (A: Central Asia, B: SE Europe/SW Asia, C: Europe and Northern Asia) and at least five lineages, with diversification in Lynx lynx commencing at least 28kyr earlier than hitherto estimated. The subspecies L. l. isabellinus harbors the most basal matriline, consistent with the origin of Lynx lynx in this subspecies’ current range. L. l. dinniki harbors the second most basal matriline, which is related to, and may be the source of, the mtDNA diversity of the critically endangered Balkan lynx L. l. balcanicus. Our results suggest that the Anatolian peninsula was a glacial refugium for Eurasian lynx, with previously unconsidered implications for the colonization of Europe by this species.

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.

A small neighborhood well-organized: seasonal and daily activity patterns of the community of large and mid-sized mammals around waterholes in the Gobi Desert, Mongolia

Background

Animal communities have complex patterns of ecological segregation at different levels according to food resources, habitats, behavior, and activity patterns. Understanding these patterns among the community is essential for the conservation of the whole ecosystem. However, these networks are difficult to study nowadays, due to anthropic disturbances and local extinctions, making it difficult to conclude if segregation patterns are natural or human-induced. We studied ecological segregation in a community of large and mid-sized mammals in the Great Gobi Desert, a remote arid area free from recent extinctions and human disturbances. Activity patterns of 10 sympatric mammal species were monitored around 6 waterholes through camera-trapping over a two-year period, and analyzed them primarily through circular statistics.

Results

Complex patterns of spatial, seasonal, and daily segregation were found. Overlap in seasonal activity was detected in only 3 of the 45 possible pairs of species. Four species used the waterholes all-year-round, while others peaked their activity during different periods. The Bactrian camel showed continuous daily activity, the grey wolf had bimodal activity, and the argali and Siberian ibex were diurnal, while the others had nocturnal peaks during different hours. Daily and spatial overlap were both detected in only 6 of the 45 pairs. Only one species pair (snow leopard and Eurasian lynx) showed an overlap at two levels: seasonal and daily. Climate and moon phase significantly affected the activity of certain species.

Conclusions

Altogether, the results showed complex patterns of ecological segregation at different levels in the use of the key resource in arid environments: waterholes. These results are important for understanding the biology of these species under natural conditions, as well as potential changes in altered ecosystems, and may help to design conservation strategies.

An early dog from southeast Alaska supports a coastal route for the first dog migration into the Americas

The oldest confirmed remains of domestic dogs in North America are from mid-continent archaeological sites dated approximately 9900 calibrated years before present (cal BP). Although this date suggests that dogs may not have arrived alongside the first Native Americans, the timing and routes for the entrance of New World dogs remain uncertain. Here, we present a complete mitochondrial genome of a dog from southeast Alaska, dated to 10 150 ± 260 cal BP. We compared this high-coverage genome with data from modern dog breeds, historical Arctic dogs and American precontact dogs (PCDs) from before European arrival. Our analyses demonstrate that the ancient dog belongs to the PCD lineage, which diverged from Siberian dogs around 16 700 years ago. This timing roughly coincides with the minimum suggested date for the opening of the North Pacific coastal (NPC) route along the Cordilleran Ice Sheet and genetic evidence for the initial peopling of the Americas. This ancient southeast Alaskan dog occupies an early branching position within the PCD clade, indicating it represents a close relative of the earliest PCDs that were brought alongside people migrating from eastern Beringia southward along the NPC to the rest of the Americas. The stable isotope δ¹³C value of this early dog indicates a marine diet, different from the younger mid-continent PCDs' terrestrial diet. Although PCDs were largely replaced by modern European dog breeds, our results indicate that their population decline started approximately 2000 years BP, coinciding with the expansion of Inuit peoples, who are associated with traditional sled-dog culture. Our findings suggest that dogs formed part of the initial human habitation of the New World, and provide insights into their replacement by both Arctic and European lineages.

Assembling mitogenome of Himalayan Black Bear (U. t. laniger) from low depth reads and its application in drawing phylogenetic inferences

The complete mitogenome of Himalayan black bear (Ursus thibetanus laniger) from Indian Himalayan region was assembled following the modified approach of mitochondrial baiting and mapping using the next-generation sequencing reads. The complete mitogenome was of 16,556 bp long, consisted of 37 genes that contained 13 protein-coding genes, 22 tRNAs, 2 rRNAs and 1 control region. The complete base composition was 31.33% A, 15.24% G, 25.45%C, and 27.98%T and gene arrangement was similar to the other sub-species of Asiatic black bear. The relative synonymous codon usage analysis revealed the maximum abundance of Isoleucine, Tyrosine, Leucine and Threonine. The assembled mitogenome of U. t. laniger exhibited 99% similarity with the mitogenomes of Himalayan black bear available from Nepal and Tibetan Plateau-Himalaya region. The findings of the present study has proven low depth sequencing data, adequate and highly efficient in rapid recovering the mitochondrial genome by overcoming the conventional strategies of obtaining long-range PCR and subsequently drawing phylogenetic inferences.

Disentangling the interplay of positive and negative selection forces that shaped mitochondrial genomes of Gammarus pisinnus and Gammarus lacustris

We hypothesized that the mitogenome of Gammarus lacustris (GL), native to the Qinghai-Tibet Plateau, might exhibit genetic adaptations to the extreme environmental conditions associated with high altitudes (greater than 3000 m). To test this, we also sequenced the mitogenome of Gammarus pisinnus (GP), whose native range is close to the Tibetan plateau, but at a much lower altitude (200-1500 m). The two mitogenomes exhibited conserved mitochondrial architecture, but low identity between genes (55% atp8 to 76.1% cox1). Standard (homogeneous) phylogenetic models resolved Gammaridae as paraphyletic, but 'heterogeneous' CAT-GTR model as monophyletic. In indirect support of our working hypothesis, GL, GP and Gammarus fossarum exhibit evidence of episodic diversifying selection within the studied Gammaroidea dataset. The mitogenome of GL generally evolves under a strong purifying selection, whereas GP evolves under directional (especially pronounced in atp8) and/or relaxed selection. This is surprising, as GP does not inhabit a unique ecological niche compared to other gammarids. We propose that this rapid evolution of the GP mitogenome may be a reflection of its relatively recent speciation and heightened non-adaptive (putatively metabolic rate-driven) mutational pressures. To test these hypotheses, we urge sequencing mitogenomes of remaining Gammarus species populating the same geographical range as GP.

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.

Mitochondrial Architecture Rearrangements Produce Asymmetrical Nonadaptive Mutational Pressures That Subvert the Phylogenetic Reconstruction in Isopoda

The phylogeny of Isopoda, a speciose order of crustaceans, remains unresolved, with different data sets (morphological, nuclear, mitochondrial) often producing starkly incongruent phylogenetic hypotheses. We hypothesized that extreme diversity in their life histories might be causing compositional heterogeneity/heterotachy in their mitochondrial genomes, and compromising the phylogenetic reconstruction. We tested the effects of different data sets (mitochondrial, nuclear, nucleotides, amino acids, concatenated genes, individual genes, gene orders), phylogenetic algorithms (assuming data homogeneity, heterogeneity, and heterotachy), and partitioning; and found that almost all of them produced unique topologies. As we also found that mitogenomes of Asellota and two Cymothoida families (Cymothoidae and Corallanidae) possess inversed base (GC) skew patterns in comparison to other isopods, we concluded that inverted skews cause long-branch attraction phylogenetic artifacts between these taxa. These asymmetrical skews are most likely driven by multiple independent inversions of origin of replication (i.e., nonadaptive mutational pressures). Although the PhyloBayes CAT-GTR algorithm managed to attenuate some of these artifacts (and outperform partitioning), mitochondrial data have limited applicability for reconstructing the phylogeny of Isopoda. Regardless of this, our analyses allowed us to propose solutions to some unresolved phylogenetic debates, and support Asellota are the most likely candidate for the basal isopod branch. As our findings show that architectural rearrangements might produce major compositional biases even on relatively short evolutionary timescales, the implications are that proving the suitability of data via composition skew analyses should be a prerequisite for every study that aims to use mitochondrial data for phylogenetic reconstruction, even among closely related taxa.

Evolutionary history and palaeoecology of brown bear in North-East Siberia re-examined using ancient DNA and stable isotopes from skeletal remains

Over 60% of the modern distribution range of brown bears falls within Russia, yet palaeoecological data from the region remain scarce. Complete modern Russian brown bear mitogenomes are abundant in the published literature, yet examples of their ancient counterparts are absent. Similarly, there is only limited stable isotopic data of prehistoric brown bears from the region. We used ancient DNA and stable carbon (δ¹³C) and nitrogen (δ¹⁵N) isotopes retrieved from five Pleistocene Yakutian brown bears (one Middle Pleistocene and four Late Pleistocene), to elucidate the evolutionary history and palaeoecology of the species in the region. We were able to reconstruct the complete mitogenome of one of the Late Pleistocene specimens, but we were unable to assign it to any of the previously published brown bear mitogenome clades. A subsequent analysis of published mtDNA control region sequences, which included sequences of extinct clades from other geographic regions, assigned the ancient Yakutian bear to the extinct clade 3c; a clade previously identified from Late Quaternary specimens from Eastern Beringia and Northern Spain. Our analyses of stable isotopes showed relatively high δ¹⁵N values in the Pleistocene Yakutian brown bears, suggesting a more carnivorous diet than contemporary brown bears from Eastern Beringia.

High genetic diversity and distinct ancient lineage of Asiatic black bears revealed by non-invasive surveys in the Annapurna Conservation Area, Nepal

Asiatic black bears (Ursus thibetanus) have a widespread distribution in mountain landscapes, and are considered vulnerable globally, but are low-priority species for conservation in Nepal. Habitat fragmentation, illegal hunting, and human-bear conflict are the major threats to Asiatic black bears across their global range. Having an adequate level of genetic variation in a population helps with adapting to rapidly changing environments, and thus is important for the long-term health of bear populations. Accordingly, we conducted non-invasive surveys of bear populations in the Annapurna Conservation Area (ACA) to elucidate genetic diversity, genetic structure, and the phylogenetic relationship of Asiatic black bears from this region of Nepal to other subspecies. To assess levels of genetic diversity and population genetic structure, we genotyped eight microsatellite loci using 147 samples, identifying 60 individuals in an area of approximately 525 km². We found that the Asiatic black bear population in the ACA has maintained high levels of genetic diversity (H_E = 0.76) as compared to other bear populations from range countries. We did not detect a signature of population substructure among sampling localities and this suggests that animals are moving freely across the landscape within the ACA. We also detected a moderate population size that may increase with the availability of suitable habitat in the ACA, so bear-related conflict should be addressed to ensure the long-term viability of this expanding bear populations. Primers specific to bears were designed to amplify a 675 bp fragment of the mitochondrial control region from the collected samples. Three haplotypes were observed from the entire conservation area. The complete mitochondrial genome (16,771 bp), the first obtained from wild populations of the Himalayan black bear (U. t. laniger), was also sequenced to resolve the phylogenetic relationships of closely related subspecies of Asiatic black bears. The resulting phylogeny indicated that Himalayan black bear populations in Nepal are evolutionary distinct from other known subspecies of Asiatic black bears.

The complete mitochondrial genome of Cymothoa indica has a highly rearranged gene order and clusters at the very base of the Isopoda clade

As a result of great diversity in life histories and a large number of described species, taxonomic and phylogenetic uncertainty permeates the entire crustacean order of Isopoda. Large molecular datasets capable of providing sufficiently high phylogenetic resolution, such as mitochondrial genomes (mitogenomes), are needed to infer their evolutionary history with confidence, but isopod mitogenomes remain remarkably poorly represented in public databases. We sequenced the complete mitogenome of Cymothoa indica, a species belonging to a family from which no mitochondrial genome was sequenced yet, Cymothoidae. The mitogenome (circular, 14484 bp, A+T = 63.8%) is highly compact, appears to be missing two tRNA genes (trnI and trnE), and exhibits a unique gene order with a large number of rearrangements. High compactness and the existence of palindromes indicate that the mechanism behind these rearrangements might be associated with linearization events in its evolutionary history, similar to those proposed for isopods from the Armadillidium genus (Oniscidea). Isopods might present an important model system to study the proposed discontinuity in the dynamics of mitochondrial genomic architecture evolution. Phylogenetic analyses (Bayesian Inference and Maximum Likelihood) conducted using nucleotide sequences of all mitochondrial genes resolved Oniscidea and Cymothoida suborders as paraphyletic. Cymothoa indica was resolved as a sister group (basal) to all remaining isopods, which challenges the accepted isopod phylogeny, where Cymothoida are the most derived, and Phreatoicidea the most basal isopod group. There is growing evidence that Cymothoida suborder might be split into two evolutionary distant clades, with parasitic species being the most basal split in the Isopoda clade, but a much larger amount of molecular resources carrying a high phylogenetic resolution will be needed to infer the remarkably complex evolutionary history of this group of animals with confidence.