• <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
    第三講 棉花肥害診斷及其防治
    此物有八面人人有两片| 日韩欧美在线二视频| 欧美色欧美亚洲另类二区| 国产日本99.免费观看| 久久这里只有精品19| 1024手机看黄色片| 99久久国产精品久久久| 美女 人体艺术 gogo| 国产精品亚洲av一区麻豆| 亚洲欧美精品综合久久99| 久久草成人影院| 日韩欧美在线二视频| 日日爽夜夜爽网站| 女性被躁到高潮视频| 午夜免费成人在线视频| 亚洲无线在线观看| 两人在一起打扑克的视频| 一级毛片精品| 天天躁夜夜躁狠狠躁躁| 黄色视频不卡| 久久久精品国产亚洲av高清涩受| 午夜日韩欧美国产| 国产野战对白在线观看| 一夜夜www| 99久久精品国产亚洲精品| 一级毛片精品| 午夜影院日韩av| 久久久国产精品麻豆| 麻豆成人午夜福利视频| 看免费av毛片| 亚洲精品国产区一区二| 国产一区二区在线av高清观看| 欧美成狂野欧美在线观看| 国产久久久一区二区三区| 女性被躁到高潮视频| 巨乳人妻的诱惑在线观看| 丝袜美腿诱惑在线| 精品午夜福利视频在线观看一区| 成年版毛片免费区| 免费av毛片视频| 久久久久国内视频| 淫妇啪啪啪对白视频| 亚洲中文字幕日韩| АⅤ资源中文在线天堂| 午夜老司机福利片| 欧美乱色亚洲激情| 久久亚洲精品不卡| 一区二区三区精品91| 最近最新中文字幕大全电影3 | 视频区欧美日本亚洲| 麻豆国产av国片精品| 精品一区二区三区av网在线观看| 精品电影一区二区在线| 色综合站精品国产| 亚洲免费av在线视频| 91在线观看av| 18禁黄网站禁片免费观看直播| 国产一区二区三区在线臀色熟女| 美女 人体艺术 gogo| 热re99久久国产66热| 91字幕亚洲| av有码第一页| 国产色视频综合| 最近最新中文字幕大全电影3 | 大型av网站在线播放| 成人亚洲精品av一区二区| 欧美国产精品va在线观看不卡| 久久国产乱子伦精品免费另类| 欧美性猛交╳xxx乱大交人| 亚洲人成伊人成综合网2020| 老司机深夜福利视频在线观看| 最好的美女福利视频网| 99国产精品一区二区蜜桃av| 午夜久久久久精精品| 国产亚洲精品久久久久久毛片| 国产蜜桃级精品一区二区三区| 久久精品国产亚洲av高清一级| 成人三级做爰电影| 中文亚洲av片在线观看爽| 国产私拍福利视频在线观看| 在线观看日韩欧美| 国产亚洲精品一区二区www| 亚洲美女黄片视频| 国产精品 国内视频| 精品久久久久久久久久久久久 | 久久精品国产亚洲av高清一级| 特大巨黑吊av在线直播 | 精品欧美国产一区二区三| 看免费av毛片| 91成人精品电影| 国内毛片毛片毛片毛片毛片| 精品久久久久久成人av| 一边摸一边做爽爽视频免费| 亚洲精品美女久久久久99蜜臀| 1024手机看黄色片| 香蕉久久夜色| 欧美在线黄色| 午夜免费激情av| 国产成人欧美在线观看| 婷婷亚洲欧美| 男女床上黄色一级片免费看| 亚洲人成伊人成综合网2020| 国产精品国产高清国产av| 亚洲片人在线观看| 好男人在线观看高清免费视频 | 欧美zozozo另类| 国产久久久一区二区三区| 一边摸一边做爽爽视频免费| 日韩三级视频一区二区三区| 欧美色欧美亚洲另类二区| 两性午夜刺激爽爽歪歪视频在线观看 | 国产一级毛片七仙女欲春2 | 少妇熟女aⅴ在线视频| 国产高清激情床上av| 日韩高清综合在线| 真人做人爱边吃奶动态| 久久人人精品亚洲av| 亚洲av第一区精品v没综合| 亚洲男人的天堂狠狠| 亚洲五月色婷婷综合| 天堂动漫精品| 国产单亲对白刺激| 一区二区三区激情视频| 听说在线观看完整版免费高清| 成人国语在线视频| 国产精品国产高清国产av| 亚洲激情在线av| 色播亚洲综合网| 国产精品 欧美亚洲| 久久香蕉精品热| 国产欧美日韩精品亚洲av| 日本一区二区免费在线视频| 欧美激情久久久久久爽电影| 中文字幕人妻丝袜一区二区| 欧美一级毛片孕妇| 久久婷婷成人综合色麻豆| 亚洲av日韩精品久久久久久密| 首页视频小说图片口味搜索| 亚洲人成网站在线播放欧美日韩| 欧美日韩中文字幕国产精品一区二区三区| 国产欧美日韩一区二区精品| 香蕉久久夜色| 免费av毛片视频| 麻豆国产av国片精品| 9191精品国产免费久久| 他把我摸到了高潮在线观看| 久久久精品国产亚洲av高清涩受| 婷婷亚洲欧美| 国产亚洲欧美在线一区二区| 成人18禁在线播放| 黄色 视频免费看| 久久久精品欧美日韩精品| 特大巨黑吊av在线直播 | 亚洲人成伊人成综合网2020| 女同久久另类99精品国产91| 国产一区二区三区在线臀色熟女| 最新美女视频免费是黄的| 少妇粗大呻吟视频| 国产高清有码在线观看视频 | 成人一区二区视频在线观看| 久久狼人影院| 欧美乱妇无乱码| 色综合婷婷激情| www日本黄色视频网| 国产精品久久视频播放| 日日爽夜夜爽网站| 日韩欧美一区二区三区在线观看| 久久这里只有精品19| 国产精品国产高清国产av| 深夜精品福利| 精品午夜福利视频在线观看一区| 日韩免费av在线播放| 国产伦人伦偷精品视频| 国产aⅴ精品一区二区三区波| 91字幕亚洲| 最近最新中文字幕大全免费视频| 窝窝影院91人妻| 国产三级在线视频| 真人一进一出gif抽搐免费| 国产亚洲精品综合一区在线观看 | 一个人观看的视频www高清免费观看 | 欧美性猛交黑人性爽| 99久久综合精品五月天人人| ponron亚洲| 最近最新中文字幕大全免费视频| 男男h啪啪无遮挡| 99re在线观看精品视频| 最近在线观看免费完整版| av电影中文网址| 他把我摸到了高潮在线观看| 狂野欧美激情性xxxx| 午夜两性在线视频| 日日干狠狠操夜夜爽| 亚洲在线自拍视频| 国产高清激情床上av| 亚洲人成电影免费在线| 精品少妇一区二区三区视频日本电影| 俄罗斯特黄特色一大片| 欧美一级a爱片免费观看看 | 欧美成人一区二区免费高清观看 | 国产又黄又爽又无遮挡在线| 日韩精品青青久久久久久| a级毛片在线看网站| 国产99久久九九免费精品| 国产成人系列免费观看| 免费在线观看黄色视频的| 曰老女人黄片| 久久人妻福利社区极品人妻图片| 国产成人av激情在线播放| 亚洲美女黄片视频| 国产极品粉嫩免费观看在线| 国产亚洲精品久久久久久毛片| 黄色视频,在线免费观看| 级片在线观看| 最好的美女福利视频网| 欧美日韩福利视频一区二区| 又紧又爽又黄一区二区| 亚洲精品中文字幕一二三四区| 黄片播放在线免费| av欧美777| 国产精品美女特级片免费视频播放器 | 成人国语在线视频| 亚洲全国av大片| 中亚洲国语对白在线视频| 看片在线看免费视频| 男女下面进入的视频免费午夜 | 国产亚洲精品久久久久久毛片| 成在线人永久免费视频| 高潮久久久久久久久久久不卡| 国产又爽黄色视频| 亚洲成a人片在线一区二区| 中文亚洲av片在线观看爽| 国产亚洲精品久久久久5区| 91老司机精品| 午夜福利一区二区在线看| 亚洲成人久久爱视频| 国产精品免费视频内射| 亚洲欧美一区二区三区黑人| 亚洲av片天天在线观看| 18禁国产床啪视频网站| 色av中文字幕| 国产三级在线视频| 美女 人体艺术 gogo| 亚洲人成网站在线播放欧美日韩| 免费观看人在逋| 人妻丰满熟妇av一区二区三区| 18禁黄网站禁片免费观看直播| 一级作爱视频免费观看| 久久精品91蜜桃| 日日干狠狠操夜夜爽| 久久 成人 亚洲| 少妇被粗大的猛进出69影院| 一卡2卡三卡四卡精品乱码亚洲| 好男人电影高清在线观看| 18美女黄网站色大片免费观看| 亚洲国产日韩欧美精品在线观看 | 久久久久久久久中文| 亚洲第一青青草原| 法律面前人人平等表现在哪些方面| 嫩草影视91久久| 99riav亚洲国产免费| 国产精品亚洲美女久久久| 亚洲精品美女久久av网站| 精品日产1卡2卡| 日本免费一区二区三区高清不卡| 国产精华一区二区三区| 精品少妇一区二区三区视频日本电影| 国产成人精品久久二区二区91| 一区二区三区国产精品乱码| 一二三四社区在线视频社区8| 午夜福利成人在线免费观看| 国产在线观看jvid| 97超级碰碰碰精品色视频在线观看| 久久久久久免费高清国产稀缺| 久久精品国产亚洲av高清一级| 亚洲中文日韩欧美视频| 国产av不卡久久| 成人国产综合亚洲| 视频区欧美日本亚洲| 久久国产精品人妻蜜桃| 18禁黄网站禁片免费观看直播| 国产成人一区二区三区免费视频网站| 亚洲国产高清在线一区二区三 | 老司机午夜福利在线观看视频| 熟女电影av网| 亚洲精华国产精华精| 国产免费男女视频| 国产精品爽爽va在线观看网站 | av在线播放免费不卡| 欧美色欧美亚洲另类二区| 亚洲成人久久性| 久久精品成人免费网站| 激情在线观看视频在线高清| 国产黄色小视频在线观看| 亚洲中文日韩欧美视频| 女性生殖器流出的白浆| 免费一级毛片在线播放高清视频| 国产黄a三级三级三级人| 黄色丝袜av网址大全| 国产成人系列免费观看| 好男人在线观看高清免费视频 | 黄色片一级片一级黄色片| 欧美三级亚洲精品| 免费看十八禁软件| 亚洲专区国产一区二区| 中亚洲国语对白在线视频| 亚洲va日本ⅴa欧美va伊人久久| a级毛片a级免费在线| 高潮久久久久久久久久久不卡| 精品国产美女av久久久久小说| 人妻久久中文字幕网| 人人妻,人人澡人人爽秒播| 1024手机看黄色片| 色哟哟哟哟哟哟| 在线观看午夜福利视频| 一级a爱视频在线免费观看| 精品免费久久久久久久清纯| 99国产精品一区二区三区| 女性生殖器流出的白浆| 99国产精品99久久久久| 亚洲av日韩精品久久久久久密| 欧美一级a爱片免费观看看 | 一级作爱视频免费观看| 99久久久亚洲精品蜜臀av| 国产色视频综合| 女同久久另类99精品国产91| 啪啪无遮挡十八禁网站| 国产精品一区二区精品视频观看| 母亲3免费完整高清在线观看| 欧美乱色亚洲激情| 精品免费久久久久久久清纯| 欧美国产精品va在线观看不卡| 国产精品精品国产色婷婷| 国产精品野战在线观看| av有码第一页| 热re99久久国产66热| 变态另类丝袜制服| 久久国产亚洲av麻豆专区| 国产成年人精品一区二区| 91大片在线观看| 丰满的人妻完整版| 国产高清videossex| 久久精品夜夜夜夜夜久久蜜豆 | 欧美久久黑人一区二区| 18禁国产床啪视频网站| 黄片播放在线免费| 天堂√8在线中文| 国产成人av激情在线播放| 亚洲第一青青草原| 亚洲欧美日韩高清在线视频| 在线天堂中文资源库| 国产精品一区二区免费欧美| www.999成人在线观看| 中文字幕最新亚洲高清| 欧美激情 高清一区二区三区| 正在播放国产对白刺激| 国产私拍福利视频在线观看| 亚洲在线自拍视频| 黄片小视频在线播放| 久久伊人香网站| 日韩精品免费视频一区二区三区| 欧美国产日韩亚洲一区| 免费电影在线观看免费观看| 两个人看的免费小视频| 国产91精品成人一区二区三区| 久久久久精品国产欧美久久久| 国产精品香港三级国产av潘金莲| 欧美国产日韩亚洲一区| 91麻豆精品激情在线观看国产| 亚洲欧洲精品一区二区精品久久久| 草草在线视频免费看| 日韩欧美在线二视频| 久久亚洲精品不卡| 一本精品99久久精品77| 天堂动漫精品| 久久国产亚洲av麻豆专区| 美女高潮到喷水免费观看| 啦啦啦韩国在线观看视频| 99在线人妻在线中文字幕| 成人一区二区视频在线观看| 欧美性猛交黑人性爽| 人人妻人人看人人澡| 国产精品爽爽va在线观看网站 | 国产精品久久久久久精品电影 | 91成年电影在线观看| 男男h啪啪无遮挡| 国产精品爽爽va在线观看网站 | 长腿黑丝高跟| 久久 成人 亚洲| 91成年电影在线观看| 99久久99久久久精品蜜桃| 欧美性猛交╳xxx乱大交人| 亚洲一区二区三区不卡视频| 欧美成人性av电影在线观看| 黄色成人免费大全| 成人手机av| 又紧又爽又黄一区二区| 国产黄色小视频在线观看| 又黄又粗又硬又大视频| 欧美性猛交╳xxx乱大交人| 在线观看舔阴道视频| 欧美性猛交╳xxx乱大交人| 午夜日韩欧美国产| 久久人妻福利社区极品人妻图片| 婷婷亚洲欧美| 桃红色精品国产亚洲av| 每晚都被弄得嗷嗷叫到高潮| 日韩一卡2卡3卡4卡2021年| 神马国产精品三级电影在线观看 | 精品国产国语对白av| 免费看美女性在线毛片视频| 熟妇人妻久久中文字幕3abv| 狠狠狠狠99中文字幕| 国产精品国产高清国产av| 母亲3免费完整高清在线观看| 少妇 在线观看| 午夜免费激情av| 老汉色∧v一级毛片| 久99久视频精品免费| 97碰自拍视频| 国产v大片淫在线免费观看| 久久狼人影院| 欧美久久黑人一区二区| 国产精品日韩av在线免费观看| 看黄色毛片网站| 黄色视频,在线免费观看| 欧美 亚洲 国产 日韩一| 精品久久久久久久毛片微露脸| 成人永久免费在线观看视频| 夜夜爽天天搞| 亚洲一区二区三区色噜噜| 国产精品美女特级片免费视频播放器 | 黄色片一级片一级黄色片| 麻豆久久精品国产亚洲av| 可以在线观看的亚洲视频| 人人妻,人人澡人人爽秒播| 在线av久久热| 黄色女人牲交| 99精品在免费线老司机午夜| 久久中文字幕人妻熟女| 色老头精品视频在线观看| 夜夜爽天天搞| 亚洲精品中文字幕在线视频| xxxwww97欧美| 久久狼人影院| 淫妇啪啪啪对白视频| 久久久国产精品麻豆| 欧美乱色亚洲激情| 视频在线观看一区二区三区| 不卡一级毛片| e午夜精品久久久久久久| 一a级毛片在线观看| 中文字幕最新亚洲高清| 久久久久久国产a免费观看| 亚洲一区二区三区不卡视频| 亚洲国产欧美一区二区综合| 久久久久久免费高清国产稀缺| 人人妻人人澡人人看| 精品久久久久久,| 熟女少妇亚洲综合色aaa.| 一进一出好大好爽视频| 看片在线看免费视频| 日本免费a在线| 国产精品一区二区精品视频观看| 最近在线观看免费完整版| 亚洲精品在线观看二区| 欧美另类亚洲清纯唯美| 亚洲成av片中文字幕在线观看| 色综合婷婷激情| 最近最新中文字幕大全免费视频| 亚洲三区欧美一区| 在线观看免费午夜福利视频| 日韩欧美国产一区二区入口| 99热这里只有精品一区 | 一区二区三区激情视频| 精品无人区乱码1区二区| 美女大奶头视频| 日韩欧美国产在线观看| videosex国产| 激情在线观看视频在线高清| 制服诱惑二区| 男女午夜视频在线观看| 亚洲第一电影网av| 女性生殖器流出的白浆| 激情在线观看视频在线高清| 久9热在线精品视频| 国产男靠女视频免费网站| www.精华液| 十八禁网站免费在线| 亚洲午夜理论影院| 老鸭窝网址在线观看| 国产99白浆流出| xxxwww97欧美| 精品国产美女av久久久久小说| 成人欧美大片| 欧美zozozo另类| 国产一区二区三区视频了| 很黄的视频免费| 午夜福利欧美成人| 亚洲欧美日韩高清在线视频| 欧美日韩乱码在线| 精品国产亚洲在线| 在线av久久热| 亚洲av美国av| av电影中文网址| www.www免费av| 老熟妇仑乱视频hdxx| 久久久久国产一级毛片高清牌| 欧美激情 高清一区二区三区| 在线国产一区二区在线| 国产精品久久久久久人妻精品电影| 动漫黄色视频在线观看| 国产精品亚洲av一区麻豆| 午夜免费观看网址| 久久精品国产综合久久久| 最近最新中文字幕大全免费视频| 亚洲一区二区三区色噜噜| 国产精品亚洲av一区麻豆| 久久中文看片网| 在线国产一区二区在线| 1024视频免费在线观看| 国产精品久久久久久亚洲av鲁大| 亚洲第一电影网av| a级毛片a级免费在线| 亚洲精品一卡2卡三卡4卡5卡| 国产极品粉嫩免费观看在线| 久久亚洲精品不卡| 中文字幕高清在线视频| 国产精品久久视频播放| 久久欧美精品欧美久久欧美| 国产片内射在线| 欧美性猛交黑人性爽| av免费在线观看网站| 人人妻人人澡人人看| 久久精品aⅴ一区二区三区四区| 日日摸夜夜添夜夜添小说| 亚洲一码二码三码区别大吗| 国产精品综合久久久久久久免费| 校园春色视频在线观看| 琪琪午夜伦伦电影理论片6080| 色老头精品视频在线观看| 一级作爱视频免费观看| 一区福利在线观看| 中文字幕高清在线视频| 视频在线观看一区二区三区| 国产男靠女视频免费网站| 午夜久久久久精精品| 叶爱在线成人免费视频播放| 在线观看www视频免费| 国产单亲对白刺激| 久久精品影院6| 婷婷亚洲欧美| 人人妻人人澡人人看| 少妇被粗大的猛进出69影院| 亚洲国产日韩欧美精品在线观看 | 波多野结衣高清无吗| 黄色丝袜av网址大全| 国产精品自产拍在线观看55亚洲| 亚洲中文字幕日韩| 精品久久蜜臀av无| 在线天堂中文资源库| 欧美日韩精品网址| 免费在线观看亚洲国产| 香蕉av资源在线| 国产成+人综合+亚洲专区| 亚洲成av人片免费观看| 97人妻精品一区二区三区麻豆 | 国产一卡二卡三卡精品| 欧美乱色亚洲激情| 麻豆成人午夜福利视频| 久久精品aⅴ一区二区三区四区| 国产成人精品久久二区二区免费| 国产成人影院久久av| 亚洲美女黄片视频| 久久青草综合色| 免费搜索国产男女视频| 神马国产精品三级电影在线观看 | 国产午夜福利久久久久久| 精品欧美一区二区三区在线| 久久香蕉国产精品| 一二三四在线观看免费中文在| 两人在一起打扑克的视频| 天堂动漫精品| 国产乱人伦免费视频| 99久久精品国产亚洲精品| 国产成年人精品一区二区| 成人国产综合亚洲| 国内精品久久久久久久电影| 国产91精品成人一区二区三区| 搞女人的毛片| 99精品久久久久人妻精品| 国产主播在线观看一区二区| 亚洲午夜理论影院| 午夜精品在线福利| 国产一区在线观看成人免费| 热re99久久国产66热| 在线观看免费日韩欧美大片| 中文字幕人成人乱码亚洲影| 国产成人精品久久二区二区免费| 在线观看免费日韩欧美大片| 午夜福利欧美成人| 9191精品国产免费久久| 一本一本综合久久| 久热这里只有精品99| 91在线观看av| 麻豆国产av国片精品| 国产精品久久久人人做人人爽| 欧美乱码精品一区二区三区| 国产色视频综合| 午夜a级毛片| 韩国av一区二区三区四区| 麻豆国产av国片精品| 日韩欧美一区二区三区在线观看| 国产免费男女视频| 亚洲成人久久爱视频|