山羊環狀RNA circ_ZCCHC24的鑒定和表達分析

陶 虎，何 雷，熊 琪，張 年，劉 洋，陳明新*

(1.湖北省農業科學院畜牧獸醫研究所/動物胚胎工程及分子育種湖北省重點實驗室，湖北 武漢 430064；2.華中農業大學動物科技學院，湖北 武漢 430070)

【研究意義】動物組織中存在各種RNA，包括信使RNA、核糖體RNA、非編碼RNA等，它們結構不同，功能各異，但都在生命活動中發揮重要作用。【前人研究進展】環狀RNA(Circular RNA, circRNA)是一類新型的內源性非編碼RNA，具有環狀穩定結構和組織表達特異性等特性。目前山羊卵泡circRNA的研究還處于初始階段，研究其調控山羊卵泡發育的調控作用，可以進一步闡明circRNA的功能和機制，為家畜育種提供理論基礎和新見解。CircRNA是在mRNA前體剪切過程中，外顯子和(或)內含子的5’端與3’端以反向剪切(Back Splicing)形式連接，最終形成共價閉合環狀結構[1-2]。circRNA廣泛存在于真核細胞內，具有物種保守性、組織特異性和穩定性等特性[3]。主要通過3種機制發揮調控作用：①CircRNA作為miRNA海綿調控基因的表達；②CircRNA結合RNA結合蛋白(RBP)形成RNA-蛋白復合物，調控線性親本基因的轉錄[4]；③CircRNA可以編碼蛋白質，發揮生物學功能[5-6]。CircRNA調控卵泡發育的研究較為有限。Capel等人在小鼠精子決定基因SRY中發現了circRNA[7]，并且SRY也被證實可以吸附miR-138，發揮“miRNA海綿”作用[8]。【本研究切入點】繼miRNA和lncRNA后，circRNA已成為非編碼RNA研究領域的新熱點。近年來對于circRNA的研究日益增多，但是對circ_ZCCHC24在山羊卵泡中的作用機制還未見詳細研究。【擬解決的科學問題】本研究以高低繁殖力山羊排卵前卵泡中差異表達的circ_ ZCCHC24為研究對象。對circ_ ZCCHC24是否正確成環進行鑒定，分析它的組織表達特異性，為揭示circ_ ZCCHC24的遺傳調控機理提供理論依據。

1 材料與方法

1.1 樣品采集

山羊組織(腎、心、肌肉、肝、脾、卵泡和肺)樣品采集自湖北省農業科學院種羊場和湖北波爾山羊保種場，山羊選擇18月齡的成年母羊；屠宰后10 min內采集腎、心、肌肉、肝、脾、卵泡和肺組織各3份，經PBS清洗后放入去RNA酶的樣品采集管內，置于液氮帶回實驗室進行總RNA的提取，剩余樣本放置在-80℃冰箱保存。

1.2 總RNA的提取及反轉錄

利用天根RNA simple總RNA提取試劑盒來山羊不同組織樣品中總RNA的提取，操作按照試劑盒說明書進行，使用NanoDrop 2000型DNA/RNA濃度測定儀，測定DNA濃度和OD值，利用1%瓊脂糖凝膠電泳檢測RNA的完整性。采用Takara的PrimeScript?RT reagent Kit with gDNA Eraser反轉錄試劑盒獲得cDNA模板。操作按照試劑盒說明書進行，反轉錄后的cDNA放置于-80 ℃冰箱保存備用。

1.3 生物信息學分析

采用轉錄組測序技術分析了波爾山羊和麻城黑山羊排卵前卵泡組織中的circRNA，利用生物信息學手段分析轉錄組測序數據，獲得山羊circ_ZCCHC24的全長序列、接頭序列和形成結構。

1.4 引物設計

以circ_ZCCHC24的全長序列為模板，利用Primer 5軟件設計鑒定山羊circ_ZCCHC24的引物，預期擴增片斷長度為78 bp，相關引物信息詳見表1，引物由北京奧科鼎盛生物科技有限公司合成。

1.5 PCR擴增

PCR反應總體系50 μl，其中山羊cDNA模板2 μl，I-5TM2×High Fidelity Master Mix試劑 25 μl，上、下游引物(0.4 μM)各2 μl，去離子水加至50 μl。PCR擴增程序為98 ℃預變性2 min，98 ℃變性10 s，60 ℃退火10 s，72 ℃延伸10 s，40個循環；72 ℃延伸5 min，16 ℃保存。PCR儀為Mini CyclerTM基因擴增儀，PCR擴增產物用1.5 %瓊脂糖凝膠電泳檢測。利用DNA快速純化回收試劑盒回收目的片段，具體操作按產品說明書進行。

1.6 目的片段測序驗證

將PCR產物膠回收產物連接pMD18-T 線性載體質粒，連接反應體系：8 μl純化后的PCR產物，10×Ligation Buffer 1 μl，pMD18-T質粒1 μl，反應體系于16 ℃金屬浴條件下連接1 h。隨后加入50 mL DH5α感受態細胞，均勻涂布于LB固體培養基上于37℃培養箱中過夜，陽性菌液送北京奧科生物科技公司測序。

1.7 實時熒光定量PCR(qRT-PCR)

qRT-PCR反應總體系20 μl，其中山羊不同組織cDNA模板各1 μl，SYBR?Green Supermix 10 μl，上、下游引物(0.4 μM)各0.25 μl，去離子水加至20 μl。PCR擴增程序為94 ℃預變性4 min，94 ℃變性15 s，59 ℃退火15 s，72 ℃延伸15 s，40個循環；72℃延伸45 s。儀器采用Mini CyclerTM基因擴增儀，利用2-ΔΔCT方法進行分析。

1.8 熒光原位雜交(FISH)

將麻城黑山羊卵泡組織放入固定液中固定6 h，經梯度酒精脫水后浸蠟包埋，石蠟經切片機切片，62 ℃烤片2 h；切片于修復液中煮沸15 min，滴加蛋白酶K(20 μg/mL)37 ℃消化30 min；滴加預雜交液，37 ℃孵育1 h；滴加含探針circ_ZCCHC24-Probe的雜交液(8 ng/μl)，37 ℃雜交過夜后進行洗滌，探針序列詳見表1；切片滴加DAPI(4',6-diamidino-2-phenylindole)染液，避光孵育8 min，沖洗后滴加抗熒光淬滅劑封片；處理后的切片利用尼康ECLIPSE CI正置熒光顯微鏡下觀察并采集圖像。

表1 引物與探針信息

2 結果與分析

2.1 山羊circRNA circ_ZCCHC24的生物信息學分析

前期利用轉錄組測序技術分析了繁殖力顯著差異的波爾山羊和麻城黑山羊排卵前卵泡組織中的circRNA，篩選出circ_ZCCHC24在2個山羊品種中的表達存在顯著差異。根據轉錄組測序得出circ_ZCCHC24堿基全長為2045 bp(圖1A)，接頭序列信息如圖1B所示。將測序所得circ_ZCCHC24的序列信息在NCBI數據庫中進行檢索，circ_ZCCHC24的全長序列是由山羊ZCCHC24基因(GenBank登錄號：NC_030835.1)的部分第1內含子和全部第2外顯子堿基序列環化而成(圖1C)。表明circ_ZCCHC24是由外顯子和內含子共同組成的circRNA。

圖1 山羊circ_ZCCHC24序列信息(A)、接頭序列信息(B)及結構示意圖(C)Fig.1 Information of goat circ_ZCCHC24 sequence (A), junction sequence (B) and schematic of goat circRNA circ_ZCCHC24 (C)

圖2 利用分散引物檢測circRNA的示意圖(A)，分散引物擴增的circ_ZCCHC24 PCR產物電泳圖(B)及采用分散引物PCR擴增circ_ZCCHC24接頭序列的測序結果(C)Fig.2 Schematic representation of the detection circRNAs using divergent primer (A), electrophoretogram of the circ_ZCCHC24 PCR products obtained with divergent primers (B), sanger sequencing of circ_ZCCHC24 PCR products resulting from amplification using divergent primers confirms head-to-tail splicing (C)

圖3 山羊circ_ZCCHC24在不同山羊品種中的表達水平(A)及在波爾母羊不同組織中的表達水平(B)Fig.3 Expression level of circ_ZCCHC24 in two different goat breeds (A) and seven tissues of Boer ewes (B)

2.2 山羊circRNA circ_ZCCHC24的鑒定

為了確定山羊circRNA circ_ZCCHC24是否正確成環，設計了分散引物(Divergent Primer)用于擴增circ_ZCCHC24的接頭序列，從而進行成環鑒定(圖2A)。利用分散引物的PCR擴增獲得了circ_ZCCHC24接頭序列產物，經凝膠電泳結果顯示circ_ZCCHC24的接頭序列順利擴增，擴增產物長度為78 bp(圖2B)，再利用Sanger測序進一步證實了接頭序列中存在環化位點(圖2C)。上述結果表明circ_ZCCHC24在山羊卵泡組織中能夠正確成環。

2.3 山羊circRNA circ_ZCCHC24的表達模式分析

通過qRT-PCR分析了circ_ZCCHC24在2個山羊品種中的表達水平，發現circ_ZCCHC24在麻城黑山羊中的表達量遠高于波爾山羊(P≤0.01，圖3A)。此外，研究了circ_ZCCHC24在波爾母羊腎、心、肌肉、肝、脾和卵泡等不同組織中的相對表達水平，結果發現，circ_ZCCHC24在上述6種組織中均有表達，其中在腎中表達量最高，卵泡中表達量最低(圖3B)。表明circ_ZCCHC24的表達水平具有組織特異性。

2.4 山羊circRNA circ_ZCCHC24的細胞定位

為了檢測山羊circ_ZCCHC24在山羊卵泡顆粒細胞中的細胞定位情況，針對接頭序列設計了CY3標記的探針，利用熒光原位雜交(FISH)技術進行了細胞定位分析，結果顯示，circ_ZCCHC24定位于山羊卵泡顆粒細胞的胞質中(圖4)。表明circ_ZCCHC24可能發揮作用的形式是通過吸附miRNA來完成的。

3 討 論

哺乳動物卵泡的產生是一種由原始卵泡發育成排卵前卵泡的連續復雜過程[9]。其中，卵泡顆粒細胞在卵泡發育過程中發揮重要作用，為卵母細胞提供必需的生長因子和特定蛋白質，卵泡顆粒細胞的增殖可誘導卵泡的生長和卵母細胞的成熟[10]。卵泡發生的復雜性表現在卵母細胞的發育需要眾多基因組成的調控網絡來發揮作用[11]。然而，動物卵泡發育和顆粒細胞增殖的決定性調控機制仍需要進一步研究。

長久以來，非編碼RNA被認為是一種轉錄垃圾，不發揮任何作用。然而越來越多的研究表明，包括miRNAs、circRNAs、piRNA和lncRNA等在各種生物學過程中發揮重要的調控功能[1, 12-14]。尤其circRNA更是成為現在基礎研究領域的熱點，因此，我們以高低繁殖力山羊排卵前卵泡中差異表達的circRNA circ_ ZCCHC24為主要研究對象，對其進行成環鑒定，分析其結構和表達分布情況。

通過生物信息學分析發現circ_ ZCCHC24由山羊ZCCHC24基因的部分第1內含子和第2外顯子環化而成。circRNA的成環來源類型分類，包括外顯子來源的circ RNA、內含子來源的circRNA和外顯子和內含子共同組成的circ RNA[15]。本研究證實的由部分內含子和完整外顯子組成的circ_ ZCCHC24，則可能是一種新的circRNA來源類型。利用sanger測序技術明確了circ_ ZCCHC24的接頭序列，從而證明circ_ ZCCHC24是真正的circRNA。利用測序技術明確接頭序列是檢驗circRNA是否正確成環的金標準[1]。分析了circ_ ZCCHC24的表達模式，發現其具有顯著的組織特異性，其中在波爾山羊的腎臟中表達最高，在脾和卵泡中表達最低，由于circ_ ZCCHC24本身實在波爾山羊卵泡中低表達，因此上述表達模式結果符合預期。

CircRNA在生物體內發揮重要的調控作用，CircRNA在細胞核內調控來源基因(Host gene)的表達[16]，在細胞質內則發揮內源性競爭RNA(competing endogenous RNAs, ceRNA)的作用[17]。為了檢測circ_ZCCHC24在山羊卵泡顆粒細胞中的細胞定位情況，利用FISH技術對circ_ZCCHC24進行了顆粒細胞定位分析，結果顯示circ_ZCCHC24定位于山羊卵泡顆粒細胞的胞質中(圖4)，表明circ_ZCCHC24可能發揮“miRNA海綿”的功能。后續將開展與circ_ZCCHC24具有結合關系的miRNA的篩選和驗證工作。

4 結 論

本研究首次對山羊卵泡差異表達的circRNA circ_ZCCHC24進行了初步研究，對其進行了成環驗證、表達分析和細胞定位。為后期進一步研究circ_ZCCHC24調控母羊卵泡發育的分子機理，挖掘出高產母羊的特有調控通路奠定了基礎，也為高繁殖力山羊品種的培育提供新靶點。