1、綿羊是羊肉的重要來源，隨著羊肉價格的逐漸盤升，如何提高羊只的產肉量、改善羊肉品質成為亟待解決的問題。研究不同生長性能綿羊骨骼肌間的全轉錄組和差異表達將有助于揭示肌肉生長和發(fā)育的分子機制。
　　近年來，基于第二代高通量測序平臺的轉錄組測序(RNA-Seq)可以實現(xiàn)快速、全面地獲得特定組織于某個狀態(tài)下幾乎所有的轉錄本。本試驗以生長性能顯著不同的兩只小尾寒羊和杜泊羊的臂二頭肌為試驗材料，分別構建了它們的轉錄組文庫，并采用Illumina

2、 HiSeq2000高通量測序平臺進行測序;獲得的測序序列與綿羊參考基因組和參考基因比對后全面揭示兩個轉錄組文庫的特征，包括基因表達、基因注釋、差異基因、可變剪接、新轉錄本、cSNP和SSR等;以qRT-PCR法對從中隨機選取的12基因進行了相對定量分析對高通量測序數(shù)據(jù)的可靠性進行驗證;采用RACE等方法克隆了小尾寒羊和杜泊羊MYL1、MYL2、MYL3和MYL4基因的全長cDNA序列，利用GENSCAN等生物軟件對cDNA序列以及編碼

3、蛋白的結構特征等方面加以分析，并對這幾個基因在綿羊的不同組織間mRNA表達譜進行分析。主要研究結果如下:
　　(1)經RNA-Seq分析，分別在小尾寒羊和杜泊羊的兩個轉錄組文庫中得到50264608和52794216條長為90 bp高質量的測序序列。這些序列中各有約三分之二以上的序列可以比對到綿羊的參考基因組中，而有42.77％和33.10％的序列被可被比對到綿羊的20236個參考基因上。在綿羊臂二頭肌中，只有約1％的參考基因屬于

4、高表達基因(RPKM≥1000)，而93％以上的基因為低表達基因(RPKM≤1000)。
　　(2)兩個轉錄組文庫間發(fā)現(xiàn)有1300個差異表達基因(FDR≤0.001且| log2Ratio|≥1)，其中有554個基因在小尾寒羊中表達上調而746個表達下調。對它們進行GO功能注釋時發(fā)現(xiàn)分別有1066、1137和1067個基因可被注釋到GO的生物學過程、細胞組分和分子功能條目中。而1152個差異表達基因被注釋到了KEGG數(shù)據(jù)庫的240

5、條通路中，其中以新陳代謝通路（含有114個基因）和肌動蛋白細胞骨架通路（21個基因）包含基因最多。綜合兩方面的注釋結果共發(fā)現(xiàn)有31個差異表達基因與肌細胞發(fā)育、分化和骨骼肌的生長相關。蛋白互作網(wǎng)絡分析發(fā)現(xiàn)差異表達基因中的細胞外基質結構組分——整合蛋白ITGA8和膠原蛋白COL4A1與其他蛋白間互作最多，而肌球蛋白(myosin)、肌鈣蛋白(troponin)以及一些與生長和發(fā)育密切相關的調節(jié)因子（IGFBP5、TGFBR3和VEGF）與其

6、他基因的互作較少。
　　(3)在小尾寒羊和杜泊羊文庫中分別存在著40481和38851個潛在的cSNPs以及4721個SSRs，并分別以A→G和(AC)n為主要的分子標記。其中cSNPs主要分布于綿羊的1、2和3號染色體上。
　　(4)小尾寒羊和杜泊羊文庫中分別有37.72％和39.03％表達的基因發(fā)生了可變剪接，并以A3SS為主要的剪接方式。依照測序序列對6989個參考基因的兩端或一端進行了延伸和優(yōu)化。共發(fā)現(xiàn)了123678

7、個平均長為343 bp的新轉錄本單元。
　　(5)克隆得到了綿羊MYL1、MYL2、MYL3和MYL4基因的全長cDNA序列。其中MYL1基因包括長度不同的兩條cDNA序列，二者在5'端有所不同而3'端的序列相同。將這些序列提交至GenBank，登錄號分別為KJ700419、KJ710701、KJ710702、KJ710703、KJ710704、KJ710705、KJ710706、KJ710707和KJ768855。同一個基因的c

8、DNA序列在兩個品種間的同源性均高于與NCBI中預測序列的同源性，主要表現(xiàn)為5'端和3'端長度和序列的不同。
　　(6)綿羊MYL1、MYL2、MYL3和MYL4基因均包括5個外顯子，它們的起始密碼子上游5'側翼序列中含有大量肌細胞增加因子MEF2以及生肌調節(jié)因子MyoD和MyoG的結合位點。
　　(7)綿羊MYL1a、MYL1b、MYL2、MYL3和MYL4蛋白中只有1～2個氨基酸與NCBI中的預測序列不同。序列分析結果表

9、明它們均為偏酸性蛋白、親水性較好、無信號肽、有多個糖基化和磷酸化位點;二級結構以α-螺旋和無規(guī)卷曲為主;三級結構以人肌球蛋白輕鏈Ⅱ為最佳模板，兩端呈對稱性桶狀;聚類結果顯示脊椎動物中綿羊與山羊、人與大熊貓、灰鼠猴子與獼猴、小鼠與褐家鼠、牛和牦牛、貓與狗相近，而斑馬魚和爪蟾與以上脊椎動物相距最遠。
　　(8)綿羊MYL1基因主要在背最長肌中表達，兩種可變剪接體MYL1a在杜泊羊中表達量高于小尾寒羊，而MYL1b在小尾寒羊中表達量高于