Whole-genome sequencing revealed local genetic differences and multi-layered genetic structure in Asian brown bear...

A multi-layered genetic structure shaped by multiple demographic events across different time periods has been well documented in some species. However, evaluating this pattern remains challenging in many wide-ranging wildlife because genome-wide datasets with broad geographic coverage are still scarce. Brown bears are genetically differentiated among regional populations, but the genomic data are still limited and the demographic events remain incompletely understood. Here, we generated whole-genome resequencing data from brown bears in Asian local populations previously reported to show significant variation and combined these data with ancient samples for population genomic analyses. Our results support previously reported demographic events, such as the early divergence of Central Asian populations. Furthermore, we found genetic differences among Asian brown bear populations, suggesting that population-specific demographic histories occurring at different times have contributed to their unique genetic characteristics. In the Central and West Asian populations, divergence times estimated from autosomal genomes differed by hundreds of thousands of years from those inferred from mitochondrial DNA and Y-chromosomal DNA haplotypes, suggesting that gene flow continued after regional dispersal. In addition, the Y-chromosomal phylogeny and ƒ4 statistics indicated genetic differentiation between the Kazakhstan and Tibetan populations, likely resulting from male-biased introgression into the Tibetan population. Finaly, the Hokkaido and Etorofu (Iturup) populations retained alleles derived from hybridization with polar bears, likely due to the limited impact of post-LGM population expansion. These findings indicate that the genetic variation among brown bear populations is shaped by multiple demographic events occurring at different times, resulting in a multi-layered genetic structure.