Am. J. Bot. Cross-Journal Searching
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a response
Right arrow Alert me when this article is cited
Right arrow Alert me when eLetters are posted
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Similar articles in PubMed
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (22)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Griffin, S. R.
Right arrow Articles by Eckert, C. G.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Griffin, S. R.
Right arrow Articles by Eckert, C. G.
Agricola
Right arrow Articles by Griffin, S. R.
Right arrow Articles by Eckert, C. G.
(American Journal of Botany. 2000;87:1246-1256.)
© 2000 Botanical Society of America, Inc.

Experimental analysis of protogyny in Aquilegia canadensis (Ranunculaceae)1

Steven R. Griffin2,0, Katherine Mavraganis0 and Christopher G. Eckert3,0

0 Department of Biology, Queen's University, Kingston, Ontario, Canada K7L 3N6

Received for publication July 27, 1999. Accepted for publication October 29, 1999.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Dichogamy is very common in flowering plants and is widely thought to reduce pollen-pistil interference, especially self-pollination. Yet, the functional significance of dichogamy has rarely been studied. We investigated the nature and functioning of dichogamy in eastern Ontario populations of Aquilegia canadensis, a highly selfing columbine previously described as protogynous. We then manipulated flowers to determine whether increased protogyny would reduce self-fertilization. Contrary to previous reports, A. canadensis is not dichogamous. Controlled pollinations in a greenhouse showed that pollen tubes generally begin to develop after anther dehiscence. Although stigmas can collect pollen early in floral development, naturally pollinated flowers collected from four populations had few pollen grains on stigmas and almost no pollen tubes in styles until after anther dehiscence. Limited pollen deposition before anther dehiscence was also associated with low nectar availability and limited sepal expansion. Because inbreeding depression is strong in this species, selection may favor increased protogyny if it reduces selfing. We tested this hypothesis by comparing the level of selfing in flowers rendered protogynous by the removal of the first 19 (of 39) anthers to develop, with nonprotogynous control flowers. Contrary to expectations, protogyny did not reduce selfing. Our results emphasize the importance of detailed field observations and manipulative experiments in understanding the nature and functional significance of dichogamy.

Key Words: Aquilegia canadensis • columbine • dichogamy • floral morphology • nectar • plant mating systems • pollen tubes • pollination • protogyny • self-fertilization


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Because plants are sessile, they rely on a variety of biotic and abiotic pollen vectors (animals, wind, and water) to mediate gamete transfer and use an array of physiological and morphological mechanisms to regulate mating. Temporal separation of male and female function within flowers (intrafloral dichogamy) is one of the most widespread morphological mechanisms and is found in >75% of cosexual angiosperm species (Bertin and Newman, 1993 ). There are two forms of dichogamy: protandry, where anthers dehisce before stigmas become receptive, and protogyny, where receptive stigmas are presented before anther dehiscence (Lloyd and Webb, 1986 ; Bertin and Newman, 1993 ). There is wide variation in the degree to which the two sexual phases are separated among species and sometimes within species (Bertin and Newman, 1993 ).

In general, dichogamy is thought to have evolved to reduce interference between pollen import and export (Wyatt, 1983 ; Lloyd and Webb, 1986 ; Bertin and Newman, 1993 ). A particularly important form of interference is self-pollination. In species that lack physiological self-incompatibility, the contamination of stigmas with self-pollen can lead to self-fertilization, which may reduce fitness if selfed offspring suffer inbreeding depression. Dichogamy is widely expected to limit self-fertilization in self-compatible taxa (Proctor and Yeo, 1972 ; Faegri and van der Pijl, 1979 ; Richards, 1986 ; Cruden, 1988 ; Barrett and Eckert, 1990 ). Protogyny, in particular, is viewed as an anti-selfing mechanism because it provides opportunities for the receipt of outcross pollen before self-pollen is shed, and it is more common in self-compatible than self-incompatible taxa (Lloyd and Webb, 1986 ; Bertin, 1993 ). However, the effect of either form of dichogamy on the mating system has rarely been investigated. A few correlative studies have detected positive covariation between the level of self-fertilization and degree of dichogamy among populations (Schoen, 1982 ; Holtsford and Ellstrand, 1992 ) or between individuals within populations (Brunet and Eckert, 1998 ). Yet, there has been little functional analysis of dichogamy in natural populations (e.g., Jonsson, Rosquist, and Widén, 1991 ; Preston, 1991 ; Dinnétz, 1997 ) and no experimental manipulation of the degree of dichogamy to determine the impact on self-fertilization.

The extent to which dichogamy reduces self-fertilization depends on (1) the timing of stigma receptivity in relation to anther dehiscence, and (2) the rate at which outcross pollen is deposited on stigmas during female phase and removed from anthers during male phase (Preston, 1991 ). In protandrous species, autogamous selfing may occur if viable pollen remains in anthers when receptive stigmas are presented. For example, Brunet and Eckert (1998) detected substantial (50%) selfing in strongly protandrous Aquilegia caerulea despite an average of 2.9 d separating the dehiscence of the last anther and the beginning of stigma receptivity. Likewise, in protogynous species, autogamous selfing may occur if stigmas are receptive and unfertilized ovules are still available when anthers shed self pollen. Further opportunities for self-fertilization arise in protandrous species if stigmas are presented and can receive self pollen before they become receptive. Consequently, the functioning of dichogamy depends on the rates of pollen deposition and removal, which, in turn, depend on pollinator visitation rates and foraging behavior.

In this study, we investigated the functional significance of protogyny in reducing self-fertilization in Aquilegia canadensis L. (Ranunculaceae), a spring flowering, short-lived perennial herb native to eastern central North America. This species usually occurs in small populations (<200 flowering plants) located on rocky outcrops and in dry woods (Payson, 1918 ; Munz, 1946 ). Individual plants bear one to several large, pendant, red flowers with deep nectar spurs that hold dilute nectar. Although floral morphology and coloration suggest adaptation to hummingbird pollination (Grant and Grant, 1968 ; Macior, 1978 ; Faegri and van der Pijl, 1979 ; Grant and Temeles, 1992 ), solitary bees and especially bumble bees may also be important pollinators (Schneck, 1901 ; Macior, 1966 ).

The flowers of A. canadensis were first described as protogynous by Schneck (1901) who observed that the styles "protrude and are ready to receive pollen from other flowers before the anthers in the same flower are ready to shed pollen" and that "close-pollination is avoided by the protogynous condition of the flower." The occurrence of protogyny in these large, rewarding flowers suggest that this species practices high levels of outcrossing. However, there are 3 observations that suggest that protogyny is either ineffective at reducing selfing and/or that there is extensive overlap between stigma receptivity and anther dehiscence. First, like other species of Aquilegia for which data are available, the flowers of A. canadensis are highly autofertile; they set almost a full complement of seeds via automatic self-pollination when isolated from pollinators (Macior, 1978 ; Eckert and Schaefer, 1998 ; Routley, Mavraganis, and Eckert, 1999 ; for other species see Miller, 1978, 1985 ; Miller and Willard, 1983 ; Brunet and Eckert, 1998 ). This suggests that protogyny is, to some extent, incomplete; stigmas are still receptive when self-pollen is shed. Second, Eckert and Schaefer (1998) observed very low levels of pollen deposition before anther dehiscence in one population, suggesting that even if flowers are protogynous, the opportunity for cross-pollination before anther dehiscence provided by protogyny is not realized. Finally, allozyme analysis of progeny arrays sampled from ten natural populations from eastern Ontario, Canada, indicated very high levels of self-fertilization (average = 75%; range = 17–100%; Routley, Mavraganis, and Eckert, 1999 ). Comparing the levels of selfing in these populations to the inbreeding coefficients of mature plants further suggested that inbreeding depression is very strong. Outcrossed offspring are more than ten times as fit as selfed offspring (Routley, Mavraganis, and Eckert, 1999 ).

Here, we investigate the nature of protogyny in A. canadensis and determine the extent to which it is, or could be, effective at reducing self-fertilization. We begin by describing floral development and quantifying variation in pollen receipt and pollen tube growth across floral development in both the greenhouse and field to assess the extent to which cross-fertilization may precede opportunities for self-fertilization. We then quantify temporal variation in floral attractiveness to see whether there are floral cues associated with the onset of stigma receptivity and anther dehiscence that may affect the timing of cross-pollination relative to self-pollination. Our results suggest that contrary to previous observations, flowers of A. canadensis are not protogynous; there is almost complete overlap between stigma receptivity and anther dehiscence, and self- and cross-pollination appear to occur simultaneously.

High levels of self-fertilization combined with very strong inbreeding depression suggest that protogyny would be advantageous in populations of A. canadensis if it promoted cross-fertilization. That protogyny does not occur in eastern Ontario populations of A. canadensis is particularly puzzling because the genus Aquilegia contains closely related species that appear to be strongly protogynous (Chase and Raven, 1975 ; Miller, 1978 ) as well as species that are strongly protandrous (Miller, 1978, 1985 ; Miller and Willard, 1983 ; Brunet and Eckert, 1998 ). Moreover, the species in this genus seem to have differentiated rapidly and very recently from each other (Hodges and Arnold, 1994 ), suggesting that the nature and degree of dichogamy are evolutionarily labile and can be rapidly altered by selection. Because the pattern of floral development in A. canadensis enables the creation of artificially protogynous flowers, we can investigate whether protogyny might be an effective mechanism for reducing self-fertilization in A. canadensis. We used floral manipulations in conjunction with marker-gene analysis to compare the level of self-fertilization in experimentally protogynous flowers vs. nondichogamous control flowers.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Flower development
To determine the temporal separation of stigma presentation and anther dehiscence, we monitored flower development in May and June of 1995 in a large population from eastern Ontario, Canada (CCR1; population descriptions and locations in Mavraganis, 1998 ). We randomly chose 39 plants, each with at least three unopened buds, and excluded them from pollinators with a wire mesh cylinder covered with fine bridal veil. Flowers were either cross- or self-pollinated when anthers began shedding pollen (for details, see Routley, Mavraganis, and Eckert, 1999 ). Every day throughout flower and fruit development we recorded flower color, the degree to which stigmas were exserted, and the position and dehiscence of anthers. Based on these observations, we divided flower and fruit development into 13 morphologically distinct stages that will be described in detail in the results (Table 1).


View this table:
[in this window]
[in a new window]
 
Table 1. Flower and fruit development in Aquilegia canadensis. Spur color changed as flowers developed from yellow-green (YG) through green (G), green-red (GR), and red (R). Means ± 1 SD for duration were calculated with the individual flower as the unit of observation. Morphological definitions of stages 10–13 are given in the text.

 
Variation in pollen receipt and pollen tube growth across floral development in hand-pollinated flowers
We examined variation in stigma receptivity and pollen tube growth across floral development in a greenhouse experiment involving two populations from eastern Ontario. Before the onset of flowering in early May 1997, 19 and 20 plants were randomly collected from populations CCR1 and QSP1, respectively, and grown in a randomized array in a common greenhouse environment. Each plant was housed in a 7.5-cm pot filled with soilless potting medium (Sunshine MixTM Number 1) and grown to flowering under ambient light and temperatures of 15°–25°C.

Each plant produced at least four flowers, and one flower on each plant was hand-pollinated at stages 3, 5, 6, or 8 (Table 1). For 3 of 39 plants, a flower was not pollinated at stage 3. Hand-pollinations involved cross pollen from a single freshly dehisced anther taken from another plant in the same population. To prevent contamination of stigmas with self-pollen, undehisced anthers were removed from flowers each day as they moved towards the stigmas on unfurling filaments. For each flower, a stigma from one of the n carpels (n = 3–5, median = 5) was collected before pollination to examine variation in the morphology of the stigma surface across floral development. Unfortunately, most of these samples were accidentally destroyed and only fragmentary information could be obtained (see Discussion). The remaining n - 1 carpels were pollinated and collected 24 h later. The pollen grains on the stigmas of one or two of these carpels (median = 2) and the pollen tubes in all n - 1 styles were visualized and counted as follows.

Stigmas were squashed in fuschin jelly under a glass coverslip, and the pollen grains on and close to each stigma were counted (Kearns and Inouye, 1993 , pp. 115–119). If more than one stigma was assayed, the number of pollen grains per stigma (Np) was averaged to yield one value per flower. Pollen tubes were visualized using fluorescence microscopy following Kearns and Inouye (1993 , pp. 124–129). The carpels from each flower were fixed in FAA (formalin-acetic acid-alcohol) for 24 h and then stored in 70% ethanol. Individual carpels were then separated using forceps, softened in 5 mol/L NaOH overnight (from 8 to 24 h), rinsed three times with tap water, and stained in 0.1% (w/v) aniline blue solution decolorized overnight in 0.1 mol/L K3PO4. Each carpel was mounted on a microscope slide with 10% glycerol and squashed under a glass cover slip. Pollen tubes were viewed under a Leica DMRBTM stereo microscope at 200x using 350-nm excitation and 420-nm barrier filters.

Each style was visually divided into k 1-mm sections (k {approx} 14)and the number of pollen tubes in each section (xi) was recorded (i = 1 - k). The mean length of pollen tubes (L) in each style was calculated as:

where di is the distance from the beginning of the ith segment to the stigma (d1 = 0 mm). We also recorded the maximum number of pollen tubes in any segment of the style (xmax, usually xmax = x1). These two variables (L and xmax) were averaged across carpels (usually three) to yield one value for each variable per flower.

We assessed variation in these measures of pollen receipt (Np) and pollen tube growth (L and xmax) among floral developmental stages and between populations using repeated-measures ANOVA with plant as subject, population as a random, between-subject effect and developmental stage as a fixed, within-subject effect (Neter, Wasserman, and Kutner, 1990 , pp. 1057–1066). The F test for population used the MS for plant[population] as denominator. All other F tests used the residual MS as denominator.

Number of pollen grains/stigma (Np) was square root transformed and L and xmax were log10 transformed to meet the assumptions of ANOVA. However, even after transformation the residuals from the analysis of Np were not entirely independent of the predicted values (Pearson r = +0.18, N = 153, P = 0.025). In this case, analysis of the raw data, transformed data, and even ranked data yielded the same qualitative results, suggesting that the analysis presented below is probably robust to violation of ANOVA assumptions. These analyses and others described below were performed using JMP (version 3.2.1; SAS, 1997 ).

Variation in pollen deposition and pollen tube growth across floral development in naturally pollinated flowers
To determine the schedule of pollen deposition and pollen tube growth across floral development in naturally pollinated flowers, we compared stigmatic pollen loads and pollen tube growth across seven stages of floral development (stages 3–9, see Table 1) in four eastern Ontario populations (AHB4, AKA1, CCP3, and CCR4) during May 1995. Ten flowers per stage (one per plant) were haphazardly sampled on a single day in each population. We counted the number of pollen grains on each stigma of each flower (as above) and calculated the average number of pollen grains per stigma for each flower. We assayed pollen tube numbers in all styles of each flower (as above) and recorded the maximum number of pollen tubes in each style as well as the number of tubes at the base of each style.

We assessed variation in pollen grain and pollen tube numbers among developmental stages in each population using one-way ANOVA with Tukey-Kramer multiple contrasts. We assessed overall variation among populations and developmental stages and whether differences between stages varied among populations using two-way mixed-model ANOVA with floral stage as a fixed effect and population as a random effect. F tests for both population and floral stage used the interaction MS as denominator. Because early floral stages had very few pollen grains on stigmas and almost no pollen tubes growing in the styles, the data were not normally distributed, there was strong heteroscedasticity in sample variances, and the data could not be successfully transformed to meet assumptions of ANOVA. Log10 transformation best improved normality of residuals and somewhat reduced the correlation between residuals and predicted values. Transformation did not qualitatively change the effect of developmental stage, which was very strong and highly significant for all variables. However, log10 transformation did affect the significance of the interaction between population and developmental stage. Hence, analyses based on log10 transformed data are reported below and should be interpreted with caution.

Variation in nectar availability and flower size across floral development in naturally pollinated flowers
To assess the potential attractiveness of flowers to pollinators in relation to periods of stigma receptivity and anther dehiscence, we measured nectar availability and components of size for flowers in stages 4, 5, 6, 7, and 8 (see Table 1) in two eastern Ontario populations (CCR1 and QBS1) in May 1997. At least ten flowers per stage (one per plant) were haphazardly sampled in each population. Nectar availability was measured as the volume of nectar in the spur tips (standing crop) by inserting a 1 x 50-mm strip of chromatography paper (Whatman 3 MM) into each of two randomly chosen nectar spurs. The strips were left in the spurs for 1 min to soak up all the nectar. Paper strips were then removed and the length of nectar migration on the chromatography paper was measured to 0.1 mm using Vernier calipers. Migration lengths of the two spurs for each flower were strongly correlated (Pearson r = +0.83, N = 124, P < 0.00001 for QBS1; and r = +0.91, N = 100, P < 0.00001 for CCR1), hence they were averaged to yield one value per flower.

A standard curve of migration length over known volume of artificial nectar with the same sucrose concentration as nectar taken from flowers in natural populations (25% sucrose; see Macior, 1978 ) was calculated and used to convert the average migration length for each flower into a nectar volume. The standard curve used the migration lengths of ten replicate strips for each of five known volumes (2, 4, 6, 8, and 10 µL). A regression of length (in millimetres) over volume (in microlitres) forced through the origin yielded the following conversion: volume = 0.262 x length (r2 = 0.99).

We measured three components of flower size to 0.1 mm using vernier calipers: (1) the maximum length of one randomly chosen spur; (2) the length of one randomly chosen sepal; and (3) sepal spread, measured as the distance between the tips of two randomly chosen opposite sepals.

We assessed variation in nectar availability and components of flower size among developmental stages in each population using one-way ANOVA with Tukey-Kramer multiple contrasts. We assessed overall variation in both flower size and per-flower nectar availability among populations and stages and whether differences among stages varied between populations using two-way ANOVA with developmental stage as a fixed effect and population as a random effect. F tests for both main effects used the interaction MS as denominator. Because flowers in early developmental stages (4 and 5) produced very little nectar, average nectar volumes were not normally distributed, and there was strong heteroscedasticity in sample variances. Thus, the data were inverse transformed (Y' = 1/[Y + 1]) to meet assumptions of ANOVA. Likewise, sepal length was log10 transformed to meet assumptions of ANOVA. The other two analyses used untransformed data.

Experimental manipulation of protogyny
In May and June 1998, we compared the level of self-fertilization in artificially protogynous flowers with non protogynous control flowers in three large populations (>150 flowering plants) in eastern Ontario (QBC1, QCA1, QWM1). We randomly chose plants with one or more flowers in which no anthers had dehisced (<stage 6) and used one flower per plant. A flower on each plant was assigned to one of four treatments. Protogynous flowers (treatment "P") were created by removing the 19 anthers that would normally dehisce first (out of an average total of 39 anthers per flower). Stamens are inserted in a column of whorls around the style and anthers most distal to the pedicel reflex and dehisce before proximal anthers. Thus the order in which anthers dehisce can be accurately predicted. These protogynous flowers may self-fertilize at a lower level than unmanipulated flowers for two reasons: (1) because artificial protogyny allows cross-fertilization before opportunities for self-fertilization; and (2) because they contain much less self-pollen. Accordingly, the protogynous flowers were compared with control flowers from which 19 randomly chosen anthers were removed before dehiscence (treatment "C"). Compared to C flowers, P flowers experienced about one day of stigma receptivity before anthers started shedding self-pollen. We also performed a third treatment in which the last 19 anthers to dehisce were removed before they started shedding pollen ("L"). In unmanipulated flowers, the 39 anthers dehisce over a 3–4 d period, however, stigmas usually acquire enough pollen to fertilize all ovules in the first day of anther dehiscence (see Results). Consequently, removing the last 19 anthers should have no effect on the level of self-fertilization.

To minimize the effect of other factors that might influence the level of self-fertilization (e.g., display size, position in the flowering sequence, time of flowering), we treated quartets of plants matched for inflorescence number and flower number. One flower on each plant received one of the three treatments (P, C, or L) or was left unmanipulated ("U"). The treated flowers were further matched for their time of flowering and position within the flowering sequence. Floral manipulations were performed between 5 May and 1 June, during which time at least ten plants in each population were flowering on any given day. We treated between 25 and 40 quartets for each population, but attrition due to herbivory, seed predation, and drought greatly reduced sample sizes. The number of fruits collected from P, C, L, and U flowers, respectively, was 20, 18, 13, and 15 in QBC1; 12, 15, 11, and 10 in QCA1; and 18, 9, 15, and 10 in QWM1. We collected most fruits just before the follicles split open and counted the plump, black glossy seeds in each complete fruit. There was an average (±1 SE) of 10.0 ± 0.5 seeds/follicle (N = 160), which did not vary significantly among populations or treatments (two-way ANOVA, whole-model F11 148 = 1.34, P = 0.21). We also collected fruits from as many additional unmanipulated plants ("X") to increase the sample of fruits used to obtain a baseline estimate of self-fertilization for each population (N for U and X combined = 35 in QBC1, 32 in QCA1 and 29 in QWM1).

We assayed ten seeds per fruit (total = 2283) for two polymorphic allozyme loci (PER and IDH) following Routley, Mavraganis, and Eckert (1999) . There were three alleles for each locus, although one allele for each was rare and not found in all populations. To obtain multilocus and single-locus estimates of the level of outcrossing (tm and ts, respectively) for each treatment (P, C, and L), we used Ritland's maximum-likelihood program MLTR, which estimates mating parameters for different classes of plants that experience a common set of pollen allele frequencies (Ritland, 1986, 1990 ). Newton-Raphson iteration was used to find maximum-likelihood estimates, and standard errors were calculated as the SD of 1000 bootstrap values generated with the progeny array as the unit of resampling. The multilocus estimate of the proportion of seeds produced through selfing was calculated as s = 1 - tm. We obtained an estimate of the baseline level of self-fertilization for each population in a similar fashion using only the progeny arrays from unmanipulated flowers (U and X). For all analyses, iterations from any set of starting values converged on only one set of maximum-likelihood estimates.

To determine whether experimental protogyny reduced the level of self-fertilization, we contrasted estimates of s between treatments within populations with a priori, nonorthogonal, pairwise comparisons (P vs. C, and C vs. L), and a per-comparison {alpha} of 2.5%. Because we expected that protogyny should reduce selfing, the P vs. C comparison was one-tailed (e.g., HO: sP = sC; HA: sP < sC). Because we expected that the removal of late anthers should have no effect on selfing, the C vs. L comparison was two-tailed (e.g., HO: sC = sL; HA: sC != sL). The one-tailed contrast was considered significant if <2.5% of the differences (sP - sC) between randomly paired bootstrap estimates were >0 (Eckert and Barrett, 1994 ). The two-tailed contrast was considered significant if <1.25% of the differences (sC - sL), were >0 or <0. Maximum-likelihood estimates of s and individual bootstrap estimates were averaged across populations to yield the average level of self-fertilization for the three emasculation treatments (P, C, and L). The standard error of the average was calculated as the SD of the 1000 average bootstrap values (see Eckert and Barrett, 1994 , for details). These average estimates were then contrasted with the same procedure used for the single population estimates.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Floral development
We divided flower and fruit development into 13 relatively distinct morphological stages (Table 1). Stage-1 flowers were small yellow-green buds. The development of clearly defined nectar spurs and a hint of a reddish color marked the onset of stage 2. At stage 3, nectar spurs were well defined, buds were a deeper red, and styles elongated to present stigmas outside closed buds. At stage 4, the sepals began to open and expose the lamellae (openings of the nectar spurs). The stamens were curled back inside the corolla and surrounded the carpels. At stage 5, the sepals had spread farther and some stamens had straightened out to bring anthers towards the stigmas, though none had dehisced. Collectively, stages 3–5 are called "stigma presentation."

Stage 6 was marked by the beginning of anther dehiscence. Fewer than ten stamens were fully unfurled to place dehiscing anthers close to (0–10 mm) the stigmas. Sepals were usually fully open at this stage. Stage-6 flowers were deep red with bright yellow lamellae. Stage-7 flowers were much like stage-6 flowers except that more anthers (11–29) had dehisced. At stage 8, almost all (>30) anthers had dehisced. At stage 9, all anthers had dehisced and the stigmas had started to curl. Collectively, stages 6–9 are called "anther dehiscence." Flowers spent more time in stigma presentation (mean ± 1 SE = 5.2 ± 1.9 d) than anther dehiscence (3.2 ± 0.9 d, paired t = 3.7, P = 0.0004).

At stage 10, nectar spurs, sepals and stamens abscised, thereby exposing the carpels. At stage 11, the unfused, thin, green carpels (follicles) had begun to swell. At stage 12, follicles were noticeably swollen. At stage 13, follicles were a darker green and were separating from one another. Stage 13 ended when mature follicles split to release seeds. On average, fruit development (stages 10–13) averaged 22.7 ± 3.1 d.

Variation in pollen receipt and pollen tube growth across floral development in hand-pollinated flowers
The results of the greenhouse pollination experiment suggest that although stigmas are presented well outside flowers several days before anther dehiscence, they are not very receptive until anthers start to dehisce (Fig. 1). The number of pollen grains deposited on stigmas via hand-pollination varied somewhat with developmental stage, although the pattern of variation differed among the two populations (Table 2). In CCR1, pollen receipt increased more or less monotonically from an average (±1 SE) of 60.3 ± 11.1 grains on stage-3 stigmas to 159.4 ± 23.7 grains on stage-8 stigmas (Fig. 1). In QSP1, pollen loads were highest on stage-3 stigmas (166.1 ± 23.3) with only a moderate increase from stage 5 (99.6 ± 17.3) to stage 8 (145.9 ± 19.3). For both populations, the number of pollen grains deposited on stigmas at all stages was much greater than the average number of ovules per carpel (average = 25 ovules; Routley, Mavraganis, and Eckert, 1999 ).



View larger version (20K):
[in this window]
[in a new window]
 
Fig. 1. Variation in pollen receipt and pollen tube growth among stages of floral development in hand-pollinated flowers from two populations of Aquilegia canadensis grown under greenhouse conditions. Points are means (±1 SE) based on 20 flowers per stage per population. Letters next to points show the results of multiple contrasts among developmental stages performed for each population separately; those not sharing a letter are significantly different. Arrows beneath the x-axis mark the onset of anther dehiscence. Analysis of these data is in Table 2

 

View this table:
[in this window]
[in a new window]
 
Table 2. Analysis of variation in pollen receipt and pollen tube growth among stages of floral development in hand-pollinated flowers from two populations of Aquilegia canadensis. The overall ANOVA model was significant for all variables (all P's < 0.05). The residual df = 108 in all analyses. Means are in Fig. 1

 
Variation in pollen germination and tube growth across floral development was much more striking than the pattern of variation in pollen receipt. ANOVA revealed a strong increase in pollen tube growth from stage 3 to stage 8 for both measures of pollen tube development (Table 2). In contrast, there was no significant variation in tube growth between populations, and the pattern of variation across developmental stages did not differ significantly between populations. Twenty-four hours after hand-pollination, there were very few pollen tubes in the styles of flowers pollinated before anther dehiscence (mean ± 1 SE of stages 3 and 5 = 4.6 ± 1.3 tubes) compared to those pollinated after anther dehiscence (stages 6 and 8 = 33.2 ± 4.8 tubes). Similarly, the mean length of pollen tubes was much lower in flowers pollinated before (0.78 ± 0.14 mm) vs. after (2.04 ± 0.13 mm) anther dehiscence.

Variation in pollen deposition and pollen tube number across floral development in naturally pollinated flowers
In all four populations, the number of pollen grains on naturally pollinated stigmas increased greatly across developmental stages (Fig. 2; P < 0.0001 for all one-way ANOVAs). Overall, there was significant variation in pollen load among developmental stages but not populations (Table 3). On average (±1 SE), only 7.9 ± 1.9 pollen grains were deposited on each stigma by just before anther dehiscence compared to 125.1 ± 21.6 grains just after the onset of anther dehiscence and 239.4 ± 19.1 grains by the end of anther dehiscence. Given that there was an average of 25 ovules in each carpel, it seems that sufficient pollen to fertilize all ovules was received by stigmas within the first day of anther dehiscence. The relation between developmental stage and pollen load varied significantly among populations. However, this interaction effect was small (F = 2.1) compared to the large effect of developmental stage (F = 91.3; Table 3). The largest increase in pollen load occurred between stages 5 and 6 in all populations (increase: range = 5–55-fold, average = 21-fold).



View larger version (31K):
[in this window]
[in a new window]
 
Fig. 2. Variation in pollen deposition and pollen tube growth among stages of floral development in naturally pollinated flowers from four populations of Aquilegia canadensis. Points are means based on ten flowers per stage per population. Error bars have been ommitted to improve clarity. = the number of pollen grains per stigma, {block} = the maximum number of pollen tubes in each style, and {square} = the number of pollen tubes reaching the base of each style. Letters next to points show the results of multiple contrasts among developmental stages; those not sharing a letter are significantly different. Arrows beneath the x-axes mark the onset of anther dehiscence. Analysis of these data is in Table 3

 

View this table:
[in this window]
[in a new window]
 
Table 3. Analysis of variation in pollen grains/stigma, pollen tubes/style, and pollen tubes at the base of the style among stages of floral developmental in naturally pollinated flowers from four populations of Aquilegia canadensis. The overall ANOVA model was significant for all variables (all P's < 0.0001). The residual df = 257 in all analyses. Means are in Fig. 2

 
Variation in the number of pollen tubes in the style closely reflected the pattern of pollen deposition (Fig. 2, P < 0.0001 for all one-way ANOVAs). Again, there was significant variation in pollen tube numbers among developmental stages but not populations (Table 3). However for this variable, the increase in tube numbers across floral development did not vary among populations. By the end of stigma presentation, styles contained an average of only 0.39 ± 0.24 pollen tubes compared to 69.5 ± 5.5 by the end of anther dehiscence. The increase in pollen tube numbers between stages generally peaked between stages 6 and 7 and strongly declined between later stages.

Variation in pollen deposition and tube growth across floral development was also reflected in the number of pollen tubes reaching the base of the style (Fig. 2; P < 0.0001 for all one-way ANOVAs). Again, there was significant variation in pollen tube numbers among developmental stages but not among populations; and the interaction between these main effects was not significant (Table 3). No pollen tubes had reached the base of the style in any of the flowers collected before anther dehiscence. The number of pollen tubes at the base continued to increase throughout the period of anther dehiscence and showed little sign of leveling off by stage 9. By the end of anther dehiscence 20.6 ± 1.9 pollen tubes had reached the base of the style. Although pollen grains accumulate very quickly on stigmas, fertilization of all ovules is probably not complete until 4 d after the beginning of anther dehiscence, when flowers have just begun to fall apart.

Variation in nectar availability and components of flower size across floral development in naturally pollinated flowers
In both populations, there was a large increase in nectar availability as flowers developed (Fig. 3; P < 0.001 for both one-way ANOVAs). Overall, nectar volume varied significantly among developmental stages and between populations (Table 4). On average (±1 SE), each spur contained only 0.53 ± 0.12 µL of nectar before anther dehiscence, compared to 3.00 ± 0.17 µL after the onset of anther dehiscence. Average nectar volume was 1.9-fold higher at QBS1 (2.5 ± 0.3 µL) than at CCR1 (1.4 ± 0.3 µL). The relationship between developmental stage and nectar volume also varied among populations. At QBS1 the largest increase in nectar volume between stages occurred between just before and just after anther dehiscence (2.8-fold increase between stages 5 and 6). At CCR1, the largest increase in nectar volume occurred between stages 6 and 7 (2.3-fold increase).



View larger version (27K):
[in this window]
[in a new window]
 
Fig. 3. Variation in nectar availability and components of flower size among stages of floral development in naturally pollinated flowers from two populations of Aquilegia canadensis. Points are means (±1 SE) based on at least ten flowers per stage per population. Letters next to points show the results of multiple contrasts among developmental stages performed for each population separately; those not sharing a letter are significantly different. Arrows beneath the x-axes mark the onset of anther dehiscence. Analysis of these data is in Table 4

 

View this table:
[in this window]
[in a new window]
 
Table 4. Analysis of variation in per-flower nectar availability (µL) and components of flower size (mm) among stages of floral development in naturally pollinated flowers from two populations of Aquilegia canadensis. The overall ANOVA model was significant for all variables (all P's < 0.0001). The residual df = 213 in all analyses. Means are in Fig. 3

 
Developmental variation in flower size reflected the pattern of nectar availability, except that the components of flower size increased more regularly across floral development (Fig. 3). All three measures of flower size showed roughly the same pattern. There was significant variation in all three components of flower size among developmental stages and between populations (Table 4). The most striking temporal variation occurred for sepal spread, which more than doubled between stages 4 and 8, though the increase was more pronounced in QBS1 than in CCR1 (Table 4). Although flower size increased steadily throughout floral development, the onset of anther dehiscence or nectar availability was not associated with a particularly large increase in any component of flower size.

The effect of experimental protogyny on the level of self-fertilization
The three populations chosen for this experiment were predominantly self-fertilizing in 1998. The average s ± 1 SE of unmanipulated (U and X) flowers was 0.771 ± 0.094. None of the differences in s between populations were significant (0.822 ± 0.063, 0.703 ± 0.084, and 0.788 ± 0.100 for QBC1, QCA1, and QWM1, respectively; P = proportion of bootstrap differences above or below zero >0.16 for all comparisons). These population-level estimates of s were well within the range of those obtained for ten other eastern Ontario populations in 1998 (B. Ozimec, C. R. Herlihy, and C. G. Eckert, Queen's University, unpublished data; see also Routley, Mavraganis, and Eckert, 1999 ).

Contrary to expectations, the removal of early-dehiscing anthers did not reduce the level of self-fertilization in any of the three populations (Fig. 4). In QBC1 and QWM1, protogynous flowers actually self-fertilized a little more than control flowers. In QCA1, the maximum-likelihood estimate for P flowers was slightly lower than the estimate for C flowers, but the difference was far from significant. None of the minor differences in selfing between C and L flowers neared significance. The level of selfing averaged across populations was 0.821 ± 0.043 for P flowers, 0.713 ± 0.070 for C flowers, and 0.783 ± 0.063 for L flowers. There was no difference in selfing between control (C) and unmanipulated (U and X) flowers in any of the three populations (Fig. 4; all P's > 0.14).



View larger version (28K):
[in this window]
[in a new window]
 
Fig. 4. The effect of experimental protogyny on the level of self-fertilization in three populations of Aquilegia canadensis. Points are maximum-likelihood estimates (±1 SE) of the proportions of seeds produced through self-fertilization for each emasculation treatment. Protogynous flowers had the first 19 anthers removed before dehiscence; late flowers had the last 19 anthers removed; control flowers had a random 19 anthers removed. The horizontal line is the average level of selfing in the population based on progeny arrays from unmanipulated flowers. The hatched box shows ±1 SE around the population mean. Sample sizes are in the Materials and Methods

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Morphological vs. functional dichogamy
Although dichogamy is one of the most common floral mechanisms in angiosperms (Bertin and Newman, 1993 ), the extent to which dichogamy separates sexual functions in time and the effect of such temporal separation on mating patterns have rarely been investigated. The extent of dichogamy is usually inferred from rough morphological cues such as the relative timing of style elongation, stigmatic swelling, and anther dehiscence (Lloyd and Schoen, 1992 ; Bertin and Newman, 1993 ). However, these cues may often be misleading (e.g., Preston, 1991 ). Such is the case in A. canadensis.

Schneck (1901) described A. canadensis as protogynous because styles elongate to present stigmas well outside the corolla several days before the filaments unfurl and anthers dehisce. The sepals and nectar spurs are almost fully developed during the period of stigma presentation before anther dehiscence. However, results from our greenhouse experiment suggest that stigmas are relatively unreceptive to pollen germination and tube growth until after anthers start to dehisce (Fig. 1b–d). In some species, including other species of Aquilegia, the onset of stigma receptivity is associated with swelling of the stigmatic papillae or some other obvious morphological change to the stigma (Bertin and Newman, 1993 ; Brunet and Eckert, 1998 ). In A. canadensis, stigmas sometimes appear increasingly swollen as floral development progresses, but the change is subtle and inconsistent (see Materials and Methods; C. G. Eckert and C. R. Herlihy, Queen's University, personal observations).

Even though stigmas become receptive at the same time as anthers start shedding self-pollen, floral development in A. canadensis may afford some opportunity for cross-fertilization to occur before self-fertilization. Styles elongate to bring stigmas to their ultimate position with respect to the rest of the flower well before anther dehiscence, and stigmas can receive substantial quantities of cross-pollen (at least via hand-pollination) before they become receptive. If pollen deposited on stigmas before the onset of receptivity remains viable, prior cross-pollination may increase the probability of cross-fertilization. However, our data on the timing of pollen deposition in four natural populations (Fig. 2) suggest that the potential for prior cross-pollination is not realized (see also Eckert and Schaefer, 1998 ). Very little pollen is deposited on stigmas before anther dehiscence.

That stigmas receive little cross-pollen before anther dehiscence may simply be because flowers are not rewarding and/or attractive to pollinators at this time. In one of the populations we examined (QBS1), nectar availability increased substantially after anther dehiscence (Fig. 3a). In the other population examined (CCR1) nectar reward steadily increased throughout floral development but there was no major increase in nectar coincident with anther dehiscence. Because nectar availability (i.e., standing crop) represents the sum of nectar production minus nectar removal by pollinators, it may not reliably indicate the rate of nectar production. However, in our populations of A. canadensis, low pollen loads on stigmas before anther dehiscence suggest very low rates of pollinator visitation. Accordingly, low nectar levels early in floral development are very likely a consequence of low nectar production.

Developmental variation in nectar production has been observed between floral sex phases in several dichogamous taxa, and in most cases variation in nectar production is associated with significant variation in pollinator visitation (Bell et. al., 1984 ; Devlin and Stephenson, 1985 ; Klinkhamer and de Jong, 1990 ; Gonzalez et al., 1995 ; Aizen and Basilio, 1998 ; but see McDade, 1986 ). In A. canadensis, the steady increase in components of flower size, especially sepal spread, throughout floral development (Fig. 3b–d) could be used by pollinators to roughly discriminate between rewarding and unrewarding flowers. Our field observations also indicate that some pollinator taxa may also visit A. canadensis to collect pollen (see also Macior, 1966 ). This too may result in flowers being visited infrequently before anthers start shedding pollen (Gonzalez et al., 1995 ; Inoue, Maki, and Masuda, 1995 ).

Another explanation for low pollen deposition before anther dehiscence is that cross-pollination is simply infrequent throughout floral development regardless of rewards or flower size. In other words, the large increase in pollen deposition following anther dehiscence observed in natural populations was due almost entirely to self-pollination. This is especially likely because A. canadensis has a well-developed capacity for automatic self-pollination (Macior, 1978 ; Eckert and Schaefer, 1998 ; Routley, Mavraganis, and Eckert, 1999 ). However, we do not have data on the relative deposition of cross- vs. self-pollen required to directly examine this possibility. However, results from a concurrent emasculation experiment in QBS1 indicated that wholly emasculated flowers receive enough allogamous pollen to achieve the same level of seed set as intact flowers (Eckert and Schaefer, 1998 ). This suggests that, at least in this population, there is a major increase in the rate of cross-pollination after anther dehiscence.

Would protogyny reduce self-fertilization if it occurred?
The simultaneous onset of stigma receptivity and anther dehiscence combined with close proximity of dehiscing anthers and receptive stigmas would seem to account for the very high levels of self-fertilization observed in eastern Ontario populations of A. canadensis. Self-fertilization estimated from progeny arrays in ten populations varied widely among populations (17–100%) and averaged 75% (Routley, Mavraganis, and Eckert, 1999 ). Comparison of the inbreeding coefficients of adult plants with those expected at inbreeding equilibrium (Ritland, 1990 ) further indicated that inbreeding depression is probably very strong ({delta} = 1 - [fitness of selfed progeny/fitness of outcrossed progeny] = 0.89; Routley, Mavraganis, and Eckert, 1999 ). Recent experimental results also suggest that although automatic self-fertilization may be favored to some extent as a mechanism of reproductive assurance, it results in severe seed discounting and a substantial net reduction in fitness (C. R. Herlihy and C. G. Eckert, Queen's University, unpublished data).

Given that high levels of selfing amount to a major fitness cost, it would seem that any mechanism that reduces the level of selfing would be strongly selected in A. canadensis. This led us to ask whether protogyny could potentially be an effective anti-selfing mechanism in this species. This "what if" approach seemed relevant for two reasons. First, many researchers view protogyny as a mechanism that might have evolved to reduce self-fertilization because it provides opportunities for cross-pollination before any possibility of self-pollination (Lloyd and Webb, 1986 ; Bertin and Newman, 1993 ). Second, the extent and nature of dichogamy appear to be evolutionarily labile in Aquilegia (see Introduction). For instance, A. elegantula resembles A. canadensis in terms of both floral morphology and pollination ecology but is strongly protogynous (Miller, 1978 ).

We created protogynous flowers in three populations by removing the first 19 anthers that would have dehisced and compared them to nonprotogynous control flowers from which 19 randomly chosen anthers had been removed. Contrary to our expectations, protogynous flowers did not self-fertilize at a lower level than control flowers in any of the populations, even though our experimental manipulation provided ~1 d during which cross-pollen could arrive on stigmas before any opportunity for self-pollination. We also did not detect any difference in self-fertilization between control flowers and unmanipulated flowers. In what follows, we consider three explanations for this unexpected result.

Protogyny, as we have attempted to simulate it, is primarily effective at reducing self-pollination occurring between anthers and stigmas in the same flower (autogamy). If most self-fertilization involved between-flower self-pollination (geitonogamy) protogyny would have little effect on the mating system. Geitonogamy may be a common form of self-fertilization in plants (de Jong, Waser, and Klinkhamer, 1993 ; Harder and Barrett, 1996 ; Snow et al., 1996 ), and the few studies that have quantified the contribution of geitonogamy to self-fertilization in natural populations have all revealed substantial geitonogamous selfing (Schoen and Lloyd, 1992 ; Leclerc-Potvin and Ritland, 1993 ; Eckert, 2000 ). Preliminary data from an emasculation experiment conducted in nine populations of A. canadensis within several kilometres of the ones involved in this study suggest that geitonogamy causes about one-third of the overall level of selfing (C. R. Herlihy and C. G. Eckert, Queen's University, unpublished data). This suggests that most selfing occurs via autogamy, thus the expectation of lower selfing in protogynous flowers compared to control flowers still holds.

The incomplete dichogamy created by our floral manipulations does not necessarily reduce the overall level of self-pollination; it simply gives cross-pollen a head start in competition for fertilizations. Consequently, protogyny will have little effect on the mating system if the rate at which cross-pollen arrives on stigmas is very low. We currently have no data on the rate of cross-pollination in natural populations. The extremely high level of self-fertilization observed in eastern Ontario populations of A. canadensis suggests that cross-pollination is infrequent. However, Eckert and Schaefer (1998) showed that, in QBS1, emasculated flowers eventually receive enough cross-pollen for full seed set. Detailed investigation of the relative timing of self- and cross-pollination will be required to further evaluate this explanation.

Finally, protogyny will be relatively ineffective at reducing selfing if pollinators tend to avoid flowers in female phase (Klinkhamer and de Jong, 1990 ; Inoue, Maki, and Masuda, 1995 ). This is particularly likely to occur if pollinators collect pollen (Gonzalez et al., 1995 ) or flowers produce less nectar during the female than male phase (Bell et al., 1994 ; Aizen and Basilio, 1998 ). Pollinators visit flowers of A. canadensis infrequently in our study populations. However, the observations we have made to date indicate that flowers are primarily visited by a few bumble bee (Bombus) species that either collect pollen, probe for nectar, or rob nectar by biting holes in the spur tips. Bees collecting pollen may avoid our protogynous flowers during their enforced female phase. Likewise, nectar-collecting bees may avoid female-phase flowers if they associate conspicuous, dehiscing anthers with the presence of nectar in the spurs (see Fig. 3a; but see Bell at al., 1994 ). Currently, we do not have the data required to rigorously address this explanation, although we have observed both pollen- and nectar-collecting bees visiting emasculated flowers. It is worth noting that none of our potential explanations for the ineffectiveness of experimental protogyny at reducing selfing involves methodological artifacts. We also expect that naturally occurring protogynous mutants would fail to enjoy the increased progeny quality associated with reduced self-fertilization. Perhaps this is why protogyny has not evolved in eastern Ontario populations of A. canadensis.

Although the functioning of dichogamy depends on how pollinators interact with flowers, much of what we know about this widespread floral mechanism is based on relatively crude observations of the timing of stigma presentation vs. anther dehiscence (Lloyd and Webb, 1986 ; Bertin and Newman, 1993 ). Few studies have combined these morphological observations with quantitative assessment of important components of stigma receptivity such as pollen adhesion, germination and pollen tube growth (e.g., Preston, 1991 ; Navarro, 1997 ). Fewer still have examined whether temporal separation of anther dehiscence and stigma receptivity reduces the potential for self-fertilization or other forms of pollen-pistil interference in natural populations (e.g., Jonsson, Rosquist, and Widén, 1991 ; Dinnétz, 1997 ). Taken together, our results emphasize the importance of field observations and experiments in evaluating the functional significance of dichogamy and other floral traits in natural populations of flowering plants.


    FOOTNOTES
 
1 The authors thank Greg Grabas, Agnes Kliber, Grace Lo, and Matthew Routley for help with field work; Celine Muis, Amy Schaefer, and Matthew Routley for help in the lab; the Queen's University Biological Station (QUBS) and the Saint Lawrence Islands National Park for logistic support in the field and access to study sites; Bob Montgomerie and Dennis Jelinski for helpful suggestions; Chris Herlihy, Agnes Kliber, and Matthew Routley for comments on the manuscript; Queen's University School of Graduate Studies for a scholarship to KM; and the National Sciences and Engineering Research Council of Canada (NSERC) for a research grant to CGE. Back

2 Current address: Department of Botany, University of Toronto, Toronto, Ontario, Canada M5S 3B2; e-mail: griffin{at}botany.utoronto.ca Back

3 Author for correspondence (phone: 613-533-6158, FAX: 613-533-6617, e-mail: eckertc{at}biology.queensu.ca ) Back


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Aizen, M. A., and A. Basilio. 1998 Sex differential nectar secretion in protandrous Alstroemeria aurea (Alstroemeriaceae): is production altered by pollen removal and receipt? American Journal of Botany 85: 245–252[Abstract]

Barrett, S. C. H., and C. G. Eckert. 1990 Variation and evolution of mating systems in seed plants. In S. Kawano [ed.], Biological approaches and evolutionary trends in plants, 229–254. Academic Press, Tokyo, Japan

Bell, G., L. Lefebvre, L.-A. Giraldeau, and D. Weary. 1984 Partial preference of insects for the male flowers of an annual herb. Oecologia 64: 287–294[CrossRef][ISI]

Bertin, R. I. 1993 Incidence of monoecy and dichogamy in relation to self-fertilization in angiosperms. American Journal of Botany 80: 557–560[CrossRef][ISI]

———, and C. M. Newman. 1993 Dichogamy in the angiosperms. Botanical Review 59: 112–152[CrossRef][ISI]

Brunet, J., and C. G. Eckert. 1998 Effect of floral morphology and display on outcrossing in blue columbine, Aquilegia caerulea (Ranunculaceae). Functional Ecology 12: 596–606[CrossRef][ISI]

Chase, V. C., and P. H. Raven. 1975 Evolutionary and ecological relationships between Aquilegia formosa and A. pubescens (Ranunculaceae), two perennial plants. Evolution 29: 474–486[CrossRef][ISI]

Cruden, R. W. 1988 Temporal dioecism: systematic breadth, associated traits, and temporal patterns. Botanical Gazette 149: 1–15[CrossRef]

de Jong, T. J., N. M. Waser, and P. G. L. Klinkhamer. 1993 Geitonogamy: the neglected side of selfing. Trends in Ecology and Evolution 8: 321–325

Devlin, B., and A. G. Stephenson. 1985 Sex differential floral longevity, nectar secretion, and pollinator foraging in a protandrous species. American Journal of Botany 72: 303–310[CrossRef][ISI]

Dinnétz, P. 1997 Male Sterility, protogyny, and pollen-pistil interference in Plantago maritima (Plantaginaceae), a wind-pollinated, self-incompatible perennial. American Journal of Botany 84: 1588–1594[Abstract]

Eckert, C. G. 2000 Contributions of autogamy and geitonogamy to self-fertilization in a mass-flowering, clonal plant. Ecology 81: 532–542[CrossRef][ISI]

———, and S. C. H. Barrett. 1994 Post-pollination mechanisms and the maintenance of outcrossing in self-compatible Decodon verticillatus (Lythraceae). Heredity 72: 396–411[ISI]

———, and A. Schaefer. 1998 Does self-pollination provide reproductive assurance in wild columbine, Aquilegia canadensis (Ranunculaceae)? American Journal of Botany 85: 919–924[Abstract]

Faegri, K., and L. van der Pijl. 1979 The principles of pollination ecology, 3rd ed. Pergamon, New York, New York, USA

Gonzalez, A., C. L. Rowe, P. J. Weeks, D. Whittle, F. S. Gilbert, and C. J. Barnard. 1995 Flower choice by honey bees (Apis mellifera L.): sex-phase of flowers and preferences among nectar and pollen foragers. Oecologia 101: 258–264[CrossRef][ISI]

Grant, K. A., and V. Grant. 1968 Hummingbirds and their flowers. Columbia University Press, New York, New York, USA

Grant, V., and E. J. Temeles. 1992 Foraging ability of rufous hummingbirds on hummingbird flowers and hawkmoth flowers. Proceedings of the National Academy of Sciences, USA 89: 9400–9404[Abstract/Free Full Text]

Harder, L. D., and S. C. H. Barrett. 1996 Pollen dispersal and mating patterns in animal pollinated plants. In D. G. Lloyd and S. C. H. Barrett [eds.], Floral biology: studies on floral evolution in animal-pollinated plants, 140–190. Chapman and Hall, New York, New York, USA

Hodges, S. A., and M. L. Arnold. 1994 Columbines: A geographically widespread species flock. Proceedings of the National Academy of Sciences, USA 91: 5129–5132[Abstract/Free Full Text]

Holtsford, T. P., and N. C. Ellstrand. 1992 Genetic and environmental variation in floral traits affecting outcrossing rate in Clarkia tembloriensis (Onagraceae). Evolution 46: 216–225[CrossRef][ISI]

Inoue, K., M. Maki, and M. Masuda. 1995 Different responses of pollinating bees to size variation and sexual phases in flowers of Campanula. Ecological Research 10: 267–273

Jonsson, O., G. Rosquist, and N. Widén. 1991 Operation of dichogamy and herkogamy in five taxa of Pulsatilla. Holarctic Ecology 14: 206–271

Kearns, C. A., and D. W. Inouye. 1993 Techniques for pollination biologists. University Press of Colorado, Niwot, Colorado, USA

Klinkhamer, P. G., and T. J. de Jong. 1990 Effects of plant size, plant density and sex differential nectar reward on pollinator visitation in the protandrous Echium vulgare. Oikos 57: 399–405

Leclerc-Potvin, C., and K. Ritland. 1994 Modes of self-fertilization in Mimulus guttatus (Scrophulariaceae): a field experiment. American Journal of Botany 81: 199–205[CrossRef][ISI]

Lloyd, D. G., and D. J. Schoen. 1992 Self- and cross-fertilization in plants. I. Functional dimensions. International Journal of Plant Sciences 153: 358–369[CrossRef]

———, and C. J. Webb. 1986 The avoidance of interference between the presentation of pollen and stigmas in angiosperms. I. Dichogamy. New Zealand Journal of Botany 24: 135–162[ISI]

Macior, L. W. 1966 Foraging behaviour of Bombus (Hymenoptera: Apidae) in relation to Aquilegia pollination. American Journal of Botany 53: 302–309[CrossRef][ISI]

———. 1978 The pollination of vernal angiosperms. Oikos 30: 452–460[CrossRef][ISI]

Mavraganis, K. 1998 The reproductive ecology and population biology of Aquilegia canadensis (Ranunculaceae). M.Sc. thesis, Queen's University, Kingston, ON, Canada

McDade, L. A. 1986 Protandry, synchronized flowering and sequential phenotypic unisexuality in neotropical Pentagonia macrophylla (Rubiaceae). Oecologia 68: 218–223[CrossRef][ISI]

Miller, R. B. 1978 The pollination ecology of Aquilegia elegantula and caerulea (Ranunculaceae) in Colorado. American Journal of Botany 65: 406–414[CrossRef][ISI]

———. 1985 Hawkmoth pollination of Aquilegia chrysantha (Ranunculaceae) in southern Arizona. Southwest Naturalist 30: 69–76[CrossRef]

———, and C. L. Willard. 1983 The pollination ecology of Aquilegia micrantha (Ranunculaceae) in Colorado. Southwest Naturalist 28: 157–164[CrossRef]

Munz, P. A. 1946 Aquilegia, the cultivated and wild columbines. Gentes Herbarum 7: 1–150

Navarro, L. 1997 Is the dichogamy of Salvia verbenaca (Lamiacae) an effective barrier to self-fertilization? Plant Systematics and Evolution 207: 111–117[CrossRef][ISI]

Neter, J., W. Wasserman, and M. H. Kutner. 1990 Applied linear statistical models, 3rd ed. Richard D. Irwin, Homewood, New Jersey, USA

Payson, E. B. 1918 The North American species of Aquilegia. Contributions from the National Herbarium 20: 133–157

Preston, R. E. 1991 The intrafloral phenology of Streptanthus tortuosus (Brassicaceae). American Journal of Botany 78: 1044–1053[CrossRef][ISI]

Proctor, M., and P. Yeo. 1972 The pollination of flowers. William Collins Sons & Co. Ltd., Toronto, Ontario, Canada

Richards, A. J. 1986 Plant breeding systems. George Allen & Unwin, London, UK

Ritland, K. 1986 Joint maximum likelihood estimation of genetic and mating structure using open-pollinated progenies. Biometrics 42: 25–43[CrossRef][ISI]

———. 1990 A series of FORTRAN computer programs for estimating plant mating systems. Journal of Heredity 81: 235–237[ISI]

Routley, M. B., K. Mavraganis, and C. G. Eckert. 1999 Effect of population size on the mating system in a self-compatible, autogamous plant, Aquilegia canadensis (Ranunculaceae). Heredity, 82: 518–528

SAS. 1997 JMP® user's guide, version 3.2.1. SAS Institute, Cary, North Carolina, USA

Schneck, J. 1901 Notes on Aquilegia canadensis (L.) and A. vulgaris (L.). Botanical Gazette 32: 304–305[CrossRef]

Schoen, D. J. 1982 The breeding