哈薩克羊和藏綿羊尾部脂肪組織差異表達(dá)基因和可變剪接的鑒定與分析

陶海溪，李貴林，郭家中

(四川農(nóng)業(yè)大學(xué) 動物科技學(xué)院，成都 611130)

哈薩克羊和藏綿羊尾部脂肪組織差異表達(dá)基因和可變剪接的鑒定與分析

陶海溪，李貴林，郭家中

(四川農(nóng)業(yè)大學(xué) 動物科技學(xué)院，成都 611130)

為鑒定哈薩克羊和藏綿羊尾部脂肪組織的差異表達(dá)基因和可變剪接，分別提取6只(3公，3母)無親緣關(guān)系的2歲成年哈薩克羊與藏綿羊尾部脂肪組織總RNA，構(gòu)建2個高通量測序文庫，利用TopHat2軟件將短序列定位到綿羊參考基因組，采用cuffdiff軟件分析基因表達(dá)豐度和差異表達(dá)基因。結(jié)果表明，共有19個基因在 2 個品種綿羊尾脂中的表達(dá)豐度均高于1 000 FPKM，包括ADIPOQ、SPARC、VIM等。其中，2 個品種尾部脂肪組織中表達(dá)量最高的轉(zhuǎn)錄本均定位于線粒體基因組5 331～14 154 bp，該區(qū)域主要包含 COX1、COX2、ATP8等基因，其表達(dá)豐度在哈薩克羊和藏綿羊中分別高達(dá)29 724.20 和28 818.80 FPKM。而 FABP4則是表達(dá)豐度最高的核基因，在哈薩克羊和藏綿羊中的表達(dá)量分別為3 194.03和2 667.09。差異表達(dá)基因分析共鑒定出627個差異表達(dá)的基因，與藏綿羊尾部脂肪組織的表達(dá)量相比，哈薩克羊的尾部脂肪組織中共有221個上調(diào)基因和406個下調(diào)基因。其中，HBB基因表達(dá)量的變化倍數(shù)最大，下調(diào)約98倍。差異表達(dá)基因的富集分析結(jié)果表明，哈薩克羊尾脂中上調(diào)表達(dá)的基因主要參與脂肪酸代謝過程。另外，在2 個品種尾脂中鑒定出11 870個可變剪接事件，其中383個為差異剪接事件。對相關(guān)差異表達(dá)基因和可變剪接的進(jìn)一步研究將有助于揭示綿羊肥脂尾形成的分子機(jī)制。

肥脂尾；差異表達(dá)基因；哈薩克羊；藏綿羊；轉(zhuǎn)錄組

脂肪組織是動物體內(nèi)一種主要由脂肪細(xì)胞組成的能量貯存器官。一般地，脂肪主要在各種動物的皮下、腹部和內(nèi)臟器官周圍大量沉積。但在現(xiàn)有綿羊品種中，除上述部位外，許多品種綿羊的尾部具有超強(qiáng)的脂肪存儲能力，從而形成肥臀尾特征。例如，Safdarian等[1]報(bào)道伊朗肥脂尾型綿羊尾部的脂肪質(zhì)量可以達(dá)到其胴體質(zhì)量的20%。在中國的綿羊品種中，最為著名的肥臀尾類型綿羊包括大尾寒羊和哈薩克羊?？偟膩碚f，綿羊肥臀(脂)尾是綿羊長期在惡劣環(huán)境下生存，受自然選擇作用所形成的一種適應(yīng)性性狀。在面對營養(yǎng)匱乏、天氣寒冷等惡劣的自然環(huán)境時，脂尾型綿羊能夠消耗脂肪提供能量來維持自身的新陳代謝。但是，隨著養(yǎng)羊業(yè)的發(fā)展和人們需求的改善，綿羊肥脂尾逐漸喪失經(jīng)濟(jì)價(jià)值，綿羊尾部脂肪的過度沉積增加飼養(yǎng)成本，降低經(jīng)濟(jì)效益。因此，如何通過遺傳育種措施阻斷脂肪在綿羊尾部的沉積成為綿羊育種工作者關(guān)心的重要問題。

目前，研究人員從不同角度對綿羊肥臀尾形成的遺傳機(jī)制進(jìn)行探索。Moradi等[2]利用全基因組SNP芯片檢測到綿羊第5和X染色體 2 個與綿羊肥脂尾形成有顯著關(guān)聯(lián)的QTL區(qū)域。劉真等[3]針對不同尾型綿羊進(jìn)行群體間選擇信號檢測發(fā)現(xiàn)，PPARG基因在肥脂尾綿羊尾脂中的表達(dá)量顯著高于瘦尾綿羊。最近，Wang等[4]利用從頭組裝的分析策略，發(fā)現(xiàn) FMO3、 NELL1、THRSP等基因在哈薩克羊和藏綿羊尾部脂肪組織中呈現(xiàn)差異性表達(dá)。現(xiàn)有家畜遺傳資源研究表明，中國綿羊品種主要分為：藏綿羊、哈薩克羊和蒙古羊。其中，藏綿羊因長期生活在青藏高原地區(qū)，從而對高海拔低氧、寒冷的氣候環(huán)境和粗放的飼養(yǎng)方式具有良好的適應(yīng)能力。此外，藏綿羊具有典型的短尾特征。而主要分布在新疆地區(qū)的哈薩克羊則屬于肉脂兼用型綿羊品種。與其他綿羊品種相比，哈薩克綿羊最為突出的表型特征是具有巨大的肥臀尾。為深入了解不同尾型綿羊尾部脂肪沉積能力差異的遺傳機(jī)制，本研究利用基于參考基因組的短序列定位策略，鑒定哈薩克羊和藏綿羊尾部脂肪組織中基因的表達(dá)豐度特征、差異表達(dá)基因以及可變剪接，在轉(zhuǎn)錄水平上探討綿羊肥脂尾形成的分子機(jī)制。

1 材料與方法

1.1 動物組織樣本

在寧夏自治區(qū)某羊場分別隨機(jī)選擇6只(3 只公羊、3 只母羊)無親緣關(guān)系的 2 歲成年哈薩克羊與藏綿羊。分別提取12只綿羊尾部脂肪組織的總RNA，并對提取的RNA進(jìn)行質(zhì)量檢測。樣品檢測合格后，從每個品種6個個體的RNA樣品中分別選取等量體積的RNA進(jìn)行混合，形成2個RNA池(RNA-pool)；然后將混合后的 2 個RNA樣品送至深圳華大基因公司，構(gòu)建2個高通量測序文庫，采用Illumina Hiseq2000測序儀進(jìn)行高通量測序。

1.2 讀段比對和組裝

高通量測序后，共獲得117 417 578條100 bp長的原始雙末端讀段(pair-end reads)數(shù)據(jù)。為去除低質(zhì)量短序列，對原始讀段進(jìn)行質(zhì)量控制，具體參數(shù)包括：過濾含有接頭的reads； 過濾掉堿基中N>10%的reads；去除質(zhì)量值Q≤5 的堿基數(shù)占整個reads 50%以上的低質(zhì)量reads；獲得高質(zhì)量讀段(clean reads)。使用TopHat2 (v2.0.10)軟件[5]及默認(rèn)參數(shù)，將clean reads比對到綿羊參考基因組[6](ftp://ftp.ensembl.org/pub/release-77/fasta/ovis_aries/dna/)。采用MarkDuplicates程序 (http://sourceforge.net/projects/picard/)去除重復(fù)讀段，最后采用cufflinks[7](v2.2.1) 對上述讀段進(jìn)行組裝和注釋。

1.3 差異表達(dá)基因鑒定和功能富集分析

由于本研究未設(shè)置生物學(xué)重復(fù)，因此利用cuffdiff軟件[7](v1.3.0) 測定基因表達(dá)豐度并鑒定差異表達(dá)基因(differentially expressed genes，DGEs)。為校正測序深度和基因長度對基因表達(dá)量的影響，采用FPKM(fragment per kilobase of exon model per million mapped reads)作為基因表達(dá)量的單位。在本研究中，以0.000 5 FPKM作為基因表達(dá)豐度的閾值來判定基因是否在樣本組織中表達(dá)。差異表達(dá)基因的鑒定則采用以下標(biāo)準(zhǔn)： 2 個品種間基因表達(dá)倍數(shù)的對數(shù)變化大于1[(|log2Ratio| ≥1)]；假陽性率(FDR)小于0.001。采用上海博豪開發(fā)的生物學(xué)通路富集分析在線工具(http://prediction.ebioservice.com:8080/path/)對差異表達(dá)基因進(jìn)行通路富集分析，根據(jù)q-value<0.05確定顯著富集的生物學(xué)通路。

1.4 可變剪接的鑒定

目前，一般根據(jù)可變剪接(alternative splicing,AS)發(fā)生的位置將基因的可變剪接分成7類：外顯子跳躍(Exon skipping)、 5′端可變剪接位點(diǎn)(alternative 5′splice sites)、3′端可變剪接位點(diǎn)(alternative 3′splice sites)、外顯子互斥(mutually exclusive exons)、內(nèi)含子保留(retained introns)、可變的起始外顯子(alternative first exon)和可變的末端外顯子(alternative last exon)。總的來說，基于轉(zhuǎn)錄組數(shù)據(jù)的可變剪接分析軟件的分析原理是：利用跨越外顯子讀段(junction reads)提供的信息，將每組樣本數(shù)據(jù)分別與參考基因組注釋文件進(jìn)行比較，鑒定每組樣本中的可變剪接。與大部分可變剪接鑒定軟件[7-11]不同的是，ASD軟件[12]的分析策略是將所有待分析的轉(zhuǎn)錄組樣本數(shù)據(jù)整合起來，首先鑒定待分析樣本與參考基因組的可變剪接，最終著重于鑒定不同組數(shù)據(jù)之間發(fā)生的差異可變剪接。該軟件更適合于不同品種的相同組織，或不同處理下相同組織的轉(zhuǎn)錄組數(shù)據(jù)。因此，本研究采用ASD軟件(http://erp.novelbio.com/ASD/)鑒定可變剪接。

2 結(jié)果與分析

2.1 讀段的基因組比對

由表1可知，哈薩克羊尾部脂肪組織轉(zhuǎn)錄組文庫共包含25 766 605對(51 533 210條)clean reads，約4.6 G序列數(shù)據(jù)，其平均GC含量50.39%；藏綿羊尾部脂肪組織轉(zhuǎn)錄組文庫共包含27 747 038對(55 494 076 條)clean reads，約5.0 G序列，其平均GC含量49.19%。clean reads的基因組比對結(jié)果表明，在哈薩克羊和藏綿羊的文庫中，分別有21 554 938和23 442 091對reads比對到綿羊參考基因組，分別占總clean reads對的83.65%和84.49%。其中，唯一匹配讀段數(shù)量分別為18 715 395和19 953 180對clean reads，分別占總clean reads的72.63% 和71.91%。

2.2 尾部脂肪組織高表達(dá)基因

根據(jù)表達(dá)量≥0.000 5 FPKM的閾值設(shè)定，本研究分別在哈薩克羊和藏綿羊尾部脂肪組織中檢測到26 173個和26 345個表達(dá)的基因。依據(jù)FPKM值，本研究進(jìn)一步將樣本組織中表達(dá)的基因分為低豐度表達(dá)基因(< 10 FPKM)，中等豐度表達(dá)基因(10 FPKM～500 FPKM)和高豐度表達(dá)基因(≥500 FPKM)。在哈薩克羊和藏綿羊的尾部脂肪組織中，表達(dá)量小于10 FPKM的基因數(shù)目分別占所有表達(dá)基因數(shù)目的70.02% 和67.77%，表明尾部脂肪組織中大部分基因都處于低豐度表達(dá)。而哈薩克羊和藏綿羊尾部脂肪組織中高表達(dá)(≥500 FPKM)的基因數(shù)目分別為150和165。如表2所示，在 2 個品種脂肪組織中表達(dá)量均大于1 000 FPKM的基因共19個，包括ADIPOQ、SPARC、VIM、CYTB、 FTH1、 GPX3、 RPS11、 RPS3A和GSN等基因。其中，在 2 個品種中表達(dá)量最高的2 個轉(zhuǎn)錄本均定位于綿羊線粒體基因組5 331～14 154 bp和0～5 018。上述 2 個區(qū)域在哈薩克羊和藏綿羊中的表達(dá)豐度為29 724.20 和28 818.80 FPKM，5 476.07和5 244.22 FPKM；這 2 個區(qū)域主要包括 COX1、 COX2、 ATP8等基因和 ND1、ND2等基因。在定位到綿羊基因組的基因中，脂肪酸結(jié)合蛋白-4( FABP4)的表達(dá)豐度最高，其在哈薩克羊和藏綿羊中的表達(dá)量分別為3 194.03和2 667.09 FPKM，而FABPs家族的 FABP5也表現(xiàn)出中等豐度的表達(dá)量(244.02 和373.41 FPKM)，其他成員均呈現(xiàn)低豐度表達(dá)或不表達(dá)。另外，脂聯(lián)素基因(ADIPOQ)在 2 個品種的尾部脂肪組織中表達(dá)量也分別高達(dá)2 645.49和1 600.51 FPKM。4個核糖體蛋白基因( RPS11、 RPS3A、 RPS17和 RPS14)在 2 個品種中的表達(dá)量大于1 000 FPKM。

表1 2 個文庫clean reads的比對Table 1 Mapping results of clean reads in both libraries

注：百分比表示各種類型讀段占高質(zhì)量讀段(clean reads)的百分比。

Note:Percentages in the table present percentages of diferent reads in clean read pairs.

表2 2 個綿羊品種表達(dá)量高于1 000 FPKM的基因Table 2 Expressed genes showing larger than 1 000 FPKM in both breeds

2.3 哈薩克羊和藏綿羊差異表達(dá)基因的鑒定及功能注釋

差異表達(dá)基因的分析結(jié)果表明，相對于藏綿羊尾部脂肪中的表達(dá)豐度，共有627個基因在哈薩克羊尾部脂肪組織中的表達(dá)量呈現(xiàn)顯著差異，包括221個上調(diào)和406個下調(diào)表達(dá)的基因。如表3所示，在哈薩克羊中，下調(diào)的高豐度表達(dá)基因包括HBB(血紅蛋白β)、 ENSOARG00000020789、FNDC1、TGFBR3、CXCL12等基因；而上調(diào)的高豐度表達(dá)基因則包括Carboxylicesterhydrolase(羧基酯水解酶)和CAT(過氧化氫酶)、THRSP、HSL、 FMO3等基因。在所有差異表達(dá)的基因中，HBB基因表達(dá)量變化倍數(shù)最大；與藏綿羊中的表達(dá)量(791.43 FPKM)相比，哈薩克羊尾部脂肪組織中HBB的表達(dá)量(8.08 FPKM)下降約 98倍。 差異表達(dá)基因的GO生物學(xué)過程(biological process,BP)富集分析結(jié)果表明，在哈薩克羊中上調(diào)表達(dá)的基因顯著富集在20個生物學(xué)過程。如圖1所示，上調(diào)基因顯著富集的過程主要包括9個組織水平或細(xì)胞水平的甘油三酯或脂類物質(zhì)代謝過程，例如，脂類代謝過程(lipid metabolic process， GO:0006629)、甘油三酯催化過程(triglyceride catabolic process，GO:0019433)、細(xì)胞脂類物質(zhì)代謝過程(cellular lipid metabolic process， GO:0044255)等。這些過程主要涉及的差異表達(dá)基因包括：MGLL、AADAC、GPLD1、LIPE、PLIN5等。差異表達(dá)基因的GO富集分析結(jié)果表明，在哈薩克羊中下調(diào)表達(dá)的基因顯著富集在139個生物學(xué)過程。如圖2所示，最顯著富集的20個生物學(xué)過程中主要涉及細(xì)胞水平的信號傳遞與轉(zhuǎn)導(dǎo)、免疫應(yīng)答過程。

2.4 可變剪接事件的檢測

如表4所示，通過與綿羊參考基因組比較，本研究在 2 個品種的尾部脂肪組織中共鑒定到11 870 個可變剪接事件。其中， 3′端可變剪接位點(diǎn)和外顯子缺失分別發(fā)生3 622次和3 232次，所占比例最高；而外顯子互斥、可變的起始外顯子和可變的末端外顯子分別發(fā)生621、1 157和689次，發(fā)生比例較低。根據(jù)卡方檢驗(yàn)(FDR<0.05)，共檢測到383個在 2 個組織間存在的差異可變剪接，其中類型最多的是外顯子跳躍(157次)，涉及132個基因。

圖1 哈薩克羊中上調(diào)表達(dá)基因最顯著富集的20個GO 生物學(xué)過程Fig.1 Top 20 significantly enriched GO biological processes for up-regulated genes in Kazak sheep

圖2 哈薩克羊中下調(diào)表達(dá)基因最顯著富集的20個GO 生物學(xué)過程Fig.2 Top 20 significantly enriched GO biological processes for up-regulated genes in Kazak sheep

表4 2 個品種尾部脂肪組織中不同類型可變剪接的數(shù)目Table 4 Summary of different alternative splicing (AS) detected in tails of two breeds

3 討 論

已有研究充分證實(shí)[13-14]，線粒體在前脂肪細(xì)胞分化和成熟脂肪細(xì)胞甘油三酯的合成與水解中扮演著重要作用。例如，細(xì)胞內(nèi)活性氧含量的升高將影響線粒體ATP合成酶的生成，從而抑制3T3-L1細(xì)胞的增殖[15]。線粒體拷貝數(shù)則與人類白色脂肪組織中脂肪生成過程顯著關(guān)聯(lián)[16]。本研究發(fā)現(xiàn)，包括 COX1、COX2、ATP8和CYTB在內(nèi)的線粒體基因，在哈薩克羊和藏綿羊尾部脂肪組織均為表達(dá)豐度最高的轉(zhuǎn)錄本，反映出綿羊尾部脂肪組織的代謝非常活躍。脂肪酸結(jié)合蛋白家族是脂肪合成與水解過程中重要的調(diào)控分子，但是各成員之間的時空表達(dá)特征與功能有差異。例如， FABP1基因主要在肝臟中高表達(dá)， FABP2基因主要在小腸中高表達(dá)；而 FABP4基因則是在脂肪組織中特異性高度表達(dá)，促進(jìn)脂肪組織中脂肪水解，同時抑制脂肪合成[17]。本研究的結(jié)果再次證明上述觀點(diǎn)： FABP4在 2 個綿羊品種的脂肪組織中均高水平表達(dá)，同時 FABP5也表現(xiàn)出中等豐度的表達(dá)，表達(dá)量分別為244.02 和373.41 FPKM，而FABPs家族其他成員均呈現(xiàn)低水平表達(dá)或不表達(dá)。另外，本研究中 4 個核糖核酸蛋白基因的高豐度表達(dá)表明，綿羊尾部脂肪組織在翻譯水平的基因表達(dá)是非常活躍的。

差異表達(dá)基因的鑒定能夠幫助揭示不同類型綿羊尾部脂肪沉積能力差異的分子機(jī)制。血紅蛋白β基因(HBB)與動物體攜帶和運(yùn)輸氧的能力密切相關(guān)，被認(rèn)為是高原人類和動物的高寒適應(yīng)性進(jìn)化的一個重要候選基因[18-19]。本研究發(fā)現(xiàn)，HBB基因在藏綿羊尾部脂肪組織中的表達(dá)量較高，但在哈薩克羊尾部脂肪組織中的表達(dá)量則顯著下調(diào)。該結(jié)果表明，HBB基因可能在藏綿羊脂肪組織的基礎(chǔ)代謝中具有重要功能。已有研究[20]表明，動物脂肪組織中脂類物質(zhì)代謝過程會產(chǎn)生大量的活性氧類物質(zhì)。而過氧化氫酶的主要作用是分解過氧化氫，從而阻止后者在生物體內(nèi)進(jìn)一步分解成有害的自由基。在本研究中，過氧化氫酶基因在哈薩克羊尾部脂肪中表達(dá)量升高，這可能是由于哈薩克羊尾部脂肪組織需要清除高度活躍的代謝過程所產(chǎn)生的大量過氧化氫。激素敏感酯酶(HSL)是動物體脂肪分解過程中的一個關(guān)鍵酶，其主要負(fù)責(zé)水解甘油三酯、甘油二酯以及膽固醇酯，HSL基因的缺失能導(dǎo)致甘油二酯在小鼠脂肪和肌肉組織中的堆積[21]，從而影響脂肪水解。本研究發(fā)現(xiàn)，HSL在哈薩克羊尾部脂肪組織中的表達(dá)量顯著高于其在藏綿羊中的表達(dá)豐度，暗示著哈薩克羊肥脂尾中的脂肪分解比藏綿羊尾部脂肪中的代謝更為活躍?？偟膩碚f，與藏綿羊相比，在哈薩克羊尾部脂肪組織中顯著上調(diào)的基因主要富集在脂肪代謝相關(guān)的生物學(xué)過程中，表明哈薩克羊尾部組織中脂肪代謝比藏綿羊中的脂肪代謝顯著活躍。而下調(diào)表達(dá)的基因所富集的一些通路則與免疫過程相關(guān)。與Wang等[4]利用從頭組裝分析策略檢測到的差異表達(dá)基因相比，本研究共檢測到60個相同的上調(diào)表達(dá)基因，46個相同的下調(diào)表達(dá)基因。相同的上調(diào)表達(dá)基因包括THRSP、FMO3等，這些基因在肥臀尾綿羊的尾部脂肪沉積中的作用值得進(jìn)一步研究。

可變剪接廣泛存在于動物各種組織的基因表達(dá)過程中[22-23]，并被認(rèn)為是導(dǎo)致動植物基因轉(zhuǎn)錄調(diào)控復(fù)雜性的一個重要原因。已有研究表明，動物脂肪組織中存在大量可變剪接并且是脂肪細(xì)胞分化的必要條件[24]，剪接因子 SRSF10基因的缺失將阻礙前脂肪細(xì)胞的分化，從而影響皮下白色脂肪的發(fā)育過程?？勺兗艚拥陌l(fā)生與動物的肥胖和胰島素抗性有密切關(guān)系[25]。例如， LPIN1基因的一種剪接體 LPIN1 α 主要在前脂肪細(xì)胞的分化中發(fā)揮作用，而另一種剪接體 LPIN1 β則會促進(jìn)前脂肪細(xì)胞增殖[26]。本研究發(fā)現(xiàn)，與綿羊參考基因組相比， LPIN1 基因也存在可變剪接，其剪接類型為5′端可變剪接。本研究的結(jié)果表明，綿羊脂肪組織中存在著大量的可變剪接。在所檢測到的不同類型可變剪接中，3′端可變剪接和外顯子跳躍發(fā)生的比例最高，這與已有研究[24]相一致，另外，有多個基因(例如 C4)的可變剪接在 2 個品種脂肪組織中存在差異，表明可變剪接可能在不同尾型綿羊尾部脂肪沉積中發(fā)揮重要作用。但相關(guān)基因的可變剪接準(zhǔn)確位置及形成的不同轉(zhuǎn)錄本作用還有待于通過試驗(yàn)進(jìn)一步證實(shí)。

Reference：

[1] SAFDARIAN M,ZAMIRI M J,HASHEMI M,etal.Relationships of fat-tail dimensions with fat-tail weight and carcass characteristics at different slaughter weights of Torki-Ghashghaii sheep[J].MeatScience,2008,80(3):686-689.

[2] MORADI M H,NEJATI-JAVAREMI A,MORADI-SHA-HRBABAK M,etal.Genomic scan of selective sweeps in thin and fat tail sheep breeds for identifying of candidate regions associated with fat deposition[J].BMCGenetics,2012,13(1):1.

[3] 劉 真,王慧華,劉瑞鑿,等.不同尾型綿羊全基因組選擇信號檢測[J].畜牧獸醫(yī)學(xué)報(bào),2015,46(10):1721-1732.

LIU ZH,WANG H H,LIU R Z,etal.Genome-wide detection of selection signatures of distinct tail types in sheep populations[J].ChineseJournalofAnimalandVeterinarySciences,2015,46(10):1721-1732.(in Chinese with English abstract).

[4] WANG X L,ZHOU G X,XU X CH,etal.Transcriptome profile analysis of adipose tissues from fat and short-tailed sheep[J].Gene,2014,549(2):252-257.

[5] KIM D,PERTEA G,TRAPNELL C,etal.TopHat2:accurate alignment of transcriptomes in the presence of insertions,deletions and gene fusions[J].GenomeBioogyl,2013,14(4):R36.

[6] JIANG Y,XIE M,CHEN W B,etal.The sheep genome illuminates biology of the rumen and lipid metabolism[J].Science,2014,344(6188):1168-1173.

[7] TRAPNELL C,WILLIAMS B A,PERTEA G,etal.Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation[J].NatureBiotechnology,2010,28(5):511-515.

[8] BROOKS A N,YANG L,DUFF M O,etal.Conservation of an RNA regulatory map between Drosophila and mammals[J].GenomeResearch,2011,21(2):193-202.

[9] KATZ Y,WANG E T,AIROLDI E M,etal.Analysis and design of RNA sequencing experiments for identifying isoform regulation[J].NatureMethods,2010,7(12):1009-1015.

[10] SHEN S,PARK J W,LU Z X,etal.rMATS:robust and flexible detection of differential alternative splicing from replicate RNA-Seq data[J].ProceedingsoftheNationalAcademyofSciences,2014,111(51):E5593-5601.

[11] LIU R,LORAINE A E,DICKERSON J A.Comparisons of computational methods for differential alternative splicing detection using RNA-seq in plant systems[J].BMCBioinformatics,2014,15(1):364.

[12] ZHOU X X,WU W W,LI H,etal.Transcriptome analysis of alternative splicing events regulated by SRSF10 reveals position-dependent splicing modulation[J].NucleicAcidsResearch,2014,42(6):4019-4030.

[13] DE PAUW A,TEJERINA S,RAES M,etal.Mitochondrial (dys) function in adipocyte (de) differentiation and systemic metabolic alterations[J].TheAmericanJournalofPathology,2009,175(3):927-939.

[14] MEDINA-GMEZ G.Mitochondria and endocrine function of adipose tissue[J].BestPractice&ResearchClinicalEndocrinology&Metabolism,2012,26(6):791-804.

[15] CARRIRE A,FERMANDEZ Y,RIGOULET M,etal.Inhibition of preadipocyte proliferation by mitochondrial reactive oxygen species[J].FEBSLetters,2003,550(1/3):163-167.

[16] KAAMAN M,SPARKS L,VAN HARMELEN V,etal.Strong association between mitochondrial DNA copy number and lipogenesis in human white adipose tissue[J].Diabetologia,2007,50(12):2526-2533.

[17] HOTAMISLIGIL G S,BERNLOHR D A.Metabolic functions of FABPs—mechanisms and therapeutic implications[J].NatureReviewsEndocrinology,2015,11:592-605.

[18] WANG G D,FAN R X,ZHAI W,etal.Genetic convergence in the adaptation of dogs and humans to the high-altitude environment of the tibetan plateau[J].GenomeBiologyandEvolution,2014,6(8):2122-2128.

[19] KATZ Y,WANG E T,AIROLDI E M,etal.Analysis and design of RNA sequencing experiments for identifying isoform regulation[J].NatureMethods,2010,7(12):1009-1015.

[20] MATSUDA M,SHIMOMURA I.Increased oxidative stress in obesity:implications for metabolic syndrome,diabetes,hypertension,dyslipidemia,atherosclerosis,and cancer[J].ObesityResearch&ClinicalPractice,2013,7(5):e330-e341.

[21] HAEMMERLE G,ZIMMERMANN R,HAYN M,etal.Hormone-sensitive lipase deficiency in mice causes diglyceride accumulation in adipose tissue,muscle,and testis[J].JournalofBiologicalChemistry,2002,277(7):4806-4815.

[22] PAN Q,SHAI O,LEE L J,etal.Deep surveying of alternative splicing complexity in the human transcriptome by high-throughput sequencing[J].NatureGenetics,2008,40(12):1413-1415.

[23] THATCHER S R,ZHOU W,LEONARD A,etal.Genome-wide analysis of alternative splicing in Zea mays:landscape and genetic regulation[J].PlantCell,2014,26(9):3472-3487.

[24] LI H,CHENG Y M,WU W W,etal.SRSF10 regulates alternative splicing and is required for adipocyte differentiation[J].MolecularandCellularBiology,2014,34(12):2198-2207.

[25] KAMINSKA D,PIHLAJAMKI J.Regulation of alternative splicing in obesity and weight loss[J].Adipocyte,2013,2(3):143-147.

[26] PTERGY M,PHAN J,REUE K.Alternatively spliced lipin isoforms exhibit distinct expression pattern,subcellular localization,and role in adipogenesis[J].JournalofBiologicalChemistry,2005,280(38):32883-32889.

(責(zé)任編輯：顧玉蘭 Responsible editor:GU Yulan)

Identification and Analysis of Differentially Expressed Genes and Alternative Splicing in Adipose Tissue from Tails between Kazak and Tibetan Sheep

TAO Haixi,LI Guilin and GUO Jiazhong

(College of Animal Science and Technology,Sichuan Agricultural University,Chengdu 611130,China)

This study aims to investigate the expression profiles and alternative splicing patterns of the genes expressed in the adipose from tails of Kazak (KS) and Tibetan (TS) sheep.Six two-year old unrelated adult individuals (three males and three females) for KS and TS breeds were randomly selected from a sheep farm located in Ningxia Autonomous Region,respectively.Total RNA from six sampled individuals of each breeds were extracted and mixed into a RNA-pool which was then used to construct the high-throughput sequencing library.A total of 19 genes showed high expression levels with more than 1 000 FPKM,such asADIPOQ,SPARC,andVIM.The most abundantly expressed transcripts which mainly included COX1,COX2.The ATP8 was mapped to 5 331-14 154 bp of the mitochondrial genome in both breeds,and the expression abundance of this transcript were 29 724.20 and 28 818.80 FPKM in KS and TS,respectively.The expression levels of FABP4 most highly expressed nuclear gene were 3 194.03 FPKM in KS breeds and 2 667.09 in TS breeds,respectively.Comparing TS group with KS group,a total of 627 differentially expressed genes (DEGs) were identified,including 221 up-regulated and 406 down-regulated genes.Notably,HBBshowed largest down-regulation expression change in the fat tail of Kazak sheep.GO enrichment analysis revealed that the up-regulated expressed genes in KS sheep were significantly related with fatty acid metabolism processes.In addition,we detected a total of 11 870 alternative splicing events in the two breed samples,of which 383 splicing events were significantly discrepant between the two breeds.The further studies for the differentially expressed genes and alternative splicing will better understand the molecular mechanisms underlying sheep fat-tail.

Fat-tail; Differentially expressed genes; Kazak sheep; Tibetan sheep; Transcriptome

TAO Haixi,male,master student.Research area:sheep and goat genetics and breeding.E-mail:thx0218@sina@qq.com

GUO Jiazhong,male,Ph.D,lecturer.Research area:sheep and goat genetics and breeding.E-mail:jiazhong.guo@sicau.edu.cn

日期：2017-03-03

網(wǎng)絡(luò)出版地址：http://kns.cnki.net/kcms/detail/61.1220.S.20170303.0836.074.html

2016-04-26

2016-06-07

四川省教育廳項(xiàng)目(15ZB0005)。

陶海溪，男，碩士研究生，從事動物遺傳育種研究。E-mail:thx0218@sina.com

郭家中，男，博士，講師，主要從事草食家畜遺傳育種研究。E-mail:jiazhong.guo@sicau.edu.cn

S826

A

1004-1389(2017)03-0342-08

Received 2016-04-26 Returned 2016-06-07

Foundation item Programs supported by Department of Education of Sichuan Povince(No.15ZB0005).