Complementarity of statistical treatments to reconstruct worldwide routes of invasion: the case of the Asian ladybird Harmonia axyridis

Inferences about introduction histories of invasive species remain challenging because of the stochastic demographic processes involved. Approximate Bayesian computation (ABC) can help to overcome these problems, but such method requires a prior understanding of population structure over the study area, necessitating the use of alternative methods and an intense sampling design. In this study, we made inferences about the worldwide invasion history of the ladybird Harmonia axyridis by various population genetics statistical methods, using a large set of sampling sites distributed over most of the species’ native and invaded areas. We evaluated the complementarity of the statistical methods and the consequences of using different sets of site samples for ABC inferences. We found that the H. axyridis invasion has involved two bridgehead invasive populations in North America, which have served as the source populations for at least six independent introductions into other continents. We al...