陳嘉敏，李 瑩，，吳會(huì)會(huì)，劉 鵬，鄭 燕，，蘇國(guó)海

1山東大學(xué)齊魯醫(yī)學(xué)院 濟(jì)南市中心醫(yī)院專(zhuān)科轉(zhuǎn)化研究中心，濟(jì)南 250013

2山東第一醫(yī)科大學(xué)附屬中心醫(yī)院專(zhuān)科轉(zhuǎn)化研究中心，濟(jì)南 250013

糖尿病心肌病是指在沒(méi)有冠狀動(dòng)脈疾病、高血壓和心臟瓣膜病的情況下發(fā)生心力衰竭的病理生理疾病[1]。歐洲及美國(guó)的多項(xiàng)回顧性研究顯示糖尿病增加了心力衰竭患者的住院率和死亡風(fēng)險(xiǎn)[2]，在我國(guó)，糖尿病導(dǎo)致的心力衰竭發(fā)病率在1993至2007 年從12.3%增加到22.1%，糖尿病是除冠心病、高血壓外導(dǎo)致心力衰竭最主要的原因[3]。2017至2018年中國(guó)心力衰竭住院患者的橫斷面研究顯示心力衰竭住院患者合并糖尿病的比例高達(dá)29.2%[4]，在包括中國(guó)人在內(nèi)的多種族的研究顯示，糖尿病增加了心力衰竭患者的發(fā)病率和死亡率[5]，近30年來(lái)我國(guó)糖尿病的患病人數(shù)逐年上漲[6]，使糖尿病心肌病的防治工作面臨更加嚴(yán)峻的挑戰(zhàn)。研究表明單純將糖尿病患者的血糖恢復(fù)到基線水平并不能降低心力衰竭的發(fā)病率和死亡率[7]，但是新型降糖藥鈉葡萄糖協(xié)同轉(zhuǎn)運(yùn)蛋白2抑制劑如達(dá)格列凈可以使心力衰竭合并糖尿病患者的心力衰竭惡化風(fēng)險(xiǎn)及因心血管事件死亡的風(fēng)險(xiǎn)降低[8]，心血管藥物沙庫(kù)巴曲-纈沙坦也可使糖尿病合并心力衰竭的患者獲益[9]。多種機(jī)制參與了糖尿病心肌病的進(jìn)展：(1)心肌細(xì)胞胰島素信號(hào)通路受損：正常情況下，磷酸肌醇3激酶-蛋白激酶B信號(hào)通路激活，刺激葡萄糖轉(zhuǎn)運(yùn)體4(glucose transporter type 4，GLUT4)向質(zhì)膜募集，促進(jìn)葡萄糖進(jìn)入心臟細(xì)胞。胰島素信號(hào)通路受損時(shí)，GLUT4表達(dá)下降，心臟攝取葡萄糖的能力下降。(2)高血糖和糖毒性：高血糖可導(dǎo)致糖基化終末產(chǎn)物增加，糖基化終末產(chǎn)物是長(zhǎng)壽命蛋白被糖化的產(chǎn)物，糖基化終末產(chǎn)物沉積使結(jié)締組織交聯(lián)、纖維化，導(dǎo)致心臟順應(yīng)性降低、舒張功能障礙。(3)脂肪酸利用增加和脂毒性：糖尿病患者葡萄糖利用下降，參與脂肪酸β氧化的信號(hào)轉(zhuǎn)導(dǎo)分子過(guò)氧化物酶體增殖物激活受體(peroxisome proliferator-activated receptor，PPAR)α表達(dá)增加促進(jìn)脂肪酸的攝取和利用，脂質(zhì)和脂質(zhì)代謝產(chǎn)物的增加導(dǎo)致胰島素抵抗和心臟纖維化。(4)氧化應(yīng)激：脂肪酸的β氧化促進(jìn)活性氧的產(chǎn)生和心臟的氧化應(yīng)激。(5)Ca2+處理受損：高血糖狀態(tài)下，內(nèi)質(zhì)網(wǎng)應(yīng)激使Ca2+處理異常，導(dǎo)致心肌細(xì)胞興奮性收縮偶聯(lián)障礙，心臟收縮能力下降[1，10- 12]。本研究旨在通過(guò)對(duì)現(xiàn)有數(shù)據(jù)庫(kù)數(shù)據(jù)的再分析，尋找導(dǎo)致糖尿病心肌病的關(guān)鍵基因，對(duì)糖尿病心肌病的發(fā)病機(jī)制有更深一步的研究。

材料和方法

GEO數(shù)據(jù)庫(kù)選取基因集進(jìn)行差異基因分析從GEO數(shù)據(jù)庫(kù)中選取與糖尿病心肌病研究相關(guān)的基因芯片。納入標(biāo)準(zhǔn)：(1)種屬：Rattus norvegicus；(2)同時(shí)具備正常心臟組織和糖尿病心肌病心臟組織。經(jīng)過(guò)檢索GSE4745[13](https：//www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE4745)和GSE5606[14](https：//www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE5606)兩個(gè)基因集符合條件。GSE4745基因集的芯片來(lái)自GPL85(Affymetrix Rat Genome U34 Array)平臺(tái)，數(shù)據(jù)集中包括正常組雄性大鼠心室組織和注射鏈脲佐菌素后3、28、42 d后糖尿病組雄性大鼠心室組織共24 份，注射鏈脲佐菌素42 d后的糖尿病組大鼠心室組織以及其對(duì)照組樣本共8 份用于差異表達(dá)基因的分析。GSE5606基因集的芯片來(lái)自GPL1355(Affymetrix Rat Genome 230 2.0 Array)平臺(tái)，數(shù)據(jù)集中包括注射鏈脲佐菌素16周的雄性大鼠正常組和糖尿病組心室組織樣本共14 份，全部用于差異基因的分析。GEO數(shù)據(jù)庫(kù)的在線分析工具GEO2R用于差異基因的分析，以P≤0.05，|Log|差異倍數(shù)(fold change，F(xiàn)C)|≥1為標(biāo)準(zhǔn)進(jìn)行篩選，在R[15](3.6.3版本)中利用ggplot2[16]包(3.3.2版本)繪制差異基因的火山圖，在線韋恩圖繪制工具(http：//bioinformatics.psb.ugent.be/webtools/Venn/)繪制韋恩圖，尋找兩個(gè)數(shù)據(jù)集中共同表達(dá)的上調(diào)和下調(diào)差異基因。

富集分析將差異表達(dá)基因利用R中的AnnotationHub包[17]、org.Rn.eg.db包[18]進(jìn)行基因ID的轉(zhuǎn)換，clusterProfile包[19]進(jìn)行基因本體論(gene ontology，GO)和京都基因和基因組百科全書(shū)(Kyoto encyclopedia of genes and genome，KEGG)富集分析，ggplot2繪圖，矯正后P≤0.05為有統(tǒng)計(jì)學(xué)意義?；蚣患治?gene set enrichment analysis，GSEA)是一種計(jì)算方法，用于確定一組預(yù)先確定的基因集在兩種生物狀態(tài)之間是否顯示出差異有統(tǒng)計(jì)學(xué)意義，GSEA軟件[20- 21](版本4.1.0，https：//www.gsea-msigdb.org/gsea/datasets.jsp)用于基因集富集分析，依據(jù)GSEA官網(wǎng)提供的文件格式[Data formats-GeneSetEnrichmentAnalysisWiki (broadinstitute.org)]構(gòu)建文件進(jìn)行GSEA，基因表達(dá)數(shù)據(jù)被分為正常和糖尿心肌病(diabetic cardiomyopathy，DCM)兩種生物狀態(tài)，PPAR信號(hào)通路基因集和脂肪酸代謝基因集[GSEA|MSigDB(gsea-msigdb.org)]作為參考基因集用于富集分析，名義P<0.05、錯(cuò)誤發(fā)現(xiàn)率即q<0.25、標(biāo)準(zhǔn)化富集評(píng)分絕對(duì)值>1被認(rèn)為有統(tǒng)計(jì)學(xué)意義。

蛋白質(zhì)互作網(wǎng)絡(luò)構(gòu)建利用STRING(https：//string-db.org)進(jìn)行蛋白質(zhì)互作網(wǎng)絡(luò)的構(gòu)建，將差異表達(dá)基因?qū)隨TRING數(shù)據(jù)庫(kù)，選擇Rattus norvegicus進(jìn)行蛋白互作網(wǎng)絡(luò)的構(gòu)建，將結(jié)果以TSV格式導(dǎo)出。將互作網(wǎng)絡(luò)中節(jié)點(diǎn)出現(xiàn)的次數(shù)進(jìn)行統(tǒng)計(jì)并排序，將STRING導(dǎo)出的TSV文件以“Network”方式導(dǎo)入Cytoscape，將排序文件以“Table”方式導(dǎo)入，將節(jié)點(diǎn)大小及顏色以“Mapping Type”方式依據(jù)節(jié)點(diǎn)出現(xiàn)次數(shù)多少進(jìn)行顏色填充和節(jié)點(diǎn)大小的調(diào)整。利用插件Cytohubba的最大集團(tuán)中心算法計(jì)算出排名前10的關(guān)鍵基因，以得分排序并調(diào)整節(jié)點(diǎn)大小繪圖。

原代心肌細(xì)胞的提取取出生3 d內(nèi)的Wistar大鼠乳鼠，75%的酒精將乳鼠短暫浸泡消毒，眼科剪剪開(kāi)乳鼠前胸，鑷子取出心臟置于含有1%青鏈霉素的Hank’s平衡鹽溶液(貨號(hào)：CC106，中科邁晨科技有限公司)中，剪碎心臟后在含有200 U/ml 膠原酶(貨號(hào)：LS004176，美國(guó)Worthington)的Hank’s平衡鹽溶液中進(jìn)行消化，預(yù)消化5 min后棄去膠原酶，重新添加膠原酶后消化1 h，取上清加入等量的心肌細(xì)胞培養(yǎng)基終止消化并離心取沉淀，余下組織再次進(jìn)行消化，重復(fù)以上步驟直至無(wú)組織團(tuán)塊，將離心取得的沉淀用培養(yǎng)基吹起后置于培養(yǎng)瓶中差速貼壁80 min后除去成纖維細(xì)胞，心肌細(xì)胞種于0.3%明膠(貨號(hào)：9000708，默克公司)包被后的6孔板中。

心肌細(xì)胞培養(yǎng)原代心肌細(xì)胞貼壁18 h后換為含有100 μmol/L的5-溴脫氧尿嘧啶核苷(貨號(hào)：S7918，Selleck)的心肌細(xì)胞培養(yǎng)基中以抑制成纖維細(xì)胞的生長(zhǎng)，心肌細(xì)胞培養(yǎng)基為8%的馬血清(貨號(hào)：16050122，Gibco)、5%的胎牛血清(貨號(hào)：10270106，Gibco)、1%的青鏈霉素和64%的低糖5.5 mmol/L DMEM的低糖培養(yǎng)基(貨號(hào)：CM1004，中科邁晨科技有限公司)構(gòu)成。原代心肌細(xì)胞在含有Brdu的低糖培養(yǎng)基培養(yǎng)72 h后正常組更換為不含Brdu的低糖培養(yǎng)基，實(shí)驗(yàn)組更換為不含Brdu的25 mmol/L高糖培養(yǎng)基(貨號(hào)：C11995500BT，Gibco)，72 h后提取RNA。

實(shí)時(shí)熒光定量PCR從6孔板收集貼壁心肌細(xì)胞，用Trizol(貨號(hào)：AG21102，艾科瑞生物工程有限公司)提取總RNA，使用反轉(zhuǎn)錄試劑盒(貨號(hào)：AG11711，艾科瑞生物工程有限公司)反轉(zhuǎn)錄后利用SYBERGeeen qPCR試劑盒(貨號(hào)：AG11718，艾科瑞生物工程有限公司)進(jìn)行實(shí)時(shí)熒光定量PCR，每次實(shí)驗(yàn)設(shè)置3 個(gè)復(fù)孔，實(shí)驗(yàn)重復(fù)3 次。引物利用美國(guó)國(guó)家生物技術(shù)信息中心的Primer-BLAST工具設(shè)計(jì)，F(xiàn)ASTA序列從美國(guó)國(guó)家生物技術(shù)信息中心的基因數(shù)據(jù)庫(kù)獲取，產(chǎn)物長(zhǎng)度限制在100～300 bp，溶解溫度設(shè)定為57～63 ℃，以β-actin作為內(nèi)參基因，引物序列見(jiàn)表1。

統(tǒng)計(jì)學(xué)處理PCR以2-△△Ct法計(jì)算目的基因表達(dá)水平，利用GraphPad(8.0.1)進(jìn)行統(tǒng)計(jì)分析，采用獨(dú)立樣本t檢驗(yàn)，P≤0.05為差異有統(tǒng)計(jì)學(xué)意義。

表1 關(guān)鍵基因的引物序列

結(jié) 果

GEO2R與韋恩圖分析結(jié)果在GEO數(shù)據(jù)庫(kù)中選擇數(shù)據(jù)集GSE4745和GSE5606作為研究對(duì)象，利用GEO2R進(jìn)行分析后，以|LogFC|≥1，P≤0.05為閾值篩選差異基因并利用ggplot2繪制火山圖，其中LogFC≥1，P≤0.05為上調(diào)基因，LogFC≤-1，P≤0.05為下調(diào)基因。在GSE4745中以GSM107428、GSM107429、GSM107430、GSM107431作為對(duì)照組，GSM107432、GSM107433、GSM107434、GSM107435作為糖尿病心肌病組，經(jīng)GEO2R分析、ggplot2繪圖后顯示上調(diào)基因164個(gè)、下調(diào)基因164個(gè)(圖1A)。在GSE5606中以GSM130860、GSM130861、GSM130862、GSM130863、GSM130864、GSM130865、GSM130866作為對(duì)照組，以GSM130867、GSM130868、GSM130869、GSM130870、GSM130871、GSM130872、GSM130873作為糖尿病心肌病組，經(jīng)GEO2R分析、ggplot2繪圖顯示上調(diào)基因286個(gè)，下調(diào)基因135個(gè)(圖1B)。在線韋恩圖繪制工具分析后顯示共同表達(dá)的上調(diào)差異基因共35個(gè)(圖1C)，共同表達(dá)的下調(diào)差異基因共15個(gè)(圖1D)。GSE5606和GSE4745基因集中共同表達(dá)的差異基因的LogFC值和P值顯示，GSE5606和GSE4745基因集中共同表達(dá)的上調(diào)基因共35個(gè)，共同表達(dá)的下調(diào)基因共15個(gè)(表2)。

基因富集分析結(jié)果在R中利用clusterProfile包，將50個(gè)差異基因進(jìn)行GO和KEGG富集，矯正后P≤0.05為有統(tǒng)計(jì)學(xué)意義。這50個(gè)差異基因在GO富集中的分子功能、生物學(xué)過(guò)程有明顯富集，而在細(xì)胞組分無(wú)富集。在生物學(xué)過(guò)程中，這50個(gè)差異基因在脂肪酸代謝過(guò)程、小分子分解過(guò)程、嘌呤核苷酸代謝過(guò)程、脂質(zhì)修飾等生物過(guò)程中明顯富集(圖2A)。在分子功能中，差異基因主要在棕櫚酰輔酶A水解酶活性、酰基輔酶A水解酶活性、輔酶A水解酶活性、脂肪酸衍生物的結(jié)合、硫酯水解酶活性這些條目中明顯富集(圖2B)。在KEGG，富集主要集中在PPAR信號(hào)通路、脂肪酸的延伸、不飽和脂肪酸的生物合成(圖2C)。對(duì)脂肪酸代謝基因集和PPAR信號(hào)通路基因集進(jìn)行GSEA富集顯示，糖尿病組在脂肪酸代謝和PPAR信號(hào)通路的基因集中明顯富集(圖3)。

FC：差異倍數(shù)

STRING在線數(shù)據(jù)庫(kù)蛋白互作網(wǎng)絡(luò)分析結(jié)果利用STRING在線數(shù)據(jù)庫(kù)(https：//www.string-db.org/)將GSE4745和GSE5606數(shù)據(jù)集共同表達(dá)的差異基因構(gòu)建蛋白互作網(wǎng)絡(luò)，利用Cytoscape繪制出蛋白質(zhì)的互作關(guān)系網(wǎng)絡(luò)(圖4A)。利用插件Cytohubba的最大集團(tuán)中心算法計(jì)算出排名前10的關(guān)鍵基因并繪制關(guān)鍵基因的互作網(wǎng)絡(luò)圖(圖4B)。等級(jí)評(píng)分見(jiàn)表3。

關(guān)鍵基因在原代大鼠心肌細(xì)胞中的驗(yàn)證結(jié)果原代心肌細(xì)胞對(duì)照組在5.5 mmol/L的葡萄糖低糖環(huán)境、實(shí)驗(yàn)組在25 mmol/L的葡萄糖高糖環(huán)境中培養(yǎng)72 h后，提取RNA進(jìn)行實(shí)時(shí)熒光定量PCR。t檢驗(yàn)結(jié)果顯示，Pdk4(t=23.520，P<0.001)、Hmgcs2(t=10.700，P=0.004)、Ucp3(t=5.184，P=0.035)、Acsl6(t=6.213，P=0.003)、Slc2a4(t=19.300，P=0.005)具有統(tǒng)計(jì)學(xué)意義且與表2的分析結(jié)果相符，其中Pdk4、Hmgcs2、Ucp3在高糖刺激后表達(dá)增加，Acsl6、Slc2a4在高糖刺激后表達(dá)降低(圖5)。而Decr1(t=17.180，P<0.001)、Acot1(t=7.186，P=0.002)、Acot2(t=34.990，P<0.001)、Slc27a1(t=5.508，P=0.005)、Cpt1a(t=30.870，P<0.001)雖具有統(tǒng)計(jì)學(xué)意義但與表2基因芯片分析結(jié)果并不相符(圖5)。

討 論

通過(guò)對(duì)GEO數(shù)據(jù)庫(kù)中的兩個(gè)糖尿病心肌病數(shù)據(jù)集進(jìn)行分析后，共篩選出50個(gè)共同表達(dá)差異基因，對(duì)差異基因進(jìn)行GO富集后顯示，這些差異基因在脂肪酸代謝和脂肪酸修飾等生物功能中明顯富集，KEGG富集分析提示差異基因在脂肪酸延伸和不飽和脂肪酸生物合成通路中明顯富集。底物代謝紊亂是糖尿病心肌病的發(fā)病機(jī)制之一，在生理情況下，心臟可使用脂肪酸和葡萄糖作為主要的代謝底物，胰島素抵抗會(huì)導(dǎo)致葡萄糖攝取減少而游離脂肪酸的攝取和氧化增加[22]，基因富集的分析結(jié)果提示這些差異基因可能參與了糖尿病心肌病的底物代謝紊亂。在差異基因的KEGG富集分析中，PPAR信號(hào)通路在富集第1位，PPAR在心臟糖脂代謝以及能量穩(wěn)態(tài)中發(fā)揮關(guān)鍵作用，PPAR有多種亞型，其中PPARα在心臟中高表達(dá)，影響脂肪酸的攝取和線粒體脂肪酸的氧化[23]。PPARα不僅參與脂肪酸和葡萄糖的代謝，還參與支鏈氨基酸、酮體在糖尿病心肌病中的代謝[24]。差異基因在PPAR信號(hào)通路中的明顯富集表明這些差異基因可能通過(guò)PPAR信號(hào)通路參與脂肪酸的代謝，進(jìn)而參與心臟的底物代謝紊亂。GSEA富集分析提示糖尿病組在脂肪酸代謝和PPAR信號(hào)通路明顯富集。

表2 基因集中差異基因的表達(dá)情況

GO：基因本體論；KEGG：京都基因和基因組百科全書(shū)；BP：生物學(xué)過(guò)程；MF：分子功能；PPAR：過(guò)氧化物酶體增殖物激活受體

表3 MCC算法計(jì)算出的排名前10的關(guān)鍵基因的分值

實(shí)時(shí)熒光定量PCR的驗(yàn)證結(jié)果顯示，在高糖刺激下，Pdk4、Ucp3、Hmgcs2基因表達(dá)增加，Slc2a4、Acsl6基因表達(dá)下降。Pdk4是細(xì)胞能量代謝的關(guān)鍵激酶[25]，Pdk4特異性過(guò)表達(dá)的轉(zhuǎn)基因小鼠表現(xiàn)出胰島素抵抗和心肌葡萄糖氧化降低[26]，Ucp3是線粒體解偶聯(lián)蛋白家族成員，Ucp3與脂肪酸的β氧化直接相關(guān)[27]。在一項(xiàng)利用胰島β細(xì)胞遺傳缺陷的糖尿病心肌病的模型研究中，高糖低胰島素狀態(tài)會(huì)導(dǎo)致Pdk4和Ucp3的表達(dá)增加[28]，這與本研究結(jié)果一致，Pdk4和Ucp3在糖尿病心肌病細(xì)胞模型中表達(dá)升高，但是目前的研究尚無(wú)針對(duì)Pdk4和Ucp3在糖尿病心肌病中的功能研究。Hmgcs2是酮生成過(guò)程中的限速酶[29]，與本研究一致，在大鼠糖尿病模型中，心臟Hmgcs2表達(dá)明顯增加而肝臟中Hmgcs2無(wú)變化[30]。Slc2a4又被稱(chēng)作Glut4，Glut4是葡萄糖轉(zhuǎn)運(yùn)體最主要的亞型，占葡萄糖轉(zhuǎn)運(yùn)體的70%，在糖尿病心肌病中Glut4表達(dá)下降[31]。但是Glut4特異性過(guò)表達(dá)的糖尿病心肌病小鼠，非但不能挽救糖尿病心肌病的心臟代謝障礙，反而可能導(dǎo)致糖尿病心肌病的進(jìn)一步發(fā)展[32]，Glut4在糖尿病心肌病中的作用有賴(lài)于進(jìn)一步的研究。

Acsl6在糖尿病心肌病中尚無(wú)研究，在人類(lèi)和大鼠骨骼肌中Acsl6調(diào)節(jié)脂質(zhì)合成和線粒體的氧化能力[33]，Acsl6與小鼠大腦ω- 3脂肪酸DHA的富集相關(guān)[34]，也是小鼠神經(jīng)保護(hù)所必需的[35]。在Acsl的各種亞型中，Acsl3、Acsl4敲減的細(xì)胞系中，葡萄糖刺激胰島素釋放降低了50%，但是這項(xiàng)研究中Acsl6對(duì)胰島素的釋放并未產(chǎn)生影響[36]。在糖尿病心肌病的細(xì)胞模型中Acsl6基因表達(dá)降低，在本研究中是首次報(bào)道，Acsl6基因在糖尿病心肌病模型中的功能可作為糖尿病心肌病的一個(gè)研究方向。

ES：富集分?jǐn)?shù)；CTR：對(duì)照組；DCM：糖尿病心肌病組；NES：標(biāo)準(zhǔn)化富集分?jǐn)?shù)；FDR：錯(cuò)誤發(fā)現(xiàn)率

與現(xiàn)有的基于基因集GSE4745和GSE5606進(jìn)行差異基因分析的文章相比，Dai等[37]通過(guò)對(duì)基因集GSE5606進(jìn)行加權(quán)基因共表達(dá)網(wǎng)絡(luò)分析顯示，Angptl4、Acot1、Decr1、Hmgcs2、Pdk4在糖尿病心肌病的小鼠心臟中表達(dá)增加，這些基因同樣是本研究中的關(guān)鍵差異表達(dá)基因，但是與Dai等[37]的研究結(jié)果不同，Decr1和Acot1基因在本研究中表達(dá)下降，考慮為實(shí)驗(yàn)選擇動(dòng)物種屬不同。另一項(xiàng)關(guān)于GSE4745基因集的研究?jī)H對(duì)GSE4745基因集的差異表達(dá)基因進(jìn)行了分析，未在動(dòng)物或細(xì)胞水平檢測(cè)差異基因的表達(dá)水平[38]。與之前的研究相比，本研究同時(shí)使用了兩個(gè)基因集的數(shù)據(jù)，關(guān)鍵差異表達(dá)基因與Dai等[37]的研究雖有重合但仍有新的發(fā)現(xiàn)，Acsl6基因不僅在現(xiàn)有的基因芯片的分析中報(bào)道較少，在糖尿病心肌病的研究中報(bào)道也較少，且與Acsl家族中的Acsl3和Acsl4胰島素的釋放相關(guān)，這為以后的研究提供了很好的理論依據(jù)。另外，在現(xiàn)有的與GSE4745和GSE5606基因集研究中僅對(duì)差異基因進(jìn)行GO和KEGG富集不同，本研究同時(shí)進(jìn)行了基因芯片整體基因代謝通路的研究，是對(duì)GO和KEGG僅對(duì)差異基因進(jìn)行富集分析的補(bǔ)充。

MCC：最大集團(tuán)中心

CTR：對(duì)照組；HG：高糖組

本研究雖篩選出了糖尿病心肌病中的關(guān)鍵基因并篩選出從未報(bào)道過(guò)的基因Acsl6，但是對(duì)篩選出基因的功能以及與相關(guān)通路的聯(lián)系并未進(jìn)行進(jìn)一步的研究，這是本研究的局限性。Pdk4、Ucp3、Hmgcs2、Slc2a4、Acsl6這些在本研究中篩選出的關(guān)鍵基因，可為以后糖尿病心肌病的研究提供方向。