Multiple origins of self-fertilization in tristylous Eichhornia paniculata (Pontederiaceae): Inferences from style morph and isozyme variation

Populations of Eichhornia paniculata (Pontederiaceae) exhibit a wide range of mating systems, from predominant outcrossing to high levels of self-fertilization. The origin of self-fertilization in this tristylous species is associated with the loss of style-length morphs from populations and the spread of self-pollinating, floral variants. We examined geographic variation in style morph and allozyme frequencies to determine whether the loss of style morphs and transition to selfing could have multiple origins in E. paniculata. Surveys of floral variation in 167 populations from six states in northeastern Brazil revealed that at least one style morph was absent from 29.3%. Non-trimorphic populations occurred in all states and ranged in frequency from 9% in Ceará to 68% in Alagoas. Selfing variants occurred in 8.5% and 55% of trimorphic and non-trimorphic populations, respectively, and were distributed among five of six states with primary concentrations in Alagoas and Pernambuco. A comparison of electrophoretic variation at 24 isozyme loci in 28 trimorphic, 13 dimorphic and 3 monomorphic populations indicated that non-trimorphic populations contained 84% of the allelic variation present in trimorphic populations and were markedly differentiated from one another. Analyses of genetic distance and the distribution of rare alleles indicated that non-trimorphic populations were often more similar to neighbouring trimorphic populations than to one another. Populations with selfing variants occurred at low frequency in three genetically distinct parts of the range. These results, in combination with genetic and morphological evidence suggest that style morphs are lost repeatedly from populations of E. paniculata and that selfing variants may have originated on at least three separate occasions in northeastern Brazil.     

THE DISSOLUTION OF A COMPLEX GENETIC POLYMORPHISM: THE EVOLUTION OF SELF-FERTILIZATION IN TRISTYLOUS EICHHORNIA PANICULATA (PONTEDERIACEAE)

Eichhornia paniculata (Pontederiaceae) displays a wide range of outcrossing levels as a result of the dissolution of the tristylous genetic polymorphism and the evolution of semihomostyly. Population surveys, comparison of fitness components of the style morphs, and computer simulations were used to investigate the breakdown of tristyly and the selective mechanisms responsible for the evolution of self-fertilization. Of 110 populations surveyed in northeast Brazil and Jamaica, 53% were trimorphic, 25% were dimorphic, and 22% were monomorphic for style morph. The short (S) morph was underrepresented in trimorphic populations and absent from nontrimorphic populations. The mid (M) morph predominated in dimorphic populations and was the only morph in monomorphic populations. Stamen modifications promoting selfing, associated with semihomostyle evolution, were largely confined to the M morph. They were rare in trimorphic populations, common in dimorphic populations, and often fixed in monomorphic populations. Stochastic simulations and comparisons of fruit set in natural populations indicate that founder events, population bottlenecks, and lowered fertility of the S morph due to an absence of long-tongued pollinators can each account for loss of the S morph from trimorphic populations. A reduced level of disassortative mating can accentuate the rate at which the S morph is lost by both random and deterministic processes. Nontrimorphic populations occur at the geographical margins of the region surveyed and tend to be smaller and less dense than trimorphic populations. These observations and the higher fruit set of the M morph relative to the L morph in dimorphic populations suggest that reproductive assurance, favoring selfing variants of the M morph under conditions of low pollinator service, has been of primary importance in the origin of most monomorphic populations. Where pollinator service is reliable, however, automatic selection of selfing genes, aided by mating asymmetries between the morphs, can cause the M morph to spread to fixation in dimorphic populations.     

POLLEN DISCOUNTING AND THE SPREAD OF A SELFING VARIANT IN TRISTYLOUS EICHHORNIA PANICULATA: EVIDENCE FROM EXPERIMENTAL POPULATIONS

Floral traits that increase self-fertilization are expected to spread unless countered by the effects of inbreeding depression, pollen discounting (reduced outcross pollen success by individuals with increased rates of self-fertilization), or both. Few studies have attempted to measure pollen discounting because to do so requires estimating the male outcrossing success of plants that differ in selfing rate. In natural populations of tristylous Eichhornia paniculata, selfing variants of the mid-styled morph are usually absent from populations containing all three style morphs but often predominate in nontrimorphic populations. We used experimental garden populations of genetically marked plants to investigate whether the effects of population morph structure on relative gamete transmission by unmodified (M) and selfing variants (M‘) of the mid-styled morph could explain their observed distribution. Transmission through ovules and self and outcross pollen by plants of the M and M’ morphs were compared under trimorphic, dimorphic (S morph absent), and monomorphic (L and S morphs absent) population structures. Neither population structure nor floral morphology affected female reproductive success, but both had strong effects on the relative transmission of male gametes. The frequency of self-fertilization in the M' morph was consistently higher than that of the M morph under all morph structures, and the frequency of self-fertilization by both morphs increased as morph diversity of experimental populations declined. In trimorphic populations, total transmission by the M and M' morphs did not differ. The small, nonsignificant increase in selfing by the M' relative to the M morph was balanced by decreased outcross siring success, particularly on the S morph. In populations lacking the S morph, male gamete transmission by the M' morph was approximately 1.5 times greater than that by the M morph because of both increased selfing and increased success through outcross pollen donation. Therefore, gamete transmission strongly favored the M' morph only in the absence of the S morph, a result consistent with the distribution of the M' morph in nature. This study indicates that floral traits that alter the selfing rate can have large and context-dependent influences on outcross pollen donation.     

Effects of population size and metapopulation dynamics on a mating-system polymorphism

The evolutionary dynamics of neutral alleles under the Wright-Fisher model are well understood. Similarly, the effect of population turnover on neutral genetic diversity in a metapopulation has attracted recent attention in theoretical studies. Here we present the results of computer simulations of a simple model that considers the effects of finite population size and metapopulation dynamics on a mating-system polymorphism involving selfing and outcrossing morphs. The details of the model are based on empirical data from dimorphic populations of the annual plant Eichhornia paniculata, but the results are also of relevance to species with density-dependent selfing rates in general. In our model, the prior selfing rate is determined by two alleles segregating at a single diploid locus. After prior selfing occurs, some remaining ovules are selfed through competing self-fertilisation in finite populations as a result of random mating among gametes. Fitness differences between the mating-system morphs were determined by inbreeding depression and pollen discounting in a context-dependent manner. Simulation results showed evidence of frequency dependence in the action of pollen discounting and inbreeding depression in finite populations. In particular, as a result of selfing in outcrossers through random mating among gametes, selfers experienced a "fixation bias" through drift, even when the mating-system locus was selectively neutral. In a metapopulation, high colony turnover generally favoured the fixation of the outcrossing morph, because inbreeding depression reduced opportunities for colony establishment by selfers through seed dispersal. Our results thus demonstrate that population size and metapopulation processes can lead to evolutionary dynamics involving pollen and seed dispersal that are not predicted for large populations with stable demography.     

VARIATION IN THE MATING SYSTEM OF EICHHORNIA PANICULATA (SPRENG.) SOLMS. (PONTEDERIACEAE)

A multilocus procedure was used to estimate outcrossing rates from allozyme data in nine populations of Eichhornia paniculata from NE Brazil and Jamaica. The populations were chosen to represent stages in a proposed model of the evolutionary breakdown of tristyly to semi-homostyly; they differed in morph structure (trimorphic, dimorphic, or monomorphic) and floral traits likely to influence the mating system. The interpopulation range in outcrossing rate, t, was 0.96–0.29. Two additional populations from Jamaica, composed exclusively of self-pollinating, semi-homostylous, mid-styled plants, were invariant at 21 isozyme loci, precluding estimates of outcrossing frequency. Trimorphic populations from Brazil had uniformly high outcrossing rates of 0.96–0.88. Values for the floral morphs within populations were not significantly different. A controlled pollination experiment, comparing the competitive ability of self and cross pollen using the Got-3 marker locus, provided evidence that the maintenance of high outcrossing rates in trimorphic populations results from the prepotency of cross pollen and/or the selective abortion of selfed zygotes. Morph-dependent variation in t was detected within a dimorphic population with the L morph outcrossing with a frequency of 0.76 in comparison with 0.36 in the M morph. The difference in the mating system of floral morphs results from modifications in position of short-level stamens in flowers of the M morph resulting in automatic self-pollination. The occurrence of E. paniculata populations composed exclusively of self-pollinating, mid-styled variants is thought to be associated with the spread of genes modifying stamen position. The high level of self-fertilization demonstrated in the M morph would allow automatic selection of these genes, augmented by fertility assurance in the absence of specialized pollinators.     

Evolutionary pathways to self-fertilization in a tristylous plant species

Evolutionary transitions from outcrossing to selfing occur commonly in heterostylous genera. The morphological polymorphisms that characterize heterostyly provide opportunities for different pathways for selfing to evolve. Here, we investigate the origins and pathways by which selfing has evolved in tristylous Eichhornia paniculata by providing new evidence based on morphology, DNA sequences and genetic analysis. The primary pathway from outcrossing to selfing involves the stochastic loss of the short-styled morph (S-morph) from trimorphic populations, followed by the spread of selfing variants of the mid-styled morph (M-morph). However, the discovery of selfing variants of the long-styled morph (L-morph) in Central America indicates a secondary pathway and distinct origin for selfing. Comparisons of multi-locus nucleotide sequences from 27 populations sampled from throughout the geographical range suggest multiple transitions to selfing. Genetic analysis of selfing variants of the L- and M-morphs demonstrates recessive control of the loss of herkogamy, although the number of factors appears to differ between the forms. Early stages in the establishment of selfing involve developmental instability in the formation of flowers capable of autonomous self-pollination. The relatively simple genetic control of herkogamy reduction and frequent colonizing episodes may often create demographic conditions favouring transitions to selfing in E. paniculata .     

Correlated evolution of floral morphology and mating-type frequencies in a sexually polymorphic plant

In sexually polymorphic species, reproductive morphology governs mating patterns and the character of negative frequency-dependent selection. If local environmental conditions cause sexual morphs to differ between populations, then frequency-dependent selection should create corresponding geographic variation in morph frequencies. We investigate this relation with a model of morph-ratio evolution and analysis of geographic variation in the heterostylous plant Narcissus triandrus. Unlike other tristylous species, N. triandrus possesses both imperfect reciprocity among morphs in sex-organ position and a self-incompatibility system that permits outcrossing within and between morphs. We sampled 137 populations throughout the Iberian Peninsula for floral-morph ratios, and measured floral morphology in 31 populations. Morph ratios exhibited three atypical features: (1) predominance of the long-styled (L) morph; (2) absence of the mid-styled (M) morph from 17.5% of populations; and (3) a negative relation between the frequencies of the L and M morphs among populations. Morph ratios varied geographically, with decreasing frequency of the M morph from the southeast to the northwest of the species' range. Much of this variation accompanied allometric change in the positions of sex organs, especially the mid-level organs, with the M morph declining in frequency and ultimately being lost in large-flowered populations. Using multivariate multiple regression, we demonstrate that variation in floral morphology among populations predicts this geographic variation in morph frequencies. Our theoretical analysis illustrates that patterns of pollen transfer governed by imperfect sex-organ reciprocity can select for unequal equilibrium morph ratios like those observed for N. triandrus. We interpret the L-biased morph ratios and the unusual morphology of N. triandrus as a consequence of its atypical intramorph compatibility system.     

Characterization of 24 microsatellite markers in Primula chungensis(Primulaceae),a distylous-homostylous species,using MiSeq sequencing

Primula chungensis is a species with considerable floral and mating-system variation,including distylous(outcrossing),homostylous(selfing) and mixed populations that contain both outcrossing and selfing forms.We isolated 24 microsatellite markers from P.chungensis using Illumina Mi Seq sequencing.Polymorphism and genetic diversity were then measured based on a sample of 24 individuals from a natural population in southern Tibet.All loci were polymorphic with the number of alleles per locus ranging from 2 to 4.The observed and expected heterozygosity ranged from 0 to 1 and 0.219 to 0.708,respectively.The microsatellite markers we have identified will serve as valuable tools for the investigation of the population genetic structure and phylogeography of P.chungensis and will inform models of the evolutionary history of mating systems in the species.     

Evolutionary maintenance of stigma-height dimorphism in Narcissus papyraceus (Amaryllidaceae)

Stigma-height dimorphism is a sexual polymorphism in which plant populations are composed of two floral morphs that differ significantly in style length but not anther position. The morphs exhibit approach and reverse herkogamy, floral designs that in most species typically occur as monomorphic conditions. We investigated the floral biology of stigma-height dimorphism in the Mediterranean geophyte Narcissus papyraceus (Amaryllidaceae) in an effort to understand the evolutionary forces maintaining stylar polymorphism. Our survey of 66 populations in Spain, Portugal, and Morocco indicated that 56% were dimorphic with the long-styled morph at an average frequency of 0.79. The remaining 44% of populations sampled were monomorphic for the long-styled morph. In dimorphic populations there was a significant positive relation between population size and the frequency of the short-styled morph. Controlled pollinations demonstrated that N. papyraceus is self-sterile with no significant differences in female fertility between intra- and intermorph crosses. Prior self-pollination reduced seed set in flowers that were subsequently cross-pollinated. Estimates of mating patterns using allozyme markers in eight populations indicated that N. papyraceus is largely outcrossing (mean t(m) = 0.81) with no significant differences between monomorphic and dimorphic populations or style morphs. Stigma-height dimorphism in N. papyraceus is maintained in populations by insect-mediated cross-pollination with biased morph ratios and stylar monomorphism likely resulting from the combined influence of the inheritance of the polymorphism, morph-specific differences in assortative mating and founder effects.     

Characterization of 24 microsatellite markers in Primula chungensis (Primulaceae), a distylous homostylous species,using MiSeq sequencing

Primula chungensis is a species with considerable floral and mating-system variation, including distylous (outcrossing), homostylous (selfing) and mixed populations that contain both outcrossing and selfing forms. We isolated 24 microsatellite markers from P.chungensis using Illumina MiSeq sequencing. Polymorphism and genetic diversity were then measured based on a sample of 24 individuals from a natural population in southern Tibet. All loci were polymorphic with the number of alleles per locus ranging from 2 to4. The observed and expected heterozygosity ranged from 0 to 1 and 0.219 to 0.708, respectively. The microsatellite markers we have identified will serve as valuable tools for the investigation of the population genetic structure and phylogeography of P.chungensis and will inform models of the evolutionary history of mating systems in the species.     

Female reproductive success and the evolution of mating-type frequencies in tristylous populations

In tristylous populations, mating-type frequencies are governed by negative frequency-dependent selection typically resulting in equal morph ratios at equilibrium. However, Narcissus triandrus generally exhibits long-styled (L)-biased populations with a deficiency of the mid-styled (M)-morph. Here we used a pollen-transfer model and measurements of female fertility in natural populations to investigate whether these uneven morph ratios were associated with variation in female reproductive success. Our theoretical analysis demonstrated that morph ratio bias can result from maternal fitness differences among the morphs, and that these effects were magnified by asymmetrical mating. In nine out of 15 populations of N. triandrus, seed set differed significantly among the morphs, but pollen limitation occurred in only two of 11 populations investigated. Average seed set of the M-morph was positively associated with its frequency in populations. Flower size was negatively correlated with the seed set of the M-morph. Our results suggest that interactions between mating patterns and female fertility are responsible for variation in morph frequencies and loss of the M-morph from tristylous populations of N. triandrus.     

Patterns of Floral Variation between Dimorphic and Monomorphic Populations in Distylous Luculia pinceana (Rubiaceae)

Flower morphology plays an important role in the evolution and maintenance of plant mating systems, including disassortative mating of heterostylous species. The transition of mating patterns may be associated with the remodification of intraspecific flower morphology. To determine the functional relationship between floral variation and transition of mating patterns, we conducted a series of morphometric analyses in a distylous species Luculia pinceana, which possesses dimorphic and monomorphic populations. Our results indicate that floral variation was higher between different types of populations than between populations of the same type. Compared to dimorphic populations, some floral characters, reduced stigma anther separation within flowers and increased overlap of stigmas and anthers (illegitimate spatial matching of sexual organs) among individuals in populations containing only the long styled morph may have been modified to increase both selfing and intra morph crossing. The observed patterns of floral variation between dimorphic and monomorphic populations coincide with the transition of mating patterns from disassortative mating to selfing and/or intra morph crossing.     

EFFECTIVE POPULATION SIZE AND GENETIC DRIFT IN TRISTYLOUS EICHHORNIA PANICULATA (PONTEDERIACEAE)

Populations of the tristylous, annual Eichhornia paniculata are markedly differentiated with respect to frequency of mating types. This variation is associated with evolutionary changes in mating system, from predominant outcrossing to high self-fertilization. To assess the potential influence of genetic drift acting on this variation, we estimated effective population size in 10 populations from northeastern Brazil using genetic and demographic methods. Effective size (N_e) was inferred from temporal changes in allele frequency at two to eight isozyme loci and also calculated using five demographic variables: 1) the number of flowering individuals (N); 2) temporal fluctuations in N; 3) variance in flower number; 4) frequency of mating types; and 5) selfing rate. Average N_e based on isozyme data was 15.8, range 3.4–70.6, and represented a fraction (mean N_e/N = 0.106) of the census number of individuals (mean N = 762.8; range: 30.5–5,040). Temporal variation in N and variance in flower number each reduced N_e to about a half of N whereas mating type frequencies and selfing rate caused only small reductions in N_e relative to N. All estimates of N_e based on demographic variables were considerably larger than those obtained from genetic data. The two kinds of estimates were in general agreement, however, when all demographic variables were combined into a single measure. Monte Carlo simulations indicated that effective size must be fewer than about 40 for drift to overcome the frequency-dependent selection that maintains the polymorphism for mating type. Applying the average N_e/N value to 167 populations censused in northeastern Brazil indicated that 72% had effective sizes below this number. This suggests that genetic drift is likely to play a dominant role in natural populations of E. paniculata.     

The relative importance of factors determining genetic drift: mating system, spatial genetic structure, habitat and census size in Arabidopsis lyrata

? The mating system, dispersal and census size are predicted to determine the magnitude of genetic drift, but little is known about their relative importance in nature. ? We estimated the contributions of several population-level features to genetic drift in 18 populations of Arabidopsis lyrata. The factors were outcrossing rate, within-population spatial genetic structure, census size and substrate type. The expected heterozygosity (H(E)) at 10 microsatellite loci was taken to reflect the effective population size (N(e)) and the strength of genetic drift. ? The mating system explained most of the variation in H(E) (60%), followed by substrate (10%), genetic structure (9%) and census size (6%). The most outcrossing population had a +0.32 higher predicted H(E) than the most selfing population; the estimated N(e) of selfing populations was less than half that of outcrossing populations. Rocky outcrops supported populations with a +0.14 higher H(E) than did sandy substrates. The most structured population had a +0.24 higher H(E) than the least structured population, and the largest population had a +0.18 higher H(E) than the smallest population. ? This study illustrates the importance of outcrossing, genetic structure and the physical environment--together with census size--in maintaining H(E), and suggests that multiple population-level characteristics influence N(e) and the action of genetic drift.     

Floral variation and environmental heterogeneity in a tristylous clonal aquatic of the Pantanal wetlands of Brazil

MethodsStyle morph frequencies were sampled from 73 populations distributed across four flooding regimes differing in depth and duration. Measurements of flower size, sex-organ dimension, pollen size and pollen production were made in selected populations, and pollinator assemblages and their functional traits were recorded.ConclusionsThe similar morph structure and evenness of populations, regardless of flooding regime, suggest that sexual reproduction and clonal dispersal are sufficiently common to prevent the signature of founder events from dominating in a region. However, the pervasive occurrence of biased morph ratios in most populations suggests that many are in a non-equilibrium state. The reduced frequency of the mid-styled morph in trimorphic and dimorphic populations may be associated with the weak self-incompatibility of this morph resulting in selfing and inbreeding depression. Clonality in E. azurea and the weak self-incompatibility of the mid-styled morph may make it more vulnerable to geitonogamous selfing.     

Pollinator visitation in populations of tristylous Eichhornia paniculata in northeastern Brazil

Summary The frequencies of floral morphs in populations of tristylous Eichhornia paniculata often deviate from the theoretical expectation of equality. This variation is associated with the breakdown of tristyly and the evolution of self-fertilization. Differences in morph frequencies could result from selection pressures due to variable levels of insect visitation to populations and contrasting foraging behavior among the floral morphs. We estimated pollinator densities in 16 populations and quantified visitation sequences to morphs in five populations of E. paniculata in northeastern Brazil. Foraging behavior among floral morphs was measured as the frequency of visits to morphs relative to their frequency in the population (preference) and number of flights between inflorescences of the same versus different morphs (constancy). Pollinator density (number/m²/minute) was not correlated with population size, plant density or morph diversity. Pollinator densities varied most among populations of less than 200 plants. Whether pollinators discriminated among the morphs, depended on whether they primarily collected nectar or pollen. In four populations, nectar-feeding bees (Ancyloscelis and Florilegus spp.) and butterflies showed no consistent preference or constancy among the morphs. In contrast, pollen-collecting bees (Trigona sp.) visited a lower proportion of longstyled inflorescences than expected and tended to visit more mid-and short-styled inflorescences in succession, once they were encountered. Pollinator constancy for morphs did not result from differences in inflorescence production or spatial patchiness among the morphs. Although non-random pollinator visitation to morphs in heterostylous populations could potentially affect mating and hence morph frequencies, the observed visitation patterns in this study do not provide evidence that pollinators play a major role in influencing floral morph frequencies.     

Inferring recent outcrossing rates using multilocus individual heterozygosity: application to evolving wheat populations

Using multilocus individual heterozygosity, a method is developed to estimate the outcrossing rates of a population over a few previous generations. Considering that individuals originate either from outcrossing or from n successive selfing generations from an outbred ancestor, a maximum-likelihood (ML) estimator is described that gives estimates of past outcrossing rates in terms of proportions of individuals with different n values. Heterozygosities at several unlinked codominant loci are used to assign n values to each individual. This method also allows a test of whether populations are in inbreeding equilibrium. The estimator's reliability was checked using simulations for different mating histories. We show that this ML estimator can provide estimates of outcrossing rates for the final generation outcrossing rate (t(0)) and a mean of the preceding rates (t(p)) and can detect major temporal variation in the mating system. The method is most efficient for low to intermediate outcrossing levels. Applied to nine populations of wheat, this method gave estimates of t(0) and t(p). These estimates confirmed the absence of outcrossing t(0) = 0 in the two populations subjected to manual selfing. For free-mating wheat populations, it detected lower final generation outcrossing rates t(0) = 0-0.06 than those expected from global heterozygosity t = 0.02-0.09. This estimator appears to be a new and efficient way to describe the multilocus heterozygosity of a population, complementary to Fis and progeny analysis approaches.     

Geographic variation in floral morphology and style-morph ratios in a sexually polymorphic daffodil

Narcissus triandrus is a tristylous daffodil from the Iberian Peninsula that has striking geographical variation in floral morphology and style-morph ratios. Here, we investigate the relation between this variation and ecological factors to understand mechanisms governing morph ratios. We estimated morph ratios in 124 populations throughout the range of N. triandrus and measured 13 morphological traits in 35-78 populations. Sampling of morph ratios conducted over 2-10-yr intervals demonstrated strong temporal stability. Variation in floral and vegetative traits enabled statistical prediction of morph-frequency variation among populations. Latitudinal gradients in precipitation and temperature were correlated with plant and flower size, with larger flowers in northern populations associated with bumblebee visitation and stylar dimorphism. Flowers of the L- and M-morphs differed significantly from the S-morph in several size-related characters, unlike other tristylous species. This pattern and the similarity in anther positions of the L- and M-morphs suggest that the M-morph of N. triandrus originated through genetic modifiers that shortened styles of the L-morph. Our findings support the hypothesis that geographical variation in style-morph ratios is largely governed by climatic gradients in the Iberian Peninsula, which influence the floral morphology and pollination biology of N. triandrus populations.     

Mating patterns and demography in the tristylous daffodil Narcissus triandrus

Mating patterns in plant populations are influenced by interactions between reproductive traits and ecological conditions, both factors that are likely to vary geographically. Narcissus triandrus, a wide-ranging heterostylous herb, exhibits populations with either two (dimorphic) or three (trimorphic) style morphs and displays substantial geographical variation in demographic attributes and floral morphology. Here, we investigate this variation to determine if demography, morphology, and mating system differ between the two sexual systems. Our surveys in Portugal and NW Spain indicated that dimorphic populations were less dense, of smaller size, and had larger plants and flowers compared to trimorphic populations. Outcrossing rates estimated using allozyme markers revealed similar outcrossing rates in dimorphic and trimorphic populations (t(m) dimorphic=0.759; t(m) trimorphic=0.710). All populations experienced significant inbreeding in progeny (mean F=0.143). In contrast, parental estimates of inbreeding were not significantly different from zero (mean F=0.062), implying that few inbred offspring survive to reproductive maturity due to inbreeding depression. Although the majority of inbreeding results from selfing, significant levels of biparental inbreeding were also detected in eight of the nine populations (mean s(s)-s(m)=0.081). Density was negatively associated with levels of selfing but positively associated with biparental inbreeding. Population size was positively associated with outcrossing but not biparental inbreeding. There were no consistent differences among the style morphs in outcrossing or biparental inbreeding indicating that the maintenance of trimorphism vs dimorphism is unlikely to be associated with inbreeding of maternal parents.     

THE DISTRIBUTION OF SELFING RATES IN HOMOSPOROUS FERNS

The models of Lande and Schemske predict that among species in which the selfing rate is largely under genetic control and not subject to tremendous environmental variation, the distribution of selfing rates should be bimodal. When this prediction was tested empirically using data from the literature for species of angiosperms and gymnosperms, the distribution of outcrossing rates for all species was clearly bimodal. To provide another empirical test of the prediction, we analyzed mating-system data for 20 species of Pteridophyta (ferns). Homosporous ferns and their allies are unique among vascular plants because three types of mating are possible: intragametophytic selfing (selfing of an individual gametophyte); intergametophytic selfing (analogous to selfing in seed plants); and intergametophytic crossing (analogous to outcrossing in seed plants). The distribution of intragametophytic selfing rates among species of homosporous ferns is clearly uneven. Most species of homosporous ferns would be classified as extreme outcrossers. In contrast, a few species are nearly exclusively inbreeding. In only a few populations of Dryopteris expansa and Hemionitis palmata and a single population of Blechnum spicant do we see convincing evidence of a mixed mating system. The uneven distribution of selfing rates we observed for homosporous ferns, coupled with a corresponding bimodality of the magnitude of genetic load, strongly supports the model.