Plant Mating Systems Research Articles

Pollination insights for the conservation of a rare threatened plant species, Astragalus tragacantha (Fabaceae)

The increase in habitat fragmentation impacts plant-pollinator interactions and threatens the sustainability of plant species. Astragalus tragacantha (Fabaceae), is a rare endangered plant species along the coastal habitats where the plant populations have undergone considerable fragmentation and decline of size. Controlled pollination treatments, the observation of pollinator activity, and pollinator captures, have been conducted to study: (1) the mating system of A. tragacantha and the potential for inbreeding depression and/or outbreeding depression based on controlled pollination treatments, (2) the pollinator composition among populations using a correspondence analysis and a hierarchical clustering, and (3) the link between pollinators and the plant reproductive success using a path-analysis model. In this study, we demonstrated that this plant was not autogamous self-pollinating and depended on pollinators for its reproduction. The absence of difference between manual and open pollinations regarding the reproductive success showed an absence of pollen limitation in our populations. We showed that populations differed in the composition of their pollinator guilds. Some pollinator species were predominant in certain populations. The pollination treatments revealed the existence of a mixed mating system in A. tragacantha populations. We showed an inbreeding depression potentially linked to a predominant pollinator-facilitated selfing, and the existence of outbreeding depression between some distant populations. These differences in pollinator guild and plant mating systems among populations must be considered during the restoration of populations along the Mediterranean coastal habitats in order to enhance the reproductive success and sustainability of A. tragacantha.

The impact of floral morphology on genetic differentiation in two closely related biennial plant species.

The genetic diversity and structure of plant populations are determined by the interaction of three distinct processes: gene flow, genetic drift and natural selection. These processes are to some extent dependent on the mating system of plants, which in turn is largely determined by floral morphology and the level of herkogamy in particular. In this study, we used molecular markers to investigate the impact of floral morphology on genetic differentiation and structure in two closely related Centaurium species that display large variation in floral morphology across two distinct geographic regions in Europe (mainland Europe and the UK). Our results showed that genetic differences between regions and populations within regions were similar for both species, but that patterns of genetic structure largely depended on the observed variation in floral morphology. Populations of Centaurium erythraea showed higher genetic differentiation and clear isolation by distance (IBD) in mainland Europe, but limited IBD in the UK. Opposite patterns were found in Centaurium littorale, with higher genetic differentiation and significant IBD in populations sampled in the UK and lower genetic differentiation in Continental populations with no pattern of IBD. Overall, these results indicate that variation in floral morphology has a profound impact on structuring of genetic diversity, with populations displaying low levels of herkogamy showing the strongest patterns of genetic structuring and the reverse pattern in populations showing high levels of herkogamy.

Does the evolution of self-fertilization rescue populations or increase the risk of extinction?

As a major evolutionary transition in seed plants, the evolution of plant mating systems has been much debated in evolutionary ecology. Over the last 10 years, well-established patterns of evolution have emerged. On the one hand, experimental studies have shown that self-fertilization is likely to evolve in a few generations (microevolution) as a response to rapid environmental change (e.g. pollinator decline), eventually rescuing a population. On the other, phylogenetic studies have demonstrated that repeated evolution towards self-fertilization (macroevolution) leads to a higher risk of lineage extinction and is thus likely to be disadvantageous in the long term. In either case - the short-term or long-term evolution of self-fertilization (selfing) - these findings indicate that a mating system is not neutral with respect to population or lineage persistence. They also suggest that selfing can have contrasting effects depending on time scale. This raises the question of whether mating system evolution can rescue populations facing environmental change. In this review, empirical and theoretical evidence of the direct and indirect effects of mating systems on population demography and lineage persistence were analysed. A simple theoretical evolutionary rescue model was also developed to investigate the potential for evolutionary rescue through selfing. Demographic studies consistently show a short-term advantage of selfing provided by reproductive assurance, but a long-term disadvantage for selfing lineages, suggesting indirect genomic consequences of selfing (e.g. mutation load and lower adaptability). However, our theoretical evolutionary rescue model found that even in the short term, while mating system evolution can lead to evolutionary rescue, it can also lead to evolutionary suicide, due to the inherent frequency-dependent selection of mating system traits. These findings point to the importance of analysing the demographic consequences of self-fertilization in order to predict the effect of selfing on population persistence as well as take into account the indirect genomic consequences of selfing. The pace at which processes such as inbreeding depression, purging, reproductive assurance and genomic rearrangements occur after the selfing transition is the key to clarifying whether or not selfing will result in evolutionary rescue.

Population density and size influence pollen dispersal pattern and mating system of the predominantly outcrossed Banksia nivea (Proteaceae) in a threatened ecological community

Gene flow is a critical component of plant mating systems and influences population fitness. However, pollen dispersal can be highly variable and influenced by natural and anthropogenic fragmentation. Gene flow through pollen dispersal was investigated in two populations of contrasting size and habitat context in Banksia nivea ssp. uliginosa, a rare species in the Busselton ironstone threatened ecological community of Western Australia with a naturally fragmented distribution. Paternity analysis was conducted with seven microsatellite loci to determine mating system parameters and patterns of pollen dispersal. Outcrossing was high in both populations with a similar level of selfing for both populations despite differences in population size, density and vegetation matrix. Most mating occurred within 10 m of a mother plant in the small, clumped population, while more dispersed mating, up to 50 m from a mother plant, was recorded in the large, less dense population. Our results show that population density and size are important influences on mating system parameters and level of pollen dispersal.

Demography and mating system shape the genome-wide impact of purifying selection in Arabis alpina

Plant mating systems have profound effects on levels and structuring of genetic variation and can affect the impact of natural selection. Although theory predicts that intermediate outcrossing rates may allow plants to prevent accumulation of deleterious alleles, few studies have empirically tested this prediction using genomic data. Here, we study the effect of mating system on purifying selection by conducting population-genomic analyses on whole-genome resequencing data from 38 European individuals of the arctic-alpine crucifer Arabis alpina We find that outcrossing and mixed-mating populations maintain genetic diversity at similar levels, whereas highly self-fertilizing Scandinavian A. alpina show a strong reduction in genetic diversity, most likely as a result of a postglacial colonization bottleneck. We further find evidence for accumulation of genetic load in highly self-fertilizing populations, whereas the genome-wide impact of purifying selection does not differ greatly between mixed-mating and outcrossing populations. Our results demonstrate that intermediate levels of outcrossing may allow efficient selection against harmful alleles, whereas demographic effects can be important for relaxed purifying selection in highly selfing populations. Thus, mating system and demography shape the impact of purifying selection on genomic variation in A. alpina These results are important for an improved understanding of the evolutionary consequences of mating system variation and the maintenance of mixed-mating strategies.

Sequence diversification and phylogenetic analysis of selfincompatibility specific determinant genes in Brassicaceae

It is a long history of theoretical interest in understanding the genetic causes and consequences of shifts between outcrossing and inbreeding in plant mating systems. Selfincompatibility (SI) systems genetically promote outcrossing, and its transitions to inbreeding potentially changes genetic diversity in the evolution of species. SI in the Brassicaceae family is sporophytically controlled by a single locus (S-locus), which locates the female determinant gene SRK (S-locus receptor kinase) and the male determinant gene SCR (S-locus cysteine rich). Operation of diversifying selection maintains the diversity of SRK and SCR during the evolutionary course of S haplotypes. Here we used deduced amino acid sequences of SRK and/or SCR publicly available from seven species in three genera, Brassica, Arabidopsis and Raphanus, to investigate the sequence diversification and phylogenetic distribution of S haplotypes in the Brassicaceae family. The results reveal that SCR shows higher sequence distance than SRK within or between genera, and also in the interspecific pairs of S haplotypes. Phylogenetic analysis supports the conclusion that SI has a single origin but evolves differently within the Brassicaceae. The genetic mechanisms and consequences of the loss of SI in the Brassicaceae family are discussed.

Research Article Isolation and cross-amplification of fifteen microsatellites for oil-rewarding Calceolaria species (Calceolariaceae)

Calceolaria is a diversified American genus whose species are mainly visited by oil-collecting bees, representing a highly specialized plant-pollinator system. However, although this is a unique plant system, different aspects of its ecology and evolution remain unexplored. In this study we characterized fifteen polymorphic microsatellites for four Andean Calceolaria species using next-generation sequencing. The number of alleles per locus (Na) for the study species ranged from 2 to 11, and the observed (Ho) and expected (He) heterozygosity ranged from 0 to 0.85 and 0.3 to 0.87, respectively, while FIS ranged from 0.03 to 1.0. Among the four species, C. cana showed the greatest genetic diversity, followed by C. filicaulis. C. arachnoidea and C. lanigera showed similar, lower genetic diversity. These makers are a useful tool to conduct futures studies on population structure and gene flow. Specifically, this set of markers will be helpful to estimate selfing rates as a first approximation to the understanding of the role of specialization in plant mating system evolution

Aggregate enantiostyly: Floral visitor interactions with a previously unreported form of floral display

Floral herkogamy, or the spatial separation of reproductive structures, plays an important role in the evolution of plant mating systems. One form of herkogamy is enantiostyly, or the alternate presentation of stigmas to the left and right of the floral axis. Theoretically, enantiostyly increases outcrossing rate by promoting reciprocal pollen deposition and export, but empirical evidence of this is lacking. Understanding how pollinators interact with enantiostylous flowers will illuminate the mechanisms by which enantiostyly encourages outcross pollination. We quantified style deployment within and among plants, of Allionia incarnata L. (Nyctaginaceae), a species with a new kind of enantiostyly we call “aggregate enantiostyly”. A. incarnata is the first species known to have “compound blossoms”, with different flowers each containing a left- or right-handed style. We quantify the autonomous selfing rate and describe floral visitor interactions with this unique floral structure. We found equal frequency of left- and right-handed blossoms within and among plants. Common floral visitors, Hemiargus ceraunus and Halictus sp., interacted with either anthers or stigmas, but not both, indicating effective outcrossing at the flower level. Outcrossing produced more fruits than selfing. The order in which floral visitors interacted with reproductive structures indicates that enantiostyly may increase pollen export. However, we did not find evidence that enantiostyly affected pollen deposition. Our study provides an intriguing first look at the reproductive ecology of a species with a new unique form of herkogamy, and provides evidence on the adaptive significance of enantiostyly.

Changes in floral biology and inbreeding depression in native and invaded regions of Datura stramonium.

Plant populations invading new environments might compromise their fitness contribution to the next generation, because of the lack of native specialist pollinators and/or potential mates. Thus, changes in plant mating system and traits linked to it are expected in populations colonising new environments where selection would favour selfing and floral traits that maximise reproductive output. To test this, we studied native (Mexico) and non-native (Spain) populations of the obligate sexual reproducing annual weed Datura stramonium. Flower size, herkogamy, total number of seeds per plant, number of visits by and type of pollinators, and inbreeding depression were assessed in native and non-native populations. Finally, we measured phenotypic selection on corolla size and herkogamy in each population. Flower size and herkogamy showed wide and similar variation in both ranges. However, the largest average flower size was found in one non-native population whereas the highest average positive herkogamy was detected in one native population. On average, flowers in the native range received more visits by pollinators. Hawkmoths were the main visitors in the native populations while only bees were observed visiting flowers in Spain's populations. Only in the native range was inbreeding depression detected. Selection to reduce herkogamy was found only in one native population. Absence of both inbreeding depression and selection on floral traits suggest a change in mating system of D.stramonium in a new range where generalist pollinators may be promoting high reproductive success. Selection against deleterious alleles might explain the reduction of inbreeding depression, promoting the evolution of selfing.

The Ecology of Mating and Its Evolutionary Consequences in Seed Plants

Mating in seed plants arises from interactions between plant traits and the environmental and demographic context in which individuals reside. These interactions commonly cause nonrandom mating, including selfing and promiscuous outcrossing within local neighborhoods. Shared features of seed plants, specifically immobility, hermaphroditism, and modularity, shape the essential character of mating mediated by animals, wind, and water. In addition, diverse floral strategies promote cross- and self-mating, depending on environmental circumstances. Extrinsic ecological factors influence all stages of the mating process—pollination, pollen-tube growth, ovule fertilization—as well as seed development, determining offspring quantity and quality. Traditionally, measures of plant mating systems have focused on a single axis of variation, the maternal outcrossing rate. Instead, we argue for an expanded perspective encompassing mating portfolios, which include all offspring to which individuals contribute genetically as maternal or paternal parents. This approach should expose key ecological determinants of mating-system variation and their evolutionary consequences.

NMπ-improved re-implementation of NM+, a software for estimating gene dispersal and mating patterns.

This study introduces the NMπ computer program designed for estimation of plant mating system and seed and pollen dispersal kernels. NMπ is a re-implementation of the NM+ program and provides new features such as support for multicore processors, explicit treatment of dioecy, the possibility of incorporating uniparentally cytoplasmic markers, the possibility of assessing assortative mating due to phenotypic similarity and inference about offspring genealogies. The probability model of parentage (the neighbourhood model) accounts for missing data and genotyping errors, which can be estimated along with regular parameters of the mating system. The program has virtually no restrictions with respect to a number of individuals, markers or phenotypic characters. A console version of NMπ can be run under a wide variety of operating systems, including Windows, Linux or Mac OS. For Windows users, a graphical user interface is provided to facilitate operating the software. The program, user manual and example data are available on http://www.ukw.edu.pl/pracownicy/plik/igor_chybicki/3694/.

Flower-level developmental plasticity to nutrient availability in Datura stramonium: implications for the mating system.

Studies of phenotypic plasticity in plants have mainly focused on (1) the effect of environmental variation on whole-plant traits related to the number of modules rather than on (2) the phenotypic consequences of environmental variation in traits of individual modules. Since environmental and developmental factors can produce changes in traits related to the mating system, this study used the second approach to investigate whether within-individual variation in herkogamy-related traits is affected by the environment during plant development in two populations of Datura stramonium , an annual herb with a hypothesized persistent mixed mating system, and to determine which morphological traits may promote self-fertilization. Full-sib families of two Mexican populations of D. stramonium , with contrasting ecological histories, were grown under low, mid and high nutrient availability to investigate the effects of genetic, environmental and within-plant flower position on flower size, corolla, stamen and pistil lengths, and herkogamy. Populations showed differences in familial variation, plasticity and familial differences in plasticity in most floral traits analysed. In one population (Ticumán), the effect of flower position on trait variation varied among families, whereas in the other (Pedregal) the effect of flower position interacted with the nutrient environment. Flower size varied with the position of flowers, but in the opposite direction between populations in low nutrients; a systematic within-plant trend of reduction in flower size, pistil length and herkogamy with flower position increased the probability of self-fertilization in the Pedregal population. Besides genetic variation in floral traits between and within populations, environmental variation affects phenotypic floral trait values at the whole-plant level, as well as among flower positions. The interaction between flower position and nutrient environment can affect the plant's mating system, and this differs between populations. Thus, reductions in herkogamy with flower positions may be expected in environments with either low pollinator abundance or low nutrients.

Novel Consequences of Bird Pollination for Plant Mating.

Pollinator behaviour has profound effects on plant mating. Pollinators are predicted to minimise energetic costs during foraging bouts by moving between nearby flowers. However, a review of plant mating system studies reveals a mismatch between behavioural predictions and pollen-mediated gene dispersal in bird-pollinated plants. Paternal diversity of these plants is twice that of plants pollinated solely by insects. Comparison with the behaviour of other pollinator groups suggests that birds promote pollen dispersal through a combination of high mobility, limited grooming, and intra- and interspecies aggression. Future opportunities to test these predictions include seed paternity assignment following pollinator exclusion experiments, single pollen grain genotyping, new tracking technologies for small pollinators, and motion-triggered cameras and ethological experimentation for quantifying pollinator behaviour.

Floral traits and mating system of Hibiscus trionum (Malvaceae)

Variations in floral traits and floral structures influence plant mating systems. Hibiscus trionum produces large, showy flowers typical of an outcrossing species, yet flowers are autonomously self-pollinated. In this study, we measured floral morphology, breeding system and outcrossing rate estimated by ISSR markers. Results indicate that two types of flowers were observed in H. trionum, and the type I with bigger petals appears to be much more visible to pollinators, demonstrated by than type II flowers with smaller petals. The flowers with hand pollination were closed 1h earlier than intact flowers, whether they were type I or II. The relationship between the amount of pollen deposited on the stigma and the number of seeds per capsule was highly significant, and 80 or more pollens per flower can make the mean number of seeds (mean=37) in H. trionum. Delayed selfing in H. trionum did not provide a large contribution to seed production, since reproductive assurance were only 0.025. However, successful reproduction of 72.5% flowers in the absence of pollinators suggested that selfing provides reproductive assurance during seasons, in which pollinators were limiting. The multilocus outcrossing rates in different populations varied from 0.982 to 1.200, with a mean of 1.116. Our data provide an empirical demonstration of a predominantly outcrossing species with potential delayed selfing when pollinators are absent or scarce.

Quantification of the zygotic barrier between interbreeding taxa using gene flow data

Hybridization and introgression via interspecific gene flow are common processes in the plant kingdom. The effectiveness of these processes is governed by the strengths of multiple zygotic barriers. These barriers have often been quantified in artificial settings using laborious and time‐consuming hand‐pollination experiments, but their quantification is nonexistent in the landscape. In this study, we utilized gene flow data within a spatially explicit simulation to assess the strengths of zygotic barriers. Our model system consisted of Populus nigra and its hybrid, P. × canadensis, which interbreed under natural conditions. The study population was located in the floodplain of the Eder River in central Germany. Pollen‐mediated introgression rates from hybrid males into the seeds of individual female trees were used as the target pattern using an inverse modeling approach. Simulations that treated pollen from both taxa equally revealed a large discrepancy between the observed and modeled rates of introgression for both taxa. The discrepancy was reduced by introducing a zygotic barrier against the pollen from the hybrid males. The best model outcome values indicated comparably strong zygotic barriers acting against pollen‐mediated introgressive gene flow into the two parental taxa, P. nigra and P. × canadensis. The sensitivity of our model was tested by applying different dispersal functions. Four common probability density functions were used along with a pollen dispersal function that had previously been fitted to gene flow data from the same dataset. The best barrier value was almost independent of the dispersal functions used here. Moreover, it was within the range previously determined in hand‐pollination‐based investigations, validating our model. These data indicate that the inverse modeling approach is a powerful method for quantifying hidden processes, and we discuss its use as a valuable tool for generating new insights into plant mating systems that are relevant to evolutionary biology and risk analyses in conservation efforts.

MYB transcription factor gene involved in sex determination in Asparagus officinalis.

Dioecy is a plant mating system in which individuals of a species are either male or female. Although many flowering plants evolved independently from hermaphroditism to dioecy, the molecular mechanism underlying this transition remains largely unknown. Sex determination in the dioecious plant Asparagus officinalis is controlled by X and Y chromosomes; the male and female karyotypes are XY and XX, respectively. Transcriptome analysis of A.officinalis buds showed that a MYB-like gene, Male Specific Expression 1 (MSE1), is specifically expressed in males. MSE1 exhibits tight linkage with the Y chromosome, specific expression in early anther development and loss of function on the X chromosome. Knockout of the MSE1 orthologue in Arabidopsis induces male sterility. Thus, MSE1 acts in sex determination in A.officinalis.

High Correlated Paternity Leads to Negative Effects on Progeny Performance in Two Mediterranean Shrub Species.

Anthropogenic habitat deterioration can promote changes in plant mating systems that subsequently may affect progeny performance, thereby conditioning plant recruitment for the next generation. However, very few studies yet tested mating system parameters other than outcrossing rates; and the direct effects of the genetic diversity of the pollen received by maternal plants (i.e. correlated paternity) has often been overlooked. In this study, we investigated the relation between correlated paternity and progeny performance in two common Mediterranean shrubs, Myrtus communis and Pistacia lentiscus. To do so, we collected open-pollinated progeny from selected maternal plants, calculated mating system parameters using microsatellite genotyping and conducted sowing experiments under greenhouse and field conditions. Our results showed that some progeny fitness components were negatively affected by the high correlated paternity of maternal plants. In Myrtus communis, high correlated paternity had a negative effect on the proportion and timing of seedling emergence in the natural field conditions and in the greenhouse sowing experiment, respectively. In Pistacia lentiscus, seedling emergence time under field conditions was also negatively influenced by high correlated paternity and a progeny survival analysis in the field experiment showed greater mortality of seedlings from maternal plants with high correlated paternity. Overall, we found effects of correlated paternity on the progeny performance of Myrtus communis, a self-compatible species. Further, we also detected effects of correlated paternity on the progeny emergence time and survival in Pistacia lentiscus, an obligate outcrossed species. This study represents one of the few existing empirical examples which highlight the influence that correlated paternity may exert on progeny performance in multiple stages during early seedling growth.

The evolution of optimal resource allocation and mating systems in hermaphroditic perennial plants.

By incorporating the effects of inbreeding depression (ID) on both juveniles and adults survivorship, we developed a new theoretical model for hermaphroditic perennial plants. Our model showed that the effect of the selfing rate on the evolutionarily stable strategy (ESS) reproductive allocation depends on three parameters: (1) the self-fertilized juvenile relative survivorship (SFJRS), (2) the self-fertilized adult relative survivorship (SFARS) and (3) the growth rate of self-fertilized adult, where the SFJRS is the survivorship of self-fertilized juveniles divided by the survivorship of outcrossed juveniles, and likewise for the SFARS. However, the ESS sex allocation decreases as the selfing rate increases. This relationship seems independent of the SFJRS, the SFARS, and the growth rate of self-fertilized adults. Additionally, our model showed that the complete outcrossing is an ESS when the fraction of juvenile inbreeding depression (FJID) is less than 1/2 − τ, where τ is the self-fertilized adults mortality rate caused by ID. In contrast, the complete selfing also acts as an ESS when the FJID is greater than 1/2 − τ. These results could explain the diversity of mating strategies and related resource allocations for plants.

Partial genotyping at polymorphic markers can improve heritability estimates in sibling groups

Accurate estimates of heritability (h2) are necessary to assess adaptive responses of populations and evolution of fitness-related traits in changing environments. For plants, h2 estimates generally rely on maternal progeny designs, assuming that offspring are either half-sibs or unrelated. However, plant mating systems often depart from half-sib assumptions, this can bias h2 estimates. Here, we investigate how to accurately estimate h2 in nonmodel species through the analysis of sibling designs with a moderate genotyping effort. We performed simulations to investigate how microsatellite marker information available for only a subset of offspring can improve h2 estimates based on maternal progeny designs in the presence of nonrandom mating, inbreeding in the parental population or maternal effects. We compared the basic family method, considering or not adjustments based on average relatedness coefficients, and methods based on the animal model. The animal model was used with average relatedness information, or with hybrid relatedness information: associating one-generation pedigree and family assumptions, or associating one-generation pedigree and average relatedness coefficients. Our results highlighted that methods using marker-based relatedness coefficients performed as well as pedigree-based methods in the presence of nonrandom mating (i.e. unequal male reproductive contributions, selfing), offering promising prospects to investigate insitu heritabilities in natural populations. In the presence of maternal effects, only the use of pairwise relatednesses through pedigree information improved the accuracy of h2 estimates. In that case, the amount of father-related offspring in the sibling design is the most critical. Overall, we showed that the method using both one-generation pedigree and average relatedness coefficients was the most robust to various ecological scenarios.

Inbreeding depression under mixed outcrossing, self-fertilization and sib-mating.

BackgroundBiparental inbreeding, mating between two relatives, occurs at a low frequency in many natural plant populations, which also often have substantial rates of self-fertilization. Although biparental inbreeding is likely to influence the dynamics of inbreeding depression and the evolution of selfing rates, it has received limited theoretical attention in comparison to selfing. The only previous model suggested that biparental inbreeding can favour the maintenance of stable intermediate selfing rates, but made unrealistic assumptions about the genetic basis of inbreeding depression. Here we extend a genetic model of inbreeding depression, describing nearly recessive lethal mutations at a very large number of loci, to incorporate sib-mating. We also include a constant component of inbreeding depression modelling the effects of mildly deleterious, nearly additive alleles. We analyze how observed rates of sib-mating influence the mean number of heterozygous lethals alleles and inbreeding depression in a population reproducing by a mixture of self-fertilization, sib-mating and outcrossing. We finally use the ensuing relationship between equilibrium inbreeding depression and population selfing rate to infer the evolutionarily stable selfing rates expected under such a mixed mating system.ResultsWe show that for a given rate of inbreeding, sib-mating is more efficient at purging inbreeding depression than selfing, because homozygosity of lethals increases more gradually through sib-mating than through selfing. Because sib-mating promotes the purging of inbreeding depression and the evolution of selfing, our genetic model of inbreeding depression also predicts that sib-mating is unlikely to maintain stable intermediate selfing rates.ConclusionsOur results imply that even low rates of sib-mating affect plant mating system evolution, by facilitating the evolution of selfing via more efficient purging of inbreeding depression. Alternative mechanisms, such as pollination ecology, are necessary to explain stable mixed selfing and outcrossing.