Identifying signatures of selection at the enhancer of split neurogenic gene complex in Drosophila

The Enhancer of split gene complex (E(spl)-C) is one of the more highly annotated gene regions in Drosophila, and the 12 genes within the complex help determine the spacing and patterning of adult bristles. Any E(spl)-C coding, transcribed, or cis-regulatory regions experiencing nonneutral evolution are strong candidates to harbor polymorphisms contributing to naturally occurring variation in bristle number. We confirm that the E(spl)-C is strongly conserved and show that 74% of regulatory elements previously identified in D. melanogaster are conserved in D. pseudoobscura. Regulatory elements in enhancer regions show lower nucleotide diversity and more rare polymorphisms compared with adjacent nonregulatory DNA, suggesting they are under purifying selection, and these effects are particularly pronounced when considering only conserved regulatory elements. The ratio of polymorphism to divergence was significantly different between binding sites and nonbinding sites for transcription factors within enhancer regions, suggesting the action of some form of selection. Too few polymorphisms in regions of the 3' UTR harboring regulatory motifs prevents adequate comparison of diversity and the polymorphism frequency spectrum between 3' UTR motif and nonmotif sequence. We identified at least two broad regions of the gene complex showing strong population subdivision among four populations, which is suggestive of local adaptation or background selection. Finally, two regions of the E(spl)-C exhibit low nucleotide diversity, a high level of rare polymorphisms, and an increase in linkage disequilibrium, which together suggest the action of positive selection. Notably, the gene m2 shows a significant deviation from neutrality by the McDonald-Kreitman test and resides in one of the two regions putatively experiencing a selective sweep. All sites in regions apparently visible to various selective forces are candidates for future work to determine their phenotypic effects.