25Recently, increasing attempts have been made to understand how plant genes function in 26 natura studies. To determine whether plant defense genes are activated under multiple biotic 27 stimuli, we combined a high-throughput RNA-Seq with insect survey data on 19 accessions 28 of Arabidopsis thaliana growing on the field site of Switzerland. We found that genes with 29 GO annotations "glucosinolate biosynthetic process" and "response to insects" were the most 30 significantly enriched, exhibiting largely variable expression among plant accessions. Nearly 31 half of the total expression variation in glucosinolate biosynthetic genes, AOPs, ESM1, ESP, 32 and TGG1, was explained by among-accession variance. Combined with the field RNA-Seq 33 data, bioassays confirmed that AOP3 was up-regulated in response to the mustard aphid 34 Lipaphis erysimi. In addition, we also found that the expression of a major cis-jasmone 35 activated gene CYP81D11 was positively correlated with the number of the flea beetles 36 Phyllotreta spp. The combined results from RNA-Seq and insect surveys suggested that 37 plants can activate their defenses even when they are exposed to multiple biotic stimuli in 38 natura. 39