A sustainable approach to address climate change and increasing water demand in agriculture is breeding for plant functional traits that conserve water and enhance climate resilience. Stomata regulate plant-water relations and are promising targets for crop improvement. Here, we investigate the variation in stomatal density (SD) in a diverse apple population (Malus domestica Borkh.) consisting of 269 accessions. Genome-wide association studies identified robust associations with SD on chromosomes 2, 9, and 10 (classified as SNPs with p value higher than adjusted Bonferroni threshold of -log10(p) > 8.78 that were consistently identified across datasets). On chromosome 9, a candidate gene that negatively regulates stomatal development, EPIDERMAL PATTERNING FACTOR 1 (EPF1), was identified inside a genomic region of 241 kb determined by six robust associations. On chromosome 10, a positive regulator candidate gene, EPIDERMAL PATTERNING FACTOR LIKE 9 (STOMAGEN), was identified 1680 kb from the robust association. Identification of positive (STOMAGEN) and negative (EPF1) regulators of SD suggest potential antagonistic roles at the population scale in determining SD. On chromosome 2, a gene co-expression analysis identified a gene cluster containing both EPF1 and STOMAGEN together with a novel candidate gene, CYTOCHROME P450 (CYP77A4), that was located 544 kb from the robust association. The percentage of SD phenotypic variance explained by each robust association was between 7% and 10%. These findings provide a foundation for understanding SD variation at the population scale and opportunities to modulate SD by genomics-assisted breeding strategies.