Egg production is an important economic trait in layer ducks and understanding the genetics basis is important for their breeding. In this study, a genome-wide association study ( GWAS) for egg production traits in 303 female Longyan Shan-ma ducks was performed based on a genotyping-by-sequencing strategy. Sixty-two single nucleotide polymorphisms ( SNPs) associated with egg weight traits were identified ( P < 9.48 × 10 -5), including 8 SNPs at 5% linkage disequilibrium ( LD)-based Bonferroni-corrected genome-wide significance level ( P < 4.74 × 10 -6). One hundred and nineteen SNPs were associated with egg number traits ( P < 9.48 × 10 -5), including 13 SNPs with 5% LD-based Bonferroni-corrected genome-wide significance ( P < 4.74 × 10 -6). These SNPs annotated 146 target genes which contained known candidate genes for egg production traits, such as prolactin and prolactin releasing hormone receptor. This study identified that these associated genes were significantly enriched in egg production-related pathways ( P < 0.05), such as the oxytocin signaling, MAPK signaling, and calcium signaling pathways. It was notable that 18 genes were differentially expressed in ovarian tissues between higher and lower egg production in Shan-ma ducks. The identified potential candidate genes and pathways provide insight into the genetic basis underlying the egg production trait of layer ducks.