邢欣 趙大慶 王思明 于士婷 李懿軒 劉美辰

中圖分類號 R285;R966 文獻標志碼 A 文章編號 1001-0408（2021）23-2859-10

DOI 10.6039/j.issn.1001-0408.2021.23.08

摘 要 目的：研究山藥蛋白（DOT）防治糖尿病性勃起功能障礙（DIED）的潛在作用機制。方法：采用高糖高脂飼料+腹腔注射鏈脲佐菌素（40 mg/kg）的方法復(fù)制DIED大鼠模型。實驗設(shè)置正常對照組（生理鹽水），模型組（生理鹽水），DOT低、中、高劑量組（0.3、0.6、0.9 mg/kg）和西地那非組（陽性對照，4.4 mg/kg），每組9只大鼠。在糖尿病模型已建立成功、開始誘導(dǎo)DIED發(fā)生的階段，各組大鼠開始灌胃相應(yīng)溶液，每天1次，連續(xù)11周。測定大鼠體質(zhì)量、空腹血糖（FPG）值、勃起次數(shù)、勃起率、空腹胰島素（FINS）、胰島素抵抗指數(shù)（IR）以及陰莖海綿體組織中內(nèi)皮型一氧化氮合酶（eNOS）、環(huán)磷酸鳥苷（cGMP）含量，以評價DOT對DIED模型大鼠的干預(yù)作用。選擇30 mmol/L葡萄糖誘導(dǎo)48 h建立高糖損傷小鼠海綿體內(nèi)皮細胞（MCECs）模型，并給予125、250、500 μg/mL DOT進行干預(yù)，然后檢測細胞活力，以評價DOT對高糖損傷MCECs的影響。繼而采用轉(zhuǎn)錄組測序技術(shù)分別篩選正常MCECs與高糖損傷MCECs、高糖損傷MCECs與250 μg/mL DOT處理的MCECs中的差異表達基因，對差異表達基因進行基因本體（GO）功能富集分析和京都基因與基因組百科全書（KEGG）通路富集分析;分析兩組細胞間共同的差異表達基因，并選擇6個關(guān)鍵基因在mRNA水平上進行驗證。結(jié)果：不同劑量DOT對DIED模型大鼠體質(zhì)量降低、FINS和IR升高、勃起次數(shù)減少、勃起率降低以及陰莖海綿體組織中eNOS、cGMP含量降低均有不同程度的逆轉(zhuǎn)作用，其中高劑量DOT作用后DIED模型大鼠的上述指標變化均被顯著逆轉(zhuǎn)（P<0.05或P<0.01）。125、250、500 μg/mL DOT均可顯著提高高糖損傷MCECs的活力（P<0.05或P<0.01）。轉(zhuǎn)錄組測序結(jié)果表明，與正常MCECs比較，高糖損傷MCECs中共有48個差異表達基因;與高糖損傷MCECs比較，DOT處理的MCECs中共有779個差異表達基因;兩組差異表達基因在生物學(xué)過程注釋中均以細胞過程為主、在細胞組分注釋中均以細胞部分為主、在分子功能注釋中均以綁定分子功能為主，并主要富集在細胞外基質(zhì)受體相互作用通路、錯配修復(fù)通路、磷脂酰肌醇3激酶-蛋白激酶B（PI3K-Akt）信號通路等。在兩組差異表達基因中，共得到Aldh1a1、Abcc5、Tac1等13個共同的差異表達基因，DOT均可顯著回調(diào)高糖損傷MCECs中上述共同差異表達基因的表達。經(jīng)驗證，DOT可以顯著逆轉(zhuǎn)高糖損傷MCECs中關(guān)鍵基因TGF-β3、Txnip、Aldh1a1、Loxl1、Mt1、Mt2的mRNA表達。結(jié)論：DOT對DIED模型大鼠的癥狀具有改善作用;其作用機制可能與抑制纖維化、降低氧化應(yīng)激等生物學(xué)途徑，從而改善海綿體內(nèi)皮功能有關(guān)。

關(guān)鍵詞 山藥蛋白;糖尿病性勃起功能障礙;轉(zhuǎn)錄組測序技術(shù);纖維化;氧化應(yīng)激;差異表達基因

Analysis of the Mechanism of Yam Protein in the Prevention and Treatment of Diabetes-induced Erectile Dysfunction Based on RNA-Seq Technology

XING Xin，ZHAO Daqing，WANG Siming，YU Shiting，LI Yixuan，LIU Meichen（Jilin Institute for Ginseng Science Research， Changchun University of Chinese Medicine， Changchun 130117， China）

ABSTRACT? ?OBJECTIVE： To study the potential mechanism of yam protein （DOT） in the prevention and treatment of diabetes-induced erectile dysfunction （DIED）. METHODS： DIED model was induced by high-glucose and high-fat diet and intraperitoneal injection of streptozotocin （40 mg/kg）. The experiment was set up in the normal control group （normal saline）， model group （normal saline）， DOT low-dose， medium-dose and high-dose groups （0.3， 0.6， 0.9 mg/kg）， sildenafil group （positive control， 4.4 mg/kg）， with 9 rats in each group. In the stage of successful establishment of diabetes model and initiation of inducing DIED， rats in each group were given relevant solution intragastrically， once a day， for consecutive 11 weeks. Body weight， fasting plasma glucose （FPG）， the times and rate of penile erection， fasting insulin （FINS）， insulin resistance index （IR）， the contents of endothelial nitric oxide synthase （eNOS） and cyclic guanosine monophosphate （cGMP） in penile cavernous tissue were determined so as to evaluate the intervention effects of DOT on DIED model rats. High-glucose damaged mice cavernous endothelial cells （MCECs） model was induced by 30 mmol/L glucose for 48 h， and then give DOT 125， 250， 500 μg/mL. The cell viability was detected so as to evaluate the effects of DOT on high-glucose damaged MCECs model. RNA-Seq technology was adopted to screen the differentially expressed genes between normal MCECs and high-glucose damaged MCECs， high-glucose damaged MCECs and MCECs treated with 250 μg/mL DOT. Gene ontology （GO） function enrichment analysis and KEGG pathway enrichment analysis were performed for differentially expressed genes. The common differentially expressed genes between 2 groups were analyzed， and mRNA expressions of six key genes were validated. RESULTS： Different doses of DOT could reverse the reduction of body weight， the increase of FINS and IR， the reduction of the times and rate of penile erection， the decrease of eNOS and cGMP contents in penile cavernous tissue of DIED model rats; above indexes of DIED model rats were reversed significantly after treated with high-dose of DOT （P<0.05 or P<0.01）. 125， 250， 500 μg/L DOT could significantly improve the activity of high-glucose damaged MCECs （P<0.05 or P<0.01）. RNA-Seq technology showed that compared with normal MCECs， a total of 48 differentially expressed genes were found in high-glucose damaged MCECs. Compared with high-glucose damaged MCECs， a total of 779 differentially expressed genes were found in MCECs treated with DOT. The differentially expressed genes of 2 groups were mainly cellular process in biological process annotation， cellular part in cell component annotation and binding molecular function in molecular function annotation， which were mainly enriched in extracellular matrix receptor interaction pathway， mismatch repair pathway， phosphatidylinositol 3 kinase-protein kinase B （PI3K-Akt） signal pathway and so on. Among differentially expressed genes of 2 groups， 13 common differentially expressed genes such as Aldh1a1， Abcc5， Tac1 were found. DOT could significantly reverse the expression of the above common differentially expressed genes in high-glucose damaged MCECs. After validation， DOT could significantly reverse the mRNA expression of TGF-β3， Txnip， Aldh1a1， Loxl1， Mt1 and Mt2 in high-glucose damaged MCECs. CONCLUSIONS： DOT could improve the symptom of DIED model rats， the mechanism of which may be related to biological pathway of inhibiting fibrosis and reducing oxidative stress， so as to improve the endothelial function of cavernous body.

KEYWORDS? ?Yam protein; Diabetes-induced erectile dysfunction; RNA-Seq technology; Fibrosis; Oxidative stress; Differentially expressed genes

糖尿病性勃起功能障礙（diabetic erectile dysfunction，DIED）是糖尿病的常見并發(fā)癥，發(fā)病率高達50%～75%，嚴重影響了患者的生活質(zhì)量[1-3]。高血糖引起的機體氧化應(yīng)激、纖維化是致使內(nèi)皮功能障礙，進而誘發(fā)DIED的重要原因[4]。然而，由于內(nèi)皮功能障礙患者的內(nèi)源性一氧化氮（NO）生成率較低，致使治療勃起功能障礙的“黃金藥物”——5型磷酸二酯酶抑制劑對DIED的療效較差[5-7]。因此，研發(fā)可有效防治DIED的新型藥物具有重要意義。

山藥是薯蕷科植物薯蕷Dioscorea opposita Thunb.的根莖，入腎經(jīng)，是《金匱要略》《圣濟總錄》等醫(yī)學(xué)經(jīng)典著作中治療陽痿、消渴等病癥的臨床常用藥[8-10]。本課題組前期研究發(fā)現(xiàn)，山藥總蛋白能夠通過改善腎陽虛模型大鼠睪丸組織纖維化和氧化應(yīng)激，發(fā)揮治療腎陽虛型勃起功能障礙的作用[11]。然而，目前關(guān)于山藥蛋白防治DIED、改善內(nèi)皮功能的研究幾近空白，且發(fā)揮作用的有效物質(zhì)基礎(chǔ)尚不明確。在以往的研究中，多數(shù)是以人臍靜脈內(nèi)皮細胞作為研究高糖所致內(nèi)皮功能損傷的體外模型，但其并不能準確地模擬海綿體內(nèi)皮細胞的微血管環(huán)境。因此，本研究以從山藥總蛋白中分離、純化得到的高純度山藥蛋白——DOT為研究對象，通過建立大鼠DIED模型，考察DOT對模型大鼠的干預(yù)作用，以及DOT對模型大鼠內(nèi)皮功能的影響;在此基礎(chǔ)上，采用Matrigel 3D技術(shù)分離、培養(yǎng)小鼠海綿體內(nèi)皮細胞（MCECs），然后利用高糖誘導(dǎo)致MCECs損傷，在體外模擬糖尿病誘導(dǎo)的內(nèi)皮功能障礙模型，探討DOT對高糖損傷MCECs的恢復(fù)作用;然后，進一步采用轉(zhuǎn)錄組測序（RNA-Seq）技術(shù)，從分子水平分析DOT對內(nèi)皮功能的保護作用機制，從改善內(nèi)皮功能角度為DOT防治DIED的新治療靶點研究提供參考。

1 材料

1.1 主要儀器

AL204型電子天平、FE20型pH計均購自德國? ?Mettler-Toledo公司;Mini Protean型蛋白電泳儀、CFX型實時熒光定量聚合酶鏈式反應(yīng)（qPCR）儀購自美國Bio- Rad公司;iBrightTMFL1000型凝膠成像系統(tǒng)、ACCU- CHEK型血糖儀均購自上海羅氏制藥有限公司;3111型CO2細胞培養(yǎng)箱購自美國Thermo Fisher Scientific公司;EVOS型數(shù)碼倒置顯微鏡購自美國AMG公司;Infnite 200 PRO型酶標儀購自美國Life Sciences公司;HiSeq 2500型基因測序儀購自美國Illumina公司;EVOS FL Auto型全自動熒光顯微鏡購自美國Life Technologies公司;WP-UO-YJ-40型超純水機購自四川沃特爾水處理設(shè)備有限公司。

1.2 主要藥品與試劑

DOT（批號CM202010003，純度99.9%）由長春中醫(yī)藥大學(xué)中醫(yī)藥大健康創(chuàng)新中心制備;乙二胺四乙酸（EDTA，批號20191221）、磷酸鹽緩沖液（PBS，批號20200428）均購自北京索萊寶科技有限公司;一氧化氮合酶（eNOS）、環(huán)磷酸鳥苷（cGMP）、胰島素酶聯(lián)免疫吸附檢測（ELISA）試劑盒（批號均為202104）均購自上海優(yōu)選生物科技有限公司;MCECs完全培養(yǎng)基（批號WH01112005XP）、Hank’s平衡鹽溶液（批號WH0111- 2005XP）均購自武漢普諾賽生命科技有限公司;CCK-8檢測試劑盒（批號I02020141G14）購自美國Invigentech公司;M199培養(yǎng)基（批號RNBJ1753）、青霉素-鏈霉素雙抗（批號A5955）、D-葡萄糖（批號SLCC4951）、鏈脲佐菌素（STZ，批號WXBD4533V）、阿撲嗎啡（APO，批號SLCB4167，純度≥98.5%）均購自美國Sigma公司;Matrigel基質(zhì)膠（批號0048009）購自美國Corning公司;兔源血小板內(nèi)皮細胞黏附分子1（PECAM-1）多克隆抗體（批號ab263899）、大鼠結(jié)蛋白（Desmin）多克隆抗體（批號ab178537）均購自英國Abcam公司;Alexa Fluor 488標記的山羊抗兔免疫球蛋白G（IgG）二抗、Alexa Fluor 594標記的山羊抗大鼠IgG二抗（批號分別為146644、133283）均購自美國Jackson Immuno Research公司;總RNA提取試劑盒（批號W9113）購自天根生化科技（北京）有限公司;RNA逆轉(zhuǎn)錄試劑盒（批號AK71647A）、TB Green? Premix Ex TaqTM試劑盒（批號AJG1875A）均購自日本Takara公司;其余試劑均為分析純或?qū)嶒炇页Ｓ靡?guī)格，水為蒸餾水。

1.3 動物

本研究所用動物為SPF級雄性SD大鼠[鼠齡為4～6周齡，體質(zhì)量為（200±20） g]和SPF級C57BL/6J雄性小鼠[鼠齡為4～6周齡，體質(zhì)量為（18±2） g]，均由長春市億斯實驗動物技術(shù)有限責(zé)任公司提供，實驗動物生產(chǎn)許可證號為SCXK（吉）2018-0007。動物購入后飼養(yǎng)于溫度（23±2） ℃、相對濕度（55±10）%的SPF級動物房中。本研究動物實驗方案經(jīng)長春中醫(yī)藥大學(xué)動物研究倫理委員會批準后實施（倫理批號為2020203）。

2 方法

2.1 DOT對DIED模型大鼠的防治作用考察

2.1.1 造模、分組與給藥 將60只SD大鼠適應(yīng)性飼養(yǎng)3 d后，分為正常對照組（n=10）和造模組（n=50）。其中，正常對照組大鼠飼以正常飼料，造模組大鼠飼以高糖高脂飼料。飼養(yǎng)4周后，正常對照組大鼠腹腔注射檸檬酸鈉緩沖液（pH=4.5），造模組大鼠腹腔注射40 mg/kg STZ誘導(dǎo)建立糖尿病模型（建模方案是在文獻[12-13]的基礎(chǔ)上加以適當修改后確定的）。注射前大鼠均先禁食禁水12 h。腹腔注射7 d后，通過大鼠尾靜脈取血，測定其空腹血糖（FPG）值，以FPG>11.1 mmol/L為糖尿病建模成功標準[13]。取45只糖尿病造模成功的大鼠，按隨機數(shù)字表法隨機分為模型組，DOT低、中、高劑量組（0.3、0.6、0.9 mg/kg，劑量依據(jù)前期預(yù)實驗結(jié)果設(shè)置）和西地那非組（陽性對照組，4.4 mg/kg，劑量依據(jù)文獻[14]設(shè)置），每組9只。正常對照組和模型組大鼠灌胃生理鹽水，各給藥組大鼠灌胃相應(yīng)藥物（均以生理鹽水為溶劑溶解），灌胃體積均為2 mL，每天1次，連續(xù)11周。在實驗期間，正常對照組大鼠繼續(xù)飼以正常飼料;其余各組大鼠繼續(xù)飼以高糖高脂飼料，以誘導(dǎo)DIED的發(fā)生。

2.1.2 指標檢測 實驗期間，每周記錄大鼠體質(zhì)量和FPG值。末次給藥后，將大鼠放置在透明籠中適應(yīng)環(huán)境20 min，然后皮下注射100 μg/mL APO，觀察30 min內(nèi)大鼠陰莖的勃起次數(shù)，并計算勃起率[勃起率（%）=某組APO實驗陽性（即陰莖有勃起）大鼠只數(shù)/同期該組大鼠只數(shù)×100%]。隨后，利用3%戊巴比妥鈉麻醉大鼠，于腹主動脈取血;將血液在4 ℃下靜置20 min，然后以? ? ? 2 500 r/min離心10 min，收集上層血清。按胰島素ELISA試劑盒說明書方法測定大鼠血清中空腹胰島素（FINS）含量，并計算胰島素抵抗指數(shù)（IR）：IR=FPG×FINS/22.5[15]。取大鼠陰莖海綿體組織，在PBS中研碎，然后以3 000 r/min離心5 min，取上清液，分別按對應(yīng)試劑盒說明書方法測定大鼠陰莖海綿體組織中eNOS、cGMP含量。

2.2 MCECs的分離、培養(yǎng)

MCECs的分離、培養(yǎng)方案是在文獻[16]的基礎(chǔ)上加以適當修改后確定的。將C57BL/6J小鼠脫頸處死，分離其陰莖海綿體組織，鏡下除去尿道及神經(jīng)血管束。將干凈的陰莖海綿體組織剪成1～2 mm3的小塊，置于預(yù)冷的24孔板底，每孔中補加200 μL含50 ng/mL血管內(nèi)皮生長因子A（VEGF-A）的Matrigel基質(zhì)膠，包被后切塊;將24孔板置于37 ℃、5%CO2培養(yǎng)箱中培養(yǎng)14 d，在顯微鏡下觀察，直至細胞長滿孔底;然后加入Dispase酶消化1 h，用10 mmol/L EDTA終止消化，以1 100 r/min離心5 min，收集細胞沉淀，置于MCECs完全培養(yǎng)基中培養(yǎng)，取第2～3代細胞進行后續(xù)實驗。

2.3 MCECs的純度鑒定

采用免疫熒光法進行鑒定。將MCECs按2×104? mL-1的密度接種在置有載玻片的6孔板中進行細胞爬片，于37 ℃、5%CO2培養(yǎng)箱中培養(yǎng)2 d后，將爬有MCECs的載玻片置于4%多聚甲醛中，室溫固定15 min;加入含5%山羊血清的0.5% Triton X-100封閉30 min;加入PECAM-1一抗（內(nèi)皮細胞標記物，稀釋比例為1 ∶ 1 000）和Desmin一抗（平滑肌細胞標記物，稀釋比例為1 ∶? ? ? ? ? ? 1 000），于4 ℃孵育過夜;分別加入Alexa Fluor 488標記的山羊抗兔IgG二抗、Alexa Fluor 594標記的山羊抗大鼠IgG二抗（稀釋比例均為1 ∶ 1 000），于室溫孵育1 h;滴加DAPI試劑，于室溫避光染色5 min;用含熒光淬滅劑的封片液進行封片，使用熒光顯微鏡觀察細胞染色情況（MCECs陽性染色呈綠色熒光，陰性染色呈紅色熒光）。統(tǒng)計陽性細胞數(shù)和細胞總數(shù)，并分析MCECs陽性率：陽性率（%）=陽性細胞數(shù)/細胞總數(shù)×100%。

2.4 DOT對高糖損傷MCECs活力的影響

采用CCK-8法進行檢測。取第2～3代MCECs，以MCECs完全培養(yǎng)基制成密度為2×104 mL-1的細胞懸液，然后按每孔100 μL接種至96孔板中，常規(guī)培養(yǎng)。待細胞貼壁后，加入無血清M199培養(yǎng)基饑餓培養(yǎng)24 h，然后將細胞分為正常對照組（MCECs完全培養(yǎng)基）、高糖組（30 mmol/L D-葡萄糖）和DOT低、中、高濃度組（30 mmol/L D-葡萄糖+125、250、500 μg/mL DOT，濃度根據(jù)預(yù)實驗結(jié)果設(shè)置），每組設(shè)置5個復(fù)孔;并設(shè)置空白組（不加細胞）。各組細胞加藥或MCECs完全培養(yǎng)基常規(guī)培養(yǎng)48 h后，在各孔中分別加入10 μL CCK-8試劑，于室溫孵育1 h，然后用酶標儀在450 nm波長處測定各孔的光密度（OD）值，并計算細胞活力：細胞活力（%）=（實驗組OD值-空白組OD值）/（正常對照組OD值-空白組OD值）×100%。實驗重復(fù)3次。

2.5 MCECs中總RNA提取及cDNA文庫構(gòu)建

取第2～3代MCECs，以MCECs完全培養(yǎng)基制成密度為2×104 mL-1的細胞懸液，然后按每孔2 mL接種至6孔板中，常規(guī)培養(yǎng)24 h。棄去原培養(yǎng)基，每孔中加入2 mL無血清M199培養(yǎng)基饑餓培養(yǎng)24 h，隨后將細胞分為正常對照組（MCECs完全培養(yǎng)基）、高糖組（30 mmol/L D-葡萄糖）、DOT組（30 mmol/L D-葡萄糖+250 μg/mL DOT，DOT濃度依據(jù)細胞活力實驗結(jié)果設(shè)置），每組設(shè)置3個復(fù)孔。各組細胞加藥或MCECs完全培養(yǎng)基常規(guī)培養(yǎng)48 h后，收集細胞，按照總RNA提取試劑盒說明書方法提取細胞中總RNA。用帶有Oligo（dT）的磁珠富集真核生物的mRNA，以片段化的mRNA為模板合成雙鏈cDNA，并通過PCR實驗富集得到各組MCECs的cDNA文庫（cDNA文庫的構(gòu)建委托諾禾致源生物科技有限公司完成）。使用基因測序儀對建好的文庫進行測序，對測序得到的原始數(shù)據(jù)（raw data）進行過濾，去除帶接頭（adapter）和低質(zhì)量讀段（reads）的數(shù)據(jù)，得到待分析數(shù)據(jù)（clean reads），并統(tǒng)計測序質(zhì)量值大于20、30的堿基占總體堿基的百分比。

2.6 MCECs差異表達基因篩選及功能富集分析

使用HISAT v2.0.4軟件將clean reads比對到小鼠參考基因組，使用RSEM v1.2.12軟件計算基因和轉(zhuǎn)錄本表達水平。采用possion Dis算法進行差異表達基因檢測，以差異倍數(shù)|fold change（FC）|≥1.2且P≤0.001為標準，分別篩選正常對照組與高糖組、高糖組與DOT組MCECs的差異表達基因，利用R語言包中的ggplot2包繪制差異基因火山圖。通過Cluster Profiler R 4.0.2軟件中的phyper函數(shù)實現(xiàn)差異表達基因的基因本體（GO）功能富集分析和京都基因與基因組百科全書（KEGG）通路富集分析，以錯誤發(fā)現(xiàn)率（FDR）≤0.01為顯著富集[17]。在得到的兩組差異表達基因及顯著富集通路中，篩選出DOT作用后能夠顯著逆轉(zhuǎn)的基因及信號通路，以探討DOT的可能作用靶點及通路。

2.7 關(guān)鍵基因的表達情況驗證

采用qPCR法對與纖維化、氧化應(yīng)激相關(guān)的6個關(guān)鍵基因（TGF-β3、Txnip、Aldh1a1、Loxl1、Mt1、Mt2）在mRNA水平上進行驗證。取“2.6”項下總RNA，按照RNA逆轉(zhuǎn)錄試劑盒方法將總RNA反轉(zhuǎn)錄合成cDNA。以cDNA為模板，按照TB Green? Premix Ex TaqTM試劑盒方法進行PCR擴增（引物序列及擴增產(chǎn)物大小如表1所示）。PCR反應(yīng)條件為：95 ℃預(yù)變性30 s;95 ℃變性5 s，60 ℃退火30 s，72 ℃延伸15 s，共40個循環(huán)。反應(yīng)體系（共25 μL）為：SYBR 12.5 μL，DEPC 8 μL，cDNA模版 2.5 μL，上、下游引物各1 μL。以GAPDH為內(nèi)參，采用2-ΔΔCt法計算各基因的表達水平（式中，Ct表示每個反應(yīng)管內(nèi)熒光信號強度達到設(shè)定閾值時所經(jīng)歷的循環(huán)次數(shù)）。實驗重復(fù)3次。

2.8 統(tǒng)計學(xué)方法

使用GraphPad Prism v 6.0軟件進行統(tǒng)計分析。計量資料以x±s表示，組間比較采用單因素方差分析和Tukey-Kramer檢驗。檢驗水準α=0.05。

3 結(jié)果

3.1 DOT對DIED模型大鼠的防治作用考察結(jié)果

3.1.1 體質(zhì)量、FPG值、FINS和IR測定結(jié)果 造模前和灌胃前，各組大鼠的體質(zhì)量、FPG值差異均無統(tǒng)計學(xué)意義（P>0.05）。灌胃結(jié)束后，與正常對照組比較，模型組大鼠的體質(zhì)量顯著減少（P<0.01），F(xiàn)PG值、FINS和IR均顯著升高（P<0.01）。與模型組比較，DOT低劑量組大鼠的FINS顯著降低（P<0.01）;DOT中、高劑量組和西地那非組大鼠的體質(zhì)量均顯著增加（P<0.01）;DOT中、高劑量組大鼠的FINS、IR均顯著降低（P<0.01）。各組大鼠的體質(zhì)量、FPG值、FINS和IR測定結(jié)果見表2。

3.1.2 勃起情況及陰莖海綿體組織中eNOS、cGMP含量測定結(jié)果 與正常對照組比較，模型組大鼠勃起次數(shù)顯著減少（P<0.01），勃起率和陰莖海綿體組織中eNOS、cGMP含量均顯著降低（P<0.01）。與模型組比較，DOT中劑量組大鼠勃起率和陰莖海綿體組織中eNOS含量均顯著升高（P<0.05或P<0.01）;DOT高劑量組大鼠勃起次數(shù)顯著增加（P<0.01），勃起率和陰莖海綿體組織中eNOS、cGMP含量均顯著升高（P<0.05或P<0.01）;西地那非組大鼠勃起次數(shù)顯著增加（P<0.05），勃起率和陰莖海綿體組織中cGMP含量均顯著升高（P<0.01）。各組大鼠勃起情況和陰莖海綿體組織中eNOS、cGMP含量測定結(jié)果見表3。

3.2 MCECs的純度測定結(jié)果

熒光顯微鏡下可見，細胞核呈現(xiàn)出藍色卵石狀的橢圓形輪廓，絕大多數(shù)細胞呈PECAM-1染色陽性。經(jīng)分析，MCECs的陽性率在90%以上，說明細胞純度較高，可用于后續(xù)實驗。MCECs的免疫熒光染色圖見圖1。

3.3 DOT對高糖損傷MCECs活力的影響考察結(jié)果

與正常對照組比較，高糖組MCECs的活力顯著降低（P<0.01）;與高糖組比較，DOT低、中、高濃度組MCECs的活力均顯著升高（P<0.05或P<0.01），其中以DOT低濃度組MCECs的活力最高。各組MCECs的活力測定結(jié)果見圖2。

3.4 MCECs的cDNA文庫測序數(shù)據(jù)評估結(jié)果

從正常對照組、高糖組和DOT組MCECs中分別得到1.52×108、1.46×108、1.51×108個clean reads，對應(yīng)的總大?。╟lean reads的個數(shù)×長度）分別為7.6 205、7.3 052、7.5 404 Gb。各組clean reads中，測序質(zhì)量值大于20的堿基在所有堿基中的比例≥98%，質(zhì)量值大于30 的堿基在所有堿基中的比例≥94%，GC堿基對的含有量在各樣本中基本保持不變。利用HISAT v2.0.4軟件比對小鼠參考基因組發(fā)現(xiàn)，各組clean reads比對到基因組中的比率在95.52%～96.25%范圍內(nèi)。以上結(jié)果說明，本研究測序質(zhì)量較好，測序數(shù)據(jù)符合后續(xù)分析要求[18]。

3.5 MCECs的差異表達基因篩選及功能富集分析結(jié)果

3.5.1 差異表達基因篩選結(jié)果 與正常對照組比較，高糖組MCECs中共有48個差異表達基因，其中34個基因表達上調(diào)、14個基因表達下調(diào);與高糖組比較，DOT組MCECs中共有779個差異表達基因，其中336個基因表達上調(diào)、443個基因表達下調(diào)。差異表達基因篩選的火山圖見圖3。

3.5.2 GO功能富集分析結(jié)果 通過對正常對照組與高糖組MCECs中48個差異表達基因進行GO功能富集分析，共得到47條GO功能注釋，其中53.19%（25條）的注釋為生物學(xué)過程、29.79%（14條）的注釋為細胞組分、17.02%（8條）的注釋為分子功能。通過對高糖組與DOT組MCECs中779個差異表達基因進行GO功能富集分析，共得到57條GO功能注釋，其中49.12%（28條）的注釋為生物學(xué)過程、31.58%（18條）的注釋為細胞組分、19.30%（11條）的注釋為分子功能。兩組差異表達基因在生物學(xué)過程注釋中均以細胞過程（cellular process）為主，在細胞組分注釋中均以細胞部分（cell part）為主，在分子功能注釋中均以綁定分子功能（binding）為主。差異表達基因的GO功能富集分析結(jié)果見圖4。

3.5.3 KEGG通路富集分析結(jié)果 通過對正常對照組與高糖組MCECs中差異表達基因進行KEGG通路富集分析發(fā)現(xiàn)，48個差異表達基因可以映射到102條不同的生物途徑。其中，顯著富集的途徑有8條，包括蛋白質(zhì)消化吸收（protein digestion and absorption）、甘油脂類新陳代謝（glycerolipid metabolism）、甲狀腺激素合成（thyroid hormone synthesis）等。通過對高糖組與DOT組MCECs中差異表達基因進行KEGG通路富集分析發(fā)現(xiàn)，779個差異表達基因可以映射到308條不同的生物途徑。其中，顯著富集的途徑有84條，包括黏著斑通路（focal adhesion）、細胞外基質(zhì)受體相互作用通路（ECM-receptor interaction）、癌癥通路（pathways in cancer）等。MCECs差異表達基因的KEGG通路富集分析結(jié)果見圖5。

3.5.4 DOT對高糖損傷MCECs的多靶點調(diào)控分析結(jié)果 從正常對照組與高糖組、高糖組與DOT組MCECs中共篩選得到Aldh1a1、Abcc5、Tac1等13個共同的差異表達基因。這些差異表達基因產(chǎn)生的調(diào)控作用涉及細胞增殖、細胞遷移、氧化還原、應(yīng)激反應(yīng)、代謝等生物過程，并涉及細胞外基質(zhì)受體相互作用通路、肌動蛋白細胞骨架調(diào)控通路（regulation of actin cytoskeleton）、錯配修復(fù)通路（mismatch repair）、磷脂酰肌醇3激酶-蛋白激酶B（PI3K-Akt）信號通路等4條顯著富集通路。與模型組比較，DOT組MCECs中上述差異表達基因的表達均被顯著回調(diào)，結(jié)果見表4。

3.6 MCECs差異表達基因的驗證結(jié)果

與正常對照組比較，高糖組MCECs中Txnip、? ? ? ?Aldh1a1、Loxl1、TGF-β3等基因的mRNA表達均顯著上調(diào)（P<0.05或P<0.01），Mt1、Mt2等基因的mRNA表達均顯著下調(diào)（P<0.05或P<0.01）;與高糖組比較，DOT組MCECs中上述基因的mRNA表達均被顯著逆轉(zhuǎn)（P<0.01）。該驗證結(jié)果與RNA-Seq分析結(jié)果一致，說明測序結(jié)果真實可信。DOT對高糖損傷MCECs中6個關(guān)鍵基因mRNA表達的驗證結(jié)果見表5。

4 討論

隨著糖尿病病程的延長，糖尿病患者的勃起功能進行性下降，造成陰莖組織出現(xiàn)嚴重病理損傷，而高血糖引發(fā)的內(nèi)皮功能障礙是誘發(fā)、加重糖尿病患者DIED的關(guān)鍵因素[18]。因此，發(fā)現(xiàn)可有效改善內(nèi)皮功能障礙的活性物質(zhì)已成為防治DIED的重中之重。為了解決上述問題，本研究首先建立糖尿病大鼠模型，并飼以高糖高脂飼料來誘發(fā)DIED，初步探究DOT對DIED的防治作用，并選擇臨床治療勃起功能障礙的一線藥物西地那非作為陽性對照藥物。

體質(zhì)量減輕、血糖和IR升高、勃起功能異常是DIED的主要表現(xiàn)。控制血糖是DIED的治療手段之一，可防止高血糖對血管功能的進一步損害;IR升高可致使胰島素不能發(fā)揮其生物學(xué)效應(yīng)，使機體靶組織器官對葡萄糖的攝取和利用能力降低，從而導(dǎo)致內(nèi)皮功能障礙[19-20]。本研究結(jié)果顯示，與模型組比較，DOT各劑量組大鼠的體質(zhì)量、勃起率均不同程度地升高，IR均不同程度地降低。這提示，DOT對DIED模型大鼠的癥狀具有一定的改善作用。相關(guān)勃起生理機制研究表明，內(nèi)皮細胞可在eNOS催化下合成并分泌活性物質(zhì)NO，并進入海綿體平滑肌抑制cGMP水解，引起海綿體平滑肌松弛，從而促使陰莖勃起[21]。本研究結(jié)果顯示，高劑量DOT可顯著升高DIED模型大鼠陰莖海綿體組織中eNOS、cGMP含量，表明DOT對DIED的防治作用可能是通過改善內(nèi)皮功能實現(xiàn)的。

為進一步探究DOT改善高糖所致內(nèi)皮功能紊亂的作用與機制，本研究選擇位于海綿體內(nèi)表面的原代MCECs建立體外細胞模型來進行相關(guān)研究。MCECs是維持海綿體功能最重要的細胞之一，是研究DIED內(nèi)皮功能的最佳選擇[22-23]。在原代MCECs的分離方式上，本研究選擇了一種新穎的非酶分離方法——Matrigel 3D技術(shù)。該技術(shù)能模擬體內(nèi)細胞生長的三維環(huán)境，使MCECs與生長因子直接接觸而誘導(dǎo)MCECs從組織中爬出，可保證MCECs的原有形態(tài)及功能特性，且方法操作省時、分離純度高、重復(fù)性好[24]。在以上研究基礎(chǔ)上，本研究又利用高糖誘導(dǎo)模擬糖尿病環(huán)境下MCECs損傷狀態(tài)，并證明了DOT能夠通過提高高糖損傷MCECs的活力來發(fā)揮其改善內(nèi)皮功能障礙的作用。

轉(zhuǎn)錄組學(xué)可從整體水平上研究基因的表達情況，從“多靶點-多通路-多途徑”闡明藥物的作用機制[25-26]。為深入分析DOT改善內(nèi)皮功能的可能作用途徑，本研究對MCECs的轉(zhuǎn)錄組差異表達基因進行重點分析。結(jié)果顯示，與正常對照組比較，高糖組MCECs有48個差異表達基因;與高糖組比較，DOT組MCECs有779個差異表達基因。進一步通過GO功能富集分析和KEGG通路富集分析發(fā)現(xiàn)，DOT對內(nèi)皮功能障礙的恢復(fù)作用是通過調(diào)控多種功能基因和信號通路實現(xiàn)的。為明確DOT發(fā)揮作用的靶點及生理過程，本研究對獲得的兩組細胞（正常細胞與高糖損傷細胞、高糖損傷細胞與DOT處理細胞）中差異表達基因數(shù)據(jù)進行深度挖掘分析，發(fā)現(xiàn)有13個重疊基因及4條顯著富集信號通路，這些基因及信號通路可能是DOT改善內(nèi)皮功能障礙、防治DIED的關(guān)鍵靶點和途徑。綜合以上分析及本課題組前期研究成果，本研究篩選出了能被DOT顯著調(diào)控的6個關(guān)鍵基因進行深入分析，包括與纖維化、氧化應(yīng)激相關(guān)的TGF-β3、Txnip、Aldh1a1、Loxl1、Mt1、Mt2。

膠原合成增加引起的纖維化因子表達上調(diào)和海綿狀纖維化是誘發(fā)勃起功能障礙的重要病理、生理機制[27]。在眾多纖維化因子中，轉(zhuǎn)化生長因子β（TGF-β）被認為是與纖維生成最為相關(guān)的細胞因子，其可引起組織纖維化，從而導(dǎo)致多種器官的功能損害[28]。TGF-β家族包括TGF-β1、TGF-β2、TGF-β3 3種亞基，其中TGF-β3能夠促使膠原大量沉積，進而引起纖維化的發(fā)生[29]。此外，TGF-β1可通過激活其下游蛋白表達，造成內(nèi)皮細胞纖維化，從而抑制MCECs的遷移和增殖[30]。本研究發(fā)現(xiàn)，經(jīng)DOT處理后，高糖損傷MCECs中TGF-β3基因的mRNA表達顯著下調(diào)，表示DOT可能是通過抑制TGF- β3基因的表達來抑制內(nèi)皮細胞纖維化，從而發(fā)揮其改善內(nèi)皮功能障礙的作用。

據(jù)報道，長期升高的葡萄糖水平可促使細胞產(chǎn)生活性氧（ROS）。當ROS產(chǎn)生后，一方面其可干擾NO的生物利用度，介導(dǎo)內(nèi)皮功能障礙;另一方面其可誘導(dǎo)內(nèi)皮細胞凋亡，進而破壞血管組織內(nèi)穩(wěn)態(tài)，影響血管舒張和陰莖勃起[31-32]。硫氧還蛋白相互作用蛋白（Txnip）是調(diào)節(jié)糖和脂代謝的中樞蛋白，其過表達可引起ROS水平升高、NO水平降低，進而啟動并放大氧化應(yīng)激反應(yīng)[33-34]。乙醛脫氫酶1家族成員a1（Aldh1a1）是一種多功能酶，其在氧化應(yīng)激狀態(tài)下可在細胞中聚集，從而緩解氧化應(yīng)激介導(dǎo)的細胞損傷;氧化應(yīng)激水平降低后，細胞中Aldh1a1的表達量也會隨之減少[35-36]。類賴氨酰氧化酶1（Loxl1）是保持細胞內(nèi)環(huán)境穩(wěn)定的關(guān)鍵酶，在氧化應(yīng)激狀態(tài)下細胞中Loxl1的含量會明顯升高[37]。金屬硫蛋白1（Mt1）和Mt2作為細胞內(nèi)的抗氧化因子可直接發(fā)揮清除ROS的作用，參與了生物體的多種保護性應(yīng)激反應(yīng)[38-41]。本研究結(jié)果顯示，DOT能夠顯著逆轉(zhuǎn)高糖損傷MCECs中Txnip、Aldh1a1、Loxl1等基因mRNA表達的上調(diào)，以及Mt1、Mt2等基因mRNA表達的下調(diào)。以上結(jié)果提示，DOT對內(nèi)皮功能的改善作用與降低氧化應(yīng)激反應(yīng)相關(guān)。

綜上所述，DOT對DIED模型大鼠的癥狀具有改善作用;其作用機制可能與抑制纖維化、降低氧化應(yīng)激等生物學(xué)途徑，從而改善海綿體內(nèi)皮功能有關(guān)。本研究為DOT改善高糖引起的內(nèi)皮功能障礙提供了治療靶點，也為DOT防治DIED的發(fā)生發(fā)展提供了實驗數(shù)據(jù)支撐。但本研究對于DOT防治DIED的分子機制研究尚不深入，有待進一步完善。

參考文獻

[ 1 ] MALAVIGE L S，LEVY J C. Erectile dysfunction in diabetes mellitus[J]. J Sex Med，2009，6（5）：1232-1247.

[ 2 ] WANG J S，F(xiàn)ENG J L，LI X，et al. Effect of leech-centipede medicine on improving erectile function in diabetes- induced erectile dysfunction rats via PDE5 signalling pathway-related molecules[J]. Pharm Biol，2021，59（1）：167- 174.

[ 3 ] KOUIDRAT Y，PIZZOL D，COSCO T，et al. High prevalence of erectile dysfunction in diabetes：a systematic review and meta-analysis of 145 studies[J]. Diabet Med，2017，34（9）：1185-1192.

[ 4 ] BELBA A，CORTELAZZO A，ANDREA G，et al. Erectile dysfunction and diabetes：association with the impairment of lipid metabolism and oxidative stress[J]. Clin Biochem，2016，49（1/2）：70-78.

[ 5 ] MARTíNEZ-JABALOYAS J M，GIL-SALOM M， VILLAMóN-FORT R，et al. Prognostic factors for re- sponse to sildenafil in patients with erectile dysfunction[J]. Eur Urol，2001，40（6）：641-646.

[ 6 ] LASKER G F，MALEY J H，KADOWITZ P J. A review of the pathophysiology and novel treatments for erectile dysfunction[J]. Adv Pharmacol Sci，2010，2010：730861.

[ 7 ] RYU J K，SUH J K，BURNETT A L. Research in pharmacotherapy for erectile dysfunction[J]. Transl Androl Urol，2017，6（2）：207-215.

[ 8 ] 周瀅，周萍.山藥在各歷史時期臨床運用的研究[J].世界科學(xué)技術(shù)（中醫(yī)藥現(xiàn)代化），2014，16（7）：1574-1577.

[ 9 ] 祝之友.張仲景應(yīng)用山藥情況[J].中國中醫(yī)藥現(xiàn)代遠程教育，2017，15（7）：141.

[10] 曹雯，喻嶸，王琦威，等.《古今名醫(yī)臨證金鑒·消渴卷》治療消渴病用藥規(guī)律分析[J].中華中醫(yī)藥學(xué)刊，2016，34（5）：1099-1011.

[11] YU S，HAN B，BAI X，et al. The cold-soaking extract of Chinese yam （Dioscorea opposita Thunb.） protects against erectile dysfunction by ameliorating testicular function in hydrocortisone-induced KDS-Yang rats and in oxidatively damaged TM3 cells[J]. J Ethnopharmacol，2020，263：113223.

[12] ZHANG J，QIU H，HUANG J，et al. Establishment of a? diabetic myocardial hypertrophy model in Mus musculus castaneus mouse[J]. Int J Exp Pathol，2018，99（6）：295- 303.

[13] WANG X，PAN J，LIU D，et al. Nicorandil alleviates apoptosis in diabetic cardiomyopathy through PI3K/Akt pathway[J]. J Cell Mol Med，2019，23（8）：5349-5359.

[14] ZHANG H，TONG W T，ZHANG C R，et al. Gross saponin of Tribulus terrestris improves erectile dysfunction in type 2 diabetic rats by repairing the endothelial function of the penile corpus cavernosum[J]. Diabetes Metab Syndr Obes，2019，12：1705-1716.

[15] 郭璇.輔助降糖顆粒與西格列汀聯(lián)合用藥對ZDF大鼠胰島素抵抗作用及機制研究[D].北京：北京中醫(yī)藥大學(xué)，2019.

[16] YIN G N，PARK S H，SONG K M，et al. Establishment of in vitro model of erectile dysfunction for the study of high-glucose-induced angiopathy and neuropathy[J]. Andrology，2017，5（2）：327-335.

[17] WANG L，F(xiàn)ENG Z，WANG X，et al. DEGseq：an R pa- ckage for identifying differentially expressed genes from RNA-seq data[J]. Bioinformatics，2010，26（1）：136-138.

[18] 羅永珍.基于轉(zhuǎn)錄組學(xué)和蛋白質(zhì)組學(xué)的小鼠氟中毒腦損傷研究[D].蘭州：西北師范大學(xué)，2020.

[19] VIDELA S，VILLORIA J，SUST M，et al. Calcium dobesilate for the treatment of erectile dysfunction in men with diabetes mellitus[J]. Basic Clin Pharmacol Toxicol，2008，103（4）：360-366.

[20] SENA C M，PEREIRA A M，SEI?A R. Endothelial dysfunction：a major mediator of diabetic vascular disease[J].Biochim Biophys Acta，2013，1832（12）：2216-2231.

[21] GUR S，ALZWERI L，YILMAZ-ORAL D，et al. Testo- sterone positively regulates functional responses and nitric oxide expression in the isolated human corpus cavernosum[J]. Andrology，2020，8（6）：1824-1833.

[22] KUMAR S，WEST D C，AGER A. Heterogeneity in endothelial cells from large vessels and microvessels[J]. Diffe- rentiation，1987，36（1）：57-70.

[23] BERNAS M J，CARDOSO F L，DALEY S K，et al. Establishment of primary cultures of human brain microvascular endothelial cells to provide an in vitro cellular model of the blood-brain barrier[J]. Nat Protoc，2010，5（7）：1265-1272.

[24] YIN G N，RYU J K，KWON M H，et al. Matrigel-based sprouting endothelial cell culture system from mouse corpus cavernosum is potentially useful for the study of endothelial and erectile dysfunction related to high-glucose exposure[J]. J Sex Med，2012，9（7）：1760-1772.

[25] CARR H L，TURNER J D，MAJOR T，et al. New developments in transcriptomic analysis of synovial tissue[J]. Front Med （Lausanne），2020，7：21.

[26] 張曉萌，李健春，王瓊，等.轉(zhuǎn)錄組測序技術(shù)在中醫(yī)藥領(lǐng)域的應(yīng)用[J].中國現(xiàn)代中藥，2016，18（8）：1084-1087.

[27] RYU J K，SONG S U，HAN J Y，et al. Establishment of penile fibrosis model in a rat using mouse NIH 3T3 fibroblasts expressing transforming growth factor beta1[J].? ? Biol Reprod，2005，72（4）：916-921.

[28] SUN T，HUANG Z，LIANG W C，et al. TGFβ2 and TGFβ3 isoforms drive fibrotic disease pathogenesis[J]. Sci Transl Med，2021，13（605）：0407.

[29] RYU C M，SHIN J H，YU H Y，et al. N-acetylcysteine prevents bladder tissue fibrosis in a lipopolysaccharide-induced cystitis rat model[J]. Sci Rep，2019，9（1）：8134.

[30] SHIN T Y，RYU J K，JIN H R，et al. Increased cavernous expression of transforming growth factor-β1 and activation of the Smad signaling pathway affects erectile dysfunction in men with spinal cord injury[J]. J Sex Med，2011，8（5）：1454-1462.

[31] CHEN Y，ZHOU B，YU Z，et al. Baicalein alleviates erectile dysfunction associated with streptozotocin-induced type Ⅰ diabetes by ameliorating endothelial nitric oxide synthase dysfunction，inhibiting oxidative stress and fibrosis

[J]. J Sex Med，2020，17（8）：1434-1447.

[32] CASTELA ?，COST A C. Molecular mechanisms associa- ted with diabetic endothelial-erectile dysfunction[J]. Nat Rev Urol，2016，13（5）：266-274.

[33] ZHANG B，LYU J，YANG E J，et al. Class Ⅰ histone deacetylase inhibition is synthetic lethal with BRCA1 deficiency in breast cancer cells[J]. Acta Pharm Sin B，2020，10（4）：615-627.

[34] LI X，KOVER K L，HERUTH D P，et al. Thioredoxin-interacting protein promotes high-glucose-induced macrovascular endothelial dysfunction[J]. Biochem Biophys Res Commun，2017，493（1）：291-297.

[35] JACKSON B，BROCKER C，THOMPSON D C，et al. Update on the aldehyde dehydrogenase gene （ALDH） superfamily[J]. Hum Genomics，2011，5（4）：283-303.

[36] SINGH S，BROCKER C，KOPPAKA V，et al. Aldehyde dehydrogenases in cellular responses to oxidative/electrophilic stress[J]. Free Radic Biol Med，2013，56：89-101.

[37] GREENE A G，EIVERS S B，DERVAN E W J，et al. Lysyl oxidase like 1：biological roles and regulation[J]. Exp Eye Res，2020，193：107975.

[38] COYLE P，PHILCOX J C，CAREY L C，et al. Metallothionein：the multipurpose protein[J]. Cell Mol Life Sci，2002，59（4）：627-647.

[39] RUTTKAY-NEDECKY B，NEJDL L，GUMULEC J，et al. The role of metallothionein in oxidative stress[J]. Int J Mol Sci，2013，14（3）：6044-6066.

[40] MIURA T，MURAOKA S，OGISO T. Antioxidant activity of metallothionein compared with reduced glutathione[J].Life Sci，1997，60（21）：301-309.

[41] SARASWATHI V，GANESAN M，PERRIOTTE-OLSON C，et al. Nanoformulated copper/zinc superoxide dismutase attenuates vascular cell activation and aortic inflammation in obesity[J]. Biochem Biophys Res Commun，2016，469（3）：495-500.

（收稿日期：2021-06-15 修回日期：2021-10-23）

（編輯：林 靜）