Global invasion history of the agricultural pest butterfly Pieris rapae revealed with genomics and citizen science.

The small cabbage white butterfly, Pieris rapae, is a major agricultural pest of cruciferous crops and has been introduced to every continent except South America and Antarctica as a result of human activities. In an effort to reconstruct the near-global invasion history of P. rapae, we developed a citizen science project, the "Pieris Project," and successfully amassed thousands of specimens from 32 countries worldwide. We then generated and analyzed nuclear (double-digest restriction site-associated DNA fragment procedure [ddRAD]) and mitochondrial DNA sequence data for these samples to reconstruct and compare different global invasion history scenarios. Our results bolster historical accounts of the global spread and timing of P. rapae introductions. We provide molecular evidence supporting the hypothesis that the ongoing divergence of the European and Asian subspecies of P. rapae (∼1,200 y B.P.) coincides with the diversification of brassicaceous crops and the development of human trade routes such as the Silk Route (Silk Road). The further spread of P. rapae over the last ∼160 y was facilitated by human movement and trade, resulting in an almost linear series of at least 4 founding events, with each introduced population going through a severe bottleneck and serving as the source for the next introduction. Management efforts of this agricultural pest may need to consider the current existence of multiple genetically distinct populations. Finally, the international success of the Pieris Project demonstrates the power of the public to aid scientists in collections-based research addressing important questions in invasion biology, and in ecology and evolutionary biology more broadly.

Climate-mediated hybrid zone movement revealed with genomics, museum collection, and simulation modeling.

Climate-mediated changes in hybridization will dramatically alter the genetic diversity, adaptive capacity, and evolutionary trajectory of interbreeding species. Our ability to predict the consequences of such changes will be key to future conservation and management decisions. Here we tested through simulations how recent warming (over the course of a 32-y period) is affecting the geographic extent of a climate-mediated developmental threshold implicated in maintaining a butterfly hybrid zone (Papilio glaucus and Papilio canadensis; Lepidoptera: Papilionidae). These simulations predict a 68-km shift of this hybrid zone. To empirically test this prediction, we assessed genetic and phenotypic changes using contemporary and museum collections and document a 40-km northward shift of this hybrid zone. Interactions between the two species appear relatively unchanged during hybrid zone movement. We found no change in the frequency of hybridization, and regions of the genome that experience little to no introgression moved largely in concert with the shifting hybrid zone. Model predictions based on climate scenarios predict this hybrid zone will continue to move northward, but with substantial spatial heterogeneity in the velocity (55-144 km/1 °C), shape, and contiguity of movement. Our findings suggest that the presence of nonclimatic barriers (e.g., genetic incompatibilities) and/or nonlinear responses to climatic gradients may preserve species boundaries as the species shift. Further, we show that variation in the geography of hybrid zone movement could result in evolutionary responses that differ for geographically distinct populations spanning hybrid zones, and thus have implications for the conservation and management of genetic diversity.

The role of latitudinal, genetic and temperature variation in the induction of diapause of Papilio glaucus (Lepidoptera: Papilionidae).

A key adaptation in insects for dealing with variable environmental conditions is the ability to diapause. The tiger swallowtail butterflies, Papilio glaucus and P. canadensis are ideal species to explore the genetic causes and population genetic consequences of diapause because divergence in this trait is believed to be a salient factor in maintaining a hybrid zone between these species. Yet little is known about the factors that influence diapause induction in this system. Here we explored how spatial (latitudinal), environmental (temperature) and genetic (hybridization) factors affect diapause induction in this system. Specifically, a series of growth chamber experiments using wild caught individuals from across the eastern United States were performed to: (1) evaluate how critical photoperiod varies with latitude, (2) isolate the stage in which induction occurs, (3) test whether changes in temperature affected rates of diapause induction, and (4) explore how the incidence of diapause is affected in hybrid offspring. We find that induction occurs in the larval stage, is not sensitive to a relatively broad range of temperatures, appears to have a complex genetic basis (i.e., is not simply a dominant trait following a Mendelian inheritance pattern) and that the critical photoperiod increases by 0.4 h with each increasing degree in latitude. This work deepens our understanding of how spatial, environmental and genetic variation influences a key seasonal adaptation (diapause induction) in a well-developed ecological model system and will make possible future studies that explore how climatic variation affects the population dynamics and genetics of this system.

Patterns of divergence across the geographic and genomic landscape of a butterfly hybrid zone associated with a climatic gradient.

Hybrid zones are a valuable tool for studying the process of speciation and for identifying the genomic regions undergoing divergence and the ecological (extrinsic) and nonecological (intrinsic) factors involved. Here, we explored the genomic and geographic landscape of divergence in a hybrid zone between Papilio glaucus and Papilio canadensis. Using a genome scan of 28,417 ddRAD SNPs, we identified genomic regions under possible selection and examined their distribution in the context of previously identified candidate genes for ecological adaptations. We showed that differentiation was genomewide, including multiple candidate genes for ecological adaptations, particularly those involved in seasonal adaptation and host plant detoxification. The Z chromosome and four autosomes showed a disproportionate amount of differentiation, suggesting genes on these chromosomes play a potential role in reproductive isolation. Cline analyses of significantly differentiated genomic SNPs, and of species-diagnostic genetic markers, showed a high degree of geographic coincidence (81%) and concordance (80%) and were associated with the geographic distribution of a climate-mediated developmental threshold (length of the growing season). A relatively large proportion (1.3%) of the outliers for divergent selection were not associated with candidate genes for ecological adaptations and may reflect the presence of previously unrecognized intrinsic barriers between these species. These results suggest that exogenous (climate-mediated) and endogenous (unknown) clines may have become coupled and act together to reinforce reproductive isolation. This approach of assessing divergence across both the genomic and geographic landscape can provide insight about the interplay between the genetic architecture of reproductive isolation and endogenous and exogenous selection.