雍子涵,張栩,譚天偉
?
糠醛對(duì)粘紅酵母生長(zhǎng)與油脂積累的影響
雍子涵,張栩,譚天偉
北京化工大學(xué)生命科學(xué)與技術(shù)學(xué)院國家能源生物煉制研發(fā)中心,北京 100029
雍子涵, 張栩, 譚天偉. 糠醛對(duì)粘紅酵母生長(zhǎng)與油脂積累的影響. 生物工程學(xué)報(bào), 2015, 31(10): 1484–1491.Yong ZH, Zhang X, Tan TW. Effects of furfural on the growth and lipid production of oleaginous yeast Rhodotorula glutinis. Chin J Biotech, 2015, 31(10): 1484–1491.
為了探究纖維素水解液中常見的發(fā)酵抑制物糠醛對(duì)粘紅酵母生長(zhǎng)與油脂積累的影響,對(duì)比了不同的糠醛濃度 (0.1、0.4、0.6、1.5 g/L) 下粘紅酵母的生物量和油脂積累情況,并探究了1.0 g/L的糠醛對(duì)粘紅酵母不同碳源 (葡萄糖和木糖) 利用的影響。研究表明,當(dāng)糠醛濃度達(dá)1.5 g/L時(shí),粘紅酵母的延遲期延長(zhǎng)至96 h,殘?zhí)歉哌_(dá)17.7 g/L,生物量最高6.6 g/L,僅為正常積累量的47%,油脂含量也減少了約50%;以木糖為碳源時(shí),糠醛對(duì)粘紅酵母的抑制程度小于葡萄糖為碳源時(shí)的情況;在糠醛存在的逆境中,粘紅酵母傾向于生成更多的18碳脂肪酸或18碳不飽和脂肪酸。
糠醛,粘紅酵母,葡萄糖,木糖,生物量,油脂
在面臨能源短缺與環(huán)境污染這兩個(gè)巨大挑戰(zhàn)的今天,生物質(zhì)資源受到了國際社會(huì)的廣泛關(guān)注[1]。我國作為一個(gè)傳統(tǒng)的農(nóng)業(yè)大國,纖維素原料非常豐富,僅小麥秸稈,每年產(chǎn)量就可達(dá)1.1億t[2]。纖維素水解液中含有豐富的可發(fā)酵性糖,尤其是葡萄糖和木糖,可用于微生物產(chǎn)油脂[3-4]。微生物油脂在脂肪酸組成上與植物油非常相似,是生產(chǎn)工業(yè)用油和生物柴油的潛在原料[5]。據(jù)文獻(xiàn)報(bào)道,粘紅酵母屬產(chǎn)油菌株中油脂含量較高的,其油脂含量可達(dá)72%[1,6]。
然而,纖維素原料在預(yù)處理過程中會(huì)產(chǎn)生以乙酸、糠醛、五羥甲基糠醛等為代表的多種發(fā)酵抑制物阻礙微生物的正常生長(zhǎng)[7-10]。前期的研究多集中在以釀酒酵母為代表的多種脫毒方法和一些相關(guān)機(jī)理[11-14],近年來也有學(xué)者探究了纖維素水解液培養(yǎng)發(fā)酵性絲孢酵母[15-18]、彎曲隱球菌[19]以及圓紅冬孢酵母[20]等產(chǎn)油酵母的情況,而關(guān)于纖維素水解液發(fā)酵抑制物糠醛對(duì)產(chǎn)油酵母粘紅酵母生長(zhǎng)的影響則鮮見報(bào)道。
本文在此基礎(chǔ)上著重探究了纖維素水解液發(fā)酵抑制物之一糠醛對(duì)粘紅酵母生長(zhǎng)和油脂積累的影響,以及對(duì)葡萄糖和木糖利用的影響,從而為優(yōu)化纖維素水解液培養(yǎng)基提供實(shí)驗(yàn)支持,并為深入探究糠醛對(duì)產(chǎn)油微生物的抑制機(jī)理和基因改良抗性高產(chǎn)油脂酵母提供科學(xué)依據(jù)。
1.1.1 菌種
本實(shí)驗(yàn)所用菌種粘紅酵母(CGMCC No. 2258) 由北京化工大學(xué)、國家能源生物煉制研發(fā)中心保藏,系從中國食品發(fā)酵所購買的粘紅酵母菌株經(jīng)過多次誘變 所得。
1.1.2 主要儀器
全波長(zhǎng)酶標(biāo)掃描儀 (Thermo Labsystems Multiskan Spectrum)、搖床、臺(tái)式離心機(jī)、高效液相色譜儀 (Thermo Fisher Sciencetific, Utimate 3 000, USA)、電子天平、電熱恒溫水浴鍋、烘箱、紫外超凈臺(tái)、pH計(jì)、自動(dòng)蒸汽滅菌鍋、氣相色譜儀。
1.1.3 培養(yǎng)基
葡萄糖基礎(chǔ)培養(yǎng)基 (g/L):葡萄糖40,KH2PO47,Na2SO42,(NH4)2SO42,MgSO4·7H2O 1.5,酵母粉1.5。
木糖基礎(chǔ)培養(yǎng)基 (g/L):木糖34,KH2PO47,Na2SO42,(NH4)2SO42,MgSO4·7H2O 1.5,酵母粉1.5。
混合糖基礎(chǔ)培養(yǎng)基 (g/L):葡萄糖20,木糖17,KH2PO47,Na2SO42,(NH4)2SO42,MgSO4·7H2O 1.5,酵母粉1.5。
根據(jù)實(shí)驗(yàn)設(shè)置,在相應(yīng)基礎(chǔ)培養(yǎng)基中添加糠醛至相應(yīng)濃度。
1.2.1 培養(yǎng)方法
本實(shí)驗(yàn)采用500 mL搖瓶,1/5裝液量,10%種子液接種量,搖床培養(yǎng)條件為200 r/min、30 ℃。
1.2.2 分析方法
采用高效液相色譜儀測(cè)定葡萄糖和木糖的條件為:Aminex HPX-87H柱 (300 mm′7.8 mm, Bio-Rad Corp., USA),65 ℃示差檢測(cè)器,流動(dòng)相為 0.005 mol/L H2SO4,流速為0.6 mL/min。糠醛的測(cè)定條件為:C18有機(jī)酸柱 (250 mm′4.6 mm, Thermo Fisher Corp., USA),35 ℃,紫外檢測(cè)器,220 nm, 流動(dòng)相為水與甲醇的混合液 (75∶25,/),流速為0.5 mL/min。
生物量:將培養(yǎng)液稀釋一定濃度后,在波長(zhǎng)為600 nm下測(cè)其吸光度,通過標(biāo)準(zhǔn)曲線計(jì)算其生物量。
油脂產(chǎn)量采用磷酸香草醛法測(cè)定,參照文獻(xiàn)[21]。
采用氯仿-甲醇法提取出油脂后 (參照文獻(xiàn)[22]),采用氣相測(cè)定油脂組分,參照文獻(xiàn)[23]。
添加了糠醛的培養(yǎng)基中粘紅酵母生長(zhǎng)速率明顯下降,且糠醛濃度越高,抑制現(xiàn)象越明顯 (圖1A)。結(jié)合糖耗曲線與糠醛濃度變化曲線 (圖1B和圖1C),可以發(fā)現(xiàn),當(dāng)糠醛濃度降為0時(shí),粘紅酵母才開始大量利用葡萄糖快速生長(zhǎng)。當(dāng)初始糠醛濃度達(dá)到1.5 g/L時(shí),粘紅酵母在96 h之后才開始大量生長(zhǎng),即對(duì)照組的酵母已經(jīng)由對(duì)數(shù)期轉(zhuǎn)入穩(wěn)定期時(shí),1.5 g/L糠醛處理的酵母才開始由延遲期進(jìn)入對(duì)數(shù)期。從最大生物量積累來看,添加了糠醛的培養(yǎng)基中生物量最高僅8 g/L (糠醛濃度為0.1 g/L時(shí)),約占正常培養(yǎng)基生物量 (14 g/L) 的57%。值得一提的是,不同糠醛濃度下,粘紅酵母最終的生物量積累幾乎相當(dāng)。另外,在初始糠醛濃度為0.1 g/L的培養(yǎng)基中,前48 h的生物量積累情況與對(duì)照組不相上下。這說明,糠醛濃度越低,粘紅酵母的延遲期越短。0.1 g/L的糠醛前期對(duì)粘紅酵母的延遲期影響并不大,但其生長(zhǎng)活力受到限制,主要表現(xiàn)為48 h 之后,生物量便不再積累,進(jìn)入穩(wěn)定期,而對(duì)照中的生物量仍繼續(xù)增加。
就葡萄糖利用情況而言 (圖1B),糠醛濃度越高,葡萄糖消耗速度越慢,而且所消耗糖量也越來越少。當(dāng)糠醛濃度達(dá)到1.5 g/L時(shí),直到第4天,粘紅酵母才開始大量利用葡萄糖,但持續(xù)時(shí)間僅48 h左右,發(fā)酵至216 h時(shí)殘?zhí)侨愿哌_(dá)17.7 g/L。
從糠醛濃度變化上看 (圖1C),24 h之后,0.1 g/L、0.4 g/L、0.6 g/L實(shí)驗(yàn)組的糠醛濃度幾乎降為0,120 h之后,1.5 g/L實(shí)驗(yàn)組的糠醛濃度也幾乎降為0。結(jié)合生物量積累 (圖1A) 和糖耗情況 (圖1B),不難發(fā)現(xiàn),培養(yǎng)基中粘紅酵母的生長(zhǎng)速率和糖耗速率與糠醛的濃度變化密切相關(guān)。
從油脂產(chǎn)量和油脂含量來看,糠醛濃度越高,其對(duì)油脂合成的抑制程度越大 (圖1D和1E)。當(dāng)糠醛濃度達(dá)1.5 g/L時(shí),粘紅酵母油脂產(chǎn)量遠(yuǎn)遠(yuǎn)低于對(duì)照組中的油脂產(chǎn)量,而油脂含量與未添加糠醛的基礎(chǔ)培養(yǎng)基情況相比,降低了約50%。
綜上分析,糠醛對(duì)粘紅酵母在葡萄糖培養(yǎng)基中的生長(zhǎng)和油脂積累均存在著明顯的抑制作用,并且抑制程度與糠醛濃度呈正相關(guān)??啡舛冉禐?之后粘紅酵母才開始大量利用葡萄糖進(jìn)行生長(zhǎng),但即使是0.1 g/L的初始糠醛濃度下,其殘?zhí)且哺哌_(dá)12.4 g/L。當(dāng)糠醛濃度達(dá) 1.5 g/L時(shí),粘紅酵母的延遲期延長(zhǎng)至96 h,殘?zhí)歉哌_(dá)17.7 g/L,生物量積累最高6.6 g/L,僅為正常積累量的47%,油脂含量也減少了約50%,糠醛的抑制作用已經(jīng)相當(dāng)嚴(yán)重了。
基于2.1的研究結(jié)果,選取1.0 g/L的糠醛濃度,探究了糠醛對(duì)粘紅酵母不同單一碳源利用的影響。從生物量積累來看 (圖2A),發(fā)酵前24 h,添加糠醛的不同碳源培養(yǎng)基中,酵母幾乎都沒有生長(zhǎng),處于適應(yīng)調(diào)整階段,24 h之后才開始大量生長(zhǎng),進(jìn)入對(duì)數(shù)生長(zhǎng)期,48 h后便進(jìn)入穩(wěn)定生長(zhǎng)期,并且持續(xù)相當(dāng)長(zhǎng)一段時(shí)間。整體比較,各培養(yǎng)基最大生物量表現(xiàn)為,對(duì)照組:葡萄糖(13.7 g/L) >混合糖(12.2 g/L) >木糖(11.6 g/L),處理組:混合糖~葡萄糖(5.6 g/L) >木糖(4.8 g/L)。處理組中等比例混合碳源與單一葡萄糖 (六碳糖) 培養(yǎng)基中的粘紅酵母生長(zhǎng)情況相當(dāng),這說明糠醛對(duì)木糖利用的抑制程度小于對(duì)葡萄糖利用的抑制程度。
從油脂積累來看 (圖2B),對(duì)照組在穩(wěn)定期的中后期 (96?120 h)油脂逐漸達(dá)到最高值,而處理組在前24 h油脂幾乎沒有積累,這與其在此期間生物量沒有積累的情況相一致,24?48 h油脂逐漸積累到最大值,48 h之后油脂產(chǎn)量保持穩(wěn)定,這也與其在48 h之后生物量保持穩(wěn)定的情況相一致??傮w比較,對(duì)照組培養(yǎng)基中的最大油脂積累情況為:葡萄糖(2.5 g/L) >混合糖(2.3 g/L) >木糖(1.5 g/L),處理組培養(yǎng)基中的最大油脂積累情況為:混合糖(0.7 g/L) >葡萄糖 (0.45 g/L) >木糖(0.33 g/L),這也在一定程度上說明了木糖為碳源時(shí)糠醛對(duì)粘紅酵母的抑制程度小于對(duì)葡萄糖為碳源時(shí)的抑制程度。
從葡萄糖消耗情況來看 (圖2C),添加糠醛的葡萄糖培養(yǎng)基與混合糖培養(yǎng)基中,前24 h葡萄糖幾乎都沒有消耗,24 h之后才開始大量消耗,但是72 h之后葡萄糖濃度便幾乎不再變化,葡萄糖利用率僅為30%?35%。而對(duì)照組中,葡萄糖則表現(xiàn)為持續(xù)消耗,直至最終消耗完畢。
圖2D為木糖濃度的變化情況,對(duì)照組中,木糖培養(yǎng)基中的木糖濃度表現(xiàn)為持續(xù)下降,其中前72 h下降較快,72 h之后相對(duì)緩慢。混合糖培養(yǎng)基中木糖在前48 h幾乎不變化,72 h之后才開始大量消耗,這是因?yàn)榛旌咸桥囵B(yǎng)基中的葡萄糖濃度在72 h時(shí)已消耗完畢,在葡萄糖消耗完之后粘紅酵母才開始大量利用木糖。處理組中,木糖培養(yǎng)基中的木糖濃度在前24 h幾乎沒有變化,24 h之后才開始消耗木糖,但72 h之后便幾乎不再利用木糖,木糖濃度穩(wěn)定在 22 g/L左右,木糖利用率僅21%?31%?;旌咸桥囵B(yǎng)基中的木糖幾乎沒有消耗,這與葡萄糖沒有消耗完的情況相一致,因?yàn)樵谄咸烟呛湍咎枪餐嬖跁r(shí),粘紅酵母在利用完葡萄糖之后才會(huì)開始大量利用木糖。
由圖3可知,在不同碳源的培養(yǎng)基中,粘紅酵母所產(chǎn)的主要脂肪酸種類是一致的,分別為16碳的軟脂酸 (16∶0)和18碳的硬脂酸 (18∶0)、油酸 (18∶1)、亞油酸 (18∶2),但相應(yīng)的比例有所差別。
基礎(chǔ)培養(yǎng)基中,葡萄糖為碳源時(shí)粘紅酵母傾向于積累大量的油酸 (68%),其次是軟脂酸 (16%)、亞油酸 (11%)和硬脂酸 (5%)。木糖為碳源時(shí),粘紅酵母傾向于積累大量的硬脂酸 (44%),其次是油酸 (34%)、亞油酸 (13%) 和軟脂酸 (8%)?;旌咸菫樘荚磿r(shí),相應(yīng)的油脂含量則介于二者之間。
添加糠醛后,葡萄糖為碳源時(shí)粘紅酵母傾向于生成更多的18碳脂肪酸 (其中油酸含量仍然是最高的),16碳的軟脂酸含量有所下降 (約占11%)。木糖為碳源時(shí)情況則相反,其軟脂酸含量有所上升 (約占15%),硬脂酸含量明顯下降 (約占25%),其他的成分含量均上升了,其中油酸含量上升為最大 (約占41%)?;旌咸菫樘荚磿r(shí),變化幅度介于二者之間,整體變化趨勢(shì)與葡萄糖為碳源時(shí)一致。
圖3 糠醛對(duì)不同碳源下粘紅酵母所產(chǎn)油脂的主要成分相對(duì)含量影響
總體而言,添加糠醛后,不同碳源的培養(yǎng)基中,粘紅酵母產(chǎn)生的主要脂肪酸不僅種類不變,相應(yīng)的含量分配也趨于一致,由高到低分別為:油酸、硬脂酸、亞油酸、軟脂酸,18碳的脂肪酸占85%?92%,其中油酸的比重占41%?55%。這說明,在糠醛存在的逆境中,粘紅酵母傾向于生成更多的長(zhǎng)鏈脂肪酸或長(zhǎng)鏈不飽和脂肪酸以抵御逆境。這可能與糠醛的抑制機(jī)理相關(guān),具體原因還有待進(jìn)一步探究。
綜上分析,糠醛明顯抑制粘紅酵母的生長(zhǎng)和油脂積累,并且其抑制程度與糠醛濃度呈正相關(guān)??啡?duì)粘紅酵母的抑制作用主要體現(xiàn)在延遲期的延長(zhǎng)和對(duì)數(shù)生長(zhǎng)期的縮短、生物量積累的減少、油脂積累的減少、糖利用能力下降等多方面??啡?duì)葡萄糖利用的抑制程度大于對(duì)木糖利用的抑制程度。從油脂的主要成分與比例分配上來看,糠醛存在時(shí),粘紅酵母傾向于生成更多的18碳長(zhǎng)鏈脂肪酸或18碳長(zhǎng)鏈不飽和脂肪酸以抵御逆境。關(guān)于糠醛對(duì)粘紅酵母的抑制機(jī)制還有待進(jìn)一步探究。
[1] Li Q. Study on the Lipid Synthesis and Pentose Sugar Metabolism of Oleaginous Yeast. Beijing: Beijing University of Chemical Technology, 2007 (in Chinese).李塹. 產(chǎn)油酵母油脂合成與戊糖代謝的研究. 北京: 北京化工大學(xué), 2007.
[2] Yong DY, Pang YZ, Yuan HR, et al.Enhancing biogas production from anaerobically digested wheat straw through ammonia pretreatment. Chin J Chem Eng, 2014, 22(5): 576?582.
[3] Huang C, Zong M, Wu H, et al. Microbial oil production from rice straw hydrolysate by. Bioresour Technol, 2009, 100(19): 4535?4538.
[4] Huang C, Chen X, Xiong L, et al. Microbial oil production from corncob acid hydrolysate by oleaginous yeast.2013, 49: 273?278.
[5] Xing X. Study of sugar contained wastewater for microbial lipid production[D]. Beijing: Beijing University of Chemical Technology, 2009 (in Chinese). 邢旭. 含糖廢水發(fā)酵生產(chǎn)微生物油脂[D]. 北京:北京化工大學(xué), 2009.
[6] Pan LX, Yang DF, Shao L, et al. Screening, identification and fermention characters of oleaginous yeasts by utilizing sugar cane bagasse hemicellulose hydrolysates. Microbiol China, 2008, 35(8): 1176?1181 (in Chinese).潘麗霞, 楊登峰, 邵利, 等.利用蔗渣半纖維素水解液產(chǎn)油酵母的篩選、鑒定及發(fā)酵實(shí)驗(yàn). 微生物學(xué)通報(bào), 2008, 35(8): 1176?1181.
[7] Wang XF, Zhang J, Xin XJ, et al. Furfural degradation by filamentous fungusZN1. Chin J Biotech, 2012, 28(9): 1070?1079 (in Chinese).王曉鳳, 張建, 辛秀娟, 等. 絲狀真菌ZN1的糠醛降解代謝分析. 生物工程學(xué)報(bào), 2012, 28(9): 1070?1079.
[8] Palmqvist E, Hahn-H?gerdal B. Fermentation of lignocellulosic hydrolysates. I: inhibition and detoxification.resour Technol, 2000, 74(1): 17?24.
[9] Wang D, Wang HH, Wang J, et al. Effects of furfural and 5-hydroxymethylfurfural on succinic acid production byChin J Biotech, 2013, 29(10): 1463?1472 (in Chinese).王丹, 王洪輝, 王競(jìng), 等. 糠醛和5-羥甲基糠醛對(duì)大腸桿菌產(chǎn)丁二酸的影響. 生物工程學(xué)報(bào), 2013, 29(10): 1463?1472.
[10] Wu J, Cheng KK, Li WY, et al. Effect of acetic acid, furfural and 5-hydroxymethylfurfural onproduction of 2,3-butanediol by. Chin J Biotech, 2013, 29(3): 350?357 (in Chinese).吳晶, 程可可, 李文英, 等.乙酸、糠醛和5-羥甲基糠醛對(duì)產(chǎn)酸克雷伯氏菌發(fā)酵生產(chǎn)2,3-丁二醇的影響. 生物工程學(xué)報(bào), 2013, 29(3): 350?357.
[11] Almeida JRM, Modig T, Petersson A, et al. Increased tolerance and conversion of inhibitors in lignocellulosic hydrolysates by. J Chem Technol Biot, 2007, 82(4): 340?349.
[12] Palmqvist E, Hahn-H?gerdal B. Fermentation of lignocellulosic hydrolysates II: inhibitors and mechanisms of inhibition. Bioresour Technol, 2000, 74(1): 25?33.
[13] Taherzadeh MJ, Gustafsson L, Niklasson C, et al. Conversion of furfural in aerobic and anaerobic batch fermentation of glucose by. J Biosci Bioeng, 1999, 87(2): 169?174.
[14] Delgenes JP, Moletta R, Navarro JM. Effects of lignocellulose degradation products on ethanol fermentations of glucose and xylose by,,, and. Enzyme Microb Technol, 1996, 19(3): 220?225.
[15] Huang C, Wu H, Li R, et al. Improving lipid production from bagasse hydrolysate withby response surface methodology. New Biotechnol, 2012, 29(3): 372?378.
[16] Huang C, Wu H, Smith TJ, et al.detoxification of furfural during lipid production by the oleaginous yeast. Biotechnol Lett, 2012, 34(9): 1637?1642.
[17] Huang C, Wen P, Wu H, et al. Evaluating the effect of acetic acid, furfural and catechol on the growth and lipid accumulation ofby response surface methodology. Afr J Microbiol Res, 2012, 6(24): 5266?5275.
[18] Huang C, Wu H, Liu Q, et al. Effects of aldehydes on the growth and lipid accumulation of oleaginous yeast. J Agric Food Chem, 2011, 59(9): 4606?4613.
[19] Yu X, Zeng J, Zheng Y, et al. Effect of lignocellulose degradation products on microbial biomass and lipid production by the oleaginous yeast. Process Biochem, 2014, 49(3): 457?465.
[20] Zhao X, Peng F, Du W, et al. Effects of some inhibitors on the growth and lipid accumulation of oleaginous yeastand preparation of biodiesel by enzymatic transesterification of the lipid. Bioproc Biosyst Eng, 2012, 35(6): 993?1004.
[21] Izard J, Limberger RJ. Rapid screening method for quantitation of bacterial cell lipids from whole cells. J Microbiol Meth, 2003, 55(2): 411?418.
[22] Bligh EG, Dyer WJ. A rapid method of total lipid extraction and purification. Can J Biochem Phys, 1959, 37(8): 911?917.
[23] Li Z, Jiang F, Li Y, et al. Simultaneously concentrating and pretreating of microalgaespp. by three-phase partitioning. Bioresour Technol, 2013, 149: 286?291.
(本文責(zé)編 陳宏宇)
Effects of furfural on the growth and lipid production of oleaginous yeast
Zihan Yong, Xu Zhang, and Tianwei Tan
,,,100029,
In order to illustrate the effects of furfural, one of the most common inhibitory compounds in lignocellulosic hydrolysate, on oleaginous yeast, we investigated the effects of different concentrations of furfural (0.1, 0.4, 0.6 and 1.5 g/L) on the biomass and lipid production of, as well as the effects of 1.0 g/L furfural on the utilization of glucose and xylose. Results showed that: when the furfural concentration reached 1.5 g/L, the lag phrase time was extended to 96 h, and the residual glucose was up to 17.7 g/L, with maximum biomass of only 6.6 g/L, which accounted for 47% of that in the basic medium (furfural-free), and the corresponding lipid content was reduced about 50%. Furfural showed lighter inhibitory degree onwhen xylose acted as the carbon source than glucose was the carbon source; more C18 fatty acids or unsaturated C18 fatty acids were generated in the presence of furfural.
furfural,,glucose, xylose, biomass, lipid
10.13345/j.cjb.140555
November 17, 2014; Accepted:January 30, 2015
National High Technology Research and Development Program of China (863 Program) (Nos. 2013AA050702, 2012AA022304).
Xu Zhang. Tel: +86-10-64450593; Fax: +86-10-64416428; E-mail:zhangxu@mail.buct.edu.cn
國家高技術(shù)研究發(fā)展計(jì)劃 (863計(jì)劃) (Nos. 2013AA050702, 2012AA022304) 資助。
2015-03-20
http://www.cnki.net/kcms/detail/11.1998.Q.20150320.1200.001.html