• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    豬耳成纖維細胞轉(zhuǎn)錄組異質(zhì)性及對核移植胚胎發(fā)育的潛在影響

    2020-09-24 01:17:08周俊趙成成吳霄石俊松周榮吳珍芳李紫聰
    遺傳 2020年9期
    關(guān)鍵詞:單細胞差異基因體細胞

    周俊,趙成成,吳霄,石俊松,周榮,吳珍芳,李紫聰

    研究報告

    豬耳成纖維細胞轉(zhuǎn)錄組異質(zhì)性及對核移植胚胎發(fā)育的潛在影響

    周俊1,趙成成1,吳霄1,石俊松2,周榮2,吳珍芳1,李紫聰1

    1. 華南農(nóng)業(yè)大學動物科學院,國家生豬種業(yè)工程技術(shù)研究中心,廣州 510642 2. 廣東溫氏種豬科技有限公司,新興 527400

    同一來源的供體細胞之間存在異質(zhì)性。許多研究已經(jīng)表明體細胞核移植(somatic cell nuclear transfer, SCNT)效率與供體細胞有關(guān)。然而,鮮有在單細胞水平分析供體細胞異質(zhì)性對核移植效率的潛在影響。本研究利用單細胞轉(zhuǎn)錄組測序技術(shù)對同一來源且隨機挑選的52個豬耳組織成纖維細胞進行測序分析。結(jié)果表明有48個單細胞的基因表達模式相似,4個單細胞(編號為D11_1、D12_1、DW61_2和DW99_2)的基因表達模式與其他單細胞存在較大的差異,并且不存在基因表達模式完全相同的兩個單細胞。以基因表達模式相似的48個單細胞作為對照,進一步分析了單細胞D11_1、D12_1、DW61_2和DW99_2的差異基因表達模式:首先利用R語言篩選4個單細胞的差異表達基因,并對前50差異表達基因進行匯總;然后對差異表達基因進行GO富集分析和KEGG通路分析。富集分析發(fā)現(xiàn)差異表達基因的主要分子功能包括能量代謝、蛋白質(zhì)代謝和細胞對刺激的反應(yīng)等;主要通路包括KEGG中富集的與細胞周期、細胞代謝、DNA復制相關(guān)的通路。根據(jù)以上研究結(jié)果并結(jié)合SCNT研究進展討論了4個單細胞的差異基因表達模式對核移植胚胎發(fā)育效率的潛在影響。本研究揭示了豬耳組織成纖維細胞的轉(zhuǎn)錄組異質(zhì)性,并提供了分析精英供體細胞的一種有效方法,為提高克隆效率帶來新的思路。

    豬;體細胞核移植;供體細胞;異質(zhì)性;單細胞測序

    體細胞核移植技術(shù)(somatic cell nuclear transfer, SCNT)又稱為體細胞克隆技術(shù)。1997年,英國胚胎學家Wilmut利用SCNT技術(shù)成功克隆出第一只哺乳動物—克隆羊“多莉”[1],隨后,通過該技術(shù)先后誕生了牛()[2]、豬()[3]、猴()[4]等20余種哺乳動物。SCNT技術(shù)在保護瀕危物種[5]、優(yōu)秀種繁擴繁[6]以及再生醫(yī)學[7]等領(lǐng)域表現(xiàn)出巨大的應(yīng)用價值,是生命科學研究的重要成果。但是,SCNT技術(shù)同樣面臨很多問題:以豬為例,體外培養(yǎng)條件下,囊胚率僅為20%左右,在不同品種或者卵母細胞質(zhì)量差的情況下可能更低。即便胚胎成功附植,相較于體外受精(fertilization, IVF)胚胎,SCNT胚胎發(fā)育能力差、在受孕母豬體內(nèi)發(fā)育異常的現(xiàn)象也十分普遍[8]。除此之外,大多數(shù)實驗條件下,克隆豬的出生率只有1%左右,遠低于人工授精(artificial insemination, AI)豬的出生效率(80%)[9],嚴重限制了SCNT技術(shù)的應(yīng)用和普及。SCNT技術(shù)基礎(chǔ)是將供體細胞移入去核的成熟卵母細胞中,最終發(fā)育成和供體細胞基因型相同的后代。因此供體細胞是影響體細胞核移植效率的關(guān)鍵因素。已有研究表明供體細胞會對SCNT效率產(chǎn)生影響:供體細胞的細胞周期[10]、傳代數(shù)[11]以及性別[12]的不同都會導致克隆胚胎發(fā)育效率出現(xiàn)差異。因此,選擇合適的供體細胞對提高SCNT胚胎發(fā)育效率十分重要。

    供體細胞核在成熟的去核卵母細胞中異常的表觀遺傳重編程被認為是阻礙SCNT發(fā)育的主要原因[13,14]。而更容易被正確重編程的核供體細胞很可能是具有更大發(fā)育潛力的精英供體細胞。Zhai等[15]分別以豬骨髓間充質(zhì)干細胞和豬胎兒成纖維細胞作為供體細胞進行核移植實驗,發(fā)現(xiàn)發(fā)育效率較高的骨髓間充質(zhì)干細胞含有更多的有利于重編程的表觀遺傳標記和較少的抑制重編程的表觀遺傳標記,暗示容易被正確重編程的供體細胞在SCNT過程中會有更好的發(fā)育潛能。同時,Yamanaka等[16]在鼠源誘導多能干細胞研究中就已經(jīng)提出精英供體細胞的概念:同一來源的細胞中會存在一些更容易被正確重編程、發(fā)育潛力更大的精英供體細胞;同樣,在克隆小鼠()的研究中也發(fā)現(xiàn),與其他品系的小鼠相比,129小鼠的基因組狀態(tài)更不穩(wěn)定,更容易被激活或抑制。用野生基因型和129小鼠基因型雜交得到的重組細胞做供體,可顯著提高克隆動物的出生率[17],說明精英供體細胞具有獨特的基因表達模式,只是其分子特征尚未探索清楚。另有研究顯示,從供體細胞遺傳而來的一些轉(zhuǎn)錄記憶可以導致克隆胚胎的發(fā)育缺陷[18],表明供體細胞在一定程度上決定了重構(gòu)胚胎的發(fā)育命運。供體細胞的異質(zhì)性是指供體細胞在基因組或表型水平上具有的不同特征。同時,因為細胞之間異質(zhì)性的存在,不同的供體細胞發(fā)育潛力是有差異的。即便同一來源的供體細胞也存在更有利于胚胎發(fā)育的精英供體細胞。但是在實際研究中,同一來源的供體細胞之間表型特征差異不明顯,分子層面的特征信息丟失嚴重,缺乏對常見供體細胞異質(zhì)性的深入研究。

    近年來,隨著測序技術(shù)的發(fā)展,特別是低輸入測序技術(shù)為體細胞核移植的研究提供了更多可能。而單細胞測序技術(shù)的誕生及發(fā)展[19~21]為人們進一步探究供體細胞之間的分子事件提供了便利。應(yīng)用單細胞測序技術(shù),可以以更精準的分辨率揭示供體細胞的異質(zhì)性?;诖?,本研究對來源相同的52個豬耳組織成纖維細胞進行單細胞轉(zhuǎn)錄組測序分析,發(fā)現(xiàn)了48個基因表達模式相似的“普通”細胞,4個基因表達模式互不相似且都與“普通”細胞的基因表達模式存在顯著差異的“另類”細胞,揭示了豬耳組織成纖維細胞的轉(zhuǎn)錄組異質(zhì)性,并結(jié)合已有研究結(jié)果討論了供體細胞轉(zhuǎn)錄組異質(zhì)性對核移植效率的潛在影響,為后續(xù)通過深入研究尋找精英供體細胞來提高豬體細胞核移植效率提供基礎(chǔ)。

    1 材料與方法

    1.1 材料

    杜洛克公豬(373日齡)耳組織樣品由廣東溫氏食品集團華農(nóng)溫氏股份有限公司提供;總RNA提取試劑盒SMART-SeqTMv4 UltraTMLow Input RNA Kit for Sequencing購于北京諾禾致源生物信息科技有限公司;胎牛血清(fetal bovine serum, FBS)、氨基酸葡萄糖培養(yǎng)基(dulbecco's modified eagle medium, DMEM)、0.25%胰蛋白酶和乙二胺四乙酸均購自美國Gibco公司。單細胞轉(zhuǎn)錄組測序由北京諾禾致源生物信息科技有限公司在illumina平臺完成。圖片處理使用Photoshop 7.0/ACDSee 9.0軟件;spliced reads比對使用HISAT軟件[22];聚類分析、基因差異表達分析、GO富集分析和主成分分析(principal component analysis, PCA)均采用R語言,其中聚類分析使用軟件包pheatmap,基因差異表達分析使用軟件包DEGSeq (1.12.0)[23],富集分析采用軟件包GOseq[24];KEGG富集使用軟件KOBAS (2.0)。

    1.2 供體細胞單細胞培養(yǎng)與分離

    將成年優(yōu)良杜洛克種公豬的耳組織樣品剪碎,用PBS洗滌兩次之后,在100 mm的培養(yǎng)皿上用手術(shù)刀和剪子切碎,先用DMEM重懸,再用胰蛋白酶和乙二胺四乙酸消化1~2 h。將胰蛋白酶消化的細胞洗滌一次后,以300的離心率離心10 min,并將其接種于100 mm的細胞培養(yǎng)皿中,放入15% FBS和10 mg/L的青霉素–鏈霉素溶液的DMEM。在39℃、5%CO2飽和濕度的恒溫培養(yǎng)箱中培養(yǎng)6~8 d。然后移除未附著的組織塊,再將附著的細胞培養(yǎng)直至匯合,期間每隔3~7 d更換DMEM。最后利用顯微操作法,在體視顯微鏡下吸取單個細胞,放于盛有裂解液的去DNase-RNase的離心管中,用于單細胞轉(zhuǎn)錄組測序。

    1.3 單細胞轉(zhuǎn)錄組數(shù)據(jù)獲得

    1.3.1 總RNA提取

    每個細胞樣品保存在6 μL SMART-SeqTMv4 kit裂解液(北京諾禾致源生物信息科技有限公司)中。細胞樣品經(jīng)過體積測量后,使用SMART-SeqTMv4 UltraTMLow Input RNA Kit for Sequencing試劑盒(美國Clontech公司)進行細胞裂解,提取總RNA,并保存在RNase-Free水中。

    1.3.2 單細胞cDNA文庫構(gòu)建及測序

    對提取的總RNA直接進行First-stand cDNA的合成,然后對First-stand cDNA進行全長LD-PCR的擴增,利用AMPure XP beads純化擴增后的雙鏈cDNA (double-standed DNA, ds cDNA),使用Qubit進行ds cDNA定量檢測;使用Covaris系統(tǒng)對ds cDNA進行超聲打斷,打斷后的雙鏈短片段進行末端修復、加A尾并連接測序接頭,然后用AMPure XP beads純化并選擇片段大小在200 bp左右的文庫;最后進行PCR富集得到最終的cDNA文庫。使用Qubit2.0對文庫進行初步定量,稀釋文庫至1 ng/μL,使用Agilent 2100對文庫的插入片段長度進行檢測。插入片段符合預(yù)期后,使用qPCR方法對文庫的有效濃度進行準確定量,以保證文庫質(zhì)量。庫檢合格后,將不同文庫按照有效濃度及目標下機數(shù)據(jù)量的需求合并后進行HiSeq測序。本研究最后擴增成功且完成單細胞測序的樣品為52個,細胞樣品編號見表1。

    1.4 單細胞轉(zhuǎn)錄組數(shù)據(jù)分析

    1.4.1 測序數(shù)據(jù)質(zhì)量控制

    fastq格式的原始數(shù)據(jù)先通過內(nèi)部perl腳本進行處理。在此步驟中,刪除包含適配器的reads、包含ploy-N的低質(zhì)量的原始數(shù)據(jù)來獲得干凈的數(shù)據(jù)(clean reads),同時對Q20、Q30和GC內(nèi)容進行計算,所有的下游分析都是基于高質(zhì)量的清潔數(shù)據(jù)。

    表1 擴增成功并測序的單細胞樣品

    1.4.2 差異基因表達分析

    參考基因組和基因模型注釋文件直接從基因組網(wǎng)(http://www.ensembl.org/)下載。選取HISAT軟件將過濾后的測序序列進行基因組定位分析。HISAT能夠有效的比對到RNA-Seq測序數(shù)據(jù)中的spliced reads,是目前比對率最高且最準確的比對軟件。先使用Hisat2 v2.0.4作為映射工具,它可以基于基因模型注釋文件生成一個拼接連接的數(shù)據(jù)庫,因此比其他非拼接映射工具具有更好的映射結(jié)果。然后使用軟件HTSeq v0.9.1計算映射到每個基因的讀取數(shù)字。然后根據(jù)基因的長度計算出每個基因的FPKM,并讀取到該基因的計數(shù)。在進行差異基因表達分析之前,通過edgeR程序包[25]對每一個序列庫標準化。利用R語言中的DEGSeq (1.20.0)軟件包進行微分表達式分析并繪制4個“另類”細胞差異基因火山圖、聚類軟件包pheatmap繪制4個“另類”細胞差異基因聚類圖。值用和法進行調(diào)整,修正的值為0.005和log2(fold change)為1,為顯著差異表達的閾值。

    1.4.3 GO富集分析和KEGG通路分析

    通過R語言中的GOseq軟件包對篩選得到的差異基因進行GO富集,展示差異基因在Gene Ontology (http://www.geneontology.org/)中的分布狀況,闡明本研究中不同細胞基因功能上的差異。

    利用KOBAS(2.0)軟件對差異表達基因進行KEGG (http://www.genome.jp/kegg/)通路統(tǒng)計,得到差異表達基因參與的最主要生化代謝途徑和信號轉(zhuǎn)導途徑。其原理是應(yīng)用超幾何檢驗,找出與整個基因組背景相比,在差異表達基因中顯著性富集的通路。

    2 結(jié)果與分析

    2.1 單細胞轉(zhuǎn)錄組測序揭示豬耳組織成纖維細胞間的異質(zhì)性

    按照PCA的前3個主要影響因素兩兩組合,對成功擴增且最終完成測序的52個單細胞進行主成分分析(圖1),其中Dim1、Dim2、Dim3分別代表影響單細胞樣本在圖中相對位置的3個主要因素,相對位置較近的細胞表示基因表達模式相似。結(jié)果表明,D11_1、D12_1、DW61_2和DW99_2在圖1中的位置與其他單細胞相比較“離群”,說明這4個細胞的基因表達模式與其他單細胞存在較大差異。同時將52個耳組織成纖維細胞的基因表達總體情況繪制成熱圖(圖2),結(jié)果表明,在同一來源的豬耳成纖維細胞中,同一基因在不同細胞間的表達情況并不相同,并不存在基因表達模式完全一致的兩個細胞,不同供體細胞之間存在轉(zhuǎn)錄組異質(zhì)性。

    2.2 D11_1、D12_1、DW61_2和DW99_2基因表達模式分析

    將基因表達模式差異不明顯的“普通”細胞作為對照組,使用R語言中的軟件包DEGSeq (1.12.0)對D11_1、D12_1、DW61_2和DW99_2進行差異基因表達分析。D11_1與對照組比較,有1860個差異表達顯著的基因,其中上調(diào)表達的基因有458個,下調(diào)表達的基因有1402個;D12_1與對照組比較,有376個差異表達顯著的基因,其中上調(diào)表達的基因有109個,下調(diào)表達的基因有267個;DW61_2與對照組比較,有330個差異表達顯著的基因,其中上調(diào)表達的基因有277個,下調(diào)表達的基因有53個;DW99_2與對照組比較,有2225個差異表達顯著的基因,其中上調(diào)表達的基因有316個,下調(diào)表達的基因有1909個。使用R語言中軟件包DEGSeq (1.12.0)將4個“另類”細胞的差異基因的整體分布情況可視化(圖3),使用R語言中的聚類軟件包pheatmap將4個“另類”細胞的差異基因整體表達情況可視化(圖4)。同時,將用軟件包DEGSeq (1.12.0)篩選出的D11_1、D12_1、DW61_1和DW99_2前50個差異最顯著的基因進行匯總:D11_1前50個差異最顯著的基因功能主要集中于細胞增殖分化過程中一些有機物質(zhì)的合成與能量代謝過程(表2);D12_1前50個差異最顯著的基因功能主要集中于蛋白質(zhì)及葡萄糖轉(zhuǎn)運等過程(表3);DW61_2前50個差異最顯著的基因功能比較多樣,但有個別基因涉及轉(zhuǎn)移酶活性(表4);DW99_2前50個差異最顯著的基因功能主要集中于蛋白質(zhì)編碼、核酸修復與能量代謝方面(表5)。

    2.3 D11_1、D12_1、DW61_2和DW99_2差異基因GO富集分析和KEGG通路分析

    對D11_1、D12_1、DW61_2和DW99_2測序得到的差異基因進行GO富集分析,這些差異基因在生物過程、細胞組分和分子功能方面的分布情況見圖5。D11_1上調(diào)的差異表達基因主要集中于與細胞代謝有關(guān)的線粒體等細胞器和細胞器膜,下調(diào)的差異表達基因主要涉及與細胞代謝有關(guān)的蛋白代謝和物質(zhì)運輸生物過程。D12_1上調(diào)的差異表達基因無顯著富集,而下調(diào)的差異表達基因主要富集于細胞對刺激的反應(yīng)與蛋白代謝的細胞過程與應(yīng)對刺激信號傳導的膜的變化。DW61_2上調(diào)和下調(diào)的差異表達基因均無顯著富集。DW99_2上調(diào)的差異表達基因主要富集于細胞有絲分裂有關(guān)的DNA復制、染色體分離、核酸代謝等生物過程和與之相關(guān)的一些胞內(nèi)有機物質(zhì)的變化,而下調(diào)的差異表達基因則主要涉及蛋白代謝和蛋白修飾等一些高分子修飾過程。結(jié)果顯示D11_1和DW99_2的差異基因功能富集趨勢更為明顯,基因表達模式也更為“另類”。

    圖1 52個單細胞轉(zhuǎn)錄組主成分分析

    每個點代表一個單細胞樣品;橫、縱坐標軸的刻度是相對距離,無實際意義;Dim1、Dim2、Dim3后的百分比代表橫、縱的差異可以解釋全面分析結(jié)果的百分比。紅色圈出部分分別為單細胞D11_1、DW61_2、D12_1和DW99_2。

    圖2 52個單細胞基因表達量熱圖

    紅色代表基因高表達,藍色代表基因低表達。

    圖3 D11_1、D12_1、DW61_2和DW99_2差異基因火山圖

    A~D分別為D11_1、D12_1、DW61_2和DW99_2的差異基因火山圖;有顯著性差異表達的基因用紅色點(上調(diào))和綠色點(下調(diào))表示,無顯著性差異表達的基因用藍色點表示;橫坐標代表基因在不同樣本中表達倍數(shù)變化;縱坐標代表基因表達量變化差異的統(tǒng)計學顯著性;篩選標準padj < 0.05。

    圖4 D11_1、D12_1、DW61_2和DW99_2差異基因聚類熱圖

    整體FPKM層次聚類圖,將log10(FPKM+1)值進行歸一化轉(zhuǎn)換(scale number)并進行聚類,紅色表示高表達基因,藍色表示低表達基因。顏色從紅到藍,表示log10(FPKM+1)從大到小。

    利用KEGG注釋系統(tǒng)對D11_1、D12_1、DW61_2和DW99_2進行差異基因代謝通路富集分析。結(jié)果顯示(圖6),D11_1上調(diào)的基因主要作用于細胞代謝,而下調(diào)的基因主要作用于細胞凋亡;D12_1與DW61_2 KEGG均無顯著通路富集;DW99_2 KEGG通路中上調(diào)的差異基因最顯著的富集于調(diào)節(jié)DNA復制、細胞周期通路。

    3 討論

    單細胞轉(zhuǎn)錄組測序技術(shù)能夠揭示單個細胞的基因表達動態(tài),反映細胞間的異質(zhì)性。在豬耳成纖維細胞異質(zhì)性研究過程中,將單個細胞從耳組織塊中分離出來是單細胞轉(zhuǎn)錄組測序的第一個關(guān)鍵步驟。為準確獲得成纖維細胞,排除其他類型細胞對細胞異質(zhì)性結(jié)果的干擾,本研究最終采用顯微操作法進行單細胞分離。該方法適用于樣本數(shù)量少的單細胞樣品制備,能夠在顯微鏡下觀察成纖維細胞的形態(tài),精準控制單個細胞的吸入與排出[26],確保52個單細胞樣品均為成纖維細胞。

    表2 D11_1前50個差異最顯著的基因信息

    表3 D12_1前50個差異最顯著的基因概況

    表4 DW61_2前50個差異最顯著的基因概況

    表5 DW99_2前50個差異最顯著的基因概況

    在生物界異質(zhì)性通常被解釋為3個層面:(1)不同物種或生物體之間具有異質(zhì)性;(2)同一生物體的不同器官或組織具有異質(zhì)性;(3)同一器官或組織中的不同細胞具有異質(zhì)性。事實上,早在1957年Novick等[27]就提出了細胞異質(zhì)性的概念,最初是根據(jù)形態(tài)和功能上的差異將細胞群分化成不同的細胞亞群。單細胞測序技術(shù)的出現(xiàn)為細胞異質(zhì)性的研究提供了全新的見解:每一個細胞都是獨特的個體,被獨一無二的DNA、RNA以及蛋白質(zhì)所編碼[28]。在醫(yī)學領(lǐng)域,已有研究利用單細胞測序技術(shù)構(gòu)建人類器官、組織的基因表達圖譜[29]。這為SCNT研究提供一個可能:應(yīng)用單細胞測序技術(shù)描述供體細胞核和卵母細胞質(zhì)相互作用的動態(tài)過程,以更高的分辨率解析阻礙胚胎發(fā)育的分子原因。在體細胞核移植研究中,由于供體細胞異質(zhì)性的存在,不同的供體核和卵母細胞質(zhì)相互作用表現(xiàn)出不同應(yīng)答模式,包括不同的染色質(zhì)重構(gòu)和表觀遺傳修飾重編程[15]。同樣,在誘導多能干細胞的研究中,來自不同組織的供體細胞在重編程過程中表現(xiàn)出不同的敏感性[30,31],由不同供體細胞獲得的多能干細胞擁有不同的轉(zhuǎn)錄模式[32]。先前大多數(shù)研究討論了不同類型的供體細胞對克隆胚胎發(fā)育效率的影響[33,34],而尚不清楚來源相同的供體細胞異質(zhì)性對核移植效率的影響。在實際研究中,往往把來源相同的供體細胞看成同質(zhì)的,導致細胞內(nèi)關(guān)鍵分子事件大量丟失。

    圖5 D11_1、D12_1、DW61_2和DW99_2差異表達基因的GO富集分析圖

    A~D分別為D11_1、D12_1、DW61_2和DW99_2與對照組差異表達基因的GO富集柱狀圖。縱坐標為富集的GO term,橫坐標為該term中差異基因個數(shù)。不同顏色用來區(qū)分生物過程、細胞組分和分子功能,帶“*”為顯著富集的GO term,對富集最顯著的30個GO term在圖中展示,如果不足30條,則全部展示。

    圖6 D11_1、D12_1、DW61_2和DW99_2差異表達基因的KEGG富集散點圖

    A~D分別代表D11_1、D12_1、DW61_2和DW99_2與對照組的差異表達基因KEGG富集散點圖??v軸表示pathway名稱,橫軸表示Rich factor,點的大小表示此pathway中差異表達基因個數(shù)多少,而點的顏色對應(yīng)于不同的value范圍。

    本研究將精度對準細胞水平,用單細胞轉(zhuǎn)錄組數(shù)據(jù)反映不同供體細胞的基因表達模式:52個成功測序的豬耳成纖維細胞轉(zhuǎn)錄組熱圖反映同一組織供體細胞間的基因表達模式差異;通過對上述52個細胞進行主成分分析,發(fā)現(xiàn) 4個基因表達模式明顯不同的“另類”細胞,證實了供體細胞間異質(zhì)性的存在。同時,推測“另類”細胞的基因表達模式和精英供體細胞存在聯(lián)系。D11_1與“普通”細胞相比,其上調(diào)基因GO富集于與細胞代謝相關(guān)的細胞器、細胞膜,尤其是線粒體,主要功能為促進細胞代謝。其下調(diào)基因在KEGG分析中主要集中在細胞凋亡通路。在正常受精胚胎中,精子來源的mtDNA在受精不久后被全部破壞,其mtRNA全部來源于卵母細胞[35]。但是在SCNT胚胎中卻存在供體來源的mtDNA[36]。因此,這些供體來源的線粒體DNA可能是影響克隆效率的關(guān)鍵線索。同時,D11_1差異基因表達結(jié)果也暗示了通過研究線粒體相關(guān)基因?qū)ふ揖⒐w細胞的可能。在D11_1表達上調(diào)的差異基因中,基因[37,38]可以調(diào)節(jié)精氨酸合成,而L-精氨酸又是細胞信號傳導、代謝功能分子(NO、多胺和肌酸)的主要合成前體[39],對胚胎發(fā)育十分重要。基因的表達能加快蛋白質(zhì)代謝和高爾基體轉(zhuǎn)運,增加胚胎發(fā)育能力[40,41]。在D11_1表達下調(diào)的差異基因中,基因家族蛋白負向調(diào)節(jié)鞘脂代謝[42],其下調(diào)表達會促進細胞生長分化;基因的表達可促進細胞凋亡[43,44],其下調(diào)表達起到抑制細胞凋亡的作用;是一種介導表觀遺傳抑制修飾的多蛋白復合物[45],可以推動細胞重編程障礙H3K9me3[46]的表觀抑制修飾,其下調(diào)表達有利于克隆胚胎正確重編程。因此,D11_1的基因表達模式可作為精英供體細胞的潛在參考。

    DW61_2和D12_1與“普通”細胞比較,GO富集與KEGG富集均無顯著趨勢,其差異表達基因的功能多為抑制蛋白活性和糖代謝、阻遏細胞生長。如D12_1中基因[47]的顯著下調(diào)會抑制結(jié)合蛋白mRNA的轉(zhuǎn)運,進而抑制細胞分裂;基因[48]的下調(diào)會影響高爾基體內(nèi)部囊泡介導物質(zhì)運輸,抑制高爾基體和微管間的相互作用,抑制細胞有絲分裂;在DW61_2中表達下調(diào)的基因中,基因[49]編碼二硫化物還原酶,其下調(diào)會抑制細胞生長,降低細胞活力;基因[50,51]在胚胎發(fā)生和傷口愈合過程中高度表達。該基因產(chǎn)物是一種基質(zhì)細胞蛋白,能促進基質(zhì)組裝,并能刺激內(nèi)皮細胞的增殖和遷移,以及血管生成活性,其表達下調(diào)不利于胚胎發(fā)生和后續(xù)克隆動物的生長發(fā)育。因此,DW61_2、DW12_1的基因表達模式呈現(xiàn)抑制細胞活力,促進細胞凋亡的趨勢,同時根據(jù)已有研究證據(jù)無法將這兩個“另類”細胞與精英供體細胞相關(guān)聯(lián),其基因表達模式并不具備代表性。

    DW99_2與對照組相比,GO富集上調(diào)的差異表達基因主要和DNA復制、染色體分離、核酸代謝等生物過程相關(guān),KEGG分析結(jié)果顯示顯著富集于DNA復制、細胞周期更迭等通路,基因組處于活躍狀態(tài),明顯區(qū)別于其他供體細胞。在顯著表達上調(diào)基因中,基因和胚胎干細胞轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò)相關(guān),其轉(zhuǎn)錄因子P52951可作為胚胎多能性因子[52,53],而多能性因子是影響克隆胚胎發(fā)育潛力的重要因素?;蚴且环N蛋白質(zhì)編碼基因,其功能主要是調(diào)控外胚層分化,在細胞重編程中發(fā)揮關(guān)鍵作用[54,55]。基因的相關(guān)途徑中有網(wǎng)格蛋白介導的內(nèi)吞作用,值得注意的是,在小鼠克隆胚2細胞期停滯胚胎單細胞轉(zhuǎn)錄組測序研究中,也報道了與內(nèi)吞作用相關(guān)的基因激活不足的現(xiàn)象[56],說明與內(nèi)吞途徑相關(guān)基因可能影響SCNT胚胎的發(fā)育效率。在體細胞重編程過程中,大量基因能否被成功激活直接影響克隆胚胎的發(fā)育命運。同時,已有研究表明供體細胞基因表達模式影響克隆效率[57]。因此,在供體細胞基因表達、核重編程、胚胎關(guān)鍵基因激活之間一定存在微妙的聯(lián)系。DW99_2活躍基因組狀態(tài),可能有利細胞重編程,但是需要在不同情況下進行區(qū)分:先前關(guān)于SCNT胚胎重編程異常原因的報道主要集中于抑制性的表觀遺傳修飾[58,59],表現(xiàn)為協(xié)調(diào)胚胎發(fā)育所需基因失敗,部分基因表達受到抑制,基因組不活躍。但是在克隆牛的研究中,供體細胞基因組的異常激活狀態(tài)會阻礙克隆胚胎發(fā)育[18]。在克隆小鼠中,供體細胞異常激活的基因卻會在胚胎發(fā)育中被成功重編程[60],說明供體細胞的基因組狀態(tài)不完全決定SCNT胚胎的發(fā)育命運,不同物種之間是存在差異的,需要進一步研究豬供體細胞基因組狀態(tài)與克隆胚胎發(fā)育命運的聯(lián)系。

    目前,關(guān)于供體細胞異質(zhì)性對核移植胚胎發(fā)育效率的影響還有很多問題亟待解決。就研究廣泛性而言,類似研究鮮有報道,同時,供體細胞異質(zhì)性受物種、遺傳變異、環(huán)境差異和物理刺激等多因素影響,需要更全面的研究來闡述造成這種差異的分子機制以及對核移植胚胎發(fā)育的影響。就研究深度而言,一方面需要有參考價值的精英供體細胞、克隆胚胎的遺傳信息,另一方面,需要建立一種體系或技術(shù),輔助研究人員根據(jù)供體細胞的分子特征篩選潛在精英供體細胞,并進行核移植實驗來驗證。在未來,可利用豬克隆胚胎活檢技術(shù)(在2細胞、4細胞期對每個胚胎抽取一個卵裂球用于單細胞測序,每個胚胎剩余的部分繼續(xù)培養(yǎng)觀察其能否發(fā)育至囊胚,從而把抽取的卵裂球分成發(fā)育正常和異常兩組樣本)結(jié)合單細胞測序技術(shù)確定發(fā)育正常及異常的克隆胚胎的基因表達模式,通過分析本研究所獲得的4個“另類”供體細胞與發(fā)育正常的克隆胚胎的基因表達模式的關(guān)聯(lián),從4個“另類”細胞中篩選出潛在的精英供體細胞,根據(jù)該精英供體細胞的基因表達特征設(shè)計可識別該精英供體細胞的探針或抗體,然后通過流式細胞儀分選等技術(shù)在供體細胞中篩選和富集精英供體細胞進行核移植,驗證利用該精英供體細胞制備的克隆胚胎的發(fā)育是否高于普通供體細胞。最近,Li等[56]利用單細胞測序技術(shù),通過比較核移植2細胞期正常發(fā)育胚胎、停滯發(fā)育胚胎,4細胞期正常發(fā)育胚胎、停滯發(fā)育胚胎的轉(zhuǎn)錄組數(shù)據(jù),找出影響小鼠克隆胚胎重編程的潛在通路和關(guān)鍵基因,為在細胞水平研究克隆胚胎異常發(fā)育的原因提供了參考。相信隨著單細胞測序技術(shù)的不斷成熟,該技術(shù)在克隆研究中的應(yīng)用有望取得進一步突破:一方面,單細胞多組學聯(lián)合分析的發(fā)展有助于建立“基因組–轉(zhuǎn)錄組–代謝組–表型”的系統(tǒng)研究網(wǎng)路;另一方面,更多種類、更多數(shù)量的供體細胞命運將通過單細胞測序技術(shù)被追蹤,形成“供體細胞–2細胞–4細胞–8細胞–囊胚–移植后胚胎”的完整研究體系,描繪供體細胞核和卵母細胞質(zhì)相互作用的動態(tài)過程,綜合分析影響克隆胚胎核移植效率的分子機制,進而提高克隆效率。

    [1] Wilmut I, Schnieke AE, Mcwhir J, Kind AJ, Campbell KH. Viable offspring derived from fetal and adult mammalian cells., 1997, 385(6619): 810–813.

    [2] Cibelli JB, Stice SL, Golueke PJ, Kane JJ, Jerry J, Blackwell C, de León FAP, Robl JM. Cloned transgenic calves produced from nonquiescent fetal fibroblasts., 1998, 280(5367): 1256–1258.

    [3] Onishi A, Iwamoto M, Akita T, Mikawa S, Takeda K, Awata T, Hanada H, Perry AC. Pig cloning by microinjection of fetal fibroblast nuclei., 2000, 289(5482): 1188–1190.

    [4] Liu Z, Cai YJ, Wang Y, Nie YH, Zhang CC, Xu YT, Zhang XT, Lu Y, Wang ZY, Poo M, Sun Q. Cloning of macaque monkeys by somatic cell nuclear transfer., 2018, 172(4): 881–887.

    [5] Beyhan Z, Iager AE, Cibelli JB. Interspecies nuclear transfer: implications for embryonic stem cell biology., 2007, 1(5): 502–512.

    [6] Niemann H, Lucas-Hahn A. Somatic cell nuclear transfer cloning: practical applications and current legislation., 2012, 47(Suppl. 5): 2–10.

    [7] Tachibana M, Amato P, Sparman M, Gutierrez N M, Tippner-Hedges R, Ma H, Kang E, Fulati A, Lee H S, Sritanaudomchai H, Masterson K, Larson J, Eaton D, Sadler-Fredd K, Battaglia D, Lee D, Wu D, Jensen J, Patton P, Gokhale S, Stouffer RL, Wolf D, Mitalipov S. Human embryonic stem cells derived by somatic cell nuclear transfer., 2013, 153(6): 1228–1238.

    [8] Ao Z, Wu X, Zhou J, Gu T, Wang XW, Shi JS, Zhao CF, Cai GY, Zheng EQ, Liu DW, Wu ZF, Li ZC. Cloned pig fetuses exhibit fatty acid deficiency from impaired placental transport., 2019, 86(11): 1569–1581.

    [9] Liu Y, Li J, L?vendahl P, Schmidt M, Larsen K, Callesen H. In vitro manipulation techniques of porcine embryos: a meta-analysis related to transfers, pregnancies and piglets., 2015, 27(3): 429–439.

    [10] Jang G, Park ES, Cho JK, Bhuiyan MM, Lee BC, Kang SK, Hwang WS. Preimplantational embryo development and incidence of blastomere apoptosis in bovine somatic cell nuclear transfer embryos reconstructed with long-term cultured donor cells., 2004, 62(3–4): 512–521.

    [11] Zhang DF, Liu D, Tang LL, Wang Y, Chen Y, Wang K, Wang GL, Schellander K, Cailu L. Effects of different donor cells on the development of nuclear-trans-ferred porcine embryos., 2007, 29(2): 211–217.張德福, 劉東, 湯琳琳, 王英, 陳茵, 王凱, 王根林, KARL Schellander, LIN Cailu. 不同供體細胞及其處理對豬核移植重構(gòu)胚體外發(fā)育的影響. 遺傳, 2007, 29(2): 211–217.

    [12] Ruan ZY, Zhao X, Li ZD, Qin XL, Shao QM, Ruan QY, Deng YF, Jiang JR, Huang B, Lu FH, Shi DS. Effect of sex differences in donor foetal fibroblast on the early development and DNA methylation status of buffalo (Bubalus bubalis) nuclear transfer embryos., 2019, 54(1): 11–22.

    [13] Rideout WR, Eggan K, Jaenisch R. Nuclear cloning and epigenetic reprogramming of the genome., 2001, 293(5532): 1093–1098.

    [14] Yang XQ, Wu ZF, Li ZC. Advances in epigenetic reprogramming of somatic cells nuclear transfer in mam-mals., 2019, 41(12): 1099–1109.楊旭瓊, 吳珍芳, 李紫聰. 哺乳動物體細胞核移植表觀遺傳重編程研究進展. 遺傳, 2019, 41(12): 1099–1109.

    [15] Zhai YH, Li W, Zhang ZR, Cao YQ, Wang ZZ, Zhang S, Li ZY. Epigenetic states of donor cells significantly affect the development of somatic cell nuclear transfer (SCNT) embryos in pigs., 2018, 85(1): 26–37.

    [16] Yamanaka S. Elite and stochastic models for induced pluripotent stem cell generation., 2009, 460(7251): 49–52.

    [17] Inoue K, Ogonuki N, Mochida K, Yamamoto Y, Takano K, Kohda T, Ishino F, Ogura A. Effects of donor cell type and genotype on the efficiency of mouse somatic cell cloning.,2003, 69(4): 1394–1400.

    [18] Zhou C, Zhang JC, Zhang M, Wang DB, Ma Y, Wang Y, Wang YZ, Huang YM, Zhang Y. Transcriptional memory inherited from donor cells is a developmental defect of bovine cloned embryos., 2020, 34(1): 1637–1651.

    [19] Picelli S, Faridani OR, Bj?rklund AK, Winberg G, Sagasser S, Sandberg R. Full-length RNA-seq from single cells using Smart-seq2., 2014, 9(1): 171–181.

    [20] Grubman A, Chew G, Ouyang JF, Sun GZ, Choo XY, Mclean C, Simmons RK, Buckberry S, Vargas-Landin DB, Poppe D, Pflueger J, Lister R, Rackham O, Petretto E, Polo JM. A single-cell atlas of entorhinal cortex from individuals with Alzheimer's disease reveals cell-type- specific gene expression regulation.,2019, 22(12): 2087–2097

    [21] Tang FC, Barbacioru C, Wang YZ, Nordman E, Lee C, Xu NN, Wang XH, Bodeau J, Tuch BB, Siddiqui A, Lao K, Surani MA. mRNA-Seq whole-transcriptome analysis of a single cell., 2009, 6(5): 377–382.

    [22] Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements., 2015, 12(4): 357–360.

    [23] Wang LK, Feng ZX, Wang X, Wang XW, Zhang XG. DEGseq: an R package for identifying differentially expressed genes from RNA-seq data., 2010, 26(1): 136–138.

    [24] Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias., 2010, 11(2): R14.

    [25] Robinson MD, Mccarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data., 2010, 26(1): 139–140.

    [26] Navin N, Hicks J. Future medical applications of single-cell sequencing in cancer., 2011, 3(5): 31.

    [27] Novick A, Weiner M. Enzyme induction as an all-or-none phenomenon., 1957, 43(7): 553–566.

    [28] Coskun AF, Eser U, Islam S. Cellular identity at the single-cell level., 2016, 12(10): 2965–2979.

    [29] Menon M, Mohammadi S, Davila-Velderrain J, Goods BA, Cadwell TD, Xing Y, Stemmer-Rachamimov A, Shalek AK, Love JC, Kellis M, Hafler BP. Single-cell transcriptomic atlas of the human retina identifies cell types associated with age-related macular degeneration.,2019, 10(1): 4902.

    [30] Takahashi K, Tanabe K, Ohnuki M, Narita M, Ichisaka T, Tomoda K, Yamanaka S. Induction of pluripotent stem cells from adult human fibroblasts by defined factors.,2007, 131(5): 861–872.

    [31] Lowry WE, Richter L, Yachechko R, Pyle AD, Tchieu J, Sridharan R, Clark AT, Plath K. Generation of human induced pluripotent stem cells from dermal fibroblasts., 2008, 105(8): 2883–2888.

    [32] Polo JM, Liu S, Figueroa M E, Kulalert W, Eminli S, Tan KY, Apostolou E, Stadtfeld M, Li YS, Shioda T, Natesan S, Wagers AJ, Melnick A, Evans T, Hochedlinger K. Cell type of origin influences the molecular and functional properties of mouse induced pluripotent stem cells., 2010, 28(8): 848–855.

    [33] Lai L, Tao T, Macháty Z, Kühholzer B, Sun QY, Park KW, Day BN, Prather RS. Feasibility of producing porcine nuclear transfer embryos by using G2/M-stage fetal fibroblasts as donors.,2001, 65(5): 1558– 1564.

    [34] Chesné P, Adenot PG, Viglietta C, Baratte M, Boulanger L, Renard JP. Cloned rabbits produced by nuclear transfer from adult somatic cells.,2002, 20(4): 366–369.

    [35] Sutovsky P, Moreno RD, Ramalho-Santos J, Dominko T, Simerly C, Schatten G. Ubiquitin tag for sperm mitochondria.,1999, 402(6760): 371–372.

    [36] Onishi A, Iwamoto M, Akita T, Mikawa S, Takeda K, Awata T, Hanada H, Perry AC. Pig cloning by microinjection of fetal fibroblast nuclei.,2000, 289(5482): 1188–1190.

    [37] Gupta S, Sahu D, Bomalaski JS, Frank I, Boorjian SA, Thapa P, Cheville J C, Hansel DE. Argininosuccinate Synthetase-1 (ASS1) loss in high-grade neuroendocrine carcinomas of the urinary bladder: implications for targeted therapy with ADI-PEG 20., 2018, 29(3): 236–241.

    [38] Moren L, Perryman R, Crook T, Langer JK, Oneill K, Syed N, Antti H. Metabolomic profiling identifies distinct phenotypes for ASS1 positive and negative GBM.,2018, 18(1): 167.

    [39] Bang J, Lee E, Lee AR, Il Lee J, Choi SH, Seol D, Park C, Lee DR. The effect of cell penetrating peptide-conjugated coactivator-associated arginine methyltransferase 1 (CPP-CARM1) on the cloned mouse embryonic development.,2018, 8(1): 16721.

    [40] Abdel-Hady AE, Beige J, Kreutz R, Bolbrinker J. Effect of UMOD genotype on long-term graft survival after kidney transplantation in patients treated with cyclosporine-based therapy.,2018, 18(2): 227–231.

    [41] Bailie C, Kilner J, Maxwell AP, Mcknight AJ. Development of next generation sequencing panel for UMOD and association with kidney disease.,2017, 12(6): e0178321.

    [42] Davis D, Kannan M, Wattenberg B. Orm/ORMDL proteins: Gate guardians and master regulators.,2018, 70: 3–18.

    [43] Fava RA, Browning JL, Gatumu M, Skarstein K, Bolstad AI. LTBR-pathway in Sjogren's syndrome: CXCL13 levels and B-cell-enriched ectopic lymphoid aggregates in NOD mouse lacrimal glands are dependent on LTBR., 2011, 691: 383–390.

    [44] Zhu QQ, Li N, Li F, Sang J, Deng H, Han QY, Lv Y, Li CY, Liu ZW. Association of LTBR polymorphisms with chronic hepatitis B virus infection and hepatitis B virus-related hepatocellular carcinoma., 2017, 49: 126–131.

    [45] Sia D, Losic B, Moeini A, Cabellos L, Hao K, Revill K, Bonal D, Miltiadous O, Zhang ZY, Hoshida Y, Cornella H, Castillo-Martin M, Pinyol R, Kasai Y, Roayaie S, Thung S N, Fuster J, Schwartz ME, Waxman S, Cordon-Cardo C, Schadt E, Mazzaferro V, Llovet JM. Massive parallel sequencing uncovers actionable FGFR2-PPHLN1 fusion and ARAF mutations in intrahepatic cholangiocarcinoma.,2015, 6: 6087.

    [46] Chung YG, Matoba S, Liu YT, Eum JH, Lu FL, Jiang W, Lee JE, Sepilian V, Cha KY, Lee DR, Zhang Y. Histone demethylase expression enhances human somatic cell nuclear transfer efficiency and promotes derivation of pluripotent stem cells., 2015, 17(6): 758–766.

    [47] Boone PM, Yuan B, Gu S, Ma ZW, Gambin T, Gonzaga-Jauregui C, Jain M, Murdock TJ, White JJ, Jhangiani SN, Walker K, Wang QY, Muzny DM, Gibbs RA, Hejtmancik JF, Lupski JR, Posey JE, Lewis RA. Hutterite-type cataract maps to chromosome 6p21.32- p21.31, cosegregates with a homozygous mutation in LEMD2, and is associated with sudden cardiac death., 2016, 4(1): 77–94.

    [48] Mcgee LJ, Jiang AL, Lan Y. Golga5 is dispensable for mouse embryonic development and postnatal survival., 2017, 55(7): 10.1002/dvg.23039.

    [49] O'Grady GL, Best HA, Sztal TE, Schartner V, Sanjuan-Vazquez M, Donkervoort S, Abath NO, Sutton RB, Ilkovski B, Romero NB, Stojkovic T, Dastgir J, Waddell L B, Boland A, Hu Y, Williams C, Ruparelia AA, Maisonobe T, Peduto AJ, Reddel SW, Lek M, Tukiainen T, Cummings BB, Joshi H, Nectoux J, Brammah S, Deleuze JF, Ing VO, Ramm G, Ardicli D, Nowak KJ, Talim B, Topaloglu H, Laing NG, North K N, Macarthur DG, Friant S, Clarke NF, Bryson-Richardson RJ, B?nnemann CG, Laporte J, Cooper ST. Variants in the oxidoreductase PYROXD1 cause early-onset myopathy with internalized nuclei and myofibrillar disorganization.,2016, 99(5): 1086–1105.

    [50] Al-Dabbagh N, Al-Shahrani H, Al-Dohayan N, Mustafa M, Arfin M, Al-Asmari AK. The SPARC-related modular calcium binding protein 2 (SMOC2) gene polymorphism in primary glaucoma: a case-control study., 2017, 11: 549–555.

    [51] Huang XQ, Zhou ZQ, Zhang XF, Chen CL, Tang Y, Zhu Q, Zhang JH, Xia JC. Overexpression of SMOC2 attenuates the tumorigenicity of hepatocellular carcinoma cells and is associated with a positive postoperative prognosis in human hepatocellular carcinoma.,2017, 8(18): 3812–3827.

    [52] Hagan N, Guarente J, Ellisor D, Zervas M. The temporal contribution of the Gbx2 lineage to cerebellar neurons., 2017, 11: 50.

    [53] Mallika C, Guo QX, Li JYH. Gbx2 is essential for maintaining thalamic neuron identity and repressing habenular characters in the developing thalamus., 2015, 407(1): 26–39.

    [54] Alber M, Kalscheuer VM, Marco E, Sherr E, Lesca G, Till M, Gradek G, Wiesener A, Korenke C, Mercier S, Becker F, Yamamoto T, Scherer SW, Marshall CR, Walker S, Dutta UR, Dalal A B, Suckow V, Jamali P, Kahrizi K, Najmabadi H, Minassian BA. ARHGEF9 disease: Phenotype clarification and genotype-phenotype correlation., 2017, 3(3): e148.

    [55] Klein KM, Pendziwiat M, Eilam A, Gilad R, Blatt I, Rosenow F, Kanaan M, Helbig I, Afawi Z. The phenotypic spectrum of ARHGEF9 includes intellectual disability, focal epilepsy and febrile seizures.,2017, 264(7): 1421–1425.

    [56] Li HS, Song MM, Yang W, Cao P, Zheng L, Zuo YC. A Comparative analysis of single-cell transcriptome identifies reprogramming driver factors for efficiency improvement., 2020, 19: 1053–1064.

    [57] Inoue K, Ogonuki N, Mochida K, Yamamoto Y, Takano K, Kohda T, Ishino F, Ogura A. Effects of donor cell type and genotype on the efficiency of mouse somatic cell cloning., 2003, 69(4): 1394–1400.

    [58] Xie BT, Zhang H, Wei RY, Li QN, Weng XG, Kong QR, Liu ZH. Histone H3 lysine 27 trimethylation acts as an epigenetic barrier in porcine nuclear reprogramming., 2016, 151(1): 9–16.

    [59] Matoba S, Liu YT, Lu FL, Iwabuchi KA, Shen L, Inoue A, Zhang Y. Embryonic development following somatic cell nuclear transfer impeded by persisting histone methylation., 2014, 159(4): 884–895.

    [60] Liu Y, Wu FR, Zhang L, Wu XQ, Li DK, Xin J, Xie J, Kong F, Wang WY, Wu QQ, Zhang D, Wang R, Gao SR, Li WY. Transcriptional defects and reprogramming barriers in somatic cell nuclear reprogramming as revealed by single-embryo RNA sequencing., 2018, 19(1): 734.

    Transcriptome heterogeneity of porcine ear fibroblast and its potential influence on embryo development in nuclear transplantation

    Jun Zhou1, Chengcheng Zhao1, Xiao Wu1, Junsong Shi2, Rong Zhou2, Zhenfang Wu1, Zicong Li1

    There is heterogeneity among donor cells of the same source.Many studies have shown that donor cell affects the efficiency of somatic cell nuclear transfer (SCNT).However, the potential influence of donor cell heterogeneity on the efficiency of nuclear transplantation were rarely analyzed at the single-cell level.In this study, single-cell transcriptome sequencing was performed on 52 porcine ear fibroblasts randomly selected from the same source to compare their gene expression patterns. The results showed that 48 cells had similar gene expression patterns, whereas 4 cells (D11_1, D12_1, DW61_2, DW99_2) had significantly different gene expression patterns from those of other cells. There were no two cells with identical gene expression patterns. The gene expression patterns of D11_1, D12_1, DW61_2 and DW99_2 were analyzed, using the 48 cells with similar gene expression patterns as controls. Firstly, we used the R language statistics to select the differentially expressed genes in the 4 single cells, and identified the top 50 most significant differentially expressed genes. Then GO enrichment analysis and KEGG pathway analysis were performed on the differentially expressed genes. Enrichment analysis revealed that the main molecular functions of the differentially expressed genes included energy metabolism, protein metabolism and cell response to stimulation. The main pathways from KEGG enrichment were related to cell cycle, cell metabolism, and DNA replication. Finally, based on the above results and in consideration with the SCNT research progress, we discussed the potential effects of differential gene expression patterns of the 4 single cells on the embryonic development efficiency of nuclear transplantation. This study revealed transcriptional heterogeneity of porcine ear tissue fibroblasts and provided an effective method to analyze elite donor cells, thereby providing new ideas on improving the cloning efficiency of SCNT.

    pigs; SCNT; donor cell; heterogeneity; single-cell sequencing

    2020-06-22;

    2020-08-11

    廣東省鄉(xiāng)村振興戰(zhàn)略專項(2018年)資助[Supported by Guangdong Provincial Promotion Project on Preservation and Utilization of Local Breed of Livestock and Poultry (2018)]

    周俊,在讀碩士研究生,專業(yè)方向:動物遺傳育種與繁殖。E-mail: 374949357@qq.com

    李紫聰,教授,博士生導師,研究方向:動物遺傳育種與繁殖。E-mail: lizicongcong@163.com

    10.16288/j.yczz.20-190

    2020/9/2 16:08:53

    URI: https://kns.cnki.net/kcms/detail/11.1913.r.20200831.1604.002.html

    (責任編委: 趙要風)

    猜你喜歡
    單細胞差異基因體細胞
    ICR鼠肝和腎毒性損傷生物標志物的篩選
    浙江:誕生首批體細胞克隆豬
    人工智能助力微生物單細胞鑒定
    科學(2020年4期)2020-11-26 08:27:16
    新型冠狀病毒入侵人體細胞之謎
    科學(2020年4期)2020-11-26 08:27:10
    基于RNA 測序研究人參二醇對大鼠心血管內(nèi)皮細胞基因表達的影響 (正文見第26 頁)
    內(nèi)皮前體細胞亞型與偏頭痛的相關(guān)性分析
    非洲菊花托的體細胞胚發(fā)生及植株再生
    聯(lián)合SNaPshot和單倍型分析技術(shù)建立G6PD缺乏癥單細胞基因診斷體系
    單細胞測序技術(shù)研究進展
    SSH技術(shù)在絲狀真菌功能基因篩選中的應(yīng)用
    夜夜爽天天搞| 国产亚洲欧美98| 免费看日本二区| 亚洲av中文字字幕乱码综合| 我的老师免费观看完整版| 国产伦精品一区二区三区视频9| 一级毛片久久久久久久久女| 好男人在线观看高清免费视频| 在现免费观看毛片| 简卡轻食公司| 性色avwww在线观看| www.www免费av| 午夜精品久久久久久毛片777| 国产精品美女特级片免费视频播放器| 很黄的视频免费| 久久久国产成人精品二区| 内地一区二区视频在线| 又黄又爽又刺激的免费视频.| 亚洲av成人精品一区久久| 中文字幕av在线有码专区| 网址你懂的国产日韩在线| 大型黄色视频在线免费观看| 舔av片在线| 成人av在线播放网站| 国产精品一区二区免费欧美| 欧美日韩黄片免| 美女黄网站色视频| 麻豆成人午夜福利视频| 午夜视频国产福利| 亚洲乱码一区二区免费版| 淫秽高清视频在线观看| 成人鲁丝片一二三区免费| 在线播放无遮挡| 久久热精品热| 日韩欧美免费精品| 如何舔出高潮| 在线观看66精品国产| 人人妻,人人澡人人爽秒播| 国产精品综合久久久久久久免费| 亚洲精品亚洲一区二区| av国产免费在线观看| АⅤ资源中文在线天堂| 国内毛片毛片毛片毛片毛片| 国产精品三级大全| 国产色爽女视频免费观看| 国产综合懂色| 午夜精品一区二区三区免费看| 日韩欧美在线乱码| 亚洲av二区三区四区| 国产午夜福利久久久久久| 国产野战对白在线观看| 亚洲激情在线av| a级毛片免费高清观看在线播放| 免费人成视频x8x8入口观看| 久久久色成人| 成人亚洲精品av一区二区| 国产色婷婷99| 亚洲国产精品合色在线| 国产高潮美女av| 国产高清激情床上av| 亚洲狠狠婷婷综合久久图片| 亚洲av免费高清在线观看| 亚洲人成网站高清观看| 久久精品91蜜桃| 精品乱码久久久久久99久播| 国产真实乱freesex| 国产aⅴ精品一区二区三区波| 免费av不卡在线播放| 最新中文字幕久久久久| 亚洲国产精品成人综合色| 亚洲av电影不卡..在线观看| 女同久久另类99精品国产91| 亚洲av成人av| 99久久精品国产亚洲精品| 好看av亚洲va欧美ⅴa在| 亚洲乱码一区二区免费版| 一级毛片久久久久久久久女| www.熟女人妻精品国产| 久久久久久九九精品二区国产| 琪琪午夜伦伦电影理论片6080| 一进一出抽搐动态| 亚洲最大成人中文| 欧美bdsm另类| 桃红色精品国产亚洲av| 国产高清三级在线| 好男人电影高清在线观看| 噜噜噜噜噜久久久久久91| 97超级碰碰碰精品色视频在线观看| 搡老熟女国产l中国老女人| 精品人妻视频免费看| 岛国在线免费视频观看| 久久精品国产亚洲av天美| 免费av不卡在线播放| 国产私拍福利视频在线观看| 久99久视频精品免费| 中文字幕高清在线视频| 欧美在线黄色| 青草久久国产| 看黄色毛片网站| 亚洲精品一区av在线观看| 性色av乱码一区二区三区2| 欧美极品一区二区三区四区| 日韩中文字幕欧美一区二区| 亚洲精品粉嫩美女一区| 亚洲精品456在线播放app | 老司机深夜福利视频在线观看| 亚洲精品日韩av片在线观看| 嫁个100分男人电影在线观看| 成人美女网站在线观看视频| 日本五十路高清| 一个人免费在线观看电影| 色av中文字幕| 又紧又爽又黄一区二区| 国产精品三级大全| 亚洲男人的天堂狠狠| 中文亚洲av片在线观看爽| 亚洲熟妇中文字幕五十中出| 久久久久性生活片| h日本视频在线播放| 美女黄网站色视频| 桃色一区二区三区在线观看| 久久香蕉精品热| 国产成人福利小说| 欧美极品一区二区三区四区| 国产高潮美女av| 欧美性猛交╳xxx乱大交人| 久久草成人影院| 桃色一区二区三区在线观看| 久久国产精品影院| 亚洲狠狠婷婷综合久久图片| 黄色丝袜av网址大全| 国产一区二区三区视频了| 琪琪午夜伦伦电影理论片6080| 国产69精品久久久久777片| 国产精品国产高清国产av| 宅男免费午夜| 亚洲人与动物交配视频| 网址你懂的国产日韩在线| 亚洲精品色激情综合| 久久久色成人| 亚洲不卡免费看| 日本精品一区二区三区蜜桃| 99国产综合亚洲精品| 热99re8久久精品国产| 免费av观看视频| 一进一出抽搐gif免费好疼| 久久婷婷人人爽人人干人人爱| 免费高清视频大片| 天天躁日日操中文字幕| 久久草成人影院| 91在线观看av| 亚洲专区中文字幕在线| 97超视频在线观看视频| 嫩草影院入口| 男人狂女人下面高潮的视频| 中文字幕高清在线视频| 久久久久久九九精品二区国产| 五月玫瑰六月丁香| 岛国在线免费视频观看| 又黄又爽又刺激的免费视频.| 日韩欧美国产一区二区入口| 18禁裸乳无遮挡免费网站照片| 窝窝影院91人妻| 国产午夜精品久久久久久一区二区三区 | 国产大屁股一区二区在线视频| 在线观看66精品国产| 久久精品夜夜夜夜夜久久蜜豆| 亚洲一区二区三区色噜噜| 亚洲经典国产精华液单 | 琪琪午夜伦伦电影理论片6080| 男人的好看免费观看在线视频| 亚洲精品乱码久久久v下载方式| 成人国产综合亚洲| 五月伊人婷婷丁香| 99国产综合亚洲精品| aaaaa片日本免费| 中文字幕av在线有码专区| 欧美一区二区精品小视频在线| 日本精品一区二区三区蜜桃| 亚洲五月婷婷丁香| 国产大屁股一区二区在线视频| 国产老妇女一区| 久久精品国产亚洲av香蕉五月| 老鸭窝网址在线观看| 亚洲avbb在线观看| 日韩精品中文字幕看吧| 日韩中字成人| 一个人观看的视频www高清免费观看| 99热这里只有精品一区| 精品福利观看| 久久久成人免费电影| 成人午夜高清在线视频| 蜜桃久久精品国产亚洲av| 又紧又爽又黄一区二区| 丰满人妻一区二区三区视频av| 久久久国产成人免费| 国产欧美日韩精品一区二区| av在线观看视频网站免费| 午夜亚洲福利在线播放| 免费黄网站久久成人精品 | 两个人视频免费观看高清| 欧美成人免费av一区二区三区| 激情在线观看视频在线高清| avwww免费| 国内毛片毛片毛片毛片毛片| 日韩 亚洲 欧美在线| 中文字幕高清在线视频| 成人特级黄色片久久久久久久| 宅男免费午夜| 禁无遮挡网站| 又黄又爽又刺激的免费视频.| 亚洲第一区二区三区不卡| 波多野结衣高清作品| 久久午夜亚洲精品久久| 欧美一级a爱片免费观看看| 国产午夜福利久久久久久| 一二三四社区在线视频社区8| 成年人黄色毛片网站| 亚洲在线观看片| 中文亚洲av片在线观看爽| 黄片小视频在线播放| 18禁裸乳无遮挡免费网站照片| 国产精品一区二区免费欧美| 国内精品美女久久久久久| av在线天堂中文字幕| 亚洲久久久久久中文字幕| 日韩精品中文字幕看吧| 亚洲精品粉嫩美女一区| 又黄又爽又刺激的免费视频.| 亚洲国产精品999在线| 青草久久国产| av福利片在线观看| 国产高清有码在线观看视频| 免费av观看视频| 免费电影在线观看免费观看| 一级av片app| 国产一级毛片七仙女欲春2| 性色av乱码一区二区三区2| 免费在线观看日本一区| 亚洲激情在线av| 国产成人欧美在线观看| 老熟妇仑乱视频hdxx| 欧美黑人巨大hd| 亚洲中文日韩欧美视频| 哪里可以看免费的av片| av黄色大香蕉| 亚洲成人精品中文字幕电影| 能在线免费观看的黄片| 国产极品精品免费视频能看的| 免费大片18禁| 午夜视频国产福利| www.999成人在线观看| 老司机午夜福利在线观看视频| 999久久久精品免费观看国产| 舔av片在线| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美国产一区二区入口| 亚洲国产色片| 熟女电影av网| 一个人观看的视频www高清免费观看| 日日摸夜夜添夜夜添小说| 日韩成人在线观看一区二区三区| 天天躁日日操中文字幕| 人妻制服诱惑在线中文字幕| 天堂动漫精品| 国内久久婷婷六月综合欲色啪| 精品人妻1区二区| 一本一本综合久久| 免费av不卡在线播放| 亚洲性夜色夜夜综合| 中文在线观看免费www的网站| av中文乱码字幕在线| 国产亚洲精品久久久com| 日韩国内少妇激情av| 国产一区二区在线av高清观看| 久久久久精品国产欧美久久久| 美女大奶头视频| 一级作爱视频免费观看| 波野结衣二区三区在线| 亚洲熟妇中文字幕五十中出| 成年人黄色毛片网站| 小蜜桃在线观看免费完整版高清| 亚洲欧美日韩无卡精品| 在线观看午夜福利视频| 国产精品国产高清国产av| 黄色视频,在线免费观看| 国产精品电影一区二区三区| 欧美精品国产亚洲| 久久久久久久亚洲中文字幕 | 国产综合懂色| 亚洲av.av天堂| 婷婷丁香在线五月| 丁香欧美五月| 麻豆av噜噜一区二区三区| 国内精品美女久久久久久| 欧美成人一区二区免费高清观看| 国产精品国产高清国产av| 18禁在线播放成人免费| 色综合婷婷激情| 一二三四社区在线视频社区8| 日韩欧美国产在线观看| 丰满的人妻完整版| 级片在线观看| 我的女老师完整版在线观看| 亚洲久久久久久中文字幕| 久久国产乱子伦精品免费另类| 非洲黑人性xxxx精品又粗又长| 国产极品精品免费视频能看的| 久久午夜福利片| 在线十欧美十亚洲十日本专区| 欧美日韩中文字幕国产精品一区二区三区| 91麻豆av在线| 欧美日韩黄片免| 国内毛片毛片毛片毛片毛片| 国产成人啪精品午夜网站| 在线观看免费视频日本深夜| 女生性感内裤真人,穿戴方法视频| 欧美黄色淫秽网站| 老熟妇乱子伦视频在线观看| 亚洲综合色惰| 国产亚洲av嫩草精品影院| 在线播放无遮挡| 国产欧美日韩一区二区三| 国产人妻一区二区三区在| 内射极品少妇av片p| 精品乱码久久久久久99久播| 精品无人区乱码1区二区| 婷婷精品国产亚洲av| 国产一区二区在线av高清观看| 成人亚洲精品av一区二区| 日本黄大片高清| 在线观看美女被高潮喷水网站 | 国产精品久久久久久久电影| or卡值多少钱| 十八禁国产超污无遮挡网站| 亚洲精品一区av在线观看| www.www免费av| 看黄色毛片网站| 国产一区二区在线观看日韩| 制服丝袜大香蕉在线| 中文在线观看免费www的网站| 亚洲成人久久性| 狠狠狠狠99中文字幕| 久久国产精品人妻蜜桃| 村上凉子中文字幕在线| 桃红色精品国产亚洲av| 亚洲无线观看免费| 亚洲精品色激情综合| 首页视频小说图片口味搜索| 夜夜看夜夜爽夜夜摸| 乱人视频在线观看| 亚洲天堂国产精品一区在线| 婷婷精品国产亚洲av在线| 一边摸一边抽搐一进一小说| 亚洲中文日韩欧美视频| 精品国产三级普通话版| 狂野欧美白嫩少妇大欣赏| 亚洲av二区三区四区| 欧美性猛交黑人性爽| 欧洲精品卡2卡3卡4卡5卡区| 国产伦精品一区二区三区视频9| 免费黄网站久久成人精品 | 精品久久久久久久久久免费视频| www.999成人在线观看| 夜夜躁狠狠躁天天躁| 精品午夜福利在线看| 淫妇啪啪啪对白视频| 深夜精品福利| 99精品在免费线老司机午夜| 久久久久久久午夜电影| 亚洲成av人片在线播放无| 嫩草影院精品99| 亚洲国产日韩欧美精品在线观看| 乱人视频在线观看| 99国产精品一区二区三区| 午夜亚洲福利在线播放| 99久久九九国产精品国产免费| 色播亚洲综合网| h日本视频在线播放| 亚洲精华国产精华精| 三级男女做爰猛烈吃奶摸视频| 精品久久久久久久末码| 久久精品91蜜桃| 精品欧美国产一区二区三| 我的老师免费观看完整版| 国产淫片久久久久久久久 | 国产成人影院久久av| 成人特级黄色片久久久久久久| 欧美成狂野欧美在线观看| 人妻夜夜爽99麻豆av| 99精品久久久久人妻精品| 日韩人妻高清精品专区| 亚洲成a人片在线一区二区| 村上凉子中文字幕在线| 精品久久久久久成人av| 91麻豆精品激情在线观看国产| 国产大屁股一区二区在线视频| 国产精品,欧美在线| 精品人妻一区二区三区麻豆 | 欧美绝顶高潮抽搐喷水| 久久亚洲真实| 一区二区三区激情视频| 久久久色成人| 波多野结衣高清无吗| 日本撒尿小便嘘嘘汇集6| 午夜影院日韩av| 亚洲第一欧美日韩一区二区三区| 国产一区二区在线观看日韩| 99久久成人亚洲精品观看| 美女黄网站色视频| 国产精品野战在线观看| 亚洲av一区综合| 亚洲国产精品合色在线| a在线观看视频网站| 午夜两性在线视频| 亚洲精华国产精华精| 国产精品人妻久久久久久| 中文资源天堂在线| 色噜噜av男人的天堂激情| 国产美女午夜福利| 亚洲av中文字字幕乱码综合| 我的女老师完整版在线观看| 久久伊人香网站| 亚洲成人精品中文字幕电影| 午夜激情欧美在线| 两个人的视频大全免费| 欧美zozozo另类| 国产探花在线观看一区二区| 男女做爰动态图高潮gif福利片| 成人av在线播放网站| 中文亚洲av片在线观看爽| 精品国产三级普通话版| 色综合欧美亚洲国产小说| 国产精品影院久久| 亚洲欧美日韩高清专用| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美日韩卡通动漫| 好看av亚洲va欧美ⅴa在| 成年女人永久免费观看视频| 男女那种视频在线观看| 热99在线观看视频| 又爽又黄无遮挡网站| 99国产极品粉嫩在线观看| 99精品久久久久人妻精品| 精品国产亚洲在线| 桃红色精品国产亚洲av| h日本视频在线播放| 国语自产精品视频在线第100页| 欧美一区二区亚洲| 免费av不卡在线播放| 久久午夜亚洲精品久久| 男人舔奶头视频| 女人被狂操c到高潮| 国产高清激情床上av| 日本一本二区三区精品| 亚洲狠狠婷婷综合久久图片| 精品福利观看| 亚洲国产色片| 美女大奶头视频| 黄色丝袜av网址大全| 久久人人爽人人爽人人片va | 国产真实伦视频高清在线观看 | 一级黄色大片毛片| 久9热在线精品视频| 国产三级在线视频| 久久精品国产亚洲av天美| www.999成人在线观看| 亚洲国产日韩欧美精品在线观看| 老鸭窝网址在线观看| 三级国产精品欧美在线观看| 久久久久性生活片| 午夜a级毛片| 久久久精品大字幕| av在线天堂中文字幕| 黄色日韩在线| 简卡轻食公司| 亚洲成人中文字幕在线播放| 十八禁国产超污无遮挡网站| 听说在线观看完整版免费高清| 老熟妇乱子伦视频在线观看| 亚洲av免费高清在线观看| av天堂在线播放| 亚洲人成网站在线播| 亚洲18禁久久av| 国产熟女xx| 欧美xxxx黑人xx丫x性爽| 日韩欧美 国产精品| 夜夜夜夜夜久久久久| 草草在线视频免费看| 精品久久久久久久久久免费视频| 18禁在线播放成人免费| 伦理电影大哥的女人| 桃红色精品国产亚洲av| 亚洲真实伦在线观看| 亚洲不卡免费看| 国产高清视频在线观看网站| 国产精品99久久久久久久久| www.www免费av| 此物有八面人人有两片| 精品99又大又爽又粗少妇毛片 | 亚洲最大成人手机在线| 亚洲美女黄片视频| 九色国产91popny在线| 国产精品自产拍在线观看55亚洲| 在线a可以看的网站| 成人av在线播放网站| 美女大奶头视频| 国产在线精品亚洲第一网站| 亚洲av一区综合| 日韩中文字幕欧美一区二区| 国产成人av教育| 男女下面进入的视频免费午夜| 日韩大尺度精品在线看网址| 在线观看美女被高潮喷水网站 | 国产成人福利小说| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | x7x7x7水蜜桃| 亚洲综合色惰| 真人做人爱边吃奶动态| 最近中文字幕高清免费大全6 | 嫩草影视91久久| 99热这里只有精品一区| 国产三级黄色录像| 婷婷精品国产亚洲av| 狂野欧美白嫩少妇大欣赏| 国产主播在线观看一区二区| 日日夜夜操网爽| 欧美激情在线99| 亚洲成人久久性| 色av中文字幕| 国产三级黄色录像| 国产白丝娇喘喷水9色精品| 熟女人妻精品中文字幕| av欧美777| 国产精品野战在线观看| 啪啪无遮挡十八禁网站| 欧美+亚洲+日韩+国产| 日韩欧美在线二视频| a级毛片免费高清观看在线播放| 网址你懂的国产日韩在线| 18禁黄网站禁片午夜丰满| 男女做爰动态图高潮gif福利片| 国产在线男女| 国模一区二区三区四区视频| 97碰自拍视频| 中文在线观看免费www的网站| 极品教师在线免费播放| 国产午夜福利久久久久久| 亚洲精品456在线播放app | 能在线免费观看的黄片| 国产高清激情床上av| 亚洲av中文字字幕乱码综合| 午夜福利在线观看吧| 欧美又色又爽又黄视频| 99久久成人亚洲精品观看| 最近视频中文字幕2019在线8| 熟女电影av网| 女人被狂操c到高潮| 欧美精品国产亚洲| 国产精品一区二区三区四区免费观看 | 97碰自拍视频| 久久久久性生活片| 欧美bdsm另类| 精品久久久久久久久久免费视频| 九色成人免费人妻av| 在线国产一区二区在线| 亚洲av第一区精品v没综合| 小说图片视频综合网站| 亚洲五月婷婷丁香| 日本黄大片高清| 成年免费大片在线观看| www日本黄色视频网| 成年女人毛片免费观看观看9| 欧美bdsm另类| 老女人水多毛片| 国产亚洲欧美在线一区二区| 国产成+人综合+亚洲专区| 少妇裸体淫交视频免费看高清| 国产白丝娇喘喷水9色精品| 精品久久久久久久末码| 999久久久精品免费观看国产| 午夜激情欧美在线| 18禁黄网站禁片免费观看直播| 日本在线视频免费播放| 日本精品一区二区三区蜜桃| 成人精品一区二区免费| 日本免费a在线| 国产精品98久久久久久宅男小说| 欧美乱妇无乱码| 欧美+亚洲+日韩+国产| 特大巨黑吊av在线直播| 欧美午夜高清在线| 国产伦精品一区二区三区四那| 淫秽高清视频在线观看| 又黄又爽又免费观看的视频| 美女cb高潮喷水在线观看| 成人毛片a级毛片在线播放| 亚洲成a人片在线一区二区| 色噜噜av男人的天堂激情| 成人美女网站在线观看视频| 淫秽高清视频在线观看| 国产毛片a区久久久久| 成人无遮挡网站| 999久久久精品免费观看国产| 一进一出抽搐动态| 97人妻精品一区二区三区麻豆| 日本在线视频免费播放| 麻豆国产97在线/欧美| 亚洲一区二区三区不卡视频| а√天堂www在线а√下载| 99国产精品一区二区三区| 久久久久九九精品影院| 每晚都被弄得嗷嗷叫到高潮| 美女高潮的动态| 久久国产乱子免费精品|