Disentangling the worldwide invasion process of Halyomorpha halys through approximate Bayesian computation

Halyomorpha halys is a polyphagous insect pest, which was first found outside its native regions in 1996, and since then it has rapidly spread worldwide causing damage to agriculture. Investigating the genetic diversity among H. halys populations is essential to understand the colonization history out of its native areas. Analyses based on mtDNA indicate multiple invasions from Asia to European and American continents, as well as serial invasions within invaded countries, but the colonization dynamics are still debated. Recently, genome-wide data (ddRAD) have been published to investigate the genomic structure of H. halys, proposing China as a putative source of multiple serial invasion events to Europe and the Americas. In this study we re-analyze published ddRAD sequences from worldwide populations of H. halys to better elucidate the colonization process. We assessed the genetic diversity in native populations identifying genetic differentiation between populations in China. Furthermore, we observed a complex pattern of population structure in the invaded countries, that may have originated from the occurrence of multiple independent colonization waves through time, from sub-populations present in the native range to Europe and the Americas. We tested alternative colonization hypotheses through Approximate Bayesian Computation comparison of demographic scenarios. Our results support multiple waves of migration from East China to invaded territories and the occurrence of European and American bridgehead effects. These results underline the importance of demographic inference through genome-wide data to investigate biological invasions, whose knowledge become fundamental to establish new strategies of management and control of invasive species.