|
|
||||||||
Reproductive Biology |
2Department of Botany, Norwegian University of Science and Technology, N-7491 Trondheim, Norway; 3Institute of Arctic Biology, University of Alaska, Fairbanks, Alaska 99775-7000 USA; 4Jepson Herbarium and Department of Integrative Biology, University of California, Berkeley, California 94720-2465 USA; 5Department of Biological Science, University of Pittsburgh, Pittsburgh, Pennsylvania 15260 USA
Received for publication January 30, 2001. Accepted for publication July 3, 2001.
| ABSTRACT |
|---|
|
|
|---|
Key Words: autogamy correlated evolution Collinsieae cross pollination flower development herkogamy mating system pollination Scrophulariaceae self-pollination
| INTRODUCTION |
|---|
|
|
|---|
We have identified sister genera in the Scrophulariaceae, sensu lato (s.l.), Collinsia Nutt. (1820 species) and Tonella Nutt. ex Gray (two species), both comprising self-compatible North American annuals, as a candidate "model system" for integrative study of the evolution of floral morphology, later stages of floral development, herkogamy, pollination ecology, and timing of self-pollination, using a phylogenetic framework. Results from traditional systematic (Newsom, 1929
; Munz, 1959
; Neese, 1993
), molecular-phylogenetic (B. G. Baldwin, W. S. Armbruster, and B. Wessa, unpublished data), and genetic (Garber, 1958a
; Grant, 1958
) investigations lead us to hypothesize that repeated transitions in mating system have occurred in Collinsia. Changes in flower size (Grant, 1958
), degree of herkogamy, and timing of self-pollination (Rust and Clement, 1977
; Armbruster, 1980
; Kalisz et al., 1999
) appear to have accompanied these putative changes in mating system.
Here, we combine comparative morphological and developmental measurements with molecular phylogenetic information to examine evolutionary patterns and functional consequences of stamen and style elongation during anthesis. We test the following predictions: (1) Large-flowered taxa have delayed self-pollination, and small-flowered taxa have early self-pollination. (2) Selection for prolonged and consistent anther-stigma separation (herkogamy) has occurred in large-flowered, cross-pollinating species. This selection has led to more precise positioning of dehiscing anthers ("relative developmental precision") and hence to a proportionately smaller region within the lower corolla lobe in which pollen is shed (the "pollen zone"). (3) Small-flowered, self-pollinating populations should have proportionally smaller anthers than outcrossing populations, because self-pollinators have lower optimal pollen-ovule ratios, smaller pollen, or both (Lloyd, 1965
; Baker, 1967
; Cruden, 1973, 1977, 2000
). (4) Flower size, anther size, developmental precision, and time of anther-stigma contact vary continuously instead of falling into discrete inbreeding and outcrossing "syndromes" (Vogler and Kalisz, 2001
; but see Lande and Schemske, 1985
, and Schemske and Lande, 1985
). Further, these traits should covary because of their functional relationships (e.g., Armbruster and Schwaegerle, 1996
).
| MATERIALS AND METHODS |
|---|
|
|
|---|
Flowers in both genera are zygomorphic and consist of a five-lobed calyx, five-lobed corolla, four epipetalous stamens, and one pistil, containing 2 to
16 ovules (Fig. 1). The corollas of Collinsia are unique among the Scrophulariaceae s.l. in resembling pea flowers: an upper lip of two lobes forms the "banner," and a lower lip of three corolla lobes includes a pair of wings and a folded (conduplicate) keel, enveloping the style and stamens. At the base of the banner, wings, and keel, the corolla is constricted into a narrow aperture, thus forming a constricted "mouth" at the top of a saccate tube (Fig. 1). Tonella flowers are similar but are more open and lack the banner, folded keel, and constricted mouth. The stigmas are receptive to pollen-tube growth either early in anthesis, in Tonella and most Collinsia species, or late in anthesis in a minority of Collinsia species (see below). In Collinsia flowers, nectar is secreted by a small nectary (probably derived from the staminode) at the base of the asymmetrically saccate corolla tube (Fig. 1) and is accessible to long-tongued bees that can reach the nectar by inserting their tongues through the opening in the constricted mouth at the base of the banner (Rust and Clement, 1972, 1977
; Armbruster, 1980
).
|
Pollinating bees at elevations below
1000 m include species of Bombus (Apidae), Anthophora, Emphoropsis, and Synhalonia (Anthophoridae), and Osmia (Megachilidae). At higher elevations species of Osmia are the main pollinators (Armbruster, 1980
). These bees land on the flowers and stand on the platform-like pair of wings with their heads at the base of the banner, obtaining nectar by inserting their proboscides through the constricted corolla mouth. Unlike in most other Scrophulariaceae s.l., the bees do not brush across the anthers and stigma as they move into the flower. Instead, as on a legume flower, the bees land directly in the foraging position and are excluded from entering the floral tube by the appressed lips of the corolla mouth. One exception is the oligolectic specialist Osmia glauca (Rust and Clement, 1972
), which is small enough to be able to enter the saccate corolla tube to reach the nectar of some large-flowered species (Armbruster, 1980
). However, this activity does not cause the depression of the keel or contact with anthers and stigmas; the latter only occurs when the bee is "properly" positioned on the platform (Armbruster, 1980
). Other pollinating bees are large enough that their mass immediately depresses the keel and wings, but not the stamens and style. Hence pollen deposition and/or pickup occurs when the stamens and style move upwards, relative to the bee, and anthers and/or stigma contact the bee in a consistent position on the underside of its head, thorax, or abdomen (depending on species of bee and Collinsia; Armbruster, 1980
). The clump of pollen placed on the bee is generally 13 mm in length. In addition to obtaining nectar, many bees also collect pollen by manipulating the anthers with their hind legs while the keel is depressed (Rust and Clement, 1972, 1977
; Armbruster, 1980
).
Phylogenetic methods
Leaf material was collected from
10 individuals from each of 19 populations representing 18 Collinsia species and one population each of the two Tonella species (Appendix 1: URL: http://ajbsupp.botany.org/v89/armbruster/) and kept on ice until returned to the laboratory, where samples were kept at 80°C. Total DNA of each sample was extracted from leaves using a modification of Doyle and Doyle's (1987)
CTAB (hexadecyltrimethylammonium bromide) procedure (plus a phenol extraction, RNase digestion, and two ethanol precipitations). The 18S26S nuclear ribosomal DNA internal transcribed spacer (ITS) region (i.e., ITS-1, 5.8S, ITS-2) was PCR-amplified and directly cycle-sequenced using the methods of Baldwin and Wessa (2000)
. Sequences of both strands of the ITS-region were resolved on 4.8% polyacrylamide gels using an ABI 377 automated sequencer (PE Applied Biosystems, Foster City, California, USA). Sequences of the ITS-region of Collinsia and Tonella were unambiguously aligned manually with outgroup sequences of Chelone L., Keckiella Straw, and Penstemon Schmid. generously provided by Andrea Wolfe. Representatives of tribe Chelonieae were chosen as the outgroup based on molecular phylogenetic evidence that tribe Chelonieae is sister to tribe Collinsieae, comprising Collinsia and Tonella (Wolfe et al., 1997
). Parsimony analyses were conducted using PAUP* 4.0 (Swofford, 1998
) using the heuristic option and 100 random addition sequences of the taxa. Analyses were conducted on the entire aligned sequence matrix plus indel characters (recoded using "simple indel coding"; Simmons and Ochoterena, 2000
), with all characters and character-state transformations given equal weight. We estimated reliability of clades by bootstrap and decay analyses, using heuristic searches (20 random addition sequences of the taxa) for each of the 100 bootstrap replicates and for the decay analysis. Decay values (Bremer support) for each clade were assessed using the reverse constraints approach implemented in AutoDecay 4.0 (Eriksson, 1998
).
Floral development
Flowers from individuals of 22 populations in 16 species of Collinsia and one population from each of the two Tonella species were collected between 1976 and 1999 (Appendix 1: URL: http://ajbsupp.botany.org/v89/armbruster/), fixed in FAA (1 part formalin, 1 part acetic acid, 18 parts ethanol), and stored in 70% ethanol. Some morphological and ecological information on two additional Collinsia species (C. concolor, C. parryi) was obtained from the literature (Newsom, 1929
; Munz, 1959
). Flower samples were examined under a dissecting microscope and measured with digital calipers or an ocular micrometer. Flowers were classified into five stages based on the number of anthers dehisced (zero to four anthers). For each population we measured six to ten flowers (each usually from a separate plant) at each of the five developmental stages, recording keel length, corolla-tube length (mouth to nectary), length of the style beyond the mouth of the corolla tube, and length of each stamen beyond the mouth of the corolla tube (to the tip of each anther). Anther length was measured on a sample of 12 or more anthers per population. Although curvature of organs was minimal except when very young, all measurements were taken on straightened organs to improve repeatability.
Similar measurements were conducted on flowers of plants grown in the greenhouse for comparison with the measurements of preserved flowers and to obtain information about the actual duration of anthesis and each stage. Species examined in the greenhouse were: C. childii, C. heterophylla, C. linearis, C. parviflora, C. rattanii, C. sparsiflora, C. torreyi, and C. verna.
On preserved flowers we also measured the distance from the base of the anther of the shortest stamen with a dehisced anther to the apex of the anther of the longest stamen with a dehisced anther to estimate the length of the pollen "zone" at each floral stage. Thus, the length of the pollen zone is a function of both anther size and the positions of the dehisced anthers. The size of the pollen zone affects the degree of herkogamy because when the stigma is within this region it is subject to spontaneous self-pollination. As the style elongates, it drives the stigma through the pollen zone, effecting self-pollination. The size and position of the pollen zone also influence the size and position of the clump of pollen deposited on pollinating bees. We compared the length of the pollen zone at floral stage 3 across populations because this is the last stage at which herkogamy is maintained in any species. As flower size increases, the absolute size of the pollen zone could be expected to increase through isometric/allometric scaling effects. Hence, selection for increased herkogamy, and hence closer placement of anthers (greater "precision") in large-flowered species, should reduce the slope of the allometric relationship of pollen-zone size to flower size; this hypothesis was tested using regression of log-transformed data (Sokal and Rohlf, 1981
).
The relative developmental precision index (DPI) was calculated as DPI = 1 [(S13 S33)/(S11 S31)], where (S13 S33) is the distance between the tip of the first dehisced anther (S1) and the tip of the third dehisced anther (S3) at stage 3 and (S11 S31) is the corresponding distance at stage 1. The DPI measures the degree of precision in the position of the dehiscing anthers relative to the positions of the undehisced anthers, and it varies from <0 (developmental divergence) to 0 (parallel stamen growth) to 1 (developmental convergence/precision). This index is mensurally independent of both anther size and flower size. We calculated DPI for floral stage 3 because we estimated pollen-zone lengths at this stage.
We used two approaches to estimate the stage at which the stigma contacts the dehiscing anthers in each population. The first measure of the timing of anther-stigma contact was based on the visual classification of each flower as having the stigma in contact with, or not in contact with, dehiscing anthers and used logistic regression (SAS PROBIT procedure; SAS Institute, 1996
) to estimate the stage of development at which 50% of the flowers in a sample had stigmas and dehiscing anthers in contact (ASC-50). The second approach to assessing timing of anther-stigma contact was based on graphs of style and stamen elongation. We determined the average stage at which the style entered the pollen zone (see RESULTS). These graphs were based on means of flower parts at each stage of development (five to ten flowers at each stage per population), and the metric is therefore sensitive to sample-size problems and differences in mean flower size of different floral cohorts. The two indices were highly correlated (r = 0.90, P < 0.001), so only ASC-50, which has known statistical properties (Sokal and Rohlf, 1981
), was used in statistical analyses.
Stigma receptivity
We determined stigmatic receptivity of fresh flowers at each of the five stages of anther dehiscence in two ways. First, in the field we tested for stigmatic peroxidase activity (SPA) using the method of Kearns and Inouye (1993)
. Intact styles were placed on a glass slide in a drop of 3% hydrogen peroxide and covered with a cover slip. Bubble production from the stigma within 23 min indicates the activity of stigmatic peroxidases and hence receptivity of the stigma. The stage at which 50% of the stigmas in the sample tested positive for peroxidase activity (SPA-50) was calculated using logistic regression (SAS PROBIT procedure, SAS Institute 1996
) and was our metric for population comparisons. Second, to verify the relationship of a positive SPA and actual receptivity to pollen-tube growth, generally 510 styles at each floral stage for one population each, from a subset of species and varieties, were preserved in ethyl alcohol, and, in the laboratory, were cleared, stained, and examined for pollen tubes using epifluorescence microscopy (see Kalisz et al. [1999]
for details on the method).
Comparative statistical analyses
Statistical analysis of among-population relationships among traits was conducted using correlation analyses, factor analysis (principal components analysis [PCA] and Varimax rotation), and path analysis (Li, 1975
) implemented with regression models in Statistica (StatSoft, 1994
). All length measurements were log transformed to assess proportional variation and to correct for heteroscedasticity (Sokal and Rohlf, 1981
). Correlation, factor, and path analyses were conducted on phylogenetically corrected (PC) independent contrasts (see below) and, for comparison, on the original "phylogenetically naïve" data (see Harvey and Pagel, 1991
; Armbruster, 1992
). This comparison gives insights not only into the value of using phylogenetic information in the analysis, but also into the nature of the phylogenetic "signal."
Independent contrasts were calculated using the program CAIC (Purvis and Rambaut, 1995
) for calculating all the paired contrasts (differences) of trait values on the molecular phylogenetic tree. In a few cases we used conspecific or con-subspecific populations as phylogenetic pairs even though one (or both) of the populations was not represented in the molecular tree. This procedure was justified by the obvious close relationship between conspecific populations (all conspecific populations treated this way shared numerous morphological synapomorphies). Only two populations of minimal taxonomic rank were used (the two with largest sample sizes) in order to have unambiguous sister relationships.
Contrasts of continuous traits were analyzed using path analysis (Li, 1975
) and multiple regression through the origin (Purvis and Rambaut, 1995
) with Statistica (StatSoft, 1994
). Simple paired comparisons of extant sister taxa were calculated using a modification of Ridley's (1983)
method, employing the paired-sample t test (Sokal and Rohlf, 1981
).
An index of relative expected outcrossing was calculated as the mean of the relative flower size and the relative ASC-50. The variables were relativized by dividing each observed value by the maximum observed value for that variable. Thus the index varies from near zero to one.
| RESULTS |
|---|
|
|
|---|
|
|
|
0.5 in the earliest self-pollinators, such as C. callosa, to 3.05 in the large-flowered C. tinctoria (Table 1).
|
|
Promotion of cross-pollination is less clear in species with late-receptive stigmas. In a few species (e.g., C. heterophylla) the stigmas are receptive shortly before anther-stigma contact, providing a short period during which cross-pollination can occur before self-pollination. In other species (e.g., C. multicolor, C. tinctoria), receptivity is nearly simultaneous with, or in two cases (C. corymbosa, C. verna) shortly after, anther-stigma contact. This last developmental condition, in particular, appears to preclude cross-pollination prior to self-pollination, yet C. verna is moderately outcrossing (Kalisz et al., 1999
), as is C. heterophylla (Mayer, Charlesworth, and Meyers, 1996
). Thus the length of the period between stigma receptivity and anther-stigma contact only partly accounts for the dynamics of pollen arrival and likelihood of cross-fertilization. Interestingly, stigma flaring generally occurred early in floral development, even in species with late receptivity (Table 1, see also Kalisz et al., 1999
).
Trait covariance
Inspection of the patterns of floral development during anthesis (Fig. 3, Table 1) suggests that populations and species with large flowers (Fig. 3, right column) have more closely placed (precise) anthers at dehiscence than do populations and species with smaller flowers (Fig. 3, left column). Upper (initially shorter) stamens tend to grow more than lower (initially longer) stamens prior to anther dehiscence in large-flowered species, and hence the anther positions tend to converge with filament elongation. The resulting pollen zone is smaller than would result from isometric development. The trajectories of stamen growth tend to be closer to parallel in smaller-flowered populations and species. In each of the five Collinsia sister-species/population pairs (C. torreyi var. wrightii vs. C. t. var. torreyi; C. rattanii vs. C. linearis; C. parviflora vs. C. grandiflora; C. sparsiflora var. collina vs. C. s. var. arvensis; and C. heterophylla, population 4 vs. C. heterophylla, population 2) the larger-flowered species or population had more closely positioned anthers (greater developmental-precision index) than the small flowered species or population (see Table 1; phylogenetically corrected paired comparison test; t = 7.35, P < 0.002).
However, species with larger flowers generally had larger pollen zones in absolute size than species with small flowers, as might be expected from allometric/isometric scaling (see Table 1), although this pattern was not significant (paired comparison test; t = 1.10, P > 0.1). This tendency towards increased absolute size of the pollen zone indicates that increased relative precision does not usually fully compensate for the positive allometric effect of increasing flower and anther size. Nevertheless, the slope of the logarithmic scaling of pollen zone with flower size was 0.50, which was significantly shallower than the isometric scaling of pollen zone with flower size (i.e., ß = 1.0; P < 0.001; Fig. 5). Anther length also increased allometrically with flower size, but its slope (0.74) was not significantly shallower than isometric or significantly steeper than that of the pollen zone (Fig. 5). These relationships suggest that smaller-than-expected (isometrically) pollen zones in large flowers are the result of more closely placed anthers rather than smaller anthers.
|
Correlation analysis showed a complex pattern of intercorrelation among the various morphological and developmental variables. The naïve and phylogenetically corrected (PC) correlation analyses gave some strikingly different results (compare above vs. below diagonal, Table 3). Phylogenetically naïve analyses indicated that the time of stigma receptivity (SPA-50) was significantly positively correlated with flower size, anther length, size of pollen zone, and ASC-50 and was significantly negatively correlated with relative precision (Table 3). Incomplete sampling precluded PC analysis of the covariance patterns of SPA-50, so we cannot be certain that these trends would remain after correcting for phylogeny.
|
|
Path analysis of the phylogenetically naïve data showed a strong positive correlation between the time of anther-stigma contact and flower size (Fig. 6A). Relative precision of stamen position responded positively to ASC-50 but not to flower size. Ovule number and anther size both increased with flower size, but there was no discernible effect of ovule number on anther size. The size of the pollen zone was strongly affected by anther size.
|
| DISCUSSION AND CONCLUSIONS |
|---|
|
|
|---|
2 mm in length and slightly zygomorphic in Tonella to nearly 20 mm in length and extreme zygomorphy in Collinsia). All species showed similar patterns of elongation of stamens and dehiscence of anthers, with one (or rarely two) anthers dehiscing at a time during the period of anthesis. Most species also showed a gradual elongation of the style during anthesis. Species varied in the floral stages during which stigma receptivity first occurred and during which stigmas came in contact with anthers, and in degree of convergence of anther position during stamen elongation and maturation. Much of this variation in developmental patterns is assignable to two syndromes: small-flowered species with widely spaced anthers at dehiscence (low developmental precision) and early self-pollination and large-flowered species with closely spaced anthers at dehiscence (high developmental precision) and late self-pollination. This variation is to a surprising extent independent of phylogeny: sister species often belong to opposite syndromes.
Variation in timing of stigma receptivity was incompletely assessed here, but enough data are available to reveal a pattern and raise a question. The data presented here support our interpretation that stigmatic peroxidase activity corresponds closely with stigma receptivity to pollen germination (see also Kalisz et al., 1999
), although additional data are needed. All populations and species studied with early anther-stigma contact showed early peroxidase activity (i.e., stigma receptivity) and generally early pollen-tube growth. These small-flowered populations are almost certainly predominantly self-pollinating; self-pollen arrives and germinates on the receptive stigma shortly after the flower opens.
Several studied populations with medium to late anther-stigma contact have stigmas that become receptive well before anther-stigma contact. This developmental pattern presumably creates an herkogamous period (dehiscing anthers and receptive stigma spatially separated) during which cross-pollination can occur prior to self-pollination (Rust and Clement, 1977
; Armbruster, 1980
). Self-pollination occurs with anther-stigma contact, but long after stigmas become receptive to pollen-tube growth. Hence, the likelihood of experiencing cross-fertilization should be much greater for these large-flowered populations than for populations with early anther-stigma contact.
Finally, a few large-flowered taxa studied have stigmas that become receptive at about the same time as, or shortly after, anther-stigma contact. This pattern, seen most strongly in C. verna (see also Kalisz et al., 1999
), C. corymbosa, and C. tinctoria, suggests that self-pollination (1) occurs as early as or even before cross-pollination and (2) occurs later than indicated by ASC-50. Absence of an herkogamous period (during which cross-pollination could precede self-pollination) in these large-flowered species is puzzling, especially because C. verna, at least, has moderately high rates of outcrossing (about as high as C. heterophylla, which appears to have an herkogamous period). More study is needed before we can understand the dynamics of cross- vs. self-pollen arrival in these species.
One possible mechanism that may promote outcrossing in populations that have stigmas becoming receptive at the same time as, or after, anther-stigma contact involves the flaring of the stigmatic lobes and expansion of papillae prior to physiological receptivity. Early flaring appears to occur in most Collinsia species (Table 1; Kalisz et al., 1999
). Cross pollen may arrive and lodge on these stigmas after flaring but before receptivity (and thus before anther-stigma contact), although the pollen apparently does not germinate. Upon receptivity, cross-pollen tubes may grow down the style prior to deposition of self-pollen on the stigma. Alternatively, cross and self-pollen may lodge on the receptive stigma simultaneously, but cross-pollen tubes may grow faster than the self-pollen tubes and reach the ovules first. Germination of pollen after accumulation of a pollen load has been observed in another member of the Scrophulariaceae s.l. (Stewart, Stewart, and Canne-Hilliker, 1996
) and a lily (Kingston, 1998
).
A possible advantage of receiving pollen on the stigmas over a protracted period but delaying receptivity until near the end of anthesis is that it may allow a large pollen load to accumulate on the stigma prior to pollen germination, leading to more intense pollen competition and higher offspring fitness. This possible advantage may hold both for pollen arriving from genetically diverse fathers (e.g., Snow and Spira, 1996
) and for self-pollen, with pollen competition possibly reducing the effects of inbreeding depression (W. S. Armbruster and D. G. Rogers, unpublished data).
Trait covariance
Correlation, factor, and path analyses all helped to describe the covariation between morphological and developmental traits. Generally, the results of naïve and phylogenetically corrected analyses were similar, which suggests that most of the morphological and developmental traits are evolutionary labile and not constrained by factors associated with phylogeny. This interpretation is consistent with sister taxa often having contrasting flower sizes and probably mating systems. "Phylogenetic lag" is suggested when phylogenetically corrected analyses detect relationships not detected by naïve analysis. This reflects the influence of ancestral character states on current phenotype in addition to selection. In contrast, phylogenetic inertia (stasis) and pseudoreplicatory sampling is indicated when naïve analyses yield (mistakenly) significant results not shown by the corresponding phylogenetically corrected analyses. This reflects the influence of ancestral character states alone.
While the factor analysis yielded similar results for both naïve and phylogenetically corrected data, the correlation and path analyses showed several differences between the two types of data. For example, PC path analysis indicated significant positive covariance between anther size and ovule number and a significant decrease in ovule number associated with delayed anther-stigma contact, relationships not detected by the naïve analysis. (Both are predicted from the theory of optimal pollen/ovule ratios; Cruden, 1977
.) These differences suggest that ancestral traits have influenced, but not prevented, the response to selection for optimal pollen/ovule ratios (phylogenetic lag).
The PC path analysis also indicated that relative precision in placement of dehiscing anthers varied in response to flower size rather than to ASC-50, whereas the naïve analysis indicated the reverse. However, this difference is more simply interpreted as a result of analytical instability generated by colinearity between ASC-50 and flower size (r = 0.78/0.72), rather than as the effect of phylogenetic inertia and/or pseudoreplication.
Tests of predictions
Our first prediction, that large-flowered species have delayed self-pollination was clearly supported, with highly significant PC and naïve correlations between the two traits (PC r = 0.72). This result leads us to conclude that small-flowered taxa are largely self-pollinating inbreeders and large-flowered taxa are more often outcrossing. Preliminary genetic data on mating systems support this interpretation (Weil and Allard, 1964
; Charlesworth and Mayer, 1995
; Mayer, Charlesworth, and Meyers, 1996
; Kalisz et al., 1999
; A. Kahler, University of California, Davis, personal communication, 1980)
Our second prediction, that large-flowered taxa have greater relative precision in anther position than small-flower taxa at the end of stamen elongation, was supported by both PC and naïve analyses. Thus, it appears that more consistent anther-stigma separation (herkogamy) is associated with delayed selfing, higher levels of outcrossing, and greater reliance on pollinators. We also found that the absolute size of the pollen zone increased through allometric scaling with flower and anther size, but significantly less steeply than the isometric relationship. This allometric relationship suggests also that selection for greater herkogamy, and hence a proportionally smaller pollen zone, has operated on larger flowered, outcrossing populations.
The third prediction, that self-pollinating populations will have smaller stamens and/or more ovules (i.e., lower pollen-ovule ratios) than cross-pollinating populations, was supported by the positive relationship between anther size and flower size in both phylogenetically corrected and naïve analyses. However, these results could also be interpreted as a result of scaling. In neither analysis were larger anthers associated with later anther-stigma contact after the effect of flower size was removed. Hence, the biological significance of these covariance relationships is difficult to assess.
Ovule number generally increased with flower size (strongly so in naïve analysis but only weakly so in the PC analysis). To the extent that flower size is a measure of mating system, this positive association between ovule number and flower size is the opposite of our prediction but consistent with flower-size scaling. However, in the phylogenetically corrected analysis ovule number was negatively related to the timing of anther-stigma contact, after the effect of flower size was removed. In other words, late-selfing species tended to have fewer ovules than earlier-selfing species of similar flower size. This pattern is consistent with the prediction of higher pollen/ovule ratios with greater outcrossing. Anther size was positively related to ovule number after removing the effect of flower size (P < 0.01; PC analysis), also as expected.
The fourth prediction was that flower size, anther size, precision of anther placement, and time of anther-stigma contact vary continuously and are intercorrelated. Indeed we found continuous variation in each trait and positive covariance between them. Although clusters of populations occurred at either end of the hypothesized mating-system continuum, over half the populations occurred in the middle of the distribution (Fig. 7). We suggest two possible explanations for this continuous variation: (1) that selection is driving populations towards the two extreme syndromes (two adaptive peaks), but genetic constraint or phylogenetic lag (manifested as a partial inertial pattern on the phylogeny) causes most populations to fall in between the extremes or (2) that continuous variation occurs in the fitness function (e.g., an adaptive ridge or cordillera) because correlational selection (Endler, 1976, 1995
; Armbruster and Schwaegerle, 1996
) acts on the various floral traits affecting mating system, and both extreme and intermediate mating systems represent stable optima. The limited phylogenetic inertia seen in this study system tends to support the first hypothesis, but the obvious evolutionary lability of traits that seem to influence mating system even more strongly supports the second hypothesis.
|
Various lines of evidence suggest that the above and related floral traits covary adaptively rather than as a result of genetic or developmental constraints. Most simply, time of anther-stigma contact and relative precision are size-independent traits, and yet they strongly covary positively and negatively, respectively, with flower size in the phylogenetically corrected analyses, as expected from our adaptive-covariance hypothesis. Second, by holding flower size statistically constant in path analysis we can assess the relationship between time of self-pollination (ASC-50) and ovule number more or less independently of confounding effects of flower size. Doing so, we see a significant negative relationship between time of self-pollination (ASC-50) and ovule number, as expected. Similarly, after correcting for phylogeny, the relationship between anther size and ovule number is positive even when flower size is held constant, as expected from the adaptive covariance hypothesis.
The adaptive scenarios presented above relating to covariation of mating system, flower size, anther size, ovule number, relative precision in anther position during anthesis, and size of the pollen zone are reasonable both in the context of evolution toward more cross-pollination and in the context of evolution toward more self-pollination. However, concluding that the above trait covariance reflects correlated response to selection for more faster, economical flowers makes sense (the above counterevidence notwithstanding) only if predominant selfers have evolved from cross-pollinating species rather than vice versa. Thus, resolving the historical events that have led to the present distribution of mating systems across the phylogeny of tribe Collinsieae may allow discrimination between hypotheses explaining floral trait covariation. More detailed phylogenetic information may also allow us to ascertain the order of change of functionally related traits (see Donoghue, 1989
; Armbruster, 1992
; Pagel, 1994, 1997
).
In summary, floral traits in Collinsieae tend to be integrated into covarying suites, probably associated with high, intermediate, and very low outcrossing rates in different populations. Populations with small flowers have early anther-stigma contact, early stigma receptivity, and parallel stamen elongation leading to widely spaced dehiscing anthers. These flowers are almost certainly mostly autogamous. In contrast, populations with large flowers have delayed anther-stigma contact, often delayed stigma receptivity, and convergent stamen elongation leading to closer positioning of the dehiscing anthers. These larger-flowered populations appear to experience more cross-pollination and probably have outcrossing or mixed mating systems. Thus, we found a complex web of covariation that could best be interpreted as the result of multitrait adaptation for early selfing and high levels of autogamy and/or for delayed selfing and higher levels of outcrossing. However, the variation in these traits is continuous, indicating either continuous variation in selective optima or the combined effects of divergent selection (bimodal fitness optima) and phylogenetic lag. Evaluating the importance of these two processes remains a challenge for future comparative studies of tribe Collinsieae.
| FOOTNOTES |
|---|
6 Author for reprint requests (scott.armbruster{at}chembio.ntnu.no
). ![]()
| LITERATURE CITED |
|---|
|
|
|---|
Armbruster W. S. 1992 Phylogeny and the evolution of plantanimal interactions. BioScience 42: 12-20
Armbruster W. S. K. E. Schwaegerle 1996 Causes of covariation of phenotypic traits among populations. Journal of Evolutionary Biology 9: 261-276[CrossRef][ISI]
Baker H. G. 1967 The evolution of weedy taxa in the Eupatorium microstemon species aggregate. Taxon 16: 293-300[CrossRef]
Baldwin B. G. B. L. Wessa 2000 Origin and relationships of the tarweedsilversword lineage (CompositaeMadiinae). American Journal of Botany 87: 1890-1908
Barrett S. C. H. L. D. Harder A. C. Worley 1996 The comparative biology of pollination and mating in flowering plants. Philosophical Transactions of the Royal Society of London, Series B 351: 1271-1280[CrossRef]
Charlesworth D. S. Mayer 1995 Genetic variability of plant characters in the partial inbreeder Collinsia heterophylla (Scrophulariaceae). American Journal of Botany 82: 112-120[CrossRef][ISI]
Cruden R. W. 1973 Reproductive biology of weedy and cultivated Mirabilis (Nyctaginaceae). American Journal of Botany 60: 802-809[CrossRef][ISI]
Cruden R. W. 1977 Pollen-ovule ratios: a conservative indicator of breeding systems in flowering plants. Evolution 31: 32-46
Cruden R. W. 2000 Pollen grains: why so many?. Plant Systematics and Evolution 222: 143-165[CrossRef][ISI]
Donoghue M. J. 1989 Phylogenies and the analysis of evolutionary sequences, with examples from seed plants. Evolution 43: 1137-1156[CrossRef][ISI]
Doyle J. J. J. L. Doyle 1987 A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemical Bulletin 19: 11-15
Endler J. 1976 Natural selection in the wild. Princeton University Press, Princeton, New Jersey, USA
Endler J. 1995 Multiple-trait coevolution and environmental gradients in guppies. Trends in Ecology and Evolution 10: 22-29
Eriksson T. 1998 AutoDecay 4.0, version 4.0 (program distributed by the author). Department of Botany, Stockholm University, Stockholm, Sweden
Fenster C. B. P. K. Diggle S. C. H. Barrett K. Ritland 1995 The genetics of floral development differentiating two species of Mimulus (Scrophulariaceae). Heredity 74: 258-266[ISI]
Fishman L. R. Wyatt 1999 Pollinator-mediated competition, reproductive character displacement, and the evolution of selfing in Arenaria uniflora (Caryophyllaceae). Evolution 53: 1723-1733[CrossRef][ISI]
Galen C. 1999 Why do flowers vary? The functional ecology of variation in flower size and form within natural plant populations. BioScience 49: 631-640[CrossRef][ISI]
Galen C. 2000 High and dry: drought stress, sex-allocation trade-offs, and selection on flower size in the alpine wildflower Polemonium viscosum (Polemoniaceae). American Naturalist 156: 72-83[CrossRef][Medline]
Galen C. R. A. Sherry A. B. Carroll 1999 Are flowers physiological sinks or faucets? Costs and correlates of water use by flowers of Polemonium viscosum. Oecologia 118: 461-470[CrossRef][ISI]
Ganders F. R. G. R. Krause 1986 Systematics of Collinsia parviflora and C. grandiflora (Scrophulariaceae). Madroño 33: 63-70
Garber E. D. 1956 The genus Collinsia. I. Chromosome number and chiasma frequency of species in two sections. Botanical Gazette 118: 71-73[CrossRef]
Garber E. D. 1958a The genus Collinsia. III. The significance of chiasma frequencies as a cytotaxonomic tool. Madroño 14: 71-72
Garber E. D. 1958b The genus Collinsia. VII. Additional chromosome numbers and chiasma frequencies. Botanical Gazette 120: 55-56[CrossRef]
Garber E. D. 1975 Collinsia. In R. C. King [ed.], Handbook of genetics, vol. 2, 333361. Plenum Press, New York, New York, USA
Grant V. 1958 The regulation of recombination in plants. Cold Springs Harbor Symposium on Quantitative Biology 23: 337-363
Harvey P. H. M. D. Pagel 1991 The comparative method in evolutionary biology. Oxford University Press, Oxford, UK
Jain S. K. 1976 The evolution of inbreeding in plants. Annual Review of Ecololgy and Systematics 7: 469-495
Jonas C. S. M. A. Geber 1999 Variation among populations of Clarkia unguiculata (Onagraceae) along altitudinal and latitudinal gradients. American Journal of Botany 86: 333-343
Kalisz S. D. Vogler B. Fails M. Finer E. Shepard T. Herman R. Gonzales 1999 The mechanism of delayed selfing in Collinsia verna (Scrophulariaceae). American Journal of Botany 86: 1239-1247
Kearns C. A. D. W. Inouye 1993 Techniques for pollination biologists. University Press of Colorado, Niwot, Colorado, USA
Kingston B. L. 1998 Papillate tepal hoods and delayed pollen germination in Albuca L. (Liliaceae). Annals of Botany 82: 263-266
Lande R. D. W. Schemske 1985 The evolution of self-fertilization and inbreeding depression in plants. I. Genetic models. Evolution 39: 24-40[CrossRef][ISI]
Li C. C. 1975 Path analysisa primer. Boxwood Press, Pacific Grove, California, USA
Lloyd D. G. 1965 Evolution of self-compatibility and racial differentiation in Leavenworthia (Cruciferae). Contributions of the Gray Herbarium, Harvard 195: 3-133
Mayer S. S. D. Charlesworth B. Meyers 1996 Inbreeding depression in four populations of Collinsia heterophylla Nutt (Scrophulariaceae). Evolution 50: 879-891[CrossRef][ISI]
Motten A. F. J. L. Stone 2000 Heritability of stigma position and the effect of stigma-anther separation on outcrossing in a predominantly self-fertilizing weed, Datura stramonium (Solanaceae). American Journal of Botany 87: 339-347