Defense by duplication: The relation between phenotypic glyphosate resistance and EPSPS gene copy number variation in Amaranthus palmeri

Gene copy number variation (CNV) has been increasingly associated with organismal responses to environmental stress, but we know little about the quantitative relation between CNV and phenotypic variation. In this study we quantify variation in EPSPS (5-enolpyruvylshikimate-3-phosphate synthase) copy number using digital drop PCR and variation in phenotypic glyphosate resistance in 22 populations of Amaranthus palmeri (Palmer Amaranth), a range-expanding agricultural weed. Overall, we detected a significant positive relation between population mean copy number and mean resistance. The majority of populations exhibited high glyphosate resistance, yet maintained low-resistance individuals resulting in bimodality in many populations. We investigated linear and threshold models for the relation between copy number and resistance, and found evidence for a threshold of ~15 EPSPS copies: there was a steep increase in resistance before the threshold, followed by a much shallower slope. Moreover, as copy number increases, the range of variation in resistance decreases. This suggests a working hypothesis that as EPSPS copies and dosage increases, negative epistatic interactions may be compensated. We detected a quadratic relation between mean resistance and variation (s.d.) in resistance, consistent with the prediction that as phenotypic resistance increases in populations, stabilizing selection decreases variation in the trait. Finally, patterns of variation across the landscape are consistent with less variation among populations in mean copy number / resistance in Georgia where glyphosate resistance was first detected, and wider variation among populations in resistance and copy number in a more northern state where resistance evolution may be at a younger evolutionary state.