大多数植物病原体是革兰氏阴性细菌[1], 其自身代谢产物N-酰基高丝氨酸内酯(N-Acylhomoserine lactones, AHLs)已被证明是参与群体感应(Quorum sensing, QS)系统中关键的信号分子。细菌能通过QS系统感知菌群密度并调节其生理活动[2], 当菌群中AHLs浓度达到一定阈值后, 就能激活病原菌致病基因表达, 使动植物出现致病性, 如参与病原体和植物之间相互作用的毒力因子[3-4]。因此, 科学家提出了一种植物病原体控制策略, 即破坏病原体的QS系统可能下调菌群密度, 从而抑制可能感染宿主植物的毒力因子的表达。
当前对N-酰基高丝氨酸内酯酶(N-acylhomos-erine Lactonase, AiiA蛋白)编码基因aiiA的研究大部分侧重于对其进行克隆表达。苏云金芽胞杆菌aiiA编码的AiiA蛋白, 通过水解内酯键失活AHLs, 从而破坏AHLs介导的QS系统[5-6]。有研究表明, 在微生物中, 蜡状芽孢杆菌AiiA蛋白对魔芋软腐病菌CZY表现出较强的抗病活性[7]; Álvaro等[8]分离到一株苏云金芽胞杆菌, 克隆其aiiA并转入大肠杆菌, 发现重组的大肠杆菌在含有欧文氏杆菌的培养基中共培养时, 可以干扰AHLs信号的积累或响应, 破坏其QS系统, 这与Wu等[9]从苏云金芽胞杆菌中克隆到的aiiA在毕赤酵母GS115中表达, 其提取物表现出对欧文氏杆菌的抗菌活性有相似之处。同时, 有学者将aiiA转化到巨尾桉中后, 用青枯雷尔氏菌侵染植株后, 发现转基因的巨尾桉出现萎蔫现象的时间显著延迟, 发病率降低, 表明转入aiiA能增强其抗青枯病能力[10]。同样, 将在巴斯德毕赤酵母中异源表达的AiiA蛋白与嗜水气单胞菌的同时注射进鲤鱼时, 其死亡率与单独注射嗜水气单胞菌相比, 降低了约25%[11]。因此, 使用AiiA蛋白破坏QS系统被认为是生物防治病害的有效手段。然而, aiiA在苏云金芽胞杆菌中的表达非常低, 对aiiA启动子的转录调控研究很少。因此, 研究aiiA的转录机制可能是增强AiiA蛋白表达量的有效方法, 启动子区域的寻找是研究转录机制的第一步。
本研究通过克隆aiiA的5'侧翼区域, 利用在线生物软件分析潜在的启动子区域。以gfp作为报告基因, 通过分段克隆及定点突变, 对aiiA的5'端侧翼序列的启动子功能进行研究。分析苏云金芽胞杆菌aiiA启动子的转录调控功能, 为aiiA在不同宿主中的高表达提供理论基础, 也为细菌病害的生物防治提供新的思路。
1 材料与方法 1.1 材料大肠杆菌E.coli DH5α、E.coli BL21(DE3)由福建师范大学生命科学学院工程研究中心保存, 苏云金芽胞杆菌BRC-ZLL5(Bacillus thuringiensis)由福建农林大学生命科学学院关雄教授惠赠。苏云金芽胞杆菌aiiA由福建师范大学生命科学学院工程研究中心前期克隆得到。所有菌株均用LB培养基培养。绿色荧光蛋白表达质粒pET28a-gfp以及由pET28a-gfp质粒敲除T7启动子改造而来的绿色荧光蛋白本底表达质粒pET28a-dP-gfp由福建师范大学发育生物学重点实验室保存, pMD18-T载体购自TaKaRa公司。
1.2 方法 1.2.1 苏云金芽胞杆菌aiiA的5'端侧翼序列的克隆将克隆的aiiA(GenBank登录号:DQ440581.1)在NCBI网站进行BLAST检索, 选择同源性最高的苏云金芽胞杆菌全基因组序列(GenBank登录号:CP001903.1), 在aiiA起始密码子ATG上游1 900 bp处设计特异性PCR扩增引物, 上游引物P1插入BglⅡ酶切位点:5'-GAAGATCTAGCTAGTGTAAAATGACGTTTGCTT-3';下游引物P2插入NheⅠ酶切位点:5'-CTAGCTAGCTATGTATCCACCTTTACATTTAGTT-3'(下划线处碱基为限制性内切酶识别位点)。此外, 在aiiA起始密码子ATG下游200 bp处设计另一条下游引物P3:5'-GTACCGTTAAAAAGCCCTTCATTAT-3', 可以克隆出一条既包含aiiA的5'端侧翼序列, 又包含200 bp aiiA的目的序列。本研究引物均由上海生工生物工程股份有限公司合成。
为扩增aiiA的5'端侧翼, 使用分子克隆方法提取苏云金芽胞杆菌的质粒[12], 进行PCR扩增, 扩增体系为模板DNA 1 μL、dNTP 0.5 μL、10×buffer 1.5 μL、引物各0.5μL、Ex Taq酶0.15 μL和ddH2O 10.9 μL。反应条件为94℃ 5 min; 94℃ 30 s, 55℃ 30 s, 72℃ 2 min, 30个循环; 72℃ 10 min。PCR扩增产物经1%琼脂糖凝胶电泳后纯化回收、连接和转化。筛选阳性克隆测序。插入aiiA的5'端侧翼区的质粒pMD18-T命名为pMD18-T-aP-1930-+ 200。
1.2.2 苏云金芽胞杆菌aiiA的5'端侧翼序列aP-1930-+200启动子区域分析使用在线软件BDGP Neural Network Promoter Prediction, V2.2(http://www.fruitfly.org/seq_tools/promoter.html)、Promoter 2.0 Prediction Server(http://www.cbs.dtu.dk/services/Promoter/)和NCBI中Finding Promoter数据库预测aiiA的5'端侧翼序列可能的启动子区域。
1.2.3 苏云金芽胞杆菌aiiA的5'端侧翼序列的分段克隆以克隆的侧翼序列aP-1930-+200为模板, 根据启动子潜在区域预测结果, 在aiiA的5'端侧翼序列的不同位置设计10对引物(表 1), 进行分段克隆。上游引物加入BglⅡ酶切位点, 下游引物加入NheⅠ酶切位点。
以pMD18-T-aP-1930-+200质粒为模板进行PCR扩增。扩增体系同1.2.1。其中片段aP-100--1、aP-295--1、aP-429--1、aP-429--296和aP-295--101的PCR反应条件为94℃ 5 min, 94℃ 30 s, 55℃ 30 s, 72℃ 30 s, 30个循环; 72℃ 10 min; 片段aP-531--1和aP-1930--1430的PCR反应条件为94℃ 5 min; 94℃ 30 s, 55℃ 30 s, 72℃ 1 min, 30个循环; 72℃ 10 min; 片段aP-1031--1、aP-1531--1和aP-1931-1的PCR反应条件为94℃ 5 min; 94℃ 30 s, 55℃ 30 s, 72℃ 2 min, 30个循环; 72℃ 10 min。
1.2.4 启动子功能鉴定载体pET28a-aP-x-gfp的构建将aiiA的5'端侧翼序列的分段克隆与经BglⅡ和NheⅠ双酶切的pET28a-gfp质粒用T4 DNA连接酶16℃过夜连接, 然后将连接的pET28a-aP-x-gfp载体转化到大肠杆菌E.coli BL21(DE3)感受态细胞, 并命名为E.coli BL21(DE3)-pET28a-aP-x-gfp。
1.2.5 aiiA的5'端侧翼区域的启动子活性分析将重组大肠杆菌E.coli BL21(DE3)-pET28a-aP-x-gfp在50 mL含有50 μg/mL卡那霉素的LB培养基中37℃震荡孵育过夜。当细菌浓度达到OD600=0.6-0.8时, 向培养基中加入IPTG(终浓度1.0 mmol/L)。取5 μL菌液制成玻片标本后, 直接用荧光显微镜观察(激发光为蓝光), 通过荧光强度比较启动子的活性。大肠杆菌E.coli BL21(DE3)-pET28a-gfp、E.coli BL21(DE3)-pET28a-dP-gfp分别作为阳性对照和阴性对照。取3.5 mL上述诱导过的菌液于荧光比色皿中, 利用荧光分光光度计在激发光波长(Ex)为491 nm、发色光波长(Em)为509 nm下测量菌液荧光强度。
1.2.6 苏云金芽胞杆菌aiiA基因启动子的定点突变在载体pET28a-aP-295--101-gfp基础上, 对aiiA启动子-150--142 bp的2个TATA序列进行点突变, 结合待突变位点设计4对引物(表 2):前2对为-144和-147位点处的单点突变引物, 第3对为-144位和-147位点处的双点突变引物, 第4对为2个TATA序列缺失突变引物。
定点突变PCR扩增体系为模板DNA 0.4 μL、dNTP Mixture 1.6 μL、5×Prime STAR buffer 4.0 μL、引物各0.4 μL、Prime STAR HS DNA Polymerase 0.2 μL和ddH2O 13.0 μL。PCR反应条件为94℃ 5 min; 94℃ 30 s, 51℃ 30 s, 72℃ 7 min, 30个循环; 72℃ 10 min。PCR扩增产物经琼脂糖凝胶电泳验证正确后, 冰浴2-3 min, 向反应物中加入1 μL DpnⅠ酶, 消化模板质粒pET28a-aP-295--101-gfp, 37℃水浴过夜后热击转化大肠杆菌E.coli BL21(DE3)感受态细胞。
通过DNA测序鉴定阳性诱变菌株, 通过荧光显微镜和荧光分光光度法分析含有突变的5'侧翼突变区域的重组菌株。
2 结果 2.1 aiiA的5'端侧翼序列的克隆与分析以苏云金芽胞杆菌质粒DNA为模板, 用引物P1和P2、P1和P3进行扩增, 获得约1.9和2.1 kb大小片段。将2.1 kb片段测序结果进行在线BLAST(http://www.ncbi.nlm.nih.gov), 结果显示克隆的aiiA的5’端侧翼序列与已经公布的苏云金芽胞杆菌YBT-1520(GenBank登录号:CP007607.1)aiiA的5’端侧翼序列同源性为100%, 其中3’端的200 bp与aiiA的5’端同源性为100%(GenBank登录号:DQ000642.1)。所克隆的片段在GenBank登录号为KP690195。
利用在线分析软件BDGP Neural Network Promoter Prediction对aiiA的5’端侧翼序列进行分析, 发现其含有潜在的启动子区域(图 1红色区域), 其中位于-1 395--1 346 bp和-339--292 bp处的启动子区域可靠性最高, 分值均为满分1.0;通过NCBI中Finding Promoter数据库预测, 发现其可能存在2个启动子区域(图 1下划线区域), 分别位于-429--1 bp和-1 930--1 822 bp处。此外, 在序列-150--142 bp处还存在2个连续的TATA序列, 可能对启动子活性有很重要的作用。
2.2 aiiA的5'端侧翼序列的启动子活性鉴定以pMD18-T-aP-1930-+200质粒为模板, 进行PCR扩增。琼脂糖凝胶电泳(图 2)可知在100、295、134、429、195、501、531、1 031、1 531和1 930 bp处出现目的条带, PCR扩增产物与目的片段大小一致, 产物特异性好。
由模板质粒pET28a-dP-gfp构建而来的分段克隆载体pET28a-aP-x-gfp(pET28a-aP-100--1-gfp、pET28a-aP-295--1-gfp、pET28a-aP-295--101-gfp、pET28a-aP-429--1-gfp、pET28a-aP-429--296-gfp、pET28a-aP-531--1-gfp、pET28a-aP-1031--1-gfp、pET28a-aP-1531--1-gfp、pET28a-aP-1930--1-gfp和pET28a-aP-1930--1429-gfp)经过双酶切(BglⅡ和NdeⅠ)和DNA测序鉴定, 10个克隆子均显示了正确的条带, 构建的载体经上海生工生物股份有限公司序列测定后表明插入序列与目的序列相同, 表明载体pET28a-aP-295--101-gfp、pET28a-aP-100--1-gfp、pET28a-aP-429--296-gfp、pET28a-aP-429--1-gfp、pET28a-aP-295--1-gfp、pET28a-aP-531--1-gfp、pET28a-aP-1531--1-gfp、pET28a-aP-1031--1-gfp、pET28a-aP-1930--1429-gfp和pET28a-aP-1930--1-gfp构建成功。
培养含有pET28a-aP-x-gfp载体的重组大肠杆菌菌株, 使用荧光显微镜和荧光分光光度计分析启动子活性。荧光显微镜(图 3)显示大肠杆菌菌株BL21(DE3)-pET28a-aP-1930--1429-gfp, BL21(DE3)-pET28a-aP-295--1-gfp和BL21(DE3)-pET28a-aP-295--101-gfp能够产生绿色荧光, 表明aiiA的5'侧翼区域aP-1930中的aP-1930-1429、aP-295--1和aP-295--101具有启动子活性, 其中aP-1930--1429和aP-295--101区域荧光强度最高。
荧光分光光度法进行定量分析(图 4)显示BL21(DE3)-pET28a-aP-1930--1429-gfp菌株启动子活性最高, 其次是BL21(DE3)-pET28a-aP-295--101-gfp菌株, BL21(DE3)-pET28a-aP-295--1-gfp菌株启动子活性相对较低。
2.3 定点突变对苏云金芽胞菌aiiA的5'端侧翼区域启动子活性的影响以pET28a-aP-295--101-gfp质粒为模板进行定点突变, 测序结果显示缺失突变质粒序列TATA盒被成功去除, 点突变质粒-144位和-147位碱基均由A突变为G, 但双点突变质粒除了预定的-144位点和-147位点发生突变外, -163位点也发生了突变, 由A变为T, 序列分析表明, -163不在启动子潜在区域内。
荧光显微镜观察显示(图 5), 相对于野生型菌株(BL21(DE3)-pET28a-aP-295--101-gfp), 所有的定点突变菌株的启动子活性均显著降低。荧光分光光度计测定结果(图 6)表明, 双点突变(P-295--101(-144、-147 Mu))和2个TATA盒缺失(P-295--101(-116--134dMu))比单点突变(P-295--101(-144 Mu)和P-295--101(-147 Mu))引起的启动子活性降低幅度更大, 提示苏云金芽胞杆菌aiiA启动子-150--142 bp区域的2个TATA盒对aiiA的5'端侧翼区域的启动子活性可能起着很关键的作用。
3 讨论将AiiA蛋白运用于农业生产, 提高aiiA的转录水平是至关重要的一步。有研究表明将苏云金芽胞杆菌cry3Aa的启动子用于启动aiiA, 可以提高aiiA的表达, 增加AiiA蛋白的产量[13]。一些报告基因, 如绿色荧光蛋白, 常用于研究启动子功能活性的鉴定[14-15]。本实验为研究aiiA的启动子活性, 克隆了苏云金芽胞杆菌aiiA的5'侧翼序列(aP-1930--1), 通过荧光显微镜和荧光分光光度计对aiiA的5'端侧翼序列的功能进行定性和定量分析, 根据报告基因GFP的荧光强度分析表明, 分段克隆到的侧翼序列aP-295--101、aP-295--1和aP-1930--1429均具有启动子活性, 这与通过生物信息学分析其启动子潜在区域结果一致, 但序列aP-1930--1不具备启动子功能。此外, 通过在NCBI上BLAST比对发现序列-1 821--1430 bp区域是氨基酸通透性酶家族蛋白(Amino acid permease family protein)基因的编码区, 表明aiiA的5'端-1 930--1 429 bp区不是aiiA启动子, 在苏云金芽胞杆菌中一个启动子同时启动两个基因表达的可能性不大。分段克隆分析显示侧翼序列aP-295--101的启动子活性比aP-295--1强, 而且aP-100--1序列没有启动GFP转录的能力。结果表明, aP-100--1序列可能是负调控序列, 这也可能导致AiiA蛋白表达量低。
定点突变是研究基因和启动子功能的有效方法[16-18], 前期研究显示, 定点突变能够有效的用于蛋白质热稳定性的研究[19]。aiiA定点突变后, 表达的AiiA蛋白其酶活力和温度耐受性均得到了提高[20]。本研究以定点突变技术来研究aiiA启动子活性位点, 结果表明单点突变使得该基因启动子的活性显著降低, 双点突变(P-295--101(-144、-147 Mu))和TATA盒缺失进一步降低其活性。
启动子具有启动和调节基因转录的作用[21-123]。传统原核启动子一般位于转录起始位点上游-10区(TATAAT)和-35区(TTGACA), 然而本研究发现, 位于aiiA转录起始位点上游-100 bp处的序列没有启动子活性(-295--101 bp可能是aiiA的启动子)。有研究表明, 苏云金芽胞杆菌的cryIIIA毒素基因的启动子区域至少由基因上游的-635--553 bp, -553--367 bp以及-367-18 bp 3个结构域组成[24]。可能苏云金芽胞杆菌基因的启动子存在距转录起始位点较远的机制。苏云金芽胞杆菌AiiA蛋白表达较低很可能是由aiiA启动子上游和下游序列对aiiA启动子的负调控作用所致。
当然本文使用大肠杆菌替代原始宿主苏云金芽胞杆菌进行研究也存在一些不足, 首先不同的表达载体其转录因子可能不同, 会对转录调控造成不同的影响; 其次, 用gfp代替aiiA, 可能会由于目的蛋白不同而导致其表达模式的不同, 从而影响对启动子功能的判断。尽管如此, 本研究根据破囊壶菌Δ4-脂肪酸脱饱和酶基因5'端侧翼区域的研究方法对苏云金芽胞杆菌aiiA的5'端侧翼序列克隆并在大肠杆菌中进行分段表达, 其确实可以启动报告基因gfp的表达, 具有启动子功能活性, 为aiiA基因最可能的启动子区域及其重要的转录元件[25]; 但其启动aiiA表达的具体调控机制仍需要进行同源性表达分析验证。该结果有利于后续进一步深入研究该启动子在原始宿主中的功能及其转录调控分子机制, 也有利于对aiiA在大肠杆菌等异源宿主中进行高表达的研究。
4 结论苏云金芽胞杆菌aiiA的5'端-295--1 bp为aiiA的启动子区, -429--296与-100--1序列是aiiA的负调控序列。-150--142 bp区域的2个TATA盒对aiiA的5'端侧翼区域的启动子活性起重要作用。
[1] |
de Kievit TR, Iglewski BH. Bacterial quorum sensing in pathogenic relationships[J]. Infect Immun, 2000, 68(9): 4839-4849. DOI:10.1128/IAI.68.9.4839-4849.2000 |
[2] |
Yeon KM, Cheong WS, et al. Quorum sensing:a new biofouling control paradigm in a membrane bioreactor for advanced wastewater treatment[J]. Environ Sci Technol, 2009, 43(2): 380-385. DOI:10.1021/es8019275 |
[3] |
Zhao CM, Zeng H, et al. N-Acyl homoserine lactonase promotes prevention of Erwinia virulence with Zwittermicin A-producing strain Bacillus cereus[J]. Biotechnol Bioeng, 2008, 100(3): 599-603. DOI:10.1002/(ISSN)1097-0290 |
[4] |
Barnard AML, et al. Quorum sensing, virule-nce and secondary metabolite production in plant soft-rotting bacteria[J]. Philos Trans R Soc Lond B Biol Sci, 2007, 362(1483): 1165-11683. DOI:10.1098/rstb.2007.2042 |
[5] |
Dong YH, Gusti AR, Zhang Q, et al. Identification of quorum-quenching N-acyl homoserine lactonases from Bacillus species[J]. Applied & Environmental Microbiology, 2002, 68(4): 1754-1759. |
[6] |
Lee SJ, Park SY, Lee JJ, et al. Genes encoding the N-acyl homoserine lactone-degrading enzyme are widespread in many subspecies of Bacillus thuringiensis[J]. Applied & Environmental Microbiology, 2002, 68(8): 3919-3924. |
[7] |
周燚, 孙明, 喻子牛. 苏云金芽胞杆菌AiiA蛋白对魔芋软腐病菌的抗病活性[J]. 武汉大学学报:理学版, 2004, 50(6): 761-764. |
[8] |
Álvaro MF, González A, Carmen J, et al. Identification, cloning and lactonase activity of recombinant protein of N -acyl homoserine lact-onase(AiiA)from Bacillus thuringiensis 147-115-16 strain[J]. Revista Colombiana De Biotecnología, 2014, 16(1): 153-162. DOI:10.15446/rev.colomb.biote.v16n1.40495 |
[9] |
Wu J, Jiao Z, et al. Constitutive and secretory expression of the AiiA in Pichia pastoris inhibits Amorphophallus konjac soft rot disease[J]. American Journal of Molecular Biology, 2016, 6(2): 79-87. |
[10] |
Ouyang LJ, Li LM. Effects of an inducible aiiA gene on disease resistance in Eucalyptus urophylla×Eucalyptus grandis[J]. Transgenic Research, 2016, 25(4): 441-452. DOI:10.1007/s11248-016-9940-x |
[11] |
Chen RD, Zhou ZG, Cao YN, et al.High yield expression of an AHL-lactonase from Bacillus sp.B546 in Pichia pastoris and its application to reduce Aeromonas hydrophila mortality in aquaculture[J].Microbial Cell Factories, 2010, 9(1): 39.
|
[12] |
Sambrook J, Fritsch EF, Maniatis T. Molecular cloning[M]. Vol 2. New York: Cold Spring Habor Laboratory Press, 1989.
|
[13] |
朱晨光, 孙明, 喻子牛. 带cry3Aa启动子的aiiA基因在苏云金芽胞杆菌中的表达[J]. 生物工程学报, 2003, 19(4): 397-401. DOI:10.3321/j.issn:1000-3061.2003.04.004 |
[14] |
Huang JZ, Jiang XZ, Xia XF, et al. Cloning and Functional Identification of Delta5 fatty acid desaturase gene and its 5'-upstream region from marine fungus Thraustochytrium sp.FJN-10[J]. Marine Biotechnology, 2011, 13(1): 12-21. DOI:10.1007/s10126-010-9262-6 |
[15] |
Lissemore JL, Jankowski JT, Thomas CB, et al. Green fluorescent protein as a quantitative reporter of relative promoter activity in E.coli[J]. Biotechniques, 2000, 28(1): 88-89. |
[16] |
Traboni C, Ciliberto G, Cortese R. A novel method for site-directed mutagenesis:its application to an eukaryotic tRNAPro gene promoter[J]. Embo Journal, 1982, 1(4): 415-420. DOI:10.1002/embj.1982.1.issue-4 |
[17] |
Ma D, Cook DN, O'Brien DA, et al. Analysis of the promoter and regulatory sequences of an oxygen-regulated bch operon in Rhodobacter capsulatus by site-directed mutagenesis[J]. Journal of Bacteriology, 1993, 175(7): 2037-2045. DOI:10.1128/jb.175.7.2037-2045.1993 |
[18] |
Hernandez K, Fernandez-Lafuente R. Control of protein immobilization:coupling immobilization and site-directed mutagenesis to improve biocatalyst or biosensor performance[J]. Enzyme & Microbial Technology, 2011, 48(2): 107-122. |
[19] |
Ben Mabrouk S, Aghajari N, Ben Ali M, et al. Enhancement of the thermostability of the maltogenic amylase MAUS149 by Gly312Ala and Lys436Arg substitutions[J]. Bioresource Technology, 2011, 102(2): 1740-1746. DOI:10.1016/j.biortech.2010.08.082 |
[20] |
杨梅, 谢盼盼, 简思美, 等. 定点突变提高N-酰基高丝氨酸内酯酶酶活和温度稳定性[J]. 微生物学报, 2014, 54(8): 905-912. |
[21] |
Shen-Orr SS, Milo R, Mangan S, et al. Network motifs in the transcriptional regulation network of Escherichia coli[J]. Nature Genetics, 2002, 31(1): 64-68. DOI:10.1038/ng881 |
[22] |
Elkon R, Linhart C, Sharan R, et al. Genome-wide in silico identification of transcriptional regulators controlling the cell cycle in human cells[J]. Genome Research, 2003, 13(5): 773-780. DOI:10.1101/gr.947203 |
[23] |
Lee TI, Rinaldi NJ, Robert F, et al. Transcriptional regulatory networks in Saccharomyces cerevisiae[J]. Science, 2002, 298(5594): 799-804. DOI:10.1126/science.1075090 |
[24] |
Agaisse H, Lereclus D. Structural and functional analysis of the promoter region involved in full expression of the cryⅢA toxin gene of Bacillus thuringiensis[J]. Mol Microbiol, 1994, 13(1): 97-107. DOI:10.1111/mmi.1994.13.issue-1 |
[25] |
江贤章, 陈金卿, 田宝玉, 等. 破囊壶菌Δ4-脂肪酸脱饱和酶基因5'端侧翼区域的克隆与鉴定[J]. 福建师大学报:自然科学版, 2008, 24(6): 89-94. |