1. 引言
花部性状特别是花色在吸引传粉者方面起着非常重要的作用。花色的呈现具有一定的遗传基础，且需要环境的参与。植物生理生化环境和花青素代谢途径上结构或调控基因的突变均可能带来花色表型的变异。自然界花色表型具有多样性，主要体现在花色参数和着色模式上。但许多花色变异模式的存在机理还不清楚，花色变异与自然选择之间的关系也需要深入研究。本文就花色呈现的生理生化基础、已知的花色变异遗传机制和传粉者及非传粉介质(如基因多效性、环境因子、病原菌、啃食者等)对花色进化的影响进行总结，以利于探讨花色多样性存在的机制。
2. 花色呈现的分子基础
2.1. 花色呈现的生化生理基础
花色呈现物质主要有三大类：甜菜碱(betalains)、类胡萝卜素(carotenoids)和类黄酮(flavonoids) [1] [2] [3] 。其中甜菜碱、类胡萝卜素和部分类黄酮(如黄酮、黄酮醇、橙酮)可使花色呈现黄、红等色，类黄酮中的花青素是花色多样性最主要的物质基础。不同类型的花青素可使花呈现红、粉、紫、蓝等颜色。花青素在细胞质中不稳定而主要以花色苷(配糖体)的形式存在于植物液泡中 [1] [4] 。花青素苷元主要有天竺葵素(pelargnidin)、矢车菊素(cyanidin)、飞燕草素(delphinidin)、芍药花素(peonidin)、碧冬茄素(petunidin)和锦葵素(malvidin) 6大类 [5] [6] [7] 。它们均具有C6-C3-C6的基本骨架 [3] [8] [9] ，结构上的不同主要体现在B-环上羟基化、甲基化的程度。羟基越多，花色越蓝；甲基多则花色趋红。花青素苷元还会发生不同程度的糖基化和酰基化使其更稳定。糖基化使花色发生轻微红移；被香豆酸、肉桂酸等芳香族氨基酸酰化会使花色发生蓝移 [3] [10] 。羟基化、甲基化、糖基化和酰基化等修饰使花青素种类丰富多样，赋予了花色的多样性。随着技术的进步，越来越多的花青素被分离鉴定 [11] [12] [13] 。至2007年，自然界共有约500种花青素被鉴定 [14] 。除了花青素类别，花瓣细胞形状 [15] [16] [17] 、液泡pH [18] - [25] 、细胞共色素及金属离子的存在 [3] [12] 都会影响花色的呈现。如飞燕草素在较高的Fe3+下会使郁金香(Tulipa gesneriana)呈现蓝色 [26] [27] 。
2.2. 花色呈现的遗传基础
花青素的合成是研究比较清楚的一个次生代谢通路(图1) [3] [28] [29] 。它以香豆酰辅酶A(4-coumaroyl-CoA)
CHS：查尔酮合酶(chalcone synthase)；CHI：查尔酮异构酶(chalcone isomerase)；F3H：黄烷酮-3-羟化酶(flavanone 3-hydroxylase)；F3’H：类黄酮-3’-羟化酶(flavonoid 3’-hydroxylase)；F3’5’H：类黄酮-3’5’-羟化酶(flavonoid 3’5’-hydroxylase)；DFR：二氢黄酮醇还原酶(dihydroflavonol reductase)；ANS：花青素合酶(anthocyanidin synthase)；UF3GT：尿苷二磷酸–葡萄糖–类黄酮-3-葡糖基转移酶(UDP glucose flavonoid-3-glucosyltransferase)
Figure 1.Anthocyanin biosynthetic pathway (modified from [132])
图1. 花青素代谢途径(修改自[132])
和丙二酰辅酶A(malonyl-CoA)为底物，由结构基因(CHS、CHI、F3H、F3’H、F3’5’H、DFR和ANS)编码的一系列酶催化合成，并经糖基化、甲基化和酰基化等修饰后被转运至液泡中。其中F3H、F3’H和F3’5’H为花青素代谢途径的三个分支酶，其合成产物分别对应于天竺葵素、矢车菊素和飞燕草素。这些酶基因在许多植物中均已被克隆 [6] [28] 。它们可能是以多酶复合物的形式在内质网膜上行使催化功能 [2] [3] [30] 。花青素代谢与初级代谢之间具有密切的联系。来自圆叶牵牛的数据显示其底物香豆酰辅酶A与丙二酰辅酶A分别来自于苯丙氨酸(phenylalanine)合成途径和糖酵解(glycolysis)过程 [31] 。
已知花青素代谢途径由MYB、碱性螺旋–环–螺旋(bHLH)和WD40重复蛋白(WDR)三类调控因子主要以MBW三元复合物的形式调控 [3] [22] 。其中WDR的表达具有广泛性，而MYB和bHLH的表达却仅限于特定组织 [32] - [38] 。除参与花青素代谢途径外，不同组合的MBW复合物还参与液泡pH的控制、原花青素的合成、根毛形成、表皮细胞形状的决定和着色模式等过程。其中MYB决定了其功能的特异性 [33] [34] [39] 。在MBW复合物中，MYB和bHLH可分别特异地结合结构基因启动子区的cis元件——ANCNNCC (MRE, MYB-recognizing element)和CACN(A/C/T) (G/T) (BRE, bHLH-recognizing element)，且MRE与BRE之间偏爱6 bp的间隔 [29] [40] [41] 。
但一些MYB可单独起作用。如在玉米(Zea mays)中，P1(subgroup7)可单独激活花青素合成的一部分基因(CHS, CHI, DFR, FLS)而促进玉米粒中鞣酐(phlobaphene)的合成 [25] [36] [42] [43] 。在拟南芥(Arabidopsis thaliana)中，CHS、CHI、F3H、F3’H、FLS可被一些冗余的R2R3MYB (MYB11, MYB12, MYB111) (subgroup7)单独激活 [44] 。
3. 常见花色表型变异的遗传机制
花色表型多样性不仅体现在花色参数上，多样的着色模式(图2)更加丰富了花色的多态。这些变异可
a. 色折(color break)：水仙(Narcissus pseudonarcissus) [61] 。b-f. 彩斑(variegation)：b. 圆叶牵牛(Ipomoea purpurea) [62] ；c-e. 裂叶牵牛(I. nil) [63] ；f. 大丽花(Dahlia variabilis) [64] 。g-j. 缘环(marginal picotee)：矮牵牛(Petunia hybrid) [65] 。k. 双色花(bicolor)大丽花(Dahlia variabilis) [66] 。l. 星状(star-type)：矮牵牛(Petunia hybrid) [67] 。m-n. 紫外模式(UV/bull’s eye)：黑心金光菊(Rudbeckia hirta)：m. 可见光下；n. 紫外光下 [68] 。o-q. 脉色(venation)：o. 金鱼草(Antirrhinum majus) [69] ；p. 矮牵牛(Petunia hybrid)；q. 大花蕙兰(Cymbidium hybrida) [70] 。r-s. 芽红(bud flush)：r. 苹果(Malus × domestica) [71] ；s. 岷江百合(Lilium regale) [71] 。t-v. 斑点(spot)：t. 细长山字草(Clarkia gracilis) [72] ；u. 郁金香(Tulipa gesneriana cv. Murasakizuisho) [27] ；v. 蝴蝶兰(Phalaenopsis spp. I-Hsin Sun Beauty “KHM1065”) [73] 。w. 抬升的斑点(raised spot)：东方百合(Oriental hybrid lily “Le Reve”) [74] 。x. 飞溅的斑点(splatter spot)：亚洲百合(Asiatic hybrid lily “Latvia”) [74] 。y. 束状斑点(brushmarks)：亚洲百合(Asiatic hybrid lily “Centrefold”) [74] 。z. 晕状斑点(spot-halo)：毛地黄(Digitalis purpurea) [69]
Figure 2. Reported types of pigmentation pattern
图2.已报道的花色着色模式类型
能有利于增加植物的环境适应性。花色的变异具有一定的遗传基础，花青素代谢途径上结构或调控基因的突变(包括编码区和调控区)都或多或少会影响基因的功能或表达量而影响花色，基因的组织特异性表达或功能分化也会影响花色表型。
3.1. 花色的获得
花青素代谢为类黄酮的一个代谢分支，根据类黄酮类物质和代谢途径基因类似片段的存在与否，人们推断类黄酮代谢途径应该起源于植物登陆 [45] 。一些基因功能的获得可能源于基因重复后的功能分化，如CHI被推测起源于脂肪酸合酶基因功能分化 [46] 。从产物类别来看，在原始植物角苔类(liverworts)、苔藓(mosses)和石松类(club mosses)中发现有黄烷酮、黄酮和黄酮醇物质的出现 [47] 。在蕨类(ferns)植物里出现原花青素 [47] ，而花青素则从在松柏类裸子植物中出现到广泛存在于被子植物中 [48] [49] 。在调控基因方面，由于MYB较强的功能和组织特异性，花色素从无到有或从少到多的获得转换有可能通过cis或trans的方式使激活型MYB高表达 [50] [51] [52] [53] [54] 或抑制型MYB低表达 [55] 来实现(表1)。
3.2. 从蓝/紫色到红/粉色的花色变异
花色从蓝/紫色到红/粉色的变异通常是由于花青素从飞燕草素或矢车菊素到天竺葵素的转变，即花青素羟基数目的减少，因而常伴随羟基化酶基因F3’H或F3’5’H的功能失活或表达量降低(表2)。牵牛花仅含有矢车菊素和天竺葵素 [56] [57] [58] ，对于本身就具有蓝/紫色花的物种如圆叶牵牛(Ipomoea purpurea)、
Table 1. Summary of MYB mutations for gain of pigment
表1. 花色获得的基因突变形式
COD代表编码区的突变；CIS代表cis调控区的突变；Trans代表反式作用转录因子的改变；“?”代表原因不详。
Table 2. Summary of mutations causing flower color changes from blue/purple to red/pink
表2. 花色由蓝/紫色到红/粉色变化的基因突变形式
COD代表编码区的突变；CIS代表cis调控区的突变；Trans代表反式作用转录因子的改变；“?”代表原因不详。
裂叶牵牛(I. nil)、三色牵牛(I. tricolor)而言，其红/粉色变异通常是由于F3’H的编码区突变造成 [59] [60] 。对于同属内蓝/紫色到红/粉色的转变，若蓝/紫色由飞燕草素决定，通常会伴随F3’5’H的功能失活 [78] [80] - [85] [87] 或表达量下降 [53] ，有时也伴随F3’H在花组织中表达量的下降 [78] [80] [81] ；若蓝/紫色由矢车菊素决定，通常红/粉花中F3’H的启动子序列(cis调控)与蓝/紫花相比会发生变异而导致其表达量下调，且组织特异性表达调控使其仅降低花组织而不影响营养组织中矢车菊素的含量 [75] [76] [77] [87] [88] 。这是因为矢车菊素除了贡献于花色外，在紫外防御、抵御病原菌侵染或动物啃食、生长素极性运输等方面还具有重要的作用，所以F3’H的功能失活可能会影响整个植株的生长，组织特异性表达调控可以降低其有害多效性(deleterious pleiotropy)。统计数据表明植物营养组织中的花青素主要为矢车菊素，很少有植物会缺少有功能的F3’H拷贝，而F3’5’H会在许多科属如蔷薇科(Rosaceae)、菊科(Asteraceae)、拟南芥属(Arabidopsis)、番薯属(Ipomoea)、紫罗兰属(Matthiola)、郁金香属(Tulipa)植物中缺失 [45] [87] 。
由于F3’H或F3’5’H分支被阻断，花青素代谢通路只能走F3H-DFR的路径，所以花色从蓝/紫色到红/粉色的变异有时还伴随着DFR底物特异性的改变，即对底物DHK亲和性的提高，对DHQ和DHM亲和性的降低 [77] [78] [79] 。FLS表达量下降导致的黄酮醇含量下降也会使花由蓝紫色变为红色 [86] 。
3.3. 从蓝/紫/红色到白/黄色的花色变异
花色从蓝/紫/红色到白/黄色的变异通常被认为是花青素从有到无的变异，是质的变化。这些变化可能来自花青素代谢途径一个或多个步骤的阻断(表3)。代谢通路上结构基因的转座子插入 [62] [63] [89] [90] [91] 或移码突变 [92] 导致的功能失活或cis突变导致的表达量下降 [95] [96] 均会使植物产生白花。但结构基因的突变通常会使整个植物都无法产生花青素，最上游基因的突变可以让植物连黄酮、黄酮醇、黄烷酮等类黄酮物质也无法产生，会因而产生较多的有害多效性，降低植物适合度(e.g. [93] [94] )。R2R3-MYB [50] [51] [96] - [103] 、bHLH [104] [105] 和WDR [99] 的突变均可以造成花青素代谢途径上结构基因的表达下调而阻断花青素的合成，但由于R2R3-MYB具有较高的组织表达特异性和功能分化 [33] [34] [39] [50] [103] [106] [107] ，所以其突变具有较低的有害多效性。有观点认为自然选择应该偏爱R2R3-MYB的突变 [88] 。这与R2R3-MYB较bHLH和WDR及结构基因具有较高的异义和同义突变比率(ω)是一致的 [98] [108] [109] [110] [111] 。
花青素量的变化可能与花青素代谢途径基因的表达量变化有关，但具体机制还不清楚。在亚洲百合(Asiatic hybrid lilies Lilium spp.)中，MYB12的表达量与花瓣中花青素的积累量呈正比 [112] 。R2R3-MYB的表达量可解释小天蓝绣球(P. drummondii) (深红)与其近缘种P. cuspidate (浅蓝色)杂交F2代群体中花色的深浅 [53] [113] 。
3.4. 花着色模式的变异
花的着色式样主要有断色、彩斑、缘环、星状、紫外模式、斑点、脉色和芽红等(图2)。目前已知的花色着色模式成因主要有四种：病毒感染、转座子的插入及部分回复突变、siRNA介导的基因转录后沉默、基因特别是调控基因的时空特异性表达(表4)。黄水仙花(Narcissus pseudonarcissus)由于水仙花叶病毒的感染会产生断色现象 [61] 。通常转座子在基因内部的插入和部分回复突变会使花色呈彩斑状 [62] [63] [64] [89] [90] [102] 。双色花(包括缘环和星状)的形成通常是由于白色区域CHS的表达下调(主要由siRNA介导的其在转录后mRNA的降解造成) [65] [67] [114] [115] [116] ，或有色区域FLS的表达上调(原因未知) [65] 。UV模式、斑点、脉色和芽红等表型通常与花青素代谢途径的结构或调控基因的时空表达模式和不同拷贝数有关。在圆叶牵牛中，CHS-D主要在花檐中表达而控制整个花檐的着色，而CHS-E主要在花筒、花脉和花药中表达，可能与花筒的脉色和花脉上有色斑点的形成有关 [117] [118] [119] 。黑心金光菊(Rudbeckia hirta)在花瓣基部有F3H、F3’H、FLS、F6H、F7GT的高量表达使UV吸收类的黄酮醇类物质
Table 3. Summary of known mutations for flower color changes from blue/purple/red to white/yellow
表3. 花色由蓝/紫/红色到白/黄色变化的基因突变形式
COD代表编码区的突变；CIS代表cis调控区的突变；Trans代表反式作用转录因子的改变；“?”代表原因不详。
Table 4. Summary of known mutations that cause pigmentation patterns
表4. 已报道的花色着色模式相关变异
COD代表编码区的突变；CIS代表cis调控区的突变；Trans代表反式作用转录因子的改变；“?”代表原因不详。
高积累而形成类似蜜导的UV模式 [68] ；F3’H1-A在整个花瓣和DFR2-A/B在斑点的早表达及F3’5’H1-A/B和DFR1-A/B在花瓣中的晚表达可能是形成细长山字草(Clarkia gracilis)粉色带红色斑点表型的主要原因，但决定斑点位置的因子还有待确定 [72] [120] 。同样，花青素代谢途径结构基因在牡丹(Paeonia suffruticosa Cultivar “Jinrong”) [121] 和三色堇(Viola × wittrockiana Gams.) [122] 斑点位置出现高表达，这对它们花瓣上形成斑点具有很重要的作用。因为在斑点的形成过程中通常涉及多个结构基因的协同表达，且结构基因的表达受调控基因的影响，所以结构基因的时空特异性表达很大程度上可能是由调控基因的时空表达特异性决定的。许多研究表明R2R3-MYB的时空特异性表达对不同部位的着色具有非常重要的影响(表4) [50] [69] [71] [73] [74] [101] [103] [106] [107] [112] [123] [124] [125] [126] 。周围有光环的斑点表型(图2z)可能与斑点周围激活型因子(WDR)的移除或抑制型因子(R3-MYB)的流入有关 [68] 。R3-MYB和WD40均具有细胞间可移动性 [69] [106] 。在斑点处MBW复合物可促进R3-MYB的产生 [127] [128] ，斑点周围区域R3-MYB的流入与WDR的流出使其花青素合成较少或无，故呈光环型斑点的表型。另外，细胞形状 [15] [17] [69] 和Fe3+转运和存储蛋白的组织特异性表达 [26] [27] [129] 也可能影响花色斑点的形成。
目前MYB的时空特异性表达机制并不清楚。影响营养组织和芽红着色的R2R3-MYB通常会受光的诱导而高表达，且抑制型的R3-MYB会被激活参与该过程以避免过激反应 [106] [127] ，但环境特别是光诱导花青素合成的分子机制并不清楚。R2R3-MYB启动子区甲基化程度的不同会决定其表达与否而影响花青素合成部位，但DNA甲基化的组织特异性机制也未知 [125] 。
4. 花色表型变异的自然选择机制
花色的变异可能会带来传粉者或其访问频率的变异，传粉者与花色的相互关系对花色进化具有重要的意义。长期以来，人们对传粉者在花色亿万年的选择中发挥的作用十分好奇。此外，非传粉者因素(自然环境、病原菌、啃食者、基因多效性等)也会影响花色多态性 [132] [180] [181] 。
4.1. 花色转换方向与传粉者偏好
在陆地植物中发生过多次独立的花色从蓝/紫色到红/粉色或花色从蓝/紫/红色到白/黄色的转换，如曼陀罗属(Iochroma)、番薯属(Ipomoea)、钓钟柳属(Penstemon)、耧斗菜属(Aquilegia) (图3a-d)、闭鞘姜属(Costus) [45] [75] - [81] [87] [92] [130] [131] [132] 。在这些属中蓝紫色常为原始色，花色从蓝/紫色到红/粉色和蓝/紫/红色到白/黄色的转换比较常见，而反向转换则较少 [92] [132] [133] 。已记录的反向转换包括在苦苣苔科(Gesneriaceae)的三个属(Sinningia，Paliavana和Vanhouttea)中曾发生过多次花色从红色到蓝色的转换 [134] ，在黄蕊花属(Dalechampia)和枫属(Acer)中曾分别发生过三次花色从绿色或黄色到红色或紫色的转换 [135] ，在沟酸浆属(Mimulus)中也发生过花色从黄色到红色或粉色的转换(图3e) [52] 。花色转换方向的不对称性对花色的进化方向具有很大的限制 [132] 。目前对于花色着色模式的进化转换研究较少。仅有研究认为金鱼草的脉色可能为祖先性状，在全红的花瓣背景下脉色的作用被红色掩盖进而导致脉色消失 [50] [124] 。
通常，花色从蓝/紫色到红/粉色的转换伴随传粉者从蜂类(bees)到蜂鸟(hummingbirds)的转换，从蓝/紫/红色到白/黄色的转换伴随传粉者从蜂类或蜂鸟到蛾类(moths)的转换 [45] [87] [132] [136] 。花色的质变导致的传粉者类型的改变在许多科属如番薯属 [75] 、沟酸浆属 [137] [138] [139] 、耧斗菜属 [100] 、碧冬茄属(Petunia) [140] 、假番薯属(Ipomopsis) [141] [142] 中均观察到。通常为了避免对同种传粉者的竞争，同域分布的两近缘种(异域分布时通常具有相同花色)会分化出不同的花色 [53] [113] [143] [159] [160] [161] [162] 。传粉者对花色的量变是否具有区分能力还存在许多不确定性 [144] [145] [146] 。着色模式的变化通常会影响传粉者的访问频率，具有较明显斑点或花蜜指示标志的花色类型 [147] - [152] 会受到偏爱。
4.2. 传粉者介导的自然选择对花色进化的影响
传粉者对花色表型变异的区别访问是否对花色的进化产生选择性作用，通常要看该区别访问是否造成不同花色类型植物适合度(fitness)的不同 [132] 。植物适合度组分(fitness components)包括结果率、结籽率、萌发率、存活率、可育率等生活史性状。适合度反映一个个体能够产生的有活力的后代数量。
多数研究缺乏对不同花色变异的适合度的评估(e.g. [75] [100] [124] [137] [138] [140] )。一些研究只测量了部分适合度成分但未检测到传粉者对花色的选择作用(e.g. [113] [143] [154] [156] )。还有的研究注意到不同花色类型间适合度成分有差异且传粉者有偏爱，但无确切的实验和分析证据表明该差异是由传粉者造成(e.g. [72] [139] [152] [164] [166] [169] [175] [191] )。但相当数量的研究观测到传粉者的区别访问会带来不同花色类型植物适合度组分的不同(表5)。不同传粉者对不同花色类型的偏爱可能会使得同一属不同物种或同种物种不同个体间基因流减少，导致花色的歧化选择(disruptive selection) [141] [142] [157] [158] 。在同域地区为了避免对同种传粉者的竞争，花色也会发生歧化 [159] 。同种传粉者对不同花色类型访问次数及取食时间的差异导致的花粉输出与沉积量的差异对花色多样性的维持具有重要作用 [153] [167] 。传粉者对花色某一类型的偏爱性可导致花色的方向性选择(directional selection) [145] [146] [151] [165] [168] [170] [172] 或稳定性选择(stabilizing selection) [163] 。有时，某一花色类型在受到传粉者低访问时自交率的增加会维持其基因频率在一定水平，对花色产生平衡选择(balancing selection) [168] [169] [170] [177] [178] [179] ，从而维持花色多态性。
Table 5. Effects of flower color on pollinator visits and fitness components
表5. 花色对传粉者和适合度组成分的影响
“比较的花色或着色模式”中粗体部分表示该花色类型适合度较高；“适合度成分”中粗体部分表示有显著差异的适合度成分；“-”表示无实验证据。
在进化过程中，植物为适应不同的传粉者会使花部特征特异化。如被蜂鸟传粉的花通常为红色、具有较长的花筒、较窄的花檐、外突的雄蕊和柱头、分泌大量但糖浓度较低的花蜜，被蜂类传粉的花大多为蓝紫色、具较宽的花筒和花檐、内嵌的雄蕊和柱头、分泌少量但高浓度的花蜜，而被蛾类传粉的花则多为白色、花筒较长、有芳香气味且经常在夜间开放，具低浓度高量花蜜 [132] [136] [182] [183] 。花色与其他花性状之间具有高度的相关性(e.g. [184] [185] [186] )，对其他花性状的选择也可能会影响花色的变异。性状之间的相互关联和权衡是限制花色方向性进化的一个因素 [187] 。
4.3. 非传粉者因素介导的自然选择对花色进化的影响
非传粉者因素对不同花色类型的适合度也具有一定的影响(表5)。因为不同的花色类型对光照、降水量、温度、海拔、经纬度等环境的适应能力可能不同，所以环境的异质性可导致不同花色类型在不同地区分布的差异 [152] [157] [158] [159] [168] [170] [171] [190] ，可能会介导花色的歧化或动态选择(fluctuating selection)。啃食者与传粉者偏爱的相同可能带来与传粉者作用相反的拮抗选择(antagonistic selection) [152] [165] 。环境因子、啃食者或病原菌在花色的多态性维持方面起着非常重要的作用，其通常还会限制花色的方向性进化 [165] [168] [170] [171] [172] 。
基因的多效性也会影响花色。有时，杂合体花色表型会表现出超显性(overdominance) [165] [173] [174] [175] [176] 而维持等位基因频率的平衡，对花色产生平衡选择。由于突变基因的有害多效性，一些白花个体与有色花相比，通常具有较低的萌发率、存活率、可育率或抵御啃食者和病原菌的能力较差，因而具有较低的适合度而受到纯化选择(purifying selection)，结果在群体中占有较低的频率 [93] [94] [165] [188] [189] 。
5. 展望
花色表型的丰富变异既存在于种间，也存在于种内。花色表型变异的分子机制和自然选择的研究对探究生物多样性的存在机制具有重要意义。目前对种内的数量变异研究相对较少。传粉者或非传粉者因素对花色变异所起的选择作用有可能影响花色的进化，但其作用机制还需要更多研究。已有的研究结果多数还需要生态学的调查加以确认。花色表型变异导致的植物-传粉者关系的变化如何影响物种的分化，进而影响生态系统的动态变化尚需更多关注。特别是随着气候变化的加剧，自然环境对植物性状和传粉者访花行为的影响将越来越受重视。显然，相关研究将会深化我们对未来气候下动植物互作关系的认识。
基金项目
国家自然科学基金：91331116。