Postglacial ecotype formation under outcrossing and self‐fertilization in Arabidopsis lyrata

The formation of ecotypes has been invoked as an important driver of postglacial biodiversity, because many species colonized heterogeneous habitats and experienced divergent selection. Ecotype formation has been predominantly studied in outcrossing taxa, while far less attention has been paid to the implications of mating system shifts. Here, we addressed whether substrate‐related ecotypes exist in selfing and outcrossing populations of Arabidopsis lyrata subsp. lyrata and whether the genomic footprint differs between mating systems. The North American subspecies colonized both rocky and sandy habitats during postglacial range expansion and shifted the mating system from predominantly outcrossing to predominantly selfing in a number of regions. We performed an association study on pooled whole‐genome sequence data of 20 selfing or outcrossing populations, which suggested genes involved in adaptation to substrate. Motivated by enriched gene ontology terms, we compared root growth between plants from the two substrates in a common environment and found that plants originating from sand grew roots faster and produced more side roots, independent of mating system. Furthermore, single nucleotide polymorphisms associated with substrate‐related ecotypes were more clustered among selfing populations. Our study provides evidence for substrate‐related ecotypes in A. lyrata and divergence in the genomic footprint between mating systems. The latter is the likely result of selfing populations having experienced divergent selection on larger genomic regions due to higher genome‐wide linkage disequilibrium.     

Mating system as a barrier to gene flow

Understanding mating system as one of reproductive isolating barriers remains important although this barrier is classified in a different sense from behavioral, ecological, and mechanical isolating barriers. Selfing enhances incipient speciation while outcrossing facilitates species integrity. Here, I study how mating system affects gene exchanges between genetically diverging species in a hybrid zone. Results show that a predominant selfing species has a greater barrier to selective gene flow than does a predominant outcrossing species. Barrier to neutral gene flow convexly changes with the selfing rate due to linkage disequilibrium, with a maximum at around intermediate selfing rate. Asymmetric transient or steady‐state barriers to neutral gene flow occur between two sides of a hybrid zone when the neutral gene is affected by its linked selective gene whose alternative alleles are adaptive to heterogeneous habitats. Selfing interacts with both a physical barrier and a density‐dependent ecological regulation (a logarithmic model) to strengthen the barriers to neutral and selective gene flow. This theory helps to interpret incipient speciation driven by selfing or to explain the asymmetric gene flow or unequal genomic mixtures between closely related species caused by their asymmetric mating systems in natural hybrid zones.     

Mating Systems of Psychotria tenuinervis (Rubiaceae): Distance from Anthropogenic and Natural Edges of Atlantic Forest Fragment

The aim of this study was to determine and compare the mating systems among Psychotria tenuinervis populations at anthropogenic edges, natural edges, and the forest interior using allozyme electrophoresis of naturally pollinated progeny arrays. P. tenuinervis showed low outcrossing rates, varying from 37% to 50% of the mating attributable to outcrossing and 50% to 63% attributable to self-fertilization, in the three habitats. The forest interior had the highest outcrossing rate (t _m = 0.50 and t _s = 0.43) among the three habitats. However, there were no differences in either multilocus or single-locus rates among the three habitats, indicating that the contribution of biparental inbreeding to the apparent selfing rate in these populations was very low. The multilocus (t _m) and single-locus (t _s) outcrossing rates for the P. tenuinervis in the sample plots within each habitat showed great heterogeneity. In conclusion, edge creation seems not to influence its mating systems. Additionally, although P. tenuinervis is a distylous species, the population’s inbreeding can be attributed almost entirely to self-fertilization.     

Consequences of multiple mating‐system shifts for population and range‐wide genetic structure in a coastal dune plant

Evolutionary transitions from outcrossing to selfing can strongly affect the genetic diversity and structure of species at multiple spatial scales. We investigated the genetic consequences of mating‐system shifts in the North American, Pacific coast dune endemic plant Camissoniopsis cheiranthifolia (Onagraceae) by assaying variation at 13 nuclear (n) and six chloroplast (cp) microsatellite (SSR) loci for 38 populations across the species range. As predicted from the expected reduction in effective population size (N_e) caused by selfing, small‐flowered, predominantly selfing (SF) populations had much lower nSSR diversity (but not cpSSR) than large‐flowered, predominantly outcrossing (LF) populations. The reduction in nSSR diversity was greater than expected from the effects of selfing on N_e alone, but could not be accounted for by indirect effects of selfing on population density. Although selfing should reduce gene flow, SF populations were not more genetically differentiated than LF populations. We detected five clusters of nSSR genotypes and three groups of cpSSR haplotypes across the species range consisting of parapatric groups of populations that usually (but not always) differed in mating system, suggesting that selfing may often initiate ecogeographic isolation. However, lineage‐wide genetic variation was not lower for selfing clusters, failing to support the hypothesis that selection for reproductive assurance spurred the evolution of selfing in this species. Within three populations where LF and SF plants coexist, we detected genetic differentiation among diverged floral phenotypes suggesting that reproductive isolation (probably postzygotic) may help maintain the striking mating‐system differentiation observed across the range of this species.     

THE EVOLUTION OF SELF-FERTILIZATION AND INBREEDING DEPRESSION IN PLANTS. I. GENETIC MODELS

The amounts of inbreeding depression upon selfing and of heterosis upon outcrossing determine the strength of selection on the selfing rate in a population when this evolves polygenically by small steps. Genetic models are constructed which allow inbreeding depression to change with the mean selfing rate in a population by incorporating both mutation to recessive and partially dominant lethal and sublethal alleles at many loci and mutation in quantitative characters under stabilizing selection. The models help to explain observations of high inbreeding depression (> 50%) upon selfing in primarily outcrossing populations, as well as considerable heterosis upon outcrossing in primarily selfing populations. Predominant selfing and predominant outcrossing are found to be alternative stable states of the mating system in most plant populations. Which of these stable states a species approaches depends on the history of its population structure and the magnitude of effect of genes influencing the selfing rate.     

The evolutionary response of mating system to heterosis

Isolation allows populations to diverge and to fix different alleles. Deleterious alleles that reach locally high frequencies contribute to genetic load, especially in inbred or selfing populations, in which selection is relaxed. In the event of secondary contact, the recessive portion of the genetic load is masked in the hybrid offspring, producing heterosis. This advantage, only attainable through outcrossing, should favour evolution of greater outcrossing even if inbreeding depression has been purged from the contributing populations. Why, then, are selfing‐to‐outcrossing transitions not more common? To evaluate the evolutionary response of mating system to heterosis, we model two monomorphic populations of entirely selfing individuals, introduce a modifier allele that increases the rate of outcrossing and investigate whether the heterosis among populations is sufficient for the modifier to invade and fix. We find that the outcrossing mutation invades for many parameter choices, but it rarely fixes unless populations harbour extremely large unique fixed genetic loads. Reversions to outcrossing become more likely as the load becomes more polygenic, or when the modifier appears on a rare background, such as by dispersal of an outcrossing genotype into a selfing population. More often, the outcrossing mutation instead rises to moderate frequency, which allows recombination in hybrids to produce superior haplotypes that can spread without the mutation's further assistance. The transience of heterosis can therefore explain why secondary contact does not commonly yield selfing‐to‐outcrossing transitions.     

Fine-scale estimation of outcrossing in western redcedar with microsatellite assay of bulked DNA

Western redcedar (Thuja plicata, Cupressaceae) is a self-fertile conifer with a mixed mating system and significant variation for outcrossing among populations. In this paper, we conducted a fine-scale study of mating system variation to identify correlates of outcrossing in natural populations. We examined variation for outcrossing within and among individual trees, and describe a new method to estimate outcrossing using bulked DNA samples. Bulking (assaying DNA tissues from several individuals simultaneously) increases the experimental power without increasing the experimental effort. We sampled 80 trees from four natural populations in southwestern British Columbia. From each tree, we sampled from up to six crown positions (three heights and inner vs outer branches). From each position, two samples of three seedlings each were bulked before DNA extractions. Using four microsatellite loci, we obtained outcrossing rates for each tree and for each of the six crown positions. We found individual tree selfing rates to increase with tree height in all four populations, but selfing rates did not differ among crown positions. The higher selfing rate of larger trees is probably due to their greater proportional contribution to local pollen clouds. Individual tree outcrossing rates ranged from 22 to 100% and the population outcrossing rates from 66 to 78%. Missed alleles due to bulking and the estimation method used both cause a downward bias in outcrossing rates, so that these estimates are probably lower than the actual outcrossing rates. Nevertheless, the trends we observed are not affected by systematic biases of estimation.     

Mitosis, stature and evolution of plant mating systems: low-Phi and high-Phi plants

There is a long-recognized association in plants between small stature and selfing, and large stature and outcrossing. Inbreeding depression is central to several hypotheses for this association, but differences in the evolutionary dynamics of inbreeding depression associated with differences in stature are rarely considered. Here, we propose and test the Phi model of plant mating system evolution, which assumes that the per-generation mutation rate of a plant is a function of the number of mitoses (Phi) that occur from zygote to gamete, and predicts fundamental differences between low-Phi (small-statured) and high-Phi (large-statured) plants in the outcomes of the joint evolution of outcrossing rate and inbreeding depression. Using a large dataset of published population genetic studies of angiosperms and conifers, we compute fitted values of inbreeding depression and deleterious mutation rates for small- and large-statured plants. Consistent with our Phi model, we find that populations of small-statured plants exhibit a range of mating systems, significantly lower mutation rates, and intermediate inbreeding depression, while large-statured plants exhibit very high mutation rates and the maximum inbreeding depression of unity. These results indicate that (i) inbred progeny typically observed in large-statured plant populations are completely lost prior to maturity in nearly all populations; (ii) evolutionary shifts from outcrossing to selfing are generally not possible in large-statured species, rather, large-statured species are more likely to evolve mating systems that avoid selfing such as self-incompatibility and dioecy; (iii) destabilization of the mating system-high selfing rate with high-inbreeding depression-might be a common occurrence in large-statured species; and (iv) large-statured species in fragmented populations might be at higher risk of extinction than previously thought. Our results help to unify and simplify a large and diverse field of research, and serve to emphasize the importance that developmental and genetic constraints play in the evolution of plant mating systems.     

Coevolution of self-fertilization and inbreeding depression. II. Symmetric overdominance in viability

We describe the evolutionary dynamics of a modifier of selfing coevolving with a locus subject to symmetric overdominance in viability under general levels of reduction in pollination success as a consequence of self-fertilization (pollen discounting). Simple models of the evolution of breeding systems that represent inbreeding depression as a constant parameter do not admit the possibility of stable mixed mating systems involving both inbreeding and random mating. Contrary to this expectation, we find that coevolution between a modifier of selfing and a single overdominant locus situated anywhere in the genome can generate evolutionarily attracting mixed mating systems. Two forms of association between the modifier locus and the viability locus promote the evolution of outcrossing. The favored heterozygous genotype at the viability locus develops positive associations with modifier alleles that enhance outcrossing and with the heterozygous genotype at the modifier locus. Associations between outcrossing and high viability evolve immediately upon the introduction of a rare modifier allele, even in the absence of linkage.     

Family-level covariation between parasite resistance and mating system in a hermaphroditic freshwater snail

Genetic compatibility, nonspecific defenses, and environmental effects determine parasite resistance. Host mating system (selfing vs. outcrossing) should be important for parasite resistance because it determines the segregation of alleles at the resistance loci and because inbreeding depression may hamper immune defenses. Individuals of a mixed mating hermaphroditic freshwater snail, Lymnaea ovata, are commonly infected by a digenetic trematode parasite, Echinoparyphium recurvatum. We examined covariation between quantitative resistance to novel parasites and mating system by exposing snail families from four populations that differed by their inbreeding coefficients. We found that resistance was unrelated to inbreeding coefficient of the population, suggesting that the more inbred populations did not carry higher susceptibility load than the less inbred populations. Most of the variation in resistance was expressed among the families within the populations. In the population with the lowest inbreeding coefficient, resistance increased with outcrossing rate of the family, as predicted if selfing had led to inbreeding depression. In the other three populations with higher inbreeding coefficients, resistance was unrelated to outcrossing rate. The results suggest that in populations with higher inbreeding some of the genetic load has been purged, uncoupling the predicted relationship between outcrossing rate and resistance. Snail families also displayed crossing reaction norms for resistance when tested in two environments that presented low and high immune challenge, suggesting that genotype-by-environment interactions are important for parasite resistance.     

Wild sorghum from different eco-geographic regions of Kenya display a mixed mating system

Knowledge of mating systems is required in order to understand the genetic composition and evolutionary potential of plant populations. Outcrossing in a population may co-vary with the ecological and historical factors influencing it. However, literature on the outcrossing rate is limited in terms of wild sorghum species coverage and eco-geographic reference. This study investigated the outcrossing rates in wild sorghum populations from different ecological conditions of Kenya. Twelve wild sorghum populations were collected in four sorghum growing regions. Twenty-four individuals per population were genotyped using six polymorphic simple sequence repeat (SSR) markers to compute their indirect equilibrium estimates of outcrossing rate as well as population structure. In addition, the 12 populations were planted in a field in a randomised block design with five replications. Their progeny (250 individuals per population) were genotyped with the six SSR markers to estimate multi-locus outcrossing rates. Equilibrium estimates of outcrossing rates ranged from 7.0 to 75.0%, while multi-locus outcrossing rates (t _m) ranged from 8.9 to 70.0% with a mean of 49.7%, indicating that wild sorghum exhibits a mixed mating system. The wide range of estimated outcrossing rates in wild sorghum populations indicate that environmental conditions may exist under which fitness is favoured by outcrossing and others under which selfing is more advantageous. The genetic structure of the populations studied is concordant with that expected for a species displaying mixed mating system.     

A novel mating system analysis for modes of self-oriented mating applied to diploid and polyploid arctic Easter daisies (Townsendia hookeri)

We have developed a new model for mating system analysis, which attempts to distinguish among alternative modes of self-oriented mating within populations. This model jointly estimates the rates of outcrossing, selfing, automixis and apomixis, through the use of information in the family structure given by dominant genetic marker data. The method is presented, its statistical properties evaluated, and is applied to three arctic Easter daisy populations, one consisting of diploids, the other two of tetraploids. The tetraploids are predominantly male sterile and reported to be apomictic while the diploids are male fertile. In each Easter daisy population, 10 maternal arrays of six progeny were assayed for amplified fragment length polymorphism markers. Estimates, confirmed with likelihood ratio tests of mating hypotheses, showed apomixis to be predominant in all populations (ca. 70%), but selfing or automixis was moderate (ca. 25%) in tetraploids. It was difficult to distinguish selfing from automixis, and simulations confirm that with even very large sample sizes, the estimates have a very strong negative statistical correlation, for example, they are not independent. No selfing or automixis was apparent in the diploid population, instead, moderate levels of outcrossing were detected (23%). Low but significant levels of outcrossing (2-4%) seemed to occur in the male-sterile tetraploid populations; this may be due to genotyping error of this level. Overall, this study shows apomixis can be partial, and provides evidence for higher levels of inbreeding in polyploids compared to diploids and for significant levels of apomixis in a diploid plant population.     

Mixed mating in natural populations of the chestnut blight fungus, Cryphonectria parasitica

As in plants, fungi exhibit wide variation in reproductive strategies and mating systems. Although most sexually reproducing fungi are either predominantly outcrossing or predominantly selfing, there are some notable exceptions. The haploid, ascomycete chestnut blight pathogen, Cryphonectria parasitica, has previously been shown to have a mixed mating system in one population in USA. In this report, we show that both selfing and outcrossing occur in 10 additional populations of C. parasitica sampled from Japan, Italy, Switzerland and USA. Progeny arrays from each population were assayed for segregation at vegetative incompatibility (vic) and DNA fingerprinting loci. Outcrossing rates (t(m)) were estimated as the proportion of progeny arrays showing segregation at one or more loci, corrected by the probability of nondetection of outcrossing (alpha). Estimates of t(m) varied from 0.74 to 0.97, with the lowest rates consistently detected in USA populations (0.74-0.78). Five populations (four in USA and one in Italy) had t(m) significantly less than 1, supporting the conclusion that these populations exhibit mixed mating. The underlying causes of variation in outcrossing rates among populations of C. parasitica are not known, but we speculate that--as in plants--outcrossing is a function of ecological, demographic and genetic factors.     

An explicit model for the inbreeding load in the evolutionary analysis of selfing

I present analytical predictions for the equilibrium inbreeding load expected in a population under mutation, selection, and a regular mating system for any population size and for any magnitude and recessivity of the deleterious effects. Using this prediction, I deduce the relative fitness of mutant alleles with small effect on selfing to explore the situations where selfing or outcrossing are expected to evolve. The results obtained are in agreement with previous literature, showing that natural selection is expected to lead to stable equilibria where populations show either complete outcrossing or complete selfing, and that selfing is promoted by large deleterious mutation rates. I find that the evolution of selfing is favored by a large recessivity of deleterious effects, while the magnitude of homozygous deleterious effects only becomes relevant in relatively small populations. This result contradicts the standard assumption that purging in large populations will only promote selfing when homozygous deleterious effects are large, and implies that previously published results obtained assuming lethal mutations in large populations can be extrapolated to nonlethal alleles of similar recessivity. This conclusion and the general approach used in this analysis can be useful in the study of the evolution of mating systems.     

The joint evolution of mating system, floral traits and life history in Clarkia (Onagraceae): genetic constraints vs. independent evolution

Genetic correlations caused by pleiotropy or linkage disequilibrium may influence the joint evolution of multiple traits as populations or taxa diverge. The evolutionary transition from outcrossing to selfing has occurred numerous times and is often accompanied by phenotypic and genetic changes in multiple traits such as flower size, pollen-ovule ratio, stigma and anther maturity and the age of reproductive maturity. Determining whether the recurring patterns of multitrait change are because of selection on each trait independently and/or the result of genetic correlations among traits can shed light on the mechanism that accounts for such convergence. Here, we evaluate whether floral traits are genetically correlated with each other and/or with whole-plant traits within- and between-populations and taxa. We report results from a greenhouse study conducted on two pairs of sister taxa with contrasting mating systems: the autogamously selfing Clarkia exilis and its predominantly outcrossing progenitor C. unguiculata and the autogamous Clarkia xantiana ssp. parviflora and its outcrossing progenitor C. xantiana ssp. xantiana. We examined variation within and covariation among maternal families in three populations of each taxon with respect to the age at first flower, the rate of successive flower production and the number of days between bud break and anther dehiscence and stigma receptivity within individual flowers. Based on phenotypic divergence between sister taxa, bivariate regressions, correlations among maternal family means and analysis of covariance (ancova), we did not find unilateral support indicating that genetic constraints govern the joint distribution of floral and whole-plant traits.     

Does cleistogamy variation translate into outcrossing variation in the annual species Lamium amplexicaule (Lamiaceae)?

The maintenance of cleistogamy, the ability to produce closed, obligately selfing flowers (CL), and open, potentially outcrossed flowers (CH), in different proportions, is classically explained through different morphological/physiological properties of the two floral types, but rarely as a mechanism of adjusting the outcrossing rate. We explore the link between CH proportion and overall outcrossing rate in natural populations of Lamium amplexicaule. We assessed number of calices, CH proportion, CH and overall outcrossing rate in four natural populations in two distant areas in France. In each region, we had one favorable and one unfavorable habitat population. Unfavorable habitats produce smaller plants (with fewer calices) with higher CH proportions compared to favorable habitats, regardless of the geographic origin of the populations. CH outcrossing rate did not change significantly among populations. Thus, the overall outcrossing rate in L. amplexicaule is mainly determined by the CH proportion. Contrary to the classical view, unfavorable environments in our study are associated with higher rate of chasmogamous flowers, supposedly more costly to produce. We propose that cleistogamy variation can be considered as a variation of the outcrossing rate and could be explained by classic forces driving the evolution of mating systems (inbreeding depression, pollinators’ abundance).     

Molecular Evolution of the GapC Gene Family in Amsinckia spectabilis Populations That Differ in Outcrossing Rate

Molecular evolutionary analysis of the glyceraldehyde 3-phosphate dehydrogenase (GapC) gene family was conducted in the plant genus Amsinckia (Boraginaceae), a group that exhibits marked variation in the mating system. GapC genes in this group differ from those of Arabidopsis thaliana in terms of both intron size and number. Phylogenetic and Southern hybridization analyses suggest the presence of multiple GapC loci, each defined by a set of base substitutions that are in strong linkage disequilibrium. One species of Amsinckia, A. spectabilis, was studied in some detail. This species consists of selfing (A. s. spectabilis) and outcrossing (A. s. microcarpa) varieties. Two selfing populations and one outcrossing population sample were analyzed in detail for variation at one of the members of this gene family, GapC3. A reduction in number of GapC3 haplotypes and level of genetic diversity was observed in the selfing populations of A. spectabilis. GapC3 in the outcrossing population (but not the two selfing populations) exhibited a significant departure from neutrality in the direction of an excess of singletons. These results are discussed in the context of forces acting on sequence evolution in populations with different mating systems.[Reviewing Editor: Dr. Martin Kreitman]     

Selfing versus outcrossing propensity of the fungal pathogen Microbotryumviolaceum across Silenelatifolia host plants

In the fungal pathogen Microbotryumviolaceum mating (i.e. conjugation between cells of opposite mating type) is indispensable for infection of its host plant Silenelatifolia. Since outcrossing opportunities are potentially rare, selfing may be appropriate to ensure reproduction. On the other hand, outcrossing may create genetic variability necessary in the coevolutionary arms race with its host. We investigated the propensity of M. violaceum to outcross vs. self in different host environments. We used haploid sporidia from each of three strains from five fungal populations for pairwise mixtures of opposite mating type, representing either selfing or outcrossing combinations. Mixtures were exposed to leaf extract from seven S. latifolia plants. The proportion of conjugated sporidia quantified mating propensity. The identity of both fungal strains and host influenced conjugation. First, individual strains differed in conjugation frequency by up to 30%, and strains differed in their performance across the different hosts. Second, selfing combinations produced, on average, more conjugations than did outcrossing combinations. Selfing appears to be the predominant mode of reproduction in this fungus, and selfing preference may have evolved as a mechanism of reproductive assurance. Third, individual strains varied considerably in conjugation frequency in selfing and outcrossing combinations across different hosts. This indicates that conjugation between outcrossing partners could be favoured at least in some hosts. Since the dikaryon resulting from conjugation is the infectious unit, conjugation frequency may correspond with infection probability. This assumption was supported by an inoculation experiment, where high infectious sporidial dosage resulted in higher infections success than did low dosage. We therefore predict that sexual recombination can provide this pathogen with novel genotypes able to infect local resistant hosts.     

Inbreeding depression is high in a self‐incompatible perennial herb population but absent in a self‐compatible population showing mixed mating

High inbreeding depression is thought to be one of the major factors preventing evolutionary transitions in hermaphroditic plants from self‐incompatibility (SI) and outcrossing toward self‐compatibility (SC) and selfing. However, when selfing does evolve, inbreeding depression can be quickly purged, allowing the evolution of complete self‐fertilization. In contrast, populations that show intermediate selfing rates (a mixed‐mating system) typically show levels of inbreeding depression similar to those in outcrossing species, suggesting that selection against inbreeding might be responsible for preventing the transition toward complete self‐fertilization. By implication, crosses among populations should reveal patterns of heterosis for mixed‐mating populations that are similar to those expected for outcrossing populations. Using hand‐pollination crosses, we compared levels of inbreeding depression and heterosis between populations of Linaria cavanillesii (Plantaginaceae), a perennial herb showing contrasting mating systems. The SI population showed high inbreeding depression, whereas the SC population displaying mixed mating showed no inbreeding depression. In contrast, we found that heterosis based on between‐population crosses was similar for SI and SC populations. Our results are consistent with the rapid purging of inbreeding depression in the derived SC population, despite the persistence of mixed mating. However, the maintenance of outcrossing after a transition to SC is inconsistent with the prediction that populations that have purged their inbreeding depression should evolve toward complete selfing, suggesting that the transition to SC in L. cavanillesii has been recent. SC in L. cavanillesii thus exemplifies a situation in which the mating system is likely not at an equilibrium with inbreeding depression.     

Host mating system and the prevalence of disease in a plant population

A modified susceptible-infected-recovered (SIR) host-pathogen model is used to determine the influence of plant mating system on the outcome of a host-pathogen interaction. Unlike previous models describing how interactions between mating system and pathogen infection affect individual fitness, this model considers the potential consequences of varying mating systems on the prevalence of resistance alleles and disease within the population. If a single allele for disease resistance is sufficient to confer complete resistance in an individual and if both homozygote and heterozygote resistant individuals have the same mean birth and death rates, then, for any parameter set, the selfing rate does not affect the proportions of resistant, susceptible or infected individuals at equilibrium. If homozygote and heterozygote individual birth rates differ, however, the mating system can make a difference in these proportions. In that case, depending on other parameters, increased selfing can either increase or decrease the rate of infection in the population. Results from this model also predict higher frequencies of resistance alleles in predominantly selfing compared to predominantly outcrossing populations for most model conditions. In populations that have higher selfing rates, the resistance alleles are concentrated in homozygotes, whereas in more outcrossing populations, there are more resistant heterozygotes.