Understanding the ecological and evolutionary mechanisms that shape a species’ range is an important goal in evolutionary biology. Evidence indicates that mating system is an effective predictor of the global range of native species or naturalized alien plants, but the mechanisms underlying this predictability are not elaborated. Here, we develop a theoretical model to account for the ranges of plants under different mating systems based on migration-selection processes (an idea proposed by Haldane). The model includes alternation of gametophyte and sporophyte generations in one life cycle and the dispersal of haploid pollen and diploid seeds as vectors for gene flow. We show that the interaction between selfing rates and gametophytic selection determines the role of mating system in shaping a species’ range. Selfing restricts the species’ range under gametophytic selection in non-random mating systems, but expands the species’ range under the absence of gametophytic selection in any mating system. Gametophytic selection slightly restricts the species’ range in random mating. Both logarithmic and logistic models of population demography yield similar conclusions in the case of fixed or evolving genetic variance. The theory also helps to explain a broader relationship between mating system and range size following biological invasion or plant naturalization.