基于转录组测序分析美洲大蠊提取物促进小鼠创面愈合的分子机制

陈佳松， 陈峰， 彭锐， 耿福能， 邹方东*

(1. 四川大学生命科学学院，成都610064； 2. 四川省药用动物工程技术研究中心，成都610064；3. 药用美洲大蠊四川省重点实验室，成都610064)



基于转录组测序分析美洲大蠊提取物促进小鼠创面愈合的分子机制

陈佳松1， 2， 3， 陈峰1， 2， 3， 彭锐1， 2， 3， 耿福能3*， 邹方东1， 2， 3*

(1. 四川大学生命科学学院，成都610064； 2. 四川省药用动物工程技术研究中心，成都610064；3. 药用美洲大蠊四川省重点实验室，成都610064)

美洲大蠊Periplanetaamericana是传统的中药材，其提取物已在临床上广泛用于创面治疗，但对其中的分子机制知之甚少。本文通过构建小鼠皮肤创伤模型，并基于转录组测序(RNA-seq)深入分析美洲大蠊提取物促进小鼠创面愈合的分子机制。首先构建C57小鼠皮肤全层切除模型，并将其分为3组，其中实验组分别以2种美洲大蠊提取物(康复新液、精粉)为敷料，而对照组以75%乙醇为敷料。用药处理3 d后取伤口皮肤组织送样进行RNA-seq，然后分析转录组数据找出与伤口愈合相关的差异性表达基因，并利用荧光定量PCR(Q-PCR)对伤口组织中相关基因的表达量进一步验证。实验结果显示，实验组小鼠伤口的结痂速度比对照组快，这表明了小鼠创面模型中，美洲大蠊提取物能促进伤口愈合。经转录组分析，得到显著性差异表达基因数：对照组vs.精粉组为545，对照组vs.康复新液组为938。结合生物信息分析结果和已发表的文献，发现3个可能参与调节伤口愈合的关键基因：表皮调节素(Ereg)、Gli-kruppel家族成员(Gli2)和表皮型转谷氨酰胺酶3(Tgm3)。Q-PCR实验表明，这3个基因在美洲大蠊提取物加药组中均高表达。本文研究结果表明，2种美洲大蠊提取物：康复新液、精粉可能通过诱导Ereg、Tgm3和Gli2的高表达来促进皮肤创面的愈合。

美洲大蠊；伤口修复；转录组分析；关键基因

美洲大蠊Periplanetaamericana是蜚蠊目Blattaria蜚蠊科Blattidae昆虫，主要分布于热带和亚热带，是世界性卫生害虫，除传播多种人类疾病外，其分泌物还引起人类过敏性反应。此外，美洲大蠊还是传统的中药材，常以干燥或鲜成虫入药，具有抗肿瘤、抑菌、抗氧化、增强免疫、抗炎、镇痛等作用。例如，在子宫内膜癌细胞株JEC中，美洲大蠊提取物能够通过阻滞G2/M期抑制肿瘤增殖(张晓巍，朱艳，2015)。美洲大蠊70%乙醇提取物：异香豆素类化合物，对肝癌细胞(HepG2)和人乳腺癌细胞(MCF-7)也有显著的细胞毒性(Luoetal.，2014)。此外，美洲大蠊精粉对金黄色葡萄球菌Staphylococcusaureus、大肠埃希菌Escherichiacoli、铜绿假单胞菌Pseudomonasaeruginosa均有抑菌作用(李远辉等，2014)。

康复新液由美洲大蠊醇提物配制而成，临床上广泛用于创面治疗。在对肛管难愈合的79例临床观察中，发现以康复新液进行局部换药时，与对照组利凡诺溶液局部换药相比治愈率高、伤口愈合时间明显较短(伍伟明等，2015)。另外，在大鼠皮下切割伤实验中，康复新液处理组伤口愈合率也高于对照组(张俊等，2014)。有实验表明美洲大蠊提取物可以加速全身放射损伤的愈合。作者使用美洲大蠊多元醇提取物W11-a12作为治疗药物发现，其能通过增加细胞外基质表达量、调控伤口细胞增值、使肉芽组织的形成速度加快，从而提高伤口愈合速度和质量(舒崇湘等，2001)。

临床及实验证据都表明美洲大蠊提取物有明显的促创面愈合能力，且能通过多种途径促进创面愈合。但是，目前对其中涉及的分子机制不甚明确。基于此，本文利用转录组测序(RNA-seq)这一成熟的基因组学研究技术手段，分析美洲大蠊提取物促进创伤修复过程中基因表达变化，从中筛选出与创面愈合相关的差异性表达基因，并进一步从细胞和个体水平进行验证，从而初步揭示美洲大蠊提取物促进创伤修复的分子机制，为美洲大蠊的深入开发奠定基础。

1 材料与方法

1.1材料与试剂

SPF级C57小鼠(四川大学实验动物中心)，实验动物生产许可证：SCSK(川)2013-026，实验动物使用许可证：SYSK(川)2013-185；康复新液、精粉(四川好医生攀西药业有限公司)；RNAiso Plus(TaKaRa，Code No. 9108)；PrimeScriptTMRT reagent Kit with gDNA Eraser(TaKaRa，RR047A)；EvaGreen 2X qPCR Master Mix(abm，Master Mix-S)；Tissue Protein Extraction Kit(康为，CW0891S)。

1.2方法

1.2.1动物模型实验5～7周龄的C57实验小鼠称量后选取体质量在平均值±3 g内的小鼠，随机分为3组：对照组(Control)、康复新液实验组(KFX)、精粉实验组(JF)(用75%乙醇配制为0.05 g·mL-1)。每组3只。饲养1 d后，将小鼠背部脱毛备用。首先用无水乙醚麻醉小鼠，然后参照《现代创伤修复学》所述(付小兵，王德文，1997)在背部用消毒手术剪，剪下边长为1 cm的正方形创口。用75%乙醇擦拭伤口消毒后，每天早晚擦药各1次，每次1 mL溶液，对照组抹75%乙醇。拍照记录，利用Photoshop计算伤口面积变化。

1.2.2转录组分析造模3 d后，取伤口周围约2 mm区域的皮肤送公司进行RNA-seq。原始reads经过滤后用tophat2和cufflink对reads进行组装拼接。用R语言包：EdgeR分析得到显著性差异基因，然后用blast2go对差异基因进行Gene Ontology(GO)和kyoto encyclo pedia of genes and genomes(KEGG)分析。

1.2.3荧光定量PCR实验使用试剂盒RNAiso Plus提取皮肤组织总RNA，步骤按照试剂所提供的说明书；反转试剂盒将总RNA反转为cDNA；从NCBI上查找Ereg、Gli2、Tgm3、GAPDH序列，用Primer 5设计引物如表1；荧光定量PCR(Q-PCR)步骤和参数设置按照EvaGreen 2X qPCR Master Mix说明书所示。

表1 荧光定量PCR引物Table 1 Primer used in this study

2 结果

2.1动物模型实验

构建C57小鼠创伤模型，实验组每天涂康复新液或精粉液，对照组每天涂75%乙醇，隔一天拍照记录。统计实验组和对照组伤口面积变化，结果表明美洲大蠊提取物：康复新液和精粉能加速伤口面积收缩(图1)。

图1 伤口面积变化Fig. 1 Wound area in mice

2.2转录组分析

2.2.1伤口组织测序结果由前期的观察实验发现，用药3 d左右加药组已有结痂现象。因此，药物处理3 d后分别选择了3组(对照组、精粉组和康复新液组)伤口组织进行RNA-seq。RNA提取、文库构建后，通过Illumina 高通量测序平台得到的原始数据用raw reads表示。各组样品得到的原始reads如表2所示。

表2 测序原始readsTable 2 Raw data of RNA-seq

2.2.2显著性差异表达基因分析完成原始reads过滤、拼接和组装后用R语言的软件包：edgeR，可以得到实验组和对照组基因差异表达情况，FDR与log2FC为筛选差异基因参数，筛选条件为FDR<0.05且|log2FC|>1，满足条件即为显著性差异表达基因，结果见表3。

表3 差异基因数目统计Table 3 Number of differentially expressed genes

2.2.3差异基因注释在得到差异基因后，对其进行GO富集分析。GO是通用的基因功能分类体系，用来全面描述基因和基因产物的属性。GO共有3个本体(ontology)，分别描述基因的分子功能、细胞组分、参与的生物过程(Ashburneretal.，2000)。通过GO富集分析可以得到的差异基因富集情况如图2。

富集结果表明：康复新液、精粉2个实验组，差异基因主要富集于生物学过程中的胞间信号、发育过程；细胞成分中的细胞外组分；分子功能中的结构分子激活。暗示这2种药物可能共同作用于上面4个GO term中的基因。结合已发表的文献，筛选出3个与创面愈合相关的关键基因：表皮调节素(Ereg)、Gli-Kruppel家族成员Gli2(Gli2)、表皮型转谷氨酰胺酶3(Tgm3)。

2.2.4差异基因KEGGpathway富集为了进一步了解药物作用的信号通路，将差异基因进行KEGG pathway富集。KEGG是一个综合数据库，通过它可以将基因和细胞生化过程联系起来(Kanehisa & Goto，2000)。实验组差异基因富集显著性最高的10个通路详见表4。

结果表明，康复新液影响最显著的通路为细胞因子受体相互作用。已有实验表明这条通路上基因在皮肤创伤早期表达量变化明显(Royetal.，2008)。

第二条通路为趋化因子信号途径，该通路可以调节炎症反应和伤口愈合(Martin & Leibovich，2005)。PPAR信号通路中基因缺失会导致伤口愈合延迟(Michalik & Wahli，2006)。在精粉实验组中，变化最明显的通路是花生四烯酸代谢。花生四烯酸是PGE2的前体，而PGE2能刺激伤口愈合(Lietal.，2010)。精粉也能调控细胞因子受体相互作用通路从而影响伤口的愈合。其他通路与伤口愈合的联系有待进一步研究。

2.3Q-PCR验证实验

加药3 d后提取伤口皮肤组织RNA，利用Q-PCR技术检测对照组和实验组中3个关键基因的表达情况。

Q-PCR结果如图3，这3个基因在实验组中均高表达。已有研究表明，这几个基因高表达能促进伤口的愈合。Ereg在精粉实验组中表达量最高，这和前期观察到精粉组创面愈合最快的结果相符。

图2 差异基因GO富集Fig. 2 GO enrichment analysis for DEGs

3 讨论

转录组是连接基因组遗传信息与生物功能的必要纽带，转录组研究已经成为揭示生物生长发育调控和逆境胁迫适应机制、生物进化规律、疾病发生发展的重要机制以及发现致病基因调控关键靶点等方面的最佳研究手段。通过转录组分析，找出在加药组中显著性差异(FDR≤0.05)表达的基因。精粉组和对照组相比：有545个基因差异表达，其中在精粉组中表达量升高的有421个，表达量降低的有124个；康复新液组和对照组相比：差异基因总数为938个，高表达739个，低表达199个。

结合已发表的文献，从2组共同的差异性表达基因中筛选出了3个与创面愈合相关的关键基因：Ereg、Tgm3和Gli2。Ereg是表皮调节素家族中的一员，而表皮调节素在皮肤损伤愈合中是关键的调节基因。Ereg可以通过增加细胞增殖和迁移来加快体外伤口的愈合，而且Ereg加速伤口愈合的效率比2种公认的促进伤口愈合因子：表皮生长因子、转化生长因子a更高(Draperetal.，2003)。这表明Ereg可能在美洲大蠊提取物促进伤口愈合中发挥着重要作用。此外，有实验显示Ereg还与炎症反应、肝祖细胞生长和结肠炎相关的肿瘤细胞生长等方面相关(Zhangetal.，2012；Neufertetal.，2013；Tomitaetal.，2014)。因此，Ereg的作用值得深入研究。Tgm3对皮肤功能的完整有重要贡献，在表皮角质化包膜形成过程中参与结构蛋白的交联。Tgm3缺陷小鼠中，皮肤炎症会增加(Bognaretal.，2014)。Gli2是一种转录因子，参与Hn信号通路的调节，已有实验表明Gli2能促进细胞增殖(Panetal.，2009)。

表4 差异基因KEGG pathway富集Table 4 KEGG pathway enrichment results

图3 Q-PCR检测关键基因表达量Fig. 3 Key genes expression in wound tissues

Q-PCR结果进一步证实在创伤组织中Gli2、Ereg、Tgm3这3个基因在实验组中均高表达，其中Ereg在精粉组中表达量最高，而Ereg高表达可以促进细胞增殖，暗示前期观察中精粉组结痂速度快，可能与Ereg表达量高有关。这表明了2种美洲大蠊提取物：康复新液、精粉都能提高伤口组织中Ereg、Gli2、Tgm3这3个参与伤口愈合分子通路的基因的表达量，从而加速皮肤伤口的愈合。此外，精粉为美洲大蠊干燥后磨成的粉末，擦药后会覆盖在伤口表面，还可以阻止外来细菌等有害微生物的侵袭，提供了伤口修复过程中细胞迁移的路径等。

近年来，对美洲大蠊的药物作用机制已有一些相关研究。谭巧云等(2016)的研究表明，美洲大蠊提取物可以降低机体的炎症因子(CD3+、CD4+和CD8+)水平和通过IL-2调节免疫功能，加快大鼠口腔黏膜愈合。杨胜群等(2016)认为美洲大蠊水提取物可能调控TGF-β1的释放促进伤口愈合又避免形成瘢痕。曾娟妮等(2016)实验表明美洲大蠊提取液能提高伤口VEGF表达量来促进创面愈合。目前，对美洲大蠊促进创面愈合的分子机制了解仍然很少，而本文基于小鼠创面愈合模型，初步探讨了其中的分子机制。基于RNA-seq分析技术，筛选了3个与美洲大蠊促进创面愈合相关的关键基因。本文结果表明，美洲大蠊提取物——康复新液、精粉能通过刺激伤口组织中3个关键基因：Gli2、Ereg、Tgm3加速伤口愈合。本研究结果可为美洲大蠊临床研究及深入开发奠定理论基础。

致谢：本实验得到四川好医生药业集团的大力支持，在此致以最诚挚的谢意！

付小兵， 王德文. 1997. 创伤修复基础[M]. 北京： 人民军医出版杜： 81.

李远辉， 李希， 冯建安， 等. 2014. 美洲大蠊不同提取部位的体外抑菌活性研究[J]. 中药与临床， 5(6)： 27-29.

舒崇湘， 程天民， 阎国和， 等. 2001. 6Gy全身照射对皮肤伤口几种细胞外基质成分的影响及W11-a12的促愈作用[J]. 中华创伤杂志， 17(10)： 604-607.

谭巧云， 满红霞， 那凯歌， 等. 2016. 美洲大蠊提取物对口腔溃疡大鼠模型的作用[J]. 中国临床药理学杂志， 32(11)： 1014-1016.

伍伟明， 金丽颖， 王真权， 等. 2015. 康复新液治疗肛管难愈合创面的临床研究[J]. 中医药导报， 21(22)： 71-73.

杨胜群， 熊茜， 邹俊波， 等. 2016. 美洲大蠊水提取物对小鼠烫伤模型的影响[J]. 中国实验方剂学杂志， 22(22)： 145-149.曾娟妮， 刘筱， 伍伟明， 等. 2016. 美洲大蠊提取液对大鼠难愈合创面VEGF表达影响的研究[J]. 中国民族民间医药杂志， (19)： 64-66.

张俊， 孟令贺， 单士军， 等. 2014.康复新液对实验性大鼠皮肤切割伤痂下愈合的影响[J]. 天津医科大学学报， 20(3)： 192-195.

张晓巍， 朱艳. 2015. 苦参总碱和美洲大蠊提取物对孕激素受体阴性的子宫内膜癌细胞株JEC的作用研究[J]. 中国中药杂志， 40(11)： 2210-2213.

Ashburner M， Ball CA， Blake JA，etal. 2000. Gene ontology： tool for the unification of biology[J]. Nature Genetics， 25(1)： 25-29.

Bognar P， Nemeth I， Mayer B，etal. 2014. Reduced inflammatory threshold indicates skin barrier defect in transglutaminase 3 knockout mice[J]. Journal of Investigative Dermatology， 134(1)： 105-111.

Draper BK， Komurasaki T， Davidson MK，etal. 2003. Epiregulin is more potent than EGF or TGFα in promotinginvitrowound closure due to enhanced ERK/MAPK activation[J]. Journal of Cellular Biochemistry， 89(6)： 1126-1137.

Kanehisa M， Goto S. 2000. KEGG： kyoto encyclopedia of genes and genomes[J]. Nucleic Acids Research， 28(1)： 27-30.

Li F， Huang Q， Chen J，etal. 2010. Apoptotic cells activate the “phoenix rising” pathway to promote wound healing and tissue regeneration[J]. Science Signaling， 3(110)： ra13.

Luo SL， Huang XJ， Wang Y，etal. 2014. Isocoumarins from American cockroach (Periplanetaamericana) and their cytotoxic activities[J].Fitoterapia， 95： 115-120.

Martin P， Leibovich SJ. 2005. Inflammatory cells during wound repair： the good， the bad and the ugly[J]. Trends in Cell Biology， 15(11)： 599-607.

Michalik L， Wahli W. 2006. Involvement of PPAR nuclear receptors in tissue injury and wound repair[J]. The Journal of Clinical Investigation， 116(3)： 598-606.

Neufert C， Becker C， Türeci Ö，etal. 2013. Tumor fibroblast-derived epiregulin promotes growth of colitis-associated neoplasms through ERK[J]. The Journal of Clinical Investigation， 123(4)： 1428-1443.

Pan Y， Wang C， Wang B. 2009. Phosphorylation of Gli2 by protein kinase A is required for Gli2 processing and degradation and the Sonic Hedgehog-regulated mouse development[J]. Developmental Biology， 326(1)： 177-189.

Roy S， Khanna S， Rink C，etal. 2008. Characterization of the acute temporal changes in excisional murine cutaneous wound inflammation by screening of the wound-edge transcriptome[J]. Physiological Genomics， 34(2)： 162-184.

Tomita K， Haga H， Mizuno K，etal. 2014. Epiregulin promotes the emergence and proliferation of adult liver progenitor cells[J]. American Journal of Physiology-Gastrointestinal and Liver Physiology， 307(1)： G50-G57.

Zhang Y， Kobayashi T， Hayashi Y，etal. 2012. Important role of epiregulin in inflammatory responses during corneal epithelial wound healing[J]. Investigative Ophthalmology & Visual Science， 53(4)： 2414-2423.

RegulatoryMolecularMechanismofPeriplanetaamericanaExtractsonWoundHealingUsingRNA-seq

CHEN Jiasong1，2，3， CHEN Feng1，2，3， PENG Rui1，2，3， GENG Funeng3*， ZOU Fangdong1，2，3*

(1. College of Life Sciences， Sichuan University， Chengdu 610064， China；2. Sichuan Engineering Research Center for Medicinal Animals， Chengdu 610064， China；3. Sichuan Key Laboratory for Medicinal American Cockroach， Chengdu 610064， China)

Periplanetaamericanais a traditional Chinese medicine， and its extracts have been widely used in the treatment of wound in clinic. However， the molecular mechanism is poorly understood. In this study， to investigate the regulatory mechanism ofP.americanaextracts on wound healing， a mouse model with excision wound was constructed to assess the efficacy ofP.americanaextracts on wound therapy. A total of 9 mice were randomly divided into 3 groups： control group (treated with 75% ethanol solution) and 2P.americanaextracts treated groups (Kangfuxin and Jingfen)， respectively. RNA-seq was used to screen the key genes in wound skin tissues on day 3. Furthermore， real-time quantitative PCR was performed to confirm the expression of key genes. The results showed that bothP.americanaextracts could accelerate wound healing. RNA-seq analyses showed that， comparison to the control group， a total of 545 differentially expressed genes (DEGs) were identified in Jingfen treated group and 938 DEGs in Kangfuxin treated group. Bioinformatics analysis revealed that 3 key genes， includingEreg，Gli2 andTgm3， may play important roles in wound healing during the treatment ofP.americanaextracts. Q-PCR experiments further confirmed thatEreg，Gli2 andTgm3 were up-regulated after the treatment of bothP.americanaextracts. Collectively， our results indicated that the extracts ofP.americanamay accelerate wound healing by regulatingEreg，Gli2 andTgm3 expression.

Periplanetaamericana； wound healing； RNA-seq； key genes

2017-03-31接受日期：2017-05-09

陈佳松(1991—)， 男， 硕士研究生， 研究方向为细胞生物学， E-mail：chenjiasong91@qq.com

*通信作者Corresponding author， E-mail：haoyishenggfn@126.com； fundzou@scu.edu.cn

10.11984/j.issn.1000-7083.20170106

Q95-33； Q786

： A

： 1000-7083(2017)04-0398-06