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

    基于CRISPR/Cas9系統(tǒng)的水稻光敏色素互作因子OsPIL15基因編輯

    2017-09-11 09:05:40季新李飛晏云孫紅正張靜李俊周彭廷杜彥修趙全志
    關(guān)鍵詞:靶位等位突變體

    季新,李飛,晏云,孫紅正,張靜,李俊周,彭廷,杜彥修,趙全志

    ?

    基于CRISPR/Cas9系統(tǒng)的水稻光敏色素互作因子基因編輯

    季新,李飛,晏云,孫紅正,張靜,李俊周,彭廷,杜彥修,趙全志

    (河南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院/河南糧食作物協(xié)同創(chuàng)新中心/河南省水稻生物學(xué)重點(diǎn)實(shí)驗(yàn)室,鄭州 450002)

    【目的】光作為一種環(huán)境信號(hào),可影響植物的基因表達(dá)、酶活性和形態(tài)建成。光敏色素互作因子在光信號(hào)傳導(dǎo)過程中起著重要作用。本研究旨在構(gòu)建水稻光敏色素互作因子的CRISPR/Cas9表達(dá)載體,創(chuàng)制突變體,挖掘水稻功能基因,豐富和完善水稻光信號(hào)調(diào)控分子機(jī)制?!痉椒ā恳罁?jù)CRISPR/Cas9技術(shù)原理,設(shè)計(jì)突變靶點(diǎn)。將所設(shè)計(jì)靶序列在水稻基因組中進(jìn)行比對(duì),排除非特異性靶位點(diǎn),同時(shí)使該靶序列含有常用酶切位點(diǎn),方便后期突變體鑒定?;瘜W(xué)合成靶位點(diǎn)寡核苷酸序列并與載體pBUN411連接構(gòu)建CRISPR/Cas9表達(dá)載體,利用農(nóng)桿菌介導(dǎo)法導(dǎo)入粳稻品種日本晴,以除草劑抗性標(biāo)記篩選獲得陽性轉(zhuǎn)基因植株。利用酶切法判斷T0代轉(zhuǎn)基因植株是否發(fā)生突變,結(jié)合測(cè)序結(jié)果分析突變單株的突變基因型。將靶點(diǎn)序列在水稻全基因組中進(jìn)行比對(duì)分析,選擇5個(gè)與靶序列同源性較高且錯(cuò)配在4 bp以內(nèi)的位點(diǎn)作為潛在脫靶位點(diǎn)進(jìn)行脫靶效應(yīng)評(píng)估,分析所設(shè)計(jì)靶序列特異性?!窘Y(jié)果】所構(gòu)建表達(dá)載體成功實(shí)現(xiàn)了對(duì)的定向編輯,酶切顯示在選取的25株T0代轉(zhuǎn)基因植株中獲得15株突變體,其中包括5株純合突變體、6株雙等位突變體和4株雜合突變體,共10種不同突變基因型和11個(gè)突變株系。突變類型以單堿基插入或缺失為主,同時(shí)也得到2種56和66 bp較大片段缺失株系。對(duì)部分純合突變、雙等位突變和雜合突變體的T1代植株進(jìn)行分析,結(jié)果表明,T0代產(chǎn)生的突變基因型絕大部分能穩(wěn)定遺傳給下一代。T0代純合突變體后代為純合突變單株,僅在株系14純合突變體后代中檢測(cè)到1株未突變單株;T0代雙等位突變體后代可得到2種純合突變型和1種雙等位突變型;T0代雜合突變體后代則可得到純合、雜合及未突變3種類型。對(duì)T0代未突變植株的后繼世代酶切分析顯示,62株T1代轉(zhuǎn)基因植株均未發(fā)生突變,表明CRISPR/Cas9在T1代轉(zhuǎn)基因陽性植株中未重新發(fā)揮基因編輯作用。對(duì)20株突變體的5個(gè)潛在脫靶位點(diǎn)進(jìn)行分析,5個(gè)潛在脫靶位點(diǎn)均未檢測(cè)出脫靶效應(yīng),表明所設(shè)計(jì)靶序列具有較高特異性。對(duì)選取的3組不同基因型T1代突變體表型進(jìn)行初步觀察,結(jié)果表明,突變體生育期和分蘗數(shù)未出現(xiàn)明顯變化,株高極顯著下降,籽粒粒長(zhǎng)極顯著增加,最大增幅達(dá)5.69%?!窘Y(jié)論】CRISPR/Cas9系統(tǒng)能對(duì)進(jìn)行定向編輯,獲得的10種不同突變基因型的突變體與野生型相比株高極顯著降低、籽粒粒長(zhǎng)極顯著增大。

    水稻;CRISPR/Cas9;基因編輯;;脫靶效應(yīng)

    0 引言

    【研究意義】光是作物生長(zhǎng)發(fā)育所需的最重要生態(tài)因子之一,是作物產(chǎn)量形成的基礎(chǔ),同時(shí)光也可作為一種重要的環(huán)境信號(hào),調(diào)節(jié)植物基因的表達(dá)、影響酶的活性以及植物的形態(tài)建成。在光信號(hào)途徑中,光敏色素互作因子起著承上啟下的作用,因此,利用CRISPR/Cas9基因編輯技術(shù)創(chuàng)制光敏色素互作因子的突變體,研究其在光信號(hào)調(diào)控途徑中的作用,對(duì)挖掘水稻功能基因,豐富和完善水稻光信號(hào)調(diào)控分子機(jī)制具有重要意義?!厩叭搜芯窟M(jìn)展】基因編輯技術(shù)指能夠?qū)δ繕?biāo)基因進(jìn)行“編輯”,實(shí)現(xiàn)對(duì)特定DNA片段的敲除、插入等。過去幾年,以鋅指核酸酶(zinc fnger nucleases,ZFNs)和類轉(zhuǎn)錄激活因子效應(yīng)物核酸酶(transcription activator-like nucleases,TALENs)為代表的序列特異性核酸酶技術(shù)可以高效率的定點(diǎn)編輯基因組,在基礎(chǔ)研究、基因治療和遺傳改良等方面發(fā)揮了重要作用[1-2]。近年來,一種名為CRISPR/Cas9(clustered regularly interspaced short palindromic repeats/CRISPR-associated nuclease 9,Cas9)的基因組定向編輯技術(shù)備受矚目。2013年初,《Science》兩篇文章首次報(bào)道了Cas9核酸酶用于人類和小鼠細(xì)胞基因組編輯[3-4],其作為第三代基因編輯技術(shù)受到人們廣泛關(guān)注。CRISPR/Cas9技術(shù)是利用核酸酶Cas9蛋白與單導(dǎo)向RNA(single guide RNA,sgRNA)形成復(fù)合體,sgRNA通過堿基互補(bǔ)配對(duì)決定靶序列特異性,Cas9蛋白作為核酸酶切割與sgRNA上的間隔序列(spacers)互補(bǔ)的基因組DNA,造成雙鏈DNA損傷,隨后通過體內(nèi)的NHEJ(non-homologous end joining)修復(fù)機(jī)制引入基因突變[5]。隨著CRISPR/ Cas9系統(tǒng)在人類和動(dòng)物細(xì)胞中建立與應(yīng)用,CRISPR/ Cas9系統(tǒng)在擬南芥、煙草、水稻、小麥、玉米、高粱、番茄和甜橙等植物中均實(shí)現(xiàn)了基因組的定向編輯[6-7]。光敏色素互作因子PIF(phytochrome-interacting factor)或稱PIL(phytochrome interacting factor-like)是bHLH轉(zhuǎn)錄因子家族中的一類轉(zhuǎn)錄因子,其主要特性是能與光敏色素互作直接或間接調(diào)控光響應(yīng)基因。作為bHLH蛋白的一種,所有PIF家族的蛋白都包含在N端與光敏色素互作的APB(active phytochrome B-binding)或APA(active phytochrome A-binding)結(jié)構(gòu)域和C端bHLH-DNA結(jié)合結(jié)構(gòu)域及核定位結(jié)構(gòu)域[8-9]。擬南芥PIF能夠與靶基因啟動(dòng)子區(qū)域的G-box(5′-CACGTG-3′)或PBE-box(PIF- binding E-box, 5′-CACATG-3′)相結(jié)合[10],可調(diào)控種子萌發(fā)[11]、幼苗去黃化[12]、避陰反應(yīng)[13]和晝夜節(jié)律[14]等植物生長(zhǎng)發(fā)育的各個(gè)方面。Nakamura等[15]通過同源性分析在水稻基因組中鑒定了6個(gè)PIF轉(zhuǎn)錄因子(OsPIL11—OsPIL16),相對(duì)于擬南芥而言,PIFs在水稻生長(zhǎng)發(fā)育過程中的作用研究相對(duì)較少。參與紅光誘導(dǎo)的幼苗去黃化反應(yīng)[16];在水稻中超表達(dá)后能夠促進(jìn)水稻節(jié)間的伸長(zhǎng),反之低表達(dá)則抑制其節(jié)間伸長(zhǎng)[17];超表達(dá)的水稻種子在黑暗環(huán)境中的萌發(fā)受到抑制[18];能夠負(fù)調(diào)控進(jìn)而調(diào)控籽粒大小[19]?!颈狙芯壳腥朦c(diǎn)】筆者前期通過高通量測(cè)序研究發(fā)現(xiàn),在水稻籽粒抽穗后10、15、21、27和35 d這5個(gè)灌漿時(shí)期的強(qiáng)勢(shì)粒和弱勢(shì)粒中,6個(gè)水稻光敏色素互作因子有3個(gè)被鑒定出,其中和只在弱勢(shì)粒某個(gè)時(shí)期表達(dá),則在強(qiáng)弱勢(shì)粒的5個(gè)灌漿時(shí)期均表達(dá)[20],而目前有關(guān)在水稻籽粒灌漿中的調(diào)控作用還不清楚?!緮M解決的關(guān)鍵問題】本研究擬利用CRISPR/Cas9系統(tǒng)對(duì)水稻光敏色素互作因子進(jìn)行編輯,轉(zhuǎn)化粳稻品種日本晴,對(duì)其突變位點(diǎn)和潛在脫靶位點(diǎn)進(jìn)行分析鑒定,通過對(duì)突變體表型進(jìn)行觀察,了解功能及其調(diào)控路徑。

    1 材料與方法

    1.1 試驗(yàn)材料

    使用粳稻品種日本晴(sspcv. Nipponbare)作為轉(zhuǎn)基因受體材料,CRISPR/Cas9載體pBUN411由中國(guó)農(nóng)業(yè)大學(xué)陳其軍教授惠贈(zèng)[21],試驗(yàn)所用引物(表1)和測(cè)序分析由蘇州金唯智生物科技有限公司完成。

    表1 本研究所用引物

    引物序列中小寫字母為Ⅰ限制性內(nèi)切酶的黏性末端接頭 The lowercase stand for sticky ends ofⅠ

    1.2 靶位點(diǎn)設(shè)計(jì)和CRISPR/Cas9表達(dá)載體構(gòu)建

    根據(jù)水稻(GenBank登錄號(hào)AK102252.1)外顯子序列設(shè)計(jì)sgRNA靶點(diǎn)序列。在盡量靠近外顯子編碼蛋白質(zhì)的N端,選取PAM(protospacer-adjacent motif)序列(NGG)5′端的一段堿基序列作為靶位點(diǎn),將該段序列在水稻基因組數(shù)據(jù)庫中進(jìn)行比對(duì),以排除非特異性靶位點(diǎn),同時(shí)盡可能使PAM 5′端存在常用酶切位點(diǎn),以方便后期突變體鑒定。

    將選取的靶點(diǎn)序列作為正義鏈,其互補(bǔ)鏈作為反義鏈,分別在2條DNA單鏈5′端添加ggcg和aaac序列作為Ⅰ限制性內(nèi)切酶的黏性末端接頭。使用限制性內(nèi)切酶Ⅰ酶切CRISPR/Cas9載體pBUN411質(zhì)粒使其線性化,T4 DNA連接酶連接后用熱激法轉(zhuǎn)化大腸桿菌DH5α感受態(tài)細(xì)胞,使用引物pBUN411- VF和pBUN411-VR進(jìn)行菌落PCR驗(yàn)證。CRISPR/Cas9表達(dá)載體轉(zhuǎn)化農(nóng)桿菌EHA105,用除草劑(Basta)篩選獲得再生苗。

    1.3 酶切和測(cè)序分析驗(yàn)證突變位點(diǎn)

    為了檢測(cè)靶位點(diǎn)的突變情況,在靶位點(diǎn)兩側(cè)設(shè)計(jì)引物OsPIL15-F和OsPIL15-R,以轉(zhuǎn)基因陽性單株DNA為模板擴(kuò)增靶位點(diǎn)序列。使用Ⅰ限制性內(nèi)切酶酶切PCR產(chǎn)物,2%瓊脂糖凝膠電泳檢測(cè)是否發(fā)生突變。將突變單株P(guān)CR產(chǎn)物測(cè)序,對(duì)靶位點(diǎn)附近出現(xiàn)套峰的PCR產(chǎn)物進(jìn)行TA克隆,隨機(jī)選取單克隆測(cè)序,測(cè)序結(jié)果與野生型序列比對(duì)分析確定突變基因型。

    1.4 潛在脫靶位點(diǎn)分析

    使用The Rice Annotation Project(http://rapdb.dna. affrc.go.jp/)在線BLAST工具,將靶序列在水稻全基因組中進(jìn)行比對(duì),選取與靶序列同源性較高、錯(cuò)配在4 bp以內(nèi)且具有PAM序列NGG的位點(diǎn)作為潛在脫靶序列,通過酶切法評(píng)估其脫靶效應(yīng)。

    2 結(jié)果

    2.1靶點(diǎn)設(shè)計(jì)和表達(dá)載體構(gòu)建

    根據(jù)CRISPR/Cas9技術(shù)原理,在外顯子區(qū)域靠近起始密碼ATG區(qū)域選擇一段19 bp的序列(5′-GACTTCTTCTCCGAGCTCC-3′)作為靶位點(diǎn)。該序列3′端PAM序列為AGG,Cas9蛋白將在AGG序列上游第3—4 bp處剪切DNA形成平滑末端,該剪切處有一段限制性內(nèi)切酶Ⅰ識(shí)別序列(5′-GAGCT▼C-3′),選擇能被限制性內(nèi)切酶識(shí)別并切割的位點(diǎn)進(jìn)行突變,有利于后續(xù)對(duì)轉(zhuǎn)基因陽性植株是否發(fā)生突變進(jìn)行篩選。CRISPR/Cas9表達(dá)載體的靶序列由啟動(dòng)子驅(qū)動(dòng),編碼Cas9蛋白基因由玉米泛素基因()啟動(dòng)子驅(qū)動(dòng)。

    2.2T0代突變體篩選鑒定

    使用引物OsPIL15-F和OsPIL15-R擴(kuò)增25株陽性轉(zhuǎn)基因單株靶點(diǎn)序列(圖1-A),野生型擴(kuò)增片段大小為673 bp,單株1、2、23和24均擴(kuò)增出片段明顯小于野生型的條帶,初步表明這4個(gè)單株可能發(fā)生較大片段缺失,其中單株1和24的PCR產(chǎn)物出現(xiàn)兩條帶,只有其中一條帶明顯小于野生型條帶,因此二者可能只有一條同源染色體有較大片段缺失;單株2和23的PCR產(chǎn)物只有一條帶且明顯小于673 bp,可初步判斷二者兩條同源染色體都發(fā)生較大片段缺失。其余21個(gè)單株不能根據(jù)條帶大小判斷是否突變,因此,利用酶切法繼續(xù)進(jìn)行鑒定。

    使用Ⅰ酶切PCR產(chǎn)物,若PCR產(chǎn)物能被完全切開,得到400和273 bp的兩條帶,表明靶位點(diǎn)未發(fā)生突變,為未突變單株;若PCR產(chǎn)物被部分切開,得到400、273和約673 bp三條帶,表明一條同源染色體突變,為雜合體;若PCR產(chǎn)物完全不能切開,表明兩條同源染色體均發(fā)生突變,則可能為2個(gè)等位基因發(fā)生不同突變的雙等位突變體或2個(gè)等位基因發(fā)生相同突變的純合體。如圖1-B所示,結(jié)果表明,單株4、5、6、7、10、12、15、16、18和25未發(fā)生突變;單株8、9、11和17的PCR產(chǎn)物被部分酶切,形成三條帶,為一條同源染色體突變的雜合體;單株1、2、3、13、14、19、20、21、22、23和24的PCR產(chǎn)物完全不能切開,則可能為雙等位突變或純合突變體。

    PCR產(chǎn)物片段大小和酶切PCR產(chǎn)物只能初步定性判斷轉(zhuǎn)基因植株靶點(diǎn)是否發(fā)生突變,而突變體靶點(diǎn)序列的改變還需測(cè)序進(jìn)行鑒定。由于雜合突變或雙等位突變的兩條同源染色體突變情況不一致,用PCR產(chǎn)物測(cè)序時(shí)會(huì)在某一位點(diǎn)(一般在突變位點(diǎn)附近)之后出現(xiàn)重疊峰,后續(xù)序列無法準(zhǔn)確得知,因此,對(duì)可能出現(xiàn)重疊峰的突變體進(jìn)行TA克隆,隨機(jī)選取單菌斑進(jìn)行測(cè)序。同時(shí)對(duì)突變基因型一致的單株進(jìn)行Southern blot分析(電子附圖1),檢測(cè)是否屬同一株系,結(jié)果表明,1和24、2和23、3和14、19和21分屬同一株系。綜合分析測(cè)序結(jié)果和Southern blot結(jié)果得到各單株突變情況(圖2),共獲得5株純合突變體(2、3、13、14和23)、6株雙等位突變體(1、19、20、21、22和24)和4株雜合突變體(8、9、11和17),包含10種不同突變基因型和11個(gè)突變株系,多為單堿基的插入或缺失,但也得到了2種56和66 bp較大片段缺失的株系。

    A:PCR產(chǎn)物酶切前;B:PCR產(chǎn)物酶切后。1—25:T0代轉(zhuǎn)基因單株;M:DL1000標(biāo)記;WT:野生型(日本晴);+RE:使用SacⅠ酶切;黑色三角:SacⅠ酶切后預(yù)期片段。下同

    藍(lán)色字母:靶位點(diǎn)序列;黃色高亮:PAM序列;紅色刪除線:缺失堿基;紅色小寫字母:插入堿基;-示缺失;+示插入;M1/M2:雙等位2種突變型

    2.3T1代突變體分析鑒定

    為驗(yàn)證T0代突變體產(chǎn)生的突變序列能否穩(wěn)定遺傳,選取部分純合突變、雙等位突變和雜合突變株系的下一代植株,對(duì)這些T1代植株突變情況進(jìn)行測(cè)序分析(表2),純合突變體(3和14)后代為純合突變單株,但株系14的30株T1代植株中有29株純合突變,發(fā)現(xiàn)1株未突變單株;3個(gè)雙等位突變體(21、22和24)后代既獲得了純合突變體,也獲得了雙等位突變體;2個(gè)雜合突變體(9和17)的T1代單株中,既有純合突變和雜合突變體,也有未突變的單株。上述結(jié)果表明,T0代突變體產(chǎn)生的突變基因型基本上能穩(wěn)定遺傳給下一代,純合株系后代為純合體,雙等位株系后代產(chǎn)生純合體和雙等位突變體,而雜合株系后代會(huì)得到未突變、純合和雜合3種類型。

    為了鑒定T1代單株是否會(huì)重新發(fā)生基因編輯,隨機(jī)選取了62株T1代為陽性(含有Cas9基因)且在T0代未發(fā)生突變的轉(zhuǎn)基因單株,使用限制性內(nèi)切酶Ⅰ進(jìn)行酶切鑒定(圖3),62個(gè)陽性轉(zhuǎn)基因單株P(guān)CR產(chǎn)物均被切開,得到400和273 bp的兩條帶,因此,判定這62株T1代轉(zhuǎn)基因植株均未發(fā)生突變。表明T0代未突變株系,T1代依然未出現(xiàn)突變。

    表2 T1代植株分析鑒定

    M1/M2示圖2中兩條同源染色體突變類型

    M1/M2 indicates the mutation types of two homologous chromosomes in Fig. 2

    62個(gè)編號(hào):T0代未突變(4、5、6、7、10、12、15、16、18和25)的T1代陽性轉(zhuǎn)基因單株;-RE:未進(jìn)行酶切。下同

    2.4突變體脫靶效應(yīng)評(píng)估

    CRISPR/Cas9系統(tǒng)基因編輯的特異性主要取決于sgRNA的識(shí)別序列,但由于高等生物基因組較大且較為復(fù)雜,sgRNA可能會(huì)與非靶點(diǎn)DNA序列發(fā)生局部錯(cuò)配,造成CRISPR/Cas9系統(tǒng)基因編輯的脫靶。因此,需要對(duì)的突變體脫靶效應(yīng)進(jìn)行評(píng)估(表3),選取水稻基因組中與sgRNA序列錯(cuò)配堿基數(shù)在4 bp以內(nèi)且具有PAM序列的5個(gè)位點(diǎn)作為潛在脫靶位點(diǎn),每個(gè)潛在脫靶位點(diǎn)均隨機(jī)選取10株T0代突變單株和10株T1代突變單株,通過酶切法進(jìn)行檢測(cè)(圖4)。Off-target1位點(diǎn)若被Ⅰ酶切可形成444和179 bp的兩條帶,Off-target2位點(diǎn)酶切后產(chǎn)生359和173 bp的兩條帶,Off-target3位點(diǎn)酶切可得到268和152 bp的兩條帶,Off-target4位點(diǎn)可被酶切為411和153 bp的兩條帶,Off-target5位點(diǎn)若被切開會(huì)得到360和254 bp兩條帶。所選取的20個(gè)突變單株的5個(gè)潛在脫靶位點(diǎn)PCR產(chǎn)物均能被Ⅰ酶成功切割,因此,這5個(gè)潛在脫靶位點(diǎn)在T0代和T1代均未產(chǎn)生脫靶效應(yīng),表明所設(shè)計(jì)的sgRNA序列具有較高的特異性。

    表3 CRISPR/Cas9系統(tǒng)潛在脫靶位點(diǎn)突變檢測(cè)

    灰色標(biāo)記為錯(cuò)配堿基,下劃線標(biāo)記為PAM序列

    mismatching bases are shown by a grey box, the PAM motif is marked by an underline

    20個(gè)編號(hào):10株T0代突變體和10株T1代突變體

    2.5功能初步分析

    選擇水稻中PIFs家族的6個(gè)基因(Os12g0610200)、(Os03g0639300)、(Os03g0782500)、(Os07g0143200)、(Os01g0286100)、(Os05g0139100)和擬南芥中PIFs家族的7個(gè)基因/(AT2G20180)、(AT1G09530)、(AT2G43010)、/(AT3G59060)、/(AT3G62090)、(AT5G61270)、(AT4G00050),根據(jù)它們的氨基酸序列進(jìn)行系統(tǒng)進(jìn)化樹分析(圖5-A),水稻中光敏色素互作因子OsPIL15和OsPIL16親緣關(guān)系最近,二者蛋白同源性幾乎達(dá)到了100%。有研究表明,超表達(dá)能使水稻籽粒減小[22],OsPIL15和OsPIL16十分相近的親緣關(guān)系表明也有可能調(diào)控水稻籽粒大小。因此,選取T1代不同基因型的三組突變體,對(duì)其表型進(jìn)行觀察,與野生型相比,突變體生育期沒有變化,株高極顯著降低,平均降低12.49%;有效分蘗數(shù)沒有明顯差異;而籽粒粒長(zhǎng)極顯著增加,其中,突變體22增加最高,增幅達(dá)5.69%,增幅最小的突變體21也達(dá)3.88%(圖5-B—圖5-E)。

    A:OsPIL15蛋白系統(tǒng)進(jìn)化樹分析;B:ospil15突變體株高變化;C:ospil15突變體有效分蘗數(shù)變化;D:ospil15突變體籽粒粒長(zhǎng)變化;E:ospil15突變體與野生型日本晴籽粒表型比對(duì)。**:在0.01水平上差異顯著;Bar=0.5 cm

    3 討論

    CRISPR/Cas9技術(shù)以其突變高效、操作簡(jiǎn)單、成本低廉等優(yōu)勢(shì)在動(dòng)植物中得到了廣泛應(yīng)用。本研究從靶位點(diǎn)選擇,表達(dá)載體構(gòu)建,突變位點(diǎn)分析和脫靶效應(yīng)評(píng)估等方面,詳細(xì)闡述了利用CRISPR/Cas9系統(tǒng)對(duì)水稻的編輯。本研究針對(duì)設(shè)計(jì)了1個(gè)靶向位點(diǎn),有研究認(rèn)為單靶點(diǎn)在T0代形成的小片段突變,使得T1代植株中sgRNA仍有可能識(shí)別靶序列造成二次編輯,后代突變具有不穩(wěn)定性[23]。本研究在純合突變體14后代中發(fā)現(xiàn)了一株未突變單株,且突變體14僅為1個(gè)單堿基插入突變,是否是由于CRISPR/ Cas9系統(tǒng)在T1代再次對(duì)靶位點(diǎn)進(jìn)行編輯形成新的堿基改變還需進(jìn)一步驗(yàn)證。本試驗(yàn)獲得突變體多為單堿基插入或缺失,但也得到了兩種分別缺失66和56 bp的較大片段缺失突變體。針對(duì)基因編碼序列設(shè)計(jì)多個(gè)靶位點(diǎn)有利于形成大片段缺失突變而獲得穩(wěn)定遺傳,但靶點(diǎn)越多產(chǎn)生脫靶的可能性越大[24],需綜合考慮目標(biāo)基因靶點(diǎn)位置及數(shù)量,盡可能保證其特異性。本試驗(yàn)采用酶切和測(cè)序相結(jié)合的方法對(duì)突變位點(diǎn)進(jìn)行分析鑒定,但由于雜合突變和雙等位突變的PCR產(chǎn)物直接測(cè)序會(huì)出現(xiàn)重疊峰而無法準(zhǔn)確讀取,需將PCR產(chǎn)物進(jìn)行TA克隆后挑選單菌斑測(cè)序,工作量較大,成本較高。為此華南農(nóng)業(yè)大學(xué)劉耀光教授課題組基于DSD(degenerate sequence decoding)法[25]開發(fā)了在線工具DSDecode(http://dsdecode.scgene.com/),可對(duì)各突變類型測(cè)序文件自動(dòng)解碼[26],但限于不同公司和不同PCR產(chǎn)物測(cè)序質(zhì)量,部分測(cè)序結(jié)果不能很好地解碼,因此,此方法可為測(cè)序分析提供參考,具體仍需驗(yàn)證。

    CRISPR/Cas9系統(tǒng)成功對(duì)植物基因組編輯后,所獲得突變基因序列能否穩(wěn)定遺傳給下一代決定了該突變體是否具有應(yīng)用前景。由于擬南芥多采用蘸花法進(jìn)行轉(zhuǎn)基因,CRISPR/Cas9系統(tǒng)誘導(dǎo)的突變?cè)赥1代中主要發(fā)生在體細(xì)胞,從而形成嵌合體,較難穩(wěn)定遺傳給下一代,而發(fā)生在生殖細(xì)胞中的突變可穩(wěn)定遺傳[27]。水稻多采用農(nóng)桿菌介導(dǎo)愈傷組織,轉(zhuǎn)基因苗多來源于單個(gè)細(xì)胞或細(xì)胞團(tuán),因此本試驗(yàn)在T0代即獲得了純合突變體,對(duì)純合突變、雙等位突變和雜合突變的T1代植株分析顯示這些突變絕大多數(shù)均能穩(wěn)定遺傳,有研究證明這些穩(wěn)定遺傳突變體自交后代分離均符合孟德爾遺傳定律[28-29]。有報(bào)道認(rèn)為CRISPR/Cas9系統(tǒng)在后代植株當(dāng)中對(duì)野生型依然可行使功能[30],T0代轉(zhuǎn)基因植株靶位點(diǎn)表現(xiàn)為野生型的植株后代仍有可能發(fā)生一定比例的突變[27]。本研究對(duì)T0代未突變植株的T1代(陽性)靶位點(diǎn)突變情況進(jìn)行鑒定,未獲得新的突變單株。Zhang等[31]認(rèn)為這可能與CRISPR/Cas9系統(tǒng)插入到基因組的位置有一定關(guān)系。本研究?jī)H對(duì)部分T1代植株進(jìn)行了鑒定,CRISPR/Cas9系統(tǒng)是否會(huì)在其他后繼世代中重新發(fā)揮作用還需作進(jìn)一步驗(yàn)證。

    CRISPR/Cas9系統(tǒng)在植物基因編輯中的特異性對(duì)于其安全高效的利用具有重要意義,因此有必要對(duì)可能出現(xiàn)的脫靶效應(yīng)進(jìn)行全面準(zhǔn)確地評(píng)估。有研究表明,針對(duì)人源CLTA基因的4個(gè)靶位點(diǎn)均存在脫靶現(xiàn)象,脫靶效率最高可達(dá)84%[32];針對(duì)鼠源Crygc基因的2個(gè)靶點(diǎn)存在40%的脫靶率[33]。而在植物中CRISPR/Cas9系統(tǒng)脫靶效率較低,朱健康教授實(shí)驗(yàn)室通過對(duì)擬南芥全基因組測(cè)序和潛在脫靶位點(diǎn)分析未檢測(cè)到脫靶效應(yīng)[27],對(duì)水稻中多個(gè)與靶位點(diǎn)高度同源序列進(jìn)行檢測(cè),僅在一個(gè)堿基不匹配處發(fā)現(xiàn)脫靶效應(yīng)[31],杜彥修等[34]對(duì)水稻基因兩個(gè)潛在脫靶位點(diǎn)進(jìn)行檢測(cè)并未發(fā)現(xiàn)脫靶。本研究選擇了PAM序列5′端同源性較高的5個(gè)位點(diǎn)作為潛在脫靶位點(diǎn),經(jīng)檢測(cè)這5個(gè)位點(diǎn)均未發(fā)生脫靶效應(yīng),因此設(shè)計(jì)的sgRNA序列具有較高特異性。有研究表明,sgRNA的特異性主要由靠近PAM序列的10—12 bp堿基所決定,而遠(yuǎn)離PAM區(qū)的8—10 bp堿基對(duì)sgRNA的特異性識(shí)別影響不大[3, 35]。本試驗(yàn)只是對(duì)最有可能出現(xiàn)脫靶效應(yīng)的5個(gè)位點(diǎn)進(jìn)行了檢測(cè),不排除其他位點(diǎn)出現(xiàn)脫靶的可能性,更有效的方法則是對(duì)獲得的突變體進(jìn)行全基因組深度測(cè)序分析,但其成本較高且費(fèi)時(shí)費(fèi)力。為此科學(xué)家開發(fā)了新的基于Cas9核酸酶消化基因組(Digenome- seq)技術(shù)和基于整合酶缺陷型慢病毒載體(IDLV)技術(shù)來精確檢測(cè)CRISPR/Cas9系統(tǒng)脫靶效應(yīng)[36-37]。通過優(yōu)化sgRNA序列(長(zhǎng)度和GC含量等)、應(yīng)用雙Cas9切口酶[38]和應(yīng)用fCas9復(fù)合體[39]等方法可降低CRISPR/Cas9系統(tǒng)脫靶效應(yīng)。

    本試驗(yàn)利用CRISPR/Cas9技術(shù)創(chuàng)制了突變體,并對(duì)其T1代表型進(jìn)行了初步觀察,發(fā)現(xiàn)突變體株高極顯著降低,籽粒粒長(zhǎng)極顯著增大,與OsPIL15蛋白系統(tǒng)進(jìn)化分析結(jié)果是一致的。后續(xù)仍需進(jìn)行多代鑒定,獲得穩(wěn)定純合的突變體,同時(shí)對(duì)不同突變類型的株系進(jìn)行觀察,以確定不同突變類型造成的氨基酸改變而產(chǎn)生的表型是否一致。

    4 結(jié)論

    依據(jù)CRISPR/Cas9系統(tǒng)原理創(chuàng)制突變體,成功獲得5株純合突變體、6株雙等位突變體和4株雜合突變體。對(duì)部分突變體后代突變序列進(jìn)行分析,發(fā)現(xiàn)絕大多數(shù)突變均能穩(wěn)定遺傳。62株T0代表現(xiàn)為未突變的T1代陽性植株經(jīng)酶切鑒定均未發(fā)現(xiàn)新的突變,表明CRISPR/Cas9系統(tǒng)在T1代轉(zhuǎn)基因陽性植株中未重新發(fā)揮基因編輯作用。5個(gè)與靶序列同源性較高且錯(cuò)配在4 bp以內(nèi)的潛在脫靶位點(diǎn)經(jīng)評(píng)估均未出現(xiàn)脫靶效應(yīng),表明所選取sgRNA有較高的特異性。對(duì)獲得的突變體表型進(jìn)行初步觀察發(fā)現(xiàn),與野生型相比T1代突變體株高極顯著降低,分蘗數(shù)未發(fā)生顯著變化,籽粒粒長(zhǎng)極顯著增大。

    References

    [1] Carroll D, Morton J J, Beumer K J, Segal D J. Design, construction and in vitro testing of zinc finger nucleases., 2006, 1(3): 1329-1341.

    [2] Li T, Liu B, Spalding M H, Weeks D P, Yang B. High- efficiency TALEN-based gene editing produces disease-resistant rice., 2012, 30(5): 390-392.

    [3] Cong L, Ran F A, Cox D, Lin S, Barretto R, Habib N, Hsu P D, Wu X, Jiang W, Marraffini L A, Zhang F. Multiplex genome engineering using CRISPR/Cas systems., 2013, 339(6121): 819-823.

    [4] Mali P, Yang L, Esvelt K M, Aach J, Guell M, Dicarlo J E, Norville J E, Church G M. RNA-guided human genome engineering via Cas9., 2013, 339(6121): 823-826.

    [5] Wiedenheft B, Sternberg S H, Doudna J A. RNA-guided genetic silencing systems in bacteria and archaea., 2012, 482(7385): 331-338.

    [6] Doudna J A, Charpentier E. The new frontier of genome engineering with CRISPR-Cas9., 2014, 346(6213): 1258096.

    [7] Belhaj K, Chaparro-Garcia A, Kamoun S, Patron N J, Nekrasov V. Editing plant genomes with CRISPR/Cas9., 2015, 32: 76-84.

    [8] Khanna R, Huq E, Kikis E A, Al-Sady B, Lanzatella C, Quail P H. A novel molecular recognition motif necessary for targeting photoactivated phytochrome signaling to specific basic helix-loop-helix transcription factors., 2004, 16(11): 3033-3044.

    [9] Shen H, Zhu L, Castillon A, Majee M, Downie B, Hu q E. Light-induced phosphorylation and degradation of the negative regulator PHYTOCHROME-INTERACTING FACTOR1 fromdepend upon its direct physical interactions with photoactivated phytochromes., 2008, 20(6): 1586-1602.

    [10] Zhang Y, Mayba O, Pfeiffer A, Shi H, Tepperman J M, Speed T P, Quail P H. A quartet of PIF bHLH factors provides a transcriptionally centered signaling hub that regulates seedling morphogenesis through differential expression-patterning of shared target genes in., 2013, 9(1): e1003244.

    [11] Oh E, Yamaguchi S, Hu J, Yusuke J, Jung B, Paik I, Lee H, Sun T, Kamiya Y, Choi G. PIL5, a phytochrome-interacting bHLH protein, regulates gibberellin responsiveness by binding directly to theandpromoters inseeds., 2007, 19(4): 1192-1208.

    [12] Leivar P, Monte E, Al-Sady B, Carle C, Storer A, Alonso J M, Ecker J R, Quail P H. Thephytochrome-interacting factor PIF7, together with PIF3 and PIF4, regulates responses to prolonged red light by modulating phyB levels., 2008, 20(2): 337-352.

    [13] Lorrain S, Allen T, Duek P D, Whitelam G C, Fankhauser C. Phytochrome-mediated inhibition of shade avoidance involves degradation of growth-promoting bHLH transcription factors., 2008, 53(2): 312-323.

    [14] Shin J, Anwer M U, Davis S J. Phytochrome-interacting factors (PIFs) as bridges between environmental signals and the circadian clock: diurnal regulation of growth and development., 2013, 6(3): 592-595.

    [15] Nakamura Y, Kato T, Yamashino T, Murakami M, Mizuno T. Characterization of a set of phytochrome-interacting factor-like bHLH proteins in., 2007, 71(5): 1183-1191.

    [16] 梁衛(wèi)紅, 李莉, 彭威風(fēng), 劉婧, 楊慧, 謝先芝. 水稻光敏色素相互作用因子基因的表達(dá)模式及其在光信號(hào)傳導(dǎo)中作用的初步分析. 中國(guó)水稻科學(xué), 2012, 26(3): 255-260.

    Liang W H, Li L, PENG W F, LIU J, YANG H, XIE X Z. Expression patterns of, a phytochrome interacting factor in rice, and preliminary analysis of its roles in light signal transduction., 2012, 26(3): 255-260. (in Chinese)

    [17] Todaka D, Nakashima K, Maruyama K, Kidokoro S, Osakabe Y, Ito Y, Matsukura S, Fujita Y, Yoshiwara K, Ohme-Takagi M. Rice phytochrome-interacting factor-like protein OsPIL1 functions as a key regulator of internode elongation and induces a morphological response to drought stress., 2012, 109(39): 15947-15952.

    [18] Zhou J, Liu Q, Zhang F, Wang Y, Zhang S, Cheng H, Yan L, Li L, Chen F, Xie X. Overexpression of, a phytochrome‐interacting factor-like protein gene, represses etiolated seedling growth in rice., 2014, 56(4): 373-387.

    [19] Heang D, Sassa H. Antagonistic actions of HLH/bHLH proteins are involved in grain length and weight in rice., 2012, 7(2): e31325.

    [20] 杜彥修, 季新, 彭廷, 孫紅正, 張靜, 李俊周, 趙全志. 水稻灌漿期籽粒中光信號(hào)相關(guān)基因鑒定與表達(dá)分析. 分子植物育種, 2016, 14(7): 1637-1647.

    DU Y X, JI X, PENG T, SUN H Z, ZHANG J, LI J Z, ZHAO Q Z. Identification and expressing analysis of light signal genes involved in endosperm development during rice grain filling., 2016, 14(7): 1637-1647. (in Chinese)

    [21] Xing H, Dong L, Wang Z, Zhang H, Han C, Liu B, Wang X, Chen Q. A CRISPR/Cas9 toolkit for multiplex genome editing in plants., 2014, 14(1): 327.

    [22] Heang D, Sassa H. Overexpression of a basic helix-loop-helix gene(APG) decreases grain length of rice., 2012, 29(1): 65-69.

    [23] Xu R, Li H, Qin R, Li J, Qiu C, Yang Y, Ma H, Li L, Wei P, Yang J. Generation of inheritable and “transgene clean” targeted genome-modified rice in later generations using the CRISPR/Cas9 system., 2015, 5: 11491.

    [24] Endo M, Mikami M, Toki S. Multigene knockout utilizing off-target mutations of the CRISPR/Cas9 system in rice., 2015, 56(1): 41-47.

    [25] Ma X, Chen L, Zhu Q, Chen Y, Liu Y. Rapid decoding of sequence-specific nuclease-induced heterozygous and biallelic mutations by direct sequencing of PCR products., 2015, 8(8): 1285-1287.

    [26] Liu W, Xie X, Ma X, Li J, Chen J, Liu Y. DSDecode: A web- based tool for decoding of sequencing chromatograms for genotyping of targeted mutations., 2015, 8(9): 1431-1433.

    [27] Feng Z, Mao Y, Xu N, Zhang B, Wei P, Yang D, Wang Z, Zhang Z, Zheng R, Yang L, Zeng L, Liu X, Zhu J. Multigeneration analysis reveals the inheritance, specificity, and patterns of CRISPR/Cas-induced gene modifications in., 2014, 111(12): 4632-4637.

    [28] Ishizaki T. CRISPR/Cas9 in rice can induce new mutations in later generations, leading to chimerism and unpredicted segregation of the targeted mutation., 2016, 36(12): 165.

    [29] Mao Y, Zhang Z, Feng Z, Wei P, Zhang H, Botella J R, Zhu J. Development of germ-line-specific CRISPR-Cas9 systems to improve the production of heritable gene modifications in., 2016, 14(2): 519-532.

    [30] Ma X, Zhang Q, Zhu Q, Liu W, Chen Y, Qiu R, Wang B, Yang Z, Li H, Lin Y, Xie Y, Shen R, Chen S, Wang Z, Chen Y, Guo J, Chen L, Zhao X, Dong Z, Liu Y. A robust CRISPR/Cas9 system for convenient, high-efficiency multiplex genome editing in monocot and dicot plants., 2015, 8(8): 1274-1284.

    [31] Zhang H, Zhang J, Wei P, Zhang B, Gou F, Feng Z, Mao Y, Yang L, Zhang H, Xu N, Zhu J. The CRISPR/Cas9 system produces specific and homozygous targeted gene editing in rice in one generation., 2014, 12(6): 797-807.

    [32] Pattanayak V, Lin S, Guilinger J P, Ma E, Doudna J A, Liu D R. High-throughput profiling of off-target DNA cleavage reveals RNA-programmed Cas9 nuclease specificity., 2013, 31(9): 839-843.

    [33] Wu Y, Liang D, Wang Y, Bai M, Tang W, Bao S, Yan Z, Li D, Li J. Correction of a genetic disease in mouse via use of CRISPR-Cas9., 2013, 13(6): 659-662.

    [34] 杜彥修, 季新, 陳會(huì)杰, 彭廷, 張靜, 李俊周, 孫紅正, 趙全志. 基于CRISPR/Cas9系統(tǒng)的基因編輯及其脫靶效應(yīng)分析. 中國(guó)水稻科學(xué), 2016, 30(6): 577-586.

    DU Y X, JI X, CHEN H J, PENG T, ZHANG J, LI J Z, SUN H Z, ZHAO Q Z. CRISPR/Cas9 system-based editing ofgene and its off-target effect analysis., 2016, 30(6): 577-586. (in Chinese)

    [35] Hsu P D, Scott D A, Weinstein J A, Ran F A, Konermann S, Agarwala V, Li Y, Fine E J, Wu X, Shalem O, Cradick T J, Marraffini L A, Bao G, Zhang F. DNA targeting specificity of RNA-guided Cas9 nucleases., 2013, 31(9): 827-832.

    [36] Kim D, Bae S, Park J, Kim E, Kim S, Yu H R, Hwang J, Kim J, Kim J. Digenome-seq: genome-wide profiling of CRISPR-Cas9 off-target effects in human cells., 2015, 12(3): 237-243.

    [37] Wang X, Wang Y, Wu X, Wang J, Wang Y, Qiu Z, Chang T, Huang H, Lin R, Yee J. Unbiased detection of off-target cleavage by CRISPR-Cas9 and TALENs using integrase-defective lentiviral vectors., 2015, 33(2): 175-178.

    [38] Ran F A, Hsu P D, Lin C, Gootenberg J S, Konermann S, Trevino A E, Scott D A, Inoue A, Matoba S, Zhang Y, Zhang F. Double nicking by RNA-guided CRISPR Cas9 for enhanced genome editing specificity., 2013, 154(6): 1380-1389.

    [39] Guilinger J P, Thompson D B, Liu D R. Fusion of catalytically inactive Cas9 to FokI nuclease improves the specificity of genome modification., 2014, 32(6): 577-582.

    (責(zé)任編輯 李莉)

    M:DNA Molecular-Weight Marker;B:Basta探針;1、24、2、23、3、14、19、21:突變單株

    M: DNA Molecular-Weight Marker; B: Basta probe; 1, 24, 2, 23, 3, 14, 19, 21: mutant plants

    附圖1 Southern blot分析突變單株的T-DNA插入

    Supplemental Fig. 1 Analysis of T -DNA region by southern blot

    CRISPR/Cas9 System-Based Editing of Phytochrome-Interacting Factor

    JI Xin, LI Fei, YAN Yun, SUN HongZheng, ZHANG Jing, LI JunZhou, PENG Ting, DU YanXiu, ZHAO QuanZhi

    (Agronomy College of Henan Agricultural University/Collaborative Innovation Center of Henan Grain Crops/Henan Key Laboratory of Rice Biology, Zhengzhou 450002)

    【Objective】 As an important environmental signal, light can regulate gene expression, affect the activity of enzymes and plant morphogenesis. Phytochrome interacting factors play an important role in the signal transduction of light. Therefore, constructing the expression vector of CRISPR/Cas9 containing rice phytochrome-interacting factorsand creating themutants have important significance to exploit functional genes and enrich light signal regulation of the molecular mechanism of rice. 【Method】 According to the principle of CRISPR/Cas9 technology, the sgRNA was designed. To exclude non-specific target sites, the sgRNA was analyzed by sequence alignment in the rice genome database. At the same time, the target sequence contained the common restriction site to identify mutants. The oligonucleotides of sgRNA were chemically synthesized and inserted into linearized plasmid pBUN411 to construct the expression vector. Transgenic rice plants harboring sgRNA:Cas9 were obtained by-mediated stable transformation and the positive transgenic plants were screened by herbicide resistance. The PCR products of T0transgenic plants were digested by restriction enzyme to judge whether they were mutants. Then, the mutated genotypes of these mutants were analyzed by DNA sequencing. After searching rice genome using sgRNA sequence, five highly identical sites with less 4 mismatching bases were selected to assess off-target efficiency and specificity of sgRNA.【Result】The recombinant vector succeeded in oriented editing of. The restriction enzyme analysis results indicated that 15 mutants from the 25 randomly selected T0transgenic lines were obtained. They included 5 homozygous mutants, 6 biallelic mutants and 4 heterozygous mutants, and a total of 10 different genotypes and 11 mutant lines. The mutant types were mainly insertions or deletions of single base, besides, two mutant types of large deletions with 56 and 66 bp were obtained. Analysis of T1transgenic plants from some T0mutants indicated that the genotypes in T0mutants could descend stably into the next generation. The progeny of homozygous mutants in T0were homozygous mutants. However, only one wild type with no mutation was detected in a homozygous mutant progenies of line 14. Two homozygous mutations and one biallelic mutations were obtained in the progeny of biallelic T0lines. Three mutation types including homozygous mutants, heterozygous mutants and wild type with no mutation were detected in the progenies of heterozygous T0mutants. Restriction enzyme analysis was used to detect the engineered target site of T1positive transgenic plants which had no mutation in T0. No mutation in 62 T1plants was detected. The results showed that CRISPR/Cas9 system played no role in gene editing in T1positive transgenic plants. After searching the rice genome using the target sequence with PAM, five highly identical sites were found. However, any mutations at these sites in T0and T1generations were not observed by restriction enzyme, which indicated the sgRNA was highly specific. Three groups of different genotypes were selected representativemutants in T1generation to observe phenotypes, the investigation result showed that the growth stage and tiller number were not obviously changed, the plant height decreased significantly (<0.01), and the grain length was significantly increased compared with the wild type, up 5.69%. 【Conclusion】CRISPR/Cas9 system succeeded in oriented editing of, themutants with 10 different genotypes were obtained and observed that the plant height decreased significantly and the grain length increased significantly compared with those of wild type.

    rice; CRISPR/Cas9; gene editing;; off-target efficiency

    2017-01-22;接受日期:2017-04-05

    河南省重大科技專項(xiàng)(141100110600)、河南省高等學(xué)校重點(diǎn)科研項(xiàng)目(16A210026)

    季新,E-mail:jixin1@126.com。李飛,E-mail:LFight123@163.com。季新和李飛為同等貢獻(xiàn)作者。通信作者杜彥修,E-mail:dumatt@163.com。通信作者趙全志,E-mail:qzzhaoh@126.com

    猜你喜歡
    靶位等位突變體
    利用CRISPR/Cas9 技術(shù)靶向編輯青花菜BoZDS
    芥菜種子顏色調(diào)控基因TT8的等位變異及其地理分布分析
    ·術(shù)語解析·
    外引大麥農(nóng)藝性狀SSR關(guān)聯(lián)位點(diǎn)及等位變異表型效應(yīng)分析
    花時(shí)間在餐廳門口等位值嗎
    CLIC1及其點(diǎn)突變體與Sedlin蛋白的共定位研究
    擬南芥干旱敏感突變體篩選及其干旱脅迫響應(yīng)機(jī)制探究
    Survivin D53A突變體對(duì)宮頸癌細(xì)胞增殖和凋亡的影響
    CT引導(dǎo)下靶位注射膠原酶治療腰椎間盤突出癥36例
    耐喹諾酮銅綠假單胞菌藥物作用靶位改變的研究
    亚洲国产看品久久| 欧美黄色片欧美黄色片| 母亲3免费完整高清在线观看| 欧美成狂野欧美在线观看| 国产精品亚洲av一区麻豆| 国内精品久久久久精免费| www日本在线高清视频| 身体一侧抽搐| 男女做爰动态图高潮gif福利片 | 欧美日本亚洲视频在线播放| 国产精品秋霞免费鲁丝片| 国产色视频综合| 麻豆av在线久日| 免费女性裸体啪啪无遮挡网站| 女人高潮潮喷娇喘18禁视频| 巨乳人妻的诱惑在线观看| 久久国产精品男人的天堂亚洲| 精品久久久久久,| 97碰自拍视频| 人人妻,人人澡人人爽秒播| 国产精品综合久久久久久久免费 | 国产成人系列免费观看| av在线播放免费不卡| 夜夜夜夜夜久久久久| 精品国产乱码久久久久久男人| 欧美黑人欧美精品刺激| 国产精品九九99| 精品国产乱码久久久久久男人| 51午夜福利影视在线观看| 国产精品爽爽va在线观看网站 | 69av精品久久久久久| e午夜精品久久久久久久| 啦啦啦韩国在线观看视频| 一级a爱视频在线免费观看| 在线视频色国产色| 欧美激情高清一区二区三区| 国产成人欧美在线观看| 黑人操中国人逼视频| 99国产精品一区二区三区| 精品卡一卡二卡四卡免费| 高潮久久久久久久久久久不卡| 国产精品免费视频内射| 色综合亚洲欧美另类图片| 69av精品久久久久久| 成熟少妇高潮喷水视频| 少妇 在线观看| 一a级毛片在线观看| 搡老妇女老女人老熟妇| 亚洲一区高清亚洲精品| 好看av亚洲va欧美ⅴa在| 欧美日韩黄片免| 男女午夜视频在线观看| 在线永久观看黄色视频| 亚洲片人在线观看| 亚洲熟妇中文字幕五十中出| 丝袜美腿诱惑在线| 母亲3免费完整高清在线观看| 国产色视频综合| 夜夜夜夜夜久久久久| 成人18禁在线播放| 精品久久蜜臀av无| 91大片在线观看| 啦啦啦观看免费观看视频高清 | av有码第一页| 日本免费a在线| 51午夜福利影视在线观看| 国产成人精品久久二区二区91| 大香蕉久久成人网| x7x7x7水蜜桃| 午夜久久久久精精品| 国产精品1区2区在线观看.| av在线播放免费不卡| 亚洲黑人精品在线| 国产日韩一区二区三区精品不卡| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩亚洲综合一区二区三区_| 国内精品久久久久久久电影| 久久精品亚洲精品国产色婷小说| 亚洲欧美日韩无卡精品| 亚洲成av片中文字幕在线观看| 99国产精品免费福利视频| 操美女的视频在线观看| 欧美中文综合在线视频| 桃色一区二区三区在线观看| 欧美一级毛片孕妇| 91麻豆av在线| 黄色丝袜av网址大全| 色尼玛亚洲综合影院| 成人欧美大片| 国产精品乱码一区二三区的特点 | 制服丝袜大香蕉在线| av片东京热男人的天堂| 成人国产综合亚洲| 午夜免费观看网址| 精品一区二区三区四区五区乱码| 99国产精品一区二区三区| 日韩 欧美 亚洲 中文字幕| 男女下面插进去视频免费观看| 色综合欧美亚洲国产小说| 日韩精品中文字幕看吧| 黄片小视频在线播放| 在线观看免费日韩欧美大片| 久久国产亚洲av麻豆专区| ponron亚洲| 琪琪午夜伦伦电影理论片6080| 中文字幕人妻熟女乱码| 中文字幕精品免费在线观看视频| av在线播放免费不卡| 久久久久久久久久久久大奶| 国产片内射在线| 久久国产精品影院| 人人妻人人爽人人添夜夜欢视频| 99riav亚洲国产免费| 日韩大尺度精品在线看网址 | 国产欧美日韩综合在线一区二区| 丝袜美足系列| 亚洲国产毛片av蜜桃av| 母亲3免费完整高清在线观看| av天堂久久9| 国产成人欧美| 嫩草影院精品99| 久久久久久亚洲精品国产蜜桃av| 亚洲最大成人中文| 国产高清激情床上av| 18美女黄网站色大片免费观看| 国内久久婷婷六月综合欲色啪| 免费无遮挡裸体视频| 两个人免费观看高清视频| 99riav亚洲国产免费| 欧美最黄视频在线播放免费| 精品久久久久久,| 国产成人精品无人区| 色av中文字幕| 亚洲 欧美 日韩 在线 免费| 久久久久久久午夜电影| 国产亚洲欧美精品永久| 国产午夜福利久久久久久| 久久久久久久午夜电影| www国产在线视频色| 久久久久久久久中文| 日日爽夜夜爽网站| 欧美日本中文国产一区发布| 欧美日韩中文字幕国产精品一区二区三区 | 精品国产超薄肉色丝袜足j| 亚洲全国av大片| 国产欧美日韩一区二区三| 久久久国产精品麻豆| 欧美丝袜亚洲另类 | 国产真人三级小视频在线观看| 国产成人精品无人区| 国产日韩一区二区三区精品不卡| 韩国精品一区二区三区| 夜夜爽天天搞| 看免费av毛片| 欧美日本亚洲视频在线播放| 久久久久久久久久久久大奶| 我的亚洲天堂| 国产不卡一卡二| 中国美女看黄片| 可以在线观看毛片的网站| 日本五十路高清| 国产精品精品国产色婷婷| 美女扒开内裤让男人捅视频| 999久久久精品免费观看国产| 女性生殖器流出的白浆| 首页视频小说图片口味搜索| 国产精品亚洲一级av第二区| 亚洲性夜色夜夜综合| 亚洲av电影不卡..在线观看| 激情视频va一区二区三区| 女性生殖器流出的白浆| 精品久久久久久成人av| 99久久久亚洲精品蜜臀av| 国产精品爽爽va在线观看网站 | 高潮久久久久久久久久久不卡| 亚洲av成人不卡在线观看播放网| 免费搜索国产男女视频| 亚洲中文日韩欧美视频| 国产主播在线观看一区二区| 久99久视频精品免费| 国产熟女xx| 免费少妇av软件| 99riav亚洲国产免费| 精品卡一卡二卡四卡免费| 亚洲国产看品久久| 国产成人av教育| 亚洲自偷自拍图片 自拍| 国产精品久久久久久精品电影 | 国产aⅴ精品一区二区三区波| 亚洲一区二区三区色噜噜| 美女高潮到喷水免费观看| 亚洲免费av在线视频| 亚洲国产精品合色在线| 亚洲久久久国产精品| 亚洲成a人片在线一区二区| 777久久人妻少妇嫩草av网站| 久久中文看片网| 日韩有码中文字幕| 成人特级黄色片久久久久久久| 成人免费观看视频高清| 精品第一国产精品| 久久人人精品亚洲av| 这个男人来自地球电影免费观看| 午夜福利影视在线免费观看| 十八禁网站免费在线| 性少妇av在线| 动漫黄色视频在线观看| 嫩草影院精品99| www.精华液| 激情视频va一区二区三区| 日韩欧美免费精品| 亚洲精品久久国产高清桃花| av在线播放免费不卡| 成人亚洲精品一区在线观看| 后天国语完整版免费观看| 成熟少妇高潮喷水视频| 国产av一区在线观看免费| 9191精品国产免费久久| 亚洲全国av大片| 欧美午夜高清在线| 电影成人av| 国内精品久久久久久久电影| 欧美乱码精品一区二区三区| 国产欧美日韩一区二区三区在线| 国产一区二区激情短视频| 99在线视频只有这里精品首页| 欧美久久黑人一区二区| 精品不卡国产一区二区三区| 黄色女人牲交| 亚洲精品中文字幕一二三四区| 亚洲欧美精品综合久久99| 国产精品亚洲美女久久久| 国产精品久久久人人做人人爽| 免费在线观看亚洲国产| 女同久久另类99精品国产91| 精品免费久久久久久久清纯| 91在线观看av| 欧美激情 高清一区二区三区| 亚洲人成伊人成综合网2020| 国产一区二区三区视频了| 性色av乱码一区二区三区2| 午夜久久久久精精品| 无遮挡黄片免费观看| 两性夫妻黄色片| 欧美乱色亚洲激情| 亚洲人成77777在线视频| 国产日韩一区二区三区精品不卡| 国产不卡一卡二| 美女 人体艺术 gogo| 麻豆成人av在线观看| 亚洲av电影在线进入| 精品第一国产精品| 国产精品久久久久久精品电影 | 国产精品一区二区在线不卡| 亚洲精品久久成人aⅴ小说| x7x7x7水蜜桃| 91字幕亚洲| 一本久久中文字幕| 老熟妇乱子伦视频在线观看| 91成年电影在线观看| 中文字幕最新亚洲高清| 精品国产一区二区三区四区第35| 操美女的视频在线观看| 日韩精品青青久久久久久| www.www免费av| 精品久久久精品久久久| 国产亚洲欧美在线一区二区| 脱女人内裤的视频| 黄片小视频在线播放| 成熟少妇高潮喷水视频| 久久久久亚洲av毛片大全| 亚洲色图 男人天堂 中文字幕| 日韩精品中文字幕看吧| 免费观看人在逋| 色尼玛亚洲综合影院| 日韩精品中文字幕看吧| 制服人妻中文乱码| www日本在线高清视频| 波多野结衣高清无吗| 日韩三级视频一区二区三区| 国产日韩一区二区三区精品不卡| 国产免费av片在线观看野外av| 欧美中文日本在线观看视频| 久久人人精品亚洲av| 操出白浆在线播放| 99久久久亚洲精品蜜臀av| 日韩大码丰满熟妇| 777久久人妻少妇嫩草av网站| 久久人妻福利社区极品人妻图片| 9191精品国产免费久久| 露出奶头的视频| 看免费av毛片| 亚洲av成人av| 欧美久久黑人一区二区| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 中文字幕人成人乱码亚洲影| 首页视频小说图片口味搜索| 午夜免费成人在线视频| 国产成人av教育| 999精品在线视频| a级毛片在线看网站| 久久久久久久久中文| 99在线人妻在线中文字幕| 国产aⅴ精品一区二区三区波| 日韩精品免费视频一区二区三区| 亚洲第一欧美日韩一区二区三区| 91字幕亚洲| 一级,二级,三级黄色视频| 欧美精品啪啪一区二区三区| 亚洲七黄色美女视频| 欧美乱码精品一区二区三区| 成在线人永久免费视频| 韩国av一区二区三区四区| 国产亚洲欧美在线一区二区| 我的亚洲天堂| 国产欧美日韩精品亚洲av| 最近最新中文字幕大全电影3 | 在线视频色国产色| 老司机午夜福利在线观看视频| 日本五十路高清| 日本一区二区免费在线视频| 欧美+亚洲+日韩+国产| 欧美乱妇无乱码| 一级a爱视频在线免费观看| 男女下面插进去视频免费观看| 欧美乱色亚洲激情| 精品国产乱码久久久久久男人| 久久国产亚洲av麻豆专区| 免费高清视频大片| 夜夜躁狠狠躁天天躁| 他把我摸到了高潮在线观看| 亚洲国产精品久久男人天堂| 色综合亚洲欧美另类图片| www.精华液| 男人舔女人下体高潮全视频| 欧美中文日本在线观看视频| 欧美在线黄色| 精品日产1卡2卡| 午夜影院日韩av| 热99re8久久精品国产| 午夜福利高清视频| 无人区码免费观看不卡| 嫩草影视91久久| 国产极品粉嫩免费观看在线| 亚洲专区中文字幕在线| 国产精品美女特级片免费视频播放器 | 一个人免费在线观看的高清视频| 久久 成人 亚洲| 亚洲熟妇中文字幕五十中出| 日日爽夜夜爽网站| 亚洲无线在线观看| 亚洲午夜理论影院| 大陆偷拍与自拍| 成人精品一区二区免费| 国产成人精品在线电影| 午夜日韩欧美国产| 亚洲性夜色夜夜综合| 国产亚洲精品综合一区在线观看 | 亚洲,欧美精品.| 国产aⅴ精品一区二区三区波| 精品国产美女av久久久久小说| 免费久久久久久久精品成人欧美视频| 免费少妇av软件| 大码成人一级视频| 久久亚洲精品不卡| 欧美大码av| 亚洲第一青青草原| 久久人人爽av亚洲精品天堂| 国产91精品成人一区二区三区| 中文字幕人妻丝袜一区二区| 9色porny在线观看| 变态另类丝袜制服| 色播在线永久视频| 亚洲九九香蕉| 岛国在线观看网站| 视频在线观看一区二区三区| 亚洲熟女毛片儿| 国产亚洲av嫩草精品影院| 国产亚洲精品久久久久久毛片| 亚洲人成电影观看| 精品一区二区三区视频在线观看免费| 亚洲色图 男人天堂 中文字幕| 久久婷婷成人综合色麻豆| 欧美乱妇无乱码| 不卡一级毛片| 亚洲人成电影免费在线| 亚洲 国产 在线| 国产国语露脸激情在线看| 日韩精品青青久久久久久| 日韩一卡2卡3卡4卡2021年| 男女午夜视频在线观看| 老司机福利观看| 在线观看免费视频网站a站| 在线av久久热| ponron亚洲| 这个男人来自地球电影免费观看| 国产亚洲欧美精品永久| 精品久久久久久成人av| 搡老熟女国产l中国老女人| 无遮挡黄片免费观看| 啪啪无遮挡十八禁网站| 国产精品一区二区免费欧美| 99re在线观看精品视频| 免费在线观看黄色视频的| 手机成人av网站| 亚洲狠狠婷婷综合久久图片| 9色porny在线观看| 国产又色又爽无遮挡免费看| 国产极品粉嫩免费观看在线| 在线观看免费视频日本深夜| 欧美激情高清一区二区三区| 成人欧美大片| 99热只有精品国产| 美女扒开内裤让男人捅视频| 精品人妻1区二区| av福利片在线| 一边摸一边做爽爽视频免费| 免费女性裸体啪啪无遮挡网站| 窝窝影院91人妻| 91麻豆精品激情在线观看国产| 欧美一级毛片孕妇| 99国产极品粉嫩在线观看| 夜夜爽天天搞| 国产精品亚洲美女久久久| 天堂动漫精品| 日本三级黄在线观看| 午夜福利18| 中文字幕色久视频| 精品一品国产午夜福利视频| 国产精品 欧美亚洲| 免费在线观看完整版高清| 久久影院123| 久久久水蜜桃国产精品网| 真人做人爱边吃奶动态| 啦啦啦韩国在线观看视频| 久久中文看片网| 久久人人爽av亚洲精品天堂| 国产精品香港三级国产av潘金莲| 久久伊人香网站| 国产精品免费一区二区三区在线| 怎么达到女性高潮| 久久中文字幕一级| 一级a爱视频在线免费观看| 久久九九热精品免费| 亚洲 国产 在线| 日本撒尿小便嘘嘘汇集6| 日本欧美视频一区| 在线播放国产精品三级| 国产欧美日韩精品亚洲av| 9191精品国产免费久久| av网站免费在线观看视频| 一边摸一边做爽爽视频免费| 黄色片一级片一级黄色片| 亚洲欧美日韩无卡精品| 亚洲天堂国产精品一区在线| 啦啦啦 在线观看视频| 好看av亚洲va欧美ⅴa在| 亚洲黑人精品在线| 丝袜在线中文字幕| 欧美中文日本在线观看视频| 1024香蕉在线观看| 一区二区三区国产精品乱码| 欧美激情高清一区二区三区| 精品少妇一区二区三区视频日本电影| 国产欧美日韩一区二区精品| 波多野结衣高清无吗| 成人三级做爰电影| 黄片小视频在线播放| 99国产精品99久久久久| 久久久久久久精品吃奶| 精品国产亚洲在线| 非洲黑人性xxxx精品又粗又长| 午夜福利18| 在线播放国产精品三级| 色综合站精品国产| 一边摸一边抽搐一进一出视频| 青草久久国产| 涩涩av久久男人的天堂| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美三级三区| 国产熟女xx| 久久精品91无色码中文字幕| 男人舔女人的私密视频| 婷婷丁香在线五月| 亚洲第一青青草原| 在线观看舔阴道视频| 亚洲国产欧美一区二区综合| 久久婷婷成人综合色麻豆| 最新在线观看一区二区三区| 黄色女人牲交| 精品久久久久久久人妻蜜臀av | 91麻豆av在线| 天堂动漫精品| 又黄又粗又硬又大视频| 午夜福利视频1000在线观看 | 熟女少妇亚洲综合色aaa.| 男女床上黄色一级片免费看| 国产精品亚洲美女久久久| 一边摸一边抽搐一进一小说| 国产精品综合久久久久久久免费 | 神马国产精品三级电影在线观看 | 十分钟在线观看高清视频www| 亚洲片人在线观看| 免费一级毛片在线播放高清视频 | 亚洲 欧美 日韩 在线 免费| e午夜精品久久久久久久| 丰满的人妻完整版| 精品卡一卡二卡四卡免费| 老鸭窝网址在线观看| 老汉色∧v一级毛片| 夜夜夜夜夜久久久久| 久久天堂一区二区三区四区| 99久久99久久久精品蜜桃| 日本 欧美在线| 日韩免费av在线播放| 一卡2卡三卡四卡精品乱码亚洲| 亚洲熟妇中文字幕五十中出| 国产精品乱码一区二三区的特点 | 久久久国产欧美日韩av| 亚洲情色 制服丝袜| 免费看美女性在线毛片视频| 亚洲 欧美 日韩 在线 免费| 国产97色在线日韩免费| 午夜福利,免费看| 9191精品国产免费久久| 日韩大码丰满熟妇| 91成人精品电影| 精品第一国产精品| 国产亚洲精品久久久久5区| 激情视频va一区二区三区| 久久久久久亚洲精品国产蜜桃av| 精品乱码久久久久久99久播| 成人永久免费在线观看视频| 亚洲成av人片免费观看| 国产在线观看jvid| 国产99久久九九免费精品| 亚洲av熟女| 老熟妇仑乱视频hdxx| 国产精品久久电影中文字幕| 久久久国产精品麻豆| 久热爱精品视频在线9| 精品一区二区三区四区五区乱码| 一进一出好大好爽视频| 国产亚洲精品综合一区在线观看 | 最近最新中文字幕大全免费视频| 妹子高潮喷水视频| 欧美中文日本在线观看视频| 一二三四社区在线视频社区8| 中文字幕精品免费在线观看视频| 国产av又大| 美女大奶头视频| 国产成人av激情在线播放| 亚洲七黄色美女视频| 少妇粗大呻吟视频| 欧美另类亚洲清纯唯美| 国产一区二区在线av高清观看| 97人妻天天添夜夜摸| 国产欧美日韩一区二区精品| avwww免费| 久久 成人 亚洲| 中亚洲国语对白在线视频| 好男人电影高清在线观看| www.精华液| 国产一区二区三区在线臀色熟女| 亚洲中文字幕日韩| 欧美一级毛片孕妇| 久久久久久久午夜电影| 一边摸一边做爽爽视频免费| 亚洲成国产人片在线观看| 亚洲熟妇熟女久久| 国产精品免费视频内射| 女人精品久久久久毛片| 一二三四社区在线视频社区8| 欧美不卡视频在线免费观看 | 夜夜夜夜夜久久久久| 欧美黄色片欧美黄色片| 岛国在线观看网站| 国产精品秋霞免费鲁丝片| 黄片小视频在线播放| 欧美国产精品va在线观看不卡| 可以在线观看的亚洲视频| 免费av毛片视频| 熟妇人妻久久中文字幕3abv| 99精品久久久久人妻精品| 国产免费av片在线观看野外av| 国产成人免费无遮挡视频| 啪啪无遮挡十八禁网站| 免费看a级黄色片| 亚洲欧美日韩高清在线视频| 高清黄色对白视频在线免费看| 嫩草影视91久久| 欧美久久黑人一区二区| 日韩av在线大香蕉| 午夜免费成人在线视频| 我的亚洲天堂| 电影成人av| 91成人精品电影| 成人欧美大片| 国产成人影院久久av| 亚洲自拍偷在线| 丁香六月欧美| 麻豆av在线久日| 97人妻精品一区二区三区麻豆 | 亚洲av片天天在线观看| 国产色视频综合| 一a级毛片在线观看| 精品国产国语对白av| 波多野结衣一区麻豆| 久久亚洲精品不卡| 午夜福利影视在线免费观看| 别揉我奶头~嗯~啊~动态视频| 欧美午夜高清在线| 色婷婷久久久亚洲欧美| 国内精品久久久久久久电影|