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

    MNase-seq與核小體定占位研究

    2020-12-18 10:43:20鄧瑋杭李鑫輝
    遺傳 2020年12期
    關(guān)鍵詞:染色質(zhì)基因組位點

    鄧瑋杭,李鑫輝

    MNase-seq與核小體定占位研究

    鄧瑋杭,李鑫輝

    上海交通大學(xué)生物醫(yī)學(xué)工程學(xué)院,上海 200240

    核小體是染色質(zhì)復(fù)雜三維結(jié)構(gòu)的基本單位,它在基因組上的定位及占位在DNA轉(zhuǎn)錄、復(fù)制和修復(fù)等基礎(chǔ)生物過程中發(fā)揮重要功能。在眾多核小體定占位研究技術(shù)中,微球菌核酸酶測序(micrococcal nuclease sequencing, MNase-seq)被認(rèn)為是目前最為高效的方法,因此應(yīng)用十分廣泛。研究人員利用該技術(shù)繪制了多種生物的核小體圖譜,并揭示了核小體組織特點的共性和差異。本文介紹了MNase-seq的技術(shù)原理以及在解析核小體組織及其功能中的應(yīng)用,總結(jié)了在染色質(zhì)構(gòu)象這一快速發(fā)展領(lǐng)域中的研究進展,并展望了染色質(zhì)生物學(xué)的未來發(fā)展方向。由MNase-seq揭示的核小體組織結(jié)構(gòu)為基因表達(dá)和發(fā)育調(diào)控提供了新的見解,也有助于人們理解疾病的發(fā)生過程。

    核小體;染色質(zhì)結(jié)構(gòu);染色質(zhì)重塑;下一代測序技術(shù);微球菌核酸酶

    構(gòu)建人類細(xì)胞的一個關(guān)鍵步驟是將近2 m的DNA組裝成染色質(zhì)并包裝進直徑10 μm的細(xì)胞核內(nèi),而核小體作為染色質(zhì)組裝的基本單位,一般是由長度為145~147 bp的DNA以左螺旋的方式纏繞組蛋白八聚體接近兩圈構(gòu)成的[1,2]。其中組蛋白八聚體由2個拷貝的4種核心組蛋白H2A、H2B、H3和H4組成,H1作為連接組蛋白。在基因組的大部分區(qū)域都有核小體覆蓋,約75%~90%的基因組DNA包裹在核小體中[3]。核小體定位、占位以及組蛋白修飾等影響著啟動子、轉(zhuǎn)錄起始位點的建立以及染色質(zhì)高級結(jié)構(gòu)的組裝,參與基因轉(zhuǎn)錄、DNA復(fù)制和修復(fù)過程。核小體動態(tài)組織通過改變?nèi)旧|(zhì)結(jié)構(gòu)的方式影響細(xì)胞功能,因此與生物體的基因調(diào)控、發(fā)育分化以及細(xì)胞應(yīng)激過程密切相關(guān)。

    應(yīng)用MNase酶切配合下一代測序技術(shù),人們加快了對核小體組織及染色質(zhì)結(jié)構(gòu)的研究,有望更深入了解基因調(diào)控模式以及染色質(zhì)功能。本文主要介紹了微球菌核酸酶測序(micrococcal nuclease seque-ncing, MNase-seq)的技術(shù)原理和相關(guān)衍生技術(shù),總結(jié)了MNase-seq等技術(shù)在揭示基因組上核小體的組織特點中的作用,綜述了近年來多種生物細(xì)胞內(nèi)核小體圖譜的研究進展,以期為今后以核小體組織為核心的染色質(zhì)構(gòu)象等研究提供參考。

    1 MNase-seq技術(shù)原理及其發(fā)展

    1.1 MNase-seq技術(shù)原理

    利用微球菌核酸酶切割染色質(zhì)纖維,回收DNA并配合下一代測序技術(shù)來繪制核小體定位圖譜,稱作MNase-seq。盡管MNase-seq在近10年來才得以飛速發(fā)展,但早在20世紀(jì)70年代,研究人員就開始利用MNase消化染色質(zhì)并研究其結(jié)構(gòu)[4,5]。MNase來源于金黃色葡萄球菌(),同時具有核酸內(nèi)切酶與外切酶活性。核小體間的連接DNA(linker DNA)對MNase的敏感性要比核小體DNA高25倍[6],根據(jù)酶的這一特性,在對染色質(zhì)DNA進行充分消化后,大部分核小體間的裸露DNA都被消化,由此可以富集到單核小體DNA片段。

    MNase-seq的主要步驟包括:(1)提取細(xì)胞核;(2)用MNase對染色質(zhì)進行酶切;(3)終止酶切、去除RNA及蛋白雜質(zhì);(4)分離提純DNA片段;(5)文庫構(gòu)建;(6)上機測序。在MNase-seq的實驗設(shè)計中,染色質(zhì)交聯(lián)程度以及MNase的酶切水平很關(guān)鍵。不少實驗利用甲醛做交聯(lián)劑(crosslinking agent),在細(xì)胞活體條件下對蛋白質(zhì)-DNA相互作用進行固定,以免在染色質(zhì)制備過程中丟失這些相互作用及破壞核小體組織結(jié)構(gòu)[7]。也有研究者用MNase直接處理未固定的染色質(zhì)[8~10],該方法減少了組蛋白以外的蛋白質(zhì)-DNA相互作用帶來的實驗誤差。盡管該方法可能導(dǎo)致核小體結(jié)構(gòu)的改變,但研究發(fā)現(xiàn)是否使用甲醛固定對核小體的組織改變不大[11]。使用中等或較高水平的酶切使得酶切產(chǎn)物中80%~100%為單核小體片段,有利于減少建庫過程中片段篩選帶來的額外技術(shù)誤差。MNase酶切后進行電泳分析,150 bp左右的條帶則是單個的完整核小體DNA。

    MNase-seq的優(yōu)點在于技術(shù)難度較低,具有較高的分辨率,且數(shù)據(jù)處理相對簡單。較高的分辨率得益于MNase的酶切特性,MNase處理染色質(zhì)可以高效去除連接DNA,得到的DNA片段末端正是包裹組蛋白的DNA兩端。相比用超聲打斷DNA的常規(guī)方法,MNase處理DNA可以獲得長度較均一的DNA片段,從而得到核小體更為精確的位置坐標(biāo)。然而,由于不同DNA片段對MNase酶切的敏感性不同,這一酶切效率的差異使得MNase-seq測序片段的末端不能準(zhǔn)確反映核小體邊緣的位置,因此測序數(shù)據(jù)的處理分析尤為關(guān)鍵。MNase-seq數(shù)據(jù)處理過程主要為數(shù)據(jù)預(yù)處理及質(zhì)量控制、序列比對、核小體定位與占位分析以及數(shù)據(jù)的可視化[12]。在序列比對后,研究者們先后運用多種算法(如iNPS[13]、DiNuP[14]等)分析MNase-seq測序數(shù)據(jù),解析全基因組核小體定位精確圖譜或?qū)Σ町惢嘉缓诵◇w進行分析。其中,Chen等[13]建立了iNPS算法,該算法在NPS(nucleosome positioning from sequencing)的基礎(chǔ)上增加了“核小體邊界信號調(diào)整”與“相鄰核小體合并或分離”步驟,它比通用的NPS算法識別核小體邊界信號的能力更強,因此可以多檢測到約60%的核小體。該算法具有更高的檢測準(zhǔn)確性和穩(wěn)健性,因此有利于下游數(shù)據(jù)分析。核小體定位的理論預(yù)測同樣具有較好的研究前景,研究者們利用核小體DNA/連接DNA的序列特性、堿基二聯(lián)體周期信號等建立數(shù)學(xué)模型(如Segal模型[3]、N-Score模型[15]等)對核小體定位進行預(yù)測。近年來,MNase-seq測序數(shù)據(jù)分析方法的進步使得該技術(shù)日趨成熟,模型的改進與優(yōu)化使得人們對核小體定位預(yù)測的準(zhǔn)確性不斷提高,MNase-seq技術(shù)目前已廣泛應(yīng)用于各種研究場景。

    MNase-seq具有的這些優(yōu)點使其成為檢測核小體分布的優(yōu)越的方法。目前,MNase-seq已被應(yīng)用于釀酒酵母()、果蠅()、人類()以及多種模式生物體內(nèi)的染色質(zhì)結(jié)構(gòu)研究。另一方面,通過MNase-seq測量染色質(zhì)被核小體及其他調(diào)控因子的占有水平,可以間接揭示染色質(zhì)可及性、發(fā)現(xiàn)潛在基因調(diào)控位點[7,12]。

    1.2 MNase-seq衍生技術(shù)及相關(guān)技術(shù)

    目前,以MNase-seq技術(shù)為基礎(chǔ),發(fā)展出了一系列衍生技術(shù)(圖1),例如(1) MNase-ChIP-seq[16],被用于對特定調(diào)控因子、組蛋白修飾或變體的檢測;(2) MNase-Exo-seq,在酶切體系中加入核酸外切酶Ⅲ來彌補MNase的外切酶活性,從而更高程度地切割到核小體核心區(qū)域(core particle),以獲得更精準(zhǔn)的核小體定位[17];(3) MACC-seq (MNase accessibility sequencing),同時對整個基因組上的核小體位置及

    其可及性進行測量,探究核小體占位與染色質(zhì)可及性的關(guān)系[18];(4) MH-seq,識別MNase高敏位點(MNase hypersensitive sites, MHSs),利用MNase可以檢測到DNase I或Tn5無法訪問到的開放染色質(zhì)區(qū)域[19];(5) Array-seq,使用低濃度MNase酶切,用于檢測核小體陣列的規(guī)律性(regularity)及核小體間距(nucleosome spacing)[20];(6) CUT&RUN (cleavage under targets and release using nuclease)技術(shù),MNase被靶向作用于特定的組蛋白修飾位點或轉(zhuǎn)錄因子等蛋白結(jié)合位點[21]。

    除MNase-seq及其衍生技術(shù)以外,染色質(zhì)免疫共沉淀測序技術(shù)(chromatin immunoprecipitation se-quencing, ChIP-seq)、染色質(zhì)開放性測序技術(shù)(assay for transposase-accessible chromatin with high throu-ghput sequencing, ATAC-seq)、DNase I超敏感位點測序(DNase I hypersensitive site sequencing, DNase-seq)以及核小體占位及甲基化測序(nucleosome occupancy and methylome sequencing, NOMe-seq)等技術(shù)也在解析核小體定位及染色質(zhì)結(jié)構(gòu)及功能的研究中起到重要作用。表1對這些技術(shù)進行了總結(jié)比較。

    1.3 少量細(xì)胞及單細(xì)胞MNase-seq技術(shù)

    多年來,大量細(xì)胞MNase-seq技術(shù)研究得較為充分并得到了廣泛應(yīng)用,但該方法存在兩個明顯的缺陷:首先,在研究臨床樣本等珍貴樣本時,收集足夠多的細(xì)胞存在困難;另外,該方法只能得到細(xì)胞群體中核小體定位的平均水平,而盡管在同源細(xì)胞群體中,細(xì)胞間的染色質(zhì)狀態(tài)仍存在明顯的異質(zhì)性,因此不能反映單個細(xì)胞內(nèi)真實的核小體定位狀態(tài)。隨著下一代測序技術(shù)的不斷進步以及單細(xì)胞測序技術(shù)的發(fā)展,MNase-seq技術(shù)也更新迭代,少量細(xì)胞起始的MNase-seq及其衍生技術(shù)、單細(xì)胞MNase-seq技術(shù)[35](single-cell MNase-seq, scMNase- seq)日趨成熟,在近年來得到重要突破。少量細(xì)胞甚至單細(xì)胞起始的MNase-seq技術(shù)有助于研究胚胎發(fā)育或疾病發(fā)展中至關(guān)重要但稀少珍貴的細(xì)胞。

    圖1 MNase-seq及其衍生技術(shù)

    A: MNase-seq; B: Array-seq; C: MNase-ChIP-seq; D: CUT&RUN。

    表1 研究核小體、染色質(zhì)結(jié)構(gòu)的常用技術(shù)

    Mattia等[36]發(fā)展了少量細(xì)胞MACC-seq,最低50個細(xì)胞的投入量即可達(dá)到和大量細(xì)胞MACC相似的可信結(jié)果,擴大了MACC技術(shù)的應(yīng)用范圍。2018年,Peter等[37]改進了CUT&RUN技術(shù),利用該技術(shù)方法對組蛋白修飾進行測定,最低只需要投入約100個細(xì)胞。2019年,Sarah等[38]實現(xiàn)了在單個細(xì)胞中利用uliCUT&RUN (Ultra-low input CUT&RUN)技術(shù)對轉(zhuǎn)錄因子結(jié)合位點進行分析,揭示了轉(zhuǎn)錄因子在胚胎干細(xì)胞中結(jié)合位點的多樣性。2018年,Lai等[39]開發(fā)了scMNase-seq技術(shù),實驗流程主要包括:(1)流式分選收集單個細(xì)胞;(2)裂解細(xì)胞、MNase酶切消化染色質(zhì);(3)分離提純DNA;(4)酶切末端補平、連接測序接頭、PCR擴增約25個循環(huán);(6)篩選合適長度的片段(150~320 bp)上機測序。單個細(xì)胞含有的DNA量極少,因此scMNase-seq的技術(shù)要點在于最小化DNA損失。為達(dá)到這一目標(biāo),實驗細(xì)節(jié)的優(yōu)化包括在DNA提純過程中加入環(huán)狀質(zhì)粒使得DNA顆粒容易可見,以及在建庫流程中不預(yù)先分離單核小體長度的DNA片段。他們利用該技術(shù)同時對NIH3T3細(xì)胞、小鼠()胚胎干細(xì)胞及小鼠幼稚CD4+幼稚T細(xì)胞的全基因組核小體定位及染色質(zhì)可及性進行研究。相比于大量細(xì)胞MNase-seq對核小體間距的研究只限于核小體定位精準(zhǔn)的區(qū)域(例如TSS附近),scMNase-seq技術(shù)可以對單個細(xì)胞全基因組范圍的核小體間距模式進行測量[40]。該研究表明,在異染色質(zhì)區(qū)域,相鄰核小體以約180 bp的間距規(guī)律排列,但在不同細(xì)胞中核小體的定位不一致;在活躍染色質(zhì)區(qū)域,相鄰核小體間的間隔存在差異,然而核小體的定位在不同細(xì)胞中趨于一致;在TH1增強子處,不同CD4+幼稚T細(xì)胞顯示不同的核小體缺失程度,在TH1增強子處缺失程度高的細(xì)胞具有更高TH1分化的潛能。

    2 MNase-seq技術(shù)解析核小體組織特點及內(nèi)在功能

    2.1 核小體組織特點及影響因素

    核小體覆蓋了可測序基因組的大部分區(qū)域,整體而言,核小體重復(fù)地、有規(guī)律地出現(xiàn)在基因組上。釀酒酵母基因組中,核小體以約165 bp的間隔重復(fù)出現(xiàn)[41],相鄰核小體由10~50 bp的連接DNA(linker DNA)相連[42~45]。然而,核小體在基因組上的組織具有異質(zhì)性,導(dǎo)致不同的染色質(zhì)折疊方式。核小體在異染色質(zhì)區(qū)域覆蓋率高,形成了緊湊封閉的染色質(zhì)結(jié)構(gòu);然而在增強子、絕緣子等調(diào)控區(qū)域常表現(xiàn)為核小體缺失,形成開放可及的染色質(zhì)[41]。并且,核小體在基因組上的組織不是固定不變的,具有定位和結(jié)構(gòu)上的動態(tài)性[46]。核小體定位的動態(tài)性體現(xiàn)在核小體可以在DNA上滑動,且容易自發(fā)發(fā)生全部或部分解聚的動態(tài)變化。在生理條件下,核小體易受染色質(zhì)重塑復(fù)合物(chromatin remodelers)、分子伴侶、聚合酶、轉(zhuǎn)錄因子等調(diào)控因子的影響而發(fā)生重定位;全局轉(zhuǎn)錄水平的變化以及外界環(huán)境刺激也會影響核小體的定位。核小體結(jié)構(gòu)的動態(tài)性體現(xiàn)在組蛋白的翻譯后修飾(post-translational modifications, PTM)[47,48],以及組蛋白變體[49,50]和非標(biāo)準(zhǔn)核小體的形成[51]等。

    研究人員使用核小體定位(nucleosome position-ning)及占位(nucleosome occupancy)來描述核小體的組織狀態(tài)。核小體定位指核小體出現(xiàn)在基因組特定位置相對于其周邊的概率,反映核小體對特定DNA序列選擇的特性;核小體占位是指在基因組特定區(qū)域出現(xiàn)的核小體平均數(shù)目,體現(xiàn)了核小體密度。核小體的定位和占位具有堿基偏好性,通常,富含G、C堿基的DNA序列更有利于DNA與核小體的緊密結(jié)合,而高A、T含量DNA序列削弱了DNA與組蛋白的相互作用,連續(xù)出現(xiàn)的A堿基區(qū)域(AAAAA)在體內(nèi)和體外實驗中都表現(xiàn)出最低的核小體占有率[52]。Albert等[53]利用MNase-ChIP-seq研究釀酒酵母核小體組織時發(fā)現(xiàn)并定義了DNA序列影響核小體定位的旋轉(zhuǎn)特性(rotational setting)和平移特性(translational setting)。其中旋轉(zhuǎn)特性與規(guī)律排列的二核苷酸有關(guān),以10 bp為周期連續(xù)出現(xiàn)的AA/AT/ TA/TT二核苷酸以及相位相差5 bp的GG/CC/GC/ CG二核苷酸序列交錯出現(xiàn),塑造了DNA的急劇彎曲特性,有利于DNA對組蛋白八聚體的纏繞,從而具有強烈的核小體定位特性。

    目前普遍認(rèn)為DNA序列對于核小體在基因組上的組織起到主要作用,僅通過核小體對DNA的序列偏好可以解釋體內(nèi)核小體組織形態(tài)的50%~60%[3]。Kaplan等[52]利用MNase-seq對培養(yǎng)于3種不同培養(yǎng)基的釀酒酵母體內(nèi)的核小體進行測定,它們的核小體圖譜顯示出較高的相似性,表明盡管存在環(huán)境差異,DNA序列仍是影響核小體組織的最關(guān)鍵因素。根據(jù)核小體定位的序列偏好性,研究人員根據(jù)DNA序列及體外合成的核小體數(shù)據(jù)建立計算模型,對生理條件下核小體組織形式進行預(yù)測[3,54]。在基因組上不同區(qū)域,DNA序列與核小體占位的相關(guān)性不同。在啟動子區(qū)域,DNA序列與體內(nèi)核小體占位的相關(guān)系數(shù)較非啟動子區(qū)域更低。由此說明,除了DNA序列的作用,核小體在生物體內(nèi)的組織形態(tài)受到眾多細(xì)胞內(nèi)調(diào)控因子及轉(zhuǎn)錄水平的影響。

    以染色質(zhì)重塑復(fù)合物為主的眾多染色質(zhì)調(diào)節(jié)因子也影響著核小體在基因組上的定位。染色質(zhì)重塑復(fù)合物是一種依賴ATP的酶類,它們可以越過核小體對DNA序列的內(nèi)在偏好性,利用ATP水解的能量移除、移動或并入組蛋白來改變核小體的定位及構(gòu)象。在促進轉(zhuǎn)錄因子與DNA的結(jié)合[55]、順式作用元件構(gòu)象建成及DNA復(fù)制、轉(zhuǎn)錄激活等過程中發(fā)揮重要作用[56]。SWI/SNF (switch/sucrose non- fermenting)復(fù)合物是研究得較為全面的染色質(zhì)重塑復(fù)合物,它們通常富集在轉(zhuǎn)錄起始位點、復(fù)制起點中核小體缺失位點(nucleosome-free region, NFR)的–1核小體處[57],主要通過促進核小體在DNA上的移動或移除來發(fā)揮功能,參與釀酒酵母的應(yīng)激反應(yīng)[58]。近年來,SWI/SNF的作用機制被研究得更為透徹,研究表明,在黑色素細(xì)胞的分化中,SWI/SNF亞基BAF60A促進色素基因及啟動子對SWI/ SNF亞基BRG1的募集,從而促進染色質(zhì)重塑與細(xì)胞分化[59]。ISWI(imitation SWI)復(fù)合物屬于SWI蛋白家族,它能協(xié)助染色質(zhì)的組裝和組織[60,61],幫助具有規(guī)律間隔的核小體序列的建成[62]。CHD (chro-modomain-helicase-DNA-binding)也屬于SWI蛋白家族,它與ISWI共同決定了基因組上核小體間距(nucleosome spacing)的全局特征[63]。

    2.2 核小體缺失位點的結(jié)構(gòu)和功能

    在利用下一代測序技術(shù)研究核小體組織之前,研究者們利用MNase切割染色質(zhì)纖維并配合基因芯片技術(shù)(MNase-chip),在釀酒酵母的核小體組織上有了開創(chuàng)性的發(fā)現(xiàn),即核小體在啟動子區(qū)域普遍顯示出缺失特性[64~67],這些區(qū)域被命名為核小體缺失位點。NFR也被發(fā)現(xiàn)存在于活躍增強子區(qū)域、復(fù)制起點以及轉(zhuǎn)錄因子結(jié)合位點[67]。核小體的缺失增加了該位點的染色質(zhì)可及性、具有更強的調(diào)控潛力,利于眾多反式作用因子正確發(fā)揮功能,包括染色質(zhì)調(diào)控因子、轉(zhuǎn)錄因子、復(fù)制和轉(zhuǎn)錄需要的酶類等[68]。

    根據(jù)染色質(zhì)的可及性可將啟動子分為兩種類型:開放啟動子(open promoter)和封閉啟動子(covered promoter)[62]。開放啟動子具有開放的染色質(zhì)狀態(tài),在起始密碼子上游200 bp左右具有一段NFR,NFR內(nèi)具有暴露的轉(zhuǎn)錄激活因子結(jié)合位點;而封閉啟動子具有較高的核小體占位,轉(zhuǎn)錄激活因子需要在染色質(zhì)重塑復(fù)合物的幫助下與核小體競爭結(jié)合位點才能開始轉(zhuǎn)錄。NFR區(qū)域上下游出現(xiàn)的第一個核小體被命名為–1/+1核小體,分別標(biāo)定了NFR的上游/下游邊界[46]。+1核小體定位在轉(zhuǎn)錄起始位點(trans-cription start site, TSS)下游的固定距離,它所結(jié)合的DNA序列具有很強的核小體定位特性。利用MNase- seq對不同物種啟動子的核小體結(jié)構(gòu)進行測量時發(fā)現(xiàn),+1核小體的位置具有物種特異性。在酵母中,+1核小體與TSS在位置上存在重疊[69],而在果蠅的基因啟動子處,+1核小體通常位于TSS下游約135 bp處,因此具有比酵母更長的NFR[42](圖2)。在表達(dá)基因TSS的+1核小體下游方向,形成了規(guī)律間隔、定位精準(zhǔn)且相位統(tǒng)一的核小體陣列(nucleosome phasing),這一陣列大約向基因內(nèi)部延伸1000 bp左右,規(guī)律排列的特性隨著與TSS間距離的增大而減弱[46]。Dustin等[44]在人類CD4+T淋巴細(xì)胞中發(fā)現(xiàn),+1核小體下游相位統(tǒng)一的核小體陣列在不表達(dá)的基因中不存在。并且,將啟動子區(qū)域核小體排列與RNA聚合酶Ⅱ的ChIP-seq數(shù)據(jù)比對發(fā)現(xiàn),啟動子區(qū)域RNA聚合酶II水平越高,+1核小體及其下游的核小體定相現(xiàn)象越明顯;不同種類的RNA聚合酶II影響著+1核小體的位置。

    圖2 酵母和果蠅+1/?1核小體的位置

    A:釀酒酵母+1/?1核小體位置;B:果蠅+1/?1核小體位置。

    啟動子區(qū)域NFR的功能主要體現(xiàn)在兩個方面:首先,轉(zhuǎn)錄起始位點的NFR有助于轉(zhuǎn)錄起始前復(fù)合物的組裝從而開始轉(zhuǎn)錄,而轉(zhuǎn)錄終止位點附近的NFR有助于轉(zhuǎn)錄復(fù)合物的解聚;另外,NFR下游+1核小體與DNA的解聚或結(jié)構(gòu)變化促進了RNA聚合酶Ⅱ與DNA的結(jié)合,促進了轉(zhuǎn)錄的進行。+1核小體的位置在轉(zhuǎn)錄活動中發(fā)揮重要作用,其位置向NFR上游偏移將會影響轉(zhuǎn)錄元件的組裝,降低轉(zhuǎn)錄效率[70,71]。在啟動子區(qū)域以外,也發(fā)現(xiàn)有NFR及相位統(tǒng)一的核小體陣列的存在,例如絕緣子復(fù)合物及轉(zhuǎn)錄因子的結(jié)合位點。絕緣子CTCF(CCCTC-binding factor) 蛋白結(jié)合位點處核小體占位很低,表現(xiàn)出核小體缺失,并且在CTCF結(jié)合位點上下游發(fā)現(xiàn)具有對稱且規(guī)律的相位統(tǒng)一的核小體陣列[72]。與+1核小體類似,這些蛋白質(zhì)與DNA結(jié)合很緊密,具有強定位特性,但基因組上核小體定相排布的機制還未完全破解[41]。

    2.3 核小體組織與基因表達(dá)調(diào)控

    核小體對DNA的包裝,一方面壓縮了DNA,阻礙了眾多DNA結(jié)合蛋白與DNA的相互作用;另一方面,在這一過程中核小體能調(diào)整染色質(zhì)的包裝方式,促進了細(xì)胞內(nèi)基因的正確表達(dá)[62]。核小體的定位與轉(zhuǎn)錄因子的結(jié)合以及基因的轉(zhuǎn)錄水平密切相關(guān),且相互影響。

    轉(zhuǎn)錄因子與核小體競爭DNA結(jié)合位點會導(dǎo)致核小體組織的變化,進而重塑該位點的染色質(zhì)可及性狀態(tài)從而調(diào)控基因表達(dá)。對于大多數(shù)轉(zhuǎn)錄因子而言,其對應(yīng)的結(jié)合位點上有核小體存在時,由于空間位阻、電荷相斥等影響,轉(zhuǎn)錄因子與DNA的親和性通常比裸露DNA低10倍以上。然而結(jié)構(gòu)蛋白(architectural proteins, AP)和染色質(zhì)重塑復(fù)合物能幫助轉(zhuǎn)錄因子與DNA的結(jié)合,這一過程通常促進核小體從DNA上分離。Daniel等[10]利用ChIP-seq結(jié)合MNase-seq等技術(shù)對35種轉(zhuǎn)錄因子結(jié)合位點附近的核小體組織形態(tài)進行分析,他們發(fā)現(xiàn)在轉(zhuǎn)錄因子結(jié)合位點峰值處通常形成NFR,且轉(zhuǎn)錄因子占位水平與其結(jié)合位點上下游核小體定位強度負(fù)相關(guān)。轉(zhuǎn)錄因子結(jié)合使得核小體解聚,進一步促進開放染色質(zhì)的形成從而激活轉(zhuǎn)錄。釀酒酵母基因啟動子是研究染色質(zhì)結(jié)構(gòu)對基因表達(dá)影響的重要模型。研究發(fā)現(xiàn),在轉(zhuǎn)錄抑制狀態(tài)的啟動子處,核小體呈規(guī)律排列形式,而轉(zhuǎn)錄因子Pho4p與UASp2位點的結(jié)合能導(dǎo)致核小體組織重構(gòu),表現(xiàn)為–2和–3核小體解聚,核小體占位水平很低。這一核小體重構(gòu)導(dǎo)致了染色質(zhì)開放位點的形成,進而激活了PHO5的表達(dá)[73,74]。轉(zhuǎn)錄因子還可以通過改變局部染色質(zhì)的空間結(jié)構(gòu)來影響轉(zhuǎn)錄。基因具有兩個重要的調(diào)控元件boxA與boxB,轉(zhuǎn)錄因子TFⅢC與這兩個調(diào)控元件的結(jié)合使boxA與boxB間的核小體向上游平移約40 bp。而后,轉(zhuǎn)錄因子TFⅢB的結(jié)合導(dǎo)致–1核小體重定位從而改變TATA上游染色質(zhì)空間結(jié)構(gòu),即在TATA框處形成一段對核酸酶高敏的活躍染色質(zhì),從而為轉(zhuǎn)錄起始前復(fù)合物(pre-initiation complex, PIC)的組裝和轉(zhuǎn)錄做好了準(zhǔn)備[75]。

    轉(zhuǎn)錄水平同樣也會對核小體在基因組上的定位和占位水平產(chǎn)生影響。Sushma等[76]利用MNase-seq對熱激(heat shock)前后釀酒酵母全基因組核小體定位及占位進行測量,發(fā)現(xiàn)細(xì)胞在經(jīng)歷轉(zhuǎn)錄干擾(transcriptional perturbation)前后,不發(fā)生全局范圍的核小體定位變化,大部分核小體的位置保持穩(wěn)定,染色質(zhì)重塑活動通常只與基因啟動子區(qū)域單個或兩個核小體的缺失或變換相關(guān)。在人CD4+T細(xì)胞中,Anton等[72]發(fā)現(xiàn)在不同轉(zhuǎn)錄頻率的基因啟動子處,其NFR的缺失程度不同,高表達(dá)基因(>8 RPKM)較低表達(dá)基因(<1 RPKM)的核小體缺失程度更高,因此轉(zhuǎn)錄活性可能與NFR區(qū)域核小體移除相關(guān)。轉(zhuǎn)錄活動還會影響核小體間隔,Lai等[35]利用單細(xì)胞MNase-seq技術(shù),對小鼠CD4+T淋巴細(xì)胞和胚胎干細(xì)胞的核小體組織及染色質(zhì)可及性進行研究,發(fā)現(xiàn)活躍轉(zhuǎn)錄基因內(nèi)部核小體間的間隔不均勻;而在沉默基因或異染色質(zhì)區(qū)域,由于沒有轉(zhuǎn)錄活動的影響,核小體間隔高度均勻,連接DNA長度變化不大。

    核小體在基因組上的排列多起到抑制轉(zhuǎn)錄的作用,通常與影響RNA聚合酶II的延伸及PIC的組裝有關(guān)。轉(zhuǎn)錄過程中,核小體的存在成為了RNA聚合酶Ⅱ的延伸的物理屏障,在體內(nèi)和體外的實驗中都發(fā)現(xiàn)了核小體參與RNA聚合酶II延伸暫停(pausing process),從而影響轉(zhuǎn)錄效率[42,77];基因啟動子區(qū)域的核小體組織特別是+1核小體的定位直接決定了PIC是否可以成功組裝。Reja等[78]利用ChIP-exo技術(shù)研究酵母核糖體蛋白啟動子時發(fā)現(xiàn),經(jīng)熱激處理后,+1核小體向上游啟動子區(qū)域移動導(dǎo)致了轉(zhuǎn)錄抑制,因為+1核小體對上游DNA序列的占領(lǐng)阻礙了PIC與足夠DNA底物的結(jié)合。然而,核小體對轉(zhuǎn)錄的抑制作用可以被細(xì)胞內(nèi)多種調(diào)控因子逆轉(zhuǎn),在人IFN-β基因啟動子處,+1核小體對TATA框(TATA box)及TSS的占位阻礙了PIC的組裝,抑制基因轉(zhuǎn)錄。當(dāng)機體受到病毒感染,在增強體(enhanceosome)及其招募的SWI/SNF、GCN5乙酰轉(zhuǎn)移酶等多種復(fù)合物的共同作用下,+1核小體向下游移動約36 bp,使得轉(zhuǎn)錄因子TFⅡD與核心啟動子區(qū)域結(jié)合,從而起始轉(zhuǎn)錄過程[79~81]。

    2.4 組蛋白修飾圖譜與表觀遺傳學(xué)

    核小體不僅是染色質(zhì)組裝的基本單位,也是生物體表觀遺傳學(xué)修飾的主要載體。核小體的定位及組蛋白翻譯后修飾是表觀遺傳學(xué)在各領(lǐng)域研究中的重要內(nèi)容。其中,MNase-ChIP-seq技術(shù)發(fā)揮了重要功能,該技術(shù)目前已經(jīng)可以對全基因組范圍內(nèi)的組蛋白PTM進行捕獲并測量。

    組蛋白PTM被研究者們稱作“組蛋白密碼”(histone code),因為PTM可以直接改變?nèi)旧|(zhì)的物理結(jié)構(gòu),或被細(xì)胞內(nèi)特定的蛋白質(zhì)識別,進而激活或抑制下游染色質(zhì)功能[82]。組蛋白翻譯后修飾中較普遍的有賴氨酸的乙?;?、甲基化及泛素化,精氨酸的甲基化以及絲氨酸的磷酸化[83,84],這些修飾多出現(xiàn)在轉(zhuǎn)錄調(diào)控元件如啟動子、增強子處。作為一種常見的修飾形式,組蛋白乙?;魅趿私M蛋白-DNA的相互作用、提高核小體周轉(zhuǎn)率和染色質(zhì)可及性,利于轉(zhuǎn)錄因子的結(jié)合從而促進基因表達(dá),在高轉(zhuǎn)錄活性的基因啟動子區(qū)域豐度較高[85~87]。研究表明,啟動子區(qū)域H3K9ac修飾水平與基因表達(dá)水平正相關(guān),且H3K9ac被證實在促進人胚胎干細(xì)胞分化[88]、胰島β細(xì)胞增殖[89]等過程中起關(guān)鍵作用。相比于乙?;揎?,組蛋白甲基化修飾形式更加穩(wěn)定。一些甲基化形式能促進轉(zhuǎn)錄,如H3K4me1/3、H3K36me3等,其中H3K4me3主要分布于活躍轉(zhuǎn)錄基因啟動子處,它有助于啟動子對轉(zhuǎn)錄因子TFⅡD、RNA聚合酶Ⅱ的招募;在增強子處,多出現(xiàn)明顯的H3K4me1修飾,而不出現(xiàn)H3K4me3修飾[90]。另一些組蛋白甲基化與基因沉默相關(guān),H3K9me3、H3K27me3在哺乳動物細(xì)胞中促進染色質(zhì)凝集[91]。

    早期的研究方向多集中于H2A.Z核小體變體定位[53]與組蛋白甲基化分布形式[92]。H2A.Z作為常見的核小體變體,H2A.Z參與轉(zhuǎn)錄調(diào)控、DNA修復(fù)等眾多生物過程[93]。H2A.Z通常出現(xiàn)在基因啟動子處+1/?1核小體處,專一定位在+1核小體處的H2A.Z與轉(zhuǎn)錄起始相關(guān)[94]。H2A.Z的占位水平由它的并入水平與解離水平共同決定。其中,SWR1促進了H2A.Z的并入,而RNA聚合酶II以及Kin28/Cdk7 激酶參與H2A.Z的解離,目前,已經(jīng)可以在活細(xì)胞中以單分子的分辨率追蹤核小體的變化過程[95]。另一方面,H2A.Z占位水平通常與DNA甲基化水平負(fù)相關(guān),Murphy等[96]發(fā)現(xiàn)在斑馬魚()的配子和胚胎發(fā)育時期,H2A.Z促使DNA甲基化模式“重編程”。其中,定位在編碼胚胎早期轉(zhuǎn)錄因子基因啟動子區(qū)的H2A.Z使得啟動子處于低甲基化水平從而促進基因表達(dá),這一過程在斑馬魚胚胎發(fā)育中發(fā)揮重要作用。

    組蛋白修飾狀態(tài)的改變通常與基因功能改變、細(xì)胞惡性轉(zhuǎn)化相關(guān)。目前,研究人員可以利用下一代測序技術(shù)探究染色質(zhì)表觀遺傳改變,從而分析細(xì)胞健康狀態(tài)。癌細(xì)胞通常具有全局范圍的表觀遺傳異常,其中,特定的組蛋白修飾與癌癥的起始和增殖緊密相關(guān)[97]。Wang等[98]發(fā)現(xiàn)了H3乙?;c腫瘤細(xì)胞耐藥性的關(guān)系。他們通過抑制組蛋白去乙?;窰DAC3的表達(dá),促進了癌細(xì)胞基因啟動子區(qū)域的H3乙?;?,從而增強基因的表達(dá)。而PD-L1有抑制T細(xì)胞活性的作用,使得癌細(xì)胞具有更強的耐藥性。這一表觀遺傳學(xué)機制的發(fā)現(xiàn)有利于發(fā)現(xiàn)治療腫瘤抗藥能力的潛在靶點。Aman等[99]發(fā)現(xiàn)在腎細(xì)胞癌中,低水平的H3K4甲基化通常與晚期癌癥和腫瘤轉(zhuǎn)移相關(guān),其中H3K4me3 可能是一個較好的腫瘤轉(zhuǎn)移預(yù)測標(biāo)志物。Liana等[100]發(fā)現(xiàn)H3K9ac促進了口腔粘膜癌細(xì)胞增殖的上皮間質(zhì)轉(zhuǎn)化(epithelial–mesenchymal transition, EMT),因此H3K9ac可以作為口腔粘膜癌預(yù)后評估的關(guān)鍵標(biāo)志物。表觀遺傳通常是動態(tài)可逆的,因此根據(jù)組蛋白表觀遺傳相關(guān)酶類設(shè)計藥物,可以恢復(fù)“正常的表觀狀態(tài)”從而起到治療效果,以組蛋白甲基轉(zhuǎn)移酶抑制劑、組蛋白去乙酰化酶抑制劑為主的表觀遺傳學(xué)藥物發(fā)展前景較為樂觀[97]。

    3 結(jié)語與展望

    21世紀(jì)以來,下一代測序技術(shù)的高速發(fā)展加快了研究人員對染色質(zhì)結(jié)構(gòu)的研究,使人們在時間和空間上對核小體組織有了更清晰的認(rèn)知。MNase-seq及其衍生技術(shù)揭示了全基因組范圍核小體的組織方式,即核小體在基因組上的分布是動態(tài)且不均勻的,在特定區(qū)域會出現(xiàn)NFR。人們已經(jīng)認(rèn)識到核小體定位主要受到DNA序列的影響,染色質(zhì)重塑因子、轉(zhuǎn)錄因子、聚合酶等同樣參與了核小體組織及染色質(zhì)結(jié)構(gòu)的形成。染色質(zhì)構(gòu)象的建立及動態(tài)平衡使得調(diào)控元件高效地對轉(zhuǎn)錄因子進行招募、調(diào)控遠(yuǎn)端序列與靶基因的互作,因此與染色質(zhì)功能特別是基因表達(dá)調(diào)控密切相關(guān)。

    近年來,對MNase-seq技術(shù)改良的重點在于減少必要的細(xì)胞投入量及降低背景信號,目前該技術(shù)已經(jīng)可以在單個細(xì)胞上實現(xiàn)[35],有助于人們從單個細(xì)胞的核小體組織結(jié)構(gòu)層面理解細(xì)胞間的染色質(zhì)異質(zhì)性。在疾病研究中,將scMNase-seq與scDNase-seq、scChIP-seq以及多種表觀遺傳學(xué)分析方法相結(jié)合,有助于人們理解腫瘤異質(zhì)性、剖析腫瘤形成過程內(nèi)在的染色質(zhì)結(jié)構(gòu)基礎(chǔ)、分析不同疾病的表觀遺傳標(biāo)記,從而解決更多疾病機理問題[101]。利用MNase-seq及相關(guān)技術(shù)對全基因組核小體定位和可及性的測量初步解析了染色質(zhì)結(jié)構(gòu)和基因調(diào)控機制,最近發(fā)現(xiàn)MNase-seq還可以預(yù)測高階染色質(zhì)結(jié)構(gòu),然而其在染色質(zhì)三維相互作用的分析中存在困難。將MNase- seq技術(shù)與近年來發(fā)展迅速的Hi-C[102,103]等染色質(zhì)捕獲技術(shù)相結(jié)合,有助于人們更好地理解染色質(zhì)結(jié)構(gòu)與其功能的關(guān)系。隨著技術(shù)方法的不斷進步以及多學(xué)科的交叉融合,人們將對復(fù)雜動態(tài)的染色質(zhì)結(jié)構(gòu)及功能研究得更加深入。

    [1] Luger K, M?der AW, Richmond RK, Sargent DF, Richmond TJ. Crystal structure of the nucleosome core particle at 2.8 ? resolution.,1997, 389(6648): 251–260.

    [2] Zhou KD, Gaullier G, Luger K. Nucleosome structure and dynamics are coming of age.,2019, 26(1): 3–13.

    [3] Segal E, Fondufe-Mittendorf Y, Chen L, Th?str?m A, Field Y, Moore IK, Wang JZ, Widom J. A genomic code for nucleosome positioning.,2006, 442(7104): 772–778.

    [4] Noll M. Subunit structure of chromatin., 1974, 251(5472): 249–251.

    [5] Lohr D, Kovacic RT, Van Holde KE. Quantitative analysis of the digestion of yeast chromatin by staphylococcal nuclease., 1977, 16(3): 463–471.

    [6] Fan XC, Moqtaderi Z, Jin Y, Zhang Y, Liu XS, Struhl K. Nucleosome depletion at yeast terminators is not intrinsic and can occur by a transcriptional mechanism linked to 3’-end formation.,2010, 107(42): 17945–17950.

    [7] Tsompana M, Buck MJ. Chromatin accessibility: a window into the genome., 2014, 7(1): 33.

    [8] Zhang WL, Jiang JM. Application of MNase-Seq in the global mapping of nucleosome positioning in plants., 2018, 1830: 353–366.

    [9] Pajoro A, Mui?o JM, Angenent GC, Kaufmann K. Profiling nucleosome occupancy by MNase-seq: experi-mental protocol and computational analysis., 2018, 1675: 167–181.

    [10] Gaffney DJ, McVicker G, Pai AA, Fondufe-Mittendorf YN, Lewellen N, Michelini K, Widom J, Gilad Y, Pitchard JKControls of nucleosome positioning in the human genome., 2012, 8(11): e1003036.

    [11] Zhang ZH, Pugh BF. High-resolution genome-wide mapping of the primary structure of chromatin., 2011, 144(2): 175–186.

    [12] Klein DC, Hainer SJ. Genomic methods in profiling DNA accessibility and factor localization., 2020, 28(1): 69–85.

    [13] Chen WZ, Liu Y, Zhu SS, Green CD, Wei G, Han JDJ. Improved nucleosome-positioning algorithm iNPS for accurate nucleosome positioning from sequencing data., 2014, 5: 4909.

    [14] Fu K, Tang QZ, Feng JX, Liu XS, Zhang Y. DiNuP: a systematic approach to identify regions of differential nucleosome positioning., 2012, 28(15): 1965–1971.

    [15] Yuan GC, Liu JS. Genomic sequence is highly predictive of local nucleosome depletion., 2008, 4(1): e13.

    [16] Wal M, Pugh BF. Genome-wide mapping of nucleosomepositions in yeast using high-resolution MNase ChIP-Seq., 2012, 513: 233–250.

    [17] Ocampo J, Cui F, Zhurkin VB, Clark DJ. The proto-chromatosome: A fundamental subunit of chro-matin?, 2016, 7(4): 382–387.

    [18] Mieczkowski J, Cook A, Bowman SK, Mueller B, Alver BH, Kundu S, Deaton AM, Urban JA, Larschan E, Park PJ, Kingston RE, Tolstorukov MYMNase titration reveals differences between nucleosome occupancy and chromatin accessibility., 2016, 7: 11485.

    [19] Zhao HN, Zhang WL, Zhang T, Lin Y, Hu YD, Fang C, Jiang JMGenome-wide MNase hypersensitivity assay unveils distinct classes of open chromatin associated with H3K27me3 and DNA methylation in., 2020, 21(1): 24.

    [20] Baldi S, Krebs S, Blum H, Becker PB. Genome-wide measurement of local nucleosome array regularity and spacing by nanopore sequencing., 2018, 25(9): 894–901.

    [21] Skene PJ, Henikoff S. An efficient targeted nuclease strategy for high-resolution mapping of DNA binding sites., 2017, 6: e21856.

    [22] Rizzo JM, Bard JE, Buck MJ. Standardized collection of MNase-seq experiments enables unbiased dataset comparisons., 2012, 13: 15.

    [23] Rizzo JM, Sinha S. Analyzing the global chromatin structure of keratinocytes by MNase-seq., 2014, 1195: 49–59.

    [24] Cui KR, Zhao KJ. Genome-wide approaches to determining nucleosome occupancy in metazoans using MNase-Seq., 2012, 833: 413–419.

    [25] Meyer CA, Liu XS. Identifying and mitigating bias in next-generation sequencing methods for chromatin biology., 2014, 15(11): 709–721.

    [26] Park PJ. ChIP-seq: advantages and challenges of a maturing technology., 2009, 10(10): 669–680.

    [27] Mardis ER. ChIP-seq: welcome to the new frontier., 2007, 4(8): 613–614.

    [28] Buenrostro JD, Wu BJ, Chang HY, Greenleaf WJ. ATAC-seq: A method for assaying chromatin accessibility genome-Wide., 2015, 109: 21.29.1–21.29.9.

    [29] Buenrostro JD, Giresi PG, Zaba LC, Chang HY, Greenleaf WJ. Transposition of native chromatin for fast and sensitive epigenomic profiling of open chromatin, DNA-binding proteins and nucleosome position., 2013, 10(12): 1213–1218.

    [30] Schep AN, Buenrostro JD, Denny SK, Schwartz K, Sherlock G, Greenleaf WJ. Structured nucleosome fingerprints enable high-resolution mapping of chromatin architecture within regulatory regions., 2015, 25(11): 1757–1770.

    [31] Thurman RE, Rynes E, Humbert R, Vierstra J, Maurano MT, Haugen E,Sheffield NC, Stergachis AB, Wang H, Vernot B, Garg K, John S, Sandstrom R, Bates D, Boatman L, Canfield TK, Diegel M, Dunn D, Ebersol AK, Frum T, Giste E, Johnson AK, Johnson EM, Kutyavin T, Lajoie B, Lee BK, Lee K, London D, Lotakis D, Neph S, Neri F, Nguyen ED, Qu H, Reynolds AP, Roach V, Safi A, Sanchez ME, Sanyal A, Shafer A, Simon JM, Song LY, Vong S, Weaver M, Yan YQ, Zhang ZC, Zhang ZZ, Lenhard B, Tewari M, Dorschner MO, Hansen RS, Navas PA, Stamatoyannopoulos G, Iyer VR, Lieb JD, Sunyaev SR, Akey JM, Sabo PJ, Kaul R, Furey TS, Dekker J, Crawford GE, Stamatoyan-nopoulos JA. The accessible chromatin landscape of the human genome., 2012, 489(7414): 75–82.

    [32] Zhong JL, Luo KX, Winter PS, Crawford GE, Iversen ES, Hartemink AJ. Mapping nucleosome positions using DNase-seq., 2016, 26(3): 351–364.

    [33] Kelly TK, Liu YP, Lay FD, Liang GN, Berman BP, Jones PA. Genome-wide mapping of nucleosome positioning and DNA methylation within individual DNA molecules., 2012, 22(12): 2497– 2506.

    [34] Krebs AR, Imanci D, Hoerner L, Gaidatzis D, Burger L, Schübeler D. Genome-wide single-molecule footprinting reveals high RNA polymerase II turnover at paused promoters., 2017, 67(3): 411–422.e4.

    [35] Lai BB, Gao WW, Cui KR, Xie WL, Tang QS, Jin WF, Hu GQ, Ni B, Zhao KJPrinciples of nucleosome organization revealed by single-cell micrococcal nuclease sequencing., 2018, 562(7726): 281–285.

    [36] Lion M, Tolstorukov MY, Oettinger MA. Low-Input MNase accessibility of chromatin (low-input MACC)., 2019, 127(1): e91.

    [37] Skene PJ, Henikoff JG, Henikoff S. Targeted in situ genome-wide profiling with high efficiency for low cell numbers., 2018, 13(5): 1006–1019.

    [38] Hainer SJ, Bo?kovi? A, McCannell KN, Rando OJ, Fazzio TG. Profiling of pluripotency factors in single cells and early embryos., 2019, 177(5): 1319– 1329.e11.

    [39] Gao WW, Lai BB, Ni B, Zhao KJ. Genome-wide profiling of nucleosome position and chromatin accessibility in single cells using scMNase-seq., 2020, 15(1): 68–85.

    [40] Baldi S. Nucleosome positioning and spacing: from genome-wide maps to single arrays., 2019, 63(1): 5–14.

    [41] Baldi S, Korber P, Becker PB. Beads on a string- nucleosome array arrangements and folding of the chromatin fiber., 2020, 27(2): 109–118.

    [42] Mavrich TN, Jiang CZ, Ioshikhes IP, Li XY, Venters BJ, Zanton SJ, Tomsho LP, Qi J, Glaser RL, Schuster SC, Gilmour DS, Albert I, Pugh BFNucleosome organization in thegenome., 2008, 453(7193): 358–362.

    [43] Lee W, Tillo D, Bray N, Morse RH, Davis RW, Hughes TR, Nislow CA high-resolution atlas of nucleosome occupancy in yeast., 2007, 39(10): 1235– 1244.

    [44] Schones DE, Cui KR, Cuddapah S, Roh TY, Barski A, Wang ZB, Wei G, Zhao KJDynamic regulation of nucleosome positioning in the human genome., 2008, 132(5): 887–898.

    [45] Valouev A, Ichikawa J, Tonthat T, Stuart J, Ranade S, Peckham H, Zeng K, Malek JA, Costa G, McKernan K, Sidow A, Fire A, Johnson SMA high-resolution, nucleosome position map ofreveals a lack of universal sequence-dictated positioning., 2008, 18(7): 1051–1063.

    [46] Lai WKM, Pugh BF. Understanding nucleosome dynamics and their links to gene expression and DNA replication., 2017, 18(9): 548– 562.

    [47] Martin C, Zhang Y. Mechanisms of epigenetic inheritance., 2007, 19(3): 266–272.

    [48] Ruthenburg AJ, Li HT, Patel DJ, Allis CD. Multivalent engagement of chromatin modifications by linked binding modules., 2007, 8(12): 983–994.

    [49] Ahmad K, Henikoff S. The histone variant H3.3 marks active chromatin by replication-independent nucleosome assembly., 2002, 9(6): 1191–1200.

    [50] Sarma K, Reinberg D. Histone variants meet their match., 2005, 6(2): 139–149.

    [51] Ramachandran S, Zentner GE, Henikoff S. Asymmetric nucleosomes flank promoters in the budding yeast genome., 2015, 25(3): 381–390.

    [52] Kaplan N, Moore IK, Fondufe-Mittendorf Y, Gossett AJ, Tillo D, Field Y, LeProust EM, Hughes TR, Lieb JD, Widom J, Segal EThe DNA-encoded nucleosome organization of a eukaryotic genome., 2009, 458(7236): 362–366.

    [53] Albert I, Mavrich TN, Tomsho LP, Qi J, Zanton SJ, Schuster SC, Pugh BFTranslational and rotational settings of H2A.Z nucleosomes across the Saccharomyces cerevisiae genome., 2007, 446(7135): 572–576.

    [54] Cui F, Chen LL, LoVerso PR, Zhurkin VB. Prediction of nucleosome rotational positioning in yeast and human genomes based on sequence-dependent DNA anisotropy., 2014, 15(1): 313.

    [55] de Dieuleveult M, Yen KY, Hmitou I, Depaux A, Boussouar F, Bou Dargham D, Jounier S, Humbertclaude H, Ribierre F, Baulard C, Farrell NP, Park B, Keime C, Carrière L, Berlivet S, Gut M, Gut L, Werner M, Deleuze JF, Olaso R, Aude JC, Chantalat S, Pugh BF, Gérard MGenome-wide nucleosome specificity and function of chromatin remodellers in ES cells., 2016, 530(7588): 113–116.

    [56] Ho L, Crabtree GR. Chromatin remodelling during development., 2010, 463(7280): 474–484.

    [57] Dutta A, Gogol M, Kim JH, Smolle M, Venkatesh S, Gilmore J, Florens L, Washburn MP, Workman JLSwi/Snf dynamics on stress-responsive genes is governed by competitive bromodomain interactions., 2014, 28(20): 2314–2330.

    [58] Ribeiro-Silva C, Vermeulen W, Lans H. SWI/SNF: Complex complexes in genome stability and cancer., 2019, 77: 87–95.

    [59] Aras S, Saladi SV, Basuroy T, Marathe HG, Lorès P, de la Serna IL. BAF60A mediates interactions between the microphthalmia-associated transcription factor and the BRG1-containing SWI/SNF complex during melanocyte differentiation., 2019, 234(7): 11780–11791.

    [60] Oppikofer M, Bai TY, Gan YT, Haley B, Liu P, Sandoval W, Ciferri C, Cochran AGExpansion of the ISWI chromatin remodeler family with new active complexes., 2017, 18(10): 1697–1706.

    [61] Levendosky RF, Bowman GD. Asymmetry between the two acidic patches dictates the direction of nucleosome sliding by the ISWI chromatin remodeler., 2019, 8: e45472.

    [62] Cairns BR. The logic of chromatin architecture and remodelling at promoters., 2009, 461(7261): 193–198.

    [63] Ocampo J, Chereji RV, Eriksson PR, Clark DJ. The ISW1 and CHD1 ATP-dependent chromatin remodelers compete to set nucleosome spacing., 2016, 44(10): 4625–4635.

    [64] Lee CK, Shibata Y, Rao B, Strahl BD, Lieb JD. Evidence for nucleosome depletion at active regulatory regions genome-wide., 2004, 36(8): 900–905.

    [65] Sekinger EA, Moqtaderi Z, Struhl K. Intrinsic histone-DNA interactions and low nucleosome density are important for preferential accessibility of promoter regions in yeast., 2005, 18(6): 735–748.

    [66] Bernstein BE, Liu CL, Humphrey EL, Perlstein EO, Schreiber SL. Global nucleosome occupancy in yeast., 2004, 5(9): R62.

    [67] Yuan GC, Liu YJ, Dion MF, Slack MD, Wu LF, Altschuler SJ, Rando OJGenome-scale identification of nucleosome positions in., 2005, 309(5734): 626–630.

    [68] Rando OJ, Ahmad K. Rules and regulation in the primary structure of chromatin., 2007, 19(3): 250–256.

    [69] Mavrich TN, Ioshikhes IP, Venters BJ, Jiang CZ, Tomsho LP, Qi J, Schuster SC, Albert I, Pugh BFA barrier nucleosome model for statistical positioning of nucleosomes throughout the yeast genome., 2008, 18(7): 1073–1083.

    [70] Kubik S, Bruzzone MJ, Challal D, Dreos R, Mattarocci S, Bucher P, Libri D, Shore DOpposing chromatin remodelers control transcription initiation frequency and start site selection., 2019, 26(8): 744–754.

    [71] Kubik S, O’Duibhir E, de Jonge WJ, Mattarocci S, Albert B, Falcone JL, Bruzzone MJ, Holstege FCP, Shore DSequence-directed action of RSC remodeler and general regulatory factors modulates +1 nucleosome position to facilitate transcription., 2018, 71(1): 89–102.e5.

    [72] Valouev A, Johnson SM, Boyd SD, Smith CL, Fire AZ, Sidow A. Determinants of nucleosome organization in primary human cells., 2011, 474(7352): 516– 520.

    [73] Boeger H, Griesenbeck J, Strattan JS, Kornberg RD. Removal of promoter nucleosomes by disassembly rather than sliding., 2004, 14(5): 667–673.

    [74] Ertel F, Dirac-Svejstrup AB, Hertel CB, Blaschke D, Svejstrup JQ, Korber P. In vitro reconstitution of PHO5 promoter chromatin remodeling points to a role for activator-nucleosome competition., 2010, 30(16): 4060–4076.

    [75] Shivaswamy S, Bhargava P. Positioned nucleosomes due to sequential remodeling of the yeast U6 small nuclear RNA chromatin are essential for its transcriptional activation., 2006, 281(15): 10461–10472.

    [76] Shivaswamy S, Bhinge A, Zhao YJ, Jones S, Hirst M, Iyer VR. Dynamic remodeling of individual nucleosomes across a eukaryotic genome in response to transcriptional perturbation., 2008, 6(3): e65.

    [77] Kulaeva OI, Hsieh FK, Chang HW, Luse DS, Studitsky VM. Mechanism of transcription through a nucleosome by RNA polymerase II., 2013, 1829(1): 76–83.

    [78] Reja R, Vinayachandran V, Ghosh S, Pugh BF. Molecular mechanisms of ribosomal protein gene coregulation., 2015, 29(18): 1942–1954.

    [79] Lomvardas S, Thanos D. Modifying gene expression programs by altering core promoter chromatin architecture., 2002, 110(2): 261–271.

    [80] Ford E, Thanos D. The transcriptional code of human IFN-beta gene expression., 2010, 1799(3–4): 328–336.

    [81] Au-Yeung N, Horvath CM. Transcriptional and chromatin regulation in interferon and innate antiviral gene expression., 2018, 44: 11–17.

    [82] Rothbart SB, Strahl BD. Interpreting the language of histone and DNA modifications., 2014, 1839(8): 627–643.

    [83] Felsenfeld G, Groudine M. Controlling the double helix., 2003, 421(6921): 448–453.

    [84] Qi HY, Zhang ZJ, Li YJ, Fang XD. Role of chromatin conformation in eukaryotic gene regulation., 2011, 33(12): 1291–1299.亓合媛, 張昭軍, 李雅娟, 方向東. 染色質(zhì)構(gòu)象調(diào)控真核基因的表達(dá). 遺傳, 2011, 33(12): 1291–1299.

    [85] Vogelauer M, Wu J, Suka N, Grunstein M. Global histone acetylation and deacetylation in yeast., 2000, 408(6811): 495–498.

    [86] Bernstein BE, Humphrey EL, Erlich RL, Schneider R, Bouman P, Liu JS, Kouzarides T, Schreiber SLMethylation of histone H3 Lys 4 in coding regions of active genes., 2002, 99(13): 8695–8700.

    [87] Klemm SL, Shipony Z, WJ. Chromatin accessibility and the regulatory epigenome., 2019, 20(4): 207–220.

    [88] Du YH, Liu ZP, Cao XK, Chen XL, Chen ZY, Zhang XB, Zhang XQ, Jiang CZNucleosome eviction along with H3K9ac deposition enhances Sox2 binding during human neuroectodermal commitment., 2017, 24(6): 1121–1131.

    [89] Daneshpajooh M, Bacos K, Bysani M, Bagge A, Ottosson Laakso E, Vikman P, Eliasson L, Mulder H, Ling CHDAC7 is overexpressed in human diabetic islets and impairs insulin secretion in rat islets and clonal beta cells., 2017, 60(1): 116–125.

    [90] Vallianatos CN, Raines B, Porter RS, Bonefas KM, Wu MC, Garay PM, Collette KM, Seo YA, Dou Y, Keegan, CE, Tronson NC, Iwase SMutually suppressive roles of KMT2A and KDM5C in behaviour, neuronal structure, and histone H3K4 methylation., 2020, 3(1): 278.

    [91] Sasidharan Nair V, El Salhat H, Taha RZ, John A, Ali BR, Elkord E. DNA methylation and repressive H3K9 and H3K27 trimethylation in the promoter regions of PD-1, CTLA-4, TIM-3, LAG-3, TIGIT, and PD-L1 genes in human primary breast cancer., 2018, 10: 78.

    [92] Barski A, Cuddapah S, Cui KR, Roh TY, Schones DE, Wang ZB, Wei G, Chepelev L, Zhao KJHigh-resolution profiling of histone methylations in the human genome., 2007, 129(4): 823–837.

    [93] Giaimo BD, Ferrante F, Herchenr?ther A, Hake SB, Borggrefe T. The histone variant H2A.Z in gene regulation., 2019, 12(1): 37.

    [94] Bagchi DN, Battenhouse AM, Park D, Iyer VR. The histone variant H2A.Z in yeast is almost exclusively incorporated into the +1 nucleosome in the direction of transcription., 2020, 48(1): 157–170.

    [95] Ranjan A, Nguyen VQ, Liu S, Wisniewski J, Kim JM, Tang XN, Mizuguchi G, Elalaoui E, Nickels TJ, Jou V, English BP, Zheng QS, Luk E, Lavis, LD, Lionnet T, Wu CLive-cell single particle imaging reveals the role of RNA polymerase II in histone H2A.Z eviction., 2020, 9: e55667.

    [96] Murphy PJ, Wu SF, James CR, Wike CL, Cairns BR. Placeholder nucleosomes underlie Germline-to-Embryo DNA methylation reprogramming., 2018, 172(5): 993–1006.e13.

    [97] Sharma S, Kelly TK, Jones PA. Epigenetics in cancer., 2010, 31(1): 27–36.

    [98] Wang HF, Fu C, Du J, Wang HS, He R, Yin XF, Li HX, Li X, Wang HX, Li K, Zheng L, Liu ZC, Qiu YREnhanced histone H3 acetylation of the PD-L1 promoter via the COP1/c-Jun/HDAC3 axis is required for PD-L1 expression in drug-resistant cancer cells., 2020, 39(1): 29.

    [99] Kumar A, Kumari N, Sharma U, Ram S, Singh SK, Kakkar N, Kaushal K, Prasad RReduction in H3K4me patterns due to aberrant expression of methyltransferases and demethylases in renal cell carcinoma: prognostic and therapeutic implications., 2019, 9(1): 8189.

    [100] Webber LP, Wagner VP, Curra M, Vargas PA, Meurer L, Carrard VC, Squarize CH, Castilho RM, Martins MDHypoacetylation of acetyl-histone H3 (H3K9ac) as marker of poor prognosis in oral cancer., 2017, 71(2): 278–286.

    [101] Torres CM, Biran A, Burney MJ, Patel H, Henser- Brownhill T, Cohen AHS, Li YL, Ben-Hamo R, Nye E, Spencer-Dene B, Chakravarty P, Efroni S, Matthews N, Misteli T, Meshorer E, Scaffidi PThe linker histone H1.0 generates epigenetic and functional intratumor heterogeneity., 2016, 353(6307): aaf1644.

    [102] Huang QT, Li Q, Zhang YB. Linking chromatin conformation to gene function., 2020, 42(1): 1–17.黃其通, 李清, 張玉波. 染色質(zhì)構(gòu)象與基因功能. 遺傳, 2020, 42(1): 1–17.

    [103] Denker A, de Laat W. The second decade of 3C technologies: detailed insights into nuclear organization., 2016, 30(12): 1357–1382.

    Resolving nucleosomal positioning and occupancy with MNase-seq

    Weihang Deng, Xinhui Li

    Nucleosomes are the basic unit of the three-dimensional structure of chromatin. It is now widely accepted that the positioning and occupancy of nucleosomes play important roles in fundamental genomic processes such as DNA transcription, replication and repair. Among the methods used to provide genome-wide nucleosomal positions and occupancy levels, MNase-seq has proven to be highly effective. Indeed, with this method, the nucleosomal landscapes of a variety of organisms have now been investigated, revealing both commonalities and differences. In this review, we first introduce the technical principles underlying MNase-seq, focusing on details essential to precisely resolve nucleosome positioning and occupancy. We then describe recent advances with this method, as well as future perspectives of its role in chromatin biology, with a particular focus of uncovering mechanistic insights of many disease process.

    nucleosome; chromatin structure; chromatin remodeling; next-generation sequencing (NGS); micrococcal nuclease

    2020-09-04;

    2020-10-18

    國家自然科學(xué)基金項目(編號:81972909)資助[Supported by the National Natural Science Foundation of China(No. 81972909)]

    鄧瑋杭,在讀碩士研究生,專業(yè)方向:系統(tǒng)生物醫(yī)學(xué)。E-mail: weihangdeng@sjtu.edu.cn

    李鑫輝,博士,助理研究員,研究方向:系統(tǒng)生物學(xué)與分子生物學(xué)。E-mail: xhli@sjtu.edu.cn

    10.16288/j.yczz.20-178

    2020/11/9 11:17:28

    URI: https://kns.cnki.net/kcms/detail/11.1913.R.20201106.1053.003.html

    (責(zé)任編委: 李海濤)

    猜你喜歡
    染色質(zhì)基因組位點
    染色質(zhì)開放性與動物胚胎發(fā)育關(guān)系的研究進展
    哺乳動物合子基因組激活過程中的染色質(zhì)重塑
    鎳基單晶高溫合金多組元置換的第一性原理研究
    上海金屬(2021年6期)2021-12-02 10:47:20
    牛參考基因組中發(fā)現(xiàn)被忽視基因
    CLOCK基因rs4580704多態(tài)性位點與2型糖尿病和睡眠質(zhì)量的相關(guān)性
    染色質(zhì)可接近性在前列腺癌研究中的作用
    “哺乳動物卵母細(xì)胞生發(fā)泡染色質(zhì)構(gòu)型的研究進展”一文附圖
    二項式通項公式在遺傳學(xué)計算中的運用*
    基因組DNA甲基化及組蛋白甲基化
    遺傳(2014年3期)2014-02-28 20:58:49
    有趣的植物基因組
    亚洲欧美日韩另类电影网站| 国产成人精品福利久久| 成人黄色视频免费在线看| 精品午夜福利在线看| 日本爱情动作片www.在线观看| 国产精品一区www在线观看| 精品人妻一区二区三区麻豆| 夫妻性生交免费视频一级片| 99久久中文字幕三级久久日本| 天美传媒精品一区二区| 丝袜在线中文字幕| 日韩,欧美,国产一区二区三区| 男的添女的下面高潮视频| 桃花免费在线播放| 国产精品熟女久久久久浪| 极品少妇高潮喷水抽搐| av.在线天堂| 少妇人妻久久综合中文| 大片免费播放器 马上看| 亚洲精品456在线播放app| 一级毛片aaaaaa免费看小| 精品熟女少妇av免费看| 有码 亚洲区| 99久国产av精品国产电影| 精品人妻一区二区三区麻豆| 看十八女毛片水多多多| 久久国产乱子免费精品| 亚洲av电影在线观看一区二区三区| 麻豆成人午夜福利视频| 亚洲国产毛片av蜜桃av| 国产高清不卡午夜福利| 国内揄拍国产精品人妻在线| 日本av手机在线免费观看| 人妻系列 视频| 两个人免费观看高清视频 | 国产黄片视频在线免费观看| a级毛色黄片| 午夜视频国产福利| 天天操日日干夜夜撸| 黑人高潮一二区| 性色avwww在线观看| 国产深夜福利视频在线观看| 黄色一级大片看看| 美女主播在线视频| 视频中文字幕在线观看| 婷婷色av中文字幕| 伦精品一区二区三区| freevideosex欧美| 丰满乱子伦码专区| 人妻少妇偷人精品九色| 精品久久久噜噜| 永久网站在线| 在线观看三级黄色| 国产亚洲精品久久久com| 91成人精品电影| 色网站视频免费| 桃花免费在线播放| 人妻人人澡人人爽人人| 97在线人人人人妻| .国产精品久久| 一区二区av电影网| 毛片一级片免费看久久久久| 三上悠亚av全集在线观看 | 欧美 亚洲 国产 日韩一| av黄色大香蕉| 国产av精品麻豆| 久久久久久久久大av| 王馨瑶露胸无遮挡在线观看| 99久久人妻综合| 国产精品一区二区在线观看99| 国内揄拍国产精品人妻在线| 99久久综合免费| 黑人巨大精品欧美一区二区蜜桃 | 国产精品不卡视频一区二区| 国产黄色免费在线视频| 狂野欧美激情性xxxx在线观看| 狂野欧美激情性bbbbbb| 国产成人午夜福利电影在线观看| 亚洲精品日韩在线中文字幕| a级毛片免费高清观看在线播放| 在线 av 中文字幕| 亚洲国产毛片av蜜桃av| 亚洲国产精品一区二区三区在线| 国产亚洲91精品色在线| 精品一区二区三区视频在线| 亚洲丝袜综合中文字幕| 中文天堂在线官网| 中文精品一卡2卡3卡4更新| 老司机亚洲免费影院| 大又大粗又爽又黄少妇毛片口| 国产永久视频网站| 久久久精品94久久精品| 久久人人爽人人片av| 熟女av电影| 如日韩欧美国产精品一区二区三区 | 午夜av观看不卡| 久久6这里有精品| 狠狠精品人妻久久久久久综合| 伦理电影大哥的女人| 搡女人真爽免费视频火全软件| 久久 成人 亚洲| 久久久久久久大尺度免费视频| 久久精品熟女亚洲av麻豆精品| 国产成人精品无人区| 精品久久国产蜜桃| 麻豆精品久久久久久蜜桃| 国产成人一区二区在线| 街头女战士在线观看网站| 色视频在线一区二区三区| 国产日韩欧美视频二区| 丝瓜视频免费看黄片| 人人妻人人澡人人爽人人夜夜| 精品人妻偷拍中文字幕| 在线观看三级黄色| 插阴视频在线观看视频| 国产精品蜜桃在线观看| 人人妻人人爽人人添夜夜欢视频 | 日韩制服骚丝袜av| 国产男女内射视频| 91精品国产国语对白视频| 亚洲人与动物交配视频| 国产亚洲av片在线观看秒播厂| 久久亚洲国产成人精品v| 国产免费一级a男人的天堂| 久久精品久久久久久噜噜老黄| 欧美少妇被猛烈插入视频| 亚洲综合精品二区| 晚上一个人看的免费电影| av视频免费观看在线观看| 欧美日韩综合久久久久久| 91午夜精品亚洲一区二区三区| 欧美另类一区| 亚洲av综合色区一区| 五月开心婷婷网| 女人精品久久久久毛片| 国产片特级美女逼逼视频| av视频免费观看在线观看| 国产一区二区三区av在线| 男男h啪啪无遮挡| 嘟嘟电影网在线观看| 97超碰精品成人国产| 国产永久视频网站| 免费不卡的大黄色大毛片视频在线观看| 亚洲色图综合在线观看| 国产深夜福利视频在线观看| 男女免费视频国产| 国产精品一区二区在线不卡| 国产精品久久久久久久久免| av天堂中文字幕网| 在线观看三级黄色| 免费在线观看成人毛片| 麻豆成人av视频| 国内揄拍国产精品人妻在线| 亚洲四区av| 少妇的逼水好多| 亚洲天堂av无毛| 18禁在线无遮挡免费观看视频| 精品久久久噜噜| 只有这里有精品99| 亚洲成人手机| 国产日韩欧美在线精品| 国产黄色免费在线视频| 久久国产亚洲av麻豆专区| 免费黄频网站在线观看国产| a级毛片在线看网站| 日本爱情动作片www.在线观看| 新久久久久国产一级毛片| 国产成人精品无人区| 国产精品女同一区二区软件| 蜜桃久久精品国产亚洲av| 成年美女黄网站色视频大全免费 | 日韩 亚洲 欧美在线| 在线精品无人区一区二区三| 亚洲怡红院男人天堂| 老司机亚洲免费影院| 人妻一区二区av| 九九久久精品国产亚洲av麻豆| 久久久久网色| 国产又色又爽无遮挡免| 久久精品国产亚洲av涩爱| 国产精品伦人一区二区| 亚洲国产毛片av蜜桃av| 成人亚洲精品一区在线观看| 亚洲天堂av无毛| 99九九在线精品视频 | 一级片'在线观看视频| 性色av一级| 黑丝袜美女国产一区| 国产国拍精品亚洲av在线观看| 丰满饥渴人妻一区二区三| 黑丝袜美女国产一区| 丰满乱子伦码专区| 自线自在国产av| 欧美精品人与动牲交sv欧美| 日日爽夜夜爽网站| 中文字幕久久专区| 成人黄色视频免费在线看| 国产精品久久久久成人av| 亚洲精品一二三| 草草在线视频免费看| 曰老女人黄片| 国产成人精品福利久久| 一边亲一边摸免费视频| 三上悠亚av全集在线观看 | 91精品国产九色| 亚洲欧美精品自产自拍| 综合色丁香网| 观看免费一级毛片| 深夜a级毛片| 高清欧美精品videossex| 高清视频免费观看一区二区| 国产成人91sexporn| 97精品久久久久久久久久精品| 欧美激情极品国产一区二区三区 | 桃花免费在线播放| 国产亚洲一区二区精品| 夜夜看夜夜爽夜夜摸| 青青草视频在线视频观看| 国产亚洲av片在线观看秒播厂| 亚洲精品,欧美精品| 午夜老司机福利剧场| 国产精品不卡视频一区二区| 日韩大片免费观看网站| 亚洲av综合色区一区| 大话2 男鬼变身卡| 啦啦啦中文免费视频观看日本| 久久狼人影院| 26uuu在线亚洲综合色| 曰老女人黄片| 国产黄片美女视频| 国产视频首页在线观看| 老女人水多毛片| 国产精品福利在线免费观看| 春色校园在线视频观看| 国产成人freesex在线| 日韩av免费高清视频| 中文字幕亚洲精品专区| 日韩一本色道免费dvd| 大片免费播放器 马上看| 国产成人午夜福利电影在线观看| 国产成人免费观看mmmm| 久久精品久久久久久久性| 婷婷色av中文字幕| 国产有黄有色有爽视频| 9色porny在线观看| 欧美精品高潮呻吟av久久| 成人黄色视频免费在线看| 下体分泌物呈黄色| 国产精品一区二区在线观看99| 我要看日韩黄色一级片| 乱码一卡2卡4卡精品| 欧美日韩视频精品一区| 女性被躁到高潮视频| 三级经典国产精品| 国产精品一区二区在线不卡| 在线观看一区二区三区激情| 精品久久久久久久久亚洲| a 毛片基地| 99精国产麻豆久久婷婷| 黄色怎么调成土黄色| 日韩精品免费视频一区二区三区 | 肉色欧美久久久久久久蜜桃| 国精品久久久久久国模美| 国产午夜精品一二区理论片| 97超视频在线观看视频| 成年av动漫网址| 亚洲一区二区三区欧美精品| 春色校园在线视频观看| 亚洲中文av在线| 成人18禁高潮啪啪吃奶动态图 | 人妻一区二区av| 国产国拍精品亚洲av在线观看| 乱人伦中国视频| 少妇人妻一区二区三区视频| 最近2019中文字幕mv第一页| 欧美三级亚洲精品| 国产免费视频播放在线视频| 人人妻人人看人人澡| 97超视频在线观看视频| 熟妇人妻不卡中文字幕| 国产精品无大码| 最近2019中文字幕mv第一页| 日本-黄色视频高清免费观看| 卡戴珊不雅视频在线播放| 国产毛片在线视频| 多毛熟女@视频| 久久久久久久亚洲中文字幕| 日韩不卡一区二区三区视频在线| 妹子高潮喷水视频| a级毛片在线看网站| 免费观看的影片在线观看| 久久久久久久久久久免费av| 黄色怎么调成土黄色| 久久毛片免费看一区二区三区| 亚洲一级一片aⅴ在线观看| 免费黄网站久久成人精品| 欧美日韩在线观看h| 日韩 亚洲 欧美在线| 久久人妻熟女aⅴ| 好男人视频免费观看在线| 人妻制服诱惑在线中文字幕| 久久ye,这里只有精品| 亚洲三级黄色毛片| 七月丁香在线播放| 99热这里只有精品一区| 深夜a级毛片| 国产日韩欧美视频二区| 国内少妇人妻偷人精品xxx网站| 免费观看的影片在线观看| 久久久国产精品麻豆| 亚洲精品久久午夜乱码| 97超碰精品成人国产| 精品久久国产蜜桃| 国产免费又黄又爽又色| 少妇人妻 视频| 婷婷色av中文字幕| 日产精品乱码卡一卡2卡三| 下体分泌物呈黄色| 国产免费福利视频在线观看| 国语对白做爰xxxⅹ性视频网站| 黄色配什么色好看| 久久综合国产亚洲精品| 欧美变态另类bdsm刘玥| 国产亚洲欧美精品永久| tube8黄色片| 亚洲怡红院男人天堂| 国产欧美日韩一区二区三区在线 | 一区二区三区四区激情视频| 人妻一区二区av| 搡老乐熟女国产| 精品视频人人做人人爽| 亚洲精品色激情综合| 丁香六月天网| 国产成人免费观看mmmm| 久久精品国产鲁丝片午夜精品| 日韩精品有码人妻一区| 男人狂女人下面高潮的视频| av播播在线观看一区| 51国产日韩欧美| 成人毛片60女人毛片免费| 99久国产av精品国产电影| 在线天堂最新版资源| 大片免费播放器 马上看| 女人久久www免费人成看片| av在线播放精品| 精品人妻熟女毛片av久久网站| 国产高清国产精品国产三级| 午夜免费男女啪啪视频观看| 亚洲欧洲日产国产| 久久精品熟女亚洲av麻豆精品| 国产亚洲91精品色在线| 乱码一卡2卡4卡精品| 亚洲激情五月婷婷啪啪| 有码 亚洲区| 七月丁香在线播放| 国产老妇伦熟女老妇高清| 国产精品麻豆人妻色哟哟久久| 蜜臀久久99精品久久宅男| 久久韩国三级中文字幕| 在线观看免费视频网站a站| 最新的欧美精品一区二区| 色视频在线一区二区三区| 亚洲图色成人| 汤姆久久久久久久影院中文字幕| 丝袜脚勾引网站| 欧美xxⅹ黑人| 五月天丁香电影| 国产精品一区二区性色av| 国产精品熟女久久久久浪| 制服丝袜香蕉在线| 在线 av 中文字幕| 天天操日日干夜夜撸| 亚洲欧美中文字幕日韩二区| 边亲边吃奶的免费视频| 国产美女午夜福利| 男女边摸边吃奶| 亚洲精品第二区| 男女啪啪激烈高潮av片| 国产中年淑女户外野战色| 汤姆久久久久久久影院中文字幕| 青春草国产在线视频| 国产精品福利在线免费观看| 成人二区视频| 少妇被粗大猛烈的视频| 亚洲不卡免费看| 日韩制服骚丝袜av| 中文在线观看免费www的网站| 亚洲一级一片aⅴ在线观看| 99热这里只有是精品50| 久久久国产一区二区| 国产伦精品一区二区三区四那| 丰满乱子伦码专区| 日日摸夜夜添夜夜添av毛片| 青青草视频在线视频观看| 国产在线视频一区二区| 在线观看免费高清a一片| .国产精品久久| 插阴视频在线观看视频| 制服丝袜香蕉在线| 人妻少妇偷人精品九色| 最新的欧美精品一区二区| 国产亚洲一区二区精品| 国产熟女午夜一区二区三区 | 国产淫语在线视频| 激情五月婷婷亚洲| 国产毛片在线视频| 一区二区三区四区激情视频| 人人妻人人澡人人爽人人夜夜| 能在线免费看毛片的网站| 一级毛片黄色毛片免费观看视频| 亚洲av免费高清在线观看| 中文天堂在线官网| 国产日韩欧美视频二区| 哪个播放器可以免费观看大片| 久久久久久久国产电影| 久久久久久伊人网av| 国产成人精品一,二区| 一级a做视频免费观看| 国产午夜精品一二区理论片| 天堂8中文在线网| 国产在线免费精品| 午夜91福利影院| 大陆偷拍与自拍| 嫩草影院入口| 成年美女黄网站色视频大全免费 | 日韩欧美一区视频在线观看 | 亚洲av二区三区四区| 亚洲精品国产av成人精品| 成人毛片60女人毛片免费| 在线免费观看不下载黄p国产| 国产精品成人在线| 三上悠亚av全集在线观看 | 国产又色又爽无遮挡免| 欧美日韩视频高清一区二区三区二| 九草在线视频观看| 性色avwww在线观看| 人妻夜夜爽99麻豆av| 亚洲精品第二区| 国产 一区精品| 午夜精品国产一区二区电影| 亚洲伊人久久精品综合| 观看免费一级毛片| 男女边吃奶边做爰视频| 久久久久久久久久成人| 18禁在线无遮挡免费观看视频| 久久人妻熟女aⅴ| 国产精品一区二区性色av| 丝袜脚勾引网站| 亚洲真实伦在线观看| 亚洲av不卡在线观看| 草草在线视频免费看| 国产色爽女视频免费观看| 国产精品欧美亚洲77777| 麻豆成人av视频| 男男h啪啪无遮挡| 永久网站在线| 99热这里只有是精品在线观看| 最黄视频免费看| 中文字幕人妻熟人妻熟丝袜美| 在线观看三级黄色| 人妻一区二区av| 嫩草影院入口| av网站免费在线观看视频| 性色av一级| 日韩制服骚丝袜av| 秋霞伦理黄片| 亚洲,一卡二卡三卡| 你懂的网址亚洲精品在线观看| 一本一本综合久久| 亚洲欧美精品自产自拍| av不卡在线播放| 久久久久久久久久久免费av| 国产免费视频播放在线视频| 人人妻人人澡人人爽人人夜夜| 亚洲婷婷狠狠爱综合网| 黄色一级大片看看| 我要看日韩黄色一级片| 亚洲人成网站在线观看播放| 国产精品国产三级国产av玫瑰| 国产一级毛片在线| 免费高清在线观看视频在线观看| a 毛片基地| 边亲边吃奶的免费视频| 2021少妇久久久久久久久久久| .国产精品久久| 国产免费视频播放在线视频| 亚洲在久久综合| 18禁在线播放成人免费| 少妇的逼水好多| a级毛片免费高清观看在线播放| 亚洲丝袜综合中文字幕| 欧美bdsm另类| 精品国产一区二区久久| 性色av一级| 国产亚洲5aaaaa淫片| 国产精品国产三级国产专区5o| 美女大奶头黄色视频| 一本大道久久a久久精品| 免费观看在线日韩| 国产精品久久久久久精品电影小说| 22中文网久久字幕| 国产中年淑女户外野战色| 男人舔奶头视频| 另类精品久久| 久久精品国产自在天天线| 中文字幕制服av| 美女国产视频在线观看| 精品人妻一区二区三区麻豆| 亚洲欧美成人综合另类久久久| 日日啪夜夜爽| 亚洲精品色激情综合| 欧美精品高潮呻吟av久久| 国产精品人妻久久久影院| 夜夜爽夜夜爽视频| 国产男女内射视频| 欧美变态另类bdsm刘玥| 噜噜噜噜噜久久久久久91| 国产伦在线观看视频一区| 日本av手机在线免费观看| 日韩成人伦理影院| 免费观看性生交大片5| 精品酒店卫生间| 午夜福利视频精品| 老司机影院毛片| 能在线免费看毛片的网站| 肉色欧美久久久久久久蜜桃| 在现免费观看毛片| 国产永久视频网站| 少妇熟女欧美另类| 国产精品成人在线| 极品教师在线视频| 精品国产乱码久久久久久小说| 蜜桃久久精品国产亚洲av| 色婷婷av一区二区三区视频| 涩涩av久久男人的天堂| 国产一区亚洲一区在线观看| 久久精品熟女亚洲av麻豆精品| 国产精品欧美亚洲77777| 国产成人午夜福利电影在线观看| 中国国产av一级| 国产成人freesex在线| 91精品伊人久久大香线蕉| 99精国产麻豆久久婷婷| 99久国产av精品国产电影| av网站免费在线观看视频| 久久精品国产亚洲网站| 国产精品成人在线| 三上悠亚av全集在线观看 | 亚洲精品第二区| 国产精品成人在线| 纯流量卡能插随身wifi吗| 精品一区二区免费观看| 国产精品麻豆人妻色哟哟久久| 国产国拍精品亚洲av在线观看| 少妇被粗大的猛进出69影院 | 在线精品无人区一区二区三| 国精品久久久久久国模美| 国产高清不卡午夜福利| 老司机亚洲免费影院| 亚洲精品亚洲一区二区| 天堂俺去俺来也www色官网| 成人黄色视频免费在线看| 亚洲美女搞黄在线观看| 亚州av有码| 久热这里只有精品99| 亚洲情色 制服丝袜| 精品少妇黑人巨大在线播放| 国产视频首页在线观看| 一个人看视频在线观看www免费| 一区二区三区四区激情视频| 国产男女超爽视频在线观看| 人妻人人澡人人爽人人| 你懂的网址亚洲精品在线观看| av有码第一页| 特大巨黑吊av在线直播| 少妇被粗大猛烈的视频| 97在线人人人人妻| 国产亚洲欧美精品永久| 欧美老熟妇乱子伦牲交| 亚洲国产精品一区三区| 岛国毛片在线播放| 一本久久精品| av网站免费在线观看视频| 久久精品国产亚洲av天美| 91久久精品国产一区二区三区| 高清午夜精品一区二区三区| 极品教师在线视频| 成人特级av手机在线观看| 伦理电影大哥的女人| 2021少妇久久久久久久久久久| 一级,二级,三级黄色视频| 亚洲精品456在线播放app| 国内揄拍国产精品人妻在线| 久久久久久久亚洲中文字幕| 亚洲欧洲日产国产| 欧美+日韩+精品| 能在线免费看毛片的网站| 日韩av不卡免费在线播放| 亚洲精品国产av蜜桃| 97超视频在线观看视频| 日韩成人伦理影院| 国产精品久久久久久av不卡| 你懂的网址亚洲精品在线观看| 97精品久久久久久久久久精品| 国产 一区精品| 高清在线视频一区二区三区| 亚洲欧洲日产国产| 日韩,欧美,国产一区二区三区| 亚洲va在线va天堂va国产| 91久久精品电影网| 超碰97精品在线观看| av又黄又爽大尺度在线免费看| 国产在线免费精品| 午夜福利视频精品| 国产黄频视频在线观看|