贺 羽,赵秘密,杨致邦,周 汛,3
孢子丝菌复合体的分子生物学研究进展
贺羽1,赵秘密1,杨致邦2,周汛1,3
孢子丝菌是人兽共患性真菌,感染人或动物可引起孢子丝菌病。以往形态学鉴定认为申克孢子丝菌(Sporothrixschenckii,S.schenckii)为孢子丝菌病的唯一致病菌种,而近年来分子生物学研究发现该菌是6个菌种构成的复合体,并且在基因鉴定、基因组学、发病机制等方面有了进一步研究,有利于理解孢子丝菌病的发生发展过程,对于早期诊断、确定治疗方案及判断预后至关重要。本文就孢子丝菌复合体的分子生物学研究进展进行综述。
孢子丝菌;基因鉴定;基因组;发病机制;分子生物学
孢子丝菌是人兽共患性真菌,感染人或动物可引起皮肤、皮下组织及附近淋巴系统的慢性肉芽肿性病变,偶可经血液、淋巴系统播散引起系统性损害[1],严重时可危及生命。孢子丝菌由美国人Shenk首次发现,故被命名为申克孢子丝菌[2]。以往形态学鉴定认为申克孢子丝菌为孢子丝菌病的唯一致病菌种,而近年来分子生物学研究发现该菌是多个菌种构成的复合体[3]。生物信息学的发展让我们对孢子丝菌有了更多更新的认识,本文就近年来孢子丝菌复合体在基因鉴定、基因组学、发病机制方面的分子生物学研究进展进行综述。
1.1基因分型孢子丝菌的基因分型通常基于mtDNA限制性内切酶HaeⅢ的限制性片段长度多态性(RFLP)分析,Lin等首次将不同地区S.schenckii菌株的mtDNA分成24种亚型,而后25到30种,以及31到32种亚型被相继发现[4],并将32种亚型划分为A、B 2个系统发育组,不同亚型其地理分布存在差异。研究者采用ITS-RELP分析将rDNA分为4个亚型[5],其中Ⅰ-Ⅲ型与A组一致,Ⅳ型与 B组一致。张振颖等采用RFLP 技术和Southern印迹法进行rDNA基因分型发现我国不同地区31株临床孢子丝菌显示出15种不同的DNA带型(A-O),51.6%DNA模带为A-C型,南方和北方菌株存在明显差异,但不能明确基因型与临床型别间的关系[6]。张萍等采用多聚酶链反应随机扩增多态性DNA法研究发现新疆申克孢子丝菌显示出新基因型,与北京固定型孢子丝菌病致病菌存在种内基因差异,而与北京、哈尔滨淋巴型孢子丝菌致病菌间未见种内基因差异[7]。
1.2菌种鉴定及基因诊断通过保守序列包括CAL、ITS及β-tubulin基因变异分析发现孢子丝菌复合体包含6个菌种,包括巴西孢子丝菌(Sporothrixbrasiliensis,S.brasiliensis)、球形孢子丝菌(Sporothrixglobosa,S.globosa)、申克孢子丝菌等[3, 8-9],菌种不同其毒力、地理分布、药物敏感性等存在差异。分子生物学技术的发展为准确而快速鉴定菌种提供了新方法,对于早期明确致病菌种类型以指导治疗具有重要意义。利用通用引物T38的PCR指纹图谱分析发现各隐含种显示出不同电泳条带,与利用CAL基因分析鉴定菌种的结果一致[10]。Rodrigues等运用CAL基因的PCR-RFLP分析实现了对孢子丝菌复合体各隐含种的基因鉴定[11],认为该方法准确、简单且成本低,可用于常规的菌种分类。Oliveira等首次将基质辅助激光解吸电离飞行时间质谱术(MALDI-TOF-MS)用于鉴定孢子丝菌各隐含种,并指出通过蛋白质谱分析可准确、快速鉴定6个菌种[12]。近期滚环扩增技术(RCA)被用于孢子丝菌不同菌种间的鉴定,研究表明其特异性和敏感性均较高[13]。孙田等利用聚合酶链反应-酶联免疫法(PCR-ELISA)筛选与孢子丝菌结合效率相对较高的探针,结果显示探针U26852、U26866、U26866、M85053及AF117945均位于PCR扩增产物序列上,而U26852显色最强,EI 值最高,由此推断探针U26852的结合效率相对较高[14]。
2.1结构基因组学Tateishi等对8株日本临床S.schenckii进行染色体核型分析显示S.schenckii基因组全长约28 Mbp,包含6~8条染色体,每条长度为460到6 200 kb,而Torres-Guerrero对不同地理来源菌株进行研究发现其基因组全长约45 Mbp[15],这种差异可能由实验菌株为不同隐含种或研究方法不同引起。近期Teixeira等通过比较基因组学研究发现S.schenckii和S.brasiliensis在基因组长度、转座子、线粒体基因等方面存在差异[16]。前者基因组含17条染色体,全长32.3 Mb,后者含20条染色体,全长33.2 Mb;转座子分别由0.34%和0.62%基因组构成;线粒体基因组比较显示编码归巢核酸内切酶的内含子只存在于巴西孢子丝菌。而双向基因组对比显示两者染色体序列相似性高达97.5%;基因组进化分析发现两个隐含种均大约在3.8~4.9万年前开始分化,推断为最近的一次物种形成事件;两者基因组中与植物腐烂相关的多糖裂解酶基因缺乏表明孢子丝菌从植物致病性向动物致病性的进化演变。
2.2功能基因组学袁立燕等利用BLASTx等生物学软件对S.schenckii的cDNA中编号为Locus-168-Contig-1的序列进行分析发现该基因与亚精胺合成酶基因同源,全长1 062 bp,可编码291个氨基酸;编码蛋白疏水性高,为非分泌型蛋白,无质体、线粒体等定位序列及跨膜螺旋结构;并成功克隆全长编码区基因和构建重组表达质粒[17]。采用同样方法翻译控制肿瘤蛋白(TCTP)的cDNA全长序列被获得,其结构及功能、潜在抗原表位被预测出,为进一步研究该蛋白生物学特性奠定了基础[18]。通过RNA干扰技术将孢子丝菌钙/钙调蛋白激酶1基因表达抑制后菌株在35 ℃环境下生长抑制表明该基因与菌株耐热性相关;通过基因敲入,可表达绿色荧光蛋白(GFP)的孢子丝菌已成功构建,为研究孢子丝菌蛋白质组学和基因组学提供了新技术[19]。Rodriguez-Caban等通过RAN干扰和酵母双杂交技术明确了HSP90与钙/钙调蛋白激酶1在孢子丝菌双相转换中的相互作用[20]。王晓慧等通过简并PCR结合RACE技术成功克隆S.schenckii的未知过氧化氢酶基因并将其命名为Sscat,为深入研究该基因结构及其功能、蛋白质表达奠定了基础[21]。
孢子丝菌病的发病机制与病原体毒力和宿主免疫状态密切相关。毒力和宿主免疫状态的分子生物学研究能够揭示疾病的发生发展过程,为孢子丝菌病治疗提供新靶点。
3.1病原体毒力孢子丝菌毒力主要与色素、黏附力及耐过氧化物能力等相关。孢子丝菌菌丝相和酵母相均可生成色素,包括1,8二羟基萘黑色素、真黑色素、脓褐素。研究发现分生孢子的黑化作用能够抵抗吞噬细胞的吞噬[22],此外Madrid等分别用黑化和白化孢子丝菌感染小鼠发现前者仅在接种部位形成肉芽肿样皮损,而后者形成多病灶性皮损,因此推断色素可能与病原体播散相关[23]。孢子丝菌粘附素通过与纤连蛋白、层粘连蛋白和Ⅱ型胶原蛋白结合粘附宿主,两个生长时相均可表达粘附素,但酵母相粘附细胞外基质的能力更强[24]。Ruiz-Baca等首次通过二维凝胶电泳和免疫印迹法对S.schenckii的免疫原性蛋白进行检测发现70 kDa糖蛋白(Gp70)是2个生长时相均存在的主要细胞壁抗原[25],而Gp60主要存在于酵母相[26]。研究显示Gp70是重要的粘附素,参与病原体入侵宿主的过程。转化生长因子β1可增加内皮下细胞外基质的暴露诱导孢子丝菌黏附内皮单层而经内皮迁移,在免疫功能受损患者中易致播散型孢子丝菌病[27]。近期研究发现S.brasiliensis的耐过氧化物能力较S.schenckii更强[28],计算机模拟分析推断这种差异可能与氧化应激信号通路中类似AP-1的转录因子差异表达和Hog1基因突变有关。致病性相关蛋白的研究目前甚少,机体感染孢子丝菌6周后,孢子丝菌细胞壁的糖肽复合物能够抑制宿主免疫反应,被认为与真菌致病性相关[19];从S.schenckii酵母相提取抗原过程中可水解人IgG的分泌蛋白被发现,该蛋白可能参与孢子丝菌抗宿主的免疫反应[29]。3.2宿主免疫状态病原体入侵后机体迅速作出防御反应,包括固有免疫和特异性免疫,后者分为细胞免疫和体液免疫。补体系统可被孢子丝菌激活,特别是替代途径激活后形成的膜攻击复合物可溶解破坏真菌细胞壁[16]。研究发现酵母相细胞壁模式识别受体Toll样受体4(TLR4)即CD284可促进促炎性细胞因子释放和氧化物质如NO的生成[30];Negrini等通过构建感染孢子丝菌的TLR2缺陷小鼠进行研究发现巨噬细胞吞噬功能受损,并且促炎性细胞因子分泌水平明显降低[31],可见TLRs在孢子丝菌固有免疫中扮演着重要角色[30-32]。参与细胞免疫的细胞包括巨噬细胞、树突状细胞、肥大细胞等。Uenotsuchi等运用流式细胞技术分析发现分离自皮肤的S.schenckii可快速激活单核细胞来源的树突状细胞(MoDCs)引起Th1细胞分泌大量IFN-γ参与反应,而分离自内脏的S.schenckii仅激活少量MoDCs和Th1细胞,有丝分裂原激活蛋白激酶P38和c-Jun氨基末端激酶信号通路被认为参与该反应过程[33]。研究发现孢子丝菌分生孢子可激活肥大细胞释放大量TNF-α和IL-6,而复合物48/80可诱导加强分生孢子刺激肥大细胞释放组胺,此外该细胞的功能消耗能够显著降低皮损程度[34]。Kajiwara等用孢子丝菌感染NADPH氧化酶和抗氧化酶缺陷的肉芽肿样疾病(chronic granulomatous disease,CGD)小鼠进行研究发现CGD小鼠不能清除接种部位真菌导致系统播散而死亡,由此可见NADPH氧化酶缺乏可致孢子丝菌致命性的感染,同时免疫细胞产生的超氧阴离子及其代谢产物如活性氧、活性氮在抑制和杀灭真菌中起着重要作用[35-36]。感染孢子丝菌5到6周后,Th2细胞产生IL-4启动机体体液免疫,体外实验中抗神经鞘糖脂类抗原的抗体可阻碍孢子丝菌的生长及分化[37]。Nascimento等研究发现感染孢子丝菌的小鼠可产生抗Gp70的特异性IgG3和IgG1抗体,可能与机体清除病原体有关[38],而近期研究认为被动免疫抗Gp70的抗体可保护小鼠免遭孢子丝菌感染[39]。
近年来随着分子生物学的发展,对孢子丝菌有了更深入认识。生物学技术的运用在基因层面上实现了对孢子丝菌的菌种鉴定,但因技术要求较高,仅限于实验室研究,探索临床快速鉴定菌种和早期诊断疾病的方法仍是目前研究热点。孢子丝菌基因组DNA由30万左右个碱基组成,确定各部位结构及功能已经成为未来研究方向和目标。同时明确发病机制,寻找新靶点在治疗及预后判断方面有重要指导意义。
[1] Tellez MD, Batista-Duharte A, Portuondo D, et al.Sporothrixschenckiicomplex biology: environment and fungal pathogenicity[J]. Microbiology, 2014, 160(11): 2352-2365. DOI: 10.1099/mic.0.081794-0
[2] Liu TT, Zhang K, Zhou X. Molecular identification ofSporothrixclinicalisolates in China[J]. J Zhejiang Univ Sci B, 2014, 15(1): 100-108. DOI: 10.1631/jzus.B1300136
[3] Rodrigues AM, de Hoog GS, de Camargo ZP. Molecular diagnosis of pathogenicSporothrixspecies[J]. PLoS Negl Trop Dis, 2015, 9(12): 1-22. DOI: 10.1371/journal.pntd.000419
[4] Ishizaki H, Kawasaki M, Anzawa K, et al. Mitochondrial DNA analysis ofSporothrixschenckiiin India, Thailand, Brazil, Colombia, Guatemala and Mexico[J]. Nihon Ishinkin Gakkai Zasshi, 2009, 50(1): 19-26.
[5] Watanabe S, Kawasaki M, Mochizuki T, et al. RFLP analysis of the internal transcribed spacer regions ofSporothrixschenckii[J]. Nihon Ishinkin Gakkai Zasshi, 2004, 45(3): 165-175.
[6] Zhang ZY, Liu XM, Yang GL, et al. Genotyping ofSporothrixschenckiiby analysis of ribosomal DNA regions[J]. Mycoses, 2006, 49(4): 305-310.
[7] Zhang P, Palida ABLZ, Chen Y, et al. New genotype ofSporothrixschenckiiin Xinjiang[J]. J Clin Dermatol, 2011, 40(5):262-264. (in Chinese)
张萍,帕丽达.阿布利孜,陈宴,等.发现于新疆的申克孢子丝菌新基因型[J].临床皮肤科杂志, 2011, 40(5): 262-264.
[8] Rodrigues AM, de Hoog S, de Camargo ZP. Emergence of pathogenicity in theSporothrixschenckiicomplex[J]. Med Mycol, 2013, 51(4): 405-412. DOI: 10.3109/13693786.2012.719648
[9] Ottonelli SCD, Magagnin CM, Castrillon MR, et al. Antifungal susceptibilities and identification of species of theSporothrixschenckiicomplex isolated in Brazil[J]. Med Mycol, 2014, 52(1): 56-64. DOI: 10.3109/13693786.2013.818726
[10] de Oliveira MM, Sampaio P, Almeida-Paes R, et al. Rapid identification ofSporothrixspecies by T3B fingerprinting[J]. J Clin Microbiol, 2012, 50(6): 2159-2162. DOI: 10.1093/femsle/fnv027
[11] Rodrigues AM, de Hoog GS, de Camargo ZP. Genotyping species of theSporothrixschenckiicomplex by PCR-RFLP of calmodulin[J]. Diagn Microbiol Infect Dis, 2014, 78(4): 383-387. DOI: 10.1016/j.diagmicrobio.2014.01.004
[12] Oliveira MM, Santos C, Sampaio P, et al. Development and optimization of a new MALDI-TOF protocol for identification of theSporothrixspecies complex[J]. Res Microbiol, 2015, 166(2): 102-110. DOI: 10.1128/JCM.00450-12
[13] Rodrigues AM, Najafzadeh MJ, de Hoog GS, et al. Rapid identification of emerging human-pathogenicSporothrixspecies with rolling circle amplification[J]. Front Microbiol, 2015, 6(1385): 1-16. DOI: 10.3389/fmicb.2015.01385
[14] Sun T, Liu XM, Zhang ZY. Screening of species-specific probes forSporothrixschenckiiby PCR-ELISA[J]. J Pract Dermatol, 2014, 7(6): 403-405. (in Chinese)
孙田,刘晓明,张振颖.聚合酶链反应-酶联免疫法筛选孢子丝菌特异性探针[J].实用皮肤病学杂志, 2014, 7(6): 403-405.
[15] Barros MB, de Almeida Paes R, Schubach AO.SporothrixschenckiiandSporotrichosis[J]. Clin Microbiol Rev, 2011, 24(4): 633-654. DOI: 10.1128/CMR.00007-11
[16] Teixeira MM, de Almeida LG, Kubitschek-Barreira P, et al. Comparative genomics of the major fungal agents of human and animal Sporotrichosis: Sporothrix schenckii andSporothrixbrasiliensis[J]. BMC Genomics, 2014, 15(1): 943-965. DOI: 10.1186/1471-2164-15-943
[17] Yuan LY, Huang HQ, Zhang J, et al. Bioinformatics analysis and molecular cloning of spermidine synthase gene fromSporothrixschenckii[J]. J Sun Yat-sen Univ (Med Sci), 2013, 34(6): 948-953. (in Chinese)
袁立燕,黄怀球,张静,等.申克孢子丝菌基因的生物信息学分析及分子克隆[J].中山大学学报(医学科学版), 2013, 34(6): 948-953.
[18] Huan HQ, Zhong Y, Zhao J, et al. Molecular structure analysis and function prediction onSporothrixschenckiitranslationally controlled tumor protein[J]. J Sun Yat-sen Univ (Med Sci), 2012, 33(4): 454-459. (in Chinese)
黄怀球,钟毅,赵静,等.申克孢子丝菌翻译控制肿瘤蛋白的结构分析与功能预测[J].中山大学学报(医学科学版), 2012, 33(4): 454-459.
[19] Mora-Montes HM, Dantas Ada S, Trujillo-Esquivel E, et al. Current progress in the biology of theSporothrixschenckiicomplex following the genomic era[J]. FEMS Yeast Res, 2015, 15(6): 1-10. DOI: 10.1093/femsyr/fov065
[20] Rodriguez-Caban J, Gonzalez-Velazquez W, Perez-Sanchez L, et al. Calcium/calmodulin kinase1 and its relation to thermotolerance and HSP90 inSporothrixschenckii: an RNAi and yeast two-hybrid study[J]. BMC Microbiol, 2011, 11(162): 1-16. DOI: 10.1186/1471-2180-11-162
[21] Wang XY, Liu W, Li RY. Cloning a novel catalase gene ofSporothrixschenckiiwith degenerate PCR and RACE[J]. Chin J Mycol, 2011, 6(5): 271-275. (in Chinese)
王晓慧,刘伟,李若瑜.简并结合技术克隆申克孢子丝菌未知过氧化氢酶基因[J].中国真菌学杂志, 2011, 6(5):271-275.
[22] Almeida-Paes R, Frases S, Araujo Gde S, et al. Biosynthesis and functions of melanoid pigment produced by species of theSporothrixcomplex in the presence of L-tyrosine[J]. Appl Environ Microbiol, 2012, 78(24): 8623-8630. DOI: 10.1128/AEM.02414-12
[23] Madrid H, Cano J, Gene J, et al.Sporothrixglobosa, a pathogenic fungus with widespread geographical distribution[J]. Rev Iberoam Micol, 2009, 26(3): 218-222. DOI: 10.1016/j.riam.2009.02.005
[24] Teixeira PA, de Castro RA, Nascimento RC, et al. Cell surface expression of adhesins for fibronectin correlates with virulence inSporothrixschenckii[J]. Microbiology, 2009, 155(Pt11): 3730-3738. DOI: 10.1099/mic.0.029439-0
[25] Lopes-Bezerra LM.Sporothrixschenckiicell wall peptidorhamnomannans[J]. Front Microbiol, 2011, 2(243): 1-4. DOI: 10.3389/fmicb.2011.00243
[26] Ruiz-Baca E, Hernandez-Mendoza G, Cuellar-Cruz M, et al. Detection of 2 immunoreactive antigens in the cell wall ofSporothrixbrasiliensisandSporothrixglobosa[J]. Diagn Microbiol Infect Dis, 2014, 79(3): 328-330. DOI: 10.1016/j.diagmicrobio.2014.02.020
[27] Figueiredo CC, Deccache PM, Lopes-Bezerra LM, et al. TGF-beta1 induces transendothelial migration of the pathogenic fungusSporothrixschenckiiby a paracellular route involving extracellular matrix proteins[J]. Microbiology, 2007, 153(Pt 9): 2910-2921.
[28] Ortega I, Soares FM, Vasconcelos AT, et al. Peroxide sensing and signaling in theSporothrixschenckiicomplex: an in silico analysis to uncover putative mechanisms regulating the Hog1 and AP-1 like signaling pathways[J]. Med Mycol, 2015, 53(1): 51-59. DOI: 10.1093/mmy/myu069
[29] Da RD, Gezuele E, Calegari L, et al. Excretion-secretion products and proteases from liveSporothrixschenckiiyeast phase: immunological detection and cleavage of human IgG[J]. Rev Inst Med Trop Sao Paulo, 2009, 51(1): 1-7.
[30] Sassa MF, Ferreira LS, Ribeiro LC, et al. Immune response againstSporothrixschenckiiin TLR-4-deficient mice[J]. Mycopathologia, 2012, 174(1): 21-30. DOI: 10.1007/s11046-012-9523-1
[31] Negrini TC, Ferreira LS, Alegranci P, et al. Role of TLR-2 and fungal surface antigens on innate immune response againstSporothrixschenckii[J]. Immunol Invest, 2013, 42(1): 36-48. DOI: 10.3109/08820139.2012.719982
[32] Guzman-Beltran S, Perez-Torres A, Coronel-Cruz C, et al. Phagocytic receptors on macrophages distinguish between differentSporothrixschenckiimorphotypes[J]. Microbes Infect, 2012, 14(12): 1093-1101. DOI: 10.1016/j.micinf.2012.06.001
[33] Uenotsuchi T, Takeuchi S, Matsuda T, et al. Differential induction of Th1-prone immunity by human dendritic cells activated withSporothrixschenckiiof cutaneous and visceral origins to determine their different virulence[J]. Int Immunol, 2006, 18(12): 1637-1646.
[34] Romo-Lozano Y, Hernandez-Hernandez F, Salinas E. Mast cell activation by conidia ofSporothrixschenckii: role in the severity of infection[J]. Scand J Immunol, 2012, 76(1): 11-20. DOI: 10.1111/j.1365-3083.2012.02706.x
[35] Grimm MJ, Vethanayagam RR, Almyroudis NG, et al. Monocyte-and macrophage-targeted NADPH oxidase mediates antifungal host defense and regulation of acute inflammation in mice[J]. J Immunol, 2013, 190(8): 4175-4184. DOI: 10.4049/jimmunol.1202800
[36] Kajiwara H, Saito M, Ohga S, et al. Impaired host defense againstSporothrixschenckiiin mice with chronic granulomatous disease[J]. Infect Immun, 2004, 72(9): 5073-5079.
[37] Toledo MS, Tagliari L, Suzuki E, et al. Effect of anti-glycosphingolipid monoclonal antibodies in pathogenic fungal growth and differentiation. Characterization of monoclonal antibody MEST-3 directed to Manpalpha1->3Manpalpha1->2IPC[J]. BMC Microbiol, 2010, 10(47): 1-22. DOI: 10.1186/1471-2180-10-47
[38] Nascimento RC, Almeida SR. Humoral immune response against soluble and fractionate antigens in experimental sporotrichosis[J]. FEMS Immunol Med Microbiol, 2005, 43(2): 241-247.
[39] de Almeida JR, Kaihami GH, Jannuzzi GP, et al. Therapeutic vaccine using a monoclonal antibody against a 70-kDa glycoprotein in mice infected with highly virulentSporothrixschenckiiandSporothrixbrasiliensis[J]. Med Mycol, 2015, 53(1): 42-50. DOI: 10.1093/mmy/myu049
Advances on molecular biology ofSporothrixcomplex
HE Yu1, ZHAO Mi-mi1, YANG Zhi-bang2, ZHOU Xun1,3
(1.DepartmentofDermatology,theFirstAffiliatedHospitalofChongqingMedicalUniversity,Chongqing400016,China;2.ExperimentalTeachingCenterforBasicMedicalPathogenBiologyandImmunologyLaboratoryofChongqingMedicalUniversity,Chongqing400016,China;3.DepartmentofDermatologyandCosmetology,ChongqingHospitalofTraditionalChineseMedicine,Chongqing400021,China)
Sporotrichosis is a polymorphic disease that affects both humans and animals worldwide caused by the fungal genusSporothrix. Sporotrichosis has been attributed to one single etiological agent,Sporothrixschencki, by morphological identification in the past decades. Recently, isolates received asSporothrixschenckiare regrouped with six cryptic species, and further research on molecular biology study on gene identification, genome and pathogenesis is helpful to understand the development and progression of Sporotrichosis, and has important significance on early diagnosis, treatment and prognosis. The present review will focus on recent advances on molecular biology ofSporothrixcomplex.
Sporothrix; gene identification; genome; pathogenesis; molecular biology
Supported by the National Natural Science Foundation of China (No. 31270062)
Zhou Xun, Email: zhouxun123@sina.com
周汛,Email:zhouxun123@sina.com
1.重庆医科大学附属第一医院皮肤科,重庆400016;2.重庆医科大学基础医学院基础医学实验教学中心病原生物学与免疫学实验室,重庆400016;3.重庆市中医院皮肤美容科,重庆400021
R379
A
1002-2694(2016)06-0576-05
2015-12-04;
2016-02-15
DOI:10.3969/j.issn.1002-2694.2016.06.013
国家自然科学基金(No.31270062)资助