Inbreeding depression and drift load in small populations at demographic disequilibrium

Inbreeding depression is a major driver of mating system evolution and has critical implications for population viability. Theoretical and empirical attention has been paid to predicting how inbreeding depression varies with population size. Lower inbreeding depression is predicted in small populations at equilibrium, primarily due to higher inbreeding rates facilitating purging and/or fixation of deleterious alleles (drift load), but predictions at demographic and genetic disequilibrium are less clear. In this study, we experimentally evaluate how lifetime inbreeding depression and drift load, estimated by heterosis, vary with census (N-c) and effective (estimated as genetic diversity, H-e) population size across six populations of the biennial Sabatia angularis as well as present novel models of inbreeding depression and heterosis under varying demographic scenarios at disequilibrium (fragmentation, bottlenecks, disturbances). Our experimental study reveals high average inbreeding depression and heterosis across populations. Across our small sample, heterosis declined with H-e, as predicted, whereas inbreeding depression did not vary with H-e and actually decreased with N-c. Our theoretical results demonstrate that inbreeding depression and heterosis levels can vary widely across populations at disequilibrium despite similar H-e and highlight that joint demographic and genetic dynamics are key to predicting patterns of genetic load in nonequilibrium systems.