The genus Narcissus L. (Amaryllidaceae) provides a model system to study the evolution and maintenance of sexual polymorphisms. In this study, we characterized microsatellite markers for N. dubius, N. cuatrecasasii, N. assoanus and N. rupicola for studies of genetic diversity and paternity analyses to investigate the stability of stylar dimorphism. We proved 40 new primer pairs from a genomic library of N. papyraceus and 12 microsatellite markers characterized also for N. papyraceus in a previous study (52 primer pairs overall). Twenty markers amplified, but their transferability and variability were different among species. Polymorphism was tested at least on 74 individuals and one population per species. The number of polymorphic loci per species ranged from four to eight. The number of alleles per locus ranged from two to 19 and the observed heterozygosity and gene diversity, from 0.107 to 0.729 and 0.103 to 0.894, respectively. These markers can be used for studies of genetic diversity and paternity analyses among individuals of N. dubius, N. cuatrecasasii, N. assoanus and N. rupicola to study the stability of stylar dimorphism.