熱補(bǔ)針法干預(yù)風(fēng)濕性關(guān)節(jié)炎寒證模型家兔尿液代謝組學(xué)研究

袁博，王金海，張星華，田亮，任朝展，杜小正

1.甘肅中醫(yī)藥大學(xué)，甘肅 蘭州 730000；2.蘭州大學(xué)第二醫(yī)院，甘肅 蘭州 730030；3.甘肅省中醫(yī)院，甘肅 蘭州 730050

熱補(bǔ)針法干預(yù)風(fēng)濕性關(guān)節(jié)炎寒證模型家兔尿液代謝組學(xué)研究

袁博1，王金海2，張星華1，田亮1，任朝展3，杜小正1

1.甘肅中醫(yī)藥大學(xué)，甘肅 蘭州 730000；2.蘭州大學(xué)第二醫(yī)院，甘肅 蘭州 730030；3.甘肅省中醫(yī)院，甘肅 蘭州 730050

目的通過對風(fēng)濕性關(guān)節(jié)炎（RA）寒證模型家兔尿液中內(nèi)源性代謝物的分析，探討熱補(bǔ)針法治療RA的特異性機(jī)制。方法采用卵蛋白誘導(dǎo)聯(lián)合低溫冷凍建立RA寒證家兔模型。將40只清潔級青紫藍(lán)家兔隨機(jī)分為正常組、模型組、平補(bǔ)平瀉組、捻轉(zhuǎn)補(bǔ)法組、熱補(bǔ)針法組，每組8只。分別應(yīng)用平補(bǔ)平瀉、捻轉(zhuǎn)補(bǔ)法、熱補(bǔ)針法針刺“足三里”，1次/d，留針30 min，共7 d。干預(yù)結(jié)束后收集新鮮尿液，氣相色譜-四級桿飛行時間質(zhì)譜采集各組尿液代謝物圖譜，利用主成分分析和偏最小二乘判別分析對數(shù)據(jù)進(jìn)行統(tǒng)計。結(jié)果與正常組比較，模型組尿液代謝物變化顯著，氨基酸代謝相關(guān)物質(zhì)降低，三羧酸循環(huán)、糖類代謝、脂肪酸代謝相關(guān)物質(zhì)增高；與模型組比較，各針刺組三羧酸循環(huán)、脂質(zhì)代謝相關(guān)物質(zhì)降低（P＜0.05）；熱補(bǔ)針法組三羧酸循環(huán)相關(guān)物質(zhì)明顯低于平補(bǔ)平瀉組和捻轉(zhuǎn)補(bǔ)法組（P＜0.05）。結(jié)論熱補(bǔ)針法治療RA的特異性體現(xiàn)在對三羧酸循環(huán)的調(diào)控。

熱補(bǔ)針法；風(fēng)濕性關(guān)節(jié)炎；代謝組學(xué)；系統(tǒng)生物學(xué)；氣相色譜-四級桿飛行時間質(zhì)譜；兔

風(fēng)濕性關(guān)節(jié)炎（rheumatoid arthritis，RA）是一種自身免疫性疾病。目前尚無根治方法，針灸具有抗炎、免疫調(diào)節(jié)等作用[1-3]，對RA具有確切的臨床療效[4]。熱補(bǔ)針法是鄭魁山教授簡化燒山火手法和進(jìn)火法而形成的一種復(fù)式手法，具有補(bǔ)益正氣、溫陽散寒功效。課題組前期研究表明，熱補(bǔ)針法在治療RA過程中具有相對特異性的效果[5-8]。但前期研究僅從鎮(zhèn)痛方面進(jìn)行了初步分析，缺乏對針刺后機(jī)體產(chǎn)生的全部生物信息的綜合分析。

代謝組學(xué)是一種能夠反映生物系統(tǒng)中所有低分子量（＜1000）代謝物質(zhì)的定性和定量分析手段，通過考察生物體系終末代謝產(chǎn)物的變化，研究生物體系代謝途徑[9-10]。氣相色譜-質(zhì)譜（GC-MS）聯(lián)用技術(shù)是最早應(yīng)用于代謝組學(xué)的色質(zhì)聯(lián)用技術(shù)，除具有較高的分辨率和檢測靈敏度外，還具有用來參考比較的標(biāo)準(zhǔn)代謝物譜庫。因此，該研究手段的優(yōu)勢與針灸作用特點(diǎn)甚相吻合。為了進(jìn)一步探討熱補(bǔ)針法治療RA的特異性機(jī)制，本研究以RA模型家兔為載體，采用氣相色譜-四極桿飛行時間質(zhì)譜（GC-Q/TOF-MS）從系統(tǒng)生物學(xué)角度探討熱補(bǔ)針法治療RA的特異性機(jī)制。

1　材料與方法

1.1動物

2～3月齡清潔級健康青紫藍(lán)家兔40只，雌雄各半，體質(zhì)量（2.5±0.5）kg，蘭州生物制品研究所，許可證號SCXK（甘）2012-0001。飼養(yǎng)于溫度（20±2）℃，相對濕度（50±10）%，每日光照時間12 h，自由進(jìn)食飲水。

1.2針具、主要試劑與儀器

針灸針（華成牌，蘇州東邦醫(yī)療器械有限公司）。卵蛋白干粉（批號A5253-250G）、弗氏完全佐劑（批號F5881），Sigma公司。甲醇、吡啶，國藥集團(tuán)；尿素酶、甲氧胺鹽酸吡，Sigma公司。Agilent 7890A/5975C GC/MS系統(tǒng)、毛細(xì)管色譜柱（HP-5 ms Agilent），J&W Scientific公司。

1.3分組及造模

家兔適應(yīng)性飼養(yǎng)7 d，隨機(jī)分為正常組、模型組、平補(bǔ)平瀉組、捻轉(zhuǎn)補(bǔ)法組、熱補(bǔ)針法組，每組8只。采用卵蛋白誘導(dǎo)聯(lián)合低溫冷凍進(jìn)行RA寒證模型復(fù)制[11-12]。實驗正式開始前3 d，用電推剃去家兔背部肩胛骨間的毛。將4 mg/mL卵蛋白溶液與等量弗氏完全佐劑混勻，充分震蕩成乳化劑，2%碘酒、75%乙醇消毒肩背部均勻選定的6個點(diǎn)，每點(diǎn)皮下注射0.2 mL。14 d后以相同方式、相同劑量重復(fù)皮下注射1次。第2次免疫后6 d，剃去家兔雙后腿膝關(guān)節(jié)周圍的毛，2%碘酒、75%乙醇消毒雙膝關(guān)節(jié)，關(guān)節(jié)內(nèi)注入20 mg/mL卵蛋白生理鹽水溶液0.4 mL。24 h內(nèi)膝關(guān)節(jié)出現(xiàn)紅腫觸痛，表明模型復(fù)制成功。成模后第2日，乙醇消毒兩后腿足，上、下、左、右圍置冰袋（3份冰＋1份結(jié)晶氯化鈣，粉碎混合，溫度降至-20～-25 ℃冷凍1.5 h，45 ℃溫水復(fù)溫5 min（室溫12 ℃)，共冷凍1次。冷凍后14 d，家兔表現(xiàn)為精神萎糜不振，兩后肢匍伏跛行，膝關(guān)節(jié)腫脹、周徑明顯大于對照組，局部青紫、皮溫不高，表示成功建立RA寒證模型。

1.4針刺方法

先將家兔固定，選取后腿“足三里”[13]。各針刺組分別施以相應(yīng)針法。熱補(bǔ)針法操作參照《鄭氏針灸全集》[14]：左手拇指緊按穴位，右手持針刺入穴內(nèi)8～12 mm，針刺得氣后，左手加重壓力，右手拇指向前連續(xù)捻按5次；針下沉緊后針尖拉著有感應(yīng)的部位，連續(xù)重插輕提5次；拇指再向前連續(xù)捻按5次，反復(fù)操作1 min；最后使針下沉緊，留針30 min。平補(bǔ)平瀉法、捻轉(zhuǎn)補(bǔ)法操作參照《刺法灸法學(xué)》[15]。平補(bǔ)平瀉法：針刺得氣后，前后均勻捻轉(zhuǎn)，不單向捻轉(zhuǎn)，指力均勻，角度180°～360°，頻率60～90 r/min，反復(fù)操作1 min。捻轉(zhuǎn)補(bǔ)法：針刺得氣后，小幅度前后捻轉(zhuǎn)針體（角度180°～360°，頻率60～90 r/min），拇指向前左轉(zhuǎn)時用力重、指力沉重向下，拇指向后右轉(zhuǎn)還原時用力輕，反復(fù)操作1 min。每日1次，每次留針30 min，連續(xù)7 d。

1.5樣本采集

針刺結(jié)束后，收集各組新鮮尿液，8000 r/min離心10 min，用EP管收集上層尿液，離心分裝后置于-80 ℃冰箱保存?zhèn)溆谩?/p>

1.6樣品衍生化處理

尿液常溫解凍后，混勻，取200 μL，12 000 r/min離心10 min，取上清液50 μL，加尿素酶10 μL，混勻，37 ℃反應(yīng)15 min；加200 μL甲醇，轉(zhuǎn)入進(jìn)樣瓶，置溫和氮?dú)庀麓蹈桑患尤?0 μL甲氧胺鹽酸吡啶（20 mg/mL）溶液，強(qiáng)烈震蕩30 s，37 ℃反應(yīng)90 min；取出，加入30 μL BSTFA（含1%TMCS，BSTFA∶TMCS＝99∶1）的衍生試劑，70 ℃反應(yīng)60 min。取出樣本，室溫放置30 min備檢。

1.7氣相色譜-質(zhì)譜條件

進(jìn)樣口溫度280 ℃，EI離子源溫度230 ℃，四極桿溫度150 ℃，高純氦氣（純度＞99.999%）作為載氣，不分流進(jìn)樣，進(jìn)樣量1.0 μL。升溫程序為：初始溫度80 ℃，維持2 min，10 ℃/min速度升至320 ℃，維持6 min。采用全掃描模式進(jìn)行質(zhì)譜檢測，質(zhì)譜檢測范圍為50～550 m/z。采用隨機(jī)順序進(jìn)行連續(xù)樣本分析，避免儀器信號波動所造成的影響。

1.8統(tǒng)計學(xué)方法

原始數(shù)據(jù)采用XCMS軟件進(jìn)行特征離子峰提取、峰對齊、色譜峰識別及峰匹配，然后構(gòu)建多維數(shù)據(jù)集，采用SIMCA-P13.0進(jìn)行多維數(shù)據(jù)分析，數(shù)據(jù)預(yù)處理采用Par的方式，數(shù)據(jù)分析方法包括主成分分析（PCA）和偏最小二乘法判別分析（PLS-DA），最后獲得差異物質(zhì)。P＜0.05表示差異有統(tǒng)計學(xué)意義。

2　結(jié)果

2.1原始色譜圖的可視化檢查

對所有樣本的總離子流色譜圖（TIC）進(jìn)行重疊與匹配分析（見圖1），結(jié)果表明，質(zhì)譜峰形明顯、界限清晰、重現(xiàn)性好，與標(biāo)準(zhǔn)品譜峰一致，5組樣品質(zhì)譜峰重疊性良好。

2.2模式識別

各組PCA圖見圖2。為了最大化反映各組別間代謝模式的差異，采用監(jiān)督性PLS-DA，見圖3。從圖3中可以看出各組明顯區(qū)別，針刺干預(yù)后向正常組回歸，并且熱補(bǔ)針法組更突出。據(jù)PLS-DA模型第一主成分VIP值（閾值＞1）并同時結(jié)合P值（P＜0.05）來尋找差異性表達(dá)代謝物。通過以上分析，最終獲得的生物標(biāo)記物見表1。可以得出，模型組尿液代謝物變化較正常組主要體現(xiàn)為氨基酸代謝相關(guān)物質(zhì)降低，三羧酸循環(huán)、糖類代謝、脂肪酸代謝增強(qiáng)。各針刺組較模型組主要體現(xiàn)為三羧酸循環(huán)代謝及脂肪酸代謝降低，且熱補(bǔ)針法組三羧酸循環(huán)明顯于低平補(bǔ)平瀉組和捻轉(zhuǎn)補(bǔ)法組。

圖1　樣本TIC（重疊）

圖2　各組尿液PCA圖

圖3　各組尿液PLS-DA圖

表1　各組間潛在生物標(biāo)志物及代謝途徑

3　討論

RA屬中醫(yī)學(xué)“痹證”范疇。正氣虛弱是RA發(fā)病的內(nèi)因，風(fēng)寒濕邪侵犯是外因，邪氣乘虛侵入四肢關(guān)節(jié)筋骨，陽氣虛不能驅(qū)邪外出，邪氣黏滯纏綿，閉阻經(jīng)絡(luò)，發(fā)為此病。熱補(bǔ)針法是鄭魁山教授臨床治療RA的經(jīng)驗結(jié)晶，具有補(bǔ)益正氣、溫陽散寒的功用，與RA的病機(jī)特點(diǎn)相吻合。現(xiàn)代醫(yī)學(xué)研究發(fā)現(xiàn)，RA的發(fā)生與三羧酸循環(huán)、氨基酸代謝、脂肪酸代謝、葡萄糖代謝等密切相關(guān)[16-19]，其涉及的相關(guān)代謝物有檸檬酸、琥珀酸、肉堿、亮氨酸、麥芽糖、乳糖、嘧啶核苷酸、尿酸等。說明RA的發(fā)病機(jī)制具有途徑廣、層次多等復(fù)雜的網(wǎng)絡(luò)系統(tǒng)特點(diǎn)。因此，從現(xiàn)代醫(yī)學(xué)角度而言，其發(fā)病機(jī)制與針灸多靶點(diǎn)、多途徑、多層次、類似于網(wǎng)絡(luò)結(jié)構(gòu)的作用特點(diǎn)相符。

本研究發(fā)現(xiàn)，模型組中琥珀酸、檸檬酸、異檸檬酸、草酰乙酸、α-酮戊二酸、丙酮酸的量異常升高，以上這些物質(zhì)均是三羧酸循環(huán)的中間產(chǎn)物。三羧酸循環(huán)反應(yīng)在線粒體內(nèi)發(fā)生，線粒體是細(xì)胞內(nèi)氧化磷酸化和合成三磷酸腺苷（ATP）的主要場所，它為細(xì)胞的活動提供了能量，這些代謝物的排泄增加可能會引發(fā)能量的供應(yīng)不足。有研究顯示，RA關(guān)節(jié)腔缺氧[20]，細(xì)胞適應(yīng)缺氧主要是從氧化磷酸化向糖酵解途徑轉(zhuǎn)換以產(chǎn)生ATP，模型組中乳糖參與糖代謝，是機(jī)體內(nèi)的能量物質(zhì)，其生成的半乳糖磷酸酯可以為糖酵解途徑提供能量，尿液中乳糖排泄量的異常說明RA體內(nèi)糖酵解途徑發(fā)生改變而導(dǎo)致能量代謝通路出現(xiàn)紊亂。因此，我們的研究證實了這一結(jié)論。模型組中L-亮氨酸水平下降，它屬于支鏈亮氨酸，對于促進(jìn)肌肉蛋白質(zhì)的合成、骨豁肌微細(xì)損傷的修復(fù)、糖異生以及骨豁肌葡萄糖的攝取都起著重要的作用[21]；模型組參與多不飽和脂肪酸代謝的2-羥基十六烷酸、前列腺素E2甲酯、癸二酸、壬二酸的水平異常增高，以上這些物質(zhì)可促使炎癥因子腫瘤壞死因子-α（TNF-α）、白細(xì)胞介素（IL）-6的產(chǎn)生，TNF-α是RA發(fā)病機(jī)制中關(guān)鍵的促炎因子，可刺激軟骨細(xì)胞合成IL-6與IL-8，產(chǎn)生進(jìn)一步的炎性反應(yīng)，與IL-6、IL-1β等協(xié)同共同促進(jìn)破骨細(xì)胞的生成導(dǎo)致骨質(zhì)的破壞[22]，TNF-α又可增加一氧化氮、前列腺素E的生成和骨的吸收[23]。綜上所述，這些代謝途徑的改變可以解釋RA機(jī)體全身炎性反應(yīng)狀態(tài)導(dǎo)致對能量的需求增加及供給的相對不足，經(jīng)針刺干預(yù)后，參與脂肪酸代謝的上述代謝物均顯著降低（P＜0.05），各針刺組之間比較差異無統(tǒng)計學(xué)意義（P＞0.05）。針刺干預(yù)后參與三羧酸循環(huán)的上述物質(zhì)普遍下調(diào)，尤以熱補(bǔ)針法組琥珀酸、α-酮戊二酸、草酰乙酸下調(diào)最顯著（P＜0.05），說明熱補(bǔ)針法可以調(diào)節(jié)RA異常的能量代謝通路。

平補(bǔ)平瀉法、捻轉(zhuǎn)補(bǔ)法及熱補(bǔ)針法均可使RA異常的三羧酸循環(huán)、脂肪酸代謝發(fā)生逆轉(zhuǎn)，提示這些針刺手法可調(diào)控RA體內(nèi)三羧酸循環(huán)（能量代謝）、脂肪酸代謝（抗炎效應(yīng)）通路，而熱補(bǔ)針法在調(diào)控三羧酸循環(huán)方面均明顯優(yōu)于平補(bǔ)平瀉、捻轉(zhuǎn)補(bǔ)法。基于此，我們認(rèn)為調(diào)控三羧酸循環(huán)（能量代謝）可能是熱補(bǔ)針法治療RA的特異性機(jī)制。

[1] ZHU J, CHEN X Y, LI L B, et al. Electroacupuncture attenuates collagen-induced arthritis in rats through vasoactive intestinal peptide signalling-dependent re-establishment of the regulatory T cell/T-helper 17 cell balance[J]. Acupuncture in Medicine,2015：acupmed-2014-010732.

[2] 李文迅,黃怡然,雷林丹,等.電針、手針“足三里”對慢性關(guān)節(jié)炎大鼠血清炎性細(xì)胞因子含量的影響[J].針刺研究,2012,37(4)：271-276.

[3] LI J, LI J, CHEN R, et al. Targeting NF-κΒ and TNF-α activation by electroacupuncture to suppress collagen-induced rheumatoid arthritis in model rats[J]. Altern Ther Health Med,2015,21(4)：26-34.

[4] OUYANG B S, GAO J, CHE J L, et al. Effect of electro- acupuncture on tumor necrosis factor-alpha and vascular endothelial growth factor in peripheral blood and joint synovia of patients with rheumatoid arthritis[J]. Chin J Integr Med,2011,17(7)：505-509. [5] 張星華,朱博雯,楊晨光,等.熱補(bǔ)針法治療風(fēng)寒濕阻型風(fēng)濕性關(guān)節(jié)炎的臨床療效[J].遼寧中醫(yī)雜志,2015,42(5)：1076-1078.

[6] 杜小正,方曉麗,東貴榮.熱補(bǔ)針法對關(guān)節(jié)炎兔的鎮(zhèn)痛后效應(yīng)及其腦脊液中β-EP和CCK-8含量的影響[J].中醫(yī)研究,2011,24(6)：18-20.

[7] 杜小正,秦曉光,方曉麗.熱補(bǔ)針法鎮(zhèn)痛后效應(yīng)的觀察及其對外周PGE-2,SP含量變化的影響[J].中醫(yī)研究,2010,23(1)：18-21.

[8] 杜小正,秦曉光,方曉麗.熱補(bǔ)針法鎮(zhèn)痛后效應(yīng)及其對關(guān)節(jié)局部組織β-內(nèi)啡肽和前列腺素E2的影響[J].針刺研究,2009,34(5)：319-323.

[9] HORNING E C, HORNING M G. Metabolic prfiles：gas-phase methods for analysis of metabolites[J]. Clin Chem,1971,17(8)：802-809.

[10] KANG J, ZHU L, LU J, et al. Application of metabolomics in autoimmune diseases：Insight into biomarkers and pathology[J]. J Neuroimmunol,2015,279：25-32.

[11] TOMINGA K, ALSTERGREN P, KURITA H, et al. Clinical course of an antigen-induced arthritis model in the rabbit temporomandibular joint[J]. J Oral Pathol Med,1999,28(6)：268-273.

[12] 吳墾莉,張珊珊.寒凝血瘀證動物模型的研制[J].中國中醫(yī)基礎(chǔ)醫(yī)學(xué)雜志,1996,2(2)：49-51.

[13] 李忠仁.實驗針灸學(xué)[M].北京：中國中醫(yī)藥出版社,2007：242-247.

[14] 鄭魁山,鄭俊江,鄭俊鵬,等.鄭氏針灸全集[M].北京：人民衛(wèi)生出版社,2000：289-290.

[15] 陸壽康.刺法灸法學(xué)[M].北京：中國中醫(yī)藥出版社,2003：92-95.

[16] YUE R, ZHAO L, HU Y, et al. Metabolomic study of collageninduced arthritis in rats and the interventional effects of Huang-Lian-Jie-Du-Tang, a traditional Chinese medicine[J]. J Ethnopharmacol,2013,145(2)：465-475.

[17] YANG X Y, DI Z K, LIN K, et al. Energy metabolism disorder as a contributing factor of rheumatoid arthritis：acomparative proteomic and metabolomic study[J]. PLoS one,2015,10(7)：e0132695.

[18] SPIES C M, STRAUB R H, BUTTGEREIT F. Energy metabolism and rheumatic diseases：from cell to organism [J]. Arthritis Res Ther,2012,14(3)：216-226.

[19] CHIMENTI M S, TUCCI P, CANDI E, et al. Metabolic profiling of human CD4+cells following treatment with methotrexate and anti-TNF-α infliximab[J]. Cell Cycle,2013,12(18)：3025-3036.

[20] SIVAKUMAR B, AKHAVANI M A, WINLOVE C P, et al. Synovial hypoxia as a cause of tendon rupture in rheumatoid arthritis[J]. J Hand Surg Am,2008,33(1)：49-58.

[21] 張慶,漆正堂,丁樹哲.亮氨酸及其代謝產(chǎn)物促進(jìn)蛋白質(zhì)合成的分子機(jī)制及其補(bǔ)充效果[J].生命的化學(xué),2009,29(1)：103-107.

[22] SCOTT D L, WOLFE F, HUIZINGA T W. Rheumatoid arthritis[J]. Lancet,2010,376(9746)：1094-1108.

[23] RALSTON S H, GRABOWSKI P S. Mechanisms of cytokine induced bone resorption：role of nitric oxide, cyclic guanosine monophosphate, and prostaglandins[J]. Bone,1996,19(1)：29-33.

Study on Intervention of Heat-reinforcing Needling for Metabonomics of Rabbits Urine with Rheumatoid Arthritis Cold Syndrome

YUAN Bo1, WANG Jin-hai2, ZHANG Xing-hua1, TIAN Liang1,REN Chao-zhan3, DU Xiao-zheng1
(1. Gansu University of Chinese Medicine, Lanzhou 730000, China; 2. The Second Hospital of Lanzhou University, Lanzhou 730030, China; 3. Gansu Provincial Hospital of Traditional Chinese Medicine, Lanzhou 730050, China)

Objective To analyze the endogenous metabolites in urine of rabbits with rheumatoid arthritis (RA) cold syndrome; To investigate the specificity mechanisms of heat-reinforcing needling for RA. Methods Egg albumin inducing method combined with low temperature and freezing method was used to establish RA cold syndrome rabbit models. A total of 40 healthy purple blue rabbits were randomly divided into normal control group (NC), model control (MC), reinforcing-reducing needling group (RRN), twirling-reinforcing needling group (TRN) and heat-reinforcing needling group (HRN) (n=8/group). Except NC and MC, RRN was given acupuncture of reinforcing-reducing needling at Zusanli (ST36), TRN was given acupuncture of twirling-reinforcing needling at Zusanli, and HRN was administrated acupuncture of heat-reinforcing needling at Zusanli, once a day and retaining 30 min, a total of seven days. Fresh urine sample after the intervention was collected and then gas chromatography quadrupole time-of-flight mass spectrometry technology were used to evaluate metabolic profiles. All the data were analyzed by principal component analysis and partial least squares discriminant analysis. Results Compared with NC, the urine metabolites of MC changed significantly, mainly reflecting in the increase of amino acid metabolism were decreased, metabolites of TCA cycle, carbohydrate and fatty acid metabolism were increased; Compared with MC, the urine metabolites of TCA cycle and fatty acid metabolism were decreased in all intervention of acupuncture group (P＜0.05), but the metabolites of TCA cycle in HRN were obviously less than that of RRN and TRN (P＜0.05). Conclusion The specificity of heat-reinforcing needling for RA manifests in regulation of TCA cycle.

heat-reinforcing needling; rheumatoid arthritis; metabonomics; systems biology; GC-Q/TOF-MS; rabbits

10.3969/j.issn.1005-5304.2016.12.016

R245

A

1005-5304(2016)12-0064-05

國家自然科學(xué)基金（81260558）；甘肅省自然科學(xué)基金（1208RJZA185）；甘肅省高校基本科研業(yè)務(wù)費(fèi)項目（2011年）

杜小正，E-mail：lz-duxiaozheng@163.com

（2016-01-23）

（2016-02-20；編輯：華強(qiáng)）