下载:
下载:
-
中药红花（Carthami Flos）是菊科植物红花（Carthamus tinctorius L.）的干燥花，传统本草学著作《本草纲目》记载，红花具有活血散瘀，通经止痛的功效[1]，其药材和制剂在临床上被广泛用于心脑血管疾病的预防和治疗。现代药理研究表明，其主要药效物质是以羟基红花黄色素A（hydroxysafflower yellow A，HSYA）为代表的查尔酮类化合物和以菸花苷为代表的黄酮醇类化合物，这些化合物均具有良好的心脑血管损伤保护活性[2-3]。红花药材的产量偏低，每平方千米产量仅为18.0～22.5 t[4]，其中特有的HSYA[5]、红花红色素等查尔酮类成分在不同品种间差异较大[6]。由于红花中的查尔酮类成分仅特异性地存在于花冠中[7]，加之体外组织培养再生率低[8]等原因，对其功能基因的研究工作一直进展缓慢。特别是对于HSYA等红花特有的有效成分，其生物合成相关的功能基因尚不完全清楚，合成通路也未被完全解析[9]。因此，用现代分子生物学技术手段以提高药效物质的含量，是提高红花品质，节约土地资源、降低制药成本的一条新途径。
短链脱氢酶/还原酶（short-chain dehydrogenases/reductases，SDR）在植物次生代谢物的生物合成中广泛参与各类碳-氧双键，碳-碳双键以及烯酮键的氧化还原催化反应。根据SDRs基因序列的特征结构，SDRs超家族可以被分为5个亚家族[10-14]。最早发现并且进行鉴定的两类主要短链还原酶命名为classical和extend，classical类的SDRs基因拥有长度约为250个氨基酸残基，被称为Extended类的SDRs基因在碳基末端因其含有多余的约100个氨基酸残基而得名。另外3种类型SDRs基因分别被命名为intermediate、complex和divergent。这些类型的SDRs基因基于其结合辅酶类型和结合催化位点的不同进行命名分类。此外，SDRs存在与传统类型不同的含有“rossmann-fold”保守结构域的氧化还原酶结构[15-18]。
黄酮类化合物起源于莽草酸途径和苯丙素生物合成途径，1个香豆酰辅酶A（coumaroyl CoA）和3个丙二酰辅酶A（malonyl CoA）在查尔酮合酶的作用下生成二氢查尔酮，然后经查尔酮异构酶催化为二氢黄酮，进一步在各类还原酶，聚合酶和糖基转移酶的作用下，生成终端次生代谢产物组合[19-21]。红花中所含的主要有效成分HSYA具有查尔酮式结构，本课题组前期研究认为：HSYA从前体物质到合成，中间存在必不可少的氧化还原过程。短链脱氢还原酶家族广泛参与植物体内次生代谢，这一类还原酶都带有相似的折叠结构以及催化位点，已有研究表明，其对苯丙烷代谢途径起重要作用[22-23]，但有关红花中还原酶基因相关报道较少[24]。故笔者通过对红花转录组数据库、基因表达谱数据库以及代谢组数据库进行分析，筛选在HSYA生物合成途径的关键还原酶基因，并进行功能验证，以期揭示红花次生代谢成分生物合成途径，为定向调控红花的品质提供科学依据。
-
云南巍山红花品系（ZHH0119），采自海军军医大学药学系温室，经海军军医大学郭美丽教授鉴定为菊科植物红花（Carthamus tinctorius L.）。红花种植条件：温度恒定25 ℃，16 h光照，8 h黑暗。采集相关花与组织后迅速存放于液氮或者−80 ℃冰箱中冷冻。
-
按照Trans ZOL Plant植物总RNA提取试剂盒（北京全式金公司，中国）说明书方法提取红花花冠总RNA，按照Transtart One-Step gDNA Removal and cDNA Synthesis Super Mix逆转录试剂盒（北京全式金公司，中国）说明书方法进行cDNA第一链的合成。cDNA于−20 ℃保存。
-
基于数据库中的基因注释以“黄酮还原酶”和“黄酮类化合物生物合成”作为关键词进行检索，筛选出其中可能与HSYA生物合成相关的还原酶基因，将筛选基因不同花期时间的表达量，将其与红花代谢组数据库中同花期的芦丁（rutin）、山柰酚（kaempferol）、槲皮素（quercetin）、HSYA、柚皮素（naringenin）、山柰酚-3-O-芸香糖苷（kaempferol-3-O-rutinoside）、山柰酚-3-O-葡萄糖苷（kaempferol-3-O-gluciside）、Carthamin、芹菜素（apigenin）、黄芩素（scutellarein）、木犀草素（luteolin）、苯丙氨酸（D-phenylalanine） 12个主要成分的含量[12,25]进行皮尔森相关性分析。
-
基于红花花冠EST转录组文库，结合第三代测序技术[26-29]红花花冠全长转录组数据库筛选得到目的基因序列。在其5'端、3'端分别设计特异性引物。按照2× Phanta Flash Master Mix（Dye Plus）高保真酶（南京诺唯赞公司，中国）说明书进行PCR扩增，扩增片段经EasyPure Quick Gel Extraction Kit胶回收试剂盒（北京全式金公司，中国）说明书操作回收后，连接于pEASY-Blunt Zero Cloning Kit（北京全式金公司，中国）载体上，转化至大肠杆菌T1感受态细胞（北京全式金公司，中国）后，涂布在LBA平板上，恒温培养37 ℃过夜，挑取阳性单克隆菌落[30-31]，送至上海生工生物有限公司进行菌液测序。
用ExPASyProtParam工具（http://web.expasy.org/compute/）对目的基因的理论等电点（pI），蛋白分子量（MW）和蛋白分子式进行预测。通过Simple Molecule Architecture Research Tool工具（http://smart.embl-heidelberg.de/）对目的基因编码的蛋白质结构功能域进行分析。使用ProtScale（http://us.Expasy.org/cgi-bin/protscale.pl）以及TMHMM（http://www.cbs.dtu.dk/services/TMHMM/）对蛋白质的亲/疏水性和跨膜区域做出预测。使用SignaIP 4.0（http://www.cbs.dtu.dk/services/SignalP/）预测目的蛋白是否含有信号肽。使用NCBI BLAST(https://blast.ncbi.nlm.nih.gov/Blast.cgi)对筛选出的SDRs基因进行BLAST序列比对。通过Neighbor-Joining相邻节点法构建系统发育进化树，自展分析法进行1000次重复[32-34]。使用PBILYON-GRLAND数据库预测构建蛋白质二级结构模型。蛋白三级结构由Protein Homology/analogy Recognition Engine预测。用WOLFPSORT软件（https://wolfpsort.hgc.jp/）进行亚细胞定位预测。
-
取盛花期新鲜红花根、茎、叶、花冠4个部位的新鲜组织和花期Ⅰ（开花前3 d）、花期Ⅱ（开花当天）、花期Ⅲ（开花后1 d）、花期Ⅳ（开花后3 d）4个花期的新鲜花冠，提取总RNA，合成cDNA第一链后，在靠近5'端处对各个基因设计引物，依据Transtart Top Green qPCR super Mix（北京全式金公司，中国）试剂盒推荐体系，以Ct60s（KJ634810）作为内参标记基因，进行qRT-PCR实验，结果使用2−ΔΔCt的方法进行计算分析[35]。
-
根据CtSDR3的开放阅读框和植物真核表达载体pMT-39序列信息，设计无缝克隆引物。以红花cDNA做模板，使用高保真酶进行PCR反应。产物经胶回收后依无缝克隆试剂盒说明书与经NcoI酶切线性化的pMT-39载体进行重组连接。重组载体转化大肠杆菌T1感受态细胞，挑取阳性克隆菌株扩大培养后抽提质粒，提取的pMT39-CtSDR3质粒用冷冻法转至农杆菌GV3101中。LBK+Rif平板筛选阳性克隆后，取1ml OD600 = 0.8的菌液经6 000 r/min，离心3 min后用1 ml 5%蔗糖溶液重悬，加入Silwet-L 1μl，用注射器注射于红花花柱，套袋避光[35]。
在pMT-39的35 s启动子区域设计5'端特异性引物，在目的基因CtSDR3中设计3'端引物。取T2代新鲜叶cDNA第一链作为模板，2× Easy Taq PCR Mix（北京全式金，中国）推荐体系进行PCR反应，确定是否存在目的条带。采集CtSDR3阳性植株花冠以及pMT-39空载体对照植株的花冠，按照上述的qRT-PCR反应体系评价CtSDR3基因的过表达水平，使用UPLC-Q-TOF/MS 检测CtSDR3过表达组和空载体对照组的黄酮代谢物含量，选择以HSYA为代表性成分的8个黄酮类化合物作为检测对象。
-
根据CtSDR3的开放阅读框及蛋白表达载体pGEX-6p-1以及pET-28a序列信息，设计同源重组克隆引物[34]。以红花花冠cDNA为模板，使用高保真酶进行PCR反应。PCR产物经胶回收后依无缝克隆试剂盒说明书与经XhoI、BamHI酶切线性化的载体pGEX-6p-1以及pET-28a进行重组连接。重组载体转化大肠杆菌T1感受态细胞，挑取阳性菌株克隆扩大培养后抽提质粒，提取的重组质粒用热激法转至大肠杆菌Rosseta（DE3）（上海唯地生物，中国）中。
在20 ml LBA液体培养基中培养至OD600为0.6左右，分2份10 ml菌液各加入终浓度为0.3 mmol/L的IPTG和生理盐水。恒温培养箱中16 ℃，100 r/min继续培养16 h[35-36]。菌液离心弃上清液，用1×PBS缓冲液洗涤两次后重悬。超声破碎仪中40 kW，工作时间5 s，循环间隔时间25 s，共15个循环进行破碎[16]，裂解完成后取上清与沉淀15 μl，上样检测。
-
通过分析，得到contig325、contig483、contig2863共3个与HSYA具有强相关性的基因（r>0.85），见图1。
图 1 不同花期红花还原酶基因表达量与黄酮类化合物积累量相关性分析热图
-
3个目的基因序列信息经测序验证结果如下：contig325全长共1523 bp，开放阅读框1341bp，编码446个氨基酸；contig483全长1393 bp，开放阅读框792 bp，编码263个氨基酸；contig2863全长序列1527 bp，开放阅读框1023 bp，编码340个氨基酸。PCR产物电泳结果如图2所示。
图 2 PCR产物电泳图
contig325基因编码446个氨基酸，命名为CtSDR1（GenBank登录号：MW792035）；Contig483基因编码263个氨基酸，命名为CtSDR2（GenBank登录号：MW792036）；Contig2863基因编码339个氨基酸，命名为CtSDR3（GenBank登录号：MW792037）。系统进化树表明CtSDR1与蓟Cirsium japonicum (QQH14901.1)同源性最高；CtSDR2与小蓬草Erigeron canadensis (XP_043636506.1)同源性最高；CtSDR3与小豆蔻Cynara cardunculus var. scolymus (KVI09206.1)同源性最高。Prot-param分析CtSDR1基因所编码的蛋白质分子式C2230H3346N606O639S7，相对分子量为49.2×103，理论等电点pI=9.61；CtSDR2基因所编码的蛋白质分子式C1289H2072N360O379S13，相对分子量为29×103，理论等电点pI=8.63；CtSDR1基因所编码的蛋白质分子式C1691H2614N442O481S9，相对分子量为37.1×103，理论等电点pI=6.80。Prot Scale分析预测CtSDR1、CtSDR2和CtSDR3蛋白为亲水性蛋白，无信号肽属非分泌蛋白。蛋白跨膜性分析显示CtSDR1、CtSDR2和CtSDR3不含有跨膜区域，为非跨膜蛋白。对CtSDR1、CtSDR2和CtSDR3蛋白二级结构的预测显示均属于不规则结构。对CtSDR1、CtSDR2、CtSDR3蛋白质三维结构预测如图3所示。系统进化树如图4所示。亚细胞定位预测显示，CtSDR1、CtSDR2、CtSDR3均可能定位于细胞质。
图 3 CtSDR1、CtSDR2、CtSDR3蛋白三级结构预测
图 4 CtSDR系统发育进化树
-
取红花花期的Ⅳ期的红花各个部位进行分析，发现红花花冠内的CtSDR1、CtSDR2、CtSDR3基因表达量从高到低依次均为花冠>叶>茎>根。其中CtSDR1在花冠中的相对表达量约为根中的3倍、而CtSDR2、CtSDR3在花冠中的相对表达量约为根中的4倍。将4个花期的红花花冠进行qRT-PCR分析表明，CtSDR1、CtSDR2、CtSDR3花冠中表达量均随着花冠发育逐渐升高，特别是CtSDR1、CtSDR2、CtSDR3的Ⅳ期花冠对比Ⅲ期花冠的表达量分别提高了7.2倍、2.7倍、2.3倍（图5）。
图 5 目的基因在不同部位（A）和不同花期（B）的表达水平
-
构建真和表达载体并通过PCR鉴定后，我们从19株农杆菌浸染的子代植株中得到5株pMT39-CtSDR3阳性红花植株（图6）。通过qPCR对其CtSDR3基因转录水平进行测定，结果发现阳性红花植株中CtSDR3基因的转录水平得到显著增加，约为空白组株系的2～3倍，CtSDR3的在花冠部位的高表达也证明了研究成功获取CtSDR3过表达红花植株（图7）。通过UPLC-QTOF/MS技术测定阳性转基因红花株系组和空白对照组的目标化合物含量，包括7个红花花冠主要黄酮类化合物及苯丙烷类代谢途径上游关键物质苯丙氨酸（图8），分别为：野黄芩素（scutellarein）、Carthamin、HSYA、山柰酚（kaempferol）、山柰酚-3-O-β-D-葡萄糖苷（kaempferol-3-O-β-D-glucoside）、山柰酚-3-O-β-D-芸香糖苷（kaempferol-3-O-β-rutinoside）、芦丁（rutin）和苯丙氨酸（D-Phenylalanine）。由图8可知，与空白组相比，CtSDR3过表达株系相较于空白组野黄芩素提高了3.6%～9.8%，HSYA提高了7.1%～16.6%，以及苯丙氨酸含量提高了5.5%～15.7%，具有显著性升高。其他化合物含量则有无显著性变化趋势。通过对过表达株系与空白组的含量分析，我们认为CtSDR3基因过表达会引起红花中黄酮类物质的变化，尤其是HSYA含量升高显著。同时，苯丙氨酸代谢途径属于植物重要的次生代谢途径，过表达组引起苯丙氨酸含量的显著上升，上述指标性成分的变化也进一步说明CtSDR3对红花黄酮类化合物次生代谢途径具有一定的影响，但目前我们尚难以判断CtSDR3红花中影响次生代谢产物积累的明确途径。
图 6 真核表达载体构建及阳性鉴定电泳图
图 7 过表达植株CtSDR3的相对表达量
图 8 阳性植株黄酮类化合物含量测定
-
目的片段成功扩增，将目的条带进行胶回收、纯化。CtSDR1、CtSDR1、CtSDR1构建的pGEX-6p-1、pET-28a原核表达载体均有在大肠杆菌内表达，但是CtSDR1-pGEX-6p-1、 CtSDR2-pGEX-6p-1、CtSDR3-pGEX-6p-1、CtSDR1-pET-28、CtSDR2-pET-28a、CtSDR3-pET-28a表达的目的蛋白均形成包涵体，存在于沉淀中。无法进行下一步大量纯化实验，唯有CtSDR2-pGEX-6p-1诱导表达了存在于上清液的目的蛋白，明显可以在上清液中观察到分子量约为50 000的蛋白条带（图9）。
图 9 目的片段PCR产物电泳及表达蛋白电泳分析
-
越来越多的红花药理学相关研究表明，红花的主要药效物质包括查尔酮类、黄酮醇类等多种黄酮类化合物，其中，查尔酮类HSYA对脑缺血具有保护作用，并且还能抗脑血栓形成以及抗氧化等。研究HSYA的生物合成分途径，对于HSYA的工业化生产具有重要意义。
本研究借助生物学分子技术、结合代谢组分析测定，筛选出3个参与HSYA生物合成途径的关键短链脱氢还原酶基因CtSDR1、CtSDR2和CtSDR3，这3个基因序列具有高度保守性，在不同器官的表达模式均呈现出花冠>叶>茎>根的特点，而且在花冠中的表达量随花冠发育逐渐升高，表明其很有可能参与红花中HSYA等主要药用成分的积累。进一步研究发现，转CtSDR3过表达T2代阳性植株花冠中CtSDR3基因的转录水平增加了2~3倍，次生代谢物HSYA的含量提高了7.1%～16.6%（P<0.05），验证了我们对CtSDR3在红花体内参与黄酮类化合物生物合成功能的推测。本研究中，体外表达CtSDR3蛋白，得到目的蛋白条带，但由于包涵体等原因，蛋白表达和纯化条件仍需要进一步摸索。下一步，我们将对可能起黄酮类生物合成途径的关键SDRs进行深入的生物学特性特别是酶结合位点的研究，为更好地阐释SDRs的生物学功能、利用分子生物育种技术培育高HSYA含量的红花新品种奠定基础。
Characterization and function of short-chain dehydrogenases/reductases in hydroxysafflower yellow A biosynthesis pathway
-
摘要:
目的 探究参与红花黄酮类生物合成途径特别是羟基红花黄色素A（HSYA）关键短链脱氢还原酶基因（SDRs）的功能。 方法 基于红花转录组数据库以及代谢组数据库，筛选参与HSYA生物合成途径的SDRs，用qRT-PCR法分析表达模式。采用无缝克隆技术构建过表达载体，以农杆菌GV3101介导遗传转化云南巍山红花品系，对转基因T2代植株进行阳性验证，并对花冠SDRs的基因表达量进行分析，UPLC-Q-TOF/MS法测定次生代谢物的含量。 结果 筛选出3个参与HSYA生物合成途径的关键短链脱氢还原酶基因CtSDR1、CtSDR2、CtSDR3，表达量从高到低依次为花冠>叶>茎>根。花冠中表达量随花冠发育逐渐升高。对转基因T2代植株进行阳性验证后的花冠进行qRT-PCR分析发现：与空白对照组相比，转CtSDR3过表达T2代阳性植株花冠中CtSDR3基因的转录水平增加了2~3倍，次生代谢物HSYA的含量提高了7.1%～16.6%（P＜0.05）。 结论 CtSDR3可能参与了红花中黄酮类化合物特别是HSYA的生物合成，为阐释CtSDR3在HSYA生物合成途径中的功能提供了数据支撑。 Abstract:Objective To explore the function of short-chain dehydrogenases/reductases (SDRs) in safflower flavonoid, especially hydroxysafflower yellow A (HSYA) biosynthesis. Methods SDRs involved in HSYA biosynthesis pathway were screened based on safflower transcriptome database and metabolome database. The expression pattern was analyzed by qRT-PCR. The overexpression vector was constructed by seamless cloning technology, then genetically transformed to the Yunnan Weishan safflower strain by Agrobacterium gv3101. The transgenic T2 generation plants were positively verified, and the gene expression of corolla SDRs was analyzed. The content of secondary metabolites was assayed by UPLC-Q-TOF/MS. Results Three SDRs genes named CtSDR1, CtSDR2 and CtSDR3 involved in HSYA biosynthesis pathway were screened. Their expression in safflower from high to low was corolla > leaf > stem > root. The expression level in corolla increased gradually with corolla development. qRT-PCR analysis of corolla with positive verification of genome insertion sequence showed that the transcription level of CtSDR3 in corolla of T2 positive plants increased by 2~3 times compared with the blank control group, and the content of secondary metabolite HSYA increased by 7.1%~16.6% (P< 0.05). Conclusion CtSDR3 may be involved in the biosynthesis of flavonoids, especially HSYA, in safflower. It provides the support data for explaining the function of CtSDR3 in HSYA biosynthesis pathway. -
Key words:
- short-chain dehydrogenase reductase /
- hydroxysafflor yellow A /
- safflower
-
骨质疏松性椎体压缩骨折（OVCF）属于脊柱骨折，多发于中老年人，为临床最为常见的骨质疏松性骨折之一，严重影响患者脊柱功能、日常活动能力等，且发病率、致残、致死率均较高。现阶段，临床多以经皮椎体成形术（PVP）对OVCF患者进行治疗，临床疗效确切，可迅速缓解骨折引起的疼痛，促进椎体形态和功能恢复，但其对患者骨质疏松程度及骨质疏松引起的疼痛改善效果欠佳，故临床倡导PVP术后予以适当的干预以改善手术效果[1]。西医临床常用药物为碳酸钙D3，但有研究[2]指出其对部分患者的干预效果较差，临床多与中医药物干预、锻炼等相结合对PVP术后OVCF患者进行治疗。中医学认为，OVCF属“骨痿”范畴，主要病机为肾精不足、气滞血瘀所致骨质枯槁，肝肾阴虚证为其主要证型之一，应治以滋补肝肾，填精壮骨[3]。中药膏方含阿胶、枸杞子、鸡血藤、党参等，根据临床中药机制可见其有补肾填精、益气健脾、通络止痛等功效，但其对PVP术后OVCF患者的疗效及机制尚未完全明确[4]。故设立本研究，以120例于我院行PVP治疗的OVCF患者为研究对象进行前瞻性随机对照研究，旨在进一步观察中药膏方联合八段锦治疗PVP术后OVCF患者的疗效，为其临床应用提供参考。
1. 资料与方法
1.1 一般资料
前瞻性选取120例2016年1月至2017年9月于我院行PVP治疗的OVCF患者为研究对象，患者均知情同意，按随机数字表法将其分为观察组（60例），对照组（60例）。本研究经我院医学研究伦理研究委员会审核通过。观察组：男35例，女25例；年龄50~80岁，平均（66.27±2.15）岁；体质量指数（BMI）19~25 kg/m2，平均（22.35±0.41）kg/m2；骨质疏松病程1~5年，平均（2.51±0.22）年；骨折病程1~3周，平均（1.55±0.12）周；骨压缩程度[5]：轻度16例，中度30例，重度14例。对照组：男33例，女27例；年龄51~78岁，平均（66.18±2.20）岁；BMI 19~25 kg/m2，平均（22.41±0.42）kg/m2；骨质疏松病程1~5年，平均（2.48±0.20）年；骨折病程1~3周，平均（1.52±0.13）周；骨压缩程度：轻度15例，中度29例，重度16例。两组性别、年龄、BMI、骨质疏松病程、骨折病程、压缩程度等一般资料比较差异无统计学意义（P>0.05），具有可比性。
1.2 纳入标准
①中医诊断符合《中医药防治原发性骨质疏松症专家共识（2015）》[6]中肝肾阴虚证相关辩证标准，西医诊断符合《骨质疏松性骨折诊疗指南》[7]中OVCF相关诊断标准，并经临床检查确诊者；②均为胸腰椎骨折，且是首次因OVCF就诊者；③无PVP相关禁忌，手术成功，且术后生命体征平稳者；④近3个月内未行特殊的抗骨质疏松治疗，同期未接受其他治疗者；⑤无脊髓损伤或神经根症状者等。
1.3 排除标准
①椎体爆裂性骨折、陈旧性胸腰椎压缩性骨折或肿瘤、感染等所致病理骨折者；②脊柱结核、骨肿瘤者；③继发性骨质疏松者；④精神障碍、认知能力或沟通能力障碍者；⑤伴严重器官障碍、急性腰扭伤、骨瘤骨转移、强直型脊柱炎、四肢新鲜骨折、软组织损伤者等。
1.4 方法
对照组予以碳酸钙D3咀嚼片[重庆海默尼制药有限公司，国药准字H20205039，碳酸钙1.25 g（相当于钙0.50 g）、维生素D3 200 IU]口服，每次1片，每日2次。同时指导患者每天摄入适量豆类、瘦肉、牛奶、鱼类等优质蛋白（1.20~1.40 g/kg），进行适当的日光浴及行走锻炼（30 min以上）。
观察组在对照组基础上予以中药膏方（将阿胶、枸杞子、鸡血藤、党参各200 g，鹿角胶、茯苓、桑寄生、熟地黄、牛膝、丹参、威灵仙各150 g，龟板胶、白术、骨碎补、神曲各100 g，炙甘草60 g制成膏剂后即得）口服，开水调服，每次10 g，每日2次。同时指导患者进行八段锦锻炼，所用方法为改良八段锦第八式：躯体直立，两足平行分开（与肩同宽），两手臂自然下垂，手指稍并拢，掌指向前，两膝关节屈曲135°~170°；两眼平视前方，微张口，平稳呼吸，全身放松2 min，紧接着保持原来体位尽可能屈膝下蹲，然后在双上臂伸直情况下慢慢上举到水平，慢慢随之站起（保持躯体直立下），至身体直立时顺势将两脚跟向上提起（期间慢慢匀速吸气到最满并憋气），停顿10~15 s；两脚跟下落着地，双手臂随之顺势下落，回到膝关节屈曲位（期间慢慢匀速呼气）；每次7个循环，每日3次。两组均治疗6个月，并随访3年。
1.5 观察指标
1.5.1 疗效
根据《中药新药临床研究指导原则》[8]可将两组治疗6个月后的疗效分为临床治愈（腰部功能完全恢复，临床症状、体征完全消失，X射线可见压缩椎体恢复正常状态），显效（腰部功能没有完全恢复，临床症状、体征基本消失，X射线可见压缩椎体恢复正常状态），有效（腰部功能没有完全恢复，临床症状、体征减轻，X射线可见大量骨痂形成，但未完全愈合），无效（腰部功能有改善，临床症状、体征有减轻，但X射线可见几乎无骨痂形成）。临床治愈率=临床治愈例数/总例数；显效率=显效例数/总例数；有效率=有效例数/总例数；无效率=无效例数/总例数。总有效率=（临床治愈+显效+有效）例数/总例数。
1.5.2 不同时点腰背疼痛情况
以疼痛视觉模拟评分法（VAS，0~10分）[9]评估两组治疗1、3、6个月的腰背疼痛情况，得分越高提示患者腰背疼痛程度越严重。
1.5.3 骨密度（BMD）、椎体后凸角度（Cobb角）及椎体前壁高度（AVBH）变化情况
以BMD测定仪（QDR 4500C，美国Hologic公司）检测两组治疗前、治疗6个月后的腰椎、股骨颈BMD；以X射线检查两组治疗前、治疗6个月后的Cobb角、AVBH。
1.5.4 骨代谢指标
抽取两组治疗前、治疗6个月后的晨起空腹静脉血3 ml，离心机3 000 r/min离心15 min，分离血清，以全自动电化学发光分析仪（罗氏Cobas e601，北京东方迈润医疗器械有限公司）检测血清I型原胶原降解产物（β-Cross I）、N端中段骨钙素（N-MID Ost）、甲状旁腺素（PTH）水平。
1.5.5 PVP术后推体再骨折发生率
记录两组随访时间及随访1年内、1~3年内PVP术后推体再骨折发生率，同1例患者发生多次再骨折仅算1例。
1.6 统计学方法
数据采用SPSS 21.0（IBM公司，美国）统计学软件进行分析，计量资料采用（
$\bar x $ ±s）表示，组内比较使用配对t检验，治疗前、治疗6个月后的比较使用独立样本t检验，多组间比较采用F检验；计数资料采用[n(%)]表示，使用χ2检验进行比较。2. 结果
2.1 两组患者临床疗效比较
治疗6个月后，观察组临床治愈率为73.33%，高于对照组的53.33%（P<0.05）；而两组显效率、有效率、无效率、总有效率比较差异无统计学意义（P>0.05，表1）。
表 1 两组患者临床疗效比较[n(%)]组别 例数 临床治愈率 显效率 有效率 无效率 总有效率 观察组 60 44(73.33) 13(21.67) 3(5.00) 0(0.00) 60(100.00) 对照组 60 32(53.33) 20(33.33) 8(13.33) 0(0.00) 60(100.00) χ2 5.167 2.048 2.502 P <0.05 >0.05 >0.05 2.2 两组患者不同时点腰背疼痛情况比较
治疗前，两组患者VAS评分比较差异无统计学意义（P>0.05）。与治疗前比较，治疗3、6个月后两组VAS评分均降低，且观察组低于对照组（P<0.05）；与治疗3个月后比较，治疗6个月后，两组VAS评分均降低（P<0.05，表2）。
表 2 两组患者不同时点腰背疼痛情况比较（$\bar x $ ±s，分）组别 例数 治疗前 治疗3个月后 治疗6个月后 观察组 60 7.55±0.37 2.47±0.33* 1.14±0.23** 对照组 60 7.52±0.35 3.51±0.51* 2.36±0.42** t 0.456 13.262 19.735 P >0.05 <0.05 <0.05 *P<0.05，与治疗前比较；**P<0.05，与治疗3个月后比较 2.3 两组患者BMD、Cobb角及AVBH变化情况比较
治疗前，两组患者BMD、Cobb角及AVBH比较差异无统计学意义（P>0.05）。与治疗前比较，治疗6个月后两组腰椎、股骨颈BMD及AVBH均升高，且观察组高于对照组（P<0.05）；而两组Cobb角均降低，且观察组低于对照组（P<0.05，表3）。
表 3 两组患者BMD、Cobb角及AVBH变化情况比较（$\bar x $ ±s）时间 组别 例数 腰椎BMD(g/cm3) 股骨颈BMD(g/cm3) Cobb角(°) AVBH(%) 治疗前 观察组 60 0.70±0.07 0.68±0.06 23.09±2.51 53.93±7.26 对照组 60 0.69±0.08 0.67±0.07 23.11±2.49 54.04±7.44 t 0.729 0.840 0.044 0.082 治疗6个月后 观察组 60 0.93±0.08* 0.93±0.07* 6.47±1.98* 96.56±9.41* 对照组 60 0.84±0.11* 0.85±0.08* 8.91±2.26* 85.73±8.04* t 5.125 5.829 6.290 6.778 *P<0.05，与治疗前比较 2.4 两组患者骨代谢指标比较
治疗前，两组患者血清β-Cross I、N-MID Ost、PTH水平比较差异无统计学意义（P>0.05）。与治疗前比较，治疗6个月后，两组血清β-Cross I、N-MID Ost、PTH水平均降低，且观察组低于对照组（P<0.05，表4）。
表 4 两组患者骨代谢指标比较（$\bar x $ ±s）时间 组别 例数 β-Cross I(ng/ml) N-MID Ost(ng/ml) PTH(pg/ml) 治疗前 观察组 60 0.79±0.13 42.10±9.53 62.28±9.42 对照组 60 0.78±0.14 41.98±9.52 62.25±9.39 t 0.405 0.069 0.017 治疗6个月后 观察组 60 0.67±0.10* 36.15±5.46* 52.47±4.02* 对照组 60 0.73±0.12* 38.93±6.14* 57.45±5.72* t 2.975 2.621 5.518 *P<0.05，与治疗前比较 2.5 两组患者PVP术后推体再骨折发生率比较
两组患者随访时间、随访1年内、1~3年内PVP术后推体再骨折发生率比较，均无统计学意义（P>0.05）；随访3年内，观察组PVP术后推体再骨折发生率为3.33%，低于对照组的20.00%（P<0.05，表5）。
表 5 两组患者PVP术后推体再骨折发生率比较组别 例数 随访1年内发生率
[n(%)]随访1~3年内 合计
[n(%)]随访时间
($\bar x $±s, 年)发生率
[n(%)]观察组 60 2(3.33) 1.94±0.43 0(0.00) 2(3.33) 对照组 60 7(11.67) 1.86±0.72 5(8.33) 12(20.00) t/χ2 3.003 0.739 3.339 8.086 P >0.05 >0.05 >0.05 <0.05 3. 讨论
骨质疏松为临床常见的慢性骨代谢异常疾病，可见骨钙含量减少、骨量下降，极易出现骨折，OVCF为其临床常见类型。OVCF发病后可致患者出现椎体高度下降、疼痛、功能障碍等，未得到及时有效的治疗，可致残、致死，且其再发骨折风险高[10]。现阶段，微创PVP治疗为大多保守治疗无效的OVCF患者的主要选择，可恢复椎体力学强度、稳定骨折，但骨转换率及骨量流失未能纠正，且手术造成的内外骨膜血管损伤及制动可能会加重骨质疏松而影响患者预后[11]。故PVP术后予以适当的干预治疗至关重要，临床常用药物为碳酸钙D3，其可促进骨质形成而缓解骨质疏松程度，但其多偏重于抑制骨吸收，减轻骨质流失。近年来，越来越多的研究[12-13]发现，中医药治疗PVP术后OVCF有一定成效，其既能抑制流失，又能促进成骨形成，在抗骨质疏松、促进骨折愈合、防治再骨折发生等方面优势明显，现已备受临床医生青睐，但关于何种药物更佳尚无统一定论。
中医认为OVCF属“骨枯”、“骨萎”、“痹证”等范畴，肾中精气是骨生长发育的根本，脾胃为后天之本，气血生化之源，故其发病机制主要为肾精亏虚、脾精不足所致骨骼失养；外伤导致筋脉损伤、气血运行不畅，加之患者湿热寒邪侵袭、肾脾虚亏，容易加重气滞血瘀症状，使得骨骼和筋骨失去濡养，导致骨关节萎缩、活动不利等而诱发OVCF，故治疗应补肾壮骨、活血行气、舒经通络、消肿止痛[14]。中药膏方中的阿胶、枸杞子、鹿角胶、龟板胶可补血益气、填精益髓，其中，龟鹿二仙胶中鹿角胶、龟板胶为主方，与党参、枸杞子配伍可增强益气养阴、补肾填精益髓之功；鸡血藤可活血补血、舒筋活络、疏风止痛；党参可补中益气、生血行血；茯苓可利水消肿、益气健脾；桑寄生、熟地黄、骨碎补、牛膝可补肾壮骨；丹参、威灵仙、神曲可健脾和胃、消食调中，进而减轻膏方的厚腻；炙甘草可补脾和胃，益气复脉，同时调和诸药；全方可发挥补肾填精、益气健脾、通络止痛之功效[15-16]。同时八段锦作为一种由八节不同动作组成的一套医疗、康复体操，其通过人体垂直方向的应力刺激成骨、下蹲平衡锻炼、膀胱经的经气疏通等可对OVCF患者产生舒筋活络的良好功效，进而有助于促进患者腰部功能恢复，提高患者临床疗效。本研究显示，治疗6个月后，观察组临床治愈率高于对照组，且治疗3、6个月后VAS评分低于对照组，进一步说明中药膏方联合八段锦治疗PVP术后OVCF可有效缓解患者腰背疼痛，疗效显著。
OVCF的发生发展与骨代谢及转换密切相关，其中，β-Cross I为骨吸收的重要指标，N-MID Ost为骨形成、骨转换的重要指标，可直接反映成骨细胞的活性与数量变化，PTH则可加速骨代谢活跃，促进骨吸收，血清β-Cross I、N-MID Ost、PTH水平升高提示患者骨转换成程度严重，可促进病情进展，不利于PVP术后腰部功能的恢复[17]。本研究显示，治疗6个月后，观察组腰椎、股骨颈BMD及AVBH高于对照组，Cobb角、血清β-Cross I、N-MID Ost、PTH低于对照组，且随访3年内PVP术后推体再骨折发生率低于对照组，提示中药膏方联合八段锦治疗PVP术后OVCF可调节患者机体骨代谢，提高腰椎、股骨颈BMD及AVBH，降低Cobb角，促进腰部功能的恢复，降低PVP术后推体再骨折发生率。PVP术后OVCF患者骨折愈合的本质为骨重建，该过程涉及局部微环境的构建、成骨细胞的活化、大量血管的生成等，而现代药理学研究[18-20]指出，中药膏方中的骨碎补含双氢黄酮苷、黄酮等成分，可提高BMD，刺激成骨细胞增殖分化，防止骨吸收、促进骨形成，发挥抗骨质疏松、促进骨折愈合等作用，进而促进腰部功能的恢复；牛膝含牛膝竹节参皂苷，可发挥较好的抗炎、止痛、改善微循环作用，有助于促进局部微环境的构建而促进患者骨折愈合，进而降低再骨折发生率。同时联合八段锦进行治疗，可有效疏通OVCF患者经络，促进机体气血运行通畅，进而可增强中药膏方促进患者腰部功能恢复的作用，降低患者PVP术后推体再骨折发生率。而对照组由于仅进行常规对症治疗，而未对患者采取中药膏方联合八段锦等治疗，故患者复发的风险明显增加。
综上，中药膏方联合八段锦治疗PVP术后OVCF可降低其血清β-Cross I、N-MID Ost、PTH水平，调节机体骨代谢，有助于提高腰椎、股骨颈BMD及AVBH，降低Cobb角，促进腰部功能的恢复，进而缓解患者腰背疼痛，降低PVP术后推体再骨折发生率，疗效显著。
-
图 6 真核表达载体构建及阳性鉴定电泳图
注：1. CtSDR3基因开放阅读框（ORF）区扩增产物电泳图，a、b泳道均为CtSDR3基因ORF区克隆PCR产物；2. 真核表达载体pMT-39载体酶切产物电泳图，a、b泳道为CtSDR3 PCR产物，c泳道为pMT-39载体，d、e泳道为pMT-39线性化载体；3. pMT39-CtSDR3重组载体阳性转化子鉴定电泳图，a、b泳道为阳性转化子菌液PCR产物；4. pMT39-CtSDR3质粒转化农杆菌GV3101，a、c和e泳道为空白对照组，b、d和f泳道为阳性克隆菌液PCR产物；5. 红花pMT39-CtSDR3阳性转化植株鉴定PCR产物电泳图，1～19为待鉴定植株，p为pMT39-CtSDR3质粒，k为空白组，WT为野生型红花植株
图 8 阳性植株黄酮类化合物含量测定
注：A.黄芩素；B.Carthamin；C.HSYA；D.山柰酚；E.山柰酚-3-O-葡萄糖苷；F.山柰酚-3-O-芸香糖苷；G.芦丁；H.苯丙氨酸；CK.空白组株系；OVX.阳性过表达株系
-
[1] 郭美丽, 张汉明, 张美玉. 红花本草考证[J]. 中药材, 1996(4):202. [2] TU Y H, XUE Y R, GUO D D, et al. Carthami Flos: a review of its ethnopharmacology, pharmacology and clinical applications[J]. REV BRAS FARMACOGN,2015,25(5):553-566. doi: 10.1016/j.bjp.2015.06.001 [3] WANG Y Q, ZHANG S S, NI H L, et al. Autophagy is involved in the neuroprotective effect of nicotiflorin[J]. J Ethnopharmacol,2021,278:114279. doi: 10.1016/j.jep.2021.114279 [4] 贾鑫磊, 何贝轩, 郭丹丹, 等. 红花花冠伸长相关基因CtXTH1的特征与功能研究[J]. 药学学报, 2019, 54(6):1132-1140. [5] 张戈. 红花的化学成分及品质评价[D]. 上海: 第二军医大学, 2001. [6] 张戈, 郭美丽, 张汉明, 等. 不同种质红花药材的高效液相色谱法指纹图谱研究[J]. 第二军医大学学报, 2006, 27(3):280-283. doi: 10.3321/j.issn:0258-879X.2006.03.012 [7] 丁丽丽, 段陈平, 李芳, 等. 红花不同采收期及不同部位中羟基红花黄色素A及山奈素的含量变化[J]. 沈阳药科大学学报, 2015, 32(1):65-69. [8] 杨捷威, 吴婷婷, 郭美丽. 红花组织培养的研究进展[J]. 药学服务与研究, 2012, 12(1):58-62. [9] HUANG L L, YANG X, SUN P, et al. The first Illumina-based de novo transcriptome sequencing and analysis of safflower flowers[J]. PLoS One,2012,7(6):e38653. doi: 10.1371/journal.pone.0038653 [10] COSGROVE D J. New genes and new biological roles for expansins[J]. Curr Opin Plant Biol,2000,3(1):73-78. doi: 10.1016/S1369-5266(99)00039-4 [11] CHENG H, LI L L, XU F, et al. Expression patterns of an isoflavone reductase-like gene and its possible roles in secondary metabolism in Ginkgo biloba[J]. Plant Cell Rep,2013,32(5):637-650. doi: 10.1007/s00299-013-1397-2 [12] MIN T, KASAHARA H, BEDGAR D L, et al. Crystal structures of pinoresinol-lariciresinol and phenylcoumaran benzylic ether reductases and their relationship to isoflavone reductases[J]. J Biol Chem,2003,278(50):50714-50723. doi: 10.1074/jbc.M308493200 [13] 崔扬, 冯彦辉, 陈众峰, 等. 玉米短链脱氢酶基因IDP2557的克隆及其抗旱功能鉴定[J]. 分子植物育种, 2019, 17(22):7300-7305. [14] 王化冰. 黄瓜短链脱氢酶SDR110C家族的生物信息学分析及CsSDR110C14的遗传转化[D]. 泰安: 山东农业大学, 2014. [15] BRAY J E, MARSDEN B D, OPPERMANN U. The human short-chain dehydrogenase/reductase (SDR) superfamily: a bioinformatics summary[J]. Chem Biol Interact,2009,178(1-3):99-109. doi: 10.1016/j.cbi.2008.10.058 [16] JÖRNVALL H, HEDLUND J, BERGMAN T, et al. Superfamilies SDR and MDR: from early ancestry to present forms. Emergence of three lines, a Zn-metalloenzyme, and distinct variabilities[J]. Biochem Bioph Res Commun,2010,396(1):125-130. doi: 10.1016/j.bbrc.2010.03.094 [17] KALLBERG Y, OPPERMANN U, PERSSON B. Classification of the short-chain dehydrogenase/reductase superfamily using hidden Markov models[J]. FEBS Journal,2010,277(10):2375-2386. doi: 10.1111/j.1742-4658.2010.07656.x [18] OPPERMANN U, FILLING C, HULT M, et al. Short-chain dehydrogenases/reductases (SDR): the 2002 update[J]. Chemico-Biological Interactions,2003,143-144:247-253. doi: 10.1016/S0009-2797(02)00164-3 [19] WU Y, WANG Y, MI X F, et al. The QTL GNP1 encodes GA20ox1, which increases grain number and yield by increasing cytokinin activity in rice panicle meristems[J]. PLoS Genet,2016,12(10):e1006386. doi: 10.1371/journal.pgen.1006386 [20] GUO T, CHEN K, DONG N Q, et al. GRAIN SIZE AND NUMBER1 negatively regulates the OsMKKK10-OsMKK4-OsMPK6 cascade to coordinate the trade-off between grain number per panicle and grain size in rice[J]. Plant Cell,2018,30(4):871-888. doi: 10.1105/tpc.17.00959 [21] 靳甜甜. 植物次生代谢途径中黄酮类化合物的生物合成研究[D]. 保定: 河北农业大学, 2014. [22] LEE J, BURNS T H, LIGHT G, et al. Xyloglucan endotransglycosylase/hydrolase genes in cotton and their role in fiber elongation[J]. Planta,2010,232(5):1191-1205. doi: 10.1007/s00425-010-1246-2 [23] SHAO M Y, WANG X D, NI M, et al. Regulation of cotton fiber elongation by xyloglucan endotransglycosylase/hydrolase genes[J]. Genet Mol Res,2011,10(4):3771-3782. doi: 10.4238/2011.October.27.1 [24] 谭政委, 鲁丹丹, 李磊, 等. 红花二氢黄酮醇4-还原酶基因CtDFR1的克隆与分析[J]. 分子植物育种, 2021:1-12. [25] XU F, DENG G, CHENG S, et al. Molecular cloning, characterization and expression of the phenylalanine ammonia-lyase gene from Juglans regia[J]. Molecules,2012,17(7):7810-7823. doi: 10.3390/molecules17077810 [26] SHOJI T, WINZ R, IWASE T, et al. Expression patterns of two tobacco isoflavone reductase-like genes and their possible roles in secondary metabolism in tobacco[J]. Plant Mol Biol,2002,50(3):427-440. doi: 10.1023/A:1019867732278 [27] 朱惠, 杨丽涛, 李杨瑞. 基于三代测序技术的转录组学研究[J]. 生物技术报, 2014(11):130-137. [28] PETIT P, GRANIER T, D'ESTAINTOT B L, et al. Crystal structure of grape dihydroflavonol 4-reductase, a key enzyme in flavonoid biosynthesis[J]. J Mol Biol,2007,368(5):1345-1357. doi: 10.1016/j.jmb.2007.02.088 [29] O'REILLY C, SHEPHERD N S, PEREIRA A, et al. Molecular cloning of the a1 locus of Zea mays using the transposable elements En and Mu1[J]. EMBO J,1985,4(4):877-882. doi: 10.1002/j.1460-2075.1985.tb03713.x [30] 郝爱平. 植物二氢黄酮醇4-还原酶的生物信息学分析[J]. 江苏农业科学, 2014, 42(6):30-34. doi: 10.3969/j.issn.1002-1302.2014.06.009 [31] ISHIDA T, HATTORI S, SANO R, et al. Arabidopsis transparent testa glabra2 is directly regulated by r2r3 myb transcription factors and is involved in regulation of glabra2 transcription in epidermal differentiation[J]. Plant Cell,2007,19(8):2531-2543. doi: 10.1105/tpc.107.052274 [32] ZHANG F, GONZALEZ A, ZHAO M Z, et al. A network of redundant bHLH proteins functions in all TTG1-dependent pathways of Arabidopsis[J]. Development,2003,130(20):4859-4869. doi: 10.1242/dev.00681 [33] 郭丹丹. 红花黄酮类成分生物合成途径关键基因的特征和功能研究[D]. 上海: 海军军医大学, 2019. [34] GRÄFF M, BUCHHOLZ P C F, STOCKINGER P, et al. The Short-chain Dehydrogenase/Reductase Engineering Database (SDRED): a classification and analysis system for a highly diverse enzyme family[J]. Proteins,2019,87(6):443-451. doi: 10.1002/prot.25666 [35] STAVRINIDES A K, TATSIS E C, DANG T T, et al. Discovery of a short-chain dehydrogenase from Catharanthus roseus that produces a new monoterpene indole alkaloid[J]. Chembiochem,2018,19(9):940-948. doi: 10.1002/cbic.201700621 [36] ZHOU Y, ZHANG L, GUI J D, et al. Molecular cloning and characterization of a short-chain dehydrogenase showing activity with volatile compounds isolated from Camellia sinensis[J]. Plant Molecular Biology Reporter,2015,33(2):253-263. doi: 10.1007/s11105-014-0751-z -
点击查看大图
计量
- 文章访问数: 7092
- HTML全文浏览量: 2473
- PDF下载量: 58
- 被引次数: 0
