吐爾遜阿依·牙力坤，曹 行，王 瓊,2，劉玲玲，劉武軍*

（1.新疆農(nóng)業(yè)大學(xué)動物科學(xué)學(xué)院，新疆烏魯木齊 830052；2.吐魯番職業(yè)技術(shù)學(xué)院，新疆吐魯番 838000）

脂肪不僅是一種重要的儲能物質(zhì)，而且在調(diào)節(jié)動物體能量平衡過程中起著重要作用，并與動物產(chǎn)肉量、肉品質(zhì)、肉風(fēng)味等產(chǎn)肉性狀密切相關(guān)。尾脂肪也被人類用作食物。但是，由于日常飲食中過多的脂肪被認(rèn)為對人體健康有害，而脂肪功能的紊亂可導(dǎo)致許多心血管疾病的發(fā)生，如高血脂、II 型糖尿病等，因此動物脂肪沉積的分子調(diào)控機(jī)制備受關(guān)注。

長鏈非編碼RNA（long non-coding RNA，lncRNA）是一類長度大于200 個核苷酸，具有調(diào)控細(xì)胞增殖、分化、凋亡及自噬等多種生物學(xué)過程的非編碼RNA 分子[1]。動物脂肪沉積是一個復(fù)雜的生物過程，受功能基因、非編碼RNA 和成脂相關(guān)信號通路等系列遺傳因子的級聯(lián)調(diào)控，其中l(wèi)ncRNA 起重要的調(diào)控作用[2]。近年來，越來越多的研究發(fā)現(xiàn)mRNA 和lncRNA 在動物脂肪沉積中發(fā)揮重要作用。前人發(fā)現(xiàn)了幾種調(diào)控脂肪生成的lncRNA，例如3T3-L1 細(xì)胞中的lncRNANEAT1[3]和小鼠中的Blnc1[4]。在人類中，lncRNA ADINR 通過激活CCAAT 增強(qiáng)子結(jié)合蛋白α（C/EBPα）的轉(zhuǎn)錄來促進(jìn)脂肪形成[5]。lncRNAH19 在骨髓間充質(zhì)干細(xì)胞進(jìn)入脂肪細(xì)胞的過程中抑制脂肪細(xì)胞的分化[6]。據(jù)報(bào)道，lncRNA MEG3 參與了人類脂肪干細(xì)胞在成脂和成骨分化之間的平衡[7]。最近，已經(jīng)發(fā)現(xiàn)了幾種調(diào)節(jié)綿羊脂肪形成的lncRNA。

Mohammad 等[8]對伊朗大尾羊和小尾羊進(jìn)行了轉(zhuǎn)錄組比較分析，研究了綿羊肥尾lncRNA 的表達(dá)及其在脂肪沉積中的可能作用，并鑒定了7 個大尾和小尾之間差異表達(dá)基因，通過lncRNA-mRNA 相互作用的基因調(diào)控網(wǎng)絡(luò)，發(fā)現(xiàn)了3 個與脂質(zhì)代謝相關(guān)的重要模塊。Miao等[9]在綿羊中通過深度測序已經(jīng)確定了心臟、皮膚、肌肉、乳腺和脂肪組織等的轉(zhuǎn)錄組譜，分析了小尾寒羊與多塞特羊之間皮下脂肪組織的轉(zhuǎn)錄組信息。Ma 等[10]為了評估蘭州大尾羊、小尾寒羊以及藏羊這些綿羊品種中與尾脂肪沉積和發(fā)育相關(guān)的lncRNA 和mRNA，對這3 個個體進(jìn)行了高通量RNA 測序。通過3 組的RNA 測序數(shù)據(jù)顯示長尾巴的綿羊和細(xì)尾巴的綿羊組之間有10個差異表達(dá)的基因和37 個差異表達(dá)的lncRNA。差異表達(dá)基因和差異表達(dá)的lncRNA 的靶基因的基因本體論和途徑分析發(fā)現(xiàn)，脂肪酸代謝和與脂肪酸伸長有關(guān)的途徑均富集，這些途徑有助于脂肪沉積。

迄今為止，關(guān)于lncRNA 在綿羊脂肪沉積中的調(diào)控功能及其如何調(diào)節(jié)mRNA 表達(dá)的報(bào)道較少。本研究選擇新疆地方綿羊巴什拜羊和巴什拜羊與盤羊雜交二代作為研究對象，采集肥尾和小尾綿羊群體的尾部組織，利用全轉(zhuǎn)錄組技術(shù)，篩選出在尾脂中高表達(dá)的差異MSTRG.37980、MSTRG.38164，對其靶基因進(jìn)行預(yù)測，并對該lncRNA 及其靶基因的表達(dá)規(guī)律進(jìn)行分析，初步確定二者表達(dá)量及其在調(diào)控綿羊脂肪代謝中的作用，以期為鑒定與綿羊脂肪代謝調(diào)節(jié)機(jī)制相關(guān)的功能基因提供依據(jù)。

1 材料和方法

1.1 實(shí)驗(yàn)樣品的采集 本次實(shí)驗(yàn)所需的動物來自于新疆塔城裕民縣地區(qū)隨機(jī)挑選相同飼養(yǎng)條件下的10 月齡的6 只巴什拜羊及6 只野生盤羊×巴什拜羊的雜交二代羊，公母各3 只，屠宰后立即采集尾脂、皮下脂、背最長肌、股四頭肌、臂三頭肌、肝臟、瘤胃、十二指腸、盲腸、大腸、小腸共11 種組織，并將采集好的新鮮組織裝入2 mL RNase-Free 凍存管，暫存于液氮中，帶回實(shí)驗(yàn)室后將全部樣品迅速轉(zhuǎn)移至-80℃冰箱保存?zhèn)溆谩?/p>

1.2 實(shí)驗(yàn)方法

1.2.1 差異表達(dá)lncRNA 的選擇 前期對6 只巴什拜羊和6 只野生盤羊× 巴什拜羊的雜交二代羊的尾脂組織中的RNA-seq 結(jié)果共檢測到了728 個差異表達(dá)的lncRNA；進(jìn)行差異lncRNA 的靶基因預(yù)測后進(jìn)行靶基因的功能富集分析，發(fā)現(xiàn)MSTRG.37980 的靶基因SCD（硬脂酰輔酶A 去飽和酶）參與PPAR 信號通路，MSTRG.38164 的靶基因GPAM（線粒體甘油-3-磷酸酰基轉(zhuǎn)移酶）與脂肪部位脂肪沉積有關(guān)。

1.3 RNA 的提取 將肌肉組織樣從液氮中拿出放入（取樣品約100 mg）提前預(yù)冷好的研缽中進(jìn)行磨碎，使用Trizol 法進(jìn)行總RNA 的提取，取2 μL 總RNA 利用1%的瓊脂糖凝膠進(jìn)行完整性檢驗(yàn)。

1.4 RNA 濃度的檢測 在Nanodrop 分光光度儀檢測孔內(nèi)加入1 μL 的RNA 原液進(jìn)行RNA 濃度檢測，分別記錄不同樣品的濃度和OD260/280的數(shù)值。其中D260nm/D260nm=1.8～2.0 的樣本可用于后續(xù)實(shí)驗(yàn)，于-80℃保存待用。

1.5 引物設(shè)計(jì)及合成并PCR 擴(kuò)增檢測

1.5.1 引物的設(shè)計(jì) 從2 種群體尾脂中差異表達(dá)篩選獲得的2 個lncRNA，根據(jù)測序結(jié)果中提供的序列設(shè)計(jì)引物，SCD基因（GenBank 登錄號：NM_001009254.1）和GPAM基因（GenBank 登錄號：XM_012102928.3）mRNA 序列從NCBI 數(shù)據(jù)庫中獲得，利用Primer 5.0 軟件設(shè)計(jì)引物，其中以β-actin（ID：HM067830.1）作為內(nèi)參基因。引物名稱和序列、退火溫度以及擴(kuò)增片段大小見表1。

表1 熒光定量引物信息

1.5.2 cDNA 的合成 使用上一步提取的總RNA，按照PrimeScript RT reagent Kit with gDNA Eraser 試劑盒說明書進(jìn)行反轉(zhuǎn)錄，取出后的cDNA 在-20℃冷凍保存。

以上一步反轉(zhuǎn)錄獲得的cDNA 作為模板，進(jìn)行目的基因的PCR 擴(kuò)增（總體系為20 μL）：2×PCR Master Mix 10 μL，上、下游引物（0.2 μmol/L）各0.4 μL，cDNA 1.5 μL（100 ng），ddH2O 7.7 μL。擴(kuò)增程序?yàn)椋?4℃5 min；94℃ 30 s；54℃ 30 s，39 個循環(huán)；72℃ 5 min。PCR 產(chǎn)物取3 μL 與2 μL 6×Loading Buffer 混勻，用2%瓊脂糖凝膠電泳檢測，180 V 20 min；使用凝膠成像系統(tǒng)采集圖像，若獲得的產(chǎn)物片段大小與預(yù)期一致則進(jìn)行下一步實(shí)驗(yàn)。

1.5.3 實(shí)時熒光定量PCR 熒光定量檢測采用7500 Fast Real-Time PCR System 型熒光定量PCR 儀進(jìn)行，使用TB Green Premix Ex Taq II 的反應(yīng)體系進(jìn)行擴(kuò)增（總體系為20 μL）：TB Green Premix Ex Taq II 10 μL，RNase-freeddH2O 6 μL，上、下游引物各0.8 μL，Rox Reference Dye II 0.4 μL，cDNA 2 μL，具體實(shí)驗(yàn)操作流程按照產(chǎn)品說明書進(jìn)行。每個樣品基因表達(dá)以β-actin基因作為內(nèi)參基因，每個樣品重復(fù)檢測3 次。PCR 程序?yàn)椋?5℃ 30 s，95℃ 5 s，54℃ 30 s，40 個循環(huán)。

1.6 統(tǒng)計(jì)分析 采用2-ΔΔCt法計(jì)算目的基因的相對表達(dá)量。運(yùn)用SPSS 19.0 軟件的單因素方差分析（One-way ANOVA）進(jìn)行組間比較，并用最小顯著差異法（Least Significant Difference，LSD）進(jìn)行多重比較，P<0.05表示差異顯著，P<0.01 表示差異極顯著。

2 結(jié)果與分析

2.1 總RNA 的檢測 圖1 可明顯看到28S 和18S 條帶，5S 條帶不明顯，說明提取的總RNA 完整性好且無降解；利用核酸檢測儀進(jìn)行總RNA 濃度測定，提取的總RNA濃度高且無蛋白質(zhì)、DNA 等污染，可用于后續(xù)實(shí)驗(yàn)。

圖1 總RNA 電泳圖

2.2 PCR 產(chǎn)物檢測 利用2% 的瓊脂糖凝膠電泳檢測PCR 產(chǎn)物（圖2），發(fā)現(xiàn)產(chǎn)物片段大小均符合預(yù)期，有且僅有1 條帶，說明引物設(shè)計(jì)的特異性良好，前期RNA 提取及cDNA 的反轉(zhuǎn)錄符合要求。

圖2 目的基因的qRT-PCR 檢測

2.3 lncRNA 表達(dá)量與RNA-Seq 結(jié)果比較 選取綿羊尾脂樣品，將挑選出的lncRNA 的qRT-PCR 結(jié)果與前期RNA-Seq 結(jié)果進(jìn)行比較，結(jié)果表明：MSTRG.37980 與MSTRG.38164 在巴什拜羊與雜交二代羊尾部脂肪組織中的定量和RNA-Seq 結(jié)果基本一致，表明RNA-Seq 結(jié)果較為準(zhǔn)確可靠，可用于后續(xù)實(shí)驗(yàn)。

2.4 巴什拜羊和雜交二代羊不同組織中MSTRG.37980的表達(dá)差異分析 由圖3 可以看出，MSTRG.37980 在巴什拜羊尾脂和小腸中的表達(dá)量極顯著高于雜交二代羊，在背最長肌中的表達(dá)量顯著高于雜交二代羊，在皮下脂、臂三頭肌、股四頭肌、十二指腸、盲腸中的表達(dá)量均高于雜交二代羊，但差異不顯著。

圖3 MSTRG.37980 在巴什拜羊與雜交二代羊不同組織中的相對表達(dá)差異

2.5 巴什拜羊和雜交二代羊不同組織中SCD基因的表達(dá)差異分析 由圖4 可以看出，SCD基因在巴什拜羊尾脂中表達(dá)量極顯著高于雜交二代羊，在皮下脂、背最長肌、臂三頭肌、股四頭肌、十二指腸、盲腸、小腸部位的表達(dá)量高于雜交二代羊，但差異不顯著；雜交二代羊大腸中表達(dá)量顯著高于巴什拜羊，肝臟、瘤胃中的表達(dá)量高于巴什拜羊，但差異不顯著。

圖4 SCD 基因在巴什拜羊與雜交二代羊不同組織中的相對表達(dá)差異

2.6 巴什拜羊和雜交二代羊不同組織中的MSTRG.38164表達(dá)差異分析 由圖5 可知，MSTRG.38164 在巴什拜羊尾脂中表達(dá)量極顯著高于雜交二代羊，背最長肌和肝臟中表達(dá)量顯著高于雜交二代羊，且皮下脂、臂三頭肌、股四頭肌、十二指腸、盲腸、大腸、小腸中的表達(dá)量均高于雜交二代羊，差異不顯著；雜交二代羊瘤胃中的表達(dá)量高于巴什拜羊。

圖5 MSTRG.38164 在巴什拜羊與雜交二代羊不同組織中的相對表達(dá)差異

2.7 巴什拜羊和雜交二代羊不同組織中GPAM基因的表達(dá)差異分析 由圖6 可以看出，GPAM基因在巴什拜羊肝臟中表達(dá)量顯著高于雜交二代羊；雜交二代羊瘤胃中的表達(dá)量高于巴什拜羊，但差異不顯著，其他部位表達(dá)量均低于巴什拜羊。

圖6 GPAM 基因在巴什拜羊與雜交二代羊不同組織中的相對表達(dá)差異

3 討 論

3.1 MSTRG.37980 與SCD基因在巴什拜羊和雜交二代羊各組織表達(dá)量的比較分析SCD是在脂肪生物合成中將飽和脂肪酸（SFA）轉(zhuǎn)換為單不飽和脂肪酸（MUFA）的關(guān)鍵酶[11-12]，它具有調(diào)節(jié)動物體內(nèi)脂代謝的功能。前人研究發(fā)現(xiàn)，哺乳動物中，SCD主要在肝臟、肌肉組織、脂肪組織中高表達(dá)[13]，可作為肌內(nèi)脂肪沉積的潛在生物標(biāo)志物[14]。而本研究發(fā)現(xiàn)SCD基因的mRNA 在巴什拜羊皮下脂肪表達(dá)量最高，雜交二代羊大腸、肝臟表達(dá)量最高，與前人研究一致，在組織中的表達(dá)量情況基本相同，但SCD基因在巴什拜羊與二代羊各組織中的表達(dá)規(guī)律有一定差異。Huang 等[15]在萊蕪豬和大白豬皮下脂肪組織中鑒定出了與脂代謝和成脂分化相關(guān)的lncRNA 和基因，表明了lncRNA 可以靶向mRNA，然后在PPAR 信號傳導(dǎo)通路、脂肪細(xì)胞分化和脂肪酸代謝中發(fā)揮重要的作用，并鑒定了SCD基因參與PPAR 信號通路。本實(shí)驗(yàn)結(jié)果推測SCD基因通過參與PPAR 信號通路，參加脂肪合成相關(guān)基因的表達(dá)，使得巴什拜羊形成大尾。Ma 等[10]人通過綿羊尾部脂肪組織中的差異表達(dá)基因和lncRNA 進(jìn)行網(wǎng)絡(luò)構(gòu)建，共得到了493 對共表達(dá)對，并發(fā)現(xiàn)這些共表達(dá)的對中，大部分是顯著正相關(guān)，只有一小部分對是負(fù)相關(guān)，表明了這些mRNA和lncRNA 可能主要通過正向調(diào)控發(fā)揮作用。本實(shí)驗(yàn)結(jié)果MSTRG.37980 和SCD基因的mRNA 在巴什拜羊和雜交二代羊各組織部位的表達(dá)規(guī)律基本相似，說明MSTRG.37980 可能參與脂肪沉積，并且可能通過正向調(diào)控作用使SCD基因的mRNA 表達(dá)高。

3.2 MSTRG.38164 與GPAM基因在巴什拜羊和雜交二代羊各組織表達(dá)量的比較分析GPAM基因的主要作用是通過催化三?；视秃土字锖铣蛇^程，從而促進(jìn)動物機(jī)體甘油三酯（TAG）的生成，并通過中心和外周調(diào)節(jié)作用，對哺乳動物脂肪沉積、能量消耗、胴體質(zhì)量以及整體的代謝進(jìn)行調(diào)控[16-17]。于海濱[18]為了驗(yàn)證GPAM基因?qū)χx的影響，通過脂質(zhì)體介導(dǎo)牛GPAM基因的RNA 干擾載體，沉默內(nèi)源，抑制GPAM表達(dá)后顯著降低牛前體脂肪細(xì)胞和牛胎兒成纖維細(xì)胞中的甘油三酯含量，證明了GPAM基因可能在調(diào)節(jié)動物脂肪代謝中起關(guān)鍵作用。并且有學(xué)者發(fā)現(xiàn)，GPAM基因在所有脂肪中都有所表達(dá)，其中在肝臟、脂肪組織中表達(dá)量最高，其次是肌肉[19]。而本實(shí)驗(yàn)發(fā)現(xiàn)GPAM基因在肌肉、脂肪以及肝臟部位的表達(dá)量較高，與該研究結(jié)果基本一致，推測出GPAM可能調(diào)控脂肪合成。本實(shí)驗(yàn)對比了GPAM基因在巴什拜羊與二代羊相同部位的表達(dá)差異發(fā)現(xiàn)，該基因在2 種羊肝臟部位的表達(dá)差異顯著。由此推測，GPAM基因在巴什拜羊中可能主要在肝臟、脂肪部位發(fā)生作用，從而增加脂肪的沉積。本實(shí)驗(yàn)還發(fā)現(xiàn)GPAM基因和MSTRG.38164 在巴什拜羊和雜交二代羊各部位的表達(dá)規(guī)律是相似的，MSTRG.38164 可能通過對GPAM的調(diào)控，參與脂肪合成并沉積。以上研究結(jié)果表明，MSTRG.38164 和GPAM基因功能可能是廣泛的，尤其是在脂肪沉積方面。

4 結(jié) 論

本研究利用qRT-PCR 對2 種差異表達(dá)lncRNA 靶基因進(jìn)行了表達(dá)量分析，發(fā)現(xiàn)MSTRG.37980 和MSTRG.3定量與RNA-Seq 結(jié)果相似；SCD基因和MSTRG.37980在巴什拜羊和雜交二代羊各組織部位表達(dá)規(guī)律基本相似，并且SCD基因通過參與PPAR 信號通路，參加脂肪合成相關(guān)基因的表達(dá)，使得巴什拜羊形成大尾，MSTRG.37980 可能參與脂肪沉積，并通過正向調(diào)控作用使SCD基因的mRNA 表達(dá)高；GPAM基因和MSTRG.38164 在巴什拜羊和雜交二代羊各部位的表達(dá)規(guī)律相似，GPAM基因可能在巴什拜羊肝臟、脂肪部位發(fā)生作用，使得脂肪沉積增加。MSTRG.38164 可能通過對GPAM的調(diào)控參與脂肪合成并沉積。