A Guide for Developing Demo-Genetic Models to Simulate Genetic Rescue

Genetic rescue is a conservation management strategy that reduces the negative effects of genetic drift and inbreeding in small and isolated populations. However, such populations might already be vulnerable to random fluctuations in growth rates (demographic stochasticity). Therefore, the success of genetic rescue depends not only on the genetic composition of the source and target populations but also on the emergent outcome of interacting demographic processes and other stochastic events. Developing predictive models that account for feedback between demographic and genetic processes (‘demo-genetic feedback’) is therefore necessary to guide the implementation of genetic rescue to minimize the risk of extinction of threatened populations. Here, we explain how the mutual reinforcement of genetic drift, inbreeding, and demographic stochasticity increases extinction risk in small populations. We then describe how these processes can be modelled by parameterizing underlying mechanisms, including deleterious mutations with partial dominance and demographic rates with variances that increase as abundance declines. We combine our suggestions of model parameterization with a comparison of the relevant capability and flexibility of five open-source programs designed for building genetically explicit, individual-based simulations. Using one of the programs, we provide a heuristic model to demonstrate that simulated genetic rescue can delay extinction of small virtual populations that would otherwise be exposed to greater extinction risk due to demo-genetic feedback. We then use a case study of threatened Australian marsupials to demonstrate that published genetic data can be used in one or all stages of model development and application, including parameterization, calibration, and validation. We highlight that genetic rescue can be simulated with either virtual or empirical sequence variation (or a hybrid approach) and suggest that model-based decision-making should be informed by ranking the sensitivity of predicted probability/time to extinction to variation in model parameters (e.g., translocation size, frequency, source populations) among different genetic-rescue scenarios.