Genome evolution in an agricultural pest following adoption of transgenic crops

Replacing synthetic insecticides with transgenic crops for pest management has been economically and environmentally beneficial, but these benefits erode as pests evolve resistance. It has been proposed that novel genomic approaches could track molecular signals of emerging resistance to aid in resistance management. To test this, we quantified patterns of genomic change in Helicoverpa zea, a major lepidopteran pest and target of transgenic Bacillus thuringiensis (Bt) crops, between 2002 and 2017 as both Bt crop adoption and resistance increased in North America. Genomic scans of wild H. zea were paired with quantitative trait locus (QTL) analyses and showed the genomic architecture of field-evolved Cry1Ab resistance was polygenic, likely arising from standing genetic variation. Resistance to pyramided Cry1A.105 and Cry2Ab2 toxins was controlled by fewer loci. Of the 11 previously described Bt resistance genes, 9 showed no significant change over time or major effects on resistance. We were unable to rule out a contribution of aminopeptidases (apns), as a cluster of apn genes were found within a Cry-associated QTL. Molecular signals of emerging Bt resistance were detectable as early as 2012 in our samples, and we discuss the potential and pitfalls of whole-genome analysis for resistance monitoring based on our findings. This first study of Bt resistance evolution using whole-genome analysis of field-collected specimens demonstrates the need for a more holistic approach to examining rapid adaptation to novel selection pressures in agricultural ecosystems.