- Research article
- Open Access
The effect of flower position on variation and covariation in floral traits in a wild hermaphrodite plant
BMC Plant Biologyvolume 10, Article number: 91 (2010)
Floral traits within plants can vary with flower position or flowering time. Within an inflorescence, sexual allocation of early produced basal flowers is often female-biased while later produced distal flowers are male-biased. Such temporal adjustment of floral resource has been considered one of the potential advantages of modularity (regarding a flower as a module) in hermaphrodites. However, flowers are under constraints of independent evolution of a given trait. To understand flower diversification within inflorescences, here we examine variation and covariation in floral traits within racemes at the individual and the maternal family level respectively in an alpine herb Aconitum gymnandrum (Ranunculaceae).
We found that floral traits varied significantly with flower position and among families, and position effects were family-specific. Most of the variance of floral traits was among individuals rather than among flowers within individuals or among families. Significant phenotypic correlations between traits were not affected by position, indicating trait integration under shared developmental regulation. In contrast, positive family-mean correlations in floral traits declined gradually from basal to distal flowers (nine significant correlations among floral traits in basal flowers and only three in distal flowers), showing position-specificity. Therefore, the pattern and magnitude of genetic correlations decreased with flower position.
This finding on covariation pattern in floral reproductive structures within racemes has not been revealed before, providing insights into temporal variation and position effects in floral traits within plants and the potential advantages of modularity in hermaphrodites.
The flower as a reproductive organ can be regarded as a phenotypic module [1, 2], in which floral traits are necessarily correlated with each other because the functional effectiveness of functionally linked floral parts (e.g. precise "fit" with pollinators, or optimal allocation to flower structures) depends on the ability to work together, and fitness depends on their interactions. Floral traits are thus expected to be more canalized and integrated than vegetative traits [3, 4]. Studies based on comparisons between vegetative and floral characters support this expectation [5–11], although some conflicting evidence showed that this pattern may depend on plant species [9, 11].
However, variation in floral traits occurs among populations, individuals, and flowers within one plant [12–18]. For example, Williams & Conner  studied sources of floral variation among different levels in wild radish (Raphanus raphanistrum) and found high intra-plant variance and high floral variation between flowers measured in different weeks. A common pattern observed in hermaphrodites is a reduction in the number or size of reproductive structures in sequentially blooming flowers [19, 20]. This intra-inflorescence variation in floral allocation has been attributed to the effects of resource competition [21–23], architectural effects [20, 24, 25] or mating environments [26, 27]. Although the ability to adjust sex allocation over time is one potential advantage of modularity in hermaphrodites [26, 28, 29], there are few empirical studies that examine temporal variation in floral traits [30, 31], particularly considering trait correlations among sequentially blooming flowers.
Although strong correlation between floral traits may play an important role in the independent evolution of one floral trait , trait correlations frequently appear plastic [2, 33, 34], reflecting the environmental sensitivity. For example, genetic correlations among floral and vegetative traits in Arabidopsis thaliana were significantly influenced by the light environment . Less attention has been paid to the possibility that intrinsic factors (i.e. ontogenetic or position effects) may be an important source of variation in trait correlations when there is great intra-individual variation of floral traits [31, 35]. In an early multilevel analysis in Dalechampia scandens, Armbruster  showed that covariance between two of three measured reproductive traits was evenly distributed among three levels (within genets, among genets and among populations), although one of the traits was lack of ontogenetically-induced variation. A study on plasticity of floral traits of Campanula rapunculoides with respect to genotype, environment and ontogeny suggested that correlation patterns among floral traits depended on the environment and the trade-off relationship between male and female function was confounded by ontogenetic effects . Considering ontogenetic variation in floral traits of Iris gracilipes, Ishii & Morinaga  found the pattern of correlations within individuals was basically similar to that among individuals. Comparing correlations among floral traits in closely related Nicotiana species, Bissell & Diggle  provided evidence of common developmental regulation of correlated traits but there were some independent trait variation with flower position and age. Although ontogenetic or positional effects have been considered in relation to covariation of floral traits, the importance of this aspect is still underappreciated for a full understanding of variation in floral traits. A more powerful approach involving calculation of genetic variance and covariance by measuring multiple traits in the offspring of full-sib or half-sib families has seldom been adopted in such studies. To our knowledge, only Mazer & Delesalle  have measured floral traits in half-sib families in Spergularia marina (Caryophyllaceae). They found that the statistical significance of correlations among six floral traits (number of ovules, normal anthers, abnormal anthers and petals, single petal area and total petal area) changed over sampling time.
In a previous study we observed variation in floral allocation in an alpine protandrous herb Aconitum gymnandrum (Ranunculaceae) with a racemose inflorescence, in which floral sex allocation shifted from female-biased to male-biased from basal to distal flowers . To demonstrate the genetic basis of such temporal shift of floral traits that occurs in many plants, here we estimate broad sense heritability by a common garden experiment. We examine effects of flower position on floral traits and trait correlations in the herb, respectively at the individual level and at the maternal family level. Specifically, we (1) address variation of floral traits within racemes and among families, and examine whether position effects are family-specific; (2) explore correlations among floral traits at the individual and family levels, and their relationship to flower position, to examine whether correlations of floral traits vary within racemes.
Variation of floral traits with flower position and among families
We examined the basal three, middle three and distal three flowers on each raceme in Aconitum gymnandrum. We measured seven traits in each raceme from 100 plants from 25 families that we cultivated in 100 pots, including height of the sepal galea, dry mass of floral structures (androecium, gynoecium and calyx), stamens and carpels per flower, and total flowers per raceme (see Additional file 1). Most of the floral traits varied significantly within racemes and/or among families. Effects of position and family on galea height (F2, 99 = 4.023, p = 0.011; F24, 99 = 2.812, P = 0.012, respectively), anther number (F2, 99 = 7. 654, P = 0.002; F24, 99 = 2.457, P = 0.021) and carpel number (F2, 99 = 86.715, P < 0.0001; F24, 99 = 2.703, P = 0.016) were significant (repeated measures ANOVA). For gynoecium mass, effect of position was significant (F2, 99 = 12. 349, P < 0.0001), while effect of family was not (F24, 99 = 1.714, P = 0.141). Significant position × family interactions were found in androecium mass (F44, 99 = 1.958, P = 0.022), anther number (F44, 99 = 2. 447, P = 0.003) and calyx mass (F44, 99 = 2.186, P = 0.006), showing that the intra-raceme variation was family-specific.
Significant among-family variation in galea height, anther number and carpel number was independent of flower position within racemes (Table 1). Among-family variations in androecium mass, gynoecium mass and calyx mass were not significant, and this finding is consistent across different positions.
Variance components of floral traits
The largest variation in floral traits was among individuals (variance components from 0.50 to 0.81). Among-flower variation in these traits (variance components from 0.15 to 0.40) was higher than among-family variation (variance components from 0.04 to 0.29), except for galea height and androecium mass (Figure 1).
Trait correlation and flower position
Partial correlation analysis indicated significant phenotypic correlations between most traits except for galea height (Table 2). The amounts of significant correlations were independent of flower positions, and phenotypic correlations did not vary within racemes. CPC analysis indicated that phenotypic correlations were similar across different positions with respect to equality, proportionality, and shared common principal components (P > 0.1 in all cases for the step-up and jump-up methods). The magnitude of the mean absolute correlation coefficients did not change significantly within racemes (Figure 2, basal vs. middle, P = 0.28, basal vs. distal, P = 0.253, bootstrapping estimation). Galea height was more closely correlated with gynoecium and calyx mass than it was with androecium mass.
The relative numbers of significant correlations among family means depended strongly on flower position within racemes (Table 2 and Figure 3): there were nine significant correlations among floral traits in basal flowers, but only three in distal flowers. Correlations between androecium mass, anther number, gynoecium mass, carpel number and calyx mass were positive in basal flowers, but these correlations with calyx mass were not significant in either middle or distal flowers (Figure 3). Correlations between androecium mass and anther number, gynoecium mass and carpel number occurred in both basal and middle flowers, but were lost in distal flowers. And mean absolute coefficients of family-mean correlation declined significantly from basal to distal positions within racemes (Figure 2, basal vs. middle, P = 0.046, basal vs. distal, P = 0.0012, bootstrapping estimation). CPC analysis showed that family-mean correlation matrices were not similar or not proportional across positions (P < 0.01, using either methods), but shared common principal components (P = 0.760 using the step-up method, P = 0.696 using the jump-up method).
Significantly positive correlations were observed between male (androecium mass, anther number) and female (gynoecium mass, carpel number) allocation (Figure 3 and Table 2). However, the significant correlations depended in part on flower position within racemes, e.g. anther number (but not androecium mass) was significantly correlated with gynoecium mass irrespective of flower position (Figure 3 and Table 2).
Variation in floral traits
We demonstrated that most floral traits varied significantly within racemes in A. gymnandrum, indicating that traits were affected by flower position. Positional or temporal effects on floral traits were also found in Dalechampia [12, 38], Spergularia marina [14, 15], Raphanus raphanistrum  and in Nicotiana . Temporal variation in floral trait expression could mask underlying genetic differences among individuals, families, populations or species in floral phenotype [15, 35, 39], resulting in sensitivity of heritability estimates to the position of measured flowers. In our studied herb, significant variation among families for most floral traits (except for gynoecium mass) was independent of flower position within racemes. The magnitude of differences among populations and among maternal family means in most floral traits of S. marina, however, was observed to change through the flowering period . We observed significant position × family interactions in androecium mass, anther number and calyx mass, showing that the intra-raceme pattern in A. gymnandrum was family-specific. A genotype-specific intra-inflorescence pattern was observed for pollen and ovule number per flower in Fragaria virginiana (Rosaceae) as shown by significant genotype-by-flower position interactions . In S. marina, no family-specific effects of sampling week on floral traits were detected except in the case of ovule number per flower which showed significant week × family interactions . The family-differential intra-raceme patterns in A. gymnandrum suggest that genetic sources of variation in floral traits contribute significantly to floral phenotype.
It has been suggested that larger variation in floral traits within individuals than among individuals may cushion differential selection pressures on floral evolution [14, 17]. For example, if pollinator-mediated selection acts as one force influencing floral variation, the evolutionary response to selection will be weakened if there is a wide range of floral variation within individuals rather than among individuals. But evidence supporting this expectation has been scarce [12, 41]. We found that variation in floral traits was greater among plants than among flowers within individuals or among families in bumblebee-pollinated Aconitum gymnandrum, consistent with the expectation. In contrast to plants pollinated by specialists, studies on generalist-pollinated species showed intra-plant variation in floral traits to be larger than inter-plant variation [16, 17]. For example, Williams & Conner  found that high within-plant variance of floral traits in the field accorded with the weak selection on floral traits in wild radish which was pollinated by diverse insects. In another Brassicaceae species Brassica rapa, however, the inter-plant variance of floral traits was similar to the intra-plant component . Variation in floral traits was also similar between intra- and inter-plant levels in Iris gracilipes . Given that few studies consider variation in floral traits both within and among individuals, at present it is impossible to assess whether the pattern of high inter-plant variation in floral traits is associated with the pollinator specialization level.
Variation in trait correlation
We observed that significant phenotypic correlations between most traits (except for galea height) were stable and independent of flower position. Although these traits independently varied among positions as exhibited above, it is not great enough to decouple the strong correlation between traits. Similarly, a study on Nicotiana exhibited stability of correlation patterns of length characters rather than width characters of flowers independent of cyme position, flower position and flower age, showing a pattern of shared developmental regulation . Variation in floral traits rather than trait correlations across environments has also been documented in other studies [41, 43]. Constant correlation patterns prove that these floral traits are developmental integrated; on the other hand, if there is indirect or correlated selection on one of the floral traits in A. gymnandrum, it would not be strongly influenced by flower position within racemes.
In contrast, we observed that significant correlations among family means depended strongly on flower position in A. gymnandrum, and the correlation matrices varied among different positions. Floral traits of basal flowers were more significantly correlated than traits of distal flowers. Positive family-mean correlations between floral traits weakened gradually from basal to distal flowers, i.e. the strength of correlations were position-specific. This suggests that later flowers within racemes are subject to less genetic constraint or developmental stability than earlier produced distal flowers. At the same time, it also indicates a high potential for evolutionary modification because of the low genetic correlation among traits; variation observed among maternal sibships could reflect genetic differences to some extent [38, 44, 45]. To our knowledge, the position-dependent decline in the magnitude of family-mean correlations rather than phenotypic correlations has not been reported before. Our result is mirrored in Spergularia marina in which the phenotypic and among-family correlations among traits changed over the sampling period although the trend was inconsistent . It has been suggested that genetic correlations are variable and could be affected by the developmental age of an individual or by environmental factors [46, 47]. For example, genetic correlations of life history traits in Geranium changed from generally negative in the early juvenile stage to strongly positive in the adult . Genetic correlations among floral and vegetative traits were significantly influenced by the light environment in Arabidopsis thaliana . The variability of genetic correlation is critical for understanding trait evolution under natural selection.
The greater variation of floral traits in later flowers compared to early flowers of A. gymnandrum may be in part due to weakened trait correlations of family means (i.e. weakened genetic correlation). Strong genetic correlations would restrict independent evolution of floral traits [31, 49], but intra-raceme variation in genetic correlations of floral traits would cushion selection on floral traits and trait integration within racemes, and consequently on the intra-raceme pattern of A. gymnandrum. For example, at an early developmental stage of the raceme, evolutionary changes of basal flowers may be constrained in terms of high genetic correlation; at a later stage, evolutionary changes of distal flowers may be less constrained due to low genetic correlation. It has been suggested that floral traits and their correlations could be regulated at the level of individual flowers [18, 28], and that developmental plasticity and canalization are independently controlled . Genetic factors and diverse environmental factors have been demonstrated to increase developmental plasticity in plants . In andromonoecious species, basal flowers invariably developed into perfect flowers but distal flowers were plastic, capable of developing into either a staminate flower or a perfect flower . Such temporal variation of later-produced flowers, with the potential to specialize as males , can be better understood if the intra-inflorescence pattern of trait correlations can be revealed, as in our demonstration in A. gymnandrum that later or distal flowers were more labile than early basal flowers. Thus, our findings provide insight into the plasticity of sexual expression and regulation of floral development in hermaphrodite plants.
In conclusion, floral traits in A. gymnandrum varied significantly within racemes and among maternal families, and the position effect was family-specific. Variation of floral traits among individuals was greater than that among flowers within individuals or among families. Although significant phenotypic correlations in floral traits did not change among flower positions, the pattern of family-mean correlations varied and the magnitude declined gradually from basal to distal flowers, exhibiting position-specific correlations. This shows that flowers within an individual have dissimilar evolutionary potential. Thus, data from A. gymnandrum suggest that position effect on the magnitude of covariation in floral traits would confound the evolution of different flowers within individuals and consequently influence selection on intra-raceme pattern in flowering plants. Our findings on variation and covariation in floral traits with position provide genetic basis of flower functional specialization within individuals [20, 26].
Aconitum gymnandrum Maxim. (Ranunculaceae) is an annual herb, widely distributed in alpine meadows (1600-3800 m a.s.l) in the Qinghai-Tibet Plateau, China. Individual plants generally produce one erect raceme consisting of 2-30 blue-purple zygomorphic flowers, which open sequentially from bottom to top. Each flower has 6-14 separate carpels (each with 8-14 ovules) surrounded by 30-90 stamens. The galea (or hood), formed from one of five petaloid sepals, contains two stalked petals with nectaries. The species is self-compatible, strongly protandrous like other related species in the same genus  and bumblebee-pollinated. The anthers dehisce over 4-5 days and stigmas become receptive 1-2 days later. Plants commonly bloom from June through August and single flowers last 6-10 days. Fruit maturation requires 20-30 days.
In this experiment seedlings from different maternal plants were transplanted into pots to investigate among-family variation and covariation of floral traits. In September 2003, we collected seeds from 25 maternal plants of A. gymnandrum from populations of an alpine meadow in northeastern Qinghai-Tibet Plateau near Hezuo County, Gansu Province, China. Plants were at least 10 m apart to ensure genetic differences among individuals. Seeds were stored in envelopes at room temperature.
On 7 May 2004, seeds from the 25 maternal plants were germinated in Petri dishes with distilled water. Eight ten-day-old seedlings of uniform size were selected from each family, and seedlings were transplanted in pairs into assigned positions in each of four plastic pots (diameter 26 cm) filled with mixed soil collected from a local natural site where A. gymnandrum grows naturally. The soil had previously been mixed and covered with film for four months to eliminate preexisting seeds. These plants were used to measure floral traits when they flowered. Irrespective of families, all pots were randomly arranged at the field station of Lanzhou University at Hezuo County (E102°53', N34°55') and exposed to the natural environment. This design was to simulate natural growing conditions in the field.
In July 2005, all of the experimental plants bloomed, and natural mortality was less than 10%. For the pots of 25 families, one plant per pot was picked randomly, resulting in a sample of four for each family. We collected all sequentially blooming flowers within racemes of a total of 100 plants from 25 families, when the flowers had just opened. Racemes of all plants contained more than 9 flowers. We grouped flowers into three positions within racemes: the basal three, the middle three and the distal three on each raceme. For each flower we measured the height of the sepal galea (from the base of the sepal to the top of the galea, to the nearest 0.01 mm, with a vernier caliper), weighed the dry mass of floral structures (androecium, gynoecium and calyx, after drying at 80°C for 24h), and counted stamens and carpels per flower and the total flowers per raceme.
Variation in the traits within racemes and variation of all traits among families were analyzed by repeated measures ANOVA (the GLM procedure of SPSS), with flower position as within-subject factor (three levels) and family as subject factor. We also analyzed variation of floral traits among families with one-way ANOVA (GLM model of SPSS for Windows) at different positions separately, with family as a random factor. To examine the relative magnitudes of sources of variation in floral traits, among flowers within individuals, among individuals and among families, variance components from fully nested random models were estimated using REML (mixed model of SPSS).
Family-mean correlations between traits were estimated separately at each position; each family's mean was calculated from the phenotypic means of the individuals representing it. Family-mean correlations reflect the genetic base to some extent and can be regarded as genetic correlation [38, 44, 45]. To remove the effect of difference in plant size, we controlled for differences among families in plant size by estimating partial correlation coefficients between traits, controlling for total flower number. We also estimated phenotypic correlations between traits across positions, removing the effect of difference in plant size using the same methods as we used for the among-family correlations. To determine if correlations (genetic and phenotypic) among floral traits vary with flower position, common principal component analysis (Flury hierarchical method) was used to compare these correlation matrices between positions by the step-up method and the jump-up method [52, 53]. Differences between flower positions in the mean absolute coefficients of family-mean and phenotypic correlations were estimated by a bootstrapping method (n = 10 000 permutations, using Data Pilot ver.1.03, ), because the values used in calculating correlation coefficients were not independent. To reduce the probability of spurious results caused by the simultaneous evaluation of multiple statistical tests, sequential Bonferroni corrections  were used in each correlation matrix (family correlations at each of flower positions separately), with a table-wide significance value of p < 0.05. All analyses were completed in SPSS statistical software (Version 12.0 for Windows).
Wagner GP: Homologues, natural kinds and the evolution of modularity. Am Zool. 1996, 36: 36-43.
Murren CJ: Phenotypic integration in plants. Plant Species Biol. 2002, 17: 89-99. 10.1046/j.1442-1984.2002.00079.x.
Berg RL: The ecological significance of correlation pleiades. Evolution. 1960, 14: 171-180. 10.2307/2405824.
Bradshaw AD: Evolutionary significance of phenotypic plasticity in plants. Adv Genet. 1965, 13: 115-155. full_text.
Conner JK, Via S: Patterns of phenotypic and genetic correlations among morphological and life-history traits in wild radish, Raphanus raphanistrum. Evolution. 1993, 47: 704-711. 10.2307/2410086.
Conner JK, Sterling A: Testing hypotheses of functional relationships - a comparative survey of correlation patterns among floral traits in five insect-pollinated plants. Am J Bot. 1995, 82: 1399-1406. 10.2307/2445866.
Sherry RA, Lord EM: Developmental stability in flowers of Clarkia tembloriensis (Onagraceae). J Evol Biol. 1996, 9: 911-930. 10.1046/j.1420-9101.1996.9060911.x.
Waitt DE, Levin DA: Genetic and phenotypic correlations in plants: a botanical test of Cheverud's conjecture. Heredity. 1998, 80: 310-319. 10.1046/j.1365-2540.1998.00298.x.
Armbruster WS, Stilio VSD, Tuxill JD, Flores TC, Runk JLV: Covariance and decoupling of floral and vegetative traits in nine neotropical plants: a reevaluation of Berg's correlation-pleiades concept. Am J Bot. 1999, 86: 39-55. 10.2307/2656953.
Brock MT, Weinig C: Plasticity and environment-specific covariances: an investigation of floral-vegetative and within flower correlations. Evolution. 2007, 61: 2913-2924. 10.1111/j.1558-5646.2007.00240.x.
Hansen TF, Pélabon C, Armbruster WS: Comparing variational properties of homologous floral and vegetative characters in Dalechampia scandens: testing the Berg hypothesis. Evol Biol. 2007, 34: 86-98. 10.1007/s11692-007-9006-3.
Armbruster WS: Multilevel analysis of morphometric data from natural plant populations: insights into ontogenetic, genetic, and selective correlations in Dalechampia scandens. Evolution. 1991, 45: 1229-1244. 10.2307/2409730.
Campbell DR: Variation in sex allocation and floral morphology in Ipomopsis aggregata (Polemoniaceae). Am J Bot. 1992, 79: 516-521. 10.2307/2444862.
Mazer SJ, Delesalle VA: Covariation among floral traits in Spergularia marina (Caryophyllaceae): geographic and temporal variation in phenotypic and among-family correlations. J Evol Biol. 1996, 9: 993-1015. 10.1046/j.1420-9101.1996.9060993.x.
Mazer SJ, Delesalle VA: Floral trait variation in Spergularia marina (Caryophyllaceae): ontogenetic, maternal family, and population effects. Heredity. 1996, 77: 269-281. 10.1038/hdy.1996.141.
Dominguez CA, Eguiarte LE, Nunez-Farfan J, Dirzo R: Flower morphometry of Rhizophora mangle (Rhizophoraceae): geographical variation in Mexican populations. Am J Bot. 1998, 85: 637-643. 10.2307/2446532.
Williams JL, Conner JK: Sources of phenotypic variation in floral traits in wild radish, Raphanus raphanistrum (Brassicaceae). Am J Bot. 2001, 88: 1577-1581. 10.2307/3558401.
Ishii HS, Morinaga SI: Intra- and inter-plant level correlations among floral traits in Iris gracilipes (Iridaceae). Evol Ecol. 2005, 19: 435-448. 10.1007/s10682-005-0896-1.
Stephenson AG: Flower and fruit abortion: proximate causes and ultimate functions. Annu Rev Ecol Syst. 1981, 12: 253-279. 10.1146/annurev.es.12.110181.001345.
Diggle PK: Architectural effects and the interpretation of patterns of fruit and seed development. Annu Rev Ecol Syst. 1995, 26: 531-552. 10.1146/annurev.es.26.110195.002531.
Herrera J: Allocation of reproductive resources within and among inflorescences of Lavandula stoechas (Lamiaceae). Am J Bot. 1991, 78: 789-794. 10.2307/2445069.
Ashman TL, Baker I: Variation in floral sex allocation with time of season and currency. Ecology. 1992, 73: 1237-1243. 10.2307/1940672.
Guitián J, Navarro L: Allocation of reproductive resources within the inflorescences of Petrocoptis grandiflora (Caryophyllaceae). Can J Bot. 1996, 74: 1482-1486. 10.1139/b96-178.
Wolfe LM: Why does the size of reproductive structures decline through time in Hydrophyllum appendiculatum (Hydrophyllaceae): developmental constraints vs. resource limitation?. Am J Bot. 1992, 79: 1286-1290. 10.2307/2445057.
Diggle PK: Ontogenetic contingency and floral morphology: the effects of architecture and resource limitation. Int J Plant Sci. 1997, 158: 99-107. 10.1086/297510.
Brunet J, Charlesworth D: Floral sex allocation in sequentially blooming plants. Evolution. 1995, 49: 70-79. 10.2307/2410293.
Huang SQ, Tang LL, Yu Q, Guo YH: Temporal floral sex allocation in protogynous Aquilegia yabeana contrasts with protandrous species: support for the mating environment hypothesis. Evolution. 2004, 58: 1131-1134.
Lloyd DG: Sexual strategies in plants. I. An hypothesis of serial adjustment of maternal investment during one reproductive session. New Phytol. 1980, 86: 69-79. 10.1111/j.1469-8137.1980.tb00780.x.
Charnov EL: The theory of sex allocation. New Jersey, Princeton: Princeton University Press, 1982, 355-
Delesalle VA, Mazer SJ, Paz H: Temporal variation in the pollen: ovule ratios of Clarkia (Onagraceae) taxa with contrasting mating systems: field populations. J Evol Biol. 2008, 21: 310-323.
Bissell EK, Diggle PK: Floral morphology in Nicotiana: architectural and temporal effects on phenotypic integration. Int J Plant Sci. 2008, 169: 225-240. 10.1086/523875.
Ashman TL, Majetic CJ: Genetic constraints on floral evolution: a review and evaluation of patterns. Heredity. 2006, 96: 343-352. 10.1038/sj.hdy.6800815.
Schlichting CD: Phenotypic integration and environmental change. Bioscience. 1989, 39: 460-464. 10.2307/1311138.
Waitt DE, Levin DA: Phenotypic integration and plastic correlation in Phlox drummondii (Polemoniaceae). Am J Bot. 1993, 80: 1224-1233. 10.2307/2445552.
Diggle PK: Architectural effects on floral form and function: a review. Deep Morphology: Toward a Renaissance of Morphology in Plant Systematics. Edited by: Stuessy T, Horandl E, Mayer V. 2003, Koeltz, Konigstein, 63-80.
Vogler DW, Peretz S, Stephenson AG: Floral plasticity in an iteroparous plant: the interactive effects of genotype, environment, and ontogeny in Campanula rapunculoides (Campanulaceae). Am J Bot. 1999, 86: 482-494. 10.2307/2656809.
Zhao ZG, Meng JL, Fan BL, Du GZ: Reproductive patterns within racemes in protandrous Aconitum gymnandrum (Ranunculaceae): potential mechanism and among-family variation. Pl Syst Evol. 2008, 273: 247-256. 10.1007/s00606-008-0002-6.
Pélabon C, Hansen TF, Carlson ML, Yoccoz NG, Armbruster WS: Consequences of inter-population crosses on developmental stability and canalization of floral traits in Dalechampia scandens (Euphorbiaceae). J Evol Biol. 2004, 17: 19-32. 10.1046/j.1420-9101.2003.00663.x.
Stephenson AG, Devlin B, Horton JB: The effects of seed number and prior fruit dominance on the pattern of fruit production in Cucurbita pepo (zucchini squash). Ann Bot. 1988, 62: 653-661.
Ashman TL, Hitchens MS: Dissecting the causes of variation in intra-inflorescence allocation in a sexually polymorphic species, Fragaria virginiana (Rosaceae). Am J Bot. 2000, 87: 197-204. 10.2307/2656906.
Møller AP, Shykoff JA: Morphological developmental stability in plants: patterns and causes. Int J Plant Sci. 1999, 160 (Suppl): S135-S146. 10.1086/314219.
Syafaruddin , Yoshioka Y, Horisaki A, Niikura S, Ohsawa R: Intraspecific variation in floral organs and structure in Brassica rapa L. analyzed by principal component analysis. Breeding Sci. 2006, 56: 189-194. 10.1270/jsbbs.56.189.
Caruso CM: Plasticity of inflorescence traits in Lobelia siphilitica (Lobeliaceae) in response to soil water availability. Am J Bot. 2006, 93: 531-538. 10.3732/ajb.93.4.531.
Mazer SJ, Dawson KA: Size-dependent sex allocation within flowers of the annual herb Clarkia unguiculata (Onagraceae): ontogenetic and among-plant variation. Am J Bot. 2001, 88: 819-831. 10.2307/2657034.
Parachnowitsch AL, Elle E: Variation in sex allocation and male-female trade-offs in six populations of Collinsia parviflora (Scrophulariaceae s. l.). Am J Bot. 2004, 91: 1200-1207. 10.3732/ajb.91.8.1200.
Fenster CB, Carr DE: Genetics of sex allocation in Mimulus (Scrophulariaceae). J Evol Biol. 1997, 10: 641-661. 10.1007/s000360050046.
Stearns S, de Jong G, Newman B: The effects of phenotypic plasticity on genetic correlations. Trends Ecol Evol. 1991, 6: 122-126. 10.1016/0169-5347(91)90090-K.
Roach DA: Life history variation in Geranium carolinianum. 1. Covariation between characters at different stages of the life cycle. Am Nat. 1986, 128: 47-57. 10.1086/284538.
Caruso CM: The ecological genetics of floral traits. Heredity. 2006, 97: 86-87. 10.1038/sj.hdy.6800853.
Diggle PK: The expression of andromonoecy in Solanum hirtum (Solanaceae): phenotypic plasticity and ontogenetic contingency. Am J Bot. 1994, 81: 1354-1365. 10.2307/2445411.
Bosch M, Simon J, Molero J, Blanché C: Breeding systems in the tribe Delphinieae (Ranunculaceae) in the western Mediterranean area. Flora. 2001, 196: 101-113.
Phillips PC, Arnold SJ: Hierarchical comparison of genetic variance-covariance matrices. I. using the Flury hierarchy. Evolution. 1999, 53: 1506-1515. 10.2307/2640896.
Roff DA, Mousseau T: The evolution of the phenotypic covariance matrix: evidence for selection and drift in Melanoplus. J Evol Biol. 2005, 18: 1104-1114. 10.1111/j.1420-9101.2005.00862.x.
Anonymous: Data Pilot v. 1.03 (developed by Oleg Adibekov). Two Pilots Inc: 2003.
Rice WR: Analyzing tables of statistical tests. Evolution. 1989, 43: 223-225. 10.2307/2409177.
We thank Man-Tang Wang and Xian-Hui Zhou for help in the field and lab, Lynda Delph, Charles Fenster and three anonymous reviewers for providing valuable comments on this work, and Sarah Corbet for correcting the English. This work was funded by the Natural Science Foundation of China (grant no. 30900162) to ZGZ and NSFC grants no. 30770135, 30825005 to SQH.
ZGZ, GZD and SQH designed the research. ZGZ performed the research and analyzed the data. ZGZ drafted the manuscript. SQH contributed to writing the manuscript. All authors read and approved the final version.