薛源++李巖++徐思

摘 要：miRNA調(diào)控機(jī)體對(duì)訓(xùn)練應(yīng)激的適應(yīng)能力，其調(diào)控模式影響有氧耐力可訓(xùn)練性。循環(huán)miRNA（c-miRNA）與miRNA密切相關(guān)，分析不同有氧耐力可訓(xùn)練性水平群體的運(yùn)動(dòng)適應(yīng)性c-miRNA表達(dá)變化差異，可獲得與訓(xùn)練敏感性相關(guān)的c-miRNA表達(dá)譜特征。結(jié)果顯示：高可訓(xùn)練性表型在有氧訓(xùn)練應(yīng)激誘導(dǎo)下差異表達(dá)17條c-miRNA，其中11條上調(diào)，6條下調(diào)或平穩(wěn)，其調(diào)控功能主要涉及低氧適應(yīng)通路和脂肪酸β氧化代謝的關(guān)鍵基因表達(dá)。提示：有氧耐力可訓(xùn)練性差異與機(jī)體低氧適應(yīng)能力和脂肪酸β氧化供能能力，對(duì)有氧訓(xùn)練應(yīng)激誘導(dǎo)作用的反應(yīng)性密切相關(guān)。miRNA調(diào)控并整合機(jī)體應(yīng)激適應(yīng)性基因表達(dá)。c-miRNA差異表達(dá)譜可用來評(píng)估有氧耐力可訓(xùn)練性，預(yù)測(cè)有氧能力發(fā)展?jié)摿Α?/p>

關(guān) 鍵 詞：運(yùn)動(dòng)生物化學(xué)；有氧耐力；可訓(xùn)練性；循環(huán)微核苷酸；最大攝氧量；關(guān)聯(lián)模型；

運(yùn)動(dòng)性適應(yīng)

中圖分類號(hào)：G804.7 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：1006-7116（2015）04-0139-06

Correlation between aerobic endurance trainability difference and the expression

spectrum characteristics of sport adaptive circulating microRNA

XUE Yuan1，LI Yan2，XU Si3

（1.School of Physical Education，Sichuan Normal University，Chengdu 610101，China；

2.School of Physical Education，Ludong University，Yantai 264025，China；

3.Sichuan Academy of Medical Science，Sichuan Provincial Hospital，Chengdu 610031，China）

Abstract： miRNA regulates the bodys ability to adapt to training stress； its regulating pattern affects aerobic endurance trainability. Circulating miRNA （c-miRNA） is closely related to miRNA； by analyzing differences in the changing of sport adaptive c-miRNA expression of groups of people at different aerobic endurance trainability levels， the authors can acquire the expression spectrum characteristics of c-miRNA sensitively correlative with training. Results： induced by aerobic training stress， the phenotype with high trainability differentially expressed 17 c-miRNAs， 11 of which were up-regulated， 16 of which were down-regulated or steady， its regulating function mainly involved with hypoxia adaptation pathway and the key gene expression of fatty acid β-oxidation metabolism. Hint： aerobic endurance trainability difference is close related to the reactivity of the bodys hypoxia adaptation ability and fatty acid β-oxidation energy supply ability to the induction function of aerobic training. miRNA regulates and integrates the bodys stress adaptive gene expression. C-miRNA differential expression spectrum can be used to evaluate aerobic endurance trainability， and to predict aerobic ability development potential.

Key words： sports biochemistry；aerobic endurance；trainability；c-miRNAs；maximum oxygen uptake；correlation model；sports adaption

杰出運(yùn)動(dòng)能力的形成與運(yùn)動(dòng)天賦和長(zhǎng)期刻苦訓(xùn)練密不可分。運(yùn)動(dòng)性適應(yīng)能力好的運(yùn)動(dòng)員的可訓(xùn)練性更強(qiáng)，成績(jī)提高更快，訓(xùn)練效果更好[1]。表觀遺傳學(xué)機(jī)制調(diào)控營(yíng)養(yǎng)、訓(xùn)練等環(huán)境應(yīng)激因素誘導(dǎo)的適應(yīng)過程的基因選擇性表達(dá)，其基因調(diào)控模式具有個(gè)體差別并可穩(wěn)定遺傳[2]，可能是產(chǎn)生運(yùn)動(dòng)能力及可訓(xùn)練性個(gè)體差異的原因。目前對(duì)于表觀遺傳學(xué)基因調(diào)控模式與運(yùn)動(dòng)能力的可訓(xùn)練性的具體關(guān)聯(lián)所知甚少，因而無法將其作為指標(biāo)應(yīng)用于運(yùn)動(dòng)員選材和個(gè)性化訓(xùn)練。

microRNA基因表達(dá)調(diào)控作用是環(huán)境應(yīng)激適應(yīng)性基因表達(dá)的主要表觀遺傳學(xué)調(diào)控方式[3]，其調(diào)控模式與環(huán)境應(yīng)激適應(yīng)能力的個(gè)體差異有關(guān)，應(yīng)激適應(yīng)過程中的microRNA表達(dá)譜反映了其調(diào)控模式和強(qiáng)度[4]。研究表明，microRNA分子（以下簡(jiǎn)稱miRNA）參與了運(yùn)動(dòng)性適應(yīng)基本生理過程的調(diào)控，如肌組織肥大和心肌/骨骼肌收縮力增強(qiáng)[5]、血管增生[6]、線粒體合成與酶活性提高[7]等。組織內(nèi)miRNA分泌進(jìn)入血液循環(huán)成為循環(huán)miRNAs（circulating miRNAs，以下簡(jiǎn)稱c-miRNA），二者在表達(dá)譜上具有穩(wěn)定關(guān)聯(lián)[8]。c-miRNA可能是運(yùn)動(dòng)應(yīng)激誘導(dǎo)機(jī)體應(yīng)答過程中整合多系統(tǒng)適應(yīng)性反應(yīng)的關(guān)鍵因子[9]，其表達(dá)譜特征與運(yùn)動(dòng)訓(xùn)練方式和訓(xùn)練效果存在關(guān)聯(lián)，在運(yùn)動(dòng)性適應(yīng)的不同狀態(tài)下，c-miRNAs表達(dá)譜發(fā)生特征性改變[10]，因此推測(cè)：運(yùn)動(dòng)適應(yīng)性c-miRNAs表達(dá)譜可以作為反映運(yùn)動(dòng)能力的可訓(xùn)練性的“分子指紋”，在選材階段準(zhǔn)確預(yù)估運(yùn)動(dòng)能力的發(fā)展?jié)摿Γ蛴脕碓u(píng)估訓(xùn)練效果/競(jìng)技狀態(tài)，作為制訂個(gè)性化訓(xùn)練方案的依據(jù)。

有氧耐力水平是決定運(yùn)動(dòng)能力的主要因素之一，不同個(gè)體的有氧耐力訓(xùn)練效果存在顯著差異[11]。本研究以最大攝氧量（VO2max）為主要指標(biāo)評(píng)估有氧耐力訓(xùn)練效果和訓(xùn)練敏感性，同步檢測(cè)運(yùn)動(dòng)適應(yīng)性c-miRNA表達(dá)譜，提取和分析有氧耐力訓(xùn)練的高敏感表型個(gè)體的運(yùn)動(dòng)性適應(yīng)c-miRNA表達(dá)譜特征，建立運(yùn)動(dòng)性適應(yīng)c-miRNA表達(dá)譜與有氧運(yùn)動(dòng)能力訓(xùn)練敏感性的關(guān)聯(lián)模型，作為早期預(yù)測(cè)和評(píng)估有氧運(yùn)動(dòng)能力發(fā)展?jié)摿瓦m應(yīng)特點(diǎn)的指標(biāo)應(yīng)用于運(yùn)動(dòng)員選材和個(gè)性化訓(xùn)練方案制訂。

1 研究對(duì)象與方法

1.1 研究對(duì)象

四川省某文武學(xué)校初中部學(xué)員92人，男性，平均年齡（13.6±0.5）歲；身高（171.85±6.28）cm；體重（59.60±4.55） kg。既往健康，無傷病史，參加試驗(yàn)前未受過系統(tǒng)的有氧運(yùn)動(dòng)能力訓(xùn)練。

1.2 有氧耐力訓(xùn)練方案

參加實(shí)驗(yàn)學(xué)員進(jìn)行為期15周的系統(tǒng)性有氧運(yùn)動(dòng)能力訓(xùn)練。每周三、五下午16：00-17：30進(jìn)行5 000 m長(zhǎng)跑訓(xùn)練，負(fù)荷強(qiáng)度：1～3周65%個(gè)人最大心率±3次/min；4～10周70%個(gè)人最大心率±3次/min；11～15周75%個(gè)人最大心率±3次/min。以polar表（瑞典產(chǎn)）監(jiān)控靶心率維持速度。訓(xùn)練方案由教練監(jiān)督執(zhí)行。

1.3 有氧運(yùn)動(dòng)能力相關(guān)指標(biāo)檢測(cè)

所有指標(biāo)在15周有氧耐力訓(xùn)練前后各檢測(cè)1次，分別記為基礎(chǔ)值（baseline）和訓(xùn)練值（post-ex）。

1）最大攝氧量相對(duì)值（R-VO2max）測(cè)試采用12 min跑法[12]，在400 m標(biāo)準(zhǔn)跑道上進(jìn)行，受試者佩戴polar表，按指令起跑，盡最大努力跑完12 min，讀取跑動(dòng)距離s并換算成英里，按以下公式折算成R-VO2max：

R-VO2max（mL/（min·kg））=（s-0.318）/0.0278 ①

2）個(gè)體乳酸閾（ILT）測(cè)定采用遞增負(fù)荷法[13]，功率自行車（瑞典產(chǎn)Monark894E、839E）起始負(fù)荷為50 W，每2 min遞增40 W，轉(zhuǎn)速維持在50 r/min，直至力竭（力竭標(biāo)準(zhǔn)：連續(xù)2次不能維持規(guī)定轉(zhuǎn)速）。取中指指尖血20 μL，采用YSI1500便攜式血乳酸自動(dòng)分析儀測(cè)定血乳酸值。分別采取安靜狀態(tài)下、每級(jí)負(fù)荷后即刻及恢復(fù)期第2、5、8、10、15 min血樣。測(cè)試結(jié)果記錄在表中，采用Stegmann方法判定ILT[14]。

1.4 c-miRNA表達(dá)譜檢測(cè)和實(shí)時(shí)熒光定量PCR法驗(yàn)證

1）c-miRNAs表達(dá)值測(cè)定。

15周有氧耐力訓(xùn)練開始前清晨空腹安靜狀態(tài)下肘靜脈取血5 mL檢測(cè)c-miRNAs表達(dá)值，記為基線表達(dá)值。完成VO2max基礎(chǔ)值測(cè)試后分別于1 h內(nèi)、恢復(fù)24和48 h肘靜脈取血5 mL檢測(cè)c-miRNAs表達(dá)值，記為測(cè)試后0、恢復(fù)24 h、恢復(fù)48 h表達(dá)值。

2）血漿總RNA提取。

TRIZOL法提取血漿總RNA，RNase-free ddH2O溶解提取的總RNA分子。紫外分光光度儀測(cè)定RNA純度和濃度，OD 260/280>1.9說明總RNA純度較好，②甲醛變性瓊脂糖電泳檢測(cè)總RNA完整性較好，未發(fā)生降解。

3）microarray檢測(cè)c-miRNA差異表達(dá)譜。

總RNA按照Ambions miRNA Isolation Kit說明書分離miRNA。使用miRCURYTM Array Labelling kit（Ambion）進(jìn)行miRNA標(biāo)記操作。microRNA芯片雜交按照miRCURYTM Array microarray kit （Exiqon）說明書進(jìn)行。Genepix 4000B雙通道激光圖像掃描儀（Molecular Devices公司）以635 nm單波長(zhǎng)進(jìn)行掃描，Genepix Pro 6.0軟件分析數(shù)據(jù)。數(shù)據(jù)預(yù)處理以后，根據(jù)各張芯片的global mean 進(jìn)行片間校正，用SAM （significance analysis of microarrays，SAM .version2.1）篩選差異表達(dá)miRNAs。篩選條件：FDR （false discovery rate）<5%，F(xiàn)old change>2或<0.5。

4）實(shí)時(shí)熒光定量PCR（qRT-PCR）驗(yàn)證差異表達(dá)c-miRNA相對(duì)表達(dá)量。

（1）c-miRNA莖環(huán)反轉(zhuǎn)錄合成cDNA：根據(jù)mirBASE提供的miRNA種子區(qū)序列，設(shè)計(jì)莖環(huán)反轉(zhuǎn)錄引物（北京六合華大基因公司合成），20 μL反應(yīng)體系（5×First Strand synthesis Buffer，4 μL；dNTPmi 1uL；RNase Inhibitor，40units；Reverse Transcriptase（M-MLV），1μL；Stem-loop RT Primer，1μL；Total RNA 500 ng；RNase-free ddH2O），反應(yīng)條件：65 ℃加熱5 min，放入冰中2 min。然后放入PCR儀16 ℃熱激30 min，42℃ 60 min，85 ℃，5 min停止反轉(zhuǎn)錄。

（2）qRT-PCR檢測(cè)c-miRNA表達(dá)水平：miRNA定量引物（北京六合華大基因公司）。反應(yīng)體系按照SYBR? Premix Ex TaqTMⅡ（Perfect Real Time）說明書進(jìn)行配置。使用Bio-Rad iCycler PCR System進(jìn)行Real Time PCR反應(yīng)，程序：95 ℃ 5 min；95 ℃ 10 s；60℃ 30 s，40cycle。以上循環(huán)結(jié)束后進(jìn)行65～95 ℃的融解曲線分析。每個(gè)樣品平行3次，溶解曲線為單一峰。按照2-ΔΔCt計(jì)算樣品的相對(duì)表達(dá)豐度。

1.5 c-miRNA差異表達(dá)譜的獲得及其與有氧耐力指標(biāo)的相關(guān)性

1）計(jì)算訓(xùn)練前后R-VO2max變化率η（R-VO2max）并采用K-means法進(jìn)行聚類分析，根據(jù)聚類分析結(jié)果將學(xué)員分為對(duì)有氧訓(xùn)練的高敏感表型組（High Response genotype Group，簡(jiǎn)稱HR組）和普通敏感表型組（Common Response genotype Group，簡(jiǎn)稱CR組）。比較兩組c-miRNA表達(dá)水平，參照Davidsen標(biāo)準(zhǔn)[10]提取HR組c-miRNA差異表達(dá)譜。

2）計(jì)算HR組差異表達(dá)的c-miRNAs在基線值-0 h階段（以下簡(jiǎn)稱為階段Ⅰ）和0～24 h階段（以下簡(jiǎn)稱為階段Ⅱ）的表達(dá)量變化率η（c-miRNA）：

η（c-miRNA）=[（階段終點(diǎn)表達(dá)量-階段起點(diǎn)表達(dá)量）/階段起點(diǎn)表達(dá)量]×100%

3）計(jì)算η（c-miRNA）與η（R-VO2max）和η（ILT）之間的Pearson相關(guān)系數(shù)，分析其相關(guān)性。

1.6 數(shù)據(jù)處理與分析

采用SPSS 17.0軟件對(duì)結(jié)果進(jìn)行處理。各項(xiàng)指標(biāo)數(shù)據(jù)均采用（±s）表示。有氧耐力訓(xùn)練前后各項(xiàng)指標(biāo)的變化率η（Χ）=[（訓(xùn)練值-基礎(chǔ)值）/基礎(chǔ)值]×100%。均數(shù)間比較采用兩獨(dú)立樣本的t檢驗(yàn)，P<0.05表示兩組間的差異有統(tǒng)計(jì)學(xué)意義。聚類分析采用k-means法。采用Bivariate過程進(jìn)行兩組數(shù)據(jù)之間的相關(guān)性分析，計(jì)算Pearson相關(guān)系數(shù)并用r表示。

2 研究結(jié)果及分析

2.1 訓(xùn)練前后R-VO2max與ILT變化率及其相關(guān)性

3 討論

本研究探討了運(yùn)動(dòng)適應(yīng)性c-miRNA表達(dá)譜與有氧耐力的可訓(xùn)練性之間的關(guān)聯(lián)，結(jié)果發(fā)現(xiàn)：η（R-VO2max）顯著高于普通水平的表型個(gè)體有17個(gè)c-miRNA差異表達(dá)，其功能主要富集于低氧適應(yīng)通路和脂肪酸β氧化代謝途徑，調(diào)控毛細(xì)血管增生、血紅素生成、糖酵解相關(guān)酶功能、線粒體脂肪酸β氧化功能及骨骼肌重塑等適應(yīng)性生理反應(yīng)的關(guān)鍵基因表達(dá)，影響機(jī)體對(duì)訓(xùn)練應(yīng)激的適應(yīng)性。其作用特點(diǎn)為協(xié)同性調(diào)控多系統(tǒng)對(duì)耐力訓(xùn)練的應(yīng)激適應(yīng)。本研究結(jié)果提示：運(yùn)動(dòng)能力及其可訓(xùn)練性與miRNA調(diào)控基因選擇性表達(dá)的模式有關(guān)。運(yùn)動(dòng)應(yīng)激誘導(dǎo)的特征性c-miRNA表達(dá)譜與杰出運(yùn)動(dòng)能力表型存在關(guān)聯(lián)，可以用來預(yù)測(cè)運(yùn)動(dòng)能力的發(fā)展?jié)摿Γ瑢⒕哂邢嚓P(guān)運(yùn)動(dòng)天賦的運(yùn)動(dòng)員早期選拔出來加以系統(tǒng)培養(yǎng)，或用以評(píng)估運(yùn)動(dòng)員的訓(xùn)練適應(yīng)狀態(tài)/競(jìng)技狀態(tài)，制訂個(gè)性化訓(xùn)練方案。

遺傳選材的一個(gè)主要問題是如何評(píng)估個(gè)體的運(yùn)動(dòng)能力可訓(xùn)練性，預(yù)測(cè)其運(yùn)動(dòng)能力發(fā)展?jié)摿15]。杰出運(yùn)動(dòng)能力是遺傳優(yōu)勢(shì)和營(yíng)養(yǎng)、訓(xùn)練等環(huán)境因素共同作用下的積累效應(yīng)。機(jī)體對(duì)訓(xùn)練應(yīng)激的適應(yīng)能力有條件地影響遺傳優(yōu)勢(shì)的表型化程度，決定了運(yùn)動(dòng)能力的可訓(xùn)練性。應(yīng)激誘導(dǎo)miRNA表達(dá)，調(diào)控基因表達(dá)強(qiáng)度，產(chǎn)生適應(yīng)性的個(gè)體差異[17]。比較分析不同可訓(xùn)練性群體的適應(yīng)性miRNA表達(dá)譜差異，不僅可以建立高可訓(xùn)練性群體的c-miRNA表達(dá)特征關(guān)聯(lián)模型，還提供了訓(xùn)練應(yīng)激與機(jī)體相互作用影響運(yùn)動(dòng)能力的有價(jià)值的信息。

最大攝氧量是評(píng)價(jià)有氧耐力的主要指標(biāo)。個(gè)體VO2max主要受遺傳影響，耐力訓(xùn)練可以在一定程度上改變VO2max水平，相同訓(xùn)練條件下VO2max變化率反映了心肺有氧機(jī)能可訓(xùn)練性的個(gè)體差異[18]。本研究采用VO2max變化率作為區(qū)別有氧能力可訓(xùn)練性的指標(biāo)分析運(yùn)動(dòng)應(yīng)激誘導(dǎo)的c-miRNA表達(dá)差異，采用Targetscan、miRanda和mirBASE聯(lián)合分析差異表達(dá)的miRNA的靶基因，并進(jìn)行GO分析和KEGG pathway分析，結(jié)合靶基因主要生理機(jī)能和代謝通路的共性特征研究HR組c-miRNA表達(dá)譜的調(diào)控特征。結(jié)果顯示：差異表達(dá)的17條miRNA的靶基因主要富集于低氧適應(yīng)通路和PPAR信號(hào)通路，HR組miRNA表達(dá)模式促進(jìn)相關(guān)代謝通路的應(yīng)激適應(yīng)能力，并顯示出多靶整合效應(yīng)特征。

低氧應(yīng)激激活HIF信號(hào)通路適應(yīng)是獲得有氧耐力訓(xùn)練效果的重要途徑[19]。常氧分壓條件下，有氧耐力訓(xùn)練引起的低氧應(yīng)激主要發(fā)生在心肌/骨骼肌等組織內(nèi)[20]。有氧工作引起耗氧量增加，組織內(nèi)氧分壓顯著降低，抑制HIF-1a降解，上調(diào)HIF-1a活性，調(diào)控下游EPO、VEGF等低氧適應(yīng)基因表達(dá)水平上調(diào)，產(chǎn)生心肌和骨骼肌肌纖維毛細(xì)血管密度增加、Mb及血紅蛋白水平提高等適應(yīng)效應(yīng)，改善心肌和骨骼肌供血供氧，提高心臟和骨骼肌有氧工作能力。

運(yùn)動(dòng)應(yīng)激誘導(dǎo)PPARδ和PPARα基因表達(dá)上調(diào)促進(jìn)心肌骨骼肌脂肪酸β氧化對(duì)有氧運(yùn)動(dòng)能力具有重要作用。線粒體脂肪酸β氧化供能是心肌有氧代謝的主要供能方式，占心肌耗能的70%以上，對(duì)于維持心泵功能具有重要作用[21]。PPARδ調(diào)節(jié)PPARα促進(jìn)心肌細(xì)胞脂肪酸β氧化，提高心泵功能[22]。進(jìn)行有氧耐力訓(xùn)練或在骨骼肌細(xì)胞過表達(dá)PPARδ都可以使PPARδ蛋白含量增加并獲得相似的表型改變和骨骼肌重塑：慢肌纖維收縮蛋白基因和脂肪酸β氧化相關(guān)基因表達(dá)上調(diào)，氧化型肌纖維增多，氧代謝效率提高，有氧耐力水平提高。阻斷PPARδ的作用，則不會(huì)產(chǎn)生上述改變[23]。

HR組差異表達(dá)的miRNA的功能主要是調(diào)控HIF-1a信號(hào)通路和PPAR信號(hào)通路對(duì)有氧訓(xùn)練應(yīng)激的敏感性和響應(yīng)強(qiáng)度，其綜合適應(yīng)效應(yīng)與其運(yùn)動(dòng)適應(yīng)表型特征是一致的。上調(diào)HIF-1a信號(hào)通路對(duì)低氧應(yīng)激的敏感性可以使循環(huán)系統(tǒng)攜氧能力和心肌/骨骼肌的供氧水平提高更顯著。PPARδ和PPARα對(duì)運(yùn)動(dòng)應(yīng)激誘導(dǎo)作用的反應(yīng)性更好，可以更顯著上調(diào)心肌/骨骼肌脂肪酸β氧化代謝能力，心肌/骨骼肌有氧工作能力提高更明顯。心泵功能增強(qiáng)提高心輸出量，骨骼肌有氧代謝效率提高，可以增加提高動(dòng)靜脈氧分壓差。VO2max水平取決于心輸出量和動(dòng)靜脈氧分壓差[12]，因此上調(diào)HIF-1a信號(hào)通路和PPAR信號(hào)通路都可以使VO2max的可訓(xùn)練性更好。同時(shí)，骨骼肌氧化型肌纖維增多，脂代謝供能比例增加還可以提高ILT水平[24]，這與HR組具有較高的η（ILT）是一致的。

HR組miRNA差異表達(dá)譜特征顯示出對(duì)多個(gè)應(yīng)激適應(yīng)通路的整合作用。差異表達(dá)的17條miRNA中，除miR210、miR423-5p和miR126專性調(diào)控低氧應(yīng)激適應(yīng)外，其余14條miRNA均同時(shí)作用于低氧應(yīng)激適應(yīng)通路和PPAR信號(hào)通路，在調(diào)控特點(diǎn)上顯示出整合效應(yīng)。本研究結(jié)果顯示VO2max與ILT具有高度相關(guān)性，一方面證明了骨骼肌有氧代謝能力提高對(duì)VO2max的可訓(xùn)練性具有重要的影響[18]，同時(shí)也提示我們miRNA對(duì)多種應(yīng)激適應(yīng)途徑的整合作用可能是有氧耐力的心肺機(jī)能和骨骼肌代謝的協(xié)調(diào)一致的重要調(diào)節(jié)機(jī)制。

分階段分析η（c-miRNA）與η（R-VO2max）和η（ILT）相關(guān)性結(jié)果表明：運(yùn)動(dòng)應(yīng)激誘導(dǎo)的脂肪酸β氧化供能能力的提高對(duì)有氧耐力的可訓(xùn)練性十分重要。在階段Ⅰ，η（miR183）與η（R-VO2max）和η（ILT）均呈顯著負(fù)相關(guān)，miR183[25]負(fù)性調(diào)節(jié)HIF-3a[20]和CROT[23]基因表達(dá)，促進(jìn)HIF-1a的功能，同時(shí)抑制線粒體脂肪酸β氧化，研究表明：低氧應(yīng)激提高糖酵解供能水平，抑制脂肪酸β氧化[19]。該結(jié)果提示：訓(xùn)練過程中低氧應(yīng)激對(duì)脂肪酸β氧化代謝功能的抑制作用不利于有氧耐力可訓(xùn)練性的提高。在階段Ⅱ，η（miR141）和η（miR15a）均與η（R-VO2max）和η（ILT）正相關(guān)。研究表明：低氧應(yīng)激誘導(dǎo)的HIF-2a[19]的高表達(dá)在應(yīng)激結(jié)束后仍持續(xù)抑制PPARa和CPT1的表達(dá)和功能，損害線粒體脂肪酸β氧化能力[23]，miR141[26]和miR15a[27]均可上調(diào)PPARa的表達(dá)，PPARa又上調(diào)CPT1的表達(dá)，該結(jié)果提示：抑制低氧應(yīng)激對(duì)脂肪酸β氧化能力的損害有助于提高VO2max和ILT的可訓(xùn)練性。

此外，本研究結(jié)果還提示：機(jī)體抗損傷能力對(duì)于有氧耐力可訓(xùn)練性也有重要的影響。在階段Ⅰ，η（R-VO2max）和η（c-miR133）呈負(fù)相關(guān)關(guān)系。miR133是肌肉特異性microRNA，運(yùn)動(dòng)過程中肌肉損傷漏出是其血漿水平升高主要原因[28]，該結(jié)果提示：骨骼肌抗運(yùn)動(dòng)損傷能力越強(qiáng)，訓(xùn)練效果越好。

綜上所述，本研究發(fā)現(xiàn)：在典型有氧訓(xùn)練周期的運(yùn)動(dòng)性適應(yīng)過程中，17條c-miRNAs的表達(dá)特征與有氧耐力的可訓(xùn)練性之間存在密切的聯(lián)系并初步建立了關(guān)聯(lián)模型，同時(shí)也提出了一種采用表觀遺傳學(xué)指標(biāo)進(jìn)行運(yùn)動(dòng)員選材和指導(dǎo)訓(xùn)練的方法，以c-miRNAs表達(dá)特征作為生物標(biāo)記評(píng)估運(yùn)動(dòng)能力的發(fā)展?jié)摿陀?xùn)練適應(yīng)狀態(tài)。進(jìn)一步的研究將對(duì)c-miRNAs特征表達(dá)譜在運(yùn)動(dòng)員選材、訓(xùn)練效果/競(jìng)技狀態(tài)評(píng)估等方面的應(yīng)用價(jià)值做出更為全面和深入的探索。

注釋：

① 該公式是根據(jù)同一受試者的R-VO2max直接測(cè)量值與12 min跑距離數(shù)值之間配對(duì)歸納得出的數(shù)值推算公式，反映的是R-VO2max與12 min跑距離之間的數(shù)值換算關(guān)系，如果采用公制單位該公式可以寫成

R-VO2max（mL/（kg·min））=（s/1 609.34-0.318）/0.027 8

由于原始文獻(xiàn)和相關(guān)文獻(xiàn)中引用該公式時(shí)12 min跑距離數(shù)值均采用英制，所以本研究也采用相同的表述方式。該公式是數(shù)值推算公式，僅表示R-VO2max與12 min跑距離之間的數(shù)值換算關(guān)系，不能從公式中得出單位。

② “OD260/280值”是提取核酸后檢驗(yàn)其純度和完整性的指標(biāo)，利用紫外分光光度計(jì)在260 nm和280 nm處的光密度比值判斷核酸提取物的純度和完整性，一般以O(shè)D260/280值在1.9～2.0范圍內(nèi)表示提取物的純度和完整性較好，過高說明雜蛋白含量較高，過低說明RNA降解明顯。

參考文獻(xiàn)：

[1] 馬力宏. 人類體質(zhì)與運(yùn)動(dòng)能力的遺傳學(xué)基礎(chǔ)[J]. 天津體育學(xué)院學(xué)報(bào)，2001，16（4）：1-5.

[2] 于紅. 表觀遺傳學(xué)：生物細(xì)胞非編碼RNA調(diào)控的研究進(jìn)展[J]. 遺傳，2009，31（11）：1077-1086.

[3] 鄭永霞，焦炳華. miRNA的生物形成及調(diào)控基因表達(dá)機(jī)制[J]. 生命的化學(xué)，2010，30（6）：821-826.

[4] Leung A K L，Sharp P A. MicroRNA functions in stress responses[J]. Molecular Cell，2010，40（2）：205-215.

[5] Williams A H，Liu N，Van Rooij E，et al. MicroRNA control of muscle development and disease[J]. Curr Opin Cell Biol，2009，21（4）：461–469.

[6] Zhang C. MicroRNAs in vascular biology and vascular disease[J]. J Cardiovasc Transl Res，2010，3（3）：235–240.

[7] Dang C V. Re-thinking the Warburg effect with Mycmicroma naging glutamine metabolism[J]. Cancer Res，2010，70（5）：859–862.

[8] Kosaka N，Iguchi H，Ochiya T. Circulating microRNA in body fluid： a new potential biomarker for cancer diagnosis and prognosis[J]. Cancer Sci，2010，101（10）：2087–2092.

[9] Baggish A L，Hale A，Weiner R B，et al. Dynamic regulation of circulating microRNA during acute exhaustive exercise and sustained aerobic exercise training[J]. J Physiol，2011，589（16）：3983-3994.

[10] Davidsen P K，Gallagher I J，Hartman J W，et al. High responders to resistance exercise training demonstrate differential regulation of skeletal muscle microRNA expression[J]. J Appl Physiol，2011，110（2）：309–317.

[11] 李義志，駱軍保，湯長(zhǎng)發(fā). HSP1A1基因+190G>C多態(tài)與有氧運(yùn)動(dòng)能力表型指標(biāo)的關(guān)聯(lián)性研究[J]. 天津體育學(xué)院學(xué)報(bào)，2013，28（2）：180-184.

[12] 李俊勇，任晉軍，曹峰銳. 最大攝氧量、無氧閾和最大攝氧量平臺(tái)同高校男生12 min跑和1 000 m跑的相關(guān)性分析[J]. 北京體育大學(xué)學(xué)報(bào)，2010，33（08）：65-67.

[13] 林文弢，鄭選梅，徐國(guó)琴，等. 不同氧濃度環(huán)境下個(gè)體乳酸閾的比較[J]. 廣州體育學(xué)院學(xué)報(bào)，2012，32（6）：81-84.

[14] Stegman H，Kindermann W，Schnabel A. Lactate kinetics and individual anaerobic threshold[J]. International Journal of Sports Medicine，1981，2（3）：160-165.

[15] 楊若愚，王予彬，沈勛章，等. 基因多態(tài)性與杰出運(yùn)動(dòng)能力[J]. 中國(guó)組織工程研究，2014，18（7）：1121-1128.

[16] 李江華，王智慧，朱小娟. 代謝組學(xué)運(yùn)動(dòng)選材與傳統(tǒng)生化選材方法的比較[J]. 體育學(xué)刊，2013，20（1）：123-128.

[17] Joshua T. Mendell Eric N. Olson. MicroRNAs in stress signaling and human disease[J]. Cell，2012，48（6）：1172–1187.

[18] 王國(guó)軍，溫含，武雅瓊. 跑節(jié)省化、最大攝氧量評(píng)價(jià)普通人有氧耐力的效果比較[J]. 中國(guó)組織工程研究，2013，17（7）：1265-1272.

[19] Mounier R，Pialoux V，Schmitt L，et al. Effects of acute hypoxia tests on blood markers in high-level endurance athletes[J]. Eur J Appl Physiol.，2009，106（5）：713-720.

[20] 趙琛，譚斐. 低氧誘導(dǎo)因子-1激活的調(diào)節(jié)[J]. 生命科學(xué)，2013，25（1）：40-46.

[21] Cheng L，Ding G，Qin Q，et al. Cardiomyocyte-restricted peroxisome proliferator-activated receptor-delta deletion perturbs myocardial fatty acid oxidation and leads to cardiomyopathy[J]. Nat Med，2004，10（11）：1245-1250.

[22] Barger P M，Brandt J M，Leone T C，et al. Deactivation of peroxisome proliferator-activated receptor αduring cardiac hypertorphic growth[J]. J Clin Invest，2000，105（12）：1723-1730.

[23] 胡琴，李隆貴. 心肌脂肪酸氧化酶的基因調(diào)控機(jī)制及PPARα的作用[J]. 中國(guó)病理生理雜志，2002，18（12）：1552-1555.

[24] 劉昭強(qiáng)，朱寄天. 心肌乳酸代謝及乳酸對(duì)心臟功能的影響[J]. 成都體育學(xué)院學(xué)報(bào)，1999，25（3）：11-16.

[25] Vohwinkel C U，Lecuona E，Sun H，et al. Elevated CO（2） levels cause mitochondrial dysfunction and impair cell proliferation[J]. J Biol Chem，2011，286（43）：37067-37076.

[26] Mangala L S，Zhang Y，He Z，et al. Effects of simulated microgravity on expression profile of microRNA in human lymphoblastoid cells[J]. J Biol Chem，2011，286（37）：32483-32490.

[27] Sun L L，Jiang B G，LI W T，et al. MicroRNA-15a positively regulates insulin synthesis by inhibiting uncoupling protein-2 expression[J]. Diabetes Res Clin Pract，2011，91（1）：94-100.

[28] Chen J F，Mandel E M，Thomson J M，et al. The role of microRNA-1 and microRNA-133 in skeletal muscle proliferation and differentiation[J]. Nat Genet，2006，38（2）：228-233.