Am. J. Bot. Join the BSA
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 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 (37)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Jonas, C. S.
Right arrow Articles by Geber, M. A.
Right arrow Search for Related Content
PubMed
Right arrow PubMed Citation
Right arrow Articles by Jonas, C. S.
Right arrow Articles by Geber, M. A.
Agricola
Right arrow Articles by Jonas, C. S.
Right arrow Articles by Geber, M. A.
(American Journal of Botany. 1999;86:333-343.)
© 1999 Botanical Society of America, Inc.

Variation among populations of Clarkia Unguiculata (Onagraceae) along altitudinal and latitudinal gradients1

Christine S. Jonas and Monica A. Geber2

Section of Ecology and Systematics, Corson Hall, Cornell University, Ithaca, New York 14853

Received for publication September 10, 1996. Accepted for publication June 25, 1998.


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
We investigated phenotypic variation in 15 traits in greenhouse-grown plants from 16 populations of Clarkia unguiculata from three elevational habitats and six latitudinal transects. Populations from the lowest and highest elevations were geographically and ecologically marginal within the species' range. We (1) describe patterns of trait variation with elevation and latitude; (2) compare latitutidinal variation between marginal and central areas of the species' range; and (3) compare patterns of variation within C. unguiculata to interspecific patterns within the genus. Although there was some evidence that traits varied clinally (i.e., increased/decreased monotonically) along environmental gradients, interaction effects between altitude and latitude dominated patterns of variation. For most traits, latitudinal trends at the low-elevation margin of the species' range differed from trends at mid- and high-elevation areas. Based on interspecific comparisons, populations at the hotter, more arid ends of both environmental gradients were expected to have rapid development, small flowers and vegetative size, low levels of herkogamy and protandry, and high rates of gas exchange. Instead, we found that while some traits were correlated with one gradient in the expected way (e.g., development time with elevation, gas-exchange physiology with latitude), all traits were not consistently associated with each other along both gradients, and intraspecific patterns of variation differed from interspecific patterns.

Key Words: altitude • Clarkia unguiculata • elevation • floral traits • geographic variation • marginal and central populations • morphology • Onagraceae • phenology • physiology


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
The study of natural variation in ecologically important traits has a long history in evolutionary biology (Turesson, 1922 ; Ford, 1954 ; Clausen, Keck, and Hiesey, 1948 ; see Endler, 1977 , 1986 ; Linhart and Grant, 1996 ). The primary motivation for studies of this nature is to further our understanding of the roles of natural selection, gene flow, and genetic drift in shaping levels of adaptation and in affecting the distribution and abundance of species.

In plants, studies have shown that populations are often genetically differentiated with respect to life history, morphology, and physiology, at both small and large spatial scales (Turesson, 1922 ; Clausen, Keck, and Hiesey, 1948 ; Hiesey and Milner, 1965 ; Langlet, 1971 ; Briggs and Walters, 1984 ; Linhart and Grant, 1996 ). Studies of variation at geographic scales have focused on two broad areas of inquiry. First, what are the patterns and causes of variation in ecologically important traits in relation to environmental gradients, such as latitude or elevation? Second, what are the patterns and causes of trait variation between populations at the center vs. the margin of a species' geographic range?

There is extensive documentation of intraspecific variation in ecologically important traits along environmental gradients (see references in Stebbins, 1950 ; Endler, 1977 , 1986 ; Briggs and Walters, 1984 ). In many cases, patterns of trait variation across environmental gradients have suggested hypotheses concerning their adaptive significance. Adaptive explanations become more convincing (a) when the same pattern of trait variation in relation to an environmental gradient is found in different parts of a species' range, (b) when functionally related traits vary in coordinated fashion across an environmental gradient, or (c) when intra- and interspecific patterns of trait variation in relation to an environmental gradient are similar. It is rare, however, that researchers document character variation along replicate transects paralleling the same environmental gradient throughout a species' range, as most studies are confined to variation along a single transect (but see Reinartz, 1984a , b ; Thorpe and Baez, 1987 ; Nevo, Carver, and Beiles, 1991 ; Endler, 1995 ). Furthermore, many studies of character variation within species along environmental gradients focus on one or a few traits (e.g., Lacey, 1988 ; Winn and Gross, 1993 ; but see Reinartz, 1984a , b ). Lastly, it is not often that intraspecific variation can be compared to patterns of interspecific variation.

Geographic range margins are often thought to represent ecologically marginal environments for a species (Mayr, 1963 ; Lewontin, 1974 ; Antonovics, 1976 ; Lewis and Lewis, 1955 ). As such, they provide an important arena for empirical and theoretical studies on the limits to adaptation (Mayr, 1963 ; Antonovics, 1976 ; Hoffmann and Parsons, 1991 ; Geber and Dawson, 1993 ; Lewis and Lewis, 1955 ). Thus, researchers have attempted to identify the environmental factors that demarcate range margins (Connell, 1961 ; Gill and Hawksworth, 1961 ; Woodward and Williams, 1987 ; Root, 1993 ), and to pinpoint the traits that may be important in limiting species' distributions at margins (St. Clair and Gregory, 1990 ; Shuter and Post, 1990 ; Loik and Nobel, 1993 ; see Lewis and Lewis, 1955 ). In addition, comparisons have been made of levels of genetic variation between marginal and central populations because marginal populations are predicted to be depauperate in variation as a result of genetic drift, inbreeding, or strong selection (Lewontin, 1974 ; Brussard, 1984 ; Ohlson, 1989 ; Hoffmann and Parsons, 1991 ; Rice and Mack, 1991a ; Blows and Hoffmann, 1993 ; Johansson, 1994 ; Stewart and Nilsen, 1995 ). Theoreticians have also modeled the interplay between selection and gene flow in shaping adaptation to ecological and/or geographical margins (Garcia-Ramos and Kirkpatrick, 1997 ; Kirkpatrick and Barton, 1997 ), and empirical work has suggested that gene flow between ecologically favorable and unfavorable habitats may swamp adaptation to marginal environments (Rice and Mack, 1991b ; Stanton and Galen, 1997 ). Lastly, marginal populations have been considered for their value in conservation efforts because they may harbor unusual phenotypes and genotypes (Millar and Libby, 1991 ; Lesica and Allendorf, 1994 ).

In this study, we characterized geographic variation in ecologically important traits across latitudes and between marginal and central populations in the annual plant Clarkia unguiculata Lindl. (Onagraceae), as part of a broader investigation of patterns and causes of the limits to adaptation in this species. Clarkia unguiculata is a primarily outcrossing endemic to California. It is highly variable in morphology and phenology, but typically germinates with the onset of winter rains and flowers in early summer (Lewis and Lewis, 1955 ).

The species has a U-shaped geographic distribution in the mountains that ring the San Joaquin Valley, from Mendocino, Butte, and Plumas Counties in the North, to Los Angeles County in the South (Fig. 1). It occurs along the coastal mountains to the west, and the Sierra Nevada to the east, but is not found in the San Joaquin Valley proper. Several disjunct populations also occur in San Diego County (Lewis and Lewis, 1955 ). Clarkia unguiculata typically grows in fairly discrete populations, on slopes in oak or oak–Gray pine woodland, and occasionally in chaparral. At low-elevation arid regions, populations occur in grassland and oak savanna habitat with sparse or no tree cover (Lewis and Lewis, 1955 ; Jonas, 1997 ).



View larger version (48K):
[in this window]
[in a new window]
 
Fig. 1. Map of California, showing the U-shaped geographic range of Clarkia unguiculata (light shading) and area of this study (dark shading). The locations of the elevational transects at six latitudes are shown. In the four southernmost transects, one population was sampled from low-, mid-, and high-elevation habitats (x). At the two northernmost transects, only populations from mid- and high-elevation habitats were sampled. In all transects the low- and high-elevation populations were located at the margin of the species range and the mid-elevation population was located in the central region of highest abundance.

 
In species with a U-shaped distribution, there are not only high- and low-latitudinal range margins, but there are also interior and exterior margins along the inner and outer edges of the U. For C. unguiculata, the interior margin in the Sierra Nevada coincides with low-elevation foothills at the edge the San Joaquin Valley, while the outer margin corresponds to a high-elevation limit on the western slopes of the mountains. Clarkia unguiculata reaches its highest abundance, both in terms of population numbers and size, at mid-elevations in the Sierra Nevada; at low- and high-elevation margins, populations are sparser and smaller (Jonas, 1997 ). The pattern of declining abundance from the "center" toward the periphery of the range is typical of many species (Brown, 1984 ; Lawton, 1993 ).

On the western slopes of the Sierra Nevada, precipitation increases, while temperature and potential evapotranspiration decrease, with elevation and latitude. The elevational gradient in climate has been of particular interest in studies of character evolution within the section Phaeostoma of the genus Clarkia, of which C. unguiculata is a member. Several species in Phaeostoma occur in or beyond the low-elevation margins of C. unguiculata habitat, where the growing season is shortened by the early onset of summer drought (Lewis and Lewis, 1955 ; Vasek, 1964 , 1968 ). Comparisons among species of Phaeostoma have suggested that low-elevation arid habitats are associated with earlier reproduction, smaller vegetative and floral size, and increased rates of self-fertilization due to reduced herkogamy and protandry (Vasek, 1964 ; Vasek and Sauer, 1971 ). In addition, common garden studies of two C. unguiculata populations have suggested that, within the species, populations from lower elevations flower earlier (Vasek, 1964 , 1971 ). Based on these associations, early flowering has been interpreted as an adaptation to avoid the early onset of summer drought at low elevations, both among species and within C. unguiculata. However, the extent to which the differences found between the two populations of C. unguiculata represent a consistent pattern of clinal variation within the species is unknown. It is also unclear whether the other traits that are associated with hotter, drier environments interspecifically (small vegetative and floral size, reduced herkogamy, and protandry) are associated in this way with C. unguiculata. Lastly, it is not known whether north–south gradients in latitude and high–low gradients in elevation result in similar patterns of character variation and association among populations of C. unguiculata.

We characterized variation in phenology, floral and vegetative morphology, and gas-exchange physiology among populations of C. unguiculata from three elevational habitats at six latitudes representing most of the species' range in the Sierra Nevada (Fig. 1). This design allowed us to compare patterns of trait variation along replicate environmental gradients. Using "common garden" experiments, we examined population differentiation in the context of (1) describing patterns of trait variation with elevation and latitude; (2) comparing patterns of trait variation across latitudes between marginal and central populations; and (3) comparing patterns of trait variation along environmental gradients within C. unguiculata to interspecific patterns of variation.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Study populations
During the summer of 1995, we collected seed from populations along six elevational transects in the Sierra Nevada spanning a latitudinal range from Kern County in the South, to Plumas County in the North (Fig. 1, Table 1). At the four southernmost transects (latitudes 1–4), we collected seed from one population in each of three habitat types: (1) low-elevation oak savanna sites had no or widely spaced tree cover (usually Quercus douglasii) and an understory dominated by grasses; (2) mid-elevation oak woodland sites were characterized by relatively open cover of Q. douglasii and Q. wislizenii and some scattered California buckeye (Aesculus californica) and Gray pine (Pinus sabiniana); (3) high-elevation oak woodland sites were characterized by a similar species assemblage, but the overstory contained a higher proportion of evergreen species (for a description of elevational gradients in California oak woodland, see Griffin [1977]). At all four latitudes, the low- and high-elevation populations were from the lowest and highest elevations, respectively, at which C. unguiculata could be located after an extensive search. They therefore represent the low- and high-elevation margins (both geographic and ecological) of the species' range in the Sierra Nevada. Mid-elevation sites were chosen in the central area of highest abundance. Clarkia unguiculata did not occur in the low-elevation habitat type at the two northernmost transects (latitudes 5–6). Consequently, only the mid- and high-habitat types are represented in these two transects. The distance between populations within an elevational transect was at least 5 km.


View this table:
[in this window]
[in a new window]
 
Table 1. Latitude and elevation of the 16 study populations of Clarkia unguiculata.

 
Study populations contained >300 individuals each, and collections were made when >95% of the plants in a given population had produced at least one seed capsule (dates of 95% maturity ranged from late June to early August). In each upright capsule, a few to hundreds of seeds may develop; in the field, the seeds disperse passively as the capsule dehisces in mid-late summer. We collected at least one mature but unopened capsule from the main stem of 100 individuals in each population. Individuals were 1 m apart on a transect placed through the middle of the population. If one transect did not yield capsules from 100 plants, we sampled plants from a second transect perpendicular to the first. Seeds were stored dry in coin envelopes for 2 mo until the beginning of the greenhouse study.

Climate information
There is a predictable gradient of temperature, precipitation, and potential evapotranspiration (potE) along the elevational cline through C. unguiculata's range in the Sierra Nevada. Major (1977) reported climes in long-term averages of these three abiotic parameters along West-to-East transects at five latitudes in the southern Sierra from Kern Canyon (35°30'N) in the South, to the hills above Stockton (38°N) in the North. At each latitude, temperature decreased with increasing elevation at average rates of -0.35 to -0.65° C/100 m. Precipitation increased at a rate of +94 to +208 mm/100 m, and potE decreased at -13 to -36 mm/100 m. Ten-year (1986–1995) monthly averages of precipitation from weather stations positioned along the elevational gradient at latitude 3 of this study support this trend (Fig. 2, data from California Department of Water Resources, Division of Flood Management). Vasek and Sauer (1971) also found that, at lowland C. unguiculata sites in Kern County, soil drying occurred more rapidly and earlier in the season, and average temperatures were higher than at higher elevation populations.



View larger version (22K):
[in this window]
[in a new window]
 
Figs. 2–3. Monthly precipitation averaged over the years 1986–1995 (data from California Department of Water Resources). 2. At high, mid, and low elevations along a transect at Latitude 3 of this study. 3. At five stations near high-elevation populations at five latitudes.

 
Latitudinal gradients in climatic factors in this region of California are less well documented. However, Major (1977) presented evidence of increasing annual precipitation and decreasing temperatures with latitude at grassland elevations and at 500 and 1000 m. Data from precipitation stations ~200 m below/above our high-elevation sites at latitudes 1, 2, and 4–6 suggest that the trend is similar at these elevations (Fig. 3).

Common garden study
In October 1995, 24 plants were randomly selected from the 100 seed plants sampled at each study population for use in the common garden experiment. Five seeds from each plant were planted together directly into tube-shaped pots (655.48 cm3; DeepotsTM, Stuewe & Sons, Inc., Corvallis, Oregon) filled with fritted clay (for a description of the properties of fritted clay, see van Bavel, Lascano, and Wilson [1978]). A small plug of Metro-Mix soil was placed at the top of each pot to prevent seeds from falling down between clay particles. Pots were arranged in a randomized block design, with two pots from each population in each of 12 blocks. Each pot was 6.35 cm from its nearest neighbor. This spacing was maintained for the duration of the study.

At the outset, plants were kept in a controlled environment chamber with diurnally ramped temperatures (10°C night, 22°C day) and 10 h of light provided by two 1000-W high-pressure sodium bulbs, and were watered evenly with a mist nozzle two times daily. Pots were observed daily for germination. Ninety-nine percent of pots had at least one emergent individual. After 8 d in the chamber, plants were moved into the greenhouse, where they were grown under natural light for 6 d, after which time each block was positioned directly below a 400-W high-pressure sodium bulb. Plants were fertilized once per wk with a 125-ppm solution of EXCEL (15-5-15, plus micronutrients) beginning 10 d after planting. Greenhouse temperatures ranged between 18°C and 32°C. Plants were watered at least once per d, with extra watering when needed. At 17 d after planting, seedlings were thinned to one per pot, leaving the tallest individual in all cases. The germination date for the remaining individual was assumed to be the date of first germination for its pot. During the course of the experiment, 56 plants developed striking morphological abnormalities due to an apparent viral infection. These plants were not included in the analysis.

Each individual was scored for 15 traits (Table 2). Flowering time was scored as the number of days from germination to petal opening of the first flower. All flower measurements were made on the first flower using electronic calipers (Mitutoyo Digimatic, Japan). Clarkia unguiculata has two sets of stamens. Measurements were done on the larger set, which has anthers that dehisce later (closer to the time of stigma receptivity). Stamen length was measured on the first day of flower opening, and petal and style lengths were measured on the first day of stigma receptivity. Herkogamy was calculated as the difference between style and stamen lengths. Protandry was scored as the number of days from first flower opening (when the pair of short anthers begin to dehisce) to stigma receptivity.


View this table:
[in this window]
[in a new window]
 
Table 2. Traits scored on Clarkia unguiculata individuals grown in a common garden.

 
The phyllotaxis of C. unguiculata is generally alternate, but early leaves are oppositely arranged. Data from a subset of individuals indicated that there were no consistent differences among populations in the node at which the switch to alternate leaf arrangement occurs. Therefore, for traits based on node number, each node was counted once in this study, regardless of leaf arrangement at that node. Average internode length was calculated by dividing height to first flowering node by the number of nodes to first flower. Total number of nodes and plant height were measured on the first day of stigma receptivity of the first flower. Height was measured from the base of the plant to the apical meristem of the main stem. Leaf area was measured using a LI-COR 3000 leaf area meter (LI-COR, Lincoln, Nebraska). Plants were grown for 4 mo, at which time flowering had finished and plants had begun to senesce.

Instantaneous rates of gas exchange were measured on a leaf area basis using a CIRAS-1 portable open gas-exchange system at ambient levels of CO2 (PP Systems, Herts, UK). Measurements were done on one block per day between 1000 and 1400, as this was determined to be the period of greatest gas-exchange activity based on several days of diurnal measurements (C. Jonas, Cornell University, unpublished data). All measurements were made on the newest fully expanded leaf of each individual, just before most plants began flowering. Leaves were allowed to acclimate to conditions within the cuvette until the photosynthesis curve reached a plateau. Light levels were 1500–1900 µol·m-2·s-1 photosynthetic photon flux density (PPFD) (the range was <100 µol·m-2·s-1 PPFD within blocks), provided by a supplemental source within the CIRAS-1 cuvette. Leaf temperatures ranged from 22.7° to 28.9° during measurements. Calculations of photosynthetic rate (Amax) and stomatal conductance (g) were based on equations of von Caemmerer and Farquhar (1981) . Measurements were made only on populations from the four southern transects, where all three elevational habitats were represented.

Differences among plants grown in a common environment reflect some combination of effects due to the genotype of the study individual and the genotype and growth environment of the maternal plant. In studies like this one, where the main goal is to examine genetically-based differences among populations, some care must be taken to factor out maternal effects. Seed size is the best documented route by which maternal genotype and environment can affect offspring traits in plants (Roach and Wulff, 1987 ). Maternal effects have been shown to influence germination, survivorship, seedling biomass, adult plant size, and competitive ability through effects on seed size (see references in Schaal, 1984 ; Roach and Wulff, 1987 ). To control as much as possible for maternal effects on the traits of interest, maternal seed mass was used as a covariate in all analyses. Five seeds from each maternal capsule were weighed on a electronic microbalance (XM1000 P, Sartorius Instruments, McGaw Park, Illinois) for an estimate of average seed mass.

Data analysis
Traits were analyzed using a mixed-model ANCOVA (JMP, SAS Institute, Inc., Cary, North Carolina) of the form: trait = seed mass + latitude + habitat + block + interaction terms + error. Block was a random effect. Three- and four-way and nonsignificant two-way interaction terms (at {alpha} = 0.1) were dropped from the model for the final analyses. Several traits were analyzed using transformed data to correct for heterogeneity of variances in the original data set (residuals were tested using Levene's test for homogeneity of variances). Days to germination, protandry, and node of first flower were log transformed, largest leaf size was square-root-transformed, and stomatal conductance (g) was squared for the final analysis. Homogeneity of variances was improved in all cases, however, the hypothesis of equal variances was rejected by Levene's test for transformed values of days to germination, leaf size, and g. Normal probability plots of the residuals from the ANCOVA analyses for all traits indicated normal distributions.

Main effects and interactions were tested for significance by a sequential Bonferroni procedure at an overall {alpha} = 0.05 (Rice, 1989 ). For each main effect or interaction effect, traits were ranked in order of P values obtained from the ANCOVA analysis, from smallest (P1) to largest (Pk). For each trait of rank i, the significance of the effect was tested at an {alpha} level of 0.05/(1 + k - i).

For traits showing a significant latitude by habitat effect, we used a series of contrasts designed to partition the interaction effect SS into parts due (1) to differences between mid- and high-elevation populations across latitudes, and (2) to differences between low- vs. mid- and high-elevation populations across latitudes (C. McCulloch, Cornell University, personal communication). For these contrasts, we used data from the four southernmost transects only, where all three elevations were represented. F values obtained for each trait were tested for significance by a sequential Bonferroni procedure at a "tablewide" level of {alpha} = 0.05 within each contrast (Rice, 1989 ).

To examine whether traits exhibited clinal variation (i.e., monotonic increasing or decreasing) across elevational habitats, we computed the Spearman rank correlation between the population mean values for each trait and habitat (ranked 1–3 from low- to high-elevation habitats). All 16 populations were included in the calculation of the correlation between a trait and elevational habitat, because, in spite of the often significant latitude by habitat interaction effect on a trait (see Results), there were simply too few populations to calculate a separate correlation at each latitude. To examine whether traits exhibited clinal variation across latitudes, we computed the Spearman rank correlation between the population mean values for each trait and latitude, including all 16 populations, and then separately for populations from low-, mid-, and high-elevational habitats. This method of examining clinal variation does not adjust for correlation between traits among populations. Traits that are correlated among populations will show similar patterns of variation across elevation or latitude.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Effects of seed mass on trait variation
Seed mass, which is a potential source of maternal influence on trait variation, had no significant effect on any life history, floral, vegetative, or physiological trait, based on the sequential Bonferroni analysis. Because seed mass was a covariate in all analyses, differences among populations in trait values most likely represent true genetic differentiation.

Patterns of trait variation
There were significant differences among elevational habitats in days to germination, days to first flower, and Amax and among latitudes in all traits (Table 3). Thus, traits varied more widely across latitudes than across elevations, as might be expected from the greater variation in precipitation across latitudes than across elevations (Figs. 2, 3), and from the greater distance between elevational transects.


View this table:
[in this window]
[in a new window]
 
Table 3. Analyses of covariance of phenological, morphological, and physiological traits in Clarkia unguiculata. Latitude refers to the north-south positions of the six elevational transects; habitat refers to the low-, mid-, and high-elevation locations of populations within transects. Nonsignificant interaction terms were dropped for the final analyses. Contrasts of latitudinal variation among elevations were restricted to populations from the four southernmost latitudes (see text).

 
Spearman rank correlations indicated that the following traits showed clinal (i.e., monotonically increasing/decreasing) variation with elevation across all populations combined: days to germination and flowering, and node of first flower increased with elevation, while there was a tendency for petal size to decrease with elevation (Table 4). Clinal variation with latitude across all populations combined was found for petal, style, and stamen lengths, node of first flower, total number of nodes, Amax, and g (Table 4, latitude, all populations). The morphological characters increased in value with latitude, while physiological traits decreased with latitude. These clinal patterns were often not statistically significant when correlations were calculated from the smaller number of populations at a single elevation, but they were of the same direction (see Table 4). For a few traits, significant clinal variation across latitude was found at one or more elevations but not for all populations combined (e.g., petal width, internode length).


View this table:
[in this window]
[in a new window]
 
Table 4. Spearman rank correlations between populations mean values for each trait and (a) latitude and (b) habitat (= elevation). Separate correlations were calculated for all populations (N = 16), and for populations at high- (N = 6), mid- (N = 6), and low-elevation (N = 4) habitats across latitudes. Correlations in bold type are statistically significant at P < 0.05; underlined correlations are significant at P < 0.1.

 
In general, traits showing clinal variation with latitude differed from those showing clinal variation with elevation, except for node of first flower. In the case of flower size, clinal patterns across latitude appeared to be in the opposite direction to clinal patterns across elevations. Latitudinal and elevational trends appeared to be quite different for traits related to the propensity for selfing (style, stamen lengths, herkogamy, and protandry), and for gas-exchange physiology.

For all traits except number of nodes, Amax, and g, there was a highly significant latitude by habitat interaction, indicating that the pattern of trait variation across latitudes differed among low-, mid-, and high-elevation populations. Contrasts of trait variation between mid- and high-elevation populations across latitude were significant for days to germination, days to first flower, petal width, style length, internode length, and plant height (i.e., 6 out of 12 traits; Table 3). Contrasts between low vs. mid and high transects were significant for all traits except days to germination and node of first flower (i.e., 10 out of 12 traits; Table 3). The first set of contrasts accounted for a lower fraction of the interaction variance than the second set, as judged by the coefficients of partial determination for the contrasts (Table 3). Thus, mid- and high-elevation populations varied more nearly in parallel with each other across latitudes than they did with the low-elevation populations (Figs. 4–18). The closer parallelism in trait variation between mid- and high-elevation populations was most obvious for leaf size, style length, stamen length, and herkogamy, all of which decreased at mid-latitudes (latitude 3 and/or 4) relative to neighboring latitudes (2–3, 5) at both elevations (Figs. 4, 8–10, 16).



View larger version (30K):
[in this window]
[in a new window]
 
Figs. 4–18. Least squares mean character values for high-, mid-, and low-elevation populations at six latitudinal transects. 4. Days to germination (d). 5. Days from germination to opening of first flower (d). 6. Petal length (mm). 7. Petal width (mm). 8. Style length (mm). 9. Stamen length (mm). 10. Herkogamy (mm). 11. Protandry (d). 12. Node of first flower. 13. Average internode length (mm). 14. Total number of nodes. 15. Plant height (cm). 16. Area of largest leaf (cm2). 17. Amax (µmol CO2·m-2·s-1). 18. g (mmol H2O·m-2·s-1). (Values reported for Figs. 11, 12, 16, and 18 are back-transformed least squares means.)

 

    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Environmental gradients and intra- vs. interspecific comparisons
Elevational gradients are thought to have played a large role in interspecific differentiation within the genus Clarkia. In the section Phaeostoma, early onset of flowering has been interpreted as an adaptive mechanism to avoid the heat and desiccation stress of early summer drought at low elevations. Studies of reproductive timing in these annuals have focused on comparisons between the protandrous, primarily outcrossing C. unguiculata and several closely related species that tend to occur in or beyond the low-elevation margins of C. unguiculata habitat. These other species, C. exlis, C. springvillensis, and C. tembloriensis, differ from C. unguiculata in that they flower earlier, have smaller leaves and floral parts, lack pubescence on the ovary and calyx, and have less separation between anthers and receptive stigma in time (protandry) and in space (herkogamy) (Lewis and Lewis, 1955 ; Vasek, 1964 , 1968 ). Within Phaeostoma, protandry and herkogamy are highly correlated with outcrossing rate (Moore and Lewis, 1965 ; Holtsford and Ellstrand, 1992 ). Thus, among species in the section, there is a suite of traits related to life history, morphology, and breeding system that is associated with the short season of low-elevation habitats.

The association of a self-fertilizing breeding system and early flowering with drier habitats is seen in several other genera as well and has led to the hypothesis that self-fertilization could evolve from outcrossing as a result of selection for rapid development in drought-stressed environments (Guerrant, 1984 ; Diggle, 1992 ; Eckhart, Geber, and Jonas, 1996 ). In this scenario, selection for rapid reproduction leads to truncated morphological development resulting in small flowers and, because of allometric relationships, reduced stigma–anther separation. To the extent that plants achieve rapid development by faster growth and an increase in metabolic rate, physiological rates such as photosynthesis and stomatal conductance are also expected to be higher in plants selected for early flowering (Geber and Dawson, 1990 ; Ehleringer, 1993 ; Jonas, 1997 ).

Previous studies of life history and morphology based on a smaller region of C. unguiculata's range have suggested that population differentiation within the species follows the elevational pattern seen among species in Phaeostoma. In the Caliente Hills of Kern County, a low-elevation population (305 m) flowered earlier than a high-elevation population (790 m) both in the field and in a greenhouse (Vasek, 1971 , 1977 ; Vasek and Sauer, 1971 ), and had reduced protandry and a nonsignificant trend toward smaller reproductive parts in the greenhouse (Vasek, 1971 , 1977 ). In a companion to the present study, leaf carbon isotope analysis of greenhouse-grown C. unguiculata from populations at latitudes 1 and 2 suggested that rates of physiology were negatively correlated with development time and with elevation (Jonas, 1997 ).

By contrast, the results of this study indicate that early flowering, small vegetative and floral morphology, selfing, and high physiological rates are not always associated in a clinal way in C. unguiculata. Low-elevation populations do tend to germinate and flower earlier and to flower at an earlier node than higher elevation populations (Figs. 4, 5, and 12; Table 4), but the elevational pattern is not repeated at all latitudes (as indicated by significant latitude by habitat interaction effects, Table 3). Thus, though there is evidence that some populations from the hottest, most arid habitat have faster development, this trend does not hold in all latitudinal regions of the species' range. There is also no evidence of smaller floral or vegetative size at lower elevations. In fact, petals tend to be larger at lower elevations (see Table 4, Figs. 6–7). Finally, there is no indication that low-elevation populations have reduced herkogamy or protandry, or higher rates of photosynthesis and stomatal conductance (Table 4, Figs. 10–11, 17–18). If interspecific patterns are indeed a product of selection for early summer drought avoidance, the results of our study suggest that the same evolutionary processes have not acted inter- and intraspecifically to shape this suite of traits.

Based on interspecific comparisons along gradients of temperature and aridity, we expected latitude to be negatively correlated with rates of development and photosynthesis and positively correlated with floral and vegetative size in C. unguiculata. The evidence for such correlations is equivocal. Spearman rank correlations indicate that node of first flower and flower size (petal, style, and stamen length) decrease, and that gas-exchange rates increase at lower latitudes for all populations combined (Table 4). There are also nonsignificant trends of earlier germination and flowering and reduced herkogamy and protandry at lower latitudes. However, these associations are not always reflected among populations at a single elevation.

In sum, this study found no strong evidence for a consistent association among rapid development, small floral and vegetative size, reduced herkogamy and protandry, and increased physiological rates at the hotter, more arid ends of latitudinal or altitudinal gradients. Furthermore, while latitudinal and elevational gradients present similar clines in temperature and precipitation, these are not reflected in congruent patterns of character variation in C. unguiculata. Lastly, intra- and interspecific patterns of trait variation differ across elevations.

Central and marginal populations and latitude by habitat interaction
We found it visually striking that trait means from mid- and high-elevation populations seemed to vary more nearly in parallel across latitudes compared to low-elevation populations (Figs. 4–16). Contrasts designed to test this observation indicated that, although mid- and high-elevation populations differ significantly in latitudinal patterns of variation for many traits, there is a greater difference between low-elevation populations and populations from the two higher elevations for all traits except days to germination (Table 3).

High- and mid-elevation populations share several unusual patterns of variation that account for part of their statistical divergence from low-elevation populations. At latitude 4, high- and mid-elevation populations have unexpectedly small styles, stamens, and herkogamy, together with unusually high protandry (Figs. 8–11). This pattern of trait association is surprising in that herkogamy and protandry are otherwise positively correlated within and among other populations of C. unguiculata, and among other members of the genus (Moore and Lewis, 1965 ; Holtsford and Ellstrand, 1992 ; C. Jonas and M. Geber, Cornell University, unpublished data). These populations also share several other floral traits that are unique to this area of C. unguiculata's range. Petals are darker red than in other regions, buds are globose, and sepals separate completely when they reflex at flower opening (C. Jonas, Cornell University, personal observation). At latitude 2, high- and mid-elevation populations have unexpectedly large leaves (Fig. 16). These were the only two sites that had a high frequency of fog cover during the summer of 1995 (C. Jonas, Cornell University, personal observation).

The greater overall divergence in trait values between low- vs. mid- and high-elevation populations may be related to greater divergence in habitat. The low-elevation region of C. unguiculata's range is characterized by a sparse or absent overstory and a grass-dominated understory, in contrast to well-developed woody overstories and similar understory species compositions in the upper and central regions. Thus, latitudinal environmental variation may be qualitatively different and lead to different patterns of selection on ecologically important traits between low vs. central and high regions of the species' range. At present, however, we cannot say whether differences in selection history or in levels of gene flow or genetic drift account for the differentiation of low-elevation populations.

Peripheral populations of widespread species are often targeted for conservation effort with the intention of preserving the among-population component of genetic diversity (Lesica and Allendorf, 1994 ). Despite similar degrees of isolation and population size (C. Jonas, Cornell University, personal observation) and similar elevational distance from the center of C. unguiculata's range, populations from high and low margins showed different amounts of similarity to central populations in patterns of trait variation. To the extent that phenotypic divergence reflects genetic divergence, all peripheral populations of C. unguiculata may not be equally valuable to the conservation of diversity within species, and divergence among marginal and central populations may correspond to habitat differences.


View this table:
[in this window]
[in a new window]
 
Table 3. Extended.

 


View larger version (34K):
[in this window]
[in a new window]
 
Figs. 4–18. Continued.

 


View larger version (13K):
[in this window]
[in a new window]
 
Figs. 4–18. Continued.

 

    FOOTNOTES
 
1 The authors thank Timothy Holtsford for help locating study populations; Todd Dawson, Roman Pausch, and Ed Grote for help with gas-exchange equipment; Sarah Yonder for help with measurements; Chuck McCulloch for invaluable statistical advice; and Sean Thomas and an anonymous reviewer for useful comments on the manuscript. This work was supported by a NSF graduate fellowship to C. S. J., and NSF (DEB-) and Andrew Mellon grants to M. A. G. Back

2 Author for correspondence. Back


    LITERATURE CITED
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS
 DISCUSSION
 LITERATURE CITED
 
Antonovics, J. 1976 The nature of limits to natural selection. Annals of the Missouri Botanical Garden 63: 224–247.[CrossRef][ISI]

Blows, M. W., and A. A. Hoffmann. 1993 The genetics of central and marginal populations of Drosophila serrata. I. Genetic variation for stress resistance and species borders. Evolution 47: 1255–1270.[CrossRef][ISI]

Briggs, D., and S. M. Walters. 1984 Plant variation and evolution. Cambridge University Press, Cambridge, UK.

Brown, J. H. 1984 On the relationship between abundance and distribution of species. American Naturalist 124: 255–279.[CrossRef][ISI]

Brussard, P. F. 1984 Geographic patterns and environmental gradients: the central-marginal model in Drosophila revisited. Annual Review of Ecology and Systematics 15: 25–64.[CrossRef][ISI]

Clausen, J., D. D. Keck, and W. M. Hiesey. 1948 Experimental studies on the nature of species. III. Environmental responses of climatic races of Achillea. Carnegie Institute of Washington Publication 581. Washington, DC.

Connell, J. H. 1961 The influence of interspecific competition and other factors on the distribution of the barnacle Chthamalus stellatus. Ecology 42: 710–723.

Diggle, P. K. 1992 Development and the evolution of plant reproductive characters. In R. Wyatt [ed.], Ecology and evolution of plant reproduction, 326–355. Chapman and Hall, New York, NY.

Eckhart, V., M. A. Geber, and C. S. Jonas. 1996 Developmental controls of flowering-time evolution in Clarkia (Onagraceae). American Journal of Botany 83: 81 (Abstract).

Ehleringer, J. R. 1993 Gas-exchange implications of isotopic variation in arid-land plants. In J. A. C. Smith, and H. Griffiths [eds.], Water deficits: plant responses from cell to community, 265–284. BIOS Scientific Publishers, Oxford, UK.

Endler, J. A. 1977 Geographic variation, speciation, and clines. Princeton University Press, Princeton, NJ.

———. 1986 Natural selection in the wild. Princeton University Press, Princeton, NJ.

———. 1995 Multiple-trait coevolution and environmental gradients in guppies. Trends in Ecology and Evolution 10: 22–29.

Ford, E. B. 1954 Problems in the evolution of geographic races. In J. S. Huxley, A. C. Hardey, and E. B. Ford [eds.], Evolution as a process, 99–108. George Allen and Unwin, Ltd., London, UK.

Garcia-Ramos, G., and M. Kirkpatrick. 1997 Genetic models of adaptation and gene flow in peripheral populations. Evolution 51: 21–28.[CrossRef][ISI]

Geber, M. A., and T. E. Dawson. 1990 Genetic variation in and covariation between leaf gas exchange, morphology, and development in Polygonum arenastrum, an annual plant. Oecologia 85: 153–158.[CrossRef][ISI]

———, and ———. 1993 Evolutionary responses of plants to global change. In P. M. Kareiva, J. G. Kingsolver, and R. B. Huey [eds.], Biotic interactions and global change, 179–197. Sinauer, Sunderland, MA.

Gill, L. S., and F. G. Hawksworth. 1961 The mistletoes: a literature review. USDA. Bulletin Number 1242, Washington, DC.

Griffin, J. R. 1977 Oak woodland. In M. G. Barbour and J. Major [eds.], Terrestrial vegetation of California, 383–416. John Wiley&Sons, New York, NY.

Guerrant, E. O. 1984 The role of ontogeny in the evolution and ecology of selected species of Delphinium and Limnanthes. Ph.D. Dissertation, University of California, Berkeley, CA.

Hiesey, W. M., and H. W. Milner. 1965 Physiology of ecological races and species. Annual Review of Plant Physiology 16: 203–216.

Hoffmann, A. A., and M. W. Blows. 1993 Evolutionary genetics and climate change: will animals adapt to global warming? In P. M. Kareiva, J. G. Kingsolver, and R. B. Huey [eds.], Biotic interactions and global change, 165–178. Sinauer, Sunderland, MA.

———, and ———. 1994 Species borders: ecological and evolutionary perspectives. Trends in Ecology and Evolution 9: 223–227.[CrossRef]

———, and P. A. Parsons. 1991 Evolutionary genetics and environmental stress. Oxford University Press, Oxford, UK.

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]

Johansson, M. E. 1994 Life history differences between central and marginal populations of the clonal aquatic plant Ranunculus lingua: a reciprocal transplant experiment. Oikos 70: 65–72[CrossRef][ISI]

Jonas, C. S. 1997 Variation and correlations between traits in the annual plant, Clarkia unguiculata. Master's thesis, Cornell University, Ithaca, NY.

Kirkpatrick, M., and N. H. Barton. 1997 Evolution of a species' range. American Naturalist 150: 1–23.[CrossRef][ISI]

Lacey, E. P. 1988 Latitudinal variation in reproductive timing of a short-lived monocarp, Daucus carota (Apiaceae). Ecology 69: 220–232.[CrossRef][ISI]

Langlet, O. 1971 Two hundred years of genecology. Taxon 20: 653–722.[CrossRef]

Lawton, J. H. 1993 Range, population abundance, and conservation. Trends in Ecology and Evolution 8: 409–413.

Lesica, P., and F. W. Allendorf. 1994 When are peripheral populations valuable for conservation? Conservation Biology 9: 753–760.[CrossRef][ISI]

Lewis, H., and M. E. Lewis. 1955 The genus Clarkia. University of California Press, Berkeley, CA.

Lewontin, R. C. 1974 The genetic basis of evolutionary change. Columbia University Press, New York, NY.

Linhart, Y., and M. C. Grant. 1996 Evolutionary significance of local genetic differentiation in plants. Annual Review of Ecology and Systematics 27: 237–277.[CrossRef][ISI]

Loik, M. E., and P. S. Nobel. 1993 Freezing tolerance and water relations of Opuntia fragilis from Canada and the United States. Ecology 74: 1722–1732.[CrossRef][ISI]

Major, J. 1977 California climate in relation to vegetation. In M. G. Barbour and J. Major [eds.], Terrestrial vegetation of California, 11-74. John Wiley & Sons, New York, NY.

Mayr, E. 1963 Animal species and evolution. Belknap Press, Harvard, Cambridge, MA.

Millar, C. I., and W. J. Libby. 1991 Strategies for conserving clinal, ecotypic, and disjunct population diversity in widespread species. In D. A. Falk and K. E. Holsinger [eds.], Genetics and conservation of rare plants, 149–170. Oxford University Press, Oxford, UK.

Moore, D. M., and H. Lewis. 1965 The evolution of self-pollination in Clarkia xantiana. Evolution 19: 104–114.

Nevo, E., B. F. Carver, and A. Beiles. 1991 Photosynthetic performance in wild emmer wheat, Triticum dicocoides: ecological and genetic predictability. Theoretical and Applied Genetics 81: 445–460.[ISI]

Ohlson, M. 1989 Ecotypic differentiation and phenotypic plasticity in Saxifraga hirculus populations in central and northern Sweden. Holarctic Ecology 12: 46–53.

Reinartz, J. A. 1984a Life history variation of common mullein (Verbascum thapsus). I. Latitudinal differences in population dynamics and timing of reproduction. Journal of Ecology 72: 897–912.[CrossRef][ISI]

———. 1984b Life history variation of common mullein (Verbascum thapsus). II. Plant size and biomass partitioning. Journal of Ecology 72: 913–925.

Rice, K. J., and R. N. Mack. 1991a Ecological genetics of Bromus tectorum. I. A hierarchical analysis of phenotypic variation. Oecologia 88: 77–83.[CrossRef][ISI]

———, and ———. 1991b Ecological genetics of Bromus tectorum. III. The demography of reciprocally sown populations. Oecologia 88: 91–101.[CrossRef][ISI]

Rice W. R. 1989 Analyzing tables of statistical tests. Evolution 43: 223–225.[CrossRef][ISI]

Roach, D. A., and R. D. Wulff. 1987 Maternal effects in plants. Annual Review of Ecology and Systematics 18: 209–235.

Root, T. L. 1993 Effects of global climate change on North American birds and their communities. In P. M. Kareiva, J. G. Kingsolver, and R. B. Huey [eds.], Biotic interactions and global change, 280–292. Sinauer, Sunderland, MA.

Schaal, B. A. 1984 Life-history variation, natural selection, and maternal effects in plant populations. In R. Dirzo and J. Sarukhan [eds.], Perspectives on plant population ecology, 188–206. Sinauer, Sunderland, MA.

Shuter, B. J., and J. R. Post. 1990 Climate, population viability, and the zoogeography of temperate fishes. Transactions of the American Fisheries Society 119: 314–336.[CrossRef][ISI]

St. Clair, R. C., and P. T. Gregory. 1990 Factors affecting the northern range limit of painted turtles Chrysemys picta: winter acidosis or freezing? Copeia 1990: 1083–1089.[CrossRef]

Stanton, M. L., and C. Galen. 1997 Life on the edge: adaptation versus environmentally mediated gene flow in the snow buttercup, Ranunculus adoneus. American Naturalist 150: 143–178.

Stebbins, G. L. 1950 Variation and evolution in plants. Columbia University Press, New York, NY.

Stewart, C. N., and E. T. Nilsen. 1995 Phenotypic plasticity and genetic variation in Vaccinium macrocarpon, the American cranberry. I. Reaction norms of clones from central and marginal populations in a common garden. International Journal of Plant Science 156: 687–697.[CrossRef]

Thorpe, R. S., and M. Baez. 1987 Geographic variation within an island: univariate and multivariate contouring of scalation, size, and shape of the lizard Gallotia galloti. Evolution 41: 256–268.

Turesson, G. 1922 The genotypical response of the plant species to the habitat. Hereditas 3: 211–350.

van Bavel, C. H. M., R. Lascano, and D. R. Wilson. 1978 Water relations of fritted clay. Soil Science Society of America Journal 42: 657–659.[Abstract/Free Full Text]

Vasek, F. C. 1964 The evolution of Clarkia unguiculata derivatives adapted to relatively xeric environments. Evolution 18: 26–42.[CrossRef][ISI]

———. 1968 The relationships of two ecologically marginal, sympatric Clarkia populations. American Naturalist 102: 25–40.[CrossRef][ISI]

———. 1971 Variation in marginal populations of Clarkia. Ecology 52: 1046–1051.

———. 1977 Phenotypic variation and adaptation in Clarkia section Phaeostoma. Systematic Botany 2: 251–279.

———, and R. H. Sauer. 1971 Seasonal progression of flowering in Clarkia. Ecology 52:1038–1045.

von Caemmerer, S., and G. D. Farquhar. 1981 Some relationships between the biochemistry of photosynthesis and the gas exchange of leaves. Planta 153: 386–387.

Winn, A. A., and K. L. Gross. 1992 Latitudinal variation in seed weight and flower number in Prunella vulgaris. Oecologia 93: 55–62.

Woodward, F. I., and B. G. Williams. 1987 Climate and plant distribution at global and local scales. Vegetatio 69: 189–197.[CrossRef][ISI]




This article has been cited by other articles:


Home page
Am. J. Bot.Home page
J. Huang, R. E. J. Boerner, and J. Rebbeck
Ecophysiological responses of two herbaceous species to prescribed burning, alone or in combination with overstory thinning
Am. J. Botany, May 1, 2007; 94(5): 755 - 763.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
D. E. Giblin
Variation in floral longevity between populations of Campanula rotundifolia (Campanulaceae) in response to fitness accrual rate manipulation
Am. J. Botany, October 1, 2005; 92(10): 1714 - 1722.
[Abstract] [Full Text] [PDF]


Home page
ANN BOT (LOND)Home page
J. HERRERA
Flower Size Variation in Rosmarinus officinalis: Individuals, Populations and Habitats
Ann. Bot., February 1, 2005; 95(3): 431 - 437.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
W. S. Armbruster, C. P. H. Mulder, B. G. Baldwin, S. Kalisz, B. Wessa, and H. Nute
Comparative analysis of late floral development and mating-system evolution in tribe Collinsieae (Scrophulariaceae s.l.)
Am. J. Botany, January 1, 2002; 89(1): 37 - 49.
[Abstract] [Full Text] [PDF]


Home page
Am. J. Bot.Home page
A. B. Carroll, S. G. Pallardy, and C. Galen
Drought stress, plant water status, and floral trait expression in fireweed, Epilobium angustifolium (Onagraceae)
Am. J. Botany, March 1, 2001; 88(3): 438 - 446.
[Abstract] [Full Text]


Home page
Am. J. Bot.Home page
C. J. Runions and M. A. Geber
Evolution of the self-pollinating flower in Clarkia xantiana (Onagraceae). I. Size and development of floral organs
Am. J. Botany, October 1, 2000; 87(10): 1439 - 1451.
[Abstract] [Full Text]


This Article
Right arrow Abstract Freely available
Right arrow Full Text (PDF)
Right arrow Submit a response