Genomic islands of heterozygosity maintained across caribou populations despite inbreeding

Small, isolated populations are prone to inbreeding, increasing the proportion of homozygous sites across the genome that can be quantified as runs of homozygosity (ROH). Caribou (Rangifer tarandus) are declining across their range in Canada; thus, understanding the effects of inbreeding on genetic potential is pertinent for conserving small, isolated populations. We quantified ROH in high-coverage whole genomes of boreal caribou from small, isolated populations in southern Ontario, Canada, in comparison to caribou from the continuous range of Ontario, other caribou ecotypes in Canada, and western Greenland. Sampled populations presented divergent evolutionary histories, differing population sizes, and extents of isolation. We conducted BLAST searches across regions of elevated heterozygosity to identify genes that have maintained variation despite inbreeding. We found caribou from recently isolated populations in Ontario had a large proportion of their genome in long ROH. We observed even larger proportions but shorter ROH in western Greenland, indicating that inbreeding has occurred over a longer period in comparison to other populations. We observed the least inbreeding in barren-ground and eastern migratory caribou, which occur in larger population sizes than boreal caribou. Despite vastly different inbreeding extents, we found regions of high heterozygosity maintained across all populations. Within these islands of heterozygosity, we identified genes associated with immunity, signaling regulation, nucleotide binding, toxin elimination, and feeding behaviour regulation. In this study, we confirm inbreeding in isolated populations of a species at risk, but also uncover high variation in some genes maintained across divergent populations despite inbreeding, suggesting strong balancing selection.