仿生控释BMP-2和BMP-7的PELA微球支架的研制及其修复大鼠骨缺损的实验研究

【摘要】临床上治疗粉碎骨折,肿瘤切除,或骨髓炎病灶清除后大段骨缺损仍然是骨科医生面临的一大挑战。自体骨移植一直作为骨缺损治疗的金标准,然而存在着诸如取骨区疼痛,出血,来源有限等不足。同种异体骨,作为自体骨替代品,近些年在临床上应用也较广泛,但其也存在成本较高,病毒传播潜在风险及不能完全避免的免疫排斥反应等问题,限制了其在临床大规模使用。作为自体骨及同种异体骨的替代品,组织工程骨的研究有望彻底解决这一难题。具备骨诱导活性的复合生物材料是组织工程骨的一个研究热点。一些有着较强的诱导成骨活性的生长因子,如重组人骨形态发生蛋白-2(rhBMP-2)及重组人骨形态发生蛋白-2(rhBMP-7),已结合不同的生物材料作为载体,被美国FDA批准用于临床上促进骨折愈合及脊柱融合。但BMP制剂广泛应用于临床后发现由于使用剂量远高于生理水平,活性因子的释放动力学不够科学等原因,出现了诸如局部软组织水肿,脊髓脊神经根炎,异位骨化,种植体周围骨吸收及引发癌症等一系列的并发症。理想的,可以广泛地应用于临床的有成骨诱导能力的人工骨填充材料仍需要进一步研究不断完善。在骨折的自然愈合过程中,有多种细胞因子共同参与,如BMP家族,血管内皮生长因子(VEGF),胰岛素样生长因子(IGF)等,他们的成骨机制包括：诱导间充质干细胞(MSCs)成骨分化,促进局部新生血管的生成,作为细胞激活分化的调控因子等,促进骨组织愈合。BMP家族中BMP-2和BMP-7是被研究得最多,有较强的成骨诱导活性的两种因子。之前的研究表明,BMP-2是骨折初始阶段各种级联反应的启动因子,在骨折后的第一天即大量表达至峰值,而后持续低浓度表达。而BMP-7则在骨折2周后才出现。最近一些研究表明,BMP-2和BMP-7的复合使用比单因子更能促进骨组织再生和骨折愈合。所以,如果研发一种骨修复材料可以仿生控释BMP-2和BMP-7,则可以在安全浓度范围内,以较低的浓度,较少的制备成本最大限度地促进使骨缺损修复。微型包囊技术系利用天然的或合成的高分子材料(囊材)作为囊膜壁壳,将固态药物或液态药物(囊心物)包裹而成药库型微型胶囊。药物在囊内通过扩散和渗透等形式在特定的位置以适当的速度和持续的时间释放出来,以达到更大限度的发挥药效的作用。到目前为止已有200多种药物采用了微囊化技术,尤其是近年来生物大分子作为药物在临床上广泛应用,包括蛋白质,多肽类,基因,生长因子,抗原等都能成功地被包封制成缓释微囊,从而最大限度地保护生物大分子的结构和活性,并能达到靶向给药,控释给药目的,使得药物释放更可控,更安全,更有效,从而带动了制药工业的飞速发展。近年来已有不少学者尝试将BMPs以不同的高分子聚合物包封成微囊制成缓释制剂,或与其他无机材料共同合成能缓释活性因子的人工骨支架,细胞实验及动物实验结果令人满意,如何让缓释系统释放出的BMPs最大限度地保留其诱导成骨活性,又不至于剂量过高引起副作用,是现在BMPs缓释载体研究的新的趋势。一直以来,载药微囊作为控释载体大都作为支架或水凝胶组成成分之一。2007年AnaJaklenec等试验了全新的支架合成方法,将包封小牛血清蛋白(BSA)的PLGA微囊通过二氯甲烷密闭熏蒸法直接融合成支架并获得成功。次年,这个研究团队继续用该方法制备了由包封胰岛素样生长因子-1及转化生长因子β1的微囊融合成的支架用于软骨组织工程的研究,这种支架能成功地缓释这两种生长因子并保留其生物活性。该方法制备过程简单,可通过控制微囊大小及熏蒸时间来控制支架的孔隙率,体外细胞培养可让纤维母细胞顺利附着为其扩增提供支持。这为组织工程骨及软骨支架的制备提供了一条崭新的思路。然而,其进一步的生物学实验成果并未见报道。将包封活性因子的高分子聚合物微囊直接熏蒸融合制成微球支架,它最大的优点在于制备过程简单,一步到位,避免了以往复合支架制备过程中反复添加不同有机溶剂,化学稳定剂及交联剂等,而且不需要旋蒸,冻干,高温,高压等步骤,这样最大限度地避免了制备过程中对功能蛋白质的结构和活性产生的负面影响。同时,合成的支架具有一定的粘弹性,可以根据骨缺损的不同形态予以塑性,并且电子显微镜证实其有一定的孔隙率,适合细胞贴附及营养液的渗入。本研究首先利用复乳溶剂挥发法制备了包封rhBMP-2的微囊,聚乳酸-聚乙二醇-聚乳酸三嵌段共聚物(PELA)被用来作为微囊的囊材。我们期望通过响应面法优化实验参数以得到最高包封率的微囊。然后应用并改良AnaJaklenec等的方法制备微囊支架用于促进骨愈合。制备的微囊其外粘附rhBMP-2,其内包封rhBMP-7。用这种微囊融合成的支架能序贯释放rhBMP-2及rhBMP-7,并进一步研究这种支架材料释放BMPs的控释曲线是否能模拟自然骨愈合过程中的BMPs的表达特点。最后通过大鼠的骨髓间充质干细胞体外实验及大鼠活体内骨缺损修复实验来进一步验证材料的成骨能力,降解能力及细胞毒性。目的：研发出制备过程简单,具有仿生控释BMP-2和BMP-7双因子的新型微球支架,其组成为BMP-2/PELA/BMP-7。研究这种生物活性支架的粘弹性,表征,生长因子的缓释曲线,缓释出因子的生物活性,生物相容性,诱导骨髓干细胞成骨转化能力。并进一步在活体骨缺损动物模型中评价其成骨能力,为今后骨缺损治疗的应用提供一种新方法和新思路。方法：1负载rhBMP-2的PELA微囊制备及包封率优化研究将改良复乳溶剂挥发法制备包封rhBMP-2的PELA微囊,通过单因素实验研究对微囊包封率影响最大的几个实验参数,包括PELA中PEG嵌段的分子量,投入的PELA的质量及rhBMP-2的质量,内水相乳化剂浓度,外水相中聚乙烯醇(PVA)浓度以及微囊固化的搅拌时间。并在单因素实验结果的基础上选取三个对包封率影响最大的实验参数,应用响应面法(RSM)设计三因素三水平的实验方案,优化实验参数,以得到最佳rhBMP-2包封率的PELA微囊。每个样本均设三个重复样以减少实验误差。采用DPS13.0软件行实验设计及其数据分析,包括回归分析及回归方程系数的确定。用STATISTICA7.0软件绘制的三维响应面图来描述自变量对于应变量的单个效应和交互效应。用方差分析检验该模型的统计学差异,P<0.05被认为差异有统计学意义。最后用扫描电镜观察制备的微囊的形态,大小,并计算平均直径。2BMPs/PELA微囊支架材料制备及体外缓释特点通过改良AnaJaklenec等的微囊支架制备方法,用PLA-PEG-PLA三嵌段共聚物(PELA)作为微囊的壁材,制备出其内包封rhBMP-7,其外粘附rhBMP-2的微胶囊,再通过二氯甲烷熏蒸将微囊聚合成支架材料(BMP-2/PELA/BMP-7,A组),并同时制备仅包封rhBMP-7的微囊支架(PELA/BMP-7,B组),仅表面粘附rhBMP-2的微囊支架(BMP-2/PELA,C组)以及单纯PELA微囊支架(D组)作为对照组。本部分首先研究这种支架材料的制备方法,将熏蒸不同时间的微囊支架在扫描电镜下观察其表面形态。再将微囊支架浸泡在37℃的PBS(pH7.4)缓冲液中,在固定的时间点称重,在46天内记录其干重及膨胀率的变化,以了解支架材料体外降解的性能。然后以同样的实验条件,在不同时间点检测缓释液中BMP-2及BMP-7的浓度,以获得微囊支架材料BMPs的释放曲线。每个样品均设六个重复样以减少误差。3BMPs/PELA微囊支架材料诱导大鼠骨髓间充质干细胞成骨分化实验研究采用体外培养的大鼠骨髓间充质干细胞(MSCs)来检验BMPs/PELA微囊支架材料的成骨能力及细胞毒性。先用全骨髓分离法分离培养大鼠骨髓间充质干细胞,通过观察细胞形态及流式细胞仪检测细胞表面标志物进一步鉴定；将微囊支架与MSCs共培养,在第3,7,14天时通过MTT实验来检验支架材料是否会对MSCs的活性产生影响；将支架的体外缓释液加入细胞培养液中诱导MSCs成骨分化,在第7天及14天时通过细胞碱性磷酸酶(AKP)活性测定来检验各组微囊支架缓释的BMP的生物活性及各种支架材料诱导MSCs成骨分化的能力。每组实验均设6个重复细胞样以减少误差,用SPSS19.0软件进行分析实验结果,组间数据比较用单因素方差分析法,两两比较用SNK法,P值小于0.05被认为差异有统计学意义。4BMPs/PELA微囊支架材料修复大鼠骨缺损实验研究建立大鼠股骨髁大块骨缺损的实验模型,于活体内检验各组微囊支架材料诱导成骨的能力。将24只大鼠随机分为4组,在全麻下于股骨髁部钻孔形成5×5mm圆形单层皮质的骨洞,分别用A,B,C,D四组支架材料填入检验其体内成骨效果。于术后第4周及第8周取材,观察大鼠的一般状况及股骨缺损部位软组织及骨组织的大体情况；将取下的股骨通过微型CT(mCT)扫描进行影像学分析以及新生骨各参数的测量；再通过HE染色及Masson三色染色法进一步行组织学观察评估新生骨生成及支架材料体内降解吸收情况。实验结果以SPSS19.0软件进行分析,组间数据比较用单因素方差分析,两两比较用SNK检验,P值<0.05被认为差异有统计学意义。结果：1负载rhBMP-2的PELA微囊制备及包封率优化研究本研究成功地以复乳溶剂蒸发法制备了包封rhBMP-2的PELA微胶囊。并应用了响应面法设计优化了制备过程中的关键参数以得到最高包封率的微囊。制备过程中一些参数,包括聚合物中PEG嵌段的分子量,投入原料的量(PELA和rhBMP-2),复乳液的搅拌固化时间,乳化剂的浓度及外水相中PVA浓度等均可影响最终微囊的rhBMP-2包封率,其最优值分别为：PEG嵌段分子量4000Da,PELA330mg,BMP3ìg,司盘-20浓度0.5%,聚乙烯醇浓度0.5%,在800rpm的转速下固化搅拌时间为30分钟。在单因素实验结果的基础上,我们选取PELA的质量(X1),PVA的浓度(X2)和rhBMP-2质量(X3)这三个对包封率影响最大的三个因素作为自变量,考察它们对包封率(Y)的影响。通过Box-Behnken法设计三因素三水平的响应面实验方案,得到的二次多项式回归方程模型：Y=27.33602+0.29086X1+86.69629X2-30.05764X3-0.00055X12-48.92593X22+2.79563X32+0.02088X1X3-7.36111X2X3方程回归的F值是47.8993,相应的P值为0.0001(<0.01),回归方程有统计学意义,其R2(决定系数)为0.9846。三维响应面图显示出PVA的浓度与rhBMP-2质量,及PELA的质量与rhBMP-2质量均对最后的包封率产生交互效应。以此模型推导出最佳的实验参数为：PELA的量为282.3mg,rhBMP-2量为1.0μg,PVA浓度为0.8%。在此条件下推算出微囊的包封率理论最大值为76.5%,而以此参数制备的微囊实际的rhBMP-2包封率为75.0%,仅存在微量误差,这也进一步证实该模型能很好地代表各变量值与响应值(包封率)之间的关系。电子显微镜下可以观察到大部分微囊都呈现出完整的圆球形,其直径分布基本均匀一致,微囊平均直径为40μm。2BMPs/PELA微囊支架材料制备及体外缓释特点参考第一部分得到的优化的实验参数,用改良的复乳溶剂挥发法成功制备出内包封BMP-7,其外粘附BMP-2的微胶囊。然后采用二氯甲烷熏蒸法能将其成功融合成多孔的三维支架材料(BMP-2/PELA/BMP-7)。并同时制备单因子材料组PELA/BMP-7和BMP-2/PELA,以及单纯PELA材料组微囊支架作为对照组。这种支架制备方法优点在于制备过程简单,一步到位,很大程度上避免了添加各种化学试剂及使用高温,紫外线照射交联等物理方法对功能蛋白质的结构和活性产生的负面影响,同时也节约了材料的制备成本。这种骨修复材料还具有一定的粘弹性,可以根据骨缺损的不同形态予以塑性。扫描电镜显示融合后的骨填充材料,其球状微囊的基本结构仍清晰可见。材料表面及剖面均可见一定的孔隙率,孔隙大小范围为50-200μm。体外降解实验表明：在磷酸盐缓冲液(PBS)中开始的两个星期,材料的膨胀率增加,于14日达到了峰值后开始快速。而材料干重在前两周无明显变化,2周后开始迅速下降,表明支架材料开始降解。缓释实验表明：BMP-2/PELA/BMP-7微囊支架材料能同时缓释BMP-2及BMP-7两大成骨活性因子,但释放曲线明显不同。rhBMP-2呈现出明显的初始突释现象,随后是一个较长时间的缓慢释放；而rhBMP-7初始突释也存在,但明显被抑制,而在缓慢释放2周后,其释放速度呈现出一个小幅加速。这两种活性因子的控释时间均长于42天。这些缓释特点能简单模拟自然骨折修复过程BMPs表达的特点。3BMPs/PELA微囊支架材料诱导大鼠骨髓间充质干细胞成骨分化实验研究用MTT法检测与各组支架材料共培养的大鼠MSCs的活性,结果表明,在共培养的第3天(F=2.08,P=0.135),第7天(F=O.12,P=0.954)及第14天(F=43.98,P=0.732)时,BMP-2/PELA/BMP-7组,PELA/BMP-7组和BMP-2/PELA组支架材料与纯细胞对照组的OD值差异无统计学意义。这结果表明BMP/PELA支架材料在短期内对大鼠MSC细胞活力没有负面影响,短期内未发现明显细胞毒性。再进一步通过检测细胞AKP活性来评价MSCs成骨分化的程度。于MSCs的培养液中加入各组支架材料缓释液,在培养第7天及14天的时候检测各组细胞的AKP活性,其结果表明：第7天时各组细胞AKP表达量有统计学差异(F=28.57,P=0.000),与纯PELA组相比,BMP-2/PELA/BMP-7组,PELA/BMP-7组和BMP-2/PELA组均能增强细胞的AKP活性,差异有统计学意义(P值均=0.000)而这三组之间AKP值差异无统计学意义。到第14天,各组细胞AKP值的总体差异也具有统计学意义(F=43.98,P=0.000)。BMP-2/PELA/BMP-7组中细胞的AKP活性高于PELA/BMP-7组和BMP-2/PELA组,差异有统计学意义(P=0.007及P=0.014)。PELA/BMP-7组和BMP-2/PELA组间差异无统计学意义。而纯PELA组的AKP值低于PELA/BMP-7组(P=0.000)和BMP-2/PELA组(P=0.007),其差异有统计学意义。这一结果提示：微囊支架缓释的的BMP仍保留很强的诱导干细胞成骨分化的生物学活性,而BMP-2和BMP-7这两大成骨活性因子对诱导MSCs成骨分化存在着协同效应。4BMPs/PELA微囊支架材料修复大鼠骨缺损实验研究大鼠手术后均成活,均未出现病理性骨折,伤口流脓,窦道等情况。取股骨标本,见所有标本周围软组织均未见脓肿,纤维化,骨组织未见死骨,大部分骨缺损生物材料填充处均可见新生组织覆盖良好。股骨mCT扫描可见：术后第4周时BMP-2/PELA/BMP-7组股骨缺损部位的新生骨较其他组多；PELA/BMP-7组和BMP-2/PELA组形成的新生骨多于纯PELA组。术后第8周各组骨缺损部的新生骨均多于第4周；BMP-2/PELA/BMP-7组的新生骨组织完全覆盖骨缺损部,并且可见新生骨皮质；PELA/BMP-7组和BMP-2/PELA组可见更多的新生骨；纯PELA组8周可在缺损部位发现少部分的新生骨组织.再通过mCT的几个重要参数来分析大鼠骨缺损部的新生骨,包括骨密度(BMD),新生骨骨体积分数(BV/TV),骨小梁数量(Tb.N)等,其结果表明：各组的骨密度值(BMD)术后第4周时各组差异有统计学意义(F=58.86,P=0.000),BMP-2/PELA/BMP-7组新生骨BMD均高于其他各组,差异有统计学意义(P值均=0.000),PELA/BMP-7组和BMP-2/PELA组骨密度值差异无统计学意义(P=0.54),但均高于PELA组(P=0.004；P=0.001)。而至术后第8周时,各组的BMD值总体差异也有统计学意义(F=22.57,P=0.000),BMP-2/PELA/BMP-7组新生骨BMD值仍高于PELA/BMP-7组(P=0.000),BMP-2/PELA组(P=0.039)和PELA组(P=0.000),差异有统计学意义。BMP-2/PELA组新生骨BMD均值要高于PELA/BMP-7组,差异则有统计学意义(P=0.013)。各组的BV-TV值在术后第4周时总体差异有统计学意义(F=171.34,P=0.000)。BMP-2/PELA/BMP-7组的BV/TV值高于其他各组,差异有统计学意义(P值均=0.000)。PELA/BMP-7组和BMP-2/PELA组的BV/TV值无统计学差异(P=0.134),而均高于PELA组,差异有统计学意义(P=0.044,P=0.001)。术后第8周时,各组总体差异有统计学意义(F=67.95,P=0.000),BMP-2/PELA/BMP-7的BV/TV值仍高于其他各组,差异有统计学意义(P值均=0.000)。BMP-2/PELA组BV/TV值高于PELA/BMP-7组,其差异有统计学意义(P=0.002)。而PELA组的BV/TV值在术后第8周仍低于BMP-2/PELA/BMP-7组(P=0.000),PELA/BMP-7组(P=0.035)和BMP-2/PELA(P=0.000),差异有统计学意义。各组新生骨的骨小梁数量(Tb.N)也体现出类似的趋势。术后第4周时方各组总体差异有统计学意义(F=20.35,P=0.000)。BMP-2/PELA/BMP-7组新生骨Tb.N值高于其他各组,差异有统计学意义(P值均=0.000)。PELA/BMP-7组和BMP-2/PELA组Tb.N值差异无统计学意义,而均大于PELA组(P=0.005及P=0.010),差异有统计学意义。术后第8周时各组总体差异也有统计学意义(F=12.01,P=0.000),BMP-2/PELA/BMP-7组Tb.N值仍高于PELA/BMP-7组(P=0.000)和PELA组(P=0.000),差异有统计学意义,而与BMP-2/PELA组比较无统计学差异。BMP-2/PELA组Tb.N值高于PELA/BMP-7组,差异有统计学意义(P=0.018)。这些数据结果表明：与BMP-2或BMP-7单因子缓释材料相比,BMP-2/PELA/BMP-7支架材料有更强大的活体内成骨诱导能力；而BMP-2/PELA支架材料中远期的成骨诱导能力强于PELA/BMP-7材料；所有复合BMPs的PELA微囊支架材料均比单纯的PELA聚合物显示出更好的成骨效果。大鼠股骨的组织学观察可见：BMP-2/PELA/BMP-7组在术后4周时可见新骨形成,其周围有残留的部分降解的材料和大量的血细胞,其切片上可见到共聚物材料降解并被新生骨替代的动态过程。到术后第8周,新生骨组织体积更大,而且成熟度更高,可以见到较多类似于成熟骨皮质的板层骨,在骨基质中可以发现带有血管的完整的骨单位。PELA/BMP-7组和BMP-2/PELA组在术后第4周时均能观察到部分新生的编织骨,骨基质中能发现支架材料的残留。到术后第8周,可以见到更多的新生骨组织,但没有发现成熟的板层骨及骨单位。而单纯的PELA组术后第4周无新骨形成,到术后第8周才能见少量新骨组织出现并开始取代支架材料。这些组织学结果与mCT图像结果相一致。结论：用改良的复乳溶剂挥发法成功制备出包封rhBMP-2的PELA微囊,并通过响应面法优化实验参数,微囊包封率达75%。然后用二氯甲烷熏蒸法将负载BMPs的微囊融合成多孔支架材料,可以简单模拟自然骨折修复过程BMPs表达的特点缓释BMP-2和BMP-7,其缓释时间超过42天。该材料具备较强的诱导大鼠间充质干细胞成骨分化的能力,并且生物相容性好,可完全降解,短期内未检测到细胞毒性。活体动物实验证实其能很好地修复大鼠股骨大块骨缺损,这为其今后广泛应用于临床治疗骨缺损提供了理论依据。
【作者】李夏林；
【导师】靳安民；
【作者基本信息】南方医科大学，外科学，2014，博士
【关键词】PLA-PEG-PLA；三嵌段共聚物；骨形态发生蛋白；微囊；支架；骨缺损；

【参考文献】
[1]申孟哲.摇荡声情而檠括于兴观群怨[D].北京大学,文艺学,2013,硕士.
[2]马荧晗.论宫本百合子与丁玲的女性解放意识[D].黑龙江大学,日语语言文学,2013,硕士.
[3]庞淑英,鲁明羽,陆玉昌.利用判定树归纳与面向属性归纳的组合方法对建筑业中施工企业的特征提取[J].计算机工程与应用,2003,16:208-211.
[4]张忠诚.中等职业学校体育教学信息化改革探讨[J].电子测试,2014,S2:214-216.
[5]高超,许震.氧化石墨烯液晶及宏观组装纤维[A].中国化学会.中国化学会第28届学术年会“中国—新加坡青年化学家论坛”摘要集[C].中国化学会:,2012:1.
[6]张新.应变式三维加速度传感器设计及相关理论研究[D].合肥工业大学,2008.
[7]张广建.论食品安全的刑法规制[D].苏州大学,刑法学,2013,硕士.
[8]郑娅.浓缩牛骨汤料工艺研究及营养风味物质分析[D].甘肃农业大学,农产品加工及贮藏工程,2013,硕士.
[9]汪天水.幼儿园自主游戏分享环节的教育功能研究[D].华东师范大学,学前教育学,2013,硕士.
[10]蒋朋.基于内在点强的加权网络演化模型的相关结果[D].南京航空航天大学,概率论与数理统计,2013,硕士.
[11]刘敏.延时相关的声触发与远程监控技术研究[D].西安工业大学,控制理论与控制工程,2014,硕士.
[12]王美科.面向除尘风机节能的转速闭环控制系统[D].浙江大学,控制科学与工程,2014,硕士.
[13]张泽静.寓中国传统文化于中学语文教学之中的策略探究[D].信阳师范学院,学科教学,2014,硕士.
[14]陈浩辰.基于微博的消费意图挖掘[D].哈尔滨工业大学,计算机技术,2014,硕士.
[15]林佳.类石墨相碳化氮的合成及其在乳腺癌肿瘤标记物检测中的应用[D].福建医科大学,2014.
[16]李凌俊.超高墩连续刚构桥施工期内风险评估研究[D].重庆交通大学,桥梁与隧道工程,2012,硕士.
[17]洪珊.从革命女儿到原始母亲[D].沈阳师范大学,中国现当代文学,2013,硕士.
[18]段良娟.主观幸福感、自我效能对老年人参与太极拳锻炼坚持性的影响研究[D].天津体育学院,体育教学（专业学位）,2013,硕士.
[19]任宏图.发展纸杆铅笔,保护生态环境——铅笔制造业的一场绿色革命[J].中国制笔.2001(03)
[20]王宇,杨莉.模糊k-prototypes聚类算法的一种改进算法[J].大连理工大学学报,2003,06:849-852.
[21]焦尚彬,刘丁,郑岗,张青.基于模糊逻辑方法的高压绝缘子污秽程度评定[J].电力系统自动化,2005,07:84-87+107.
[22]朱博华.基于J2EE的茂名市政协提案管理系统的设计与实现[D].电子科技大学,软件工程（专业学位）,2013,硕士.
[23]周海峰.油井射孔段管柱对爆炸冲击载荷动态响应研究[D].北京理工大学,兵器科学与技术,2014,硕士.
[24]赵军.春秋时期晋国战争史研究[D].山西师范大学,专门史,2013,硕士.
[25]郭雯.试论巴赫《康塔塔》的艺术风格与演唱技巧[D].西安音乐学院,声乐演唱,2013,硕士.
[26]莫以为,萧德云.基于进化粒子滤波器的混合系统故障诊断[J].控制与决策,2004,06:611-615.
[27]郎洪明.流量计式饲料油脂添加设备的系统设计分析[J].饲料广角.2007(09)
[28]蓝云春.《德伯家的苔丝》中的现代主义主题[D].江西师范大学,英语语言文学,2004,硕士.
[29]薛永生,翁伟,文娟,王劲波,张宇.LSNCCP——一种基于最大不相含核心点集的聚类算法[J].计算机研究与发展,2004,11:1930-1935.
[30]杨浩楠.外国人眼中的英国：试论四本非小说作品中的英国形象[D].西安外国语大学,英语语言文学,2013,硕士.
[31]陈昌根.中国电信企业大客户关系管理[D].广西大学,2004.
[32]许兆新,周双娥,郝燕玲.最优关联规则的形式和挖掘思想的研究[J].计算机工程与应用,2001,20:118-119.
[33]张成良.新型城域网MSTP的关键技术[J].电信网技术.2003(07)
[34]邓彤.铁路货车检修生产安全风险管理研究[D].西南交通大学,项目管理（专业学位）,2013,硕士.
[35]开金丹.MCM3、ki67及PCNA在非小细胞肺癌中的表达及临床意义[D].华中科技大学,胸外科,2013,硕士.
[36]古新.管壳式换热器数值模拟与斜向流换热器研究[D].郑州大学,2006.
[37]倪进东.芜湖市区计划免疫服务满意度调查研究[D].安徽医科大学,社会医学与卫生事业管理,2003,硕士.
[38]杨伟才.基于FPGA的嵌入式图形处理系统的设计[D].河北科技大学,测试计量技术及仪器,2013,硕士.
[39]李虹儒.基于居民归属感的居住区景观设计[D].中南林业科技大学,风景园林,2013,硕士.
[40]李玲玲.基于颜色Petri网的再制造系统不确定性建模与优化调度方法[D].重庆大学,工业工程（专业学位）,2014,硕士.
[41]戴培华.Pd/LaNiO_3复合薄膜电极材料在碱性溶液中的电化学行为研究[D].兰州理工大学,材料学,2014,硕士.
[42]杨卫平,黄烟波,段丹青.智能入侵检测技术述评[J].网络安全技术与应用,2005,04:58-59+57.
[43]石纯一,王家廞,宋星,万方.运输调度专家系统[J].自动化学报,1986,03:245-251.
[44]苏茹茹.齿轮泵困油问题的理论分析与仿真试验研究[D].河北工程大学,机械设计及理论,2012,硕士.
[45]王艺潼.论行政处罚裁量权的规制[D].吉林大学,宪法与行政法学,2014,硕士.
[46]董伯骞.油松人工林天然更新特征对林隙和皆伐年限的响应[D].北京林业大学,森林培育,2014,博士.
[47]王彪,孙国梁.爆炸法消除应力技术在压力容器制造中的应用[J].石油化工设备.1993(03)
[48]温钧明.嵌入式智能数字测试仪的研究[D].河北经贸大学,计算机应用技术,2013,硕士.
[49]TchisterMorrelEBISSA.氮沉降背景下枯落物对草地生产力的影响[D].东北师范大学,生态学,2012,硕士.
[50]郭以贺,谢志远,徐智坚.基于ARM的馈线自动化系统设计[J].电力系统自动化,2010,10:105-108.