徐 彤 ，馬巖團(tuán)進(jìn) ，張宇航 ，何秋月 ，黃 薇 ，呂夢欣 ，唐 健 ，何建萍 ，蔣國慶 ，錢 源 1,2,3,4,5,

（1）云南省檢驗醫(yī)學(xué)重點(diǎn)實驗室，云南 昆明 650032;2）云南省醫(yī)學(xué)檢驗臨床醫(yī)學(xué)研究中心，云南昆明 650032;3）昆明醫(yī)科大學(xué)第一附屬醫(yī)院醫(yī)學(xué)檢驗科，云南 昆明 650032;4）昆明市婦幼保健院醫(yī)學(xué)遺傳與產(chǎn)前診斷中心，云南 昆明 650032;5）昆明市婦幼保健院國家婦產(chǎn)疾病臨床醫(yī)學(xué)研究中心（云南）聯(lián)合重點(diǎn)實驗室，云南 昆明 650032;6）昆明醫(yī)科大學(xué)第一附屬醫(yī)院產(chǎn)科，云南 昆明 650032）

維生素 D 是人體內(nèi)不可缺少的類固醇激素，也是一種脂溶性維生素，與血管生成、調(diào)節(jié)組織器官免疫、細(xì)胞的凋亡和增殖等方面密切相關(guān)[1]。維生素 D 可以在蛻膜和胎盤組織中表達(dá)，在整個胚胎發(fā)育過程當(dāng)中起到非常重要的作用[2]。即使維生素 D 在人體生命過程中發(fā)揮著十分重要的作用，但維生素 D 在孕婦人群中的缺乏仍是一個普遍現(xiàn)象[3]。有學(xué)者評價了來自7 個國家的13 項研究，發(fā)現(xiàn)孕婦維生素D 缺乏的患病率在51.3%至100%之間，而亞洲、中東和非洲這些地區(qū)的孕婦的維生素D 缺乏患病率在全球范圍內(nèi)最高，其值高達(dá)70%～80%[4]。妊娠與維生素 D 的代謝有著明顯相關(guān)性，維生素D 的缺乏與產(chǎn)后抑郁、胎兒骨骼發(fā)育、胎兒骨骼發(fā)育、早產(chǎn)、流產(chǎn)、妊娠期糖尿病、小于胎齡兒、子癇前期、不孕不育等圍孕期疾病密切相關(guān)[5-8]。

在本研究中，筆者通過對維生素 D 刺激后胎盤滋養(yǎng)細(xì)胞的RNA 測序分析，尋找DEG，為維生素 D 缺乏引起的不良妊娠結(jié)局的診治提供參考。

1 材料與方法

1.1 細(xì)胞培養(yǎng)與收集

人胎盤滋養(yǎng)細(xì)胞株 HTR-8/SVneo 購自BioVector 質(zhì)粒載體菌種細(xì)胞基因保藏中心。將1，25（OH）2D3粉末用無水乙醇配制成100 μmol/L 的原液，再用1640 培養(yǎng)基分別稀釋成0.1 nmol/L、1 nmol/L、10 nmol/L、100 nmol/L 的稀釋液。將滋養(yǎng)細(xì)胞分為 5 組，包括不用維生素 D 刺激的對照組（Control 組）、維生素 D 刺激濃度為0.1 nmol/L組（VitD 0.1 組）、維生素D 刺激濃度為1 nmol/L組（VitD 1 組）、維生素 D 刺激濃度為 10 nmol/L組（VitD 10 組）和維生素 D 刺激濃度為100 nmol/L 組（VitD 100 組）。于37 ℃、5% CO2培養(yǎng)箱中，每2～3 d 對細(xì)胞進(jìn)行傳代培養(yǎng)，取處于對數(shù)生長期的細(xì)胞進(jìn)行下一步的刺激實驗。用不同濃度的維生素 D（0.1、1、10、100 nmol/L）刺激細(xì)胞，于24 h 后收取并放入-80 ℃冰箱保存，便于細(xì)胞的RNA 提取和測序分析。

1.2 總 RNA 提取、質(zhì)檢和測序分析

建庫起始提取細(xì)胞的總 RNA，總量 ＞=1 μg。使 用 Agilent 對 RNA 進(jìn)行長度檢測，樣品檢測達(dá)到合格標(biāo)準(zhǔn)后，構(gòu)建文庫，Illumina Hiseq 平臺上測序。

1.3 DEG 的篩選

使用R 語言的 limma 包篩選出DEG。設(shè)置P值和差異倍數(shù)（fold change，F(xiàn)C）進(jìn)行DEG 的篩選，當(dāng)P值 ＜ 0.05 和|log2FC| ＞ 1 被認(rèn)為差異具有統(tǒng)計學(xué)意義。使用“Venny”（http://bioinfogp.cnb.csic.es/tools/venny/index.html）在線工具繪制EDG 的交集韋恩圖。

1.4 GO 和 KEGG 對 DEG 進(jìn)行富集分析

DAVID（http://david.ncifcrf.gov/）在線工具對篩選出的DEG 進(jìn)行GO 和KEGG 富集分析。

1.5 PPI 網(wǎng)絡(luò)構(gòu)建與中樞基因篩選

利用STRING（https://string-db.org）識別DEG編碼蛋白質(zhì)之間的相互作用。用Cytoscape 軟件（V.3.8.0）可視化 PPI 網(wǎng)絡(luò)，使用CytoHubba 中的對Hub 基因進(jìn)行排序和評估。

2 結(jié)果

2.1 DEG 的篩選結(jié)果

分析了4 種不同濃度維生素D 刺激的細(xì)胞樣本和未加維生素D 刺激的正常對照樣本的基因表達(dá)譜數(shù)據(jù)（圖2）。通過對表達(dá)矩陣數(shù)據(jù)分析，VitD 0.1 vs Control 組共獲取了 354 個DEG，包括242 個上調(diào)，112 個下 調(diào)（圖1A）；VitD 1 vs Control組共獲取了 320 個DEG，包括 216 個上調(diào)，104個下調(diào)（圖1B）；VitD 10 vs Control 組共獲取了374 個DEG，包括 277 個上調(diào)，97 個下調(diào)（圖1C）；VitD 100 vs Control 組共獲取了 300 個DEG，包括151 個上調(diào)，149 個下調(diào)（圖1D）。最后對 4 組差異基因數(shù)據(jù)取交集，共獲得 9 個共同的DEG（圖3，表1）。

表1 共同差異基因表達(dá)情況Tab.1 Expression of common DEGs

圖1 滋養(yǎng)細(xì)胞差異表達(dá)基因火山圖Fig.1 Volcano map of differentially expressed genes in trophoblast cells

圖2 滋養(yǎng)細(xì)胞基因表達(dá)的熱圖Fig.2 Heat map of trophoblast gene expression

圖3 差異表達(dá)基因韋恩圖Fig.3 Venn diagram of differentially expressed genes

2.2 GO 富集分析

對 4 組DEG 進(jìn)行GO 的富集分析，包括生物學(xué)過程（BP）、細(xì)胞組成（CC）、和分子功能（MF）。VitD 0.1 組vs Control 組DEG 顯著富集的 BP 主要包括膠原原纖維組織、細(xì)胞分化和趨化作用，CC主要包括膠原三聚體復(fù)合物、蛋白質(zhì)的細(xì)胞外基質(zhì)和膠原蛋白三聚物，MF 主要包括受體配體活動、受體調(diào)節(jié)器活動和蛋白質(zhì)橋接（圖4 A）。VitD 1 組vs Control 組DEG 顯著富集的BP 包括蛋白質(zhì)分泌的正調(diào)控和蛋白質(zhì)分泌等，CC 包括膠原三聚體復(fù)合物和膠原蛋白三聚物，MF 主要包括受體調(diào)節(jié)器活動和細(xì)胞因子活性（圖4 B）。VitD10 組vs Control 組DEG 顯著富集的 BP 主要包括細(xì)胞分化的負(fù)調(diào)控、信號受體活性的調(diào)節(jié)等，CC 主要包括蛋白質(zhì)的細(xì)胞外基質(zhì)和細(xì)胞的頂端部分，MF 主要包括受體配體活動和受體調(diào)節(jié)器活動（圖4 C）。VitD 100 組vs Control 組DEG 顯著富集的BP 包括細(xì)胞外結(jié)構(gòu)組織和細(xì)胞外基質(zhì)的組織，CC 主要包括膠原蛋白三聚物、膠原三聚體復(fù)合物和細(xì)胞外基質(zhì)成分，MF 主要包括受體配體活動、受體調(diào)節(jié)器活動（圖4D）。

圖4 GO 富集分析結(jié)果Fig.4 GO enrichment analysis results

2.3 KEGG 通路分析

KEGG 分析顯示VitD 0.1 組 vs Control 組DEG 主要富集于刺激神經(jīng)組織的受體配體相互作用、細(xì)胞粘附分子和胃癌等信號通路（圖5A），VitD 1 組 vs Control 組DEG 主要富集于細(xì)胞因子受體相互作用和產(chǎn)生IgA 的腸道免疫網(wǎng)絡(luò)等信號通路（圖5 B），VitD10 組 vs Control 組DEG 主要富集于炎癥性腸病和瘧疾等相關(guān)信號通路（圖5 C），VitD 100 組 vs Control 組DEG 主要富集于炎癥性腸病、阿米巴病和 PI3K-Akt 信號通路（圖5 D）。

圖5 KEGG 通路分析結(jié)果Fig.5 Results of KEGG pathway analysis

2.4 PPI 分析

STRING 分析4 組DEG 構(gòu)建PPI 網(wǎng)絡(luò)，使用Cytoscape 軟件進(jìn)行網(wǎng)絡(luò)可視化，用Cytoscape 軟件中Cytohubba 插件中的 “Degree ”來計算DEG的重要性排名。如 圖6A 和表2 所示VitD 0.1 組vs Control 組DEG 的Degree 排名前5 的分別是COL1A2、ACTA2、S100A4、TAGLN 和CSF1R。如圖6B 和表3 所示VitD 1 組 vs Control 組DEG的Degree 排名前5 的分別是 TLR4、TNFSF13B、FTCD、S100A4 和 APOBEC3G。如圖6C 和表4 所示VitD 10 組 vs Control 組DEG 的Degree 排名前5 的分別是IL6、IGF1、PDGFRB、TGFB2 和BGLAP。如圖6D 和表5 示VitD 100 組 vs Control組DEG 的Degree 排名前5 的分別是TLR4、COL4A4、S100A4、COL8A1 和 COL11A1。

表2 VitD 0.1 vs Control top 前5 的基因Tab.2 VitD 0.1 vs Control Top 5 genes

表3 VitD 1vs Control top 前5 的基因Tab.3 VitD 1 vs Control Top 5 genes

表4 VitD 10 vs Control top 前5 的基因Tab.4 VitD 10 vs Control Top 5 genes

表5 VitD 100 vs Control top 前5 的基因Tab.5 VitD 100 vs Control Top 5 genes

圖6 PPI 網(wǎng)絡(luò)互作圖Fig.6 PPI network interworking diagram

3 討論

活性形式的維生素 D 能改善胎盤細(xì)胞的功能障礙[9]。隨著對維生素 D 在妊娠過程中作用機(jī)制研究的不斷深入，已有研究表明維生素 D 的缺乏與妊娠期的高血壓、糖尿病等不良妊娠結(jié)局顯著相關(guān)[10]。促炎性細(xì)胞因子的高表達(dá)會導(dǎo)致不良妊娠結(jié)局，維生素 D 可通過抑制核因子-?B 通路信號級聯(lián)，下調(diào)促炎細(xì)胞因子的表達(dá)，抑制免疫反應(yīng)進(jìn)而改善妊娠結(jié)局[11]。有研究表明，維生素 D可以通過抑制胎盤滋養(yǎng)層細(xì)胞凋亡從而減少不良妊娠結(jié)局的發(fā)生[12]。

本研究中，筆者獲得了9 個共同的DEG。其中CYP24A1 在維生素D 降解中具有重要作用，其表達(dá)增加可能是維生素 D 缺乏和疾病進(jìn)展的潛在原因[13]。在多種腫瘤組織中GAS5-AS1 的表達(dá)水平明顯低于正常組織，高表達(dá)的 GAS5-AS1 能夠抑制多種腫瘤性疾病如宮頸癌、膠質(zhì)瘤、肝細(xì)胞癌和肺癌的生長和轉(zhuǎn)移[14-17]。PCSK4 基因主要在睪丸生精細(xì)胞中轉(zhuǎn)錄，PCSK4 缺陷小鼠的精子表現(xiàn)出加速獲能、早熟頂體反應(yīng)、受精能力的受損，它在體內(nèi)表達(dá)量的下降會導(dǎo)致嚴(yán)重的雄性不育[18]。這些基因在維生素 D 和不良妊娠結(jié)局的研究中尚未見有報道，它們是否是維生素 D 缺乏引起不良妊娠結(jié)局的潛在分子標(biāo)志物，有待進(jìn)一步的深入研究。對維生素D 在妊娠過程中發(fā)揮的作用機(jī)制研究具有重要意義。筆者發(fā)現(xiàn)不同濃度的維生素D 可能通過作用于不同的信號通路影響著胎盤滋養(yǎng)細(xì)胞的遷移和侵襲過程。筆者構(gòu)建了PPI 網(wǎng)絡(luò)，并篩選出5 個關(guān)鍵基因編碼蛋白。對該網(wǎng)絡(luò)的進(jìn)一步研究將有利于理解DEG 之間的相互作用，為尋找維生素 D 缺乏引起的不良妊娠結(jié)局新的治療靶點(diǎn)提供思路。

已有研究發(fā)現(xiàn)，維生素D 缺乏可能通過影響機(jī)體炎性因子和脂肪因子水平異常，從而增加肥胖風(fēng)險，導(dǎo)致糖尿病發(fā)生[19].維生素D 可能通過全身炎癥反應(yīng)的間接途徑從而影響葡萄糖穩(wěn)態(tài)。母體維生素D 缺乏可能會增加機(jī)體氧化應(yīng)激增加促炎反應(yīng)，從而導(dǎo)致內(nèi)皮功能障礙[20]，這可能是發(fā)生子癇前期的原因之一。有研究發(fā)現(xiàn)，維生素D缺乏的患者在補(bǔ)充維生素D 后氧化應(yīng)激減少[21]。維生素D 的缺乏與NF-KB 信號相關(guān)的血管內(nèi)皮的相關(guān)功能抑制增加、細(xì)胞炎癥、維生素D 受體和1-a 輕化酶表達(dá)減少有關(guān)[22]。維生素D 通過減少炎癥細(xì)胞因子的產(chǎn)生和作用，從而減少全身炎癥并提高細(xì)胞的存活率[23]。而本次研究維生素 D刺激后差異基因主要富集在炎癥信號通路中，這就解釋了維生素D 缺乏可以通過炎癥通路來致病，也可以通過維生素D 來治療炎癥通路影響的相關(guān)疾病。

臨床通過檢測血清 25（OH）D 的量來判斷機(jī)體是否缺乏維生素 D，＜ 30 ng/mL 為不足，＜ 20 ng/mL即為缺乏，＜ 10 ng/mL 為嚴(yán)重缺乏[24]。炎癥性腸?。↖BD）是一種慢性炎癥性疾病包括潰瘍性結(jié)腸炎（UC）和克羅恩病（CD）。已有研究表明，維生素D 缺乏普遍存在于 IBD 患者中，UC 患者為31.6%，而 CD 患者為 38.1%[25]。筆者課題組前期研究表明4 個濃度中10 nmol/L 的濃度最接近人體維生素D 缺乏的生理濃度，而本次研究VitD10 vs Control 組DEG 主要富集在IBD 相關(guān)信號通路，這就表明了維生素D 的缺乏，參與了炎癥性腸病的發(fā)生和發(fā)展過程，有待進(jìn)一步的實驗研究。

有研究表明，與健康妊娠婦女相比，子癇前期（PE）維生素D 缺乏會導(dǎo)致胎盤miR-21 表達(dá)水平升高的現(xiàn)象。筆者課題組在前期研究中發(fā)現(xiàn)，用1，25（OH）2D3刺激HTR（人胎盤滋養(yǎng)細(xì)胞株HTR-8/SVneo）細(xì)胞，1，25（OH）2D3抑制HTR 細(xì)胞miR-21 的表達(dá)，減弱了miR-21 抑制HTR 細(xì)胞侵襲及遷移能力。這可以表明1，25（OH）2D3可能通過減弱PE 患者miR-21 的表達(dá)從而影響PE 的病理生理過程，從而改善不良妊娠結(jié)局[26]。

綜上所述，本研究通過對維生素 D 刺激后的胎盤滋養(yǎng)細(xì)胞進(jìn)行mRNA 測序分析。篩選出維生素 D 的缺乏在妊娠過程中可能影響的相關(guān)基因，并通過生物信息學(xué)分析了一些重要基因的功能和參與的信號通路。結(jié)果表明維生素 D 的缺乏可能通過相關(guān)信號通路里的基因影響胎盤發(fā)育和母體健康，進(jìn)而導(dǎo)致妊娠并發(fā)癥的產(chǎn)生。因此，筆者的發(fā)現(xiàn)將有助于了解維生素 D 在妊娠并發(fā)癥的發(fā)病機(jī)制中可能的作用，并為維生素 D 缺乏引起的不良妊娠結(jié)局的早期診斷、治療和預(yù)后提供重要依據(jù)。