林麝源蜡样芽孢杆菌SCBCM001株的分离鉴定及全基因组序列分析

【背景】细菌性疾病是林麝规模化养殖的重要制约因素,蜡样芽孢杆菌曾在林麝化脓灶中检出,但是目前对林麝源蜡样芽孢杆菌的研究报道很少。【目的】对分离自病死林麝肝脏中的一株疑似蜡样芽孢杆菌进行分离鉴定和全基因组序列分析,为林麝相关疾病的防治奠定基础。【方法】将病原菌纯化培养后,对病原菌进行生化试验、药敏试验和小鼠致病性试验;并通过第3代单分子测序技术进行全基因组测序,根据测序结果进一步评估物种间的亲缘关系,并进行基因功能注释和遗传进化分析。【结果】该病原菌经平均核苷酸相似度分类和系统发育树分析属于蜡样芽孢杆菌群,生化结果符合蜡样芽孢杆菌的一般特征,将分离菌株命名为SCBCM001。该菌株对小鼠的半数致死量为8.3×107 CFU,对大多数β-内酰胺类、四环素和磺胺异噁唑耐药,对氨基糖苷类、头孢氨苄、头孢哌酮和亚胺培南等药物敏感;全基因组测序结果表明该菌株的染色体大小为5292570bp,GC含量为35.37%,多位点序列分型显示该菌株属于ST427序列类型。在菌株SCBCM001基因组内发现hblA、hblC、hblD、nheA、nheB、clo和cytK等多种毒力因子,同时菌株携带对β-内酰...     

一株拮抗黄单胞菌的贝莱斯芽孢杆菌的分离和鉴定

【目的】为了筛选防治水稻条斑病(bacterial leaf streak,BLS)的生防细菌。【方法】以水稻条斑病菌(Xanthomonas oryzae pv. oryzicola,Xoc)的模式菌株RS105为靶标菌,采用平板稀释和抑菌圈法,从空心菜根际土壤中筛选到一株对RS105具有拮抗作用的细菌菌株504。通过形态学、生理生化特征以及16SrDNA和gyrA序列分析对菌株504进行了鉴定。利用牛津杯法测定504对植物病原黄单胞菌的拮抗活性及其无菌发酵液拮抗活性的稳定性。通过PCR扩增预测504编码合成脂肽类和聚酮类化合物的合成相关基因。采用苗期水稻注射接菌法来评价水稻组织中504对Xoc的拮抗活性。【结果】菌株鉴定结果表明504为贝莱斯芽孢杆菌,命名为Bacillusvelezensis504。抑菌实验显示,B.velezensis504对黄单胞菌属的细菌具有较好的抑菌活性,对水稻白叶枯病菌(X. oryzae pv. oryzae,Xoo)的拮抗效果最显著。基因预测结果显示,B. velezensis 504含有fenA、dhbA、sfrA、bmyA、beaS、dfnA及bacA等编码脂肽类和聚酮糖类抑菌化合物的基因簇。其无菌发酵液的活性物质耐高温和蛋白酶降解,但不耐强酸、强碱,在pH值为5.5–8.9时仍具有稳定的拮抗活性。在高感水稻品种原丰早上,B. velezensis 504对Xoc在水稻叶片中引起的水渍症状具有显著的抑制作用。【结论】B. velezensis 504能够特异性拮抗黄单胞菌,在黄单胞菌引起的细菌性病害的生物防治中将具有较大的应用潜力。     

一株产低温蛋白酶地衣芽胞杆菌NWMCC0046全基因组测序及分析

地衣芽胞杆菌(Bacillus licheniformis)NWMCC0046是从西藏日喀则地区屠宰场废弃血污放置土壤中分离得到的一株益生菌,产生的碱性蛋白酶在低温下有作为洗涤剂添加酶的潜力。深入分析菌株NWMCC0046的基因组序列信息,并挖掘该菌株功能特性基因及潜在应用价值。使用PacBio RS II平台和Illumina HiSeq 4000平台对菌株NWMCC0046的基因组进行测序,并对测序数据进行基因组组装、基因预测与功能注释、共线性分析、进化分析及次级代谢产物合成基因簇预测。菌株NWMCC0046全基因组大小为4 321 565 bp,平均GC含量为46.78%,共编码4 504个基因。基因注释揭示了其益生菌特性,如胃肠道内独特的适应性、抗氧化活性和抗菌活性。此外,菌株NWMCC0046还编码工业上许多重要的酶。进化树及共线性结果表明,菌株NWMCC0046属地衣芽胞杆菌且与地衣芽胞杆菌ATCC 14580具有较好的共线性。同时,预测到菌株NWMCC0046中有11个次级代谢产物合成基因簇,编码地衣素、丰原素、杆菌肽和丁酰苷菌素等生物活性物质。基因组信息存储于GenBa...     

生防菌解淀粉芽孢杆菌 SQ-2 全基因组测序及生物信息学分析

【背景】解淀粉芽孢杆菌SQ-2是从市售豆瓣酱中分离得到的一株益生菌,实验表明菌株SQ-2具有较强的抑制植物病原真菌的能力,说明SQ-2具有较好的生物防治能力,有作为生物农药的潜力。【目的】研究菌株SQ-2的遗传信息并揭示其抑菌机制。【方法】在Illumina MiSeq X10平台上对菌株SQ-2进行全基因组测序,使用Trim Galore V.0.4.0清理原始数据,并使用FastQC检查质量;使用SOAPdenovo2进行从头组装,使用antiSMASH鉴定负责次级代谢产物生物合成的基因。【结果】解淀粉芽孢杆菌SQ-2基因组大小为3 486 537 bp,GC含量为46.63%,共编码4298个基因,共发现11个与次级代谢产物生物合成相关的基因簇,其中6个被确定为抗真菌物质合成簇,分别编码bacillaene、bacilysin、butirosin、fengycin、bacillibactin和surfactin。基因组框架测序数据提交至NCBI获得GenBank登录号为JAHXSB000000000。超高液相色谱/质谱联用分析表明产生了儿茶酚型嗜铁素、多烯类和表面活性肽。【结论】...     

鸭源鸡杆菌的分离鉴定及其全基因组序列分析

【背景】鸭源鸡杆菌作为一种条件致病菌能引起家禽卵巢炎、输卵管炎和腹膜炎等疾病,严重威胁养殖业的发展。【目的】四川某养鸡场送检了一批疑似感染鸭源鸡杆菌的病死鸡,为探究其感染机制与防治方法,对该菌进行分离鉴定及全基因组测序分析。【方法】从病料中分离并纯化细菌,再依次进行生化试验、16S rRNA基因序列分析和药敏试验,同时通过全基因组测序对其进行物种分型与毒力、耐药等基因功能注释及遗传进化分析。【结果】该分离菌被鉴定为鸭源鸡杆菌,菌株命名为TS0001,药敏试验显示其仅对硝基呋喃类和少数β-内酰胺类药物敏感,对部分β-内酰胺类、氯霉素、部分氨基糖苷类、大环内酯类、四环素类和磺胺类药物具有耐药性。全基因组序列长度为2 626 722 bp,蛋白质编码基因功能注释显示其有较强的自我加工修饰能力,全基因组注释到83个与毒力因子和耐药性相关的基因,包含4个前噬菌体区域。序列类型分析结果显示,该菌株为ST69型,而且管家基因联合建树表明其与墨西哥普通家鸡分离株7990一致性最高。【结论】本研究为鸭源鸡杆菌的感染机制与防治研究提供了参考,丰富了后续研究的分子生物学背景。     

解淀粉芽孢杆菌生防菌BS-3全基因组测序及生物信息分析

[背景]解淀粉芽孢杆菌BS-3是从健康橡胶树树根中分离获得的一株对真菌具有较强抗菌活性的内生细菌,有作为生物农药的潜力.[目的]解析菌株BS-3的基因组序列信息,以深入研究该菌株防病促生机制及挖掘次级代谢产物基因资源.[方法]采用第二代BGISEQ与第三代PacBio平台相结合的测序技术,对生防菌BS-3进行全基因组测...     

微生物基因组全序列分析进展

目前已测定了１３种微生物的基因组全序列，为追踪生物进化网络提供了宝贵信息，本文综述了已经取得的重要成果，并提出进一步研究的意见。     

一株高度变异的中国SV40分离株的全基因组序列分析

对SV40中国云南分离株YNQD38进行了全基因组核苷酸序列测定。覆盖了整个基因组的9个重叠的基因片段被扩增和测序,与其它SV40株进行了序列比对并基于全基因序列建立了遗传进化树。结果显示:基因组全长5125bp,基因组构成与其它SV40毒株相似,均有6个开放读码框架和1个调控区。YNQD38与已被证实高度保守的其它SV40比,全基因组核苷酸同源性仅为91.0%。在SV40的保守区VP1、VP2、VP3、小t抗原(t-ag)和部分大T抗原(不包括大T抗原C末端)区,YNQD38与其它SV40之间核苷酸同源性分别为90.7%~91.1%、91.7%~92.0%、90.2%~90.8%、92.8%~93.3%、88.5%~89.7%。在SV40的可变区大T抗原C末端(T-ag-C)编码区,YNQD38同源性更低,仅为65.7%~74.3%。YNQD38发生在保守区的核苷酸变异多为无义突变,而发生在变异区的核苷酸变异多为有义突变。YNQD38的调控区缺少一个完整的72bp增强子,这种特别的调控区的结构以前未见报道。基于整个基因组构建的进化树显示该株病毒形成了一个独特的组。以上结果表明YNQD38是目前报道的SV40中变异最大的一株,而且也是第一株被完整测序的SV40中国株。这个报道不仅为SV40中国株的基础研究提供了一个完整清楚的分子生物学资料,还对这样一株高度变异的SV40能否成为人类致病因子进行了初步探讨。     

表型特征结合基因组分析鉴定不同菌落形态Bacillus velezensis ACCC 19742

在中国农业微生物菌种中心(ACCC)菌株的定期转接保藏过程中,发现解淀粉芽孢杆菌ACCC 19742在同一培养基上出现两种不同的菌落形态,将这两个不同形态的菌株编号为19742-1和19742-2。通过形态学、生理生化及基因组分析相结合鉴定不同菌落形态ACCC 19742,并进一步确定该菌株的分类地位。首先将菌株进行分离与纯化,其次将纯化后的菌株进行16S rRNA及gyrB基因扩增及序列分析,通过MEGA 7.0软件构建系统发育树;API 20NE、BIOLOG及脂肪酸等分析菌株的生理生化特性;全基因组分析菌株的ANI和DDH值。两株菌在API 20NE中,仅葡萄糖酸盐同化反应存在差异;脂肪酸检测中主要组成相同,仅是百分含量方面略有差别;两株菌的16S rRNA基因相似性为100%,gyrB基因的相似性为99.4%;全基因组测序表明,两株菌的ANI值为99.95%,DDH值为99.62%。综合遗传学特征和表型特征,证实两者为来源于同一菌株不同的形态变异型,而并非污染所致。同时,19742-1和19742-2与Bacillus velezensis NRRL_B 41580~T的ANI及DDH值最高,分别为97%和77%,且16S rRNA和gyrB系统进化分析也表明,该菌株在分类地位上属于贝莱斯芽孢杆菌(Bacillus velezensis),而非解淀粉芽孢杆菌。这为菌株的保藏提供了一定的参考价值。     

来自南极洲类诺卡氏菌(Nocardioides sp.) InS609-2的全基因组序列分析

【背景】类诺卡氏菌（Nocardioides sp.） InS609-2是一株分离自南极洲罗斯海特拉诺瓦湾的恩克斯堡岛土壤中潜在的极地放线菌新种。尚无研究报道Nocardioides sp.InS609-2的全基因组序列，缺少对其功能基因、代谢产物合成途径及比较基因组学等的研究。【目的】解析Nocardioides sp.InS609-2的基因组序列信息，以深入挖掘次级代谢产物基因资源。【方法】利用Illumina HiSeq高通量测序平台对菌株InS609-2进行全基因组完成图测序，使用相关软件对测序数据进行基因组组装、基因预测和功能注释、预测次级代谢产物合成基因簇和共线性分析等。【结果】基因组最后得到的总长度为4 524 052 bp，G+C含量为69.42%，预测到4 656个基因、56个tRNA和6个rRNA。根据Nocardioides sp.InS609-2的全基因组测序结果，分别有3 761、3 052、1 767、4 134和2 725个基因在COG、GO、KEGG、NR和Swiss-Prot数据库中提取到注释信息。同时，还预测得到19个次级代谢产物合成基因簇。基因组测序数据提交至NCBI获得GenBank登录号为CP060034。Nocardioides sp.InS609-2与N. dokdonensis CP015079、N. yefusunii CP034929、N. euryhalodurans CP038267、N. seonyuensis CP038436、N. daphniae CP038462和N. okcheonensis CP087710这6株基因组同源性比较高的类诺卡氏菌进行共线性分析和蛋白聚类分类分析，得到的结果是7个基因组间既有保守性又各自有独特性。七个基因组共有44个蛋白聚类簇。最后进行16S rRNA基因系统发育树、泛基因组、core基因组和基因组进化树分析，进一步挖掘了Nocardioides sp.InS609-2的基因组信息。【结论】从基因组层面上预测了Nocardioides sp.InS609-2的次级代谢产物的生物合成基因簇，为InS609-2的后续相关研究提供了参考信息，具有重要意义。     

一株猪葡萄球菌的分离鉴定、生物学特性研究及全基因组测序分析

【背景】2021年6月,广东省茂名市某散养户送检了一头发病仔猪,猪身上长有脓疱,四肢关节肿大,关节内可见脓液。【目的】确定引起仔猪发病的病原菌,分析其药物敏感性,为临床用药提供指导;对分离菌株进行全基因组序列分析,挖掘其毒力因子和耐药基因,揭示该菌致病和耐药的分子机制。【方法】取关节脓液分离细菌;通过革兰氏染色、16S rRNA基因和全基因组测序分析,鉴定细菌种类;通过溶血试验、血浆凝固酶试验和生长曲线测定,确定分离菌株的溶血活性、血浆凝固酶活性和生长特性;用小鼠感染模型评估分离菌株的致病性;用纸片扩散法测定分离菌株的药物敏感性;通过全基因组序列分析挖掘分离菌株的毒力因子和耐药基因。【结果】分离菌株被鉴定为猪葡萄球菌(Staphylococcus hyicus);该菌不溶血,无血浆凝固酶活性,在胰蛋白胨大豆肉汤培养基中于37℃、120r/min条件下生长良好;小鼠感染试验结果显示,该菌具有高致病性;药敏试验结果显示,该菌对苯唑西林、大观霉素等7种药物敏感,对青霉素G、红霉素等9种药物耐药;全基因组序列分析结果显示,该菌携带多个毒力因子和耐药基因。【结论】从发病仔猪的关节脓液中分离到一株猪葡萄球菌,可用苯唑西林、大观霉素等药物防控该菌感染;解析了该菌的基因组信息,为后续深入研究该菌致病和耐药的分子机制奠定了基础。     

一株花椒根腐病拮抗菌的筛选、发酵条件及其抑菌物质的初步分析

【背景】腐皮镰孢(Fusarium solani)是陇南花椒根腐病的主要病原物,是危害花椒生产的重要原因之一。【目的】筛选拮抗腐皮镰孢的生防菌株,为开发防治花椒根腐病的微生物农药提供新的生物资源。【方法】采用平板对峙法结合16S rRNA基因序列分析等方法从花椒根围中分离具有较强抑菌效果的细菌并进行筛选鉴定;通过Box-Benhnken中心组合设计试验探究该菌最佳发酵条件;采用液质联用(LC-MS)方法对该菌发酵液的有机溶剂提取相进行非靶向代谢组学分析及KEGG、Lipidmaps数据库注释。【结果】分离筛选出一株拮抗菌株贝莱斯芽胞杆菌(Bacillus velezensis) W-1,该菌株对腐皮镰孢的抑菌率为89%,对藤仓赤霉(Gibberella fujikuroi)、禾谷镰孢(Fusarium graminearum)、梅毒镰孢(Fusarium syphilis)、层生镰孢(Fusarium proliferatum)和松针刺盘孢(Colletotrichum fioriniae)的抑制率为57%-90%。优化后W-1的最佳发酵条件为pH值6.5、温度30℃、接种量5%、转速...     

一株兔源A型多杀性巴氏杆菌的分离鉴定和全基因组测序及分析

【背景】多杀性巴氏杆菌可导致猪肺疫、牛出血性败血症和兔出血性败血症等多种疾病,严重威胁多种动物畜牧养殖业的健康发展。【目的】重庆某兔场送检一批病死兔,为研究其病原和治疗方法,对病原进行了微生物分离和全基因组测序分析。【方法】从2022年重庆某兔场送检兔病料中进行细菌分离纯化、生化试验、16S rRNA基因鉴定、荚膜血清型分型、药敏试验和毒力基因检测,同时通过全基因组测序结果进行毒力、耐药基因注释和遗传进化等分子生物学信息分析。【结果】该菌为兔源A:ST74多杀性巴氏杆菌,命名为LXSS001,基因组序列上传到NCBI数据库(登录号为CP119523.1),药敏试验显示该菌对四环素、杆菌肽、复方新诺明和磺胺异恶唑耐药,对头孢噻肟、头孢哌酮和丁胺卡那等药物敏感。全基因组长度为2 480 671 bp,并注释到了58个毒力基因和9类药物的靶向抗药基因。通过联合建树表明其与3480株一致性最高。【结论】本研究完成了一株A型多杀性巴氏杆菌的分离鉴定和全基因组测序,并揭示了其与国内外其他分离株的进化关系,为多杀性巴氏杆菌的后续研究提供了参考依据。     

一株鳢源鰤诺卡氏菌致病性与全基因组分析

【背景】鰤诺卡氏菌是一种典型的条件致病菌,感染鳢、鲈等多种名优鱼类,易造成存在机体损伤或免疫机能下降的鱼持续性感染,给水产养殖业造成了巨大损失。【目的】了解临床分离鳢源鰤诺卡氏菌对乌斑杂交鳢（斑鳢♀×乌鳢♂）的致病性,并从全基因组层面了解该病原菌的基因组和致病因子信息,为鰤诺卡氏菌后续病原学及鰤诺卡氏菌病防治技术和疫苗的开发研究提供有利的数据支撑。【方法】鰤诺卡氏菌NK201610020通过回归感染试验和发病鱼靶器官组织病理分析,了解鰤诺卡氏菌的毒力和病理特征。通过对试验菌进行全基因组测序和比较基因组学分析,挖掘该菌的基因组与毒力特征。【结果】回归感染试验结果显示,除1.5×10³组外,其余5个感染组致死率高达90%,LD₅₀为1.079×10³ CFU/mL,说明试验菌毒力较强。组织病理学观察到肝、脾、肾呈现严重的病理损伤,而且有肉芽肿结构形成。试验菌全基因组测序发现,全基因组大小为8294329bp,GC含量为68.10%,共预测到编码基因7812个。比较基因组学分析发现,不同地区不同宿主来源鰤诺卡氏菌在基因组基础特...     

植物内生链霉菌Streptomyces sp. SAT1的基因组测序和比较基因组分析

【背景】一直以来,链霉菌都是活性物质的主要生产者,近年来随着抗生素滥用引起的环境和微生物抗药性问题越发严重,挖掘高效生物防治因子和新型抗生素成为了解决以上问题的重要手段。【目的】通过获得植物内生链霉菌SAT1全基因组序列和次级代谢基因簇信息,利用比较基因组学和泛基因组学分析SAT1菌株的特殊性以及与其他链霉菌的共性,为阐明SAT1抑菌和内生机制提供理论基础,为揭示链霉菌的生态功能提供可靠数据。【方法】通过三代测序平台PacBio Sequel完成SAT1基因组测序,利用生物信息学技术进行注释和功能基因分类;分别利用RAxML和PGAP软件进行系统发育树的构建和泛基因组分析;次级代谢基因簇的预测和分析通过antiSMASH网站完成。【结果】获得SAT1菌株的全基因组完成图,该菌线性染色体长度约7.47 Mb,包含有4个质粒,GC含量近73%,共预测到7 550个蛋白编码基因,含有37个次级代谢基因簇,分属29个类型,其中默诺霉素基因簇与加纳链霉菌具有较高相似性。42株代表链霉菌中,单个菌株次级代谢基因簇数量为20-55个,主要类型为PKS类、Terpene类和Nrps类,而且含有大量杂合基因簇,各个菌株中特有基因数目较为庞大。【结论】链霉菌SAT1菌株在基因组特点以及次级代谢基因簇的数量和类型上与其余41株链霉菌具有一定的共性,其中潮霉素B基因簇和默诺霉素基因簇合成的相关物质可能与SAT1抑菌活性密切相关。42株链霉菌中次级代谢基因簇数量的多少与基因组大小成正相关,同时大量杂合基因簇以及庞大的特有基因数目的存在说明链霉菌在长期进化过程中存在了很高程度的水平基因转移现象,可能具有重要的生态功能。     

出芽短梗霉菌株PA-2全基因组测序及分析

【背景】出芽短梗霉菌株PA-2是一株分离自青海省海东市平安区患病杨树叶片上的真菌,前期研究表明该菌株具有除草和抑菌能力,说明其在生物农药方面具有潜在的应用前景。【目的】了解菌株PA-2的基因组序列信息,挖掘其生防相关功能基因簇,为进一步研究解析该菌株生防机理及生防功能改造提供遗传背景信息。【方法】利用IlluminaHiSeq高通量测序平台对生防菌株PA-2进行全基因组测序,用生物信息学的方法对测序数据进行基因组组装、基因预测及功能注释、碳水化合物活性酶预测、次级代谢产物合成基因簇预测,利用刚果红染色等方法对水解酶活性进行衡量。【结果】菌株PA-2基因组序列全长28 932 793 bp,平均GC含量为50%,共编码10 839个基因,预测到该菌株具有4个已知的次级代谢产物合成基因簇,编码Melanin、Burnettramic Acid A、ACR-Toxin I、Abscisic Acid,该菌株能水解纤维素和果胶。【结论】有助于在基因组层面上解析菌株PA-2生防机制的内在原因,为深入了解出芽短梗霉菌次级代谢物合成途径提供参考,对菌株PA-2的下一步相关研究具有重要意义。     

一株嗜酸硫杆菌M4-422-6的分离及其全基因组测序和比较基因组学分析

【目的】开展具有硫氧化能力的嗜酸硫杆菌属(Acidithiobacillus)的分离及其比较基因组学分析,不仅可以丰富硫氧化细菌菌种资源,而且有助于加深理解嗜酸硫杆菌的分子进化与生态适应机制。【方法】利用以硫代硫酸钠为唯一能源的培养基分离具有硫氧化能力的细菌;利用Illumina HiSeq X和Oxford Nanopore测序平台对一株嗜酸硫杆菌M4-422-6进行全基因组测序;利用相关生物信息学分析软件对原始数据进行组装和基因组注释,并与一株亲缘关系最近的菌株Igneacidithiobacillus copahuensis VAN18-1进行比较基因组学分析。【结果】分离获得一株具有硫氧化能力的嗜酸硫杆菌M4-422-6。基因组注释结果显示,菌株M4-422-6基因组由1个染色体和2个质粒组成,基因组大小为2 917 823 bp,G+C含量为58.54%,共编码2 925个蛋白。16S rRNA基因和基因组系统发育树显示,菌株M4-422-6代表嗜酸硫杆菌属的一个潜在新种。功能基因注释结果显示,菌株Acidithiobacillus sp. M4-422-6拥有与菌株特性相关的众多基因,包括硫氧化相关基因、CO₂固定相关基因和耐酸相关基因。比较基因组学分析发现,虽然菌株M4-422-6与VAN18-1的亲缘关系最近,但两者仍拥有众多的差异基因,主要包括噬菌体抗性相关基因和移动元件编码基因。【结论】菌株M4-422-6代表嗜酸硫杆菌属的一个潜在新种,该菌株具有同种内菌株所不具有的特有基因,并据此推测嗜酸硫杆菌种内分化可归因于对特定生态位的适应。     

Comparative analysis of whole genome structure of Streptococcus suis using whole genome PCR scanning

An outbreak associated with Streptococcus suis infection in humans emerged in Sichuan province, China in 2005. The outbreak is atypical for the apparent large number of human cases, high fatality rate and geographical spread. To determine whether the bacterium has changed, we compared both human and animal isolates from the Sichuan outbreak with those collected previously within China and in other countries using whole genome PCR scanning (WGPScaning) comparative sequencing of several known virulence factor genes and multilocus sequence typing (MLST) analysis. WGPScanning analysis showed that all primer pairs yielded PCR products of the expected sizes in all four strains tested. The nucleotide sequences of all the detected virulence factor genes are identical in the four strains and MLST results showed that the four isolates studied and reference strain all belonged to the ST1 complex. No new genetic changes were found in the genome structure of the isolates from this Sichuan outbreak. Contributed equally to this work Supported by the National Key Technologies Research and Development Program (Grant No. 2005BA711A09) from the Ministry of Science and Technology of China