Diversity of mitochondrial D-loop haplotypes from ancient Thracian horses in Bulgaria.

The domestication of the horse began possibly more than 5000 years ago in the western part of the Eurasian steppe, and according to the leading hypothesis, horses first spread from the Steppe toward the region of the Thracian culture, starting in the second half of the 2nd millennium BCE and flourished from the fifth to first centuries BCE, mainly located in present-day Bulgaria. We analyzed 17 horse bone remains excavated from Thracian archaeological sites (fourth to first centuries BCE) in Bulgaria and successfully identified 17 sequences representing 14 different haplotypes of the mitochondrial D-loop. Compared with the mtDNA haplotypes of modern horses around the world, ancient Thracian horses in Bulgaria are thought to be more closely related to modern horses of Southern Europe and less related to those of Central Asia. In addition, the haplotypes we obtained represented 11 previously reported modern horse mtDNA haplogroups: A, B, D, E, G, H, I, L, N, P, and Q. All the haplogroups contain modern and regionally predominant haplotypes occurring in Europe, the Middle East, and Central Asia. Our results indicate that Thracian horses in Bulgaria have had relatively high genetic diversity and are closely related to modern horse breeds.

Diversity of the MHC class II DRB gene in the wolverine (Carnivora: Mustelidae: Gulo gulo) in Finland.

The wolverine (Gulo gulo) in Finland has undergone significant population declines in the past. Since major histocompatibility complex (MHC) genes encode proteins involved in pathogen recognition, the diversity of these genes provides insights into the immunological fitness of regional populations. We sequenced 862 amplicons (242 bp) of MHC class II DRB exon 2 from 32 Finnish wolverines and identified 11 functional alleles and three pseudogenes. A molecular phylogenetic analysis indicated trans-species polymorphism, and PAML and MEME analyses indicated positive selection, suggesting that the Finnish wolverine DRB genes have evolved under balancing and positive selection. In contrast to DRB gene analyses in other species, allele frequencies in the Finnish wolverines clearly indicated the existence of two regional subpopulations, congruent with previous studies based on neutral genetic markers. In the Finnish wolverine, rapid population declines in the past have promoted genetic drift, resulting in a lower genetic diversity of DRB loci, including fewer alleles and positively selected sites, than other mustelid species analyzed previously. Our data suggest that the MHC region in the Finnish wolverine population was likely affected by a recent bottleneck.

Population genetic structure and diversity of the East Balkan Swine (Sus scrofa) in Bulgaria, revealed by mitochondrial DNA and microsatellite analyses.

The East Balkan Swine (EBS) is the only indigenous pig breed in Bulgaria. We analyzed the mitochondrial DNA (mtDNA) control region and 21 microsatellite loci for 198 individuals from 11 farms in Bulgaria. Obtained 11 mtDNA haplotypes including three novel ones were grouped to two major clades, European clade E1 (146/198 individuals, 73.7%) and Asian clade A (52/198, 26.3%). The mixture of the two clades may have resulted from historical crossbreeding between the European and Asian pig breeds. Clade A was frequent in southeastern Bulgaria (Burgas Province), but less frequent or absent in northeastern Bulgaria (Varna and Shumen Provinces). The distribution of Europe- and Asia-specific haplotypes relative to EBS farm locations could be attributed to regional differences of breeding systems (e.g., crossbreeding with imported commercial pigs). A microsatellite analysis showed high heterozygosities for all the EBS farms, and negative inbreeding coefficients presumably due to crossing with commercial pigs or wild boars and/or efforts to reduce inbreeding by farmers. Bayesian clustering analyses showed that all farm populations are genetically well distinguishable from one another. Although diversity has been maintained by the efforts of farmers and a breeding association, the effective population size remains small, and conservation efforts should be continued.

Demographic History of the Brown Bear (Ursus arctos) on Hokkaido Island, Japan, Based on Whole-Genomic Sequence Analysis.

Previous studies of the brown bear (Ursus arctos) on Hokkaido Island, Japan, have detected three geographically distinct subpopulations representing different mitochondrial lineages and shown that gene flow between subpopulations has occurred due to male-biased dispersal. In this study, we determined whole-genomic sequences for six Hokkaido brown bears and analyzed these data along with previously published genomic sequences of 17 brown bears from other parts of the world. We found that the Hokkaido population is genetically distinct from the other populations, keeping genetic diversity higher than the endangered populations in western Europe but lower than most populations on the continents. A reconstruction of historical demography showed no increase in population size for the Hokkaido population during the Eemian interglacial period (130,000-114,000 years ago). In a phylogenetic analysis of the autosomal data, the Hokkaido population formed a clade distinct from North American and European populations, showing that it has maintained genetic diversity independently from continental populations following geographical isolation on the island. This autosomal genetic similarity contrasts with the geographically separate mitochondrial lineages on Hokkaido and indicates the occurrence of male-driven gene flow between subpopulations.

Evolution of MHC class I genes in Japanese and Russian raccoon dogs, Nyctereutes procyonoides (Carnivora: Canidae).

Major histocompatibility complex (MHC) genes have been widely studied to assess the immunological fitness and evolutionary adaptation of animal populations. Among the Canidae, the raccoon dog's adventurous nature, omnivorous behavior, and high variability of intracellular pathogens make it ideal to study selection on MHC class I in a non-model canid species. Here, we examined allelic diversity and evolutionary patterns of MHC class I genes in the raccoon dog (Nyctereutes procyonoides). We identified 48 novel MHC class I alleles from 31 raccoon dogs from Japan and Russia. Some alleles were geographically restricted, whereas others were widely distributed across the species' range. The rate of non-synonymous substitutions was greater than that of synonymous substitutions for both exon 2 and exon 3 encoding α1 and α2 domains, respectively, in the α chain of the MHC class I protein. Positively selected sites at the amino acid level were evident in both the α1 and α2 domains, and a recombination breakpoint was found in exon 3. Bayesian phylogenetic trees showed no evidence of trans-species polymorphism (TSP) with alleles from carnivoran species in other families but did detect TSP between raccoon dogs and the domestic dog, Canis familiaris, indicative of long-term balancing selection in canids. Our results indicate that the extensive allelic diversity of MHC class I in Japanese and Russian raccoon dogs has been influenced and maintained by pathogen-driven positive selection, recombination, and long-term balancing selection. Supplementary Information: The online version contains supplementary material available at 10.1007/s13364-021-00561-y.

Evolution of MHC class I genes in Eurasian badgers, genus Meles (Carnivora, Mustelidae).

Because of their role in immune defense against pathogens, major histocompatibility complex (MHC) genes are useful in evolutionary studies on how wild vertebrates adapt to their environments. We investigated the molecular evolution of MHC class I (MHCI) genes in four closely related species of Eurasian badgers, genus Meles. All four species of badgers showed similarly high variation in MHCI sequences compared to other Carnivora. We identified 7-21 putatively functional MHCI sequences in each of the badger species, and 2-7 sequences per individual, indicating the existence of 1-4 loci. MHCI exon 2 and 3 sequences encoding domains α1 and α2 exhibited different clade topologies in phylogenetic networks. Non-synonymous nucleotide substitutions at codons for antigen-binding sites exceeded synonymous substitutions for domain α1 but not for domain α2, suggesting that the domains α1 and α2 likely had different evolutionary histories in these species. Positive selection and recombination seem to have shaped the variation in domain α2, whereas positive selection was dominant in shaping the variation in domain α1. In the separate phylogenetic analyses for exon 2, exon 3, and intron 2, each showed three clades of Meles alleles, with rampant trans-species polymorphism, indicative of the long-term maintenance of ancestral MHCI polymorphism by balancing selection.

Genetic Diversity of MHC Class II DRB1 Exon 2 in the Red Fox (Vulpes vulpes) on Hokkaido, Japan.

To assess the genetic diversity of the red fox (Vulpes vulpes) population on Hokkaido Island, northern Japan, we examined genotypes of exon 2 of the major histocompatibility complex (MHC) class II DRB1 gene for 232 individuals and identified 17 novel alleles. The subpopulation in the Southern area was genetically differentiated from those in all other areas on Hokkaido, suggesting isolation in southern Hokkaido. In addition, the Southern subpopulation is lower in genetic diversity than the other subpopulations, possibly resulting from purifying selection and/or a recent bottleneck. The non-synonymous substitutions exceeded the synonymous substitutions for codons encoding antigen-binding sites (ABSs) in exon 2, indicating that the red fox DRB1 alleles have evolved under positive selection. In a Bayesian phylogenetic tree, although most of the DRB1 alleles from the Hokkaido red fox were contained within a red fox-like canid (Vulpes) clade, some belonged to another canid clade. This means trans-species polymorphism maintained by balancing selection. Our results showed the local variability and the presence of selection on the MHC gene in this population, which contributes to the understanding of the historical background and the molecular evolution.

Cranial Variation in the Siberian Weasel Mustela sibirica (Carnivora, Mustelidae) and its Possible Taxonomic Implications.

Alexei V. Abramov, Andrey Yu. Puzachenko, and Ryuichi Masuda (2018) Morphometric variation in 23 cranial characters of 555 Siberian weasels (Mustela sibirica) was studied across its whole distribution range. Most of the distribution range in Siberia and China is occupied by medium-sized weasels, whereas the eastern part of the species range - including the Russian Far East, Korea and eastern China - is occupied by the larger form. Specimens from the Pacific islands (Jeju and Tsushima) were morphologically closely related to the western form of M. sibirica than to the neighboring continental weasels. The western form can be treated as nominotypical subspecies M. s. sibirica Pallas, 1773, whereas the eastern form can be treated as M. s. manchurica Brass, 1911. Small-sized weasels from the eastern Himalayan area (Myanmar and southwestern China) form a distinct group within M. sibirica, and they were treated as a subspecies; M. s. moupinensis (Milne- Edwards, 1874). Specimens from the western Himalayas (Kashmir, Nepal and Sikkim) are morphologically distinct from all other populations of Mustela sibirica and can be treated as a separate species Mustela subhemachalana Hodgson, 1837.

Low genetic variation in the MHC class II DRB gene and MHC-linked microsatellites in endangered island populations of the leopard cat (Prionailurus bengalensis) in Japan.

Isolated populations of the leopard cat (Prionailurus bengalensis) on Tsushima and Iriomote islands in Japan are classified as subspecies P. b. euptilurus and P. b. iriomotensis, respectively. Because both populations have decreased to roughly 100, an understanding of their genetic diversity is essential for conservation. We genotyped MHC class II DRB exon 2 and MHC-linked microsatellite loci to evaluate the diversity of MHC genes in the Tsushima and Iriomote cat populations. We detected ten and four DRB alleles in these populations, respectively. A phylogenetic analysis showed DRB alleles from both populations to be closely related to those in other felid DRB lineages, indicating trans-species polymorphism. The MHC-linked microsatellites were more polymorphic in the Tsushima than in the Iriomote population. The MHC diversity of both leopard cat populations is much lower than in the domestic cat populations on these islands, probably due to inbreeding associated with founder effects, geographical isolation, or genetic drift. Our results predict low resistance of the two endangered populations to new pathogens introduced to the islands.

Paternal phylogeographic structure of the brown bear (Ursus arctos) in northeastern Asia and the effect of male-mediated gene flow to insular populations.

BACKGROUND: Sex-biased dispersal is widespread among mammals, including the brown bear (Ursus arctos). Previous phylogeographic studies of the brown bear based on maternally inherited mitochondrial DNA have shown intraspecific genetic structuring around the northern hemisphere. The brown bears on Hokkaido Island, northern Japan, comprise three distinct maternal lineages that presumably immigrated to the island from the continent in three different periods. Here, we investigate the paternal genetic structure across northeastern Asia and assess the connectivity among and within intraspecific populations in terms of male-mediated gene flow. RESULTS: We analyzed paternally inherited Y-chromosomal DNA sequence data and Y-linked microsatellite data of 124 brown bears from Hokkaido, the southern Kuril Islands (Kunashiri and Etorofu), Sakhalin, and continental Eurasia (Kamchatka Peninsula, Ural Mountains, European Russia, and Tibet). The Hokkaido brown bear population is paternally differentiated from, and lacked recent genetic connectivity with, the continental Eurasian and North American populations. We detected weak spatial genetic structuring of the paternal lineages on Hokkaido, which may have arisen through male-mediated gene flow among natal populations. In addition, our results suggest that the different dispersal patterns between male and female brown bears, combined with the founder effect and subsequent genetic drift, contributed to the makeup of the Etorofu Island population, in which the maternal and paternal lineages show different origins. CONCLUSIONS: Brown bears on Hokkaido and the adjacent southern Kuril Islands experienced different maternal and paternal evolutionary histories. Our results indicate that sex-biased dispersal has played a significant role in the evolutionary history of the brown bear in continental populations and in peripheral insular populations, such as on Hokkaido, the southern Kuril Islands, and Sakhalin.

Duplication and Variation in the Major Histocompatibility Complex Genes in Blakiston's Fish Owl, Bubo blakistoni.

The major histocompatibility complex (MHC) includes many genes that are essential for the adaptive immune system, and variation in the antigen binding site (ABS) is related to resistance against pathogens. In the present study, quantitative real-time PCR indicated a larger number of MHC gene copies in the endangered population of Blakiston's fish owl (Bubo blakistoni) than in five other owl species, and massively parallel pyrosequencing detected more MHC class IIß per individual alleles in B. blakistoni than in the other species. A chromosomal fluorescence in situ hybridization (FISH) analysis showed that the MHC class I and class IIß loci are closely linked on a single pair of microchromosomes, indicating that the MHC genes were tandemly duplicated in a limited chromosomal region. Because B. blakistoni has twice as many MHC genes as its sister species, the tawny fish owl (Bubo flavipes), the duplication of MHC genes occurred after these species diverged by speciation. A Bayesian molecular phylogenetic analysis showed that the DAB1 and DAB2 lineages of MHC class IIß alleles from various strigid species each formed a separate clade, indicating that the two allelic lineages preceded the radiation of Strigidae and evolved as paralogs. By contrast, the ABS sequences did not form distinct clades between DAB1 and DAB2 alleles but were intermixed, presumably due to gene conversion. Despite the low diversity of alleles per locus, B. blakistoni had many lineages of MHC class IIß alleles. Gene duplication increases variation in the MHC genes in this species, and could have facilitated adaptation in small populations.

Low Genetic Variation of Red-Crowned Cranes on Hokkaido Island, Japan, Over the Hundred Years.

The red-crowned crane (Grus japonensis) is recognized internationally as an endangered species. Migratory populations breed in eastern Russia and northeastern China, whereas the resident population inhabits the island of Hokkaido, Japan. Although the population inhabiting Hokkaido had experienced a severe bottleneck by the end of the 19th century, the population size has recovered to about 1500 and continues to increase now thanks to conservation efforts. A previous study reported that no marked genetic differences were seen in the island population, and that the genetic variation of the whole population on Hokkaido was lower than that of the continental population. However, the precise genetic structure of the island population in the past or near present remains unclear. To better understand the spatiotemporal changes in the genetic structure of the island population, we performed mitochondrial DNA (mtDNA) analyses using stuffed specimens (years 1878-2001) and tissue or blood samples (years 1970-2014). We found three haplotypes in the island population, one of which was a novel mtDNA haplotype in 1997 and 2007 samples. In addition, there was no clear difference in the haplotype frequency through the time span. These results suggest that the low genetic variation of the island population persisted for the last hundred years. It is thus nearly impossible for the island population to recover its genetic variation in isolation. Conservation plans for this species should therefore include the promotion of genetic exchanges between the continental and island populations, such as through artificial introduction to Hokkaido.

Genetic variation of major histocompatibility complex genes in the endangered red-crowned crane.

Populations that have drastically decreased in the past often have low genetic variation, which may increase the risk of extinction. The genes of major histocompatibility complex (MHC) play an important role in the adaptive immune response of jawed vertebrates. Maintenance of adaptive genetic diversity such as that of MHC genes is important for wildlife conservation. Here, we determined genotypes of exon 3 of MHC class IA genes (MHCIA) and exon 2 of MHC class IIB genes (MHCIIB) to evaluate genetic variation of the endangered red-crowned crane population on Hokkaido Island, Japan, which experienced severe population decline in the past. We identified 16 and 6 alleles of MHCIA and MHCIIB, respectively, from 152 individuals. We found evidence of a positive selection at the antigen-binding sites in MHCIA exon 3 and MHCIIB exon 2. The phylogenetic analyses indicated evidence of trans-species polymorphism among the crane MHC genes. The genetic variability in both classes of MHC genes at the population level was low. No geographic structure was found based on the genetic diversity of microsatellite and MHC genes. Our study provides useful data for the optimal management of the red-crowned crane population in Hokkaido and can contribute to future studies on MHC genes of the continental populations of the red-crowned crane and other crane species.

A Comparison of Visual and Genetic Techniques for Identifying Japanese Marten Scats - Enabling Diet Examination in Relation to Seasonal Food Availability in a Sub-Alpine Area of Japan.

We compared the reliability of visual diagnostic criteria to DNA diagnostic techniques, including newly designed primers, to discriminate Japanese marten (Martes melampus) feces from those of other sympatric carnivore species. Visual criteria proved > 95% reliable for fresh, odoriferous scats in good condition. Based upon this verification, we then examined if and how Japanese marten diet differs among seasons at high elevation study site (1500-2026 m). We also considered how intra-specific competition with the Japanese red fox (Vulpes vulpes japonica) may shape marten feeding ecology. From 120 Japanese marten fecal samples, high elevation diet comprised (frequency of occurrence) 30.6-66.0% mammals, 41.0-72.2% insects and 10.6-46.2% fruits, subject to seasonal variation, with a Shannon-Weaver index value of 2.77. These findings contrast substantially to seasonal marten diet reported in adjacent lowland regions (700-900 m), particularly in terms of fruit consumption, showing the trophic adaptability of the Japanese marten. We also noted a substantial dietary overlap with the red fox (n = 26 scats) with a Shannon-Weaver index of 2.61, inferring little trophic niche mutual exclusion (trophic niche overlap: 0.95), although some specific seasonal prey selection differences were likely related to relative differences in body size between foxes and martens. This additional information on the feeding ecology of the Japanese marten enables a better assessment of the specific risks populations face in mountainous regions.

Gene duplication and concerted evolution of mitochondrial DNA in crane species.

The gene duplication in mitochondrial DNA (mtDNA) has been reported in diverse bird taxa so far. Although many phylogenetic and population genetic analyses of cranes were carried out based on mtDNA diversity, whether mtDNA contains duplicated regions is unknown. To address the presence or absence of gene duplication in cranes and investigate the molecular evolutionary features of crane mtDNA, we analyzed the gene organization and the molecular phylogeny of mtDNA from 13 crane species. We found that the mtDNA in 13 crane species shared a tandem duplicated region, which consists of duplicated sequence sets including cytochrome b (Cytb), NADH6, control region (CR) and three genes of tRNA. The gene order in the duplicated region was identical among all the 13 crane species, and the nucleotide sequences found within each individual showed high similarities. In addition, phylogenetic trees based on homologous sequences of CR and Cytb indicated the possibility of concerted evolution among the duplicated genes. The results suggested that the duplication event occurred in the common ancestor of crane species or some older ancestors.

Population Structure of the Raccoon Dog on the Grounds of the Imperial Palace, Tokyo, Revealed by Microsatellite Analysis of Fecal DNA.

The raccoon dog (Nyctereutes procyonoides, Canidae, Carnivora) is highly adaptable to urban environments. Populations of carnivorans inhabiting urban areas sometimes differ ecologically and genetically from those in rural areas. However, there is little information on urban raccoon dogs. This study focused on raccoon dog populations in Tokyo, one of the most highly urbanized cities in the world. We examined the genotypes of 10 microsatellites for 101 fecal samples from raccoon dogs inhabiting the grounds of the Imperial Palace, a green space in central Tokyo. We successfully genotyped 58 samples originating from 31 individuals. We also analyzed muscle tissue samples from raccoon dogs from the grounds of the Imperial Palace, the Akasaka Imperial Grounds (a green space close to the Imperial Palace), and the surrounding urban area, and then investigated the genetic structure and diversity of these populations, and the genetic differentiation among them. The population on the grounds of the Imperial Palace was genetically differentiated from that in the Akasaka Imperial Grounds, suggesting that the roads and buildings act as barriers to gene flow. In addition, the population on the grounds of the Imperial Palace showed greater genetic difference from that in the surrounding area than that in the Akasaka Imperial Grounds. We speculate that the moats around the Imperial Palace restrict individual ranges within the palace grounds and limit migration and gene flow to other areas.

Spatial and temporal variation at major histocompatibility complex class IIB genes in the endangered Blakiston's fish owl.

INTRODUCTION: Quantifying intraspecific genetic variation in functionally important genes, such as those of the major histocompatibility complex (MHC), is important in the establishment of conservation plans for endangered species. The MHC genes play a crucial role in the vertebrate immune system and generally show high levels of diversity, which is likely due to pathogen-driven balancing selection. The endangered Blakiston's fish owl (Bubo blakistoni) has suffered marked population declines on Hokkaido Island, Japan, during the past several decades due to human-induced habitat loss and fragmentation. We investigated the spatial and temporal patterns of genetic diversity in MHC class IIß genes in Blakiston's fish owl, using massively parallel pyrosequencing. RESULTS: We found that the Blakiston's fish owl genome contains at least eight MHC class IIß loci, indicating recent gene duplications. An analysis of sequence polymorphism provided evidence that balancing selection acted in the past. The level of MHC variation, however, was low in the current fish owl populations in Hokkaido: only 19 alleles were identified from 174 individuals. We detected considerable spatial differences in MHC diversity among the geographically isolated populations. We also detected a decline of MHC diversity in some local populations during the past decades. CONCLUSIONS: Our study demonstrated that the current spatial patterns of MHC variation in Blakiston's fish owl populations have been shaped by loss of variation due to the decline and fragmentation of populations, and that the short-term effects of genetic drift have counteracted the long-term effects of balancing selection.

Recent fragmentation of the endangered Blakiston's fish owl (Bubo blakistoni) population on Hokkaido Island, Northern Japan, Revealed by Mitochondrial DNA and Microsatellite Analyses.

INTRODUCTION: Blakiston's fish owl (Bubo blakistoni) was previously widespread on Hokkaido Island, Japan, but is now distributed only in limited forest areas. The population size on Hokkaido decreased during the 20th century due to reduction and fragmentation of the owl's habitat. To elucidate temporal and spatial changes in population structure and genetic diversity, we analyzed 439 individuals collected over the last 100 years. RESULTS: We detected a population bottleneck and fragmentation event indicated by mitochondrial DNA (mtDNA) haplotype and microsatellite analyses. The lowest value for effective population size, which was estimated by moment and temporal methods from microsatellite data, occurred in the 1980s. Five haplotypes were found in the mtDNA control region; most haplotypes were previously widespread across Hokkaido, but have become fixed in separate areas after the bottleneck period. Genetic differentiation among local populations, as indicated by both mtDNA and microsatellite data, likely arose through population fragmentation. CONCLUSIONS: The owl population may have been divided into limited areas due to loss of habitats via human activities, and have lost genetic variability within the local populations through inbreeding. Our mtDNA and microsatellite data show that genetic diversity decreased in local populations, indicating the importance of individuals moving between areas for conservation of this species on Hokkaido.

Isolation and characterization of major histocompatibility complex class II B genes in cranes.

In this study, we isolated and characterized the major histocompatibility complex (MHC) class II B genes in cranes. Genomic sequences spanning exons 1 to 4 were amplified and determined in 13 crane species and three other species closely related to cranes. In all, 55 unique sequences were identified, and at least two polymorphic MHC class II B loci were found in most species. An analysis of sequence polymorphisms showed the signature of positive selection and recombination. A phylogenetic reconstruction based on exon 2 sequences indicated that trans-species polymorphism has persisted for at least 10 million years, whereas phylogenetic analyses of the sequences flanking exon 2 revealed a pattern of concerted evolution. These results suggest that both balancing selection and recombination play important roles in the crane MHC evolution.

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.