Genomic Footprints of Recovery in the European Bison

Abstract

After extinction in the wild in the beginning of the 20th century, the European bison has been successfully recovered in 2 distinct genetic lines from only 12 and 7 captive founders. We here aimed at characterizing the levels of realized inbreeding in these 2 restored lines to provide empirical insights into the genomic footprints left by population recovery from a small number of founders. To that end, we genotyped 183 European bison born over the last 40 years with the Illumina BovineHD beadchip that contained 22 602 informative autosomal single-nucleotide polymorphisms after data filtering. We then identified homozygous-by-descent (HBD) segments and classified them into different age-related classes relying on a model-based approach. As expected, we observed that the strong and recent founder effect experienced by the 2 lines resulted in very high levels of recent inbreeding and in the presence of long HBD tracks (up to 120 Mb). These long HBD tracks were associated with ancestors living approximately from 4 to 32 generations in the past, suggesting that inbreeding accumulated over multiple generations after the bottleneck. The contribution to inbreeding of the most recent groups of ancestors was however found to be decreasing in both lines. In addition, comparison of Lowland individuals born at different time periods showed that the levels of inbreeding tended to stabilize, HBD segments being shorter in animals born more recently which indicates efficient control of inbreeding. Monitoring HBD segment lengths over generations may thus be viewed as a valuable genomic diagnostic tool for populations in conservation or recovery programs.

Keywords: European bison; ROH; RZooRoH; homozygous-by-descent; inbreeding; recovery.