Landscape and climate determine patterns of spread for all colour morphs of the alien ladybird Harmonia axyridis

Aim Understanding why some areas are more vulnerable to invasion by invasive alien species than others represents a key challenge in invasion biology. We investigated the roles of landscape, climate, proximity to initial invasion sites, recording intensity and colour pattern polymorphism in explaining the spread of the invasive ladybird Harmonia axyridis across Great Britain. Location The British mainland and inshore islands (Anglesey, Isle of Wight and the Inner Hebrides). Methods Patterns in time to invasion of H. axyridis were quantified at a hectad level (10 km grid squares) from over 23,000 records. Bayesian spatial survival models were developed relating spread patterns to concurrent patterns in climate, land cover and recording intensity. Separate models were developed for melanic (f. conspicua and f. spectabilis) and non-melanic (f. succinea) colour forms that were hypothesized to differ in their response to climate factors, particularly temperature, and subsequent spread. Results Harmonia axyridis spread more rapidly into hectads containing a high proportion of urban land cover and marginally more slowly into hectads containing high conifer cover. Invasions were recorded more rapidly in hectads with high recording intensity, or that were close to initial invasion sites. Harmonia axyridis colour morphs spread to similar extents, although increased sunshine significantly enhanced the spread of the non-melanic form (f. succinea). Correspondence between observed and predicted times of invasion from survival models was substantial (Pearson's rp for all H. axyridis colour forms = 0.834). Main conclusions The ability of H. axyridis to use a wide range of habitats and prey species and overwinter inside buildings may have allowed this species to establish rapidly in urban areas. Native ladybird species in urban areas may suffer combined negative impacts of habitat deterioration and competition from H. axyridis. Spatial survival models help us to identify areas vulnerable to invasion and to understand the impacts of multiple environmental drivers on the spread of invasive alien species, given detailed data on time to invasion.

[1]  L. Reino,et al.  Using survival regression to study patterns of expansion of invasive species: will the common waxbill expand with global warming? , 2009 .

[2]  P. Chesson,et al.  Community ecology theory as a framework for biological invasions , 2002 .

[3]  Petr Pyšek,et al.  European map of alien plant invasions based on the quantitative assessment across habitats , 2009 .

[4]  M. McClure.  Role of Predators in Regulation of Endemic Populations of Matsucoccus matsumurae (Homoptera: Margarodidae) in Japan , 1986 .

[5]  Sw. Banerjee,et al.  Hierarchical Modeling and Analysis for Spatial Data , 2003 .

[6]  P. Grambsch,et al.  A Package for Survival Analysis in S , 1994 .

[7]  H. Rue,et al.  Approximate Bayesian inference for hierarchical Gaussian Markov random field models , 2007 .

[8]  Marius Gilbert,et al.  Long-distance dispersal and human population density allow the prediction of invasive patterns in the horse-chestnut leafminer Cameraria ohridella , 2004 .

[9]  P. Brakefield.  Ecological studies on the polymorphic ladybird Adalia bipunctata in The Netherlands. I. Population biology and geographical variation of melanism , 1984 .

[10]  Catherine S. Jarnevich,et al.  Challenges of predicting the potential distribution of a slow-spreading invader: a habitat suitability map for an invasive riparian tree , 2010, Biological Invasions.

[11]  P. Brakefield,et al.  The basis of thermal melanism in the ladybird Adalia bipunctata: Differences in reflectance and thermal properties between the morphs , 1985, Heredity.

[12]  A. Estoup,et al.  The harlequin ladybird, Harmonia axyridis , 2016 .

[13]  G. Muenchow.  Ecological Use of Failure Time Analysis , 1986 .

[14]  Neal Alexander,et al.  Bayesian Disease Mapping: Hierarchical Modeling in Spatial Epidemiology , 2011 .

[15]  Bradley P. Carlin,et al.  Bayesian measures of model complexity and fit , 2002 .

[16]  Helen E. Roy,et al.  Lessons from lady beetles: accuracy of monitoring data from US and UK citizen-science programs , 2012 .

[17]  Andrew Gelman,et al.  Data Analysis Using Regression and Multilevel/Hierarchical Models , 2006 .

[18]  R. Koch.  The multicolored Asian lady beetle, Harmonia axyridis: A review of its biology, uses in biological control, and non-target impacts , 2003, Journal of insect science.

[19]  D. Roy,et al.  Harmonia axyridis in Great Britain: analysis of the spread and distribution of a non-native coccinellid , 2008, BioControl.

[20]  T. Sparks,et al.  Decline in native ladybirds in response to the arrival of Harmonia axyridis: early evidence from England , 2011 .

[21]  D. Lonsdale,et al.  Melanism in Adalia bipunctata L. (Col., Coccinellidae) and its Relationship to Atmospheric Pollution , 1975 .

[22]  F. Jiggins,et al.  Melanic through nature or nurture: genetic polymorphism and phenotypic plasticity in Harmonia axyridis , 2010, Journal of evolutionary biology.

[23]  R. Meentemeyer,et al.  Equilibrium or not? Modelling potential distribution of invasive species in different stages of invasion , 2012 .

[24]  David R. Anderson,et al.  Model selection and multimodel inference : a practical information-theoretic approach , 2003 .

[25]  S. Dullinger,et al.  Modelling climate change‐driven treeline shifts: relative effects of temperature increase, dispersal and invasibility , 2004 .

[26]  T. Steenberg,et al.  The harlequin ladybird (Harmonia axyridis Pallas) in Denmark: spread and phenology during the initial phase of invasion. , 2009 .

[27]  D. Roy,et al.  Invasive alien predator causes rapid declines of native European ladybirds , 2012 .

[28]  M. Majerus,et al.  Aphid‐mediated coexistence of ladybirds (Coleoptera: Coccinellidae) and the wood ant Formica rufa: seasonal effects, interspecific variability and the evolution of a coccinellid myrmecophile. , 2000 .

[29]  M. Majerus,et al.  The ant-associations and diet of the ladybird Coccinella magnifica (Coleoptera: Coccinellidae) , 2002 .

[30]  H. Roy,et al.  The potential impacts of the arrival of the harlequin ladybird, Harmonia axyridis (Pallas) (Coleoptera: Coccinellidae), in Britain , 2006 .

[31]  Xavier Font,et al.  Habitat invasions by alien plants: a quantitative comparison among Mediterranean, subcontinental and oceanic regions of Europe , 2008 .

[32]  Remy L. Ware,et al.  Intraguild predation of immature stages of British and Japanese coccinellids by the invasive ladybird Harmonia axyridis , 2008, BioControl (Dordrecht).

[33]  P. Brakefield,et al.  A steep cline in ladybird melanism has decayed over 25 years: a genetic response to climate change? , 2011, Heredity.

[34]  Richard A. Wadsworth,et al.  Final Report for LCM2007 - the new UK land cover map. Countryside Survey Technical Report No 11/07 , 2011 .

[35]  J. Klein,et al.  Survival Analysis: Techniques for Censored and Truncated Data , 1997 .

[36]  Helen E. Roy,et al.  Using biological traits to explain ladybird distribution patterns , 2012 .

[37]  Helen E. Roy,et al.  Understanding citizen science and environmental monitoring: final report on behalf of UK Environmental Observation Framework , 2012 .

[38]  A. Estoup,et al.  The global spread of Harmonia axyridis (Coleoptera: Coccinellidae): distribution, dispersal and routes of invasion , 2011, BioControl.

[39]  A. Peterson,et al.  Evidence of climatic niche shift during biological invasion. , 2007, Ecology letters.

[40]  G. F. Peterken,et al.  Natural Woodland: Ecology and Conservation in Northern Temperate Regions. , 1996 .

[41]  A. Masson,et al.  Seasonal phenotypic plasticity: wild ladybirds are darker at cold temperatures , 2011, Evolutionary Ecology.

[42]  W. Gerrard.  Effect of Temperature , 1976 .

[43]  P. Borges,et al.  Harmonia axyridis: What will stop the invader? , 2008, BioControl.

[44]  W. Hutchison,et al.  The chemical ecology of Harmonia axyridis , 2011, BioControl.

[45]  H. Rue,et al.  Approximate Bayesian Inference for Survival Models , 2010 .

[46]  H. Rue,et al.  Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations , 2009 .

[47]  D. Richardson,et al.  Plant Invasions and Invasibility of Plant Communities , 2013 .

[48]  Challenges in predicting invasive reservoir hosts of emerging pathogens: mapping Rhododendron ponticum as a foliar host for Phytophthora ramorum and Phytophthora kernoviae in the UK , 2013, Biological Invasions.

[49]  Steven J. Phillips,et al.  The art of modelling range‐shifting species , 2010 .

[50]  Glenn Marion,et al.  Bayesian Inference for the Spatio-Temporal Invasion of Alien Species , 2007, Bulletin of mathematical biology.

[51]  J. Besag,et al.  Bayesian image restoration, with two applications in spatial statistics , 1991 .

[52]  N. Belyakova.  Polymorphism of the harlequin ladybird Harmonia axyridis (Coleoptera, Coccinellidae) Baikal population , 2013, Entomological Review.

[53]  C. Pipper,et al.  Failure Time Analysis , 2009 .

[54]  A. Estoup,et al.  Investigating the genetic load of an emblematic invasive species: the case of the invasive harlequin ladybird Harmonia axyridis , 2013, Ecology and evolution.

[55]  Andrew B Lawson,et al.  Bayesian parametric accelerated failure time spatial model and its application to prostate cancer , 2011, Journal of applied statistics.

[56]  M. Majerus,et al.  Adaptations of Coccinella magnifica, a myrmecophilous coccinellid to aggression by wood ants (Formica rufa group). II. Larval behaviour, and ladybird oviposition location1 , 2003 .

[57]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[58]  David N. Bonter,et al.  Citizen Science as an Ecological Research Tool: Challenges and Benefits , 2010 .

[59]  R. Cowling,et al.  Rainfall reliability, a neglected factor in explaining convergence and divergence of plant traits in fire-prone mediterranean-climate ecosystems , 2005 .