文章信息
- 张珊珊, 康洪梅, 杨文忠, 向振勇
- ZHANG Shangshan, KANG Hongmei, YANG Wenzhong, XIANG Zhenyong.
- 干旱胁迫下AMF对云南蓝果树幼苗生长和光合特征的影响
- Effects of arbuscular mycorrhizal fungi on growth and photosynthetic characteristics of Nyssa yunnanensis seedlings under drought stress
- 生态学报[J]. 2016, 36(21): 6850-6862
- Acta Ecologica Sinica[J]. 2016, 36(21): 6850-6862
- http://dx.doi.org/10.5846/stxb201504270862
-
文章历史
- 收稿日期: 2015-04-27
- 网络出版日期: 2016-03-03
云南蓝果树(Nyssa yunnanensis W.C.Yin)为蓝果树科(Nyssaceae)蓝果树属(Nyssa Gronov ex Linn.)植物, 国级Ⅰ级重点保护野生植物, 云南特有种和极危种[1]。云南蓝果树天然林内幼苗更新极其困难[1], 但对云南蓝果树的研究仅见于系统分类、形态修订及种子萌发特性等方面的研究[2-3], 未见相关报道。
更新幼苗与降水特点密切相关。关于木本植物的研究发现, 干旱期的长短很大程度上影响了幼苗的生长和存活, 连续的降水过程或连阴雨天气才是更新幼苗发生的有利条件[4-5]。另外, 土壤水分是幼苗生长存活的关键因子, 在光照强度大的森林中, 由于蒸发量大, 土壤含水量较低, 森林内的幼苗在冬季干旱时容易大量死亡[6]。而云南蓝果树的原生境所在地西双版纳地区38 a(1974—2003年)来平均的气候情况为:从9 月份至次年2 月份, 西双版纳均处于旱季[7-8]。而云南蓝果树种子一般于每年的8 月份之前就基本成熟, 错过了种子萌发的最适季节。这说明干旱已经严重影响了云南蓝果树的种群延续和扩大。然而, 干旱胁迫究竟对云南蓝果树幼苗有哪些影响?云南蓝果树能否表现出相应的表型可塑性?影响这种表型可塑性的因素又有哪些?这都需进一步研究。
丛枝菌根真菌(arbuscular mycorrhizal fungi, AMF)能改善植物养分吸收[9]、植株生长状况[10-11]等。研究表明, 丛枝菌根能显著促进植物的生长[12]、提高植物养分吸收的能力[13]、增强植物的抗逆性(例如:抗干旱、耐瘠薄和耐重金属污染等)和抗病性[14]。何跃军等观察到接种AM真菌的构树(Broussonetia papyrifera)生长量及净光合速率都显著高于未接种的对照株[15]。虽然, 目前还未见AMF与云南蓝果树共生关系的相关报道, 但是笔者在对云南蓝果树幼苗的根进行解剖时发现其根部可被AMF侵染, 且侵染率较高(45.66%)(文章未发表)。那么, 干旱胁迫条件下AMF对云南蓝果树的影响如何, 以及菌根共生体功能的发挥与程度又如何呢?
本试验选用苯菌灵为杀真菌剂, 用盆栽方法人工模拟土壤干旱条件, 研究正常供水和干旱胁迫条件下AMF对云南蓝果树幼苗生长和光合特征的影响, 进而揭示云南蓝果树濒危原因, 为其科学保育及种苗繁育提供理论依据。苯菌灵能很好地抑制AMF的活性, 降低其与植物形成共生体的机会, 抑制丛枝菌根的形成。
1 材料和方法 1.1 材料选择云南蓝果树1年实生幼苗为研究对象。2013年4月中旬将生长基本一致的幼苗移栽到容积为10 L的花盆中, 每盆1株, 栽培基质为云南省林业科学院苗圃红壤。试验用红壤于次氯酸钠(NaClO)中浸泡48 h, 后洗净后120℃烘8 h。
1.2 试验设计试验于云南省林业科学院温室条件下进行, 自然温度和自然光照, 水分人工控制, 此温室主要种兰花用, 需对其进行定期加湿, 因此整个实验过程中的空气湿度较高。试验为2因子(AMF×水分)试验。AMF处理为施加苯菌灵(低AMF)和不加苯菌灵(高AMF), 对于施加苯菌灵处理, 将2 g杀真菌剂苯菌灵溶于2 L自来水, 加到盆钵中, 每月处理1次, 获得低AMF处理的土壤, 同时, 在对照处理的盆钵中每次均加入相同量的自来水。水分设置6个水平, 共得到2(AMF)×6(水分)=12个处理, 每个处理10个重复, 共有120个盆钵。正式实验前, 测得栽培用土的田间持水能力(WHC), 此时的实际含水率为32.32%(干土重基础)。实验处理以相对土壤含水量(RWC), 即各盆的土壤含水率占WHC的百分比来区分(表 1):第1个处理始终保持饱和, 每天补足由于蒸散所散失的水分(100%);第2个处理保持在90%—100%(91.68%), 当由于蒸散所散失的水分达到90%时, 立即补充水分, 使之达到饱和;以同样的方法保持第3个处理在80%—100%(82.85%);第4个处理在60%—100%(60.00%);第5个处理在40%—100%(41.86%);第6个处理在20%—100%(21.28%)。文中分别用W1, W2, W3, W4, W5和W6表示6个处理。浇水时用缓慢滴加的方法, 尽量避免水分从花盆底部流失, 这样可以最大限度地保证土壤养分不损失。浇水在每天下午进行, 天气炎热时, W1处理每天中午补水100—200 mL。实验期限为6个月, 2014年10月中旬采样, 测定相应指标。
参数Parameters | 有无AMF处理(0, 低AMF;高AMF)With or without AMF(0, low AMF; 1, high AMF) | W1 | W2 | W3 | W4 | W5 | W6 | 可塑性指数Plasticity index |
叶总面积/cm2 | AMF1 | 111.66±3.90aA | 69.47±3.16bA | 50.96±4.71bA | 36.81±3.93bcA | 23.15±4.47cdA | 12.95±1.72dA | 0.88 |
Total leave area | AMF0 | 68.06±4.94aB | 34.62±4.68bcB | 19.60±3.29cdB | 23.54±3.53cdA | 11.45±1.51dA | 10.04±0.92dA | 0.85 |
叶片数目 | AMF1 | 20.60±2.38aA | 20.20±2.13abA | 16.25±2.56abcA | 11.20±1.58cdA | 6.75±1.49dA | 7.33±1.07dA | 0.67 |
Number of leaves | AMF0 | 14.60±1.61aB | 8.80±1.89bB | 7.50±1.76bB | 8.00±1.84bA | 5.50±1.32bA | 5.50±1.02bA | 0.62 |
株高/cm | AMF1 | 51.22±2.17aA | 36.18±2.58bA | 30.25±2.14bcA | 27.75±2.28cdA | 23.66±2.47cdA | 22.63±1.95dA | 0.56 |
Height | AMF0 | 33.15±2.62aB | 23.81±2.57bB | 21.18±1.99bB | 20.85±2.56bA | 20.85±2.45bA | 18.96±2.39bA | 0.43 |
相对生长速率 | AMF1 | 17.79±2.29aA | 14.64±3.55aA | 13.45±3.77aA | 10.05±2.28aA | 9.19±2.28aA | 6.48±1.08aA | 0.63 |
Relative growth ration | AMF0 | 12.77±3.07aB | 11.21±2.89aA | 10.95±1.69aA | 10. 72±2.09aA | 7.16±1.69aA | 4.93±1.36aA | 0.61 |
总生物量/(g/株) | AMF1 | 11.66±1.05aA | 6.95±1.03bA | 4.64±1.02cA | 2.49±0.59cdA | 1.54±0.51cdA | 1.03±0.27dA | 0.91 |
Total biomass | AMF0 | 6.09±1.04aB | 2.00±0.98bB | 1.26±0.95bB | 1.12±0.36bB | 0.90±0.08bB | 0.96±0.06bA | 0.85 |
地上部分生物量 /g | AMF1 | 7.61±1.03aA | 4.74±1.00bA | 3.16±1.00cA | 2.07±0.39cdA | 0.76±0.24deA | 0.60±0.21eA | 0.92 |
Aboveground biomass | AMF0 | 4.14±1.01aB | 1.36±0.95bB | 0.92±0.30bcB | 0.77±0.16bcB | 0.54±0.06cA | 0.51±0.08cA | 0.88 |
地下部分生物量/g | AMF1 | 4.05±0.69aA | 2.21±0.03bA | 1.48±0.41bcA | 0.84±0.03cdA | 1.04±0.17cdA | 0.43±0.07dA | 0.89 |
Belowground biomass | AMF0 | 1.95±0.25aB | 0.64±0.023bB | 0.34±0.03bB | 0.34±0.02bB | 0.35±0.01bB | 0.41±0.07bA | 0.83 |
地上部分生物量比 | AMF1 | 0.69±0.11abA | 0.73±0.12aA | 0.72±0.14aA | 0.69±0.05abA | 0.61±0.03abA | 0.57±0.07bA | 0.22 |
Aboveground biomass/ Total biomass | AMF0 | 0.65±0.04aA | 0.68±0.06aA | 0.67±0.07aA | 0.56±0.33bA | 0.56±0.31bB | 0.57±0.09abA | 0.18 |
地下部分生物量比 | AMF1 | 0.35±0.04bcA | 0.27±0.01cA | 0.33±0.016bcA | 0.44±0.022bA | 0.60±0.03aA | 0.43±0.008bA | 0.41 |
Belowground biomass/ Total biomass | AMF0 | 0.31±0.03bcA | 0.32±0.03bcA | 0.28±0.013cA | 0.31±0.01bcA | 0.39±0.014abA | 0.43±0.007aA | 0.28 |
根冠比Aboveground biomass/ | AMF1 | 0.55±0.11aA | 0.40±0.005aA | 0.52±0.004aA | 0.46±0.003aA | 0.53±0.006aA | 0.51±0.002aA | 0.28 |
Belowground biomass ratio | AMF0 | 0.47±0.125aA | 0.47±0.003aA | 0.43±0.001aA | 0.46±0.006aA | 0.42±0.001aA | 0.46±0.001aA | 0.24 |
叶面积比率 | AMF1 | 197.21±18.04abA | 212.88±17.56aA | 178.62±27.10abA | 165.45±19.07bA | 101.81±22.81cA | 63.23±16.17dA | 0.70 |
Leaf area ratio, LAR | AMF0 | 163.09±12.50aB | 152.52±15.01aB | 117.10±18.65bB | 168.46±14.19aA | 70.35±12.50cB | 57.50±1.00cA | 0.58 |
比叶面积 | AMF1 | 302.27±29.21aA | 295.93±25.63aA | 262.38±28.63aA | 199.18±17.56bB | 206.30±14.05bA | 157.77±12.06cA | 0.48 |
Specific leaf area, SLA | AMF0 | 239.89±22.55aB | 223.78±14.53aA | 159.75±6.08bB | 243.34±15.01aA | 117.13±18.68cB | 107.71±15.01cB | 0.34 |
表中同行不同小写字母表示不同水分条件下叶片的解剖结构指标在P<0.05水平差异显著;同列不同大写字母表示不同AMF处理下叶片的解剖结构指标在P<0.05水平差异显著 |
取部分新鲜根部样品固定于FAA溶液中(37%甲醛-冰醋酸-50%乙醇溶液, 体积比 9:0.5:0.5) 用于检测AMF侵染率。先将根部的固定液清洗干净, 然后浸泡在10%的KOH中, 90℃水浴加热5 min, 然后用1%的盐酸酸化15 min, 并用酸性品红染色过夜, 将根部剪成 2 cm长的根段, 在显微镜下10倍物镜观察, 用十字交叉法计算侵染率(gridline intersect method)(Giovannetti & Mosse, 1980) 。侵染率的计算公式为:侵染率= 侵染根段长度/根段总长度×100%。
1.3.2 生长指标测定实验结束后, 用直尺测量株高, 用数显游标卡尺测量叶长和叶宽, 并记录叶片数目。用CI-203 型便携式激光叶面积仪(CID Inc, Wash ington, USA)测定野外叶片的面积, 得到叶面积与叶长和叶宽的拟合公式:叶面积= 0.7416×叶长×叶宽-0.7477 (R2 = 0.984, P < 0.001, n=120) , 形态测量的叶面积通过拟合公式得出。采用King 等人的方法计算株高的相对生长速率(RGR) =(lnH2-lnH1)/Δt, 式中:H1和H2 为前后两次测量时的株高, Δt为测量间隔时间。
将幼苗整棵挖出清洗干净, 分为地上部分和地下部分, 105℃杀青30 min后在80℃下烘干至恒重。称量各部分干重并计算如下参数:地上部分生物量比(AMR, 地上部分重/总生物量)、地下部分生物量比(BMR, 地下部分重/总生物量)、根冠比(R/S, 地下部分生物量/地上生物量), 结合形态测量的总叶面积计算叶面积比率(LAR, 总叶面积/总生物量)和比叶面积(SLA, 总叶面积/叶重)。
1.3.3 气体交换和资源利用参数测定实验中后期用用LI-6400XT便携式光合仪于09:00—11:00测定净光合速率(Pn)、蒸腾速率(Tr)、气孔导度(Gs)、胞间二氧化碳浓度(Ci)等参数, 并记录光合有效辐射(PAR)和大气CO2 浓度(Ca)等环境参数, 并计算气孔限制值(Ls) = 1-Ci /Ca、瞬时水分利用效率(WUE) = Pn /Tr。
1.4 数据统计分析试验采用SPSS16.0对数据进行单因素方差分析、Duncan多重比较、变量间的相关性分析以及主成分分析。方差分析时, 不满足方差齐性检验的数据通过[arcsin]或[log(x+1) ]转换以满足方差分析的要求。通常采用Post-hoc Tukey方法检验变量的显著性, 如果数据不满足参数检验的条件, 采用Kruskall-Wallis方法检验。5%为显著水平, 1%为极显著水平。根据Valladares等的方法计算各参数的可塑性指数, 即6个水分梯度下某变量的最大值减去最小值再除以最大值。其中隶属函数值具体公式为:如果指标与抗旱性成正相关, X(μ)=(X-Xmin)/(Xmax-Xmin);如果指标与抗旱性成负相关, X(μ)=1-(X-Xmin)/(Xmax-Xmin)。式中, X表示各指标的测定值, Xmax和Xmin分别表示各处理条件下指标的最大和最小测定值。
2 结果与分析 2.1 丛枝菌根真菌侵染率图 1显示, 施加苯菌灵处理显著降低了AMF对云南蓝果树的侵染率, 从而形成低AMF处理(F=38.141, P<0.01) , 而不施加苯菌灵的处理则形成高AMF处理。不管是高AMF处理还是低AMF处理, AMF侵染率都随着土壤水分含量的降低显著下降(F=11.093, P<0.01;F=37.175, P<0.01) 。
2.2 不同水分条件下丛枝菌根真菌对云南蓝果树幼苗生长的影响表 1表明, 随着土壤相对含水量的降低, 不管是AMF侵染率高还是低, 云南蓝果树生长指标中除了相对生长速率和根冠比2个指标上没有发生显著变化, 叶片数目、叶面积、株高、总生物量、地下部分生物量、地上部分生物量、LAR和SLA8个指标呈递减状态。生物量的分配方式也发生了显著的变化(P<0.05) , 地上部分生物量比和地下部分生物量比都呈现先升高后降低的趋势。
在测定的10个生长指标中, AMF处理的影响也不同(表 1)。W1和W2条件下AMF处理对9个生长指标有显著影响;W3条件下AMF处理对8个生长指标有显著影响;W4条件下AMF处理对4个生长指标有显著影响;W5条件下AMF处理对5个生长指标有显著影响;W6条件下AMF处理只对1个生长指标有显著影响。因此, AMF处理从土壤含水量处理为W1时便开始显著影响生长指标的特征值。当土壤过于干旱(W6) 时, AMF处理对其幼苗生长的影响并不显著。
2.3 不同水分条件下丛枝菌根真菌对云南蓝果树幼苗光合特征的影响由表 2看出, 随着土壤相对含水量的降低, 不管是AMF侵染率高还是低, 云南蓝果树光合特征指标中Pn、Tr和Gs 3个指标呈显著递减状态。Ci表现为先降后升的趋势, Ls和WUE则为先升后降的趋势。
参数Parameters | 有无AMF处理(0, 低AMF;高AMF)With or without AMF(0, low AMF; 1, high AMF) | W1 | W2 | W3 | W4 | W5 | W6 | 可塑性指数Plasticity index |
Pn /(μmol m-2 s-1) | AMF1 | 3.71±0.51aA | 2.69±0.07bA | 2.06±0.08cA | 1.29±0.36dA | 0.81±0.11eA | 0.51±0.15eA | 0.86 |
AMF0 | 2.05±0.04aB | 1.67±0.04bB | 0.89±0.26cB | 0.53±0.07dB | 0.47±0.08dB | 0.42±0.13dA | 0.84 | |
Tr/(mmol m-2 s-1) | AMF1 | 1.61±0.04aA | 1.58±0.26aA | 0.84±0.22bA | 0.51±0.04cA | 0.38±0.05cdA | 0.19±0.11dA | 0.88 |
AMF0 | 1.07±0.07aB | 0.92±0.10bB | 0.46±0.06cB | 0.21±0.03dB | 0.22±0.04dB | 0.15±0.07dA | 0.86 | |
Gs/(mol m-2 s-1) | AMF1 | 0.07±0.01aA | 0.06±0.01aA | 0.05±0.01abA | 0.04±0.01bA | 0.03±0.007bcA | 0.02±0.003cA | 0.7 |
AMF0 | 0.05±0.01aA | 0.03±0.01abB | 0.02±0.008bcB | 0.01±0.005bcA | 0.01±0.005bcA | 0.02±0.004cA | 0.61 | |
Ci/(μmol/m2) | AMF1 | 258.00±20.90cA | 221.00±24.74dA | 263.80±39.27cA | 290.40±31.20bcA | 325.20±31.14bA | 372.40±15.87aA | 0.41 |
AMF0 | 203.60±16.49cB | 194.20±36.50cA | 233.00±30.23cA | 248.40±23.14bcB | 316.20±22.87abA | 324.80±13.01aB | 0.37 | |
Ls | AMF1 | 0.86±0.01aA | 0.88±0.07aA | 0.84±0.08aA | 0.83±0.02aA | 0.79±0.02bA | 0.78±0.01bA | 0.12 |
AMF0 | 0.83±0.01abB | 0.87±0.02aA | 0.82±0.03bA | 0.81±0.02bcB | 0.78±0.02cA | 0.75±0.01dB | 0.14 | |
WUE/ (mmol/mol) | AMF1 | 2.29±0.26aA | 3.82±0.05aA | 3.40±0.20aA | 2.38±0.45aA | 1.97±0.23aA | 1.76±0.35aA | 0.54 |
AMF0 | 1.93±0.12bcA | 2.51±0.11aA | 2.15±0.26bA | 1.89±0.38bcA | 1.85±0.22bcB | 1.74±0.24cB | 0.31 |
在测定的10个生长指标中, AMF处理的影响也不同(表 2)。W1条件下AMF处理对4个光合特征有显著影响;W2、W3、W5和W6条件下AMF处理对3个光合特征有显著影响。W4条件下AMF处理对4个光合特征有显著影响;
2.4 云南蓝果树幼苗生长指标和光合特征间的相关性分析由表 3相关分析可以看出:AMF侵染率与叶片数目、叶面积、株高和总生物量等12个指标呈显著正相关;与其他指标间虽然存在一定的或正或负的相关性, 但均未达到显著水平(表 3)。土壤相对含水量与叶片数目、叶面积和总生物量12个指标呈显著正相关。其余18个指标之间存在56个极显著正相关, 12个显著正相关, 5个极显著负相关, 3个显著负相关。
相关性Correlations | 生长指标Growth parameters | 光合指标Photosynthetic characteristics | ||||||||||||||||||
AMF侵染率AMF root colonization rate | 土壤相对含水量Relative-soilwater content | 叶片数目Number of leaves | 叶面积Leave area | 株高Height | 相对生长速率Relative growth ration | 总生物量Total biomass | 地上部分生物量Aboveground biomass | 地下部分生物量Belowground biomass | 地上部分生物量比Aboveground biomass/ Total biomass | 地下部分生物量比Belowground biomass/Total biomass | 根冠比Aboveground biomass/Belowground biomass ratio | 叶面积比率LAR | 比叶面积SLA | Pn | Tr | Gs | Ci | Ls | WUE | |
AMF侵染率AMF root colonization rate | 1 | |||||||||||||||||||
土壤相对含水量Relativesoilwater content | 0.652* | 1 | ||||||||||||||||||
叶片数目Number of leaves | 0.913** | 0.744** | 1 | |||||||||||||||||
叶面积Leave area | 0.976** | 0.688* | 0.935** | 1 | ||||||||||||||||
株高Height | 0.932** | 0.471 | 0.895** | 0.885** | 1 | |||||||||||||||
相对生长速率Relative growth ration | -0.128 | 0.482 | 0.183 | 0.023 | -0.256 | 1 | ||||||||||||||
总生物量Total biomass | 0.948** | 0.772** | 0.975** | 0.947** | 0.883** | 0.127 | 1 | |||||||||||||
地上部分生物量Aboveground biomass | 0.927** | 0.809** | 0.930** | 0.929** | 0.802** | 0.201 | 0.971** | 1 | ||||||||||||
地下部分生物量Belowground biomass | 0.954** | 0.736** | 0.943** | 0.961** | 0.883** | 0.061 | 0.955** | 0.936** | 1 | |||||||||||
地上部分生物量比Aboveground biomass/ Total biomass | -0.398 | -0.857** | -0.505 | -0.488 | -0.130 | -0.670* | -0.543 | -0.632* | -0.470 | 1 | ||||||||||
地下部分生物量比Belowground biomass/ Total biomass | 0.398 | 0.857** | 0.505 | 0.488 | 0.130 | 0.670* | 0.543 | 0.632* | 0.470 | -1.000** | 1 | |||||||||
根冠比Aboveground biomass/Belowground biomass ratio | 0.700* | 0.404 | 0.393 | 0.639* | 0.499 | -0.406 | 0.526 | 0.588* | 0.602* | -0.282 | 0.282 | 1 | ||||||||
叶面积比率LAR | 0.905** | 0.671* | 0.988** | 0.922** | 0.915** | 0.145 | 0.961** | 0.906** | 0.930** | -0.406 | 0.406 | 0.365 | 1 | |||||||
比叶面积SLA | 0.924** | 0.689* | 0.991** | 0.945** | 0.915** | 0.143 | 0.969** | 0.921** | 0.943** | -0.440 | 0.440 | 0.399 | 0.996** | 1 | ||||||
Pn | 0.862** | 0.618* | 0.964** | 0.859** | 0.925** | 0.086 | 0.925** | 0.848** | 0.889** | -0.305 | 0.305 | 0.302 | 0.983** | 0.967** | 1 | |||||
Tr | 0.380 | 0.718** | 0.409 | 0.481 | 0.099 | 0.549 | 0.449 | 0.552 | 0.403 | -0.855** | 0.855** | 0.347 | 0.356 | 0.385 | 0.226 | 1 | ||||
Gs | -0.327 | -0.552 | -0.366 | -0.453 | -0.073 | -0.546 | -0.404 | -0.474 | -0.348 | 0.747** | -0.747** | -0.203 | -0.332 | -0.370 | -0.173 | -0.908** | 1 | |||
Ci | 0.231 | -0.026 | 0.315 | 0.173 | 0.460 | -0.347 | 0.268 | 0.096 | 0.192 | 0.292 | -0.292 | -0.024 | 0.335 | 0.280 | 0.464 | -0.390 | 0.455 | 1 | ||
Ls | 0.807** | 0.793** | 0.791** | 0.846** | 0.643* | 0.142 | 0.800** | 0.804** | 0.817** | -0.732** | 0.732** | 0.594* | 0.708* | 0.747** | 0.631* | 0.625* | -0.538 | 0.031 | 1 | * |
WUE | 0.853** | 0.737** | 0.839** | 0.868** | 0.747** | 0.029 | 0.836** | 0.829** | 0.853** | -0.594* | 0.594* | 0.623* | 0.776** | 0.791** | 0.738** | 0.534 | -0.384 | 0.282 | 0.932** | 1 |
*:0.05水平上相关性显著(双侧检验); **:0.001水平上相关性显著(双侧检验) |
本文采用主成分分析法对18项指标进行分析, 并根据各主成分中每个指标载荷量及其变异系数的大小筛选出具有代表性的指标。
从基于相关矩阵的主成分分析结果(表 4和表 5)可以看出, 低AMF处理条件下, 云南蓝果树幼苗第2个主成分的特征值为3.176, 累计贡献率为87.771%, 说明干旱胁迫下低AMF对云南蓝果树的前2个主因子基本上能概括18个变量的主要信息, 所以共提取了2个主成分因子。在确定了主成分因子数量后, 在第一主成分中, 地下部分生物量、株高和总生物量7个指标得分较高, 均为0.90以上;在第二主成分中, 特征向量系数按照绝对值大小依次排列, 得分最高的地上部分生物量比和地下部分生物量比, 分别是0.908和-0.908。
成分Component | 初始特征值Initial eigenvalue | 因子提取结果Rotation Sums of squared loadings | ||||||||||
合计Total | 方差贡献率/%Proportion | 方差累积率/%Accumulation | 合计Total | 方差贡献率/%Proportion | 方差累积率/%Accumulation/% | |||||||
AMF0 | AMF1 | AMF0 | AMF1 | AMF0 | AMF1 | AMF0 | AMF1 | AMF0 | AMF1 | AMF0 | AMF1 | |
1 | 12.623 | 14.074 | 70.129 | 78.189 | 70.129 | 78.189 | 9.522 | 11.378 | 53.068 | 63.212 | 53.068 | 63.212 |
2 | 3.176 | 2.725 | 17.642 | 15.138 | 87.771 | 93.327 | 6.247 | 5.421 | 34.703 | 30.114 | 87.771 | 93.327 |
3 | 0.982 | 0.536 | 5.454 | 2.977 | 93.226 | 96.304 | ||||||
4 | 0.918 | 0.472 | 5.101 | 2.623 | 98.326 | 98.927 | ||||||
5 | 0.301 | 0.193 | 1.674 | 1.073 | 100.000 | 100.000 | ||||||
6 | 6.613×-16 | 3.626×-15 | 3.674×-15 | 2.015×-14 | 100.000 | 100.000 | ||||||
7 | 4.882×-16 | 4.941×-16 | 2.712×-15 | 2.745×-15 | 100.000 | 100.000 | ||||||
8 | 4.422×-16 | 3.984×-16 | 2.457×-15 | 2.214×-15 | 100.000 | 100.000 | ||||||
9 | 3.311×-16 | 3.045×-16 | 1.839×-15 | 1.692×-15 | 100.000 | 100.000 | ||||||
10 | 2.684×-16 | 2.875×-16 | 1.491×-15 | 1.597×-15 | 100.000 | 100.000 | ||||||
11 | 2.512×-16 | 8.763×-17 | 1.396×-15 | 4.868×-16 | 100.000 | 100.000 | ||||||
12 | 9.834×-17 | 2.067×-17 | 5.464×-16 | 1.148×-16 | 100.000 | 100.000 | ||||||
13 | 2.199×-17 | -6.063×-17 | 1.222×-16 | -3.368×-16 | 100.000 | 100.000 | ||||||
14 | -5.165×-17 | -1.136×-16 | -2.869×-16 | -6.312×-16 | 100.000 | 100.000 | ||||||
15 | -1.622×-16 | -2.030×-16 | -9.010×-16 | -1.128×-15 | 100.000 | 100.000 | ||||||
16 | -2.636×-16 | -2.639×-16 | -1.464×-15 | -1.466×-15 | 100.000 | 100.000 | ||||||
17 | -3.455×-16 | -4.063×-16 | -1.919×-15 | -2.257×-15 | 100.000 | 100.000 | ||||||
18 | -4.260×-16 | -5.056×-16 | -2.366×-15 | -2.809×-15 | 100.000 | 100.000 |
从基于相关矩阵的主成分分析结果(表 4和表 5)可以看出, 高AMF处理条件下, 云南蓝果树幼苗第2个主成分的特征值为2.725, 累计贡献率为93.327%, 说明干旱胁迫下高AMF对云南蓝果树的前2个主因子基本上能概括18个变量的主要信息。在第一、二主成分中, 将2个主成分的特征向量系数按照绝对值大小依次排列, 第一主成分居得分较高的的为地下部分生物量、株高和总生物量等7个指标, 得分均为0.90以上, 综合反映了云南蓝果树幼苗的生长情况;第二主成分得分为0.90以上的为WUE, 分值为0.931, 主要反映了叶片的水分利用效率。
指标Characteristic | 主成分Principal component | |||
1 | 2 | |||
AMF0 | AMF1 | AMF0 | AMF1 | |
地下部分生物量Belowground biomass | 0.993 | 0.995 | 0.082 | -0.039 |
株高Height | 0.983 | 0.989 | 0.125 | 0.082 |
总生物量Total biomass | 0.979 | 0.986 | 0.163 | 0.143 |
地上部分生物量Aboveground biomass | 0.969 | 0.970 | 0.207 | 0.220 |
相对生长速率Relative growth ration | 0.958 | 0.945 | 0.251 | 0.286 |
叶片数目Number of leaves | 0.923 | 0.981 | 0.383 | 0.187 |
叶面积Leave area | 0.905 | 0.821 | 0.398 | 0.554 |
Gs | 0.871 | 0.880 | 0.244 | 0.471 |
Pn | 0.796 | 0.950 | 0.508 | 0.312 |
Tr | 0.758 | 0.856 | 0.534 | 0.487 |
根冠比Aboveground biomass/Belowground biomass ratio | 0.536 | 0.238 | 0.265 | -0.824 |
地上部分生物量比Aboveground biomass/ Total biomass | 0.162 | 0.530 | 0.908 | 0.778 |
地下部分生物量比Belowground biomass/ Total biomass | -0.162 | -0.491 | -0.908 | -0.679 |
Ci | -0.450 | -0.633 | -0.886 | -0.733 |
Ls | 0.450 | 0.633 | 0.886 | 0.733 |
叶面积比率LAR | 0.415 | 0.692 | 0.820 | 0.648 |
水分利用效率WUE | -0.059 | 0.049 | 0.809 | 0.931 |
比叶面积SLA | 0.474 | 0.833 | 0.753 | 0.507 |
生长指标方面, 不管是高AMF还是低AMF条件下, 总生物量、地上部分生物量、地下部分生物量、叶片数目, Pn、Tr和Gs具有较高的可塑性指数(表 1, 表 2);高AMF条件下, 云南蓝果树幼苗各个参数及除Ls之外的其余5个光和特征指标都表现出较高的可塑性指数(表 1, 表 2)。通过对上述18个指标的综合分析, 不同处理下云南蓝果树幼苗抗旱性的隶属函数结果列于表 6。从表 6可以看出, 土壤水分含量在W1—W6范围时, 高AMF处理都增加了云南蓝果树幼苗的隶属函数值, 隶属函数值分别为0.55、0.68、0.54、0.45、0.43和0.44, 增强了其在各水分条件下的抗旱性。
指标Characteristic | W1 | W2 | W3 | W4 | W5 | W6 | ||||||
AMF0 | AMF1 | AMF0 | AMF1 | AMF0 | AMF1 | AMF0 | AMF1 | AMF0 | AMF1 | AMF0 | AMF1 | |
叶片数目Number of leaves | 0.47 | 0.44 | 0.42 | 0.46 | 0.45 | 0.72 | 0.53 | 0.40 | 0.35 | 0.29 | 0.30 | 0.44 |
株高Height | 0.53 | 0.90 | 0.49 | 0.57 | 0.57 | 0.33 | 0.58 | 0.44 | 0.57 | 0.30 | 0.52 | 0.86 |
相对生长速率Relative growth ration | 0.40 | 0.68 | 0.48 | 0.35 | 0.51 | 0.40 | 0.46 | 0.40 | 0.44 | 0.33 | 0.36 | 0.08 |
Pn | 0.60 | 0.56 | 0.53 | 0.43 | 0.59 | 0.29 | ||||||
WUE | 0.53 | 2.35 | 0.33 | 0.44 | 0.35 | 0.52 | ||||||
总生物量Total biomass | 0.43 | 0.41 | 0.63 | 0.28 | 0.51 | 0.72 | 0.40 | 0.49 | 0.44 | 0.58 | 0.37 | 0.48 |
地上部分生物量Aboveground biomass | 0.48 | 0.45 | 0.48 | 0.37 | 0.59 | 0.72 | 0.41 | 0.50 | 0.39 | 0.47 | 0.58 | 0.48 |
地下部分生物量Belowground biomass | 0.36 | 0.35 | 0.51 | 0.47 | 0.42 | 0.54 | 0.29 | 0.46 | 0.28 | 0.52 | 0.25 | 0.37 |
叶面积Leave area | 0.92 | 0.42 | 0.44 | 0.56 | 0.38 | 0.23 | ||||||
地上部分生物量比Aboveground biomass/Total biomass | 0.57 | 0.53 | 0.55 | 0.32 | 0.34 | 0.30 | ||||||
地下部分生物量比Belowground biomass/Total biomass | 0.43 | 0.47 | 0.45 | 0.04 | 0.47 | 0.47 | ||||||
平均隶属函数The average membership function | 0.51 | 0.55 | 0.49 | 0.68 | 0.50 | 0.54 | 0.40 | 0.45 | 0.41 | 0.43 | 0.38 | 0.44 |
本试验中采用苯菌灵灭菌的方法来抑制AMF对云南蓝果树实生苗根系的侵染, 得到理想的结果, 即杀真菌剂处理显著降低了云南蓝果树幼苗根侵染率, 使得本试验中AMF处理的实生苗在生长和光合特征上的差异可以用苯菌灵处理导致的侵染差异加以解释[16]。
干旱胁迫会限制植物的生长[17]。本研究表明, 干旱胁迫显著抑制了云南蓝果树幼苗的10个生长指标。土壤相对含水量较高时, 幼苗具有较多的叶片和较大的叶面积, 从而拥有较大的光合速率和较大的生物量。随着干旱胁迫程度的增加, 云南蓝果树幼苗的大部分生长指标都发生了显著的变化, 如叶片数目、叶面积和株高等8个指标都显著降低, 这与有些研究发现的轻度胁迫条件下某些植物可以通过形态的可塑性作出适应相反, 证明云南蓝果树缺乏应对干旱胁迫的策略;生物量分配方式也发生了显著变化, 水分充足时, 云南蓝果树将生物量较多地分配到地上部分的生长上, 可以提高光合产量, 满足植物的需要;而在干旱胁迫时, 大多数植物会降低地上部分生物量的分配, 将更多的资源分配到地下, 从而获得更多的水分和营养, 获得更多的地下部分生物量, 根冠比呈现先升高后降低的趋势[17-18], 然而, 云南蓝果树却随着干旱胁迫的增加, 根冠比没有发生显著的变化, 这一定程度上揭示了云南蓝果树根部表型可塑性差的特点;生物量分配方式也发生了紊乱, 例如地上部分生物量比和地下部分生物量比都呈现先升高后降低的趋势;LAR 和SLA 是调控植物功能的重要性状[19], 水分充足时, 云南蓝果树具有较大的LAR 和SLA, 从而保证了较高的光资源捕获面积, 两者共同保证了较高的生物量积累[17]。干旱胁迫导致SLA的减少是适应干旱的一种表现[20-21], 因此云南蓝果树的SLA表现出较强的可塑性。
干旱影响植物的生长在很大程度上要归结于光合作用的响应, 而光合作用对缺水特别敏感[17, 22-24]。付士磊等在研究杨树的光合作用与抗早能力的关系中指出, 在一定时间内, 随胁迫程度的增加, 叶片的净光合速率显著降低[25]。气孔可以快速感应大气湿度或根系水势的变化, 通过降低Gs可以有效防止水分散失。受旱植物为了尽量维持其Pn, 叶片Gs会降低。但是, 随着水分胁迫程度的增强, Ci和CO2浓度达到平衡, 出现先降低再升高的趋势, Ls和WUE则为先升后降[26]。本研究中, 云南蓝果树的Pn、Tr和Gs 3个指标随着干旱程度的增加呈显著递减状态, Ls和WUE则为先升后降的趋势。由此表明, 干旱胁迫条件下云南蓝果树幼苗不仅能获得的资源量较少, 而且光合机构的伤害造成了资源利用率的低下, 共同影响了云南蓝果树幼苗的生长。
3.2 不同水分条件下丛枝菌根真菌对云南蓝果树幼苗生长和光合特征的影响干旱胁迫条件下AM真菌对植物生长的影响已有一些研究[22, 27-30]。陈冬青等通过施加苯菌灵对云南蓝果树进行AMF处理, 发现水分胁迫下AMF对黄顶菊的生物量和生理指标影响显著, AMF共生能够促进黄顶菊土壤水分和矿质营养吸收, 改变植物代谢活动, 提高抗旱性[16]。贺学礼和赵丽莉发现接种摩西球囊霉(Glomus moseae)显著促进了小麦(Triticum aestivum)的营养生长和磷吸收, 提高了叶绿素含量, 增强了光合作用[29]。贺学礼和李生秀采用土培试验研究了水分胁迫条件下接种摩西球囊霉对玉米(Zea mays)营养生长和抗旱性的影响[31]。贺学礼等发现, 不同水分条件下, 接种AMF提高了民勤绢蒿的生长和抗旱性[32]。本试验中丛枝菌根的处理对干旱胁迫条件下云南蓝果树幼苗的生长产生了较大的影响, AMF侵染率与8个生长指标呈显著正相关关系。在胁迫前期, AMF对云南蓝果树幼苗的作用明显, 但是到了胁迫后期, 干旱胁迫对其影响大于AMF的影响, 导致AMF对重度干旱胁迫条件下的云南蓝果树幼苗没有作用。
AMF可以提高光合作用效率, 增强宿主植物的抗旱性[22, 27-28]。很多研究发现AMF可以提高很多植株叶片的光合速率(Pn)、蒸腾速率(Tr)、增加了气孔导度(Gs)和水分利用效率(WUE), 水分胁迫下促进作用更加显著[33]。Morte等在向日葵上也得出接种菌根能提高蒸腾速率、气孔导度和净光合速率等结论[34]。本试验中丛枝菌根的处理对云南蓝果树幼苗的光合特征产生了较大的影响, AMF处理能够提高云南蓝果树幼苗的光合速率、蒸腾速率、气孔导度及水分利用效率等参数, 增强其抗旱能力, 这无论对菌根化育苗, 还是云南蓝果树生态恢复都将具有极其重要的实践意义, 而且云南蓝果树幼苗生长指标和光合参数之间的相关性比较强, 意味着AMF处理导致的净光合速率等光合能力的提高, 使得云南蓝果树幼苗生长指标的增加。
3.3 不同水分条件下丛枝菌根真菌对云南蓝果树幼苗抗旱性的影响可塑性指数在一定程度上可以反映植物的可塑性程度[17]。不论是水分充足还是干旱胁迫条件下, 高AMF处理的云南蓝果树幼苗生长和光合参数的可塑性指数普遍较高, 可以较好地适应干旱胁迫等逆境环境。
隶属函数分析提供了在多指标测定的基础上对材料进行综合评价的方法, 更具科学性和可靠性[35]。本研究通过主成分分析简化了评价指标体系, 确定了生长和光合特征的主要结构指标, 并采用隶属函数值法进行了综合评价。结果显示, 不管土壤水分条件如何(从W1到W6) , AMF处理都显著影响了幼苗的抗旱性, 即高AMF处理下幼苗的平均隶属函数值更高, 具有更强的抗旱性。
研究证实水分显著影响丛枝菌根侵染率。Kaya等发现菌根侵染会显著提高黄瓜的产量及水分利用效率[36]。但是, 从云南有气象记录的1961年以来, 云南的年降水量则出现减少的趋势, 半个世纪以来年降水量减少了39mm, 减少速率为-8mm/10a, 其中夏季和秋季减少趋势明显于春季和冬季。最显著的例子是云南蓝果树原产地西双版纳年降水日由20世纪50年代的每年270d锐减到目前的150d;年雾日由180d减少到30d, 以往湿润的热带雨林气候已经发生明显变化。这有可能是导致云南蓝果树菌根侵染率降低的主要原因, 而AMF侵染率的降低, 影响了云南蓝果树幼苗的抗旱性。因此, 云南蓝果树抗旱性不强的特点, 加上逐渐干旱的原生境, 种群逐渐变小, 使得此物种的处境岌岌可危, 如果不加以适当保护, 本物种必将灭绝。
综合分析表明, 干旱胁迫对云南蓝果树幼苗的生长和光合特征均产生了显著的影响, 而云南蓝果树幼苗本身的抗旱性较差, 但是AMF能与云南蓝果树根系形成良好的共生关系, 一定程度上可以缓解干旱胁迫对植株所造成的伤害, 增强其抗旱性, 促进宿主植物生长。这说明自然生境下云南蓝果树幼苗数量少, 天然更新缺乏很大一部分原因可以归结于当地生境的干旱, 加上可以缓解其干旱胁迫的AMF也受当地气候导致水分缺乏的影响, 也间接影响了云南蓝果树的抗旱性。目前对云南蓝果树开展了一系列保护措施, 例如就地保护、近地保护、迁地保护和回归引种等。保护地的云南蓝果树近几年开始出现生长缓慢、部分死亡的现象。基于本研究云南蓝果树抗旱性不强的结果, 在对云南蓝果树开展保护措施的过程中, 建议将云南蓝果树种植在水源充足的生境, 以便为云南蓝果树生长提供足够的水分;基于AMF可以增强云南蓝果树抗旱性的研究结果, 建议对保护地的云南蓝果树施加AMF菌剂。
[1] | 陈伟, 史富强, 杨文忠, 周云, 陈宏伟. 云南蓝果树的种群状况及生态习性. 东北林业大学学报 , 2011, 39 (9) : 17–19. |
[2] | 孙宝玲, 张长芹, 周凤林, 史富强, 吴之坤. 极度濒危植物——云南蓝果树的种子形态和不同处理条件对种子萌发的影响. 云南植物研究 , 2007, 29 (3) : 351–354. |
[3] | 孙宝玲, 张长芹. 极度濒危植物云南蓝果树的形态修订. 云南植物研究 , 2007, 29 (2) : 173–175. |
[4] | Poorter L, Hayashida-Oliver Y. Effects of seasonal drought on gap and understorey seedlings in a Bolivian moist forest. Journal of Tropical Ecology , 2000, 16 (4) : 481–498. DOI:10.1017/S026646740000153X |
[5] | Engelbrecht B M J, Kursar T A. Comparative drought-resistance of seedlings of 28 species of co-occurring tropical woody plants. Oecologia , 2003, 136 (3) : 383–393. DOI:10.1007/s00442-003-1290-8 |
[6] | McDowell N, Pockman W T, Allen C D, Breshears D D, Cobb N, Kolb T, Plaut J, Sperry J, West A, Williams D G, Yepez E A. Mechanisms of plant survival and mortality during drought:why do some plants survive while others succumb to drought? New Phytologist, 2008, 178(4):719-739. |
[7] | 刘文杰, 李红梅. 西双版纳旅游气候资源. 自然资源 , 1997 (2) : 62–66. |
[8] | 宋富强, 赵俊斌, 张一平, 许再富, 肖来云. 西双版纳区域气候变化对植物生长趋势的影响. 云南植物研究 , 2010, 32 (6) : 547–553. |
[9] | Smith S E, Read D J. Mycorrhizal Symbiosis. 2nd ed. San Diego:Academic Press, 1997:233-289. |
[10] | Varma A. Mycorrhizae-the friendly fungi:What we know, what should we know, and how do we know?//Varma A ed. Mycorrhiza Manual. Berlin:Springer, 1998:1-24. |
[11] | Rai M, Acharya D, Singh A, Varma A. Positive growth responses of the medicinal plants Spilanthes calva and Withania somnifera to inoculation by Piriformospora indica in a field trial. Mycorrhiza , 2001, 11 (3) : 123–128. DOI:10.1007/s005720100115 |
[12] | Wang X R, Pan Q, Chen F X, Yan X L, Liao H. Effects of co-inoculation with arbuscular mycorrhizal fungi and rhizobia on soybean growth as related to root architecture and availability of N and P. Mycorrhiza , 2011, 21 (3) : 173–181. DOI:10.1007/s00572-010-0319-1 |
[13] | Veresoglou S D, Shaw L J, Sen R. Glomus intraradices and Gigaspora margarita arbuscular mycorrhizal associations differentially affect nitrogen and potassium nutrition of Plantago lanceolata in a low fertility dune soil. Plant and Soil , 2011, 340 (1/2) : 481–490. |
[14] | Klironomos J N, McCune J, Hart M, Neville J. The influence of arbuscular mycorrhizae on the relationship between plant diversity and productivity. Ecology Letters , 2000, 3 (2) : 137–141. DOI:10.1046/j.1461-0248.2000.00131.x |
[15] | 何跃军, 钟章成, 刘锦春, 刘济明. 石灰岩土壤基质上构树幼苗接种丛枝菌根(AM)真菌的光合特征. 植物研究 , 2008, 28 (4) : 452–457. |
[16] | 陈冬青, 皇甫超河, 刘红梅, 王楠楠, 杨殿林. 水分胁迫和杀真菌剂对黄顶菊生长和抗旱性的影响. 生态学报 , 2013, 33 (7) : 2113–2120. |
[17] | 徐飞, 郭卫华, 徐伟红, 王仁卿. 刺槐幼苗形态、生物量分配和光合特性对水分胁迫的响应. 北京林业大学学报 , 2010, 32 (1) : 24–30. |
[18] | Ackerly D D, Bazzaz F A. Leaf dynamics, self-shading and carbon gain in seedlings of a tropical pioneer tree. Oecologia , 1995, 101 (3) : 289–298. DOI:10.1007/BF00328814 |
[19] | 王满莲, 冯玉龙. 紫茎泽兰和飞机草的形态、生物量分配和光合特性对氮营养的响应. 植物生态学报 , 2005, 29 (5) : 697–705. |
[20] | Abrams M D. Adaptations and responses to drought in Quercus species of North America. Tree Physiology , 1990, 7 (1/4) : 227–238. |
[21] | 郑淑霞, 上官周平. 黄土高原油松和刺槐叶片光合生理适应性比较. 应用生态学报 , 2007, 18 (1) : 16–22. |
[22] | 宋会兴, 彭远英, 钟章成. 干旱生境中接种丛枝菌根真菌对三叶鬼针草(Bidens pilosa L.)光合特征的影响. 生态学报 , 2008, 28 (8) : 3744–3751. |
[23] | 孙存华, 李扬, 杜伟, 金会丽, 王东升, 陈湘玲, 徐新娜. 干旱胁迫下藜的光合特性研究. 植物研究 , 2007, 27 (6) : 715–720. |
[24] | 陈昕, 徐宜凤, 张振英. 干旱胁迫下石灰花楸幼苗叶片的解剖结构和光合生理响应. 西北植物学报 , 2012, 32 (1) : 111–116. |
[25] | 付士磊, 周永斌, 何兴元, 陈炜. 干旱胁迫对杨树光合生理指标的影响. 应用生态学报 , 2006, 17 (11) : 2016–2019. |
[26] | Lawlor D W, Tezara W. Causes of decreased photosynthetic rate and metabolic capacity in water-deficient leaf cells:a critical evaluation of mechanisms and integration of processes. Annals of Botany , 2009, 103 (4) : 561–579. |
[27] | 田帅, 刘振坤, 唐明. 不同水分条件下丛枝菌根真菌对刺槐生长和光合特性的影响. 西北林学院学报 , 2013, 28 (4) : 111–115. |
[28] | 梁宇, 高玉葆, 陈世苹, 任安芝. 干旱胁迫下内生真菌感染对黑麦草实验种群光合、蒸腾和水分利用的影响. 植物生态学报 , 2001, 25 (5) : 537–543. |
[29] | 贺学礼, 赵丽莉. 非灭菌条件下VA菌根真菌对小麦生长发育的影响. 土壤通报 , 1999 (2) : 57–59. |
[30] | 赵金莉, 贺学礼. AM真菌对油蒿生长和抗旱性的影响. 华北农学报 , 2007, 22 (5) : 184–188. |
[31] | 贺学礼, 李生秀. 不同VA菌根真菌对玉米生长及抗旱性的影响. 西北农业大学学报 , 1999, 27 (6) : 49–53. |
[32] | 贺学礼, 高露, 赵丽莉. 水分胁迫下丛枝菌根AM真菌对民勤绢蒿生长与抗旱性的影响. 生态学报 , 2011, 31 (4) : 1029–1037. |
[33] | Sánchez-Blanco M J, Ferrández T, Morales M A, Morte A, Alarcón J J. Variations in water status, gas exchange, and growth in Rosmarinus officinalis plants infected with Glomus deserticola under drought conditions. Journal of Plant Physiology , 2004, 161 (6) : 675–682. DOI:10.1078/0176-1617-01191 |
[34] | Morte A, Lovisolo C, Schubert A. Effect of drought stress on growth and water relations of the mycorrhizal association Helianthemum almeriense-Terfezia claveryi. Mycorrhiza , 2000, 10 (3) : 115–119. DOI:10.1007/s005720000066 |
[35] | 韩瑞宏, 卢欣石, 高桂娟, 杨秀娟. 紫花苜蓿抗旱性主成分及隶属函数分析. 草地学报 , 2006, 14 (2) : 142–146. |
[36] | Kaya C, Higgs D, KiInak H, Tas I. Mycorrhizal colonisation improves fruit yield and water use efficiency in watermelon (Citrullus lanatus Thunb.) grown under well-watered and water-stressed conditions. Plant and Soil, 2003, 253(2):287-292. |