盛玉珍時小東王　麗高繼海莊國慶**

（1.四川省林業(yè)科學研究院，成都　610081；2.四川大學生命科學學院，成都　610064；3.成都中醫(yī)藥大學藥學院，四川　成都　611137）

?

川產(chǎn)楨楠礦物元素因子分析及聚類分析*

盛玉珍1時小東2王麗1高繼海3莊國慶1**

（1.四川省林業(yè)科學研究院，成都610081；2.四川大學生命科學學院，成都610064；3.成都中醫(yī)藥大學藥學院，四川成都611137）

摘要：采用等離子體原子發(fā)射光譜法測定了不同苗齡楨楠幼苗和川產(chǎn)楨楠木材中9種礦物元素的含量，并運用SPSS軟件對數(shù)據(jù)進行了因子分析和聚類分析。結果表明，不同苗齡楨楠幼苗中Ca和S元素變化呈現(xiàn)出富集的變化趨勢，其他元素變化較?。徊煌a(chǎn)地楨楠元素間存在一定的關聯(lián)性，主成分分析表明，其87%以上的貢獻來源于前4個主成分，對應元素分別為Ca、Mn、Mo、K；聚類分析表明，川西北和西部等地楨楠聚集為一類，川西南和川南等地楨楠聚集為另一類。運用該方法初步分析了礦物元素和楨楠生長的關系，同時驗證了該方法在楨楠聚類分析中的可行性。

關鍵詞：川產(chǎn)楨楠；礦物元素；聚類分析

*國家自然基金青年科學基金項目（批準號：31200504）

楨楠（Phoebe zhennan）是我國特有的樟科楠[1]，其木材優(yōu)良、堅實，是我國二級重點保護瀕危樹種，屬亞熱帶常綠喬木，主要集中分布于我國長江流域及以南的地區(qū)，如云南、四川、湖北、貴州、廣西等省份[2 - 4 ]。四川是我國楨楠的歷史分布中心，但目前楨楠的分布星散，群體結構衰退趨勢明顯，使得資源日趨瀕危[5]。加之楨楠本身生長發(fā)育緩慢，嚴重制約了楨楠人工林的發(fā)展和對楨楠木材的需求。

關于楨楠元素的研究起步相對較晚，僅有賀維研究了施肥對楨楠幼苗生長和光合作用等方面的作用[6- 11]。目前，指紋圖譜技術已經(jīng)被廣泛應用于植物產(chǎn)地溯源、藥材質量控制等方面[12- 13]，基于礦物元素的指紋圖譜構建和輔助資源評價，已經(jīng)在多種物種中進行了研究，如咖啡、水稻、薯蕷、蕓豆、枸杞等[14 - 18 ]。其中，因子分析和聚類分析是處理化學計量學分析方法中獲得大量數(shù)據(jù)的常用技術和有力工具[19 ]。主成分分析作為一種常用的多元統(tǒng)計分析方法，可以分析各指標對主成分的貢獻大小，用最少的因子去解釋待觀測的事物，以揭示事物本身之間的關系。系統(tǒng)聚類是目前國內外使用最多的一種聚類方法，其中分類距離是系統(tǒng)聚類分析的依據(jù)和基礎。近年來還沒有關于楨楠礦物元素方面的研究報道，并且，基于化學物質或礦物元素的楨楠因子分析和聚類研究等方面也沒有開展。本文以四川楨楠為研究對象，利用ICP法測定了楨楠木材中9種礦物元素的含量，并對不同樹齡中元素差異和因子分析進行研究，以期為楨楠規(guī)范化栽培以及產(chǎn)業(yè)化發(fā)展提供理論依據(jù)。

1　材料與方法

1.1試驗材料

選取大邑野生楨楠樹木種子，萌發(fā)后種植于四川省林業(yè)科學研究院沙河堡實驗基地，第3年的12月（2014年12月）選取生長一致的楨楠苗6棵，自然條件下晾干，去除葉片和側枝等，分別選取1~3年不同生長年齡的主莖段，進行礦物元素測定和分析。

四川楨楠分布廣泛，對采集地進行實地考察、調研，2014年10月選取11個地區(qū)野生成材楨楠樹進行研究，截取當年生側枝（距離頂端10 cm莖段）進行礦物元素測定和聚類分析。樣品來源及具體信息見表1。

表1　楨楠樣品的來源

1.2方法

按照DB53/T 288- 2009標準方法進行樣品處理和元素測定，分別測定Mo、Zr、Sr、Fe、Mn、Cr、Ca、K、S等9種礦物元素，每個樣品測定重復3次。

1.3數(shù)據(jù)分析

Kaiser- Meyer- Olkin（KMO）檢驗和Bartletts’s球形檢驗、主成分分析（Principal component analysis，PCA）、Pearson相關性、余弦相關性和聚類分析（Hierarchical cluster analysis，HCA）等采用軟件SPSS完成，其他數(shù)據(jù)處理用EXCEL完成。

2　結果與分析

2.1礦物元素含量隨樹齡變化及分析

由圖1可知，楨楠礦物元素中Ca含量最高，其次是K和S，Sr含量最少。對于不同苗齡而言，除Fe、Cr和S元素外，其他各元素在不同苗齡之間雖然存在差異，但是相對差異較小，沒有出現(xiàn)較大的波動。這可能是由于在試驗環(huán)境中，各種營養(yǎng)元素和水分充足，完全能夠滿足楨楠的生長需要，不會出現(xiàn)因環(huán)境壓力造成的元素缺乏脅迫。Ca和S元素變化呈現(xiàn)出富集的變化趨勢，可能是在楨楠樹齡增加的過程中，需要吸收Ca和S元素進行生理活動，所以在楨楠生長過程中需要注意Ca和S元素的添加，這一研究結果與我們對楨楠不同樹齡的轉錄組差異基因分析，發(fā)現(xiàn)了與S元素代謝有關的結果一致（數(shù)據(jù)還未發(fā)表）。

2.2不同產(chǎn)地楨楠木材礦物元素測定及分析

2.2.1測定結果及適當性量數(shù)檢驗

對11種川產(chǎn)楨楠材料的9種元素進行測定，發(fā)現(xiàn)各礦物質成分指標平均含量相差較大，范圍從9.99±2.07~5309.62±1358.97ppm。Ca元素的平均含量最高；其次是K元素和S元素；含量最低的是Sr元素。同時，各樣本間元素含量變異系數(shù)較大，表明不同川產(chǎn)楨楠元素含量具有較大差異。其中各元素變異系數(shù)大小順序為：Ca>K>S>Mn>Mo> Sr>Zr>Cr>Fe，說明川產(chǎn)楨楠中Ca元素的變異程度最大，含量變化范圍較寬；而Fe元素含量在不同樣本間變化較小。

圖1　不同苗齡礦物元素含量的變化

Kaiser- Meyer- Olkin KMO檢驗和Bartletts’s球形檢驗用于分析變量間的相關系數(shù)和偏相關系數(shù)，以及變量間獨立性的指標，以便判斷是否適合進行變量主成分分析。對川產(chǎn)楨楠進行KMO和Bartlett's檢驗（表2），KMO值為0.549，大于0.5，表明變量間具有相關性，可以進行后續(xù)相關性分析。另外，Bartlett’s球形檢驗值（167.335）顯著性Sig. 為0.000，達到了顯著水平，說明楨楠元素間存在共性，而并非獨立，也進一步證明了川產(chǎn)楨楠元素間可運用主成分分析法進行分析。

表2　川產(chǎn)楨楠樣本礦物元素的KMO與Bartalett’s檢驗

2.2.2主成分分析

主成分分析可以解釋事物本質之間的聯(lián)系，其重要指標為特征根和方差貢獻率。由表3可知，占總方差87.007%的貢獻率來自前4個主成分，因此選擇前4個主成分進行進一步分析，其中整體共同度的平均值為0.817（表4）。前4個主成分對Ca元素的共同度最大，為0.929，其次是對Mo元素的0.891和K元素的0.843，對Sr (0.710)元素共同度最小。說明前4個主成分對所測定的9個元素的共同性較高，足以進行主成分分析，但是各主成分間載荷系數(shù)相差不大，不易對各成分進行明確，因此需要進行因子載荷矩陣的旋轉，旋轉后各載荷系數(shù)見表6（采用方差最大化旋轉）。

表3　各主成分的特征值和貢獻率

表4　各變量在主成分中的載荷系數(shù)

由表5可知，旋轉后各系數(shù)發(fā)生了明顯極化，第一成分中起主要作用的是Ca元素，第二主成分中表現(xiàn)出與Mn元素呈高度的正相關，第三成分中起主要作用的是Mo元素，第四主成分中表現(xiàn)出與K元素有高度的正相關。

表5　旋轉后各原始變量對主成分的載荷表

2.2.3樣本相似性和聚類分析

運用Pearson相似性對不同地區(qū)的楨楠基于礦物元素的含量進行相似性分析（表6），樣品2:5、3:5的Pearson相關系數(shù)大于0.8；樣品1:4、2:3、4:5 的Pearson相關系數(shù)大于0.6，樣品2:9、3:4、5:9的Pearson相關系數(shù)大于0.5。余弦相關系數(shù)（表7）和Pearson相關系數(shù)均相近，表明樣品可分為2類，1、2、3、4、5、9樣本構成相近歸為第一類，其余樣本可歸為第二類。

表6　 Pearson相關系數(shù)

表7　余弦相關系數(shù)

對不同產(chǎn)地的楨楠材料的礦物元素含量的平均值作為變量，采用類平均法進行系統(tǒng)聚類分析，得到川產(chǎn)楨楠基于礦物元素含量的聚類樹狀圖（圖2），不同的聚類距離下會出現(xiàn)不同的聚類結果。當聚類距離為0時，每個楨楠樣本各自聚集為1類，當聚類距離擴大時，不同樣品聚類才依次分開，當聚類距離達到25時，不同樣品被聚集為1類。

當選擇歐式距離dij=15為分級界限時，供試川產(chǎn)楨楠可分為3個類別，6和11樣本劃分為一類，7、8和10樣本劃分為一類，剩余6份樣本劃分為另一類；當選擇歐式距離dij=20為分級界限時，供試川產(chǎn)楨楠可分為2個類別，1、2、3、4、5、9樣本劃分為一類，6、7、8、10、11樣本劃分為另一類，與Pearson和夾角余弦相似性得到的結果相似?？傮w看來，四川楨楠分布中，川西北和西部等地楨楠聚集為一類，川西南和川南等地楨楠聚集為另一類。

圖2　川產(chǎn)楨楠的聚類分析圖

3　　討論

3.1微量元素在生物體內含量極少，但其與植物的生長息息相關，在植物體內具有重要的生物學功能，含量過多或過少都會導致植物不正常的生理生化變化[20 ]。目前微量元素與植物生長關系的研究是一大熱點，同時也將為植物的營養(yǎng)調控提供重要的理論基礎和指導意義。研究表明，微量元素鐵參與植物葉綠素形成，同時對植物呼吸作用和光合作用產(chǎn)生影響[21- 24]。同時有研究表明，礦物元素之間存在著相互影響現(xiàn)象[25]。楨楠前50年左右材積生長緩慢，之后才進入快速生長期[ 4 ]，其生長速度與礦物元素，特別是微量元素之間限制關系的研究至今沒有開展。通過楨楠不同礦物元素積累呈現(xiàn)的不同規(guī)律，可以探索出礦物元素與楨楠生長緩慢生長的關系，為楨楠的高效栽培提供科學指導。本研究對楨楠幼苗中9中礦物元素進行測定，初步分析表明Ca和S元素含量波動相對較大，一定程度上說明該兩種元素可能是限制楨楠幼苗生長緩慢的制約礦物元素，在實際生產(chǎn)中應注意該元素的跟蹤分析，并及時降低該元素對楨楠生長的制約。

3.2基于礦物元素因子分析和聚類分析的方法已經(jīng)在植物中得到了廣泛的應用[ 26 - 27 ]，本研究對不同產(chǎn)地楨楠木材中9種礦物元素進行測定，并利用SPSS軟件進行分析，結果表明，不同產(chǎn)地和生長環(huán)境中的楨楠木材中所含有的礦物元素存在差別，這種地域特征可能來源于生長環(huán)境、種植方式、大氣環(huán)境、人類活動等多方面[17 ]。通過聚類分析表明，川西北和西部等地楨楠聚集為一類，川西南和川南等地楨楠聚集為另一類，這為楨楠木材的鑒定和區(qū)分提供了指導。

3.3本文初步探索了礦物元素對楨楠生長的影響，并運用礦物元素指標進行了楨楠聚類分析的研究，為其指紋圖譜建立提供了參考。鑒于楨楠分布廣泛和生活環(huán)境的復雜多樣性，采集難度較大，加之生長緩慢，造成實驗存在一定的局限性?？赏ㄟ^測定楨楠中存在的其他礦物元素，以及采取更多的樣本量，來進一步確定數(shù)據(jù)的可靠性。這將更加有助于分析楨楠生長和礦物元素之間的內在關系，打破礦物元素對楨楠生長的限制，為加快楨楠生長速度和人工林發(fā)展進程提供指導。同時借助分子生物學方法，挖掘內在相關的調控機制，為楨楠品種改良和速生品種的培育奠定基礎。

參考文獻

［1］蔣紅星.木中金子——楠木［J］.國土綠化, 2000, 26(5): 41- 41.

［2］蔡春軼,黃建輝.四川都江堰地區(qū)楨楠林、杉木林和常綠闊葉林土壤N庫的季節(jié)變化［J］.生態(tài)學報, 2006, 26(8): 2 540 -2548.

［3］王琦,李因剛,柳新紅,等.湖北恩施楨楠林群落組成與結構［J］.林業(yè)科學研究, 2013, 26(1): 21 - 28.

［4］張煒,何興炳,唐森強,等.四川楨楠生長特性與分布［J］.林業(yè)科技開發(fā), 2012, 26(5): 38 - 41.

［5］龍漢利,張煒,宋鵬,等.四川楨楠生長初步分析［J］.四川林業(yè)科技, 2011, 32(4): 89 - 91.

［6］孫瑾,王曉靜,王飛,等.滇西南地區(qū)樟科17種木材解剖特征研究［J］.華南農(nóng)業(yè)大學學報, 2014, 35(5): 102 - 107.

［7］莊琳,黃群,徐燕紅.楠屬和潤楠屬4種木材的紅外光譜鑒別［J］.福建林業(yè)科技, 2014, 41(4): 21 - 25.

［8］劉勇.我國苗木培育理論與技術進展［J］.世界林業(yè)研究, 2000, 13(5): 43 - 49.

［9］張書娜,王慶成,郝龍飛,等.光照和施肥對白樺林冠下水曲柳、胡桃楸苗木生長的影響［J］.森林工程, 2015, 31(2): 51- 56.

［10］賀維,張煒,胡庭興,等.施肥對1年生楨楠幼苗生長及主要土壤肥力指標的影響［J］.應用與環(huán)境生物學報, 2015, 21(1): 129 - 138

［11］賀維,胡庭興,王銳,等.施肥對楨楠幼苗光合生理及生長特性的影響［J］.西北植物學報, 2014, 34(6): 1 187 - 1 197.

［12］Podio N S, Baroni MV, Badini R G, et al. Elemental and isotopic fingerprint of argentinean wheat. matching soil, water, and crop composition to differentiate provenance［J］. J Agr Food Chem, 2013, 61(16): 3 763 - 3 773.

［13］Deng J W, Fan CL, Jiang ZJ, et al. Qualityassessment and origin tracing of Guangdong Liangcha granules using direct mass spectrometry fingerprinting［J］. Anal Methods, 2012, 4 (11): 3 638 - 3 645.

［14］Anderson K A, Smith B W. Chemical profiling to differentiate geographic growing origins of coffee［J］. J Agr Food Chem, 2002, 50(7): 2 068 - 2 075.

［15］Laursen K H, Schjoerring J K, Olesen J E, et al. Multielemental fingerprintingas a tool for authentication of organic wheat, barley, faba bean, and potato［J］. J Agr Food Chem, 2011, 59(9): 4 385 - 4 396.

［16］沈濤,賈琳,楊亞麗,等.云南薯蕷屬植物礦質元素指紋圖譜研究［J］.熱帶作物學報, 2014, 35(12): 2 332 - 2 339.

［17］錢麗麗,李平惠,楊義杰,等.不同產(chǎn)地蕓豆中礦物元素的因子分析與聚類分析［J］.食品科學, 2015, 36(14):102- 106.

［18］史秀紅,常璇,袁毅,等.不同產(chǎn)地枸杞微量元素的因子分析與聚類分析［J］.安徽農(nóng)業(yè)科學, 2010, 38(4):1839- 1841.

［19］FARNHAMI M, JOHANNESSONKH, SINGH AK, et al. Factor analytical approaches for evaluating groundwater trace element chemistrydata［J］.Analytica Chimica Acta,2003,490:123- 138.

［20］Marschner, H., 2012. In: Marschner, P. (Ed.), Marschner's Mineral Nutrition of HigherPlants. , 3rd ed. Academic press, London, ISBN978- 0- 12- 384905- 2.

［21］Zekri, M., Obreza, T.A., 2003. Micronutrient Deficiencies in Citrus: Iron, Zinc, andManganese. University of Florida Cooperative Extension Service, Institute of Food and Agricultural Sciences, EDIS.

［22］Zekri, M., Obreza, T.A., 2009. Plant Nutrients for Citrus Trees. University of FloridaCooperative Extension Service, Institute of Food and Agricultural Sciences, EDIS.

［23］Srivastava A K, Singh S. Zinc Nutrition in 'Nagpur' Mandarin on Haplustert［J］. Journal of Plant Nutrition, 2009, 32(7): 1 065 - 1 081.

［24］He C, Wang L, Liu J, et al. Evidence for 'silicon' within the cell walls of suspension- cultured rice cells [J]. The New phytologist, 2013, 200(3): 700 - 709.

［25］Zhang Y, Hu C- X, Tan Q- L, et al. Plant nutrition status, yield and qualityof satsuma mandarin (Citrus unshiu Marc.) under soil application of Fe- EDDHA and combination with zinc and manganese in calcareous soil［J］. Scientia Horticulturae, 2014, 174(4): 46 - 53.

［26］Pilgrim T S, Watling R J, Grice K. Application of trace element and stable isotope signatures to determine the provenance of tea (Camellia sinensis) samples［J］. Food Chemistry, 2010, 118(4): 921 - 926.

［27］Chudzinska M, Baralkiewicz D. Estimation of honey authenticity by multielements characteristics using inductively coupled plasma- mass spectrometry (ICP- MS) combined with chemometrics［J］. Food and chemical toxicology : an international journal published for the British Industrial Biological Research Association, 2010, 48(1): 284 - 290.

（責任編輯：張亞楠）

Factor Analysis and Clustering Analysis of Mineral Elements in Sichuan Phoebe zhennan

SHENGYuzhen
（Sichuan Academy of Forestry, Chengdu 610081）

Abstract The content of mineral elements of Sichuan Phoebe zhennan in different seedling age and areas was determined by using inductively coupled plasmatorch-atomic emission spectrometry. The factor analysis and clustering analysis in SPSS were applied to analyze and compare the determination results. Ca and S were richened on different seedling age of Phoebe zhennan, and the content of other element were steady. The factor analysis results showed that the first four principal components accounted for above 87% of the total variance, which indicated that Ca, Mn, Mo and K were the characteristic mineral elements. The results of clustering analysis showed that Phoebe zhennan of northwest Sichuan and western Sichuan were culstered into one group, southwest Sichuan and south Sichuan as the second group. This paper preliminarily analyzes the relationship between mineral elements and seedling growth, verified the feasibility of the method and lay the foundation for the development and utilization of Phoebe zhennan resources in different areas.

Key words Sichuan Phoebe zhennan；Mineral Elements；Clustering Analysis

收稿日期：2015 - 12 - 19

通訊作者：莊國慶(1977- )，男，副研究院，主要從事林學方面研究工作。

作者簡介：第1盛玉珍（1989-），女，碩士，研究方向是林產(chǎn)品加工及檢測。

文章編號：1001 - 9499（2016）01 - 0016 - 05

中圖分類號：S792. 24, S718. 55+4. 2

文獻標識碼：A