Subtelomeric 5-enolpyruvylshikimate-3-phosphate synthase (EPSPS) copy number variation confers glyphosate resistance in Eleusine indica

Genomic structural variation (SV) can have profound effects on an organism’s evolution, often serving as a novel source of genetic variation. Gene copy number variation (CNV), a specific form of SV, has repeatedly been associated with adaptive evolution in eukaryotes, especially to biotic and abiotic stresses. Resistance to the most widely used herbicide, glyphosate, has evolved through target-site CNV in many weedy plant species, including the economically important cosmopolitan grass, Eleusine indica (goosegrass); however, the origin and mechanisms of these resistance CNVs remain elusive in many weed species due to limited genetic and genomics resources. In order to study the target site CNV in goosegrass, we generated high-quality reference genomes for both glyphosate-susceptible and -resistant individuals, fine assembled the duplication of glyphosate’s target site gene enolpyruvylshikimate-3-phosphate synthase (EPSPS), and revealed a novel rearrangement of EPSPS into the subtelomeric region of the chromosomes, ultimately leading to herbicide resistance evolution. This discovery adds to the limited knowledge of the importance of subtelomeres as rearrangement hotspots and novel variation generators as well as provides an example of yet another unique pathway for the formation of CNVs in plants.