-
从生理机制上探讨植物与环境的关系是植物生理生态学的核心内容[1]。在气候变化的大背景下,通过改变植物的生长环境来对比植物生理性状的差异是目前植物生理生态的研究热点[2]。对于许多分布较广的植物而言,其生理及其他性状都有可能随地理纬度变化表现出适应性变异[3-4]。植物的光合生理性状可反映植物对光及CO2的利用特点,其中的叶绿素含量以及叶绿素荧光参数可以作为评价植物对环境的生长适应能力的主要生理指标[5-6]。种源试验不仅是研究植物地理变异(由环境与遗传因素主导)对植物性状作用机制的重要方法,也是预测在气候变化背景下物种分布范围变迁的重要研究手段[7-8]。
栓皮栎(Quercus variabilis Bl.)是东亚天然地理分布最广泛的树种之一,其分布区域包括中国大陆、台湾和海南岛、日本列岛以及朝鲜半岛。栓皮栎在我国分布可北至辽宁,南达广西、台湾等地,是我国暖温带、亚热带落叶阔叶林中具有代表性的树种之一[9-10],也是我国最重要的造林树种之一[11-12]。因而,栓皮栎被认为是在大区域尺度上研究植物对环境适应机理的理想树种之一[13-14]。
本研究选择栓皮栎天然分布的5个不同纬度的种源,并在其分布南界以南(广东鹤山,HS)进行田间栽培试验。通过测定不同种源栓皮栎幼苗的光合荧光参数,分析各种源南移后的光合和荧光生理特性及其对新生境下的适应性差异,为将来气候变化情景下栓皮栎的引种驯化以及理解植物对水平梯度(纬度)带来的环境差异的响应方式和适应策略提供科学参考。
HTML
-
2014年9—10月,选择我国5个不同纬度的栓皮栎天然成熟林分收集种子,分别为北京平谷(PG)、河南内乡(NX)、湖北秭归(ZG)、江西永修(YX)及湖南城步(CB)。2015年2月初,将种子播种于广东鹤山(HS)林科所苗圃(112°54′ E,22°41′ N,海拔56 m)。2015年5月,选择健康、大小均一的不同种源栓皮栎幼苗分别移植至广东鹤山试验地进行观察,随机排列各种源幼苗的栽植位置。常规管理,至第2年生长旺季(2016年7月)对每个种源随机选取10株幼苗,分别测量其地径、株高,随机选择3~4株健康植株作为重复,编号,选择每株幼苗中上部有代表性的健康、成熟的叶片进行光合生理指标的测定(包括叶绿素含量、光合及叶绿素荧光的测定)。另外,收集种源地临近气象台站的气象数据(1950—2000年),使用Kriging空间插值法(考虑海拔对温度的影响)进行拟合,根据其经纬度坐标提取所需气象数据(表 1)。
种源地/试验地
Sites纬度
NLatitude经度
ELongitude海拔/m
Elevation年均气温/℃
MAT年降水量/mm
MAP气候区
Climatic zone北京平谷 PG 40°16′ 117°07′ ~300 10.3 556 暖温带 河南内乡 NX 33°29′ 111°54′ ~900 10.9 851 暖温带 湖北秭归 ZG 30°46′ 110°20′ ~950 12.6 1 241 北亚热带 江西永修 YX 29°05′ 115°35′ ~820 15.8 1 290 中亚热带 湖南城步 CB 26°18′ 110°07′ ~1 050 13.2 1 540 中亚热带 广东鹤山 HS 22°41′ 112°54′ ~50 22.4 1 462 南亚热带 Table 1. Geographical location and climate conditions of the provenance and cultivated sites
-
图 1表明:在低光合有效辐射范围内(PAR<450 μmol·m-2·s-1),不同生境下各种源栓皮栎叶片的Pn对光照强度的响应均较敏感,即随着PAR的增加,Pn也随之迅速增加,且各个种源的增长趋势基本一致;当PAR>450 μmol·m-2·s-1之后,各种源栓皮栎在同一PAR下的Pn值大小顺序为CB>NX>YX>PG>ZG,其净光合速率缓慢地逐渐增加,没表现出明显的光饱和现象, 其中,PG种源的净光合速率随光强的增加表现出相对平缓的趋势。
光响应曲线参数拟合结果(表 2)表明:YX种源栓皮栎的表观量子效率最高,显著高于ZG种源,而与其它3个种源差异不显著。CB种源的最大净光合速率(Pmax)显著高于PG种源和ZG种源,而NX种源和YX种源的Pmax与上述3个种源差异不显著。不同种源间的暗呼吸速率、光饱和点和光补偿点的差异不显著(P>0.05),但暗呼吸速率和光补偿点有随高纬度(北方)种源向低纬度(南方)种源逐级递减的大致趋势,而低纬度种源CB的光饱和点最高。
指标 Parameters 北京平谷 PG 河南内乡 NX 湖北秭归 ZG 江西永修 YX 湖南城步 CB AQY/(μmol·μmol-1) 0.059±0.003ab 0.057±0.002ab 0.045±0.009b 0.069±0.002a 0.055±0.009ab Pmax/(μmol·m-2·s-1) 19.938±0.787b 25.934±4.035ab 17.716±0.347b 22.089±3.255ab 29.653±4.446a Rd/(μmol·m-2·s-1) 3.536±0.732a 3.535±0.223a 2.659±0.298a 3.179±1.439a 2.410±0.448a LSP/(μmol·m-2·s-1) 582.124±64.83a 697.669±91.75a 591.42±105.571a 563.932±149.35a 718.241±62.139a LCP/(μmol·m-2·s-1) 78.145±24.571a 74.99±5.329a 71.474±23.884a 59.17±36.011a 47.701±4.328a 注 Notes:AQY:表观量子效率 Apparent quantum yield;Pmax:最大净光合速率 Maximum net photosynthetic rate;Rd:暗呼吸速率 Dark respiration rate;LSP:光饱和点 Light saturation point;LCP:光补偿点 Light compensation point。数据为平均值±标准偏差,同一行数值带有不同小写字母表示在0.05水平上差异显著(P<0.05);下同 Table 2. Fitted photosynthetic parameters of Q. variabilis leaves among various provenances
-
采用FvCB模型模拟的5个栓皮栎种源叶片净光合速率对CO2响应的拟合曲线(Pn-Ci)(图 2)表明:在1 200 μmol·m-2·s-1光强下,各种源Pn均随胞间CO2浓度(Ci)的增加而增加,最后逐渐平缓并达到饱和。FvCB模型模拟的Pn-Ci曲线与光响应曲线不同,具有明显的分段性。虚线之前为第一阶段,细线为CO2限制阶段的净光合速率(Ac);虚线之后为第二阶段,细线为最大电子传递速率(Jmax)限制阶段的净光合速率(Aj)。粗线代表 2个阶段的限制速率。第一阶段的Pn随Ci的增加而呈线性增加,此时光合作用受Rubisco酶活性大小的限制;第二阶段的Pn随Ci的增加呈逐渐平缓的曲线增加,此时光合作用受RuBP再生速率大小的限制。虚线处Ci值则为这2个限制阶段的临界胞间CO2浓度(Ac_Aj)。
Figure 2. Photosynthesis CO2 response curves(Pn-Ci) of Q. variabilis leaves among various provenances
从表 3看出:栓皮栎幼苗的Rubisco酶最大羧化速率(Vcmax)随纬度的升高而逐渐降低(CB>YX>ZG>NX>PG),但各种源间差异不显著(P>0.05),这表明不同种源栓皮栎叶片叶肉细胞Rubisco酶的羧化能力相差不大。另外,最大电子传递速率(Jmax)、光呼吸速率(Rp)、Ac_Aj和CO2补偿点(CCP)在不同种源间的差异均不显著,但ZG种源的Ac_Aj明显低于其他4个种源。
指标 Parameters 北京平谷 PG 河南内乡 NX 湖北秭归 ZG 江西永修 YX 湖南城步 CB Vcmax/(μmol·m-2·s-1) 39.87±17.10a 43.15±11.33a 48.96±12.72a 49.60±0.20a 58.01±6.10a Jmax/(μmol·m-2·s-1) 173.29±2.15a 190.95±77.41a 226.59±58.08a 221.37±31.11a 212.86±14.71a Rp/(μmol·m-2·s-1) 1.170±0.360a 1.020±0.080a 1.140±0.410a 1.140±0.490a 1.130±0.060a Ac_Aj/(μmol·mol-1) 434.46±200.64a 404.17±17.73a 290.06±74.50a 371.95±20.97a 374.47±57.64a CCP/(μmol·mol-1) 72.922±12.905a 66.343±4.467a 74.541±11.537a 74.110±1.774a 65.557±4.513a 注 Notes: Vcmax:Rubisco 酶最大羧化速率 Rubisco enzyme maximum carboxylation rate;Jmax:最大电子传递速率 Maximum electron transport rate;Rp:线粒体在光照下的呼吸速率 Respiration rate of mitochondria under light;Ac_Aj:Rubico 酶限制阶段与 RuBP 限制阶段临界点的胞间CO2浓度 Intercellular CO2 concentration of the transition poin tbetween Rubisco carboxylation-limited and RuBP regeneration-limited state;CCP:CO2补偿点 CO2 compensation point. Table 3. Fitted parameters from Pn-Ci curves of Q. variabilis leaves among various provenances
-
通过对5个种源栓皮栎的叶绿素荧光参数进行统计分析,结果(表 4)表明:光系统Ⅱ(PSⅡ)原初光能转换效率(Fv/Fm)、PSⅡ潜在活性(Fv/F0)和PSⅡ实际光化学反应效率ΦPSⅡ参数在不同种源间差异不显著(P>0.05),且不同种源栓皮栎在非光化学猝灭系数(NPQ)、光化学猝灭系数(qP)和相对电子传递速率(ETR)参数上差异也较小(P>0.05)。叶绿素相对含量(SPAD值)方面,NX种源和CB种源显著高于YX种源(P<0.05),而与其他2个种源的SPAD值差异不显著。表 4还表明:不同种源株高、地径之间差异不显著,但PG种源的地径、株高较大,生长较快。
指标 Parameters 北京平谷 PG 河南内乡 NX 湖北秭归 ZG 江西永修 YX 湖南城步 CB D(mm) 13.50±0.890a 10.77±3.745a 11.61±3.337a 12.25±2.865a 12.51±2.094a H(cm) 112.0±42.48a 80.00±14.54a 77.75±20.19a 90.00±26.24a 93.25±23.37a SPAD 48.00±2.50ab 50.73±4.84a 47.56±1.75ab 44.25±2.28b 49.05±2.76a ΦPSⅡ 0.228±0.034a 0.221±0.027a 0.228±0.018a 0.227±0.026a 0.223±0.009a ETR 78.25±11.85a 76.20±9.18a 78.45±6.17a 78.08±8.95a 76.60±3.09a qP 0.437±0.040a 0.410±0.070a 0.423±0.047a 0.416±0.045a 0.420±0.017a NPQ 2.738±0.519a 2.347±0.521a 2.583±0.364a 2.595±0.364a 2.642±0.165a Fv/Fm 0.800±0.029a 0.797±0.038a 0.807±0.033a 0.810±0.031a 0.805±0.007a Fv/F0 4.071±0.727a 4.057±0.948a 4.276±0.833a 4.365±0.771a 4.133±0.188a 注 Notes: D:地径 Basal diameter;H:高度 Height; SPAD 值:相对叶绿素含量;ΦPSⅡ:PSⅡ 实际光化学反应效率 Effective quantum yield of photosystem Ⅱ photochemistry;ETR:相对电子传递速率 Rel. electron transport rate;qP:光化学猝灭 Photochemical quenching;NPQ:非光化学猝灭 Non-photochemical quenching;Fv/Fm:PSⅡ原初光能转换效率 Primary conversion of light energy of photosystem Ⅱ;Fv/F0:PSⅡ潜在活性 Potential activation ofphotosystem Ⅱ;重复4次。 Table 4. Chlorophyll content and chlorophyll fluorescence parameters of Q.variabilis seedlings among different provenances