多浪羊與卡拉庫爾羊FSHR 基因啟動子區甲基化水平差異研究

劉 波，芮 雪，方 翟，黃 飛，陶維昆，趙建清，李宏健，高慶華,3*

（1.新疆塔里木大學動物科學學院，新疆阿拉爾 843300；2.新疆塔里木大學生命科學學院，新疆阿拉爾 843300；3.新疆生產建設兵團塔里木畜牧科技重點實驗室，新疆阿拉爾 843300；4.新疆烏魯木齊畜牧獸醫總站，新疆烏魯木齊 830000）

多浪羊有較高的繁殖能力，性成熟早，且為常年發情，一般兩年產三胎，而且雙羔率較高，繁殖成活率在220%左右；卡拉庫爾羊繁殖率較低，一般發情受季節影響，產羔率在110%左右[1]。生殖激素受體基因是影響繁殖率的因素之一，促卵泡素受體（FSHR）基因介導促卵泡激素（FSH）作用于卵巢，是卵泡發育、成熟以及最終引發排卵的必需物質[2]。研究FSHR基因如何加快卵泡的生長分化、促進優勢卵泡的形成以及提高低繁殖家畜的繁殖率具有重要意義。王惠娥等[3]研究多浪羊和卡拉庫爾羊卵巢FSHR基因發現多浪羊對FSH的敏感性高于卡拉庫爾羊，表明FSHR 含量越高繁殖率越高。早期的研究表明缺乏FSH 或FSH 受體的卵泡未能進展到排卵前階段，導致不孕[4]。目前認為，FSHR基因的組織特異性表達與轉錄因子的作用、DNA 甲基化、染色體的重建、染色體組蛋白的修飾及mRNA、蛋白的轉換有關[5]，其中DNA 甲基化對基因表達水平的調控尤其受到關注。Priya 等[6]對大鼠FSHR基因啟動子研究中觀察到了FSHR基因啟動子區甲基化帶，且研究DNA 甲基化可能有助于更好地理解綿羊生殖能力的表觀遺傳調控。DNA 甲基化是調節FSHR基因細胞特異性沉默的主要因素之一，但是目前DNA 甲基化與FSHR基因表達的研究尚處于早期，還不完全了解其調控機制。本實驗通過亞硫酸氫鹽測序（Bisulfite Genomic se-quencing PCR，BSP）技術檢測不同綿羊FSHR基因啟動子區DNA 甲基化水平差異來探究綿羊FSHR基因啟動子區DNA 甲基化水平與繁殖性狀的關系，為進一步研究綿羊FSHR基因啟動子區DNA 甲基化水平與FSHR基因表達量的關系提供重要實驗依據。

1 材料與方法

1.1 材料

1.1.1 主要試劑 血液/ 組織/ 細胞基因組提取試劑盒（DP304）、RNaseA（100 mg/mL）溶液（RT405-12）、DH5α感受態細胞（CB101）、紅細胞裂解液（RT122）、10×Loading Buffer、DL 2000（8000）Marker 分子量標記、快捷性瓊脂凝膠DNA 回收試劑盒、pMD19-T 載體、Hot Start Taq 聚合酶、瓊脂、瓊脂糖、乙二胺四乙酸二鈉（EDTA）、無水乙醇等均購自天根生化科技有限公司。甲基化處理試劑盒（D5005）購置于ZYMO RESEARCH 公司。

1.1.2 主要儀器 PCR 儀（BIO RAD，C100）；核酸濃度檢測儀（Thermo，ND-One）；-80℃超低溫冰箱（Thermo，ExF24086V）；低溫高速冷凍離心機（Sigma，3K30）；高壓蒸汽滅菌鍋（TOMY，SX-500）；超凈工作臺（Boxum，SHC-CT-22）；恒溫水浴鍋（上海博訊，HHS-21-4）；電泳儀；凝膠成像儀。

1.1.3 試劑配制 試劑配制均以滅菌蒸餾水為溶劑，高壓滅菌條件：高壓蒸汽滅菌30 min。其余參考《分子克隆實驗指南》第三版李載平等[7]。

1.2 實驗方法

1.2.1 樣品的采集 本實驗使用的動物分別為多浪羊和卡拉庫爾母羊各3 只，均為1 歲母羊，且同品種間無血緣關系。血液采自塔里木大學動物實驗站，使用一次性采集容器（EDTA-K2）靜脈采血后送回實驗室-20℃保存備用。

1.2.2 基因組DNA 的提取與檢測 使用試劑盒DP304按照步驟提取不同綿羊血液組織基因組DNA。提取前所用器具高溫高壓滅菌并做烘干處理；使用DNA 提取試劑盒時，應該按照使用說明提取DNA；提取完成后使用1%的瓊脂糖凝膠進行電泳檢測質量，使用凝膠成像儀觀察基因組DNA 的片段大小及亮度情況，使用核酸濃度檢測儀檢測提取的綿羊血液組織DNA 濃度和純度，并拍照記錄。

1.2.3 基因啟動區引物設計合成與預測 參考綿羊FSHR基因序列（Gene ID：443299），根據FSHR基因啟動子區，使用MethPrimer 在線網站（www.bioon.com.cn/sub/showarticle.asp）預測FSHR基因啟動子序列存在的CpG 島檢測區域是否符合實驗要求，根據符合要求的區域設計FSHR基因啟動子區的引物，上游引物：5'-GTGTTTTTGTTTGGAGAATTTTAGG-3'，下游引物：5'-TCAATCAAAAAACCATCTACTTTCA-3'，送往廣州伯信生物科技有限公司進行合成。

1.2.4 基因啟動子區的擴增及回收純化 以BSP 處理后的綿羊血液DNA 為模板，對FSHR基因設置最佳PCR擴增體系，PCR 擴增程序設置為95℃預變性5 min；94℃變性30 s，68～53℃退火，每個循環降3℃，10 個循環；94℃變性30 s，51℃退火30 s，72℃延伸30 s，循環30 次；72℃終延伸10 min；4℃保存備用。擴增完成后，2%瓊脂糖凝膠電泳檢測，舍棄擴增失敗引物，并拍照做好記錄。按瓊脂糖凝膠回收純化試劑盒說明進行產物回收純化。

1.2.5 T-A 克隆及熱激轉化 將綿羊FSHR基因啟動子區PCR 擴增后回收產物與pMD19-T 載體連接，加DH5α感受態細胞進行單克隆篩選。

1.2.6 陽性克隆鑒定 過夜培養后會產生藍白色菌斑。吸取1 mL 的LB 液體培養基加入到1.5 mL 離心管內，用滅菌處理過的牙簽或者槍頭挑取飽滿的白色菌落放于離心管中，放置于恒溫搖床（37℃，200 rpm）中振蕩培養8 h 左右。參照DNA 模板PCR 擴增體系將DNA模板換成菌液進行PCR 鑒定，2%凝膠電泳檢測，挑選優質（單一清晰，片段大小與目的片段一致的條帶）對應菌液送至測序。

1.2.7 統計分析 CpG 位點的甲基化水平=各CpG 位點的甲基化/ 克隆數；應用QUMA（http://quma.cdb.riken.jp/）在線軟件分析克隆測序結果。

2 結果與分析

2.1 CpG 島在線預測 預測結果為：CpG 島長度>100 bp，GC 含量>50%，期望值>0.6，檢測區域符合實驗要求，預測結果如圖1 所示。

圖1 FSHR 基因啟動子區CpG 島預測結果

2.2 基因組DNA 的提取 利用基因組提取試劑盒按照操作步驟提取基因組DNA，使用1% 的瓊脂糖凝膠電泳以及核酸濃度檢測儀分別檢測其濃度、純度。提取的基因組DNA 濃度均可達到200 ng/μL 以上，OD260/OD280值均在1.8～2.0 之間。因此，提取的血液組織基因組DNA 質量較好，滿足后續的實驗條件，電泳檢測結果如圖2 所示。

圖2 綿羊基因組DNA

2.3 基因組DNA 甲基化BSP 擴增 分別以重亞硫酸鹽修飾處理過的多浪羊與卡拉庫爾羊血液組織基因組DNA 為模板，對FSHR基因啟動子區進行BSP-PCR 擴增，產物用2%瓊脂糖凝膠電泳檢測，擴增片段長度為189 bp，擴增片段大小均與預期擴增片段大小相符（圖3）。

圖3 BSP 處理后FSHR 擴增產物電泳結果

2.4 DNA 甲基化BSP 擴增產物純化回收 對FSHR基因擴增產物進行切膠回收純化，回收純化的DNA 產物片段長度為189 bp，與預期結果一致（圖4）。

圖4 純化后FSHR 擴增產物電泳結果

2.5 DNA 甲基化BSP 克隆 回收純化的DNA 片段分別進行TA 克隆和菌液PCR 檢測，菌液PCR 擴增產物使用2%瓊脂糖凝膠電泳檢測，片段長度都為189 bp（圖5），表明重組克隆載體pGM-T 構建成功。

圖5 重組克隆載體菌液PCR 產物電泳結果

2.6 測序結果及分析 如圖6 所示，非甲基化位點中的非甲基化的“C”均轉化為“T”，而其余堿基完全一致。

圖6 DNAMAN 序列比對結果

運用QUMA 在線網站（http://quma.cdb.riken.jp/）對6 個測序序列進行分析后發現，目的片段中的4 個位點為甲基化位點，如圖7 所示，多浪羊中28、49、71位點均不發生甲基化，61 位點發生甲基化。卡拉庫爾羊中1-3、1-4 羊在28、71 位點不發生甲基化，在49、61 位點發生甲基化。1-5 號卡拉庫爾羊在28 位點不發生甲基化，而在49、61、71 位點均發生甲基化（圖7）。

圖7 甲基化位點QUMA 分析結果

由圖8 可知，多浪羊FSHR基因啟動子區中的目的序列甲基化位點低于卡拉庫爾羊。

圖8 甲基化位點QUMA 分析結果

如圖9 所示，多浪羊組與卡拉庫爾羊組FSHR基因在啟動子區目的片段中的28 位點都不發生甲基化；在49 位點只有卡拉庫爾羊發生甲基化，甲基化率100%，而多浪羊組不發生甲基化；61 位點多浪羊組與卡拉庫爾羊組全部發生甲基化，甲基化率100%；71 位點只有卡拉庫爾羊發生甲基化，甲基化率33.3%，多浪羊組的平均甲基化率為25.0%，卡拉庫爾羊組的平均甲基化率為58.3%。

圖9 甲基化位點QUMA 分析結果

3 討 論

本實驗通過DNA 甲基化檢測的“黃金標準”BSP檢測了多浪羊與卡拉庫爾羊FSHR基因啟動子區甲基化水平，結果顯示在多浪羊和卡拉庫爾羊FSHR基因啟動子區目的片段中發現4 個甲基化位點，且多浪羊甲基化發生率顯著低于卡拉庫爾羊，推測在多浪羊與卡拉庫爾羊中FSHR基因表達與DNA 甲基化水平及甲基化位點具有一定的相關性。

研究發現，在哺乳動物中，DNA 甲基化模式似乎在分化和分化細胞的基因調控方面起著關鍵作用[8]。DNA 甲基化通過圖式和數量的改變對生物信息學進行調節，在基因表達調控等方面發揮重要作用[9]。基因啟動子區及附近區域CpG 島的甲基化是許多基因完成表達與沉默的調控方式，通過對啟動子區CpG 島甲基化狀態檢測來反映出基因表達情況。通過對大鼠FSHR基因的核心啟動子區研究發現，7 個CpG 二核苷酸位點與小鼠FSHR基因啟動子區的4 個CpG 二核苷酸位點，經BSP 測序發現，在表達FSHR的大鼠/小鼠細胞中這些潛在甲基化位點均未被甲基化，而在不表達FSHR的大鼠/ 小鼠的組織細胞中這些位點的CpG 二核苷酸均被甲基化[10]。盡管大鼠/ 小鼠FSHR啟動子中CpG 二核苷酸并不豐富，但啟動子中非特定位點的DNA 甲基化與基因失活相關，DNA 甲基化在大鼠/ 小鼠FSHR基因的調節過程中起著主要作用[11]。蔣曹德等[12]利用甲基敏感擴增多態技術（MSAP）分析了甲基化與生長性狀的關系，在檢測到的1 274 個甲基化位點中，有252 個多態性甲基化位點，其中有81 個位點對1 個或多個生長性狀有顯著影響。本實驗在綿羊血液FSHR基因啟動子區目的片段中找到4 個甲基化潛在位點，其中在多浪羊組中只有在61 位點發生甲基化，而卡拉庫爾羊組中分別在49、61 位點發生甲基化，71 位點部分卡拉庫爾羊發生甲基化，推測在不同位點發生甲基化可能對FSHR基因的表達與繁殖力具有一定的影響。研究發現，源自睪丸腫瘤的小鼠睪丸支持細胞株（MSC-1）具有非活性的FSHR啟動子[13]，并且MSC-1 細胞中轉錄的非活性狀態與FSHR基因核心啟動子的胞嘧啶甲基化相關[14]。張燕麗等[15]對綿羊卵巢與多產性相關的DNA甲基化譜進行全基因組分析，發現一些與激素功能相關的差異甲基化區域相關基因，如FSHR和FST，高繁殖力群體組與低繁殖力群體組相比卵巢組織中的甲基化水平（FSHR和FST上調，LHCGR下調）有顯著差異，這可能會影響這些基因的mRNA 表達。本實驗發現，在高繁殖力的多浪羊與低繁殖力的卡拉庫爾羊中FSHR基因啟動子區的目的片段中49 位點差異極顯著，這可能是影響FSHR基因表達量的因素之一，這一甲基化位點與FSHR基因的表達量是否相關還需進一步研究。

研究表明，基因的啟動子區域和外顯子的低甲基化和高表達呈負相關，即甲基化密度高抑制基因表達，密度低促進基因表達[16]。梁慧慧等[17]對綿羊下丘腦GNAQ基因研究發現適量的葉酸濃度使DNMT1基因表達上調，GNAQ啟動子區CpG 位點甲基化水平上調，引起GNAQ表達下調，進而間接調控了GnRH的分泌。羅榮松等[18]對奶綿羊與蒙古羊全基因組DNA 甲基化研究中發現，奶綿羊和蒙古羊肌肉和尾脂組織具有一致的DNA 甲基化動態；尾脂組織的甲基化水平高于肌肉組織。啟動子區的GC 含量從轉錄起始位點向側翼區域逐漸降低隨后趨于穩定，并與DNA 甲基化水平呈負相關，GC 含量越高甲基化水平越低。王小莉等[19]通過對大鼠不同細胞的研究，得到大鼠胰島和INS-1細胞中INS-1基因高表達，啟動子區域的甲基化程度非常低，然而在大鼠其他組織和WB 細胞中INS-1基因不表達，啟動子區域的甲基化程度很高。本實驗研究結果為在卡拉庫爾羊血液中FSHR基因甲基化水平顯著高于多浪羊，而基因的甲基化一般對基因的表達量有負調控作用，推測多浪羊與卡拉庫爾羊DNA 甲基化水平可能與FSHR基因的表達呈一定負相關。FSHR基因的DNA 甲基化水平及位點差異可能是影響繁殖力的重要因素之一，本研究結果為進一步探究FSHR基因在綿羊卵巢中的表達與DNA 甲基化水平的關系提供依據。

4 結 論

多浪羊FSHR基因啟動子區中目的片段平均甲基化水平為25%，卡拉庫爾羊FSHR基因啟動子區目的片段平均甲基化水平為58.3%，多浪羊與卡拉庫爾羊在目的片段CpG 島上49 號位點甲基化水平差異較大，2 個品種間FSHR基因啟動子區目的片段DNA 甲基化水平差異顯著，推測FSHR基因甲基化水平與綿羊FSHR基因的表達呈負相關。