Adaptation to plant defence in an agricultural insect pest: Integrating genome scans and gene expression in the soybean aphid reveals multi‐genic pathways

In agroecosystems, intense selection pressures cause species to adapt and spread, often leading to the evolution and persistence of pests. Understanding how pests rapidly adapt can help develop sustainable strategies for their management and improve agroecosystem health. Pest adaptation involves stable variations in DNA sequence, as well as dynamic shifts in gene expression, often mediated by non-coding regulatory elements. We examined adaptation to plant defences in the soybean aphid, Aphis glycines, in which virulent aphids have overcome plant defences and avirulent aphids have not. Previous data with laboratory colonies suggested that virulent aphids have higher overall gene expression, including transposable elements, some of which influence gene regulation. However, we lack information on how genetic variation in natural populations impacts adaptation and potentially gene regulation. We integrated population genome scans of field-collected, soybean aphid populations with gene expression profiles of virulent and avirulent laboratory colonies to uncover connections between genetic differentiation and gene regulation for virulence. Genome scan methods found 2144 single nucleotide polymorphisms (SNPs) with significant genetic differentiation (i.e., outliers) in field-collected populations. These SNPs were near 1004 genes, representing 5.16% of the effective number of genes. Based on previous RNA-Seq data with laboratory colonies, we found 3160 genes and 147 long non-coding RNAs (lncRNAs) with differential expression among virulent and avirulent biotypes. By integrating both data sets, we identified 16 genes and 5 long non-coding RNAs with differential expression and that were associated with an outlier SNP (within 10 kbp). We validated SNPs with additional field collected aphids and found an aphid clone with stronger virulence than our laboratory virulent colony, surviving on 2 different aphid-resistant soybean varieties. This new virulent clone had fixed allele differences at 9 SNPs compared to our avirulent and other virulent colony. Field collected soybean aphids matching the phenotype of this new virulent clone had significant genetic differentiation with 3 outlier SNPs near genes related to zinc transport and lachesin compared to field collected avirulent aphids. Our entire data reinforced the importance of a potential multi-genetic response to overcome plant defence and generates new insights into complex genetic and regulatory mechanisms involved in insect-plant interactions.