吳思騏 劉建博

摘要 目的：通過生物信息學技術分析比較慢性阻塞性肺疾?。–OPD）患者與健康人的基因芯片數(shù)據(jù)，初步篩選COPD的差異表達基因，并借助差異表達基因預測治療COPD的潛在中藥藥物。方法：首先通過R分析GSE27597和GSE106986基因數(shù)據(jù)，得到139個差異表達基因，對差異表達基因進行基因本體（GO）富集分析，并通過String數(shù)據(jù)庫、Cytoscape軟件及其插件分析得到14個核心基因，然后對核心基因進行Reactome、KEGG、PANTHER等通路注釋，最后通過核心基因與醫(yī)學本體信息檢索平臺（Coremine Medical）相互映射，篩選治療COPD的中藥。結果：通過MCODE篩選出14個核心基因包括MMP8、PTX3、DEFA4、CHRM2、FPR2、CX3CR1、FPR1、ORM1、APLNR、PPBP、CHIT1、AGTR2、ARG1，與中性粒細胞脫顆粒、先天免疫系統(tǒng)、趨化因子信號通路等密切相關，其中FPR1和CX3CR1通過影響中性粒細胞向肺部的遷移活動成為COPD病理變化的關鍵。治療COPD的潛在中藥藥物包括黃連、黃芪、黃芩、半夏、丹參、厚樸、姜黃、郁金、瓜蔞等。結論：部分基因會促進香煙煙霧暴露誘發(fā)COPD的發(fā)生，通過miRNA進行調節(jié)可能是中藥治療COPD的潛在作用機制。

關鍵詞 差異表達基因;慢性阻塞性肺疾病;微小RNA;中藥;香煙煙霧

Abstract Objective：By analyzing and comparing the gene chip data of chronic obstructive pulmonary disease （COPD） patients and healthy people through bioinformatics technology，the differentially expressed genes of COPD are initially screened，and the potential Chinese medicinal drugs for the treatment of COPD are predicted with the help of differentially expressed genes.Methods：The gene expression profiles of lung tissues of COPD cases and healthy cases were obtained from GSE27597 and GSE106986 of the Gene Expression Omnibus.A otal of 139 differentially expressed genes were obtained，and we performed gene ontology （GO） enrichment analysis of differentially expressed genes，and obtained 14 core genes through String database，Cytoscape software and its plug-in analysis，and then we performed Reactome，KEGG，PANTHER and other pathways on the core genes.Finally through the mutual mapping between core genes and the medical ontology information retrieval platform （Coremine Medical），we screened Chinese medicinal for the treatment of COPD.Results：A total of 14 core genes screened by MCODE including MMP8，PTX3，DEFA4，CHRM2，F(xiàn)PR2，CX3CR1，F(xiàn)PR1，ORM1，APLNR，PPBP，CHIT1，AGTR2，ARG1，and neutrophil degranulation，innate immune system，chemokines signal pathways were closely related，among which FPR1 and CX3CR1 have become the key to the pathological changes of COPD by affecting the migration of neutrophils to the lungs.Potential Chinese medicinals for the treatment of COPD include Rhizoma Coptidis，Radix Astragali seu Hedysari，Radix Scutellariae，Rhizoma Pinelliae，Radix Salviae Miltiorrhizae，Cortex Magnoliae Officinalis，Rhizoma Curcumae Longae，Radix Curcumae，F(xiàn)ructus Trichosanthis.Conclusion：Some genes can promote the occurrence of COPD induced by cigarette smoke exposure，and regulation by miRNA may be a potential mechanism of Chinese medicinal in the treatment of COPD.

Keywords Differentially expressed gene; Chronic obstructive pulmonary disease; miRNA; Chinese medicinal; Cigarette smoke

中圖分類號：R282;R563文獻標識碼：Adoi：10.3969/j.issn.1673-7202.2021.14.023

慢性阻塞性肺疾?。–hronic Obstructive Pulmoriary Disease，COPD）是以持續(xù)性呼吸道癥狀和氣流受限為主要特征的一種呼吸系統(tǒng)常見病，感染是誘發(fā)COPD急性加重的主要原因，香煙煙霧則是COPD最重要的獨立危險因素[1]。COPD是全球第4大死因，根據(jù)2020年COPD全球創(chuàng)議（Global Initiative for Chronic Obstructive Lung Disease）可知，隨著發(fā)展中國家吸煙率的升高和高收入國家老齡化加劇，COPD的發(fā)病率在未來40年仍會繼續(xù)上升，至2060年可能每年有超過540萬人死于COPD及其相關疾病[2]。目前COPD的基礎治療用藥為吸入性支氣管擴張劑，而在COPD的防治過程中，戒煙的重要性越來越得到認可。

盡管吸煙是COPD的重要誘因，但并非所有吸煙者都必然發(fā)展為COPD患者。部分COPD患者無明顯的香煙煙霧暴露史，也發(fā)展為重癥COPD，可見吸煙與COPD之間存在著復雜的內在關系。有研究表明，COPD具有一定的遺傳性，但其具體機制尚不清楚。因此，本研究著眼于健康吸煙者與COPD患者、健康非吸煙者與COPD患者的差異表達基因，通過生物信息學分析，嘗試從基因水平探尋COPD的發(fā)病因素。同時通過差異表達基因預測防治COPD的中藥，并從活性單體、蛋白分子、miRNA等多個層面探討其作用機制，從發(fā)病原因到治療方法對COPD進行全面的剖析。

1 資料與方法

1.1 數(shù)據(jù)提取 從國家生物技術信息中心基因表達綜合數(shù)據(jù)庫（https：//www.ncbi.nlm.nih.gov/geo/）下載了GSE27597和GSE106986 2個數(shù)據(jù)集。其中GSE27597包含了64例COPD患者的肺組織樣本和16個來自健康供體的樣本的mRNA表達譜[3]，GSE106986包含了14例COPD患者的肺組織樣本和5個來自健康供體的樣本的mRNA表達譜[4]。二者的表達譜來自Affymetrix Human Exon 1.0 ST GeneChip和Agilent-026652 Whole Human Genome Microarray 4x44K v2平臺。

1.2 差異基因的篩選 使用R軟件的注釋包對2個數(shù)據(jù)集的探針I(yè)D進行轉換并將二者集合并為一個數(shù)據(jù)集，利用R包（sva包和limma包）對2個數(shù)據(jù)集進行批量歸一化。通過R 3.60中的limma包篩選COPD病例與對照組之間的差異表達基因，以log2 fold change（簡稱log2FC）的絕對值>1.0和P<0.05作為篩選標準，篩選出差異表達基因（Differentially Expressed Gene，DEG）。其中，log2FC為負數(shù)的基因為下調基因，log2FC為正數(shù)的基因為上調基因。利用R語言制作差異表達基因火山圖。

1.3 差異基因的功能注釋 對差異基因進行多個層面的功能注釋，包括基因本體（GO）富集分析和生物過程富集分析?；虮倔w（GO）富集分析包含以下3個方面：細胞組分（CC）、分子功能（MF）和生物過程（BP）。借助R語言的clusterProfiler包和ggplot2包進行基因本體（GO）富集分析和京都基因和基因組百科全書（KEGG）富集分析和可視化。同時，通過Funrich軟件（Funrich 2.1.2）尋找差異基因富集的生物通路。

1.4 核心基因的篩選和功能注釋 使用STRING數(shù)據(jù)庫（https：//string-db.org/）分析顯著差異表達基因的蛋白質-蛋白質相互作用（PPI）網(wǎng)絡，并將得到的PPI網(wǎng)絡導入Cytoscape軟件（version 3.8.0），利用CytoHubba[5]插件和MCODE插件進行PPI分析（篩選標準：combined score>0.4）和網(wǎng)絡分析。MCODE可以發(fā)現(xiàn)PPI網(wǎng)絡中緊密聯(lián)系的區(qū)域，同時兼顧到網(wǎng)絡邊緣的分子，避免了因高通量技術帶來假陽性影響，更有利于挖掘與特定功能相關的分子集合。通過KOBAS（http：//kobas.cbi.pku.edu.cn/）進行通路富集，尋找核心基因富集的通路，并通過string構建核心基因之間的交互作用網(wǎng)絡。

1.5 靶點對接及中藥預測 將篩選出的核心基因與Coremine Medical（http：//www.coremine.com/）相互映射，篩選出核心基因相關的治療中藥，根據(jù)P<0.05篩選出差異有統(tǒng)計學意義的中藥。

2 結果

2.1 差異基因 基于GSE27597和GSE106986的原始數(shù)據(jù)，以log2FC的絕對值>1.0和P<0.05為篩選標準，共得出139個差異表達基因，其中包括62個上調基因和77個下調基因。見表1，圖1。

2.2 差異基因的功能注釋 GO發(fā)現(xiàn)上調的差異表達基因主要與中性粒細胞活化、中性粒細胞脫顆粒、中性粒細胞活化參與免疫反應、中性粒細胞介導的免疫應答等生物過程有關，它們介導的分子則富集于以下生物過程：碳水化合物結合、肽結合、細胞因子受體活性、G蛋白偶聯(lián)肽受體活性等。見圖2。在生物通路方面，差異基因富集于肽配體結合受體、整合素細胞表面相互作用、鈣調神經(jīng)磷酸酶調節(jié)淋巴細胞中NFAT依賴性轉錄等通路。見圖3。且較之其他部位，這些差異基因在肺組織的表達程度更高。見圖4。為差異基因與COPD的相關性提供佐證。

2.3 核心基因的篩選及功能注釋 首先，以required confidence（combined score）>0.4為PPI關系的閾值，利用String數(shù)據(jù)庫進行PPI分析，并下載相關TSV格式文件，得到PPI關系對文件。接著，使用Cytoscape軟件的CytoHubba插件對其進行下一步分析，發(fā)掘復雜網(wǎng)絡中的關鍵目標和子網(wǎng)絡;然后使用Cytoscape軟件的MCODE插件對其進行模塊分析。CytoHubba分析依據(jù)degree值找到14個表達量最高的基因（CX3CR1、PPBP、PTGS2、FPR1、FPR2、VCAM1、S100A12、ARG1、EGR1、CD163、FGG、ORM1、S100A8，S100A9），MCODE分析則得到3個重要模塊（MMP8、PTX3、DEFA4、CHRM2、FPR2、CX3CR1、FPR1、ORM1、APLNR、PPBP、CHIT1、AGTR2、ARG1;S100A9、VNN2、S100A12、AQP9、CEACAM4、S100A8;NR4A3、ATF3、EGR2、EGR3、EGR1、FOSB）。Hubba篩選出的模塊1包含了Cytohubba篩選出的14個基因中表達量最高的4個基因，且富集于特定的功能，同時囊括了上調表達基因與下調表達基因，更能全面反映差異表達基因對COPD的影響，因此選擇MOCDE篩選出的模塊1所包含的14個基因作為核心基因。見圖5～7。

通過KOBAS（http：//kobas.cbi.pku.edu.cn/）對篩選出的14個核心基因進行通路富集分析，其中涵蓋了Reactome、KEGG、PANTHER等多個數(shù)據(jù)庫。根據(jù)Reactome數(shù)據(jù)庫，14個核心基因主要富集于中性粒細胞脫顆粒、先天免疫系統(tǒng)、免疫系統(tǒng)，KEGG通路方面，14個核心基因則富集于神經(jīng)活性配體-受體相互作用、金黃色葡萄球菌感染、病毒蛋白與細胞因子和細胞因子受體的相互作用、趨化因子信號通路、細胞因子與細胞因子受體的相互作用等通路。趨化因子和細胞因子信號轉導途徑介導的炎癥是PANTHER預測的核心基因富集的生物過程。見表2。

將14個核心基因映射到CoremineMedical數(shù)據(jù)庫，以P<0.05為篩選標準，篩選治療COPD潛在中藥，共篩選出35種可能具有防治COPD的中藥。見表3。

3 討論

3.1 COPD與健康人之間存在差異基因 盡管吸煙是COPD發(fā)病最重要的危險因素，但吸煙與COPD并非簡單的線性關系。本研究綜合分析了COPD與健康吸煙者以及COPD與健康非吸煙者的DEG，表明COPD的易感性與基因有關。顯然，COPD并非單基因疾病，其發(fā)病與否受到多種基因的共同調節(jié)，且部分差異基因之間存在明顯的相互作用。

3.2 核心基因與COPD存在著復雜關系 本研究篩選出14個核心基因，并提示CX3CR1、APLNR、CHIT1、AGTR2的高表達與COPD正相關，MMP8、PTX3、DEFA4、FPR2、FPR1、ORM1、PPBP、ARG1的低表達與COPD負相關。它們通過多種通路和生物過程，對COPD的誘發(fā)、病情進展、疾病預后起到重要的影響。而根據(jù)Cytohubba結果，CX3CR1、FPR1、FPR2、ORM1、PPBP、ARG1表達量最高，與COPD的發(fā)病關系也最為密切。

CX3CR1具有促進結構破壞和重塑的作用，在慢性阻塞性肺疾?。–OPD）和肺氣腫等慢性炎癥性肺部疾病的發(fā)生發(fā)展中發(fā)揮重要作用。釋放CX3CR1的趨化性炎癥細胞可能導致單個核細胞在肺實質和肺血管壁聚集，釋放介質誘導損傷，刺激增殖和趨化性炎癥細胞[6]。單核巨噬細胞通過產(chǎn)生TNF-α和IL-6對香煙煙霧產(chǎn)生先天免疫應答，從而促進肺氣腫的發(fā)生[7]，而CX3CR1作為單核巨噬細胞的產(chǎn)物，在這個過程中不可或缺。

FPR1的激活可觸發(fā)中心粒細胞的多種功能，包括趨化性、脫顆粒、產(chǎn)生ROS和吞噬作用，并參與了香煙煙霧誘導的COPD炎癥的發(fā)生。Leslie等[8]表明FPR1基因的缺乏降低了肺部發(fā)生炎癥和纖維化的風險，同時也證明了FPR1基因表達上調對中性粒細胞產(chǎn)生重要影響，使之能夠在肺纖維化中發(fā)揮特異性的關鍵作用，對于患病肺的病理組織重塑至關重要。另有研究表明，COPD患者肺中的中性粒細胞顯著增多，并且中性粒細胞相關的產(chǎn)物與COPD的發(fā)生和嚴重程度相關[9]。FPR1與香煙煙霧暴露下中性粒細胞的肺部轉移程度負相關也在動物實驗中得到驗證[10]。纖維化和肺部炎癥是COPD的重要表現(xiàn)，香煙煙霧會導致中性粒細胞向肺內遷移從而發(fā)生炎癥反應，而FPR1基因下調減少了中性粒細胞向肺的遷移，這為闡明香煙煙霧與COPD的個體差異性提供了一種可能解釋。

ORM1編碼的血漿蛋白在急性炎癥顯著增加，被歸類為急性期反應物。有研究通表明AECOPD患者的ORM1的基因和蛋白質表達顯著高于穩(wěn)定COPD或健康對照者[11]。

PPBP能刺激滑膜細胞的纖溶酶原激活物的形成和分泌，是一種強有力的中性粒細胞趨化劑和激活劑，其在支氣管黏膜的高表達可能參與了COPD的發(fā)病[11]。ARG1、FPR2與COPD的發(fā)病也有一定關系。

FPR1和CX3CR1與COPD的發(fā)病機制與病情進展密切相關，并得到深入研究，它們影響了香煙煙霧暴露情況下中性粒細胞的活動，其表達量的差異使得個體在COPD的發(fā)病和進展中呈現(xiàn)出不一樣的轉歸，為揭示COPD易感性與差異性提供了遺傳水平的解釋。

3.3 相關中藥防治COPD的潛在機制探討

通過Coremine Medical篩選出的中藥包括了臨床中治療COPD的常用藥物。實驗表明，相關中藥及其活性分子與對COPD的治療有積極的作用，同時它們的作用機制也與對應的核心基因密切相關。此外，部分中藥可能通過影響COPD相關的miRNA而對疾病的進展產(chǎn)生影響，并得到了實驗證明。

3.3.1 中藥的活性成分治療COPD的分子機制 ? 有研究表明，半夏生物堿可能通過抑制TNF-α和NO等早期炎癥介質的釋放而保護損傷的肺上皮細胞，同時抑制IL-8和ICAM-1的表達、緩解IL-8和ICAM-1過表達而觸發(fā)的中性粒細胞趨化作用，減緩中性粒細胞聚集，延緩肺損傷的發(fā)展[12]。而中性粒細胞聚集正是差異基因GO富集分析的主要通路之一。厚樸酚能通過抗氧化及抗炎作用抑制膿毒癥大鼠急性肺損傷[13]，同時可通過減輕肺纖維化大鼠的肺組織纖維化程度，降低肺組織含量，達到抗肺纖維化作用[14]。另有研究表明，厚樸揮發(fā)油可能通過降低炎癥組織PGE2、TNF-α和IL-1β水平從而發(fā)揮抗炎作用[15]。有研究表明，黃芪多糖可通過抑制COPD大鼠肺組織MMP9、TIMP1的表達、改善肺組織損傷及肺功能[16]。

根據(jù)中藥的功能分類，預測的治療COPD潛在靶點的藥物大致可分為以下5類：止咳化痰平喘藥（半夏、瓜蔞）、清熱藥（黃芩、大黃、黃連）、補虛藥（黃芪、干姜、吳茱萸）、活血化瘀藥（丹參、姜黃、郁金）、理氣藥（厚樸）。其中黃芪補益肺脾之氣，臨床研究表明，COPD患者大多有肺脾氣虛，黃芪也是臨床上最常應用于COPD治療的中藥之一[17]。肺脾腎虛是COPD患者最主要的病機，急性加重期患者往往受到痰熱、痰濕、瘀血等病理因素的影響，預測的中藥也涵蓋了多種中醫(yī)病理因素，體現(xiàn)了中醫(yī)“標本兼治”的思想。瓜蔞清熱滌痰，半夏燥濕化痰，而痰濕是COPD主要的病理因素及病理產(chǎn)物，二者聯(lián)用可有效減輕COPD的癥狀，對COPD起到治療作用。肺主一身之氣，主宣發(fā)肅降，肺氣受損，氣機升降失司，則發(fā)為喘，厚樸燥濕消痰、下氣除滿，在COPD的治療中得到廣泛應用，療效顯著。黃芩、大黃、黃連、干姜、吳茱萸、丹參、姜黃、郁金對COPD的療效也得到臨床研究的證實[18-21]。

3.3.2 miRNA可能是中藥干預COPD進程的關鍵 ?微小RNA（microRNA，miRNA）是在多種真核細胞及病毒中發(fā)現(xiàn)的一類參與基因表達調控的非編碼單鏈RNA。近期研究表明，植物來源的miRNA可以穩(wěn)定地存在于動物的血液與組織器官中、并參與調控多種靶點蛋白的表達而發(fā)揮作用。尹作靜等[22]發(fā)現(xiàn)黃連素和吳茱萸堿作用于肝癌細胞后hsa-miR-188-5p的表達量升高，而hsa-miR-188-5p與CX3CR1相關聯(lián)，由此推測，黃連素可能通過上調hsa-miR-188-5p的表達干預CX3CR1的表達，從而影響COPD的進程。研究表明，丹參酮ⅡA的抗炎功效與多種miRNAs密切相關，通過下調miR-155、miR-147、miR-184、miR-29b和miR-34c表達，從而降低細胞因子、趨化因子和急性期GM-CSF、sICAM-1、cxcl-1、MIP-1α水平，抑制核因子κB活化，降低TNF-α、IL-1β、COX-2的水平，這些因子與COPD的發(fā)病及進展有著密切的關系[23]。彭先祝[24]發(fā)現(xiàn)參芪補肺方通過介導miRNA-146a對煙霧溶液誘導COPD效應細胞，正向調控Nrf2表達，負向調控核因子κB表達。黃芪總黃酮可通過抑制miRNA-21，增加Smad7表達，使TGF-β1/Smad激活能力下降，抑制肺纖維化中上皮細胞間質轉分化的作用[25]。有研究表明，吸煙者和COPD患者的血清中miRNA-21的水平較高，并且與第1秒用力呼氣容積比用力肺活量（FEV1/FVC）的呼氣量負相關[26]。而miRNA-21能夠促進肺成纖維細胞分化為肌成纖維細胞，可能是COPD氣道重塑的重要機制[27]。

研究表明，部分miRNA與COPD密切相關，可作為診斷和預后的生物標志物，miRNA作為COPD治療靶點更具有廣泛的應用前景[28]。部分中藥可能通過調控miRNA影響COPD的病理變化，進而對COPD的發(fā)生和發(fā)展產(chǎn)生影響。從miRNA的層面對中藥進行分析，為COPD的研究提供了新的思路，為理解COPD發(fā)生過程中的復雜分子機制提供了新的視角，也為進一步推動臨床藥物研發(fā)提供了依據(jù)。然而，考慮到本研究尚有一定的局限性，未來還需要進一步開展更為深入的實驗驗證。

參考文獻

[1]Wang C，Xu J，Yang L，et al.Prevalence and risk factors of chronic obstructive pulmonary disease in China（the China Pulmonary Health[CPH]study）：a national cross-sectional study[J].Lancet，2018，391（10131）：1706-1717.

[2]Global strategy for the diagnosis，management and prevention of chronic obstructive pulmonary disease 2020 report[EB/OL].（2019-11-30）.https：//goldcopd.org/gold-reports/.

3]Campbell JD，McDonough JE，Zeskind JE，et al.A gene expression signature of emphysema-related lung destruction and its reversal by the tripeptide GHK[J].Genome Med，2012，4（8）：67.

[4]Heinbockel L，Marwitz S，Schromm AB，et al.Identification of novel target genes in human lung tissue involved in chronic obstructive pulmonary disease[J].Int J Chron Obstruct Pulmon Dis，2018，13：2255-2259.

[5]Chin CH，Chen SH，Wu HH，et al.cytoHubba：identifying hub objects and sub-networks from complex interactome[J].BMC Syst Biol，2014，8 Suppl 4：S11.

[6]Zhang J，Patel JM.Role of the CX3CL1-CX3CR1 axis in chronic inflammatory lung diseases[J].Int J Clin Exp Med，2010，3（3）：233-244.

[7]Xiong Z，Leme AS，Ray P，et al.CX3CR1+ lung mononuclear phagocytes spatially confined to the interstitium produce TNF-α and IL-6 and promote cigarette smoke-induced emphysema[J].J Immunol，2011，186（5）：3206-3214.

[8]Leslie J，Millar BJ，Del Carpio Pons A，et al.FPR-1 is an important regulator of neutrophil recruitment and a tissue-specific driver of pulmonary fibrosis[J].JCI Insight，2020，5（4）：e125937.

[9]Dicker AJ，Crichton ML，Pumphrey EG，et al.Neutrophil extracellular traps are associated with disease severity and microbiota diversity in patients with chronic obstructive pulmonary disease[J].J Allergy Clin Immunol，2018，141（1）：117-127.

[10]Cardini S，Dalli J，F(xiàn)ineschi S，et al.Genetic ablation of the fpr1 gene confers protection from smoking-induced lung emphysema in mice[J].Am J Respir Cell Mol Biol，2012，47（3）：332-339.

[11]Shi L，Zhu B，Xu M，et al.Selection of AECOPD-specific immunomodulatory biomarkers by integrating genomics and proteomics with clinical informatics[J].Cell Biol Toxicol，2018，34（2）：109-123.

[12]吳偉斌，祝春燕，羅超.半夏生物堿對肺上皮細胞炎癥損傷的保護作用研究[J].內蒙古農(nóng)業(yè)大學學報：自然科學版，2018，39（4）：1-8.

[13]段金旗，林艷，趙芹芳，等.厚樸酚對膿毒癥大鼠模型急性肺損傷的保護作用及機制[J].中國老年學雜志，2018，38（11）：2745-2747.

[14]劉波，莫鎮(zhèn)濤，李意奇，等.厚樸酚對肺纖維化大鼠肺病理學和羥脯氨酸含量的影響[J].中國老年學雜志，2017，37（6）：1325-1327.

[15]曹迪，徐照輝，王芳芳.厚樸揮發(fā)油化學成分及其抗炎作用的實驗研究[J].中國中醫(yī)藥科技，2015，22（6）：647-649.

[16]余維巍，黃驍燕，張艷，等.黃芪多糖對慢性阻塞性肺疾病大鼠肺組織內基質金屬蛋白酶-9、金屬基質蛋白酶抑制劑-1表達的影響[J].中國醫(yī)藥導報，2012，9（2）：25-27.

[17]李國棟，牛潔，吳志松，等.近20年中藥治療慢性阻塞性肺疾病穩(wěn)定期組方用藥規(guī)律文獻研究[J].中醫(yī)雜志，2019，60（7）：614-618.

[18]鄭文江，彭紫荊，陳淑婉，等.三子養(yǎng)親湯加味治療慢性阻塞性肺疾病急性加重期的系統(tǒng)評價[J].中國中藥雜志，2019，44（10）：2171-2178.

[19]鄭文江，彭紫荊，嚴倩，等.桑白皮湯加減治療慢性阻塞性肺疾病急性加重期的Meta分析和試驗序貫分析[J].中國中藥雜志，2019，44（17）：3806-3815.

[20]張葵，張櫻，陳翌江，等.參芪補肺湯對肺氣虛證慢性阻塞性肺疾病大鼠氣道重構中MMP-9和TIMP-1蛋白表達的影響[J].中國中醫(yī)藥信息雜志，2008，15（8）：40-42.

[21]謝文英，尚立芝，胡文豪，等.慢性阻塞性肺疾病的發(fā)病機制及中醫(yī)藥治療進展[J].中國實驗方劑學雜志，2015，21（9）：227-230.

[22]尹作靜，曹志偉，閆鑫淼，等.黃連素和吳茱萸堿協(xié)同抗癌的miRNA網(wǎng)絡機制研究[J].中國藥理學通報，2017，33（6）：772-780.

[23]Fan G，Jiang X，Wu X，et al.Anti-Inflammatory Activity of Tanshinone ⅡA in LPS-Stimulated RAW264.7 Macrophages via miRNAs and TLR4-NF-κB Pathway[J].Inflammation，2016，39（1）：375-384.

[24]彭先祝.參芪補肺方通過miRNA-146a調控COPD大鼠模型Nrf2和NF-κB的機理研究[D].哈爾濱：黑龍江中醫(yī)藥大學，2019.

[25]馬林納，王瑞，苗明三.基于數(shù)據(jù)挖掘的治療肺纖維化中藥用藥規(guī)律分析[J].中草藥，2020，51（6）：1406-1411.

[26]Gon Y，Shimizu T，Mizumura K，et al.Molecular techniques for respiratory diseases：MicroRNA and extracellular vesicles[J].Respirology，2020，25（2）：149-160.

[27]Li Y，Yin Z，F(xiàn)an J，et al.The roles of exosomal miRNAs and lncRNAs in lung diseases[J].Signal Transduct Target Ther，2019，4：47.

[28]Mohan A，Agarwal S，Clauss M，et al.Extracellular vesicles：novel communicators in lung diseases[J].Respir Res，2020，21（1）：175.

（2020-10-26收稿 責任編輯：張雄杰）