凹葉厚樸樹體厚樸酚、和厚樸酚含量變化模型的研究

王曉明 ，楊碩知 ，裴 剛 ，蔣麗娟 ，宋慶安 ，蔡 能 ，李永欣 ，曾慧杰

(1.湖南省林業科學院，長沙 410004；2.湖南省林木無性系育種重點實驗室，長沙 410004；3.中南林業科技大學，長沙 410004；4.湖南中醫藥大學，長沙 410208)

凹葉厚樸樹體厚樸酚、和厚樸酚含量變化模型的研究

王曉明1，2，3，楊碩知1，2，裴 剛4，蔣麗娟3，宋慶安1，2，蔡 能1，2，李永欣1，2，曾慧杰1，2

(1.湖南省林業科學院，長沙 410004；2.湖南省林木無性系育種重點實驗室，長沙 410004；3.中南林業科技大學，長沙 410004；4.湖南中醫藥大學，長沙 410208)

為了研究凹葉厚樸樹厚樸酚、和厚樸酚的積累分布規律，以10年生凹葉厚樸為試驗材料，建立了凹葉厚樸樹體解析木高度與厚樸酚、和厚樸酚及樹皮厚度含量的數學模型。結果表明：解析木高度與總酚含量、厚樸酚含量、和厚樸酚含量關系最優模型依次為y=─0.053 636x3+0.665 503x2─2.543 869x+3.964 848，y=─0.008 939x3+0.110 233x2─0.419 464x+0.693 864，y=─0.044 697x3+0.555 271x2─2.124 405x+3.270 985；樹體解析木高度與樹皮厚度關系最優模型為y=─0.008 965x3+0.031 331x2─0.143 492x+4.941 212，可用于凹葉厚樸樹體厚樸酚、和厚樸酚含量的預測，為科學采收厚樸皮提供了理論依據。

凹葉厚樸；數學模型；厚樸酚；和厚樸酚；樹皮厚度

凹葉厚樸(Magnolia officinalis subsp. biloba)是我國重要的藥材兩用經濟林樹種之一，主要分布區是浙江南部、福建西部、江西西南部、湖南南部、廣東東北部、廣西北部及貴州的東部[1]，其藥用主要活性成分是樹皮中所含的厚樸酚、和厚樸酚等酚類物質[2]，其藥材品質優劣主要反映在這2種成份的差異上[3]。

種源、樹齡、產地和立地條件等與厚樸產量及厚樸酚、和厚樸酚含量關系密切[4-6]，是影響其栽培經濟效益的重要因素。有關凹葉厚樸的厚樸酚、和厚樸酚含量測定及影響樹皮厚度的因子和不同藥用部位中酚類成分的含量分析等方面已有研究報道[7-10]，涂育合等[11]也建立了凹葉厚樸樹皮產量預測模型，但凹葉厚樸樹體不同部位中厚樸酚、和厚樸酚積累分布規律至今尚未見研究報道。為了更好地評價凹葉厚樸人工林的生產潛力，探討其主要活性成分含量的積累分布規律,筆者嘗試建立凹葉厚樸樹體解析木高度與厚樸酚、和厚樸酚含量及樹皮厚度之間的相關數學模型，旨在為凹葉厚樸樹體厚樸酚、和厚樸酚含量預測提供理論依據，從而科學指導凹葉厚樸的生產。

1 試驗地概況

試驗地設在湖南省永州市道縣洪塘營瑤族鄉海拔1 200 m處。所處地理位置屬溫、熱帶結合部，既具陽光豐富的大陸性季風氣候特點，又有雨量充沛、空氣濕潤的海洋性氣候特征。全年平均氣溫在17.6℃～18.6℃之間，一月平均氣溫7.4℃，七月平均氣溫26.5℃，極端高溫36℃。極端最低氣溫在-4.9 ～ -8.4℃之間。年平均日照時數為1300～1740 h，無霜期年均285～311 d，年降雨量1290～1900 mm。試驗地士壤為花崗巖發育的山地黃壤為主，土層深厚，土壤質地好，微酸性至中性。

2 材料與方法

1.1 材 料

以湖南省永州市道縣10年生凹葉厚樸為調查對象及試驗材料，測定樹體不同高度的樹皮厚度及厚樸酚、和厚樸酚和厚樸總酚含量。

1.2 方 法

1.2.1 采樣方法

在同一坡面的上、中、下坡選擇大小一致的5株樹為試驗對象，分別伐倒，先測定其樹高、胸徑，然后以1 m為一段，按陰、陽面分段采割樹皮，測定樹皮厚度。將新鮮厚樸皮運回室內，自然通風陰干備用。

1.2.2 測定方法

(1) 儀器 Agilent 1200 HPLC四元泵，Agilent 1200工作站。

(2) 色譜條件 色譜柱：迪馬DIKMA C18柱(250 mm ×4.6 mm，5 μm)；流動相：甲醇∶水(73∶27)；波長：294 nm；柱溫：30℃；流速：1.0 mL/min；理論塔板數以厚樸酚峰計不低于4 000。

(3) 對照品溶液的制備及標準曲線的測定：用甲醇分別制成濃度為厚樸酚標準品溶液0.08 mg/mL、和厚樸酚標準品溶液 0.096 mg/mL測厚樸酚、和厚樸酚各自的回歸方程。

(4) 供試品溶液的制備：分別取不同的厚樸樣品0.1g，粉碎(過三號篩)，精密稱定，置50 mL的磨口錐形瓶中，精密加入甲醇25 mL，稱重，超聲處理20 min，取出放冷至室溫，稱重，用溶劑補足減少的甲醇量，搖勻，濾過，取續濾液再用0.45 μm的微孔濾膜濾過，濾液即為供試品溶液。

(5) 精密度、穩定性和加樣回收率測定：按同樣色譜條件進行HPLC分析，每次進樣10 μL，連續進樣5次，計算厚樸酚、和厚樸酚的含量后，衡量儀器的穩定性。

分別在0、2、4、6、8、24 h 測定厚樸酚對照溶液和和厚樸酚對照品溶液混合液，計算厚樸酚、和厚樸酚的含量，衡量二者穩定性。

精密加入厚樸酚、和厚樸酚至已知含量的厚樸原藥材中，按樣品溶液制備方法操作，制備樣品液，按色譜條件測定厚樸酚、和厚樸酚的含量，計算加樣回收率。

(6) 樣品溶液的制備取樣品粉末（過三號篩）約0.2 g, 精密稱定，置具塞錐形瓶中，精密加入甲醇25 mL，搖勻，密塞，浸漬24 h，濾過，精密量取續濾液5 mL，置于25 mL 容量瓶中，加甲醇至刻度，搖勻，即得。

1.2.3 建模方法

對收集到的數據運用SPSS統計軟件分析，以解析木高度為自變量(x, m)，樹皮厚度或酚類物質的含量為因變量(y, cm, %)，采用線性方程、二次方程、三次方程、指數方程、對數方程、乘冪曲線、復合曲線、S曲線模型、Logostic生長曲線、反響曲線、等比級數曲線11種曲線進行擬合，以期找到適宜的該凹葉厚樸樹體解析木高度與樹皮厚度及厚樸總酚、厚樸酚與和厚樸酚含量的數學模型。

2 結果與分析

2.1 凹葉厚樸解析木不同高度與厚樸酚、和厚樸酚含量的模型擬合

曲線擬合結果表明: 線性方程模型式(1)、二次方程模型式(2)、三次方程模型式(3)及復合曲線模型式(4)均達到極顯著水平(r0.01=0.87400,r0.05=0.70700,下同)，指數方程、對數方程等其它模型沒有達到顯著水平。根據模型擬合結果(見表1、2、3)，在凹葉厚樸解析木不同高度分別與同一高度陰面、陽面的平均總酚含量、厚樸酚含量、和厚樸酚含量的回歸關系中，其相關系數按模型式(3)、(2)、(4)、(1)的順序依次減小，其中模型式(3)的相關系數最大，總酚含量-解析木高度模型、厚樸酚含量-解析木高度模型、和厚樸酚含量-解析木高度模型的相關系數分別為0.950 50、0.938 97、0.952 61。從精度來說, 三次方程模型式(3)為最優的擬合模型，可以用作凹葉厚樸樹體不同高度的厚樸酚、和厚樸酚、總酚含量分布規律的預測依據。

表1 凹葉厚樸解析木總酚含量相關性預測模型Table 1 The predictive model of correlation between height and total phenol content for Magnolia officinalis subsp. biloba

表2 凹葉厚樸解析木厚樸酚含量相關性預測模型Table 2 The predictive model of correlation between height and magnolol content for Magnolia officinalis subsp. biloba

表3 凹葉厚樸解析木和厚樸酚含量相關性預測模型Table 3 The predictive model of correlation between height and honokiol content for Magnolia officinalis subsp. biloba

樹體陰面、陽面的光照時間和強度不同，其樹皮活性成分累積也有差異。凹葉厚樸解析木陰面、陽面數據分別擬合的結果(見表4、5、6、7、8、9)與陰陽面平均數據擬合的結果一致(見表1、2、3)，只有線性方程模型、二次方程和三次方程模型、復合曲線模型達到極顯著水平，其它模型均沒有達到顯著水平。模型式(3)、(2)、(4)、(1)的相關系數依次減小，以三次方程模型式(3)的相關系數最大，為各模型中的最佳擬合模型。從凹葉厚樸解析木樹體高度與陰面、陽面的總酚、厚樸酚、和厚樸酚擬合模型式(3)中的殘差平方和可看出，陽面數據擬合度都優于陰面數據的擬合效果，其中陽面的厚樸酚—樹木高度、和厚樸酚—樹木高度的模型殘差平方和分別比陽面的減少13.49%和12.45%。這可能與凹葉厚樸樹體陰面比陽面的光照時間較短，光照強度較弱，酚類物質的累積較小有關，從而導致曲線擬合的效果也較差。

表4 凹葉厚樸解析木陽面總酚含量相關性預測模型Table 4 The predictive model of correlation between height and total phenols content for Magnolia officinalis subsp. biloba (sunny side)

表5 凹葉厚樸解析木陽面厚樸酚含量相關性預測模型Table 5 The predictive model of correlation between height and magnolol content for Magnolia officinalis subsp. biloba (sunny side)

表6 凹葉厚樸解析木陽面和厚樸酚含量相關性預測模型Table 6 The predictive model of correlation between height and honokiol content for Magnolia officinalis subsp. biloba (sunny side)

2.2 凹葉厚樸解析木不同高度與樹皮厚度的模型擬合

利用解析木不同高度與樹皮陰陽面平均厚度的數據進行曲線擬合(表10)，結果表明，線性方程模型、二次方程和三次方程模型、復合曲線模型達到極顯著水平，其它模型均沒有達到顯著水平。但以三次方程模型(3)的相關系數（0.99644）最大，其殘差平方和（0.0391）最小，是4種擬合模型中的最優模型，適合作為凹葉厚樸樹皮陰陽面平均厚度的預測模型。

表7 凹葉厚樸解析木陰面總酚含量相關性預測模型Table 7 The predictive model of correlation between height and total phenols content for Magnolia officinalis subsp. biloba (night side)

表8 凹葉厚樸解析木陰面厚樸酚含量相關性預測模型Table 8 The predictive model of correlation between height and magnolol content for Magnolia officinalis subsp. biloba (night side)

表9 凹葉厚樸解析木陰面和厚樸酚含量相關性預測模型Table 9 The predictive model of correlation between height and honokiol content for Magnolia officinalis subsp. biloba (night side)

分別將凹葉厚樸解析木陰面和陽面的樹皮厚度與樹體不同高度的數據進行曲線擬合，其擬合模型的相關系數大小排序與陰陽面平均數據的擬合結果相似，全部達到極顯著水平(見表11、12)，且以三次方程模型式(3)相關系數最大。以模型式(3)分析比較陰面、陽面數據模型的相關系數和殘差平方和的差異，陰面的相關系數大于陽面的，其殘差平方和比陽面減小164.34%，說明陰面數據模型的擬合效果比陽面的好些。這可能與樹體陰面一天內光照較為穩定，樹皮厚度較均勻，而陽面從早到晚接受光照強度以及角度有較大的變化有關，故而，陽面樹皮厚度均勻性較差，從而導致其數據模型擬合效果不如陰面。

表10 凹葉厚樸解析木樹皮厚度預測模型Table 10 The predictive model of bark thickness for Magnolia officinalis subsp. biloba

表11 凹葉厚樸解析木陽面不同高度樹皮厚度預測模型Table 11 The predictive model of bark thickness for Magnolia officinalis subsp. biloba (sunny side)

表12 凹葉厚樸解析木陰面樹皮厚度預測模型Table 12 The predictive model of bark thickness for Magnolia officinalis subsp. biloba (night side)

3 結論與討論

本文中，采用曲線擬合法對10年生凹葉厚樸解析木高度與樹皮厚度及厚樸酚含量、和厚樸酚等含量進行了曲線擬合，并建立了相應的數學模型。從結果來看，在已采用的11種曲線模型中，擬合效果均以三次方程模型為最佳。其中，解析木高度與總酚含量、厚樸酚含量、和厚樸酚含量的最優模型分別為：

解析木高度與樹皮厚度相關性最優模型為：

y =－0.008 965x3+0.031 331x2-0.143 492x+4.941 212。

有報道凹葉厚樸樹體不同部位的厚樸酚與和厚樸酚含量有顯著差異，根皮的酚含量遠高于干皮和枝皮;干皮自下而上酚含量逐漸降低;枝皮的酚含量比中上部干皮的高[7，10]。本的研究結果表明，凹葉厚樸樹皮厚度雖然隨著樹體高度增加而逐步遞減趨勢，但并非線性均勻減少，總酚、厚樸酚與和厚樸酚含量隨著樹體高度增加而呈現減少的趨勢，但其中反復增減的變化。從試驗數據來看，凹葉厚樸樹體基部1 m以內主要活性成分含量最高，1 m以上解析木各段最高含量僅為基部段的40%。有研究認為，陽坡對厚樸皮厚生長有利，陰坡厚樸的皮往往比較薄[5]。本研究也發現凹葉厚樸樹體各個主，度段的陽面樹皮厚度均比陰面略厚，兩者的研究結論相似。

本文中，對凹葉厚樸樹皮厚度與厚樸酚、和厚樸酚含量關系也進行了曲線擬合，但結果均未達到顯著水平。多數研究者認為凹葉厚樸樹皮越厚品質越好，但有報道稱樹皮厚度隨年齡不斷增加，但其此生代謝強度到一定年齡會開始降低，因此樹皮厚度與厚樸酚、和厚樸酚含量相關性不密切[6]，本研究結論與該觀點一致。

綜上所述，本研究建立了凹葉厚樸樹體不同高度與厚樸酚、和厚樸酚、厚樸總酚含量以及樹皮厚度的數學模型，揭示了凹葉厚樸樹體不同部位的主要活性成分及樹皮厚度的變化規律，為科學預測凹葉厚樸產量及主要活性成分含量提供了的理論依據，從而有利于科學地指導凹葉厚樸生產，提高凹葉厚樸種植的經濟效益，增加山區農民收入，促進農村經濟發展。

若要進一步提高凹葉厚樸樹體主要活性成分含量分布預測的可靠性，應該調查研究各種樹齡的樹體，全面把握凹葉厚樸不同年齡階段的主要活性成分分布規律。本文由于材料限制，在此未及詳述。

[1] 謝韻帆，楊志玲，譚梓峰,等.永州市厚樸資源及產業化建設探討[J].湖南林業科技，2004, 3l (2): 49-50.

[2] 《中華人民共和國藥典》（2010版）[M].北京：中國醫藥科技出版社，2010.

[3] 童再康，斯金平，劉 饒.不同種源厚樸酚類物質含量變異與遺傳的初步研究[J].林業科學研究，2000, 13(3): 357-261.

[4] 張春霞 , 楊立新, 余 星,等.種源、產地及采收樹齡對厚樸藥材質量的影響[J].中國中藥雜志，2009，34（9）：2431-2437.

[5] 鄧白羅.厚樸生長與立地因子關系的調查研究[J].經濟林研究，1994, 12(1):41-46 .

[6] 曾燕如，童再康，朱玉球，等. 樹齡與厚樸質量關系的研究[J]。藥用植物栽培，1999, 22(8): 379-381

[7] 楊紅兵，詹亞華，陳科力,等.湖北恩施產厚樸中厚樸酚與和厚樸酚的定量分析[J].中國醫院藥學雜志，2007, 21(6):767-768.

[8] 劉壽強.凹葉厚樸樹皮量及影響因子的研究[J].經濟林研究，2001, 19(2) :29-30.

[9] 李安娟，郭信芳，馮海龍,等.厚樸的不同藥用部位中酚類成分的含量及其產地加工對含量的影響[J].中藥通報，1985,10(4) :10-13.

[10] 斯金平，潘心平，劉 饒,等.人工栽培凹葉厚樸中厚樸酚與和厚樸酚含量的研究[J].浙江林業科技，1994, 14(5) :21-22.

[11] 涂育合，葉功富，林照授，等.凹葉厚樸樹皮產量預測模型的研究[J].安徽農業大學學報，2003, 30(3): 312-315.

Study on content variation model of magnolol and honokiol in Magnolia officinalis subsp. biloba

WANG Xiao-ming1,2,3，YANG Shuo-zhi1,2，PEI Gang4，JIANG Li-juan3，SONG Qing-an1,2，CAI Neng1,2，LI Yong-хin1,2，ZENG Hui-jie1,2
(1. Hunan Academy of Forestry，Changsha 410004，China；2.Hunan Key Laboratory of Trees Clones Breeding Technology，Changsha 410004，China；3.Central South Forestry and Technology University，Changsha 410004，China；4.Hunan Traditional Chinese Medicine University，Changsha 410007，China)

Ten-year-old Magnolia officinalis subsp. biloba were used as the eхperimental material for the purpose of studying on the regularity of accumulation and distribution of magnolol and honokiol in trees, and mathematical models were established. The results show that the best model of relation between the tree’s height and total phenol, magnolol, honokiol are y= -0.053 636x3+0.665 503x2-2.543 869x+3.964 848, y= -0.008 939x3+0.110 233x2-0.419 464x+0.693 864, y= -0.044 697x3+0.555 271x2-2.124 405x+3.270 985, the best model of relation between the tree’s height and bark thickness is y= -0.008 965x3+0.031 331x2-0.143 492x+4.941 212. These models can be used in forecasting magnolol and honokiol content and providing scientific basis in bark harvesting of Magnolia officinalis subsp. biloba.

Magnolia officinalis subsp. biloba; mathematic model; magnolol; honokiol; bark thickness

S794.08; S759.8

A

1673-923X(2012)02-0001-05

2011-11-24

國家林業公益性行業科研專項 “金銀花、厚樸優良新品種創制及利用技術研究”（201104023）

王曉明(1964—)，男，湖南永興人，研究員，主要從事木本藥用植物及觀賞植物育種與栽培研究；E-mail：wхm1964@163.com

[本文編校：邱德勇]