Mapping of quantitative trait loci (QTL) was used to investigate the genetic architecture of divergence in floral characters associated with the mating system, an important adaptive trait in angiosperms. Two species of Leptosiphon (Polemoniaceae), one strongly self-fertilizing (L. bicolor) and the other partially outcrossing (L. jepsonii), were crossed to produce F2 and both backcross progenies. For each crossing population, a linkage map was created using amplified fragment length polymorphism markers, and QTL were identified for several dimensions of floral size. For each of the five traits examined, three to seven QTL were detected, with independent datasets yielding congruent results in some but not all cases. The phenotypic effect of individual QTL was generally moderate. We estimated that many of the QTL were additive or showed dominance toward L. bicolor, whereas comparison of mean trait values for parental and cross progenies showed apparent overall dominance of L. jepsonii traits. Colocalization of QTL for different dimensions of floral size was consistent with high phenotypic correlations between floral traits. Substantial segregation distortion was observed in marker loci, the majority favoring alleles from the large-flowered parent. A low frequency of male sterility in the F2 population is consistent with the Dobzhansky-Muller model for the evolution of reproductive isolation.