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

    棉花中不同植物病毒介導的VIGE體系的研究

    2023-02-27 02:00:42雷建峰代培紅尤揚子賈建國曲延英劉曉東
    作物學報 2023年4期
    關鍵詞:內源基因組棉花

    雷建峰 李 月 代培紅 趙 燚 尤揚子 賈建國 趙 帥 曲延英,* 劉曉東,*

    棉花中不同植物病毒介導的VIGE體系的研究

    雷建峰1李 月2代培紅2趙 燚2尤揚子2賈建國1趙 帥1曲延英1,*劉曉東2,*

    1新疆農業(yè)大學農學院 / 教育部棉花工程研究中心, 新疆烏魯木齊 830052;2新疆農業(yè)大學生命科學學院, 新疆烏魯木齊 830052

    植物病毒介導sgRNA的傳遞與表達相比傳統轉化完整的編輯載體進行基因編輯具有巨大優(yōu)勢, 因為sgRNA的表達可以伴隨著病毒的復制和移動而快速擴增、累積, 從而產生高效的基因編輯效率。為在棉花中開發(fā)應用更多植物病毒介導的基因編輯系統(virus-induced genome editing, VIGE)。本研究以超表達Cas9 (Cas9-OE)棉花作為VIGE受體, 在棉花中分別建立了棉花葉皺縮病毒(cotton leaf crumple virus, CLCrV)和煙草脆裂病毒(tobacco rattle virus, TRV)介導的VIGE系統。首先CLCrV和TRV介導的VIGE均可以靶向敲除棉花和基因A亞組和D亞組基因組序列, 驗證了這2個系統的可行性和有效性。進一步量化這2個系統對和基因的編輯效率結果發(fā)現, CLCrV和TRV介導的VIGE均能夠產生高效的基因編輯效率, 并且2種系統基因編輯效率不存在顯著差異。本研究還探究了使用棉花內源U6啟動子和擬南芥U6啟動子驅動sgRNA對VIGE基因編輯效率的影響。結果顯示, 這2個啟動子介導的VIGE均能夠產生高效的基因編輯效率, 并且基于這2種啟動子介導的VIGE系統基因編輯效率不存在顯著差異。以上結果預示著可以開發(fā)更多的植物病毒載體應用于棉花基因編輯的研究, 進而獲得高效的sgRNA篩選體系。

    棉花; CLCrV; TRV; VIGE; 基因編輯

    棉花基因功能的鑒定和新品種的選育離不開遺傳變異突變體的獲得, 傳統的誘變技術和T-DNA插入等隨機突變遠遠不能滿足科學研究和棉花育種的需求?;蚪M編輯技術是一種強大的技術工具, 它可以利用可編程序列和特異性核酸內切酶(Cas9或Cpf1)在生物體基因組特定位點上產生靶向突變[1-3]。近年來, II型CRISPR/Cas9系統被廣泛應用于農作物的基因編輯[4-7]。該系統由single guide RNA (sgRNA)引導核酸內切酶Cas9在非常精確的目標位點上切斷雙鏈DNA (dsDNA), 然后生物體通過非同源末端連接(non-homologous end joining, NHEJ)或同源重組(homologous recombination, HR)進行修復。在缺乏同源DNA模板的情況下, NHEJ途徑在修復DSBs過程中可引起目標位點堿基的缺失、插入或替換, 從而在特定位點產生遺傳變異[8]。

    目前, 棉花基因編輯技術體系已經建立并陸續(xù)得到升級和改良(CRISPR/Cpf1和單堿基編輯系統)[9-13]。然而對于遺傳轉化困難且周期長的棉花而言, 要獲得基因編輯的突變體是困難的。即使可以實現遺傳轉化, 也只是局限于該物種的少數幾個基因型。另外, 研究發(fā)現作為CRISPR系統關鍵元件之一的sgRNA, 在很多基因組位點上是低效或者是無效的。而sgRNAs的效率會直接影響到CRISPR/Cas9在棉花中的有效應用。通過生物信息學的方法很難預測sgRNAs的活性[14-15]。因此, 在采用CRISPR系統獲得突變體之前, 必須盡量減少使用無效sgRNA的風險, 就需要快速篩選出高效的sgRNA。在棉花中將完整的基因編輯載體進行原生質體轉化[16]和農桿菌介導的瞬時轉化棉花葉片[17]是2種最常見驗證sgRNA活性的方法。然而, 這些方法檢測到的突變效率低, 甚至很多時候無法有效驗證sgRNA的活性。

    近期, 植物病毒誘導的基因編輯(virus-induced genome editing, VIGE)被廣泛應用于一些模式植物和農作物基因功能的研究[18-22]。該系統需要一個穩(wěn)定超表達Cas9 (Cas9-OE)的轉基因植株作為VIGE受體, 同時將sgRNA表達元件組裝在病毒載體上, 借助病毒可以在植物體內系統侵染、復制、擴散的能力投送sgRNA, 從而高效地實現靶基因的定向編輯。棉花葉皺縮病毒(CLCrV, DNA病毒)和煙草脆裂病毒(TRV, RNA病毒)是2種常用于棉花基因沉默的病毒載體, 可用于沉默棉花內源基因的表達[23-24]。本課題組在前期工作中, 對DNA病毒CLCrV載體進行改造分別在擬南芥和棉花中建立了CLCrV介導的VIGE系統[22,25]。在棉花中, RNA病毒TRV是否也可以作為一個理想的病毒載體用于構建VIGE系統, 能否高效地實現對棉花內源基因的靶向編輯尚缺乏相關研究。與此同時, 在前期的工作中, 本課題組使用截短的擬南芥AtU6-26啟動子(330 bp)驅動sgRNA表達, 在棉花中獲得了比較理想的基因編輯效率。使用棉花內源U6啟動子驅動sgRNA表達能否進一步提升VIGE的基因編輯效率, 相關研究也鮮有報道。因此, 在棉花中建立不同植物病毒介導的VIGE系統和探究不同啟動子驅動sgRNA表達對基因編輯效率的影響, 對篩選出高效的VIGE體系具有重要意義。

    本研究基于DNA病毒CLCrV和RNA病毒TRV載體, 構建了相應的VIGE體系。采用農桿菌注射法將不同sgRNA衍生物接種棉花Cas9-OE植株, 驗證這2種VIGE系統的有效性和高效性; 與此同時, 探究不同啟動子對VIGE基因編輯效率的影響, 以篩選出適用于棉花基因編輯高效的VIGE系統。

    1 材料與方法

    1.1 試驗材料

    試驗所用AtU6-26P::sgRNA、GhU6-5P::sgRNA[26]、棉花葉皺縮病毒CLCrV載體、煙草脆裂病毒TRV載體、棉花Cas9超表達(Cas9-OE)植株YZ-1均為新疆農業(yè)大學農業(yè)生物技術重點實驗室保存。各種限制性內切酶購于Thermo公司; T4DNA連接酶、Blunt Zero平端載體、Trans1-T1感受態(tài)細胞、瓊脂糖凝膠回收試劑盒、DNA聚合酶、植物基因組DNA提取試劑盒、反轉錄試劑盒、熒光定量試劑盒、2K Plus II DNA Marker均購自于北京全式金生物技術有限公司; Phusion超保真DNA聚合酶購自于紐英倫生物技術(北京)有限公司; 其他常規(guī)試劑均為國產分析純; 引物合成及測序均由上海杰李生物技術有限公司完成。

    1.2 qRT-PCR分析Cas9表達量

    利用Biospin Plant Total RNA提取試劑盒從野生型和Cas9-OE棉花子葉中分離總RNA。然后用反轉錄試劑盒One-Step gDNA Removal and cDNA Synthesis SuperMix合成cDNA。通過qRT-PCR擴增188 bp基因片段, 以棉花基因為內參對照(引物設計見表1), 每個樣本設置3個技術重復。反應結束后, 根據目的基因和內參基因的Ct值, 利用2–DDCt方法[27]計算基因表達量。

    1.3 不同植物病毒介導的VIGE體系的建立

    1.3.1 sgRNA中間載體的構建 使用截短的AtU6-26 (330 bp)和GhU6-5P (363 bp)啟動子驅動sgRNA表達。文中所提及的和基因的guide RNA (表1)通過退火: 95℃每5 s降低0.5~20.0℃, 插入到經I酶切的AtU6-26::sgRNA和GhU6-5P::sgRNA載體上, 測序驗證- sgRNA和-sgRNA表達載體是否正確。

    1.3.2 CLCrV-sgRNA和TRV-sgRNA載體構建 用I和I酶切回收AtU6-26::- sgRNA、AtU6-26::-sgRNA、GhU6-5P::-sgRNA和GhU6-5P::-sgRNA片段, 將其分別組裝在CLCrV-A上; 同時用I和I酶切回收AtU6-26::-sgRNA、AtU6-26::-sgRNA、GhU6-5P::- sgRNA和GhU6-5P::-sgRNA片段, 將其分別組裝在TRV-V2上。CLCrV-sgRNA和TRV-sgRNA載體示意圖見圖1。采用酶切鑒定的方法驗證這些重組質粒是否正確, 將鑒定正確的表達載體轉入農桿菌GV3101, 用于下一步侵染。

    表1 本研究中使用的引物序列

    圖1 CLCrV-sgRNA和TRV-sgRNA載體構建

    1.4 不同植物病毒投送sgRNA瞬時轉化Cas9- OE棉花

    Cas9-OE棉花品種YZ-1種子在ddH2O水中浸泡3 d, 然后在黑暗中28℃發(fā)芽48 h。發(fā)芽后, 將幼苗移栽到土壤中, 在28℃、12 h光照/12 h黑暗條件下生長。待棉花幼苗2片子葉完全展開時進行接種轉化試驗。將構建好的不同CLCrV-sgRNA和TRV-sgRNA載體的農桿菌菌液接種于2 mL LB培養(yǎng)基(含50 μg mL–1Kan和50 μg mL–1Rif)中活化, 28℃、200轉 min-1過夜培養(yǎng)。按1∶100比例分別吸取不同活化菌液接種到15 mL LB培養(yǎng)基(含50 μg mL–1Kan和50 μg mL–1Rif)中, 再次28℃、200轉 min-1搖菌, 統一搖至菌液的OD600= 1.6~1.8時, 1000×離心15 min收集農桿菌培養(yǎng)物, 重懸于轉化液(10 mmol L–1MgCl2, 10 mmol L–1MES和200 μmol L–1乙酰丁香酮)中, 調整OD600至1.0左右, 在室溫下避光孵育3~4 h。接種病毒時, 將CLCrV-B和CLCrV-A的衍生物, TRV-V1和TRV-V2的衍生物1∶1等比例混合, 利用1 mL注射器將轉化液注射到Cas9-OE棉花子葉中進行瞬時表達。

    1.5 突變檢測

    sgRNA可伴隨著病毒的復制和擴散在棉花體內系統的表達, 接種相應病毒載體的棉花子葉是sgRNA表達的起點。采用Plant Genomic DNA提取試劑盒分別從接種7~10 d后的CLCrV-B和CLCrV-A衍生物, TRV-V1和TRV-V2衍生物的棉花子葉中提取基因組DNA。利用PCR/RE[17]的方法檢測和基因目標位點是否發(fā)生突變, PCR擴增包含目標位點的基因組片段, 擴增引物見表1。用靶位點的限制性內切酶消化PCR產物。同時, 未被消化的PCR擴增子被克隆到Blunt Zero平端克隆載體上進行測序。

    1.6 不同植物病毒介導的VIGE基因編輯效率分析

    分別提取每一株接種CLCrV-AtU6-26::- sgRNA、TRV-AtU6-26::-sgRNA、CLCrV- AtU6-26::-sgRNA和CLCrV-AtU6-26::-sgRNA棉花葉片基因組DNA。以每一株轉化植株的基因組DNA為模板, 參照和基因組序列進行PCR擴增, 并用靶位點上相應的限制性內切酶消化PCR產物, 檢測每個單株CLCrV和TRV介導的VIGE系統對和基因的編輯效率。

    1.7 不同U6啟動子介導的VIGE基因編輯效率分析

    以棉花基因為靶標, 分別提取每一株接種CLCrV-AtU6-26::-sgRNA和CLCrV-GhU6-5P::-sgRNA棉花基因組DNA。以每一株轉化植株的基因組DNA為模板進行PCR擴增并用靶位點上的限制性內切酶I消化PCR產物, 檢測每個單株在CLCrV介導的VIGE系統下, 使用棉花內源U6啟動子和擬南芥U6啟動子驅動-sgRNA表達對棉花基因的編輯效率。

    2 結果與分析

    2.1 棉花Cas9-OE植株Cas9基因表達量檢測

    以野生型棉花YZ-1和Cas9-OE棉花cDNA為模板, 進行qRT-PCR擴增。成功檢測到不同Cas9-OE植株中基因穩(wěn)定表達(圖2)。qRT-PCR驗證了在轉基因棉花植株中高效表達, 本研究選擇了表達量較高的Cas9-OE: 1~6轉基因株系用于后續(xù)研究。

    圖2 qRT-PCR分析Cas9表達量

    在同一指標中不同小寫字母表示不同處理在0.05概率水平差異顯著。

    Different letters indicate significant difference among different treatments at the 0.05 probability level.

    2.2 CLCrV-sgRNA和TRV-sgRNA載體構建

    將4個AtU6-26::-sgRNA、AtU6-26::-sgRNA、GhU6-5P::-sgRNA和GhU6-5P::-sgRNA分別與CLCrV-A和TRV-V2片段連接, 對8個重組的質粒進行酶切鑒定:用I和I對CLCrV-AtU6-26::- sgRNA、CLCrV-AtU6-26::-sgRNA重組質粒進行雙酶切, 產生了516 bp和11,300 bp 2條帶, 對CLCrV-GhU6-5P::-sgRNA、CLCrV-GhU6- 5P::-sgRNA雙酶切, 產生了514 bp和11,300 bp 2條帶; 用I和I對TRV-AtU6- 26::-sgRNA、TRV-AtU6-26::- sgRNA重組質粒進行雙酶切, 產生了516 bp和11,003 bp 2條帶, 對TRV-GhU6-5P::- sgRNA、TRV-GhU6-5P::-sgRNA雙酶切, 產生了514 bp和11,003 bp 2條帶, 所有酶切結果與預期酶切片段大小相符(圖3), 表明不同CLCrV-sgRNA和TRV-sgRNA載體均已構建成功。

    2.3 CLCrV和TRV介導靶向敲除棉花GhBsrk1和GhMAPKKK2基因

    隨機選取3~4株接種CLCrV-AtU6-26::- sgRNA、CLCrV-AtU6-26::-sgRNA、TRV-AtU6-26::-sgRNA、TRV-AtU6-26::-sgRNA棉花葉片混樣, 提取總DNA。以野生型棉花基因組DNA和混樣DNA為模板, 通過PCR/RE的方法驗證和基因是否在目標位點發(fā)生突變。PCR擴增基因(777 bp), 對I酶切未消化的PCR產物進行克隆測序, 發(fā)現所有突變類型都為靶序列區(qū)域堿基的缺失或插入(圖4-A~C)。PCR擴增基因(607 bp), 對I酶切未消化的PCR產物克隆測序, 發(fā)現所有突變類型均為PAM位點上游堿基的缺失(圖4-D~F)。而在對照野生型棉花中,和基因的PCR產物被相應的限制性內切酶消化完全, 幾乎沒有擴增產物殘留。表明利用CLCrV和TRV介導的VIGE系統均能夠實現對棉花內源基因(A亞組和D亞組)的靶向編輯。

    圖3 不同CLCrV-sgRNA和TRV-sgRNA表達載體的酶切鑒定

    M: 2K Plus II DNA 標準分子量; 1: CLCrV-AtU6-26::- sgRNA; 2: CLCrV-AtU6-26::-sgRNA; 3: CLCrV- GhU6-5P::-sgRNA; 4: CLCrV-GhU6-5P::- sgRNA; 5: TRV-AtU6-26::-sgRNA; 6: TRV-AtU6-26::-sgRNA; 7: TRV-GhU6-5P::G-sgRNA; 8: TRV- GhU6-5P::-sgRNA。

    M: 2K Plus II DNA marker; 1: CLCrV-AtU6-26::-sgRNA; 2: CLCrV-AtU6-26::-sgRNA; 3: CLCrV-GhU6-5P::-sgRNA; 4: CLCrV-GhU6-5P::-sgRNA; 5: TRV-AtU6-26::-sgRNA; 6: TRV-AtU6-26::- sgRNA; 7:TRV-GhU6-5P::G-sgRNA; 8:TRV-GhU6-5P::-sgRNA.

    (圖4)

    M: 2K Plus II DNA標準分子量。A和B:-sgRNA靶向突變檢測。A: WT為對照, 1和2分別為接種CLCrV-AtU6-26::-sgRNA和TRV-AtU6-26::-sgRNA植株編號, 凝膠電泳圖顯示基因經I酶切消化的產物和未消化的PCR產物。B: 未消化的PCR產物缺乏I位點(由于存在突變), 隨后進行純化、克隆和測序分析。C:突變測序峰圖; D和E:-sgRNA靶向突變檢測。D: WT為對照, 1和2分別為接種CLCrV-AtU6-26::-sgRNA和TRV-AtU6-26::-sgRNA植株編號, 凝膠電泳圖顯示基因經I酶切消化的產物和未消化的PCR產物。E: 未消化的PCR產物缺乏I位點(由于存在突變), 隨后進行純化、克隆和測序分析。F:突變測序峰圖。PAM位點為綠色, 藍色下畫線表示靶序列上的酶切位點。‘M’表示突變序列, ‘+’堿基插入用紅色大寫字母表示, ‘-’堿基缺失用紅色線段表示。

    M: 2K Plus II DNA marker. A and B: the detection of-sgRNA targeted mutations. A: wild type served as a control, 1 and 2 are plant numbers inoculated with CLCrV-AtU6-26::-sgRNA and TRV-AtU6-26::G-sgRNA respectively. The gel image shows digested PCR products ofgene withI and undigested PCR products. B: the undigested PCR products lacking theI site (due to the presence of a mutation) that were subsequently purified, cloned, and analyzed by sequencing. C:mutation sequencing peak map; D and E: the detection of-sgRNA targeted mutations. D: wild type served as a control, 1 and 2 are plant numbers inoculated with CLCrV-AtU6-26::-sgRNA and TRV-AtU6-26::-sgRNA, respectively. The gel image shows digested PCR products ofgene withI and undigested PCR products. E: the undigested PCR products lacking theI site (due to the presence of a mutation) that were subsequently purified, cloned, and analyzed by sequencing. F:mutation sequencing peak map. Green color indicates the PAM sequence, underline in blue indicates the restriction site on the target sequence. M indicates the mutation sequence, insertions are denoted with red uppercase letters, deletions are shown as red dashes.

    2.4 CLCrV和TRV介導的VIGE基因編輯效率分析

    為進一步量化CLCrV和TRV介導的VIGE基因編輯效率, 提取每一株接種CLCrV-AtU6-26::-sgRNA (8株)、CLCrV-AtU6-26::- sgRNA (12株)、TRV-AtU6-26::-sgRNA (7株)、TRV-AtU6-26::-sgRNA (13株)載體的棉花基因組DNA。對每個單株進行突變檢測, 以野生型為對照。在使用相同啟動子(AtU6-26)驅動-sgRNA和-sgRNA表達情況下, PCR/RE檢測結果顯示: 對于基因, CLCrV介導的VIGE基因編輯效率為100% (12/12), TRV介導的VIGE基因編輯效率為92.31% (12/13); 對于基因, CLCrV介導的VIGE基因編輯效率為75% (6/8), TRV介導的VIGE基因編輯效率為85.71% (6/7) (圖5)。表明在棉花中針對同一個基因(或者)進行定向編輯, 利用CLCrV和TRV介導的VIGE系統均能夠產生高效的基因編輯效率, 并且2種系統基因編輯效率不存在顯著差異。

    2.5 不同U6啟動子介導的VIGE基因編輯效率分析

    為探究使用棉花內源GhU6-5P啟動子和擬南芥AtU6-26啟動子介導的VIGE系統的基因編輯效率是否存在差異。以基因為靶標, 提取每一株接種CLCrV-AtU6-26::-sgRNA (12株)、CLCrV-GhU6-5P::-sgRNA (11株)載體的棉花基因組DNA。對每個單株進行突變檢測, 以野生型為對照。由于CLCrV和TRV介導的VIGE基因編輯效率無明顯差異, 本研究在CLCrV介導的VIGE下, 探究不同U6啟動子介導的VIGE基因編輯效率。PCR/RE突變檢測結果顯示: 擬南芥AtU6-26啟動子介導的VIGE基因編輯效率為100% (12/12), 而棉花GhU6-5P啟動子介導的VIGE基因編輯效率同樣為100% (11/11) (圖6)。表明在棉花中使用棉花內源U6啟動子和擬南芥U6啟動子驅動sgRNA表達均能夠產生高效的基因編輯效率, 并且基于這2種啟動子介導的VIGE基因編輯效率不存在顯著差異。這暗示著在同源關系較遠的不同物種中, U6啟動子是可以替換使用的。

    圖5 CLCrV和TRV介導的VIGE基因編輯效率比較

    M: 2K Plus II DNA標準分子量。A: WT為對照, 1~12為接種CLCrV-AtU6-26::-sgRNA單株編號, 13~25為接種TRV-AtU6-26::-sgRNA單株編號。B: WT為對照, 1~8為接種CLCrV-AtU6-26::-sgRNA單株編號, 9~15為接種TRV-AtU6-26::-sgRNA單株編號。

    M: 2K Plus II DNA marker. A: wild type, 1–12 are individual plant numbers inoculated with CLCrV-AtU6-26::-sgRNA, 13–25 are individual plant numbers inoculated with TRV-AtU6-26::-sgRNA. B: wild type, 1–8 are individual plant numbers inoculated with CLCrV-AtU6-26::-sgRNA, 9–15 are individual plant numbers inoculated with TRV-AtU6-26::-sgRNA.

    圖6 不同U6啟動子介導的VIGE基因編輯效率比較

    M: 2K Plus II DNA標準分子量。A: WT為對照, 1~12為接種CLCrV-AtU6-26::-sgRNA單株編號。B: WT為對照, 1~11為接種CLCrV-GhU6-5P::-sgRNA單株編號。

    M: 2K Plus II DNA marker. A: wild type, 1–12 are individual plant numbers inoculated with CLCrV-AtU6-26::-sgRNA; B: wild type , 1–11 are individual plant numbers inoculated with CLCrV-GhU6-5P::-sgRNA.

    3 討論

    棉花全基因組測序已經完成, 然而棉花功能基因的挖掘以及全基因組序列的生物學意義還知之甚少。從根源上解析棉花基因功能, 必需要以棉花基因組變異的突變體材料作為研究基礎。對于遺傳轉化困難且周期長的棉花而言, 獲得少數幾個基因編輯的棉花是經濟可行的, 但要大規(guī)模構建棉花基因突變體庫和實現全基因組水平上的精準遺傳改良是費錢、費時, 不可想象的。鑒于獲得穩(wěn)定遺傳棉花突變體的時間和難度, 開發(fā)一種快速驗證sgRNAs活性, 進而獲得高效基因編輯效率的方法至關重要。

    利用病毒載體進行基因編輯早已被人們認知, 但是由于病毒龐大基因組的限制, 嚴重阻礙了利用病毒載體進行基因編輯的研究。對于DNA病毒CLCrV和RNA病毒TRV而言, 研究者對病毒基因組進行改造將其組裝在植物表達載體上, 使之成為以農桿菌接種的方式侵染棉花葉片, 相繼在陸地棉和海島棉中建立了病毒誘導的基因沉默(VIGS)體系[23-24,28]。雖然這2種病毒載體受限于病毒承載力的限制, 不適于表達Cas9蛋白, 但足以表達sgRNA。以Cas9-OE植株作為VIGE受體, 利用病毒載體單獨投送sgRNA可避免轉化完整編輯載體編輯效率低的弊端。本研究為了擴大這2種病毒載體在棉花中的應用范圍, 分別在棉花中建立CLCrV和TRV介導的VIGE系統進行測試。在VIGE體系內, 要想獲得高效的基因編輯效率,基因的高表達是首要關鍵因素。本研究首先檢測了棉花Cas9-OE植株中的表達量發(fā)現,基因在所有Cas9-OE植株中高效表達(圖2), 這為驗證VIGE的有效性和高效性提供了良好的開端。利用CLCrV和TRV病毒載體投送-sgRNA和-sgRNA, 突變分析結果顯示, 在使用相同啟動子(AtU6-26)驅動sgRNA表達的情況下, CLCrV和TRV介導的VIGE均可以實現棉花內源基因(A亞組和D亞組)的靶向編輯, 突變類型主要以堿基的缺失和插入為主(圖4),相比以往本實驗室通過構建完整的編輯載體轉化棉花原生質體獲得較低的編輯效率(大多數突變類型為點突變)[16], 這2種VIGE系統大大提高了基因編輯的檢出率。為進一步量化CLCrV和TRV介導的VIGE基因編輯效率, 對每一個轉化單株進行突變檢測, 發(fā)現CLCrV和TRV介導的VIGE均能夠產生高效的基因編輯效率, 并且2種系統基因編輯效率不存在顯著差異(圖5)。另外, 作為VIGE系統關鍵元件之一的sgRNA, 使用棉花內源U6啟動子驅動sgRNA表達能否進一步提升VIGE的基因編輯效率?本研究繼續(xù)探究截短的棉花GhU6-5P (363 bp)和擬南芥AtU6-26 (330 bp)啟動子介導的VIGE (CLCrV)對每一個轉化單株的基因編輯情況發(fā)現, 棉花內源U6啟動子和擬南芥U6啟動子介導的VIGE均能夠產生高效的基因編輯效率, 并且基于這2種啟動子介導的VIGE基因編輯效率不存在顯著差異(圖6)。與此同時, 本研究結果也證明了U6啟動子在同源關系較遠的不同物種之間是可以替換使用的。在今后的研究中CLCrV和TRV介導的VIGE系統或許還可以借助核糖核酸酶或者tRNA的幫助使用一個U6啟動子同時表達多個sgRNAs, 實現多基因編輯。

    另外, 本研究還通過Hi-TOM高通量測序檢測了部分突變單株的真實編輯效率, 突變效率在15.74%~42.68%之間。推測造成不同植株編輯效率存在差異的因素可能與植株自身基因表達量以及sgRNA伴隨著病毒復制、擴散能力有關, sgRNA并不是在所有病毒侵染的組織細胞中都表達產生靶向突變。盡管病毒可以高效投送sgRNAs, 實現靶基因的定向編輯。但該事件只能發(fā)生在當代植株體內, 由于植物自身存在某種特殊機制使得病毒難以進入到植物的頂端分生組織(SAM)中[29], 實現植物生殖細胞的編輯, 因此基因編輯無法遺傳給后代。近期有研究表明將可移動RNA元件(和tRNA)與sgRNA融合后克隆到TRV病毒載體中, 能夠將sgRNA長距離運輸至煙草的SAM中, 并以較高的效率產生可遺傳的基因編輯[30]。除此之外, 研究還發(fā)現不修飾的BSMV病毒傳遞sgRNA能夠直接在超表達Cas9的小麥中獲得可遺傳的基因編輯[31]。這些研究成果為利用VIGE系統在棉花中實現可遺傳的基因編輯后代提供了借鑒, 后期的工作仍需開展突破棉花遺傳轉化途徑限制的研究, 從而獲得棉花突變體進行功能驗證。

    以上結果顯示基于CLCrV和TRV介導的VIGE系統在棉花上的成功應用, 預示著通過將更多種病毒作為VIGE載體, 在棉花上進行驗證是有可能篩選到適合的VIGE系統的, 這為定向創(chuàng)制棉花突變體材料奠定了堅實的技術基礎。

    4 結論

    在棉花中分別建立了CLCrV和TRV介導的VIGE系統, 這2個系統均可以實現棉花內源基因的定向編輯, 且基因編輯效率不存在顯著差異; 使用棉花內源U6啟動子或者異源U6啟動子驅動sgRNA表達都可以獲得高效的基因編輯效率。研究結果預示著將來可以開發(fā)更多的植物病毒應用于農作物基因編輯的研究。

    [1] Chen K L, Gao C X. Targeted genome modification technologies and their applications in crop improvements., 2014, 33: 575–583.

    [2] Zhou H B, Liu B, Weeks D P,Spalding M H, Yang B. Large chromosomal deletions and heritable small genetic changes induced by CRISPR/Cas9 in rice., 2015, 42: 10903–10914.

    [3] Alok A, Sandhya D, Jogam P, Rodrigues V, Bhati K K, Sharma H, Kumar J. The rise of the CRISPR/Cpf1 system for efficient genome editing in plants., 2020, 11: 264.

    [4] Wang C, Liu Q, Shen Y, Hua Y F, Wang J J, Lin J R, Wu M G, Sun T T, Cheng Z K, Mercier R, Wang K J. Clonal seeds from hybrid rice by simultaneous genome engineering of meiosis and fertilization genes., 2019, 37: 283–286.

    [5] Zhang S J, Zhang R Z, Song G Q, Gao J, Li W, Han X D, Chen M L, Li Y L, Li G Y. Targeted mutagenesis using the-mediated CRISPR-Cas9 system in common wheat., 2018, 18: 302.

    [6] Wang D D, Samsulrizal N, Yan C, Allcock N S, Craigon J, Blanco-Ulate B, Ortega-Salazar I, Marcus S E, Bagheri H M, Fons L P, Fraser P D, Foster T, Fray R, Knox J P, Seymour G B. Characterization of CRISPR mutants targeting genes modulating pectin degradation in ripeningtomato., 2019, 179: 544–557.

    [7] Osakabe Y, Liang Z C, Ren C,Nishitani C, Osakabe K, Wada M, Komori S, Malnoy M, Velasco R, Poli M, Jung M H, Koo O J, Viola R, Kanchiswamy C N. CRISPR-Cas9-mediated genome editing inapple and grapevine., 2018, 13: 2844–2863.

    [8] Jiang W Z, Zhou H B, Bi H H, Fromm M, Yang B, Weeks D P. Demonstration of CRISPR/Cas9/sgRNA-mediated targeted gene modification in, tobacco, sorghum and rice., 2013, 41: e188.

    [9] Chen X G, Lu X K, Shu N, Wang S, Wang J J, Wang D L, Guo L X, Ye W W. Targeted mutagenesis in cotton (L.) using the CRISPR/Cas9 system., 2017, 7: 44304.

    [10] Long L, Guo D D, Gao W, Yang W W, Hou L P, Ma X N, Miao Y C, Botella J R, Song C P. Optimization of CRISPR/Cas9 genome editing in cotton by improved sgRNA expression., 2018, 14: 85.

    [11] Manghwar H, Lindsey K, Zhang X L, Jin S X. CRISPR/Cas system: recent advances and future prospects for genome editing., 2019, 24: 1102–1125.

    [12] Li B, Rui H P, Li Y J, Wang Q Q, Alariqi M, Qin L, Sun L, Ding X, Wang F Q, Zou J W, Wang Y Q, Yuan D J, Zhang X L, Jin S X. Robust CRISPR/Cpf1 (Cas12a)-mediated genome editing in allotetraploid cotton ()., 2019, 17: 1862–1864.

    [13] Qin L, Li J Y, Wang Q Q, Xu Z P, Sun L, Alariqi M, Manghwar H, Wang G Y, Li B, Ding X, Rui H P, Huang H M, Lu T L, Lindsey K, Daniell H, Zhang X L, Jin S X. High-efficient and precise base editing of C?G to T?A in the allotetraploid cotton () genome using a modified CRISPR/Cas9 system., 2019, 18: 45–56.

    [14] Zhao C Z, Wang Y L, Nie X W, Han X S, Liu H L, Li G L, Yang G J, Ruan J X, Ma Y L, Li X Y, Cheng H J, Zhao S H, Fang Y P, Xie S S. Evaluation of the effects of sequence length and microsatellite instability on single-guide RNA activity and specificity., 2019, 15: 2641–2653.

    [15] Xie S S, Shen B, Zhang C B, Huang X X, Zhang Y L. sgRNAcas9: a software package for designing CRISPR sgRNA and evaluating potential off-target cleavage sites., 2014, 9: e100448.

    [16] 李繼洋, 雷建峰, 代培紅, 姚瑞, 曲延英, 陳全家, 李月, 劉曉東. 基于棉花U6啟動子的海島棉CRISPR/Cas9基因組編輯體系的建立. 作物學報, 2018, 44: 227–235.

    Li J Y, Lei J F, Dai P H, Yao R, Qu Y Y, Chen Q J, Li Y, Liu X D. Establishment of CRISPR/Cas9 genome editing system based on GbU6 promoters in cotton (L.)., 2018, 44: 227–235 (in Chinese with English abstract).

    [17] Gao W, Long L, Tian X Q, Xu F C, Liu J, Singh P K, Botella J R, Song C P. Genome editing in cotton with the CRISPR/Cas9 system., 2017, 8: 1364.

    [18] Yin K Q, Han T, Liu G, Chen T Y, Wang Y, Alice Y Y, Liu Y L. A geminivirus-based guide RNA delivery system for CRISPR/Cas9 mediated plant genome editing., 2015, 5: 14926.

    [19] Hu J C, Li S, Li Z L, Li H Y, Song W B, Zhao H M, Lai J S, Xia L Q, Li D W, Zhang Y L. A barley stripe mosaic virus-based guide RNA delivery system for targeted mutagenesis in wheat and maize., 2019, 20: 1463–1474.

    [20] Li T D, Hu J C, Sun Y,Li B S, Zhang D L, Li W L, Liu J X, Li D W, Gao C X, Zhang Y L, Wang Y P. Highly efficient heritable genome editing inwheat using an RNA virus and bypassing tissue culture., 2021, 14: 1787–1798.

    [21] Luo Y J, Na R, Nowak J S, Qiu Y, Lu Q S, Yang C Y, Marsolais F, Tian L N. Development of a-processed guide RNA delivery system with soybean-infecting virus ALSV for genome editing., 2021, 21: 419.

    [22] Lei J F, Dai P H, Li Y, Zhang W Q, Zhou G T, Liu C, Liu X D. Heritable gene editing usingmobile guide RNAs and DNA viruses., 2021, 17: 20.

    [23] Gu Z H, Huang C J, Li F F, Zhou X P. A versatile system for functional analysis of genes and microRNAs in cotton., 2014, 12: 638–649.

    [24] Zhang J X, Wang F R, Zhang C Y, Zhang J H, Chen Y, Liu G D, Zhao Y X, Hao F S, Zhang J. A novel VIGS method by agroinoculation of cotton seeds and application for elucidating functions ofin salt-stress response., 2018, 37: 1091–1100.

    [25] 周冠彤, 雷建峰, 代培紅, 劉超, 李月, 劉曉東. 棉花CRISPR/Cas9基因編輯有效sgRNA高效篩選體系的研究. 作物學報, 2021, 47: 427–437.

    Zhou G T, Lei J F, Dai P H, Liu C, Li Y, Liu X D. Efficient screening system of effective sgRNA for cotton CRISPR/Cas9 gene editing., 2021, 47: 427–437 (in Chinese with English abstract)

    [26] 雷建峰, 李月, 徐新霞, 阿爾祖古麗·塔什, 蒲艷, 張巨松, 劉曉東. 棉花不同GbU6啟動子截短克隆及功能鑒定. 作物學報, 2016, 42: 675–683.

    Lei J F, Li Y, Xu X X, A’erzuguli T, Pu Y, Zhang J S, Liu X D. Cloning and functional analysis of different truncated GbU6 promoters in cotton., 2016, 42: 675–683 (in Chinese with English abstract).

    [27] Livak K J, Schmittgen T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2–ΔΔCTmethod., 2001, 25: 402–408.

    [28] Wang W N, Yuan Y L, Yang C, Geng S P, Sun Q, Long L, Cai C W, Chu Z Y, Liu X, Wang G H, Du X M, Miao C, Zhang X, Cai Y F.Characterization, expression, and functional analysis of a novelgene associated with resistance to Verticillium wilt and abiotic stress in cotton., 2016, 6: 3951–3961.

    [29] Wu H J, Qu X Y, Dong Z C, Luo L J, Shao C, Forner J, Lohmann J U, Su M, Xu M C, Liu X B, Zhu L, Zeng J, Liu S M, Tian Z X, Zhao Z.triggers innate antiviral immunity in plant stem cells., 2020, 370: 227–231.

    [30] Ellison E E, Nagalakshmi U, Gamo M E, Huang P J, Dinesh-Kumar S, Voytas D F. Multiplexed heritable gene editing using RNA viruses and mobile single guide RNAs., 2020, 6: 620–624.

    [31] Li T D, Hu J C, Sun Y, Li B S, Zhang D L, Li W L, Liu J X, Li D W, Gao C X, Zhang Y L, Wang Y P. Highly efficient heritable genome editing in wheat using an RNA virus and bypassing tissue culture., 2021, 14: 1787–1798.

    Study on VIGE system mediated by different plant viruses in cotton

    LEI Jian-Feng1, LI Yue2, DAI Pei-Hong2, ZHAO Yi2, YOU Yang-Zi2, JIA Jian-Guo1, ZHAO Shuai1, QU Yan-Ying1,*, and LIU Xiao-Dong2,*

    1College of Agronomy, Xinjiang Agricultural University / Research Center of Cotton Engineering, Ministry of Education, Urumqi 830052, Xinjiang, China;2College of Life Sciences, Xinjiang Agricultural University, Urumqi 830052, Xinjiang, China

    Plant virus-mediated sgRNA delivery and expression has tremendous advantage over traditional transformation of intact editing vectors for gene editing, because sgRNA expression can be rapidly amplified and be accumulated along with virus replication and movement, resulting in efficient gene editing efficiency. In this study, to develop and apply more plant virus-mediated gene editing (VIGE) systems in cotton, the overexpressing Cas9 (Cas9-OE) cotton was used as VIGE receptor and the cotton leaf crumple virus (CLCrV) and tobacco rattle virus (TRV) mediated VIGE systems were established in cotton. First of all, both CLCrV and TRV-mediated VIGE could target and knock out theandin subgroup A and subgroup D genomic sequences, which verified the feasibility and effectiveness of these two systems. Further quantificaiton of the editing efficiency ofandgenes by these two systems showed that both CLCrV and TRV-mediated VIGE could produce high-efficiency gene editing efficiency, and there was no significant difference in gene editing efficiency between the two systems. This study also explored the effect of sgRNAs driven by cotton endogenous U6 promoter andU6 promoter on VIGEgene editing efficiency. The results showed that both promoter-mediated VIGE were able to produce high-efficiency gene editing efficiency, and there was no significant difference in gene editing efficiency based on these two promoter-mediated VIGE systems. The above results indicate that more plant virus vectors can be developed for the application of cotton gene editing research, thus obtaining an efficient sgRNA screening system.

    cotton; CLCrV; TRV; VIGE; gene editing

    10.3724/SP.J.1006.2023.24071

    本研究由新疆維吾爾自治區(qū)自然科學基金項目(2022D01B23)和自治區(qū)高?;究蒲袠I(yè)務費科研項目(XJEDU2022J007)資助。

    This study was supported by the Natural Science Foundation of Xinjiang Uygur Autonomous Region (2022D01B23) and the Basic Research Funds for Universities in the Autonomous Region (XJEDU2022J007).

    劉曉東, E-mail: xiaodongliu75@aliyun.com; 曲延英, E-mail: xjyyq5322@126.com

    E-mail: kyleijianfeng@163.com

    2022-03-30;

    2022-07-21;

    2022-08-17.

    URL: https://kns.cnki.net/kcms/detail/11.1809.S.20220816.1830.008.html

    This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

    猜你喜歡
    內源基因組棉花
    棉花是花嗎?
    牛參考基因組中發(fā)現被忽視基因
    棉花
    小讀者(2020年4期)2020-06-16 03:33:54
    內源多胺含量在砂梨果實發(fā)育過程中的變化研究
    內源信號肽DSE4介導頭孢菌素C酰化酶在畢赤酵母中的分泌表達
    一次注射15N-亮氨酸示蹤法檢測雞內源氨基酸損失量適宜參數的研究
    不同層積處理滇重樓種子內源激素變化的研究
    心中的“棉花糖”
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    第三講 棉花肥害診斷及其防治
    色婷婷久久久亚洲欧美| 最新美女视频免费是黄的| www国产在线视频色| 成年女人毛片免费观看观看9| 岛国视频午夜一区免费看| 国产亚洲精品久久久久5区| 亚洲一卡2卡3卡4卡5卡精品中文| 久久午夜综合久久蜜桃| 激情在线观看视频在线高清| 午夜福利欧美成人| 最近最新中文字幕大全免费视频| 免费在线观看亚洲国产| 亚洲情色 制服丝袜| 国产精品免费视频内射| 精品国产美女av久久久久小说| 亚洲九九香蕉| 成年女人毛片免费观看观看9| 咕卡用的链子| 此物有八面人人有两片| av在线天堂中文字幕| 亚洲五月婷婷丁香| 9色porny在线观看| 天天躁夜夜躁狠狠躁躁| 免费高清在线观看日韩| 男人操女人黄网站| 高清在线国产一区| 一级a爱视频在线免费观看| 亚洲第一av免费看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲国产看品久久| 搡老妇女老女人老熟妇| 欧美中文日本在线观看视频| 国产精品亚洲美女久久久| 99国产极品粉嫩在线观看| 免费高清视频大片| 国产精品爽爽va在线观看网站 | 性少妇av在线| 成年人黄色毛片网站| 精品国产乱码久久久久久男人| 亚洲国产精品999在线| 久久狼人影院| 亚洲一区中文字幕在线| 又黄又爽又免费观看的视频| 变态另类成人亚洲欧美熟女 | 丝袜美腿诱惑在线| 夜夜爽天天搞| 久久九九热精品免费| 成年女人毛片免费观看观看9| 美女扒开内裤让男人捅视频| 级片在线观看| 成人三级做爰电影| 亚洲成人国产一区在线观看| 成人精品一区二区免费| 亚洲天堂国产精品一区在线| 亚洲精品国产区一区二| 在线天堂中文资源库| 成人手机av| 正在播放国产对白刺激| 国产精品99久久99久久久不卡| 国产野战对白在线观看| 国产成人精品久久二区二区91| 国产区一区二久久| 日韩精品中文字幕看吧| 久久狼人影院| 一区二区三区精品91| 最近最新中文字幕大全电影3 | 黄色毛片三级朝国网站| 日韩欧美三级三区| 岛国视频午夜一区免费看| 男人的好看免费观看在线视频 | 久久久久国内视频| 日本三级黄在线观看| 色老头精品视频在线观看| 亚洲国产精品sss在线观看| 中文字幕另类日韩欧美亚洲嫩草| 一卡2卡三卡四卡精品乱码亚洲| 天天一区二区日本电影三级 | 法律面前人人平等表现在哪些方面| 天天一区二区日本电影三级 | 亚洲第一电影网av| av免费在线观看网站| 国产xxxxx性猛交| 日本a在线网址| 国产精品亚洲av一区麻豆| 搡老熟女国产l中国老女人| 国产精品久久久久久亚洲av鲁大| 黄频高清免费视频| 欧美中文综合在线视频| 女性被躁到高潮视频| 精品久久久久久,| 国内毛片毛片毛片毛片毛片| 亚洲精品国产一区二区精华液| 在线十欧美十亚洲十日本专区| 黑人欧美特级aaaaaa片| 国产精品乱码一区二三区的特点 | 国产高清videossex| 在线观看66精品国产| 国产91精品成人一区二区三区| 在线观看免费日韩欧美大片| 日日摸夜夜添夜夜添小说| 纯流量卡能插随身wifi吗| 亚洲精品国产精品久久久不卡| 香蕉国产在线看| 极品人妻少妇av视频| 91在线观看av| 国产亚洲精品第一综合不卡| 亚洲精华国产精华精| 国产成人啪精品午夜网站| 久久国产精品影院| 天堂√8在线中文| 一级作爱视频免费观看| 老司机深夜福利视频在线观看| 无遮挡黄片免费观看| 黄频高清免费视频| 91大片在线观看| 动漫黄色视频在线观看| 国产一区二区三区综合在线观看| 亚洲一区高清亚洲精品| 不卡av一区二区三区| 搡老岳熟女国产| 男女之事视频高清在线观看| 狂野欧美激情性xxxx| 搞女人的毛片| 亚洲成人国产一区在线观看| 在线观看免费午夜福利视频| 91九色精品人成在线观看| 欧美日本亚洲视频在线播放| 久久中文字幕人妻熟女| 可以免费在线观看a视频的电影网站| 一本综合久久免费| 此物有八面人人有两片| 99国产精品一区二区蜜桃av| 精品久久久久久成人av| 久久国产精品影院| 久久久国产精品麻豆| 国产麻豆69| 欧美性长视频在线观看| 欧美老熟妇乱子伦牲交| 亚洲欧美激情综合另类| 亚洲欧美精品综合一区二区三区| 国产精品二区激情视频| 亚洲成a人片在线一区二区| 757午夜福利合集在线观看| 一区二区日韩欧美中文字幕| 久久久国产成人免费| 国产精品野战在线观看| 国产又爽黄色视频| 在线观看免费日韩欧美大片| 色综合婷婷激情| 色在线成人网| 午夜福利在线观看吧| 免费观看精品视频网站| 日本欧美视频一区| 99精品欧美一区二区三区四区| 国产三级黄色录像| 久久久国产成人免费| 久久亚洲真实| 日韩欧美三级三区| 国产av精品麻豆| 国产精品 欧美亚洲| 亚洲精品国产色婷婷电影| 高潮久久久久久久久久久不卡| 久久 成人 亚洲| 在线观看www视频免费| 老司机午夜十八禁免费视频| 99久久精品国产亚洲精品| 最近最新中文字幕大全电影3 | 18禁黄网站禁片午夜丰满| 可以免费在线观看a视频的电影网站| 咕卡用的链子| 精品人妻1区二区| 18禁国产床啪视频网站| 免费在线观看完整版高清| 大陆偷拍与自拍| 18美女黄网站色大片免费观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲欧美激情在线| 无遮挡黄片免费观看| 色老头精品视频在线观看| 色av中文字幕| 中出人妻视频一区二区| 国产片内射在线| 免费在线观看日本一区| 中文字幕av电影在线播放| 手机成人av网站| 每晚都被弄得嗷嗷叫到高潮| 99re在线观看精品视频| 黄片小视频在线播放| 日本免费a在线| 久久午夜亚洲精品久久| 美女 人体艺术 gogo| 91大片在线观看| 欧美性长视频在线观看| ponron亚洲| 久久天堂一区二区三区四区| 免费高清在线观看日韩| 最近最新免费中文字幕在线| 大码成人一级视频| 久久久久久免费高清国产稀缺| 亚洲久久久国产精品| 人人妻人人澡人人看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲片人在线观看| 精品无人区乱码1区二区| 欧美国产精品va在线观看不卡| 欧美亚洲日本最大视频资源| 女性生殖器流出的白浆| 国产精品久久久人人做人人爽| АⅤ资源中文在线天堂| 无限看片的www在线观看| 欧美不卡视频在线免费观看 | 国产成人影院久久av| 桃色一区二区三区在线观看| 久久青草综合色| 麻豆成人av在线观看| 亚洲少妇的诱惑av| 99在线人妻在线中文字幕| 满18在线观看网站| 大香蕉久久成人网| 国产精品亚洲av一区麻豆| 中文字幕av电影在线播放| 麻豆一二三区av精品| 无限看片的www在线观看| 久久青草综合色| 欧美色欧美亚洲另类二区 | 久久热在线av| 好看av亚洲va欧美ⅴa在| 午夜视频精品福利| 他把我摸到了高潮在线观看| 亚洲无线在线观看| 午夜福利视频1000在线观看 | 亚洲自拍偷在线| 午夜福利,免费看| 久久天躁狠狠躁夜夜2o2o| 亚洲一区高清亚洲精品| 看片在线看免费视频| 中文亚洲av片在线观看爽| 精品久久久久久久人妻蜜臀av | 欧美黑人精品巨大| 国产精品自产拍在线观看55亚洲| 高清毛片免费观看视频网站| 老司机午夜福利在线观看视频| 中文字幕精品免费在线观看视频| 黑人巨大精品欧美一区二区蜜桃| 88av欧美| a在线观看视频网站| 午夜久久久在线观看| 久久精品亚洲精品国产色婷小说| 亚洲国产毛片av蜜桃av| 免费看美女性在线毛片视频| 久久久精品国产亚洲av高清涩受| 中文字幕久久专区| 女人被躁到高潮嗷嗷叫费观| 99热只有精品国产| 一区在线观看完整版| 亚洲 欧美一区二区三区| 久久国产乱子伦精品免费另类| 久久久久国内视频| 少妇熟女aⅴ在线视频| 国产精品二区激情视频| 美女扒开内裤让男人捅视频| 午夜免费成人在线视频| 久久九九热精品免费| 757午夜福利合集在线观看| 国产一区二区三区综合在线观看| 欧美成人免费av一区二区三区| 老汉色∧v一级毛片| 最近最新中文字幕大全免费视频| 淫妇啪啪啪对白视频| 亚洲色图av天堂| 中出人妻视频一区二区| 久久伊人香网站| 天堂影院成人在线观看| √禁漫天堂资源中文www| 国产伦一二天堂av在线观看| 一区二区三区高清视频在线| 中文字幕人妻丝袜一区二区| 亚洲欧美一区二区三区黑人| 久久久精品国产亚洲av高清涩受| 精品国产国语对白av| 自线自在国产av| 啦啦啦 在线观看视频| 真人一进一出gif抽搐免费| 国产成人av激情在线播放| 亚洲精品在线美女| 不卡av一区二区三区| 欧美av亚洲av综合av国产av| 久久这里只有精品19| 麻豆成人av在线观看| 日韩欧美免费精品| 日日摸夜夜添夜夜添小说| 亚洲视频免费观看视频| 自拍欧美九色日韩亚洲蝌蚪91| www.自偷自拍.com| 一本久久中文字幕| 国产精品一区二区精品视频观看| 中文字幕另类日韩欧美亚洲嫩草| 久久久久久久久中文| 久久中文看片网| 婷婷丁香在线五月| 99国产极品粉嫩在线观看| 桃红色精品国产亚洲av| 曰老女人黄片| а√天堂www在线а√下载| 波多野结衣巨乳人妻| 一进一出好大好爽视频| 色综合欧美亚洲国产小说| 国产精品综合久久久久久久免费 | 伦理电影免费视频| 亚洲成人精品中文字幕电影| 18禁裸乳无遮挡免费网站照片 | 久久精品人人爽人人爽视色| 69精品国产乱码久久久| 级片在线观看| 麻豆成人av在线观看| 亚洲男人天堂网一区| 91国产中文字幕| 九色国产91popny在线| 亚洲国产高清在线一区二区三 | 男人舔女人的私密视频| 国产成人精品在线电影| 女人被狂操c到高潮| 免费无遮挡裸体视频| 在线观看免费视频网站a站| 久久久久久久精品吃奶| 欧美黄色淫秽网站| 两个人看的免费小视频| 9热在线视频观看99| 亚洲成人国产一区在线观看| 国产欧美日韩综合在线一区二区| 日日夜夜操网爽| 亚洲成人精品中文字幕电影| 亚洲人成伊人成综合网2020| 18禁观看日本| 真人一进一出gif抽搐免费| av免费在线观看网站| 中文字幕人妻丝袜一区二区| 午夜福利成人在线免费观看| 久久人人精品亚洲av| 99精品久久久久人妻精品| 亚洲欧美精品综合一区二区三区| 69av精品久久久久久| 久久精品国产综合久久久| 老汉色av国产亚洲站长工具| 人人澡人人妻人| 亚洲片人在线观看| 亚洲精品国产一区二区精华液| 男人舔女人的私密视频| 久99久视频精品免费| 在线观看www视频免费| 久久久久国内视频| 麻豆成人av在线观看| 久久久久久久久久久久大奶| 国产精品一区二区精品视频观看| 老鸭窝网址在线观看| 成人永久免费在线观看视频| 中文字幕人妻丝袜一区二区| 女人精品久久久久毛片| 国产精品亚洲av一区麻豆| 校园春色视频在线观看| 黄网站色视频无遮挡免费观看| 国产精品九九99| 久久久久国产一级毛片高清牌| 麻豆成人av在线观看| 88av欧美| av超薄肉色丝袜交足视频| 欧美中文综合在线视频| www日本在线高清视频| 午夜福利一区二区在线看| 午夜精品久久久久久毛片777| 午夜福利影视在线免费观看| 国产免费av片在线观看野外av| 99精品久久久久人妻精品| 不卡av一区二区三区| 亚洲,欧美精品.| 国产麻豆成人av免费视频| 久久精品亚洲熟妇少妇任你| 欧美激情 高清一区二区三区| 黄色片一级片一级黄色片| 一边摸一边抽搐一进一出视频| 久久中文字幕一级| 国产成+人综合+亚洲专区| 琪琪午夜伦伦电影理论片6080| 免费观看精品视频网站| 久久久久久国产a免费观看| 精品国内亚洲2022精品成人| 琪琪午夜伦伦电影理论片6080| 久久 成人 亚洲| 欧美激情极品国产一区二区三区| 啦啦啦韩国在线观看视频| 亚洲成av片中文字幕在线观看| 人人澡人人妻人| 国产91精品成人一区二区三区| 中文字幕最新亚洲高清| 在线免费观看的www视频| 两个人看的免费小视频| 亚洲精品国产精品久久久不卡| 窝窝影院91人妻| 亚洲精品美女久久av网站| 午夜免费鲁丝| 国产一区二区在线av高清观看| 久久青草综合色| 亚洲av成人av| 成人欧美大片| 黄色毛片三级朝国网站| 美女高潮喷水抽搐中文字幕| 窝窝影院91人妻| 香蕉丝袜av| 国产极品粉嫩免费观看在线| av视频在线观看入口| 久久久久久久精品吃奶| 国产欧美日韩综合在线一区二区| 老司机在亚洲福利影院| 国产午夜精品久久久久久| 美女扒开内裤让男人捅视频| 一进一出抽搐gif免费好疼| 国产亚洲精品久久久久5区| 午夜免费鲁丝| av在线播放免费不卡| 国产极品粉嫩免费观看在线| 国产精品亚洲美女久久久| 国产精品二区激情视频| 91麻豆av在线| 老司机靠b影院| av网站免费在线观看视频| 一级毛片高清免费大全| АⅤ资源中文在线天堂| 91精品国产国语对白视频| 久久精品人人爽人人爽视色| 好看av亚洲va欧美ⅴa在| 国产精品av久久久久免费| 精品国产国语对白av| 精品久久久久久久毛片微露脸| 一区二区三区国产精品乱码| 国产精品av久久久久免费| 国产伦一二天堂av在线观看| av网站免费在线观看视频| 国产日韩一区二区三区精品不卡| 亚洲激情在线av| 成人18禁在线播放| 亚洲av美国av| 中文字幕av电影在线播放| 亚洲av日韩精品久久久久久密| 午夜日韩欧美国产| 制服丝袜大香蕉在线| 国产亚洲精品综合一区在线观看 | 国产又爽黄色视频| 午夜福利视频1000在线观看 | 纯流量卡能插随身wifi吗| a级毛片在线看网站| 国产1区2区3区精品| 88av欧美| 精品久久久久久成人av| 91字幕亚洲| www.www免费av| 757午夜福利合集在线观看| or卡值多少钱| 最新美女视频免费是黄的| 一a级毛片在线观看| 操美女的视频在线观看| 久久婷婷成人综合色麻豆| 国产精品久久视频播放| 欧美绝顶高潮抽搐喷水| 变态另类丝袜制服| 久久国产精品人妻蜜桃| 在线天堂中文资源库| 真人做人爱边吃奶动态| 嫩草影院精品99| 免费在线观看视频国产中文字幕亚洲| 啪啪无遮挡十八禁网站| 久久精品91蜜桃| 波多野结衣一区麻豆| 99久久久亚洲精品蜜臀av| 日韩有码中文字幕| 色在线成人网| 多毛熟女@视频| 日韩欧美三级三区| 免费一级毛片在线播放高清视频 | 成人国产综合亚洲| 免费搜索国产男女视频| 亚洲一区中文字幕在线| 亚洲人成电影观看| 18禁观看日本| 久久久精品国产亚洲av高清涩受| 999久久久国产精品视频| av超薄肉色丝袜交足视频| 99国产精品一区二区三区| 嫩草影院精品99| 真人做人爱边吃奶动态| 中文字幕另类日韩欧美亚洲嫩草| 国产精品一区二区精品视频观看| 久久午夜亚洲精品久久| 午夜影院日韩av| 亚洲国产中文字幕在线视频| 国产xxxxx性猛交| 中文字幕最新亚洲高清| videosex国产| 亚洲国产精品合色在线| 香蕉国产在线看| 村上凉子中文字幕在线| 婷婷丁香在线五月| 国产精品二区激情视频| 欧美绝顶高潮抽搐喷水| 黄片大片在线免费观看| 亚洲国产欧美日韩在线播放| 狂野欧美激情性xxxx| 欧美午夜高清在线| 成人欧美大片| а√天堂www在线а√下载| 精品少妇一区二区三区视频日本电影| 男人操女人黄网站| 涩涩av久久男人的天堂| √禁漫天堂资源中文www| 亚洲熟妇熟女久久| 亚洲va日本ⅴa欧美va伊人久久| av网站免费在线观看视频| 激情视频va一区二区三区| aaaaa片日本免费| 国产真人三级小视频在线观看| 成熟少妇高潮喷水视频| 国产又色又爽无遮挡免费看| 一级毛片女人18水好多| 国产一区二区激情短视频| 搞女人的毛片| 国产成年人精品一区二区| 真人做人爱边吃奶动态| 欧美激情极品国产一区二区三区| 日日摸夜夜添夜夜添小说| 国产精品精品国产色婷婷| 50天的宝宝边吃奶边哭怎么回事| 国产精品野战在线观看| 欧美激情 高清一区二区三区| 久久影院123| 天堂√8在线中文| 两个人看的免费小视频| 日韩欧美一区视频在线观看| 69精品国产乱码久久久| 亚洲精品在线美女| 禁无遮挡网站| 黄网站色视频无遮挡免费观看| 久久精品91无色码中文字幕| 搞女人的毛片| 久久精品国产亚洲av高清一级| 一级毛片精品| 岛国视频午夜一区免费看| 人人澡人人妻人| 两个人看的免费小视频| 久久久久久免费高清国产稀缺| 亚洲五月婷婷丁香| 国产伦人伦偷精品视频| 欧美大码av| 免费av毛片视频| 狂野欧美激情性xxxx| 精品久久久久久,| 久久久久久亚洲精品国产蜜桃av| 日韩精品中文字幕看吧| 最近最新免费中文字幕在线| 一边摸一边抽搐一进一出视频| 亚洲国产精品成人综合色| 在线十欧美十亚洲十日本专区| 国产成人精品久久二区二区91| 欧美激情久久久久久爽电影 | 久久这里只有精品19| 国产精品综合久久久久久久免费 | 电影成人av| 麻豆一二三区av精品| 国产亚洲欧美在线一区二区| xxx96com| 99精品欧美一区二区三区四区| 9色porny在线观看| 男女床上黄色一级片免费看| 国产成人av激情在线播放| 啦啦啦 在线观看视频| 国产成人免费无遮挡视频| 成人亚洲精品一区在线观看| 国产亚洲精品综合一区在线观看 | 国产熟女xx| 久久久久久久午夜电影| 免费在线观看日本一区| 99国产精品免费福利视频| 日韩欧美国产在线观看| 午夜福利一区二区在线看| 成人18禁在线播放| 深夜精品福利| 国产精品久久视频播放| 99香蕉大伊视频| 18禁观看日本| 日本撒尿小便嘘嘘汇集6| 999久久久国产精品视频| 午夜老司机福利片| 国产av在哪里看| 精品无人区乱码1区二区| 村上凉子中文字幕在线| 大码成人一级视频| 精品人妻在线不人妻| 最新美女视频免费是黄的| 黄片小视频在线播放| 99国产精品99久久久久| 国产精品av久久久久免费| 午夜激情av网站| 性少妇av在线| 欧美激情 高清一区二区三区| 午夜福利影视在线免费观看| 两性夫妻黄色片| 国产亚洲精品第一综合不卡| 夜夜躁狠狠躁天天躁| 国产精品 国内视频| 午夜福利18| 丁香六月欧美| a在线观看视频网站| 午夜福利欧美成人| 久久精品国产综合久久久| 麻豆国产av国片精品| 久久欧美精品欧美久久欧美| www.精华液|