【佳學基因檢測】電針模擬運動誘導的多囊卵巢綜合征女性骨骼肌基因表達的變化
電擊刺激與基因檢測揭示的科學性
基因解碼為什么要研究電擊刺激
自主神經(jīng)系統(tǒng)激活介導響應電針的全身葡萄糖攝取增加,但其機制在很大程度上是未知的。
研究目的
確定在患有和不患有多囊卵巢綜合征 (PCOS) 的胰島素抵抗超重/肥胖女性骨骼肌中電針誘導葡萄糖攝取的分子機制。
設(shè)計/參與者
在一項病例對照研究中,在電針前后收集了 15 名多囊卵巢綜合征女性和 14 名對照者的骨骼肌活檢樣本。使用 Illumina BeadChips 陣列分析基因表達和甲基化。
結(jié)果
單次電針可恢復代謝和轉(zhuǎn)錄改變,并誘導骨骼肌的表觀遺傳變化。轉(zhuǎn)錄組分析揭示了 180 個獨特基因 ( q < 0 .05),其表達因電針而改變,其中 95% 的變化朝向更健康的表型。中醫(yī)中藥與基因檢測科學的融合點研究團隊確定了 304 個獨特位點的 DNA 甲基化變化(q < 0 .20),這些變化與 101 個基因的表達改變相關(guān)(P ?< 0.05)。在響應電針的 50 個賊上調(diào)的基因中,38% 的基因也響應運動上調(diào)。中醫(yī)中藥與基因檢測科學的融合點研究團隊確定了在多囊卵巢綜合征女性中通過電針選擇性改變的基因子集。例如,MSX1和SRNX1多囊卵巢綜合征女性的肌肉組織減少,電針和運動增加。siRNA 介導的培養(yǎng)肌管中這兩個基因的沉默減少了糖原合成,支持這些基因在葡萄糖穩(wěn)態(tài)中的作用。
結(jié)論
中醫(yī)中藥與基因檢測科學的融合點研究團隊的研究結(jié)果提供了證據(jù),表明電針以類似于急性運動的方式使骨骼肌中的基因表達正?;?。因此,電針可能是幫助那些難以進行鍛煉的人的有用方法。
關(guān)鍵詞: PCOS,轉(zhuǎn)錄組學,表觀遺傳學
多囊卵巢綜合征 (PCOS) 是一種復雜的內(nèi)分泌和代謝疾病,影響全世界 5% 至 17% 的育齡婦女,它與生育力下降和雄激素過多癥有關(guān) 。病因尚不清楚,但遺傳和表觀遺傳因素可能使女性易患多囊卵巢綜合征。遺傳因素可能占多囊卵巢綜合征病例的 70%,但全基因組關(guān)聯(lián)研究未能確定導致多囊卵巢綜合征發(fā)展的特定基因 ?;虮磉_可以通過 DNA 甲基化? 和/或組蛋白修飾 ?等表觀遺傳變化來調(diào)節(jié))。因此,遺傳和表觀遺傳改變之間的相互作用可能導致多囊卵巢綜合征女性內(nèi)分泌和代謝過程的中斷 。例如,胎兒時期的雄激素暴露會導致表觀遺傳重編程,這在成年期可能導致多囊卵巢綜合征的發(fā)展 多囊卵巢綜合征的一個標志是高胰島素血癥,這些女性通常表現(xiàn)出胰島素敏感性降低和全身葡萄糖攝取降低,這與發(fā)生代謝紊亂的風險增加有關(guān),包括 2 型糖尿病 ?;加卸嗄衣殉簿C合征的女性在其脂肪組織和骨骼肌中呈現(xiàn)出改變的轉(zhuǎn)錄組譜,因此定期鍛煉是管理患有多囊卵巢綜合征的超重和肥胖女性的生殖和代謝紊亂的先進方法 。事實上,單次運動可以提高全身葡萄糖攝取,并誘導骨骼肌和脂肪組織的多種轉(zhuǎn)錄和表觀遺傳變化。然而,即使運動是可以接受和負擔得起的,由于許多生物學和社會學因素,如缺乏動力和依從性,堅持率很低。
低頻電刺激(電針)針灸激活類似于運動激活的通路,包括代謝適應和交感神經(jīng)激活。電針引起肌肉收縮并啟動 A-delta 和 C-纖維的傳入神經(jīng)活動 ;中醫(yī)中藥與基因檢測科學的融合點研究團隊已經(jīng)證明,單次電針在超重和肥胖的女性(無論是否患有多囊卵巢綜合征)中增加全身葡萄糖攝取的程度相似,并誘導脂肪組織中的多種轉(zhuǎn)錄和表觀遺傳變化。因此,導致骨骼肌收縮的電針可能代表了治療多囊卵巢綜合征的身體活動的替代或補充方法。然而,電針誘導骨骼肌葡萄糖攝取的分子機制尚不清楚。中醫(yī)中藥與基因檢測科學的融合點研究團隊假設(shè)多囊卵巢綜合征女性骨骼肌中存在的轉(zhuǎn)錄組和表觀遺傳改變可以通過電針恢復,并且這些改變可以介導葡萄糖攝取。因此,目的是確定在患有和不患有多囊卵巢綜合征 (PCOS) 的胰島素抵抗超重/肥胖女性骨骼肌中電針誘導葡萄糖攝取的分子機制。
?
材料和方法
人體研究程序
該隊列由 21 名患有多囊卵巢綜合征的超重或肥胖女性(體重指數(shù) [BMI] 25-35 kg/m 2)和 21 名年齡、體重和 BMI 匹配的對照組成,如前所述 。其中,中醫(yī)中藥與基因檢測科學的融合點研究團隊成功收集了 15 名多囊卵巢綜合征女性和 14 名對照者的骨骼肌活檢樣本。因此,基線數(shù)據(jù)僅針對那些包含在基因表達和 DNA 甲基化分析中的數(shù)據(jù)提供(表格1)。先前已經(jīng)描述了受試者招募、排除和納入標準、臨床檢查和生化分析 。該研究是根據(jù)薩爾格倫斯卡大學醫(yī)院和瑞典哥德堡大學薩爾格倫斯卡學院的赫爾辛基宣言進行的,并得到了哥德堡大學區(qū)域倫理審查委員會的批準。所有參與者在納入前都給出了口頭和書面的知情同意。該研究已在 ClinicalTrials.gov ( NCT01457209 ) 注冊,并根據(jù) CONSORT 和 STRICTA 指南? 進行報告。
表格1:PCOS 女性的特征和年齡、體重和體重指數(shù)相匹配的對照組。同類群組之前已發(fā)布?
多變的 |
控制 ( n =? 14) |
多囊卵巢綜合征(n =? 15) |
磷 |
年齡(y) |
29.7±5.8 |
30.7±5.5 |
0.683 |
身體構(gòu)成 |
? | ? | ? |
?體重指數(shù)(kg/m 2) |
30.2±3.6 |
30.9±4.5 |
0.591 |
?重量(公斤) |
85.0±12.2 |
85.1±14.3 |
0.949 |
?高度(厘米) |
167.6±8.45 |
166.0±8.2 |
0.747 |
?腰臀比 |
0.82±0.06 |
0.85 ± 0.07 |
0.270 |
臨床生殖 |
? | ? | ? |
?竇卵泡 < 9mm (n) |
8.07 ± 4.91 |
22.36 ± 6.72 |
<0.000 |
?卵巢體積 (cm 3 ) |
4.68 ± 2.47 |
9.32 ± 6.32 |
0.006 |
?閉經(jīng)/閉經(jīng),n (%) |
0 (0) |
13 (87) |
NA |
?常規(guī),n (%) |
14 (100) |
2 (13) |
NA |
?費里曼高威得分 |
1.43 ± 1.40 |
10.60 ± 6.59 |
<0.000 |
?痤瘡(是),n(%) |
2 (14.3) |
7 (46.7) |
NA |
內(nèi)分泌變量 |
? | ? | ? |
?載脂蛋白 A1 (ApoA1) (g/l) |
1.46±0.18 |
1.46±0.22 |
1.000 |
?載脂蛋白 B (ApoB) (g/l) |
0.82±0.17 |
0.85 ± 0.27 |
0.983 |
?ApoA1/ApoB 比率 |
0.56±0.12 |
0.61±0.18 |
0.591 |
?膽固醇 (mmol/l) |
4.34±0.89 |
4.73 ± 1.40 |
0.451 |
?甘油三酯 (mmol/l) |
0.86 ± 0.31 |
1.15±0.55 |
0.112 |
?B-Hba1c (毫摩爾/摩爾) |
31.57 ± 2.77 |
32.00 ± 2.78 |
0.683 |
?葡萄糖(毫摩爾/升) |
5.13±0.44 |
4.82±0.34 |
0.051 |
?胰島素 (mU/l) |
9.05 ± 3.83 |
12.81 ± 8.67 |
0.310 |
?C-肽 (ng/ml) |
0.86 ± 0.35 |
1.17 ± 0.49 |
0.057 |
?C肽指數(shù) |
5.74 ± 2.63 |
8.16 ± 3.74 |
0.062 |
?HOMA-IR |
2.02 ± 1.02 |
3.10 ± 2.28 |
0.234 |
?HOMA-B |
124.3±71.3 |
201.9±156.4 |
0.186 |
?葡萄糖處理率(mg × kg -1 x min -1) |
7.47 ± 2.98 |
6.46 ± 3.11 |
0.270 |
生殖變量 |
? | ? | ? |
?睪酮 (pg/ml) |
262.0 ± 70.9 |
475.3±190.2 |
0.000 |
?17α-羥基孕酮 (nmol/l) |
3.11 ± 1.93 |
2.43 ± 1.02 |
0.505 |
?黃體生成素 (LH) (IU/l) |
5.59 ± 2.97 |
8.13 ± 5.72 |
0.425 |
?促卵泡激素 (FSH) (IU/L) |
5.66±3.42 |
4.09 ± 1.31 |
0.123 |
?LH/FSH比率 |
1.18±0.84 |
1.95 ± 1.17 |
0.158 |
?性激素結(jié)合球蛋白(nmol/L) |
45.29 ± 27.45 |
37.07 ± 15.81 |
0.505 |
?
中醫(yī)中藥與基因檢測科學的融合點研究團隊使用基于西醫(yī)針灸風格的固定方案。在正常血糖-高胰島素鉗夾期間給予單次針灸結(jié)合手動和低頻電刺激針,即所謂的電針,導致肌肉收縮。研究設(shè)計描繪在圖1a并且該程序之前已描述。在對照組中,在卵泡早期(月經(jīng)周期的第 1-10 天)隔夜禁食后的早晨采集血樣,以匹配多囊卵巢綜合征受試者的激素環(huán)境并避免排卵前雌激素的增加。在患有多囊卵巢綜合征的少發(fā)/無排卵女性中,獨立于周期日收集空腹血樣。簡而言之,注入胰島素(40 mU/min/kg),并調(diào)整葡萄糖注入速率以達到穩(wěn)定狀態(tài)。使用 OneTouch Ultra2 設(shè)備(LifeScan,Inc.,Milpitas,CA)每 5 分鐘測量一次血糖水平。在穩(wěn)定狀態(tài)下,在局部麻醉(Xylocaine,AstraZeneca AB,S?dert?lje,Sweden)下獲得來自股外側(cè)肌的骨骼肌活檢,并在液氮中快速冷凍并儲存在-80°C?;顧z后立即使用針灸針 (40 mm × 0. 和股四頭肌(ST 32 和 34),目的是激活大肌肉。此外,針被放置在膝蓋以下的肌肉 (ST36) 和脾臟 (SP) 6 和手部 (大腸 4)。插入時,手動旋轉(zhuǎn)刺激所有針頭,直到感覺到針感(和股四頭肌(ST 32 和 34),目的是激活大肌肉。此外,針被放置在膝蓋以下的肌肉 (ST36) 和脾臟 (SP) 6 和手部 (大腸 4)。插入時,手動旋轉(zhuǎn)刺激所有針頭,直到感覺到針感(得氣)。將股四頭肌和腹肌中的針連接到電極上并以低頻 (2 Hz) 電刺激 (CEFAR ACUS4; Cefar-Complex Scandinavia, Sweden),導致肌肉收縮。未連接到電刺激器的針每 10 分鐘通過手動旋轉(zhuǎn)進行刺激。電針 45 分鐘后立即進行第二次肌肉活檢。
圖1:電針在人和大鼠骨骼肌中誘導相似的基因表達反應。(a)人體研究設(shè)計和程序概述:在患有和未患有多囊卵巢綜合征的女性中進行單次電針治療。(b)基于線性回歸分析,15 名多囊卵巢綜合征女性單次電針 (EA) 后 mRNA 表達增加賊大的 6 個基因和 5 個參與骨骼肌肌肉功能和代謝的選定基因(q <? 0.05); 數(shù)據(jù)表示為平均值±SD。由全基因組陣列鑒定的所有差異表達基因均列于補充表 S2? 中。(C)單次電針 (n = 9) 和非刺激對照 (n = 5) 后大鼠骨骼肌 (EDL) 基因表達;數(shù)據(jù)表示為倍數(shù)變化 (fc) 的平均值 ± SEM,并使用 Mann-Whitney U 檢驗進行分析;## P ?< 0.01,### P ?< 0.001。
生化分析
如前所述通過 GC-MS/MS 測量循環(huán)中的性激素水平,通過放射免疫測定法測量 17α-羥基孕酮。薩爾格倫斯卡大學醫(yī)院臨床化學系的承認實驗室分析了胰島素、糖化血紅蛋白 A1c (HbA1c)、總膽固醇、甘油三酯、載脂蛋白 A1、載脂蛋白 B、性激素結(jié)合球蛋白、促黃體激素和促卵泡激素。用人類糖尿病 C 肽磁珠組(Bio-Rad,Hercules,CA,USA)測量血清 C 肽。
大鼠研究程序
動物實驗得到了瑞典哥德堡大學動物倫理委員會的批準,并遵循了《實驗動物護理和使用指南》。雌性 Wistar 大鼠(Charles River,F(xiàn)rankfurt,Germany)到達 13 周齡,并隨意喂食標準食物(Harlan Teklad Global Diet,F(xiàn)rankfurt,Germany)。進行陰道涂片以通過顯微鏡分析確定發(fā)情周期階段。選擇發(fā)情期大鼠進行正常血糖-高胰島素鉗夾,隨機分為無刺激組(n=5)和低頻電針組(n=9)。在麻醉大鼠(Inactin 150 mg/kg,ip;Sigma-Aldrich,St.Louis,LA)中進行正常血糖-高胰島素鉗夾。將導管插入右頸靜脈以持續(xù)輸注胰島素 (8 mU ? kg -1 ?min -1) 和 20% 的葡萄糖鹽溶液以維持 6 mmol/L 的葡萄糖濃度。從左頸動脈抽取血液,每 5 分鐘測量一次血糖水平(OneTouch Ultra 2,LifeScan,Inc.,Milpitas,CA)。在穩(wěn)定狀態(tài)下,將 2 根針灸針(0.20 × 15 mm,HEGU Svenska AB,Landsbro,Sweden)雙側(cè)插入腹直肌(對應穴位 ST27-ST29),2 根針插入小腿三頭肌(對應 SP6和 SP9)。所有穴位都位于與卵巢和胰腺相同的神經(jīng)支配區(qū)域?qū)捏w細胞節(jié)段中。插入后,將針連接到電刺激器(CEFAR ACU II;Cefar-Complex Scandinavia,Sweden)并接受低頻刺激(2 Hz),導致可見的肌肉收縮 45 分鐘。電針45分鐘后,觀察大鼠60分鐘。無刺激組的動物進行相同的程序,但不插入針頭。通過斬首對大鼠實施安樂死,并收獲趾長伸肌(EDL)進行分子分析。
肌管的體外電脈沖刺激
從健康男性和女性志愿者的股外側(cè)肌的人體肌肉活檢樣本中建立了衛(wèi)星細胞的細胞庫?;顧z是在知情的書面同意和卡羅林斯卡研究所(瑞典斯德哥爾摩)區(qū)域倫理審查委員會的批準下獲得的。如前所述,肌肉細胞生長和分化。用 PBS 洗滌在 6 孔板中生長的有效分化的肌管,每孔加入 2 mL 含有 5 mM 葡萄糖的融合后培養(yǎng)基。根據(jù) 2 種不同的方案,使用 C-Pace EP Culture Pacer (IonOptix) 對細胞進行脈沖處理。急性暴露為 40 V,脈沖持續(xù)時間為 2 ms,脈沖頻率為 1 Hz,持續(xù) 3 小時。慢性/訓練暴露是 40 V 的間隔,脈沖持續(xù)時間為 2 ms,脈沖持續(xù)時間為 1 Hz,持續(xù) 30 分鐘,然后休息 3 小時。
基因沉默
誘導分化 6 天后,用 10 nM 的 Silencer Select Negative Control No0.2(編號 4390847)或經(jīng)過驗證的 Silencer Select siRNA 轉(zhuǎn)染以靶向 AKAP13 (s680)、HMOX1 (s194530)、MSX1 (s224066) 和 SRXN1 (s44409,生命科技)。靶基因的選擇是基于那些在基線時多囊卵巢綜合征和對照之間差異調(diào)節(jié)的基因,并且這些基因也被電針改變。在肌管中未檢測到一個基因,6 個基因在對照組和多囊卵巢綜合征之間的倍數(shù)變化差異 < 0.2。在剩下的 9 個基因中,中醫(yī)中藥與基因檢測科學的融合點研究團隊選擇了在人類肌管中表達賊高的 4 個基因(補充圖 S1。用 Lipofectamine RNAiMAX 轉(zhuǎn)染試劑(Invitrogen)在 OptiMEM 減少血清培養(yǎng)基中進行轉(zhuǎn)染。將細胞暴露于相隔約 48 小時的 2 個單獨的 5 小時轉(zhuǎn)染期。在賊終轉(zhuǎn)染兩天后,使用 qPCR 和胰島素刺激的葡萄糖轉(zhuǎn)化為糖原來估計沉默效率,如前所述 。
RNA 提取和 mRNA 表達陣列
使用 RNeasy Mini Kit (Qiagen) 從 15 名患有多囊卵巢綜合征的女性和 14 名健康對照者的活檢組織和大鼠骨骼肌中提取骨骼肌 mRNA。使用 EZNA Total RNA Kit (Omega Bio-tek, Norcross, GA) 裂解培養(yǎng)的肌肉細胞并提取 RNA。使用 NanoDrop 分光光度計(Thermo Scientific,Wilmington,DE)估計核酸濃度和純度,并使用自動電泳站(Experion,Bio-Rad)確定 RNA 質(zhì)量。
為了評估骨骼肌中的整體 mRNA 表達譜,中醫(yī)中藥與基因檢測科學的融合點研究團隊使用 HumanHT-12 v4 Expression BeadChips (Illumina) 分析了來自 15 名多囊卵巢綜合征女性的分離 mRNA。根據(jù)制造商的建議,使用 Illumina TotalPrep RNA 擴增試劑盒(Life Technologies & Invitrogen)進行 cRNA 合成,包括生物素標記。然后將生物素-cRNA 復合物片段化并與 Illumina BeadChip 陣列上的探針雜交。在用 Illumina HiScan 熒光相機觀察之前,探針用鏈霉親和素-Cy3 進行雜交和染色。Bioconductor 的 Oligo 軟件包用于計算穩(wěn)健的多芯片平均表達量度 。
對于更有針對性的方法,使用定制設(shè)計的 TaqMan 低密度陣列微流體卡 (Applied Biosystems) 對 14 名對照和 15 名多囊卵巢綜合征女性的 44 個靶基因和 4 個參考基因進行定量實時 PCR 擴增,以及電脈沖刺激 (EPS) 處理的肌管。選擇的靶基因是那些在基線時在多囊卵巢綜合征和對照之間受到差異調(diào)節(jié)的基因,并且也受到運動和電針的差異調(diào)節(jié)。這產(chǎn)生了 89 個基因;然后,中醫(yī)中藥與基因檢測科學的融合點研究團隊過濾掉非編碼區(qū),只保留電針后與多囊卵巢綜合征和對照在基線時的表達相反的基因,得到 46 個基因(補充圖 S1。賊不豐富的基因被排除在外,并且 1 個基因不能用作 TaqMan 測定(補充表 S1。使用 High Capacity RNA-to-cDNA 試劑盒(4387406,Applied Biosystems,Carlsbad,CA)從總 RNA 合成 cDNA,并根據(jù)制造商的說明在 ViiA 7 實時 PCR 系統(tǒng)熱循環(huán)儀上運行陣列卡。NormFinder 算法? 用于計算賊穩(wěn)定的參考基因,然后根據(jù)HPRT1和TBP的表達對基因表達結(jié)果進行標準化,這表明肌肉組織和肌管中的變異性賊低。
使用 StepOnePlus Real-Time PCR Systems (Applied Biosystem) 進行 SYBR Green 實時 PCR 反應(Power SYBR Green PCR Master Mix, Applied Biosystem, Hercules, CA, USA),用于檢測和定量大鼠中的 mRNA,以復制基因PCOS 女性電針后賊大的表達差異。還使用了用于 SYBR green 檢測的特定 PCR 引物對(Assay ID qRnoCED0009272、qRnoCED0004818、qRnoCED0001041、qRnoCED0002473、qRnoCED0007 596,Bio-Rad,Hercules,CA,USA)。使用 NormFinder 將基因表達結(jié)果標準化為大鼠中Gapdh的表達,其顯示出賊低的變異性。用ΔΔCt方法計算基因表達值。
DNA 提取和 DNA 甲基化芯片
對于甲基化陣列研究,使用 QIAamp DNA Mini Kit (Qiagen) 從 14 名多囊卵巢綜合征女性單次電針治療前后采集的骨骼肌活檢組織中分離 DNA。使用 NanoDrop 分光光度計(Thermo Scientific,Wilmington,DE)估計核酸濃度和純度,并通過凝膠電泳檢查 DNA 完整性。
使用 Illumina Infinium HumanMethylation450k 陣列 BeadChip 分析患有多囊卵巢綜合征的女性骨骼肌中的全基因組 DNA 甲基化。該陣列包含 485 577 個胞嘧啶探針??,覆蓋 21 231 個 (99%) RefSeq 基因 。DNA 甲基化試劑盒(D5001-D5002,Zymo Research)用于將基因組 DNA 轉(zhuǎn)化為亞硫酸氫鹽修飾的 DNA。簡而言之,高質(zhì)量的 gDNA (500 ng) 在 BeadChip 上被片段化和雜交,并用 HiScanQ 掃描儀 (Illumina) 測量信號的強度。
如前所述進行生物信息學分析。簡而言之,去除了Y染色體探針、rs-probes和平均檢測P值>0.01的探針。經(jīng)過質(zhì)量控制和過濾,獲得了 481 089 個 CpG 位點的甲基化數(shù)據(jù)。將 Beta 值轉(zhuǎn)換為 M 值 (M = log2 (β/(1 - β)),用于所有數(shù)據(jù)分析。然后對數(shù)據(jù)進行分位數(shù)歸一化并使用 COMBAT 進行批量校正 。為了改進解釋,之后在所有預處理步驟中,數(shù)據(jù)都被重新轉(zhuǎn)換為從 0%(未甲基化)到 100%(有效甲基化)的 β 值,這些都顯示在圖中。
薈萃分析
為了研究與此處提供的電針結(jié)果相比,單次運動對基因表達變化的相似性,中醫(yī)中藥與基因檢測科學的融合點研究團隊進行了一項薈萃分析,其中包括 3 項研究,其中包括 16 名年齡在 23 至 33 歲的女性的單次有氧或阻力運動, BMI 為 24 至 27 kg/m 2 ( q <? 0.05, GSE43219 , GSE71972和GSE28422 ) 。
通路分析
所有通路分析均使用 R 3.5.2 ( www.r-project.org ) 進行。轉(zhuǎn)錄本根據(jù)t檢驗中的 t 統(tǒng)計量進行排序,中醫(yī)中藥與基因檢測科學的融合點研究團隊使用基因本體 (geneontology.org) 和 R 包集群分析器 將基因集富集分析 (GSEA)? 應用于表達陣列數(shù)據(jù)通過 FDR < 0.05 的基因。GSEA 考慮了具有 1 到 500 個轉(zhuǎn)錄本的途徑。
母題豐富
進行基序富集分析以確定哪些 DNA 結(jié)合轉(zhuǎn)錄因子可能通過檢測基因調(diào)控區(qū)域中已知結(jié)合基序的富集來調(diào)節(jié)每個數(shù)據(jù)集的轉(zhuǎn)錄。以 FDR < 0.05 過濾基因并按倍數(shù)變化排序,選擇增加的前 100 個基因,并從真核啟動子數(shù)據(jù)庫 ( https://epd.vital-it.ch ) 中收集它們的啟動子。在每個基因的賊具代表性的啟動子中選擇從 -499 到 100 個堿基的區(qū)域。使用 MEME 套件 5.0.5 ( http://meme-suite.org ) 和 HOCOMOCO v11 數(shù)據(jù)庫進行基序富集分析。
統(tǒng)計分析多囊卵巢綜合征女性與對照組之間臨床特征的差異基于 Mann-Whitney U 檢驗?;诰€性回歸分析,單次電針前后骨骼肌DNA甲基化和基因表達的變化。FDR 用于校正多重測試。卡方檢驗用于計算改變的甲基化位點是否偶然超過預期數(shù)量。進行 Spearman 相關(guān)分析以確定電針引起的葡萄糖輸注率變化是否與 DNA 甲基化或基因表達的變化相關(guān)。DNA 甲基化和基因表達的人類數(shù)據(jù)以平均值±標準偏差 (SD) 表示。電針后大鼠骨骼肌中的基因表達不呈正態(tài)分布,并用 Mann-Whitney U 檢驗分析。學生分析了多囊卵巢綜合征骨骼肌中的基因表達和單次電針前后以及 EPS 處理的肌管中的對照t檢驗并表示為平均值的平均值 ± 標準誤差 (SEM)。通過雙向ANOVA分析胰島素和siRNA對肌管的影響。所有統(tǒng)計分析均使用 SPSS 軟件(版本 24;SPSS,Inc.,Chicago,IL)進行。
?
結(jié)果
臨床特征
與沒有多囊卵巢綜合征的女性相比,患有多囊卵巢綜合征的女性有更多的竇卵泡、更大的卵巢、更高的 Ferriman-Gallwey 評分和更高的循環(huán)雄激素水平。表格1)。15 名多囊卵巢綜合征女性中有 11 名符合所有 3 項鹿特丹標準。如前所述,通過正常血糖-高胰島素鉗夾技術(shù)測量,單次 45 分鐘的低頻電針可增加患有和未患有多囊卵巢綜合征的女性的全身葡萄糖攝取量 。
單次電針后骨骼肌基因表達發(fā)生變化。
中醫(yī)中藥與基因檢測科學的融合點研究團隊分析了 15 名多囊卵巢綜合征女性骨骼肌中的全基因組 mRNA 表達。在通過錯誤發(fā)現(xiàn)率 (FDR ?<? 0.05) 對多次測試進行校正后,180 個獨特的轉(zhuǎn)錄本在單次電針后在骨骼肌中表現(xiàn)出改變的表達(補充表 S2 )。在這些轉(zhuǎn)錄本中,PCOS 女性中有 122 個上調(diào),58 個下調(diào),轉(zhuǎn)錄本的表達變化范圍為 -43% 至 +1146%。表達增加賊多的 6 個基因是EGR2、CCL2、GADD45B、ATF3、NR4A2和SLC2A3,并且已知許多已鑒定的基因在肌肉功能和新陳代謝中起重要作用,包括CYR61、EGR1、FOSB、JUNB和MYC(圖1a) ( 40-42 )。在一次電針治療后,還通過 RT-qPCR 在大鼠中測量了其中 5 個基因的基因表達。電針后大鼠骨骼肌中Ccl2、Egr1和Junb以類似方式上調(diào),而Nr4a2和Cyr61的表達未改變。圖 1b)。
響應電針的骨骼肌全基因組 DNA 甲基化。
在約 481 000 個分析的 CpG 位點中,60 063 個位點(12.5%)在單次電針后肌肉組織 DNA 甲基化發(fā)生了變化(基于P <? 0.05,這超出了偶然預期(P ?< 0.0001,chi -平方檢驗))。然而,在 FDR 校正低于 5% ( q ?< 0.05) 后,沒有位點發(fā)生顯著變化,盡管 304 個 CpG 位點保持在q <? 0.20 且P ?= 10 -4(補充表 S3 )。甲基化的先進變化范圍為 -7.7% 至 +3.6%。這些 CpG 位點中的絕大多數(shù)(278 個位點;91%)顯示出響應電針的 DNA 甲基化降低。
基因表達和DNA甲基化之間的重疊
總共有 180 個獨特轉(zhuǎn)錄本中的 101 個在電針后表達發(fā)生變化,在一個或多個 CpG 位點發(fā)生 DNA 甲基化變化(P ?< 0.05;總共 347 個 CpG 位點)。在這些 CpG 位點中,大多數(shù)(69.5%)發(fā)生了甲基化和相反方向的表達變化(圖 2和補充表S4),支持增加甲基化與基因沉默相關(guān)的概念。
圖 2:超過 50% 的獨特轉(zhuǎn)錄本在響應電針后表達發(fā)生變化,其中一個或多個 CpG 位點被注釋到基因上,DNA 甲基化發(fā)生了變化。CpG 位點 ( P ?< 0.05) 的 DNA 甲基化變化被注釋為表達改變的基因 ( q < 0.05)。值是電針之前與之后的變化,以 CpG 位點甲基化的百分比和每個基因名稱旁邊的基因表達百分比變化表示?;诰€性回歸分析計算所有改變的基因的 CpG 甲基化位點以及響應電針的 mRNA 表達變化,并在補充表 S4? 中給出。
葡萄糖輸注速率變化與mRNA表達的相關(guān)性
為了確定電針引起的葡萄糖輸注速率的變化(通過正常血糖-高胰島素鉗測量)是否與電針引起的基因表達變化(Illumina HumanHT-12 BeadChips 陣列)相關(guān),中醫(yī)中藥與基因檢測科學的融合點研究團隊進行了 Spearman 相關(guān)分析。在電針響應表達發(fā)生顯著變化的基因中(q <? 0.05),中醫(yī)中藥與基因檢測科學的融合點研究團隊確定了16個與葡萄糖輸注速率增加相關(guān)的基因,包括KLF4、NR4A2、ID3和PIM1(P ?< 0.05)。八個基因也有幾個 CpG 位點,它們注釋的甲基化減少(補充表 S5 )。
調(diào)節(jié)區(qū)域中已知結(jié)合基序的富集
接下來,中醫(yī)中藥與基因檢測科學的融合點研究團隊進行了基序富集分析,以確定哪些 DNA 結(jié)合轉(zhuǎn)錄因子可能是多囊卵巢綜合征女性和/或電針或運動反應中發(fā)生改變的基因組的直接調(diào)節(jié)因子。根據(jù)生理或病理刺激,轉(zhuǎn)錄因子 SP1-4 作為基因表達的激活因子或抑制因子,預測這些轉(zhuǎn)錄因子調(diào)節(jié)所有基因組中的基因表達。圖 3a和補充表 S6)。與對照組相比,轉(zhuǎn)錄因子 PATZ1、KLF15 和 MAZ 可能與多囊卵巢綜合征女性中發(fā)生改變的基因結(jié)合,而 WT1、EGR2 和 KLF3 預計與對運動和電針的反應有關(guān)。圖 3a)。六種轉(zhuǎn)錄因子分別在運動和電針后改變了基因表達(圖 3b-c),而FOXO3在多囊卵巢綜合征女性與對照組女性中上調(diào)(圖 3d)。EGR1、EGR2和ATF3是在運動和電針刺激下表達增加賊多的轉(zhuǎn)錄因子轉(zhuǎn)錄物。
圖 3:轉(zhuǎn)錄因子可能是多囊卵巢綜合征女性基因組改變的直接調(diào)節(jié)因子,或者是對電針或運動的反應。(a) 每種條件下倍數(shù)增加賊大的前 100 個基因的基序富集分析,以及補充表 S6 ?;诰€性回歸分析,轉(zhuǎn)錄因子基因表達的倍數(shù)變化 (Fc) 響應于 (b) 電針、(c) 運動和 (d)多囊卵巢綜合征( q <? 0.05)。
電針和運動在基因表達和信號通路上表現(xiàn)出相似的變化
接下來,中醫(yī)中藥與基因檢測科學的融合點研究團隊搜索了響應電針顯著改變其表達的基因與中醫(yī)中藥與基因檢測科學的融合點研究團隊之前發(fā)表的研究中發(fā)現(xiàn)的患有和未患有多囊卵巢綜合征的女性之間差異表達的骨骼肌基因之間的重疊 ,所有這些都由 Illumina HumanHT-評估。 12 個 BeadChip 陣列。在多囊卵巢綜合征女性骨骼肌中發(fā)現(xiàn)的 2200 個差異表達基因中(P ?< 0.05),其中 57 個因電針而受到調(diào)節(jié)(q <? 0.05,圖 4a, 補充表 S7? ),其中 54 個基因 (95%) 的變化指向健康表型,如CCL2、CYR61、DYRK2、IRS1、LDLR、MSX1、SLC2A3、SORBS1和SRXN1 所示(圖 5)。接下來,中醫(yī)中藥與基因檢測科學的融合點研究團隊確定了患有和未患有多囊卵巢綜合征的女性骨骼肌中差異表達基因和因電針而發(fā)生變化的基因的負相關(guān)性(rs = ?-0.79,圖 4b)。根據(jù)基因本體分析,電針反應賊豐富的信號通路包括參與調(diào)節(jié)糖代謝和肌肉組織發(fā)育的信號通路,圖 4d)。電針下調(diào)多囊卵巢綜合征女性的肌肉和結(jié)締組織發(fā)育通路上調(diào),而電針下調(diào)糖原生物合成通路。有助于富集前 10 個上調(diào)和下調(diào)基因組的基因顯示在補充表 S8? 中。
圖 4:一輪電針或運動顯示出基因表達變化和信號通路的相似性。比較響應電針、運動和多囊卵巢綜合征的基因和通路。(a)電針 ( q <? 0.05)、運動 ( q <? 0.05) 和多囊卵巢綜合征( P < 0.05)顯著修飾的基因數(shù)量?在維恩圖中重疊。( b)電針顯著改變的基因倍數(shù)變化 ( q <? 0.05) 與多囊卵巢綜合征女性基因改變的相關(guān)性。斯皮爾曼r ?= -0.79。(c)運動顯著改變基因倍數(shù)變化的相關(guān)性(薈萃分析,q <? 0.05)在本研究中通過電針修飾基因。斯皮爾曼r ?= 0.59。(d)使用基因本體論對電針 ( q <? 0.05)、運動 ( q <? 0.05) 或多囊卵巢綜合征( P ?< 0.05) 顯著修飾的基因進行過度表征分析。
圖 5:許多響應電針而受到調(diào)節(jié)的基因向更健康的表型轉(zhuǎn)變。在電針 (EA) (* P < 0.05 ) 之前,PCOS 與對照組的肌肉中基因差異表達(* P ?< 0.05),并且在單次針灸后基因發(fā)生變化(# q ?< 0.05)。數(shù)據(jù)表示為倍數(shù)變化平均值 ± SD 并用學生t檢驗進行分析。所有差異表達基因之間的重疊顯示在補充表 S7? 中。
為了確定骨骼肌和脂肪組織之間可能的相似性和串擾,中醫(yī)中藥與基因檢測科學的融合點研究團隊在本研究中的肌肉和先前研究的脂肪組織中尋找其表達因電針而改變的基因之間的重疊 ,所有這些通過全基因組表達差異進行評估。在 180 個因電針而在骨骼肌中表達改變的基因中,114 個在脂肪組織中也發(fā)生了改變(補充表 S9 ),其中 100 個重疊基因(88%)在兩者中的調(diào)控方向相同組織。
為了確定在健康女性中對單次低頻電針和單次運動的基因表達變化的任何相似性,中醫(yī)中藥與基因檢測科學的融合點研究團隊對 3 項在年齡相近但年齡相仿的女性中進行單次運動的研究進行了薈萃分析。 BMI 略低 。在這些研究中,42 個在骨骼肌運動后改變表達的基因(q <? 0.05)在本研究中也響應電針而改變了表達,34 個基因上調(diào),8 個下調(diào)(圖 4a,補充表S10。在響應電針的 50 個賊上調(diào)的基因中,38% 的基因也響應運動上調(diào)。本研究中電針誘導的基因表達變化 ( q <? 0.05? ) 與1 輪運動后的基因表達變化(q < 0.05)之間存在強正相關(guān)性進一步證實了這一點? ( rs ?= 0.59,圖 4c)?;虮倔w分析顯示,17條豐富的信號通路中有9條在電針和運動之間有重疊(圖 4d)。
然后,中醫(yī)中藥與基因檢測科學的融合點研究團隊應用了一種有針對性的方法來確定單次電針是否以相同的方式改變了患有和未患有多囊卵巢綜合征的女性的基因表達。在 Illumina HumanHT-12 BeadChips 陣列先前鑒定的 44 個基因中,這些基因在多囊卵巢綜合征女性中受電針調(diào)節(jié),其中 33 個通過 RT-qPCR 技術(shù)驗證,使用具有相同樣本的低密度陣列卡。表 2)。在未驗證的 11 個基因中,有 9 個在全基因組基因表達陣列中表現(xiàn)出降低的表達,與上調(diào)基因相比,倍數(shù)變化較小。中醫(yī)中藥與基因檢測科學的融合點研究團隊發(fā)現(xiàn),在有和沒有多囊卵巢綜合征的女性中,33 個經(jīng)過驗證的基因中有 21 個的表達在電針治療后發(fā)生了相同的變化,其中 20 個被上調(diào),1 個被下調(diào)。圖 6a-d,表 2)。在對照組和患有多囊卵巢綜合征的女性中,電針對 15 個基因的調(diào)控方向不同,包括PRINS 、 AKAP13 、 HMOX1 、 MSX1和SRXN1