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

    真核生物基因組長內含子遞歸剪接事件的分子機制

    2019-02-28 07:05:54魏金川徐添翼吳靜宋曉峰
    遺傳 2019年2期
    關鍵詞:套索內含子堿基

    魏金川,徐添翼,吳靜,2,宋曉峰

    ?

    真核生物基因組長內含子遞歸剪接事件的分子機制

    魏金川1,徐添翼1,吳靜1,2,宋曉峰1

    1. 南京航空航天大學自動化學院,南京 211106 2. 南京醫(yī)科大學生物醫(yī)學工程與信息學院,南京 211166

    在高等真核生物基因組轉錄過程中,一次剪接可完成短內含子的去除,而較長內含子(>10 kb)則需通過多次剪接方可去除。多次剪接去除長內含子的過程通常被稱為遞歸性剪接。已有研究表明,遞歸性剪接事件與諸多生物學過程及疾病的發(fā)生發(fā)展有著密切的聯(lián)系。近年來,關于遞歸性剪接的研究越來越多,研究者已經在果蠅()和多種脊椎動物基因組轉錄過程中發(fā)現(xiàn)了遞歸剪接事件,通過不同的生物信息學方法找到了多個遞歸剪接位點并進行了實驗驗證。目前國際上對遞歸性剪接的研究主要集中在遞歸剪接過程、剪接位點識別及其對生物學過程的影響等方面。本文針對真核生物基因組轉錄過程中遞歸剪接事件的分子機制和國內外研究現(xiàn)狀進行了綜述,旨在為深入理解RNA剪接過程分子機制提供參考。

    遞歸剪接;套索結構;RS位點;Total RNA-seq

    真核生物基因組轉錄過程是以DNA雙鏈中一條鏈為模板,經堿基互補配對得到前體mRNA,后通過對內含子與外顯子的剪接加工得到成熟的mRNA。內含子是基因編碼區(qū)中的非編碼序列,在生物體內不編碼蛋白,因此在前體mRNA轉錄加工過程中會被剪接體(spliceosome)識別并去除。起初研究者認為內含子不參與基因的表達調控,但近幾年隨著二代測序技術的發(fā)展和生物信息學工具的開發(fā),人們發(fā)現(xiàn)內含子在轉錄及基因表達調控中具有重要的生物學功能,比如:有些內含子含有增強子、啟動子或者其他的順式作用元件,對基因的表達起調控作用;內含子的剪接增加了mRNA的穩(wěn)定性;內含子還可能發(fā)生選擇性剪接事件[1~5]。

    前體mRNA的剪接過程中,內含子通過兩個外顯子的連接被完全切除,這個過程需要剪接體介導。前體mRNA的剪接需要兩步實現(xiàn):第一步,靠近內含子3′端分支位點的A殘基的2′-OH對5′端的磷酸基團進行親核攻擊,形成2′→5′磷酸二酯鍵連接的套索結構;第二步,被剪接的外顯子3′端核苷酸的-OH進攻內含子3′末端的磷酸基團,使得內含子3′端在剪接位點處斷開,釋放套索結構,兩個外顯子連接[6~8]。此過程是直接去除整個內含子的,屬于一步剪接[9],但是有一些較長的內含子是經過了多步剪接來去除的。1999年,Hatton等[10]在黑腹果蠅()基因中發(fā)現(xiàn),其第1個內含子分3步去除,并將這一剪接過程命名為遞歸剪接(recu-rsive splicing) (圖1)。內含子中被剪接的核酸位點稱為遞歸性剪接位點(RS-site)。近年來研究人員發(fā)現(xiàn),除果蠅以外,遞歸剪接在其他真核生物基因組中也廣泛存在,并與疾病的發(fā)生發(fā)展有重要關系。本文對真核生物基因組中長內含子遞歸剪接分子機制的國內外研究現(xiàn)狀進行了綜述,以期為RNA的剪接研究提供參考,同時為疾病的診斷和治療提供新的 思路。

    圖1 遞歸剪接示意圖

    含有兩個剪接位點的內含子分步去除過程。exon1與exon2分別表示兩個相鄰外顯子,中間內含子有兩個遞歸剪接位點,灰色箭頭指向剪接位點,一共經過了3步拼接在一起。

    1 遞歸剪接的特征

    1.1 遞歸剪接的過程

    通常內含子5′端剪接位點被稱為供體位點(donor site),3′端剪接位點被稱為受體位點(acceptor site)[11,12]。研究表明,內含子起始5′端兩對堿基和結尾3′端的兩對堿基最為保守,在所有的剪接位點中有99.24%的位點為GU-AG,0.7%的位點為GC-AG,0.05%的位點為AT-AC[13]。除了這兩對保守堿基外,它們附近的堿基雖然在不同的物種之間存在差異,但在同一物種中通常具有保守性,如在脊椎動物中5′剪接端堿基通常為AG|GUAAGU[14]。由堿基序列可以看出供體與受體間不存在堿基互補配對的可能,所以它們并不是通過堿基互補的方式拼接到一起的。外顯子轉錄是從5′端到3′端進行的,遞歸剪接就是把剪接位點處默認為含有0個核苷酸的外顯子,所以遞歸剪接位點通常為AG|GU。

    遞歸剪接的具體過程大致如下[15~17]:(1)首先內含子3′端上游10~50 bp處的特定腺苷的2′羥基攻擊5′剪接位點(圖2A),形成5′-2′磷酸二酯鍵,從而構成套索結構,進而釋放5′外顯子(圖2B);(2) 5′外顯子的3′羥基攻擊內含子3′端最近的遞歸剪接位點(RS sites),釋放套索結構,與此同時,5′外顯子的3′-羥基與下游外顯子的5′端連接形成新的磷酸二酯鍵(圖2C);(3)重復前面兩個步驟,直到全部內含子去除。

    1.2 遞歸剪接位點的特征

    針對可變剪接的研究發(fā)現(xiàn)內含子5′和3′端堿基幾乎都是GU和AG,并且內含子位點為AGGU的比例高達99.24%,所以人們對遞歸剪接位點的研究主要也主要針對AGGU型。此外,遞歸剪接在形成套索結構的過程中,在剪接位點上游第5和第6位置處的堿基基本上都為“U”,這也是尋找遞歸剪接位點的主要依據(jù)。而Kelly等[18]研究發(fā)現(xiàn)在人類()全基因組中并非所有遞歸剪接位點都是AGGU,而是AGMN(M、N分別表示A、G、C、U 4種堿基中的一種),其中AGGN的概率為44%,AGGU的概率僅為14% (表1)。因此,對于非經典(剪接位點不是AGGU)遞歸剪接位點的識別已成為當前研究的熱點問題。

    圖2 遞歸剪接套索結構形成與切除示意圖

    A:內含子3′端上游10~50 bp處的特定腺苷的2′羥基攻擊5′剪接位點;B:構成套索結構,釋放5′外顯子;C:釋放套索結構外顯子與遞歸剪接位點相互連接。

    表1 遞歸剪接位點AGMN中M、N為各堿基概率統(tǒng)計表

    遞歸性剪接位點AGMN,M與N分別對應4種堿基的概率[18]。M/N分別表示剪接位點AGMN中的M/N兩個堿基;A、C、G、U為4種堿基。

    1.3 遞歸剪接位點的識別

    Duff等[15]對黑腹果蠅進行了轉錄組測序,并將得到的Total RNA-seq數(shù)據(jù)用Tophat[19]進行了處理。從得到的reads與全基因組匹配結果中,他們發(fā)現(xiàn)基因的轉錄是多個RNA聚合酶同時進行,所以靠近5′外顯子處的轉錄酶更多(圖3A)。存在遞歸剪接現(xiàn)象的內含子在reads的分布上存在著波峰與波谷,并且整體呈現(xiàn)鋸齒狀,這種分布模式體現(xiàn)了遞歸剪接去除內含子的機理,其理想的reads密度圖為倒三角狀(圖3B)。遞歸剪接通常存在于大于10 kb的內含子中,人類基因組中此類內含子較多,長度大于24 kb的內含子有8000個以上;大于50 kb的內含子有3400多個;而超過100 kb的內含子都有1200個以上[20,21]。Sibley等[22]對人腦組織Total RNA-seq數(shù)據(jù)進行了分析,從中發(fā)現(xiàn)了8個遞歸剪接位點,這些位點所在的內含子長度都超過150 kb。Duff等[15]用Tophat將果蠅Total RNA-seq數(shù)據(jù)與果蠅全基因組(modENCODE注釋文件)相匹配,發(fā)現(xiàn)含有剪接位點的內含子長度都大于2 kb。

    近來,遞歸剪接位點的識別主要基于Total RNA-Seq數(shù)據(jù),輔以對應物種的基因組注釋數(shù)據(jù)及基因組序列信息。通常,遞歸剪接位點的識別主要分為以下4步(圖4):(1) 通過FastQC、Fastx-toolkit分別對其質量檢測、控制;(2)利用Tophat進行reads映射,得到Bam文件與Junction文件;(3)依據(jù)內含子長度、剪接位點保守性對Junction位點進行篩選;(4)將Bam文件轉換成Sam文件,并與Junction文件相結合,進而通過生物信息學方法識別遞歸剪接位點。目前,常將周邊reads分布呈鋸齒狀的位點判定為潛在的遞歸剪接位點。

    圖3 Total RNA-seq數(shù)據(jù)匹配全基因組統(tǒng)計得到的結果圖

    A:含遞歸剪接的基因組轉錄示意圖。AG/GU為遞歸剪接位點,曲線表示RNA轉錄情況,橢圓形表示RNA聚合酶Ⅱ,代表一條DNA鏈在同時進行多次轉錄。B:Total RNA-seq的reads密度圖。

    圖4 遞歸剪接位點識別流程圖

    矩形框中的是數(shù)據(jù)處理過程、用到的文件、得到的結果等。

    2 遞歸剪接研究的進展

    2.1 果蠅基因組中遞歸剪接的研究

    Hatton等[10]于1998年在黑腹果蠅基因中首次發(fā)現(xiàn)了遞歸性剪接現(xiàn)象,在基因的第1個外顯子和第2個外顯子之間發(fā)現(xiàn)了新的剪接位點,并且在3′剪接位點的下游的堿基序列為GTAAGA。研究者對基因進行了深入研究,發(fā)現(xiàn)該基因第1個內含子長73 kb,該內含子的去除經歷了3次遞歸剪接[24]。

    Conklin等[25]于2005年通過套索結構識別遞歸剪接位點,并于基因的第一個內含子中成功識別出Hatton等[10]找到的兩個遞歸剪接位點。針對套索結構RNA分子半衰期通常較短,且轉錄的過程不能準確把握等問題,他們通過控制PCR技術對分支連接處2′→5′磷酸二酯鍵的特異性擴增,提高了RNA套索結構分析的靈敏度。

    2015年,Duff 等[15]以黑腹果蠅為材料,對其進行實驗處理,通過二代測序技術得到對應的Total RNA-seq數(shù)據(jù),然后通過Tophat/Bowtie軟件比對得到Junction位點信息,結合reads在Junction位點附近的分布信息,最終篩選出197個潛在的遞歸剪接位點。這些遞歸剪接位點分布于115個基因的130個內含子中,他們挑選了其中14個基因用作實驗驗證,證實這些基因含有24個遞歸剪接位點。對存在遞歸性剪接的基因(115 個)和內含子統(tǒng)計分析發(fā)現(xiàn),其中有100個基因只存在一個遞歸性剪接內含子,剩下的15個基因含有兩個遞歸性剪接內含子。此外,這130個內含子的剪接位點個數(shù)都在1~6之間。遞歸性剪接內含子的堿基長度為11 341~132 736 bp,內含子的平均長度為45 164 bp。進一步分析發(fā)現(xiàn),存在遞歸剪接的內含子長度都相對較長,但是并不是所有的較長內含子都存在遞歸剪接,在較長的內含子中只有大約6%存在遞歸性剪接。他們還發(fā)現(xiàn)U2AF蛋白減少會導致含有遞歸剪接的Pre-mRNA加工速率降低。

    Pre-mRNA在內含子轉錄完成幾秒后就會剪切掉內含子,遞歸剪接是內含子中間片段的剪接。所以,2018年Pai等[26]利用4sU對果蠅S2細胞進行標記,分別獲取轉錄5 min、10 min、20 min的Total RNA-seq數(shù)據(jù),從中發(fā)現(xiàn)了541個遞歸剪接位點,其中實驗證明的有119個,這些位點所屬內含子長度大部分都超過40 kb。

    2.2 脊椎動物基因組中遞歸剪接的研究

    除果蠅以外,研究者還對多種脊柱動物進行了遞歸剪接的研究,包括人、大鼠()、斑馬魚()等。

    2012年,Taggart 等[27]開發(fā)了一個基于RNA套索結構識別內含子3′剪接位點的生物信息學方法,其假定套索結構都位于3′剪接位點附近。利用該方法,他們對人類12個組織的轉錄組測序數(shù)據(jù)進行了處理,共篩選出861個內含子3′剪接位點。其中,2118個reads經匹配映射到這些剪接位點附近的套索結構中,而屬于遞歸性剪接套索結構的reads有70個(3%)。

    2015年,Kelly等[18]對基因轉錄過程中的內含子遞歸剪接事件進行了探究。基因全長221 kb,其中第1個內含子長134 kb。研究人員用人臍靜脈內皮細胞(human umbilical vein endoth-elial cells, HUVEC)作為實驗材料,利用腫瘤壞死因子(TNFα)對HUVEC進行刺激,并通過二代測序技術(high-throughput sequencing)得到相關的Total RNA- seq數(shù)據(jù)。接著,將reads匹配到全基因組(hg19),發(fā)現(xiàn)reads匹配的位置主要分為3類情況:第一類為直接整體匹配到單一外顯子上;第二類匹配到exon-exon上,即分別匹配到上游外顯子的3′ 端部分堿基和下游外顯子的5′端部分堿基;第三類為exon-intron,即一部分匹配到上游外顯子上,而另一部分匹配到相鄰內含子的中間部分。其中第三類reads能反映出內含子存在遞歸剪接與可變剪接的可能。最后研究者用生物信息學的方法對套索結構進行篩選,從中發(fā)現(xiàn)基因的第1個內含子存在遞歸性剪接現(xiàn)象?;蛟谌四氺o脈內皮細胞中基本不表達,而TNFα對基因的轉錄有著很強的刺激作用,所以研究人員通過熒光原位雜交技術(fluorescencehybridization, FISH)對基因表達情況進行檢測,發(fā)現(xiàn)較長的內含子并不是轉錄結束后直接切除的,而是邊轉錄邊剪接的。對于基因的熒光標記發(fā)現(xiàn),剪切部分的水解時間占整個內含子轉錄時間的1/15;它每次剪切大約經歷了15 min,其平均半衰期大約為4.2 min。研究表明,在基因遞歸剪接過程中,內含子片段的水解過程與形成套索切除的過程是同步進行的;在較長的內含子中,轉錄的速度大約3 kb/min。研究人員使用CRISPR-Cas9技術[28]對遞歸剪接的分支點(branch point)進行突變與刪除,發(fā)現(xiàn)這樣會大大降低Pre-mRNA到成熟mRNA的加工速率。并且,對任意一個遞歸剪接位點(RS-site)進行基因突變或基因敲除處理,基因的表達量都會降低35%~50%,表明遞歸剪接的分支點和RS位點在對應基因的表達過程中有著重要的作用。同時研究人員還發(fā)現(xiàn),在人體靜脈細胞中含有多個剪接位點的基因在某一個剪接位點失效后,會通過下一個剪接位點或者其他的補償機制進行處理,雖然效率降低,但是還可以表達[18]。因此,遞歸剪接位點突變后,應同樣存在某種補償機制,使得基因仍能正常轉錄。

    2015年,Sibley等[22]對脊椎動物神經元細胞進行了RNA深度測序,分析發(fā)現(xiàn)在脊椎動物神經元中存在大量的遞歸剪接現(xiàn)象。由此,他們基于對人類腦組織中遞歸剪接的研究建立了一個雙重剪接的模型,該模型大體對AGGU為剪接位點的現(xiàn)象進行了研究,發(fā)現(xiàn)遞歸剪接位點上游內含子去除后,供體部分的外顯子與剪接后的受體首端的GU相互連接起來,形成新的外顯子-內含子-外顯子復合物,在進行下次剪接過程中,有可能本次外顯子-內含子會保持連在一起,最終的產物會有兩種,一種是外顯子-外顯子(exon-exon),另外一種為外顯子-遞歸剪接中間小段內含子(exon-intron)這兩種產物。雖然該模型還不夠完善,但也說明了人類細胞中mRNA亞型具有多樣性。Emmett[29]發(fā)現(xiàn)在人類腦組織中,存在遞歸剪接現(xiàn)象的7個基因與神經系統(tǒng)疾病密切有關。

    以上研究都是以組織的Total RNA-seq數(shù)據(jù)為材料進行分析實驗。2018年Hayashi等[30]開發(fā)了一種對單細胞Total RNA-seq進行測序的方法,即RamDA-seq。與其他方法相比,這種測序方法對non- poly(A) RNA敏感性更高。Hayashi等[30]用RamDA- seq對不同時期的小鼠胚胎干細胞進行實驗分析,在長度大于150 kb的內含子中發(fā)現(xiàn)了207個潛在的遞歸剪接位點。2018年Pai等[31]對果蠅的30 000個內含子進行研究,發(fā)現(xiàn)存在遞歸剪接的內含子也具有半衰期。此外,Srndic 等[32]在人類多個組織中檢測到12 000個遞歸剪接事件,并發(fā)現(xiàn)遞歸剪接對無義介導的RNA降解(nonsense-mediated decay)可能起著反作用。

    3 遞歸剪接研究的意義

    隨著RNA測序技術的發(fā)展,越來越多的內含子被發(fā)現(xiàn)在基因表達調控中發(fā)揮著重要的作用,因此對遞歸剪接的深入研究有助于理解內含子調控功能的分子機制[33]。全基因組研究發(fā)現(xiàn)機體可以通過改變RNA聚合酶Ⅱ的合成速度來控制轉錄的速率[34],堿基突變或缺失容易造成各類疾病[35,36]。人類基因組除了必要的編碼信息外,還有許多必要的剪接序列,其內部位點的突變或者缺失經常會導致疾病的發(fā)生[37~39]。研究人員還發(fā)現(xiàn)在剪接位點發(fā)生突變或者刪除時,成熟mRNA的生成速率會大大降低。

    在人類基因組中,遞歸剪接位點比其他位置更具有保守性,并且廣泛存在[18,22,29,40,41]。研究發(fā)現(xiàn),人體內的突變主要來自于不編碼蛋白的區(qū)域,占總突變事件的90%以上,且超過40%的突變來自于內含子區(qū)域[35]。剪接過程中的拼接錯誤是導致疾病的重要原因,遞歸剪接位點突變可能對某些疾病的研究具有重要意義[42]。目前,人們已對神經元[29]及內皮細胞[18]中的遞歸剪接事件進行了研究,發(fā)現(xiàn)剪接位點與神經系統(tǒng)疾病如帕金森綜合征、循環(huán)系統(tǒng)疾病如視網(wǎng)膜硬化或者高血壓等病變的關聯(lián)[43]。

    遞歸性剪接對生物體的生長發(fā)育具有重要的調節(jié)作用。例如,剪接位點的選擇通過順式作用元件[44]和反式作用因子共同控制,例如SR蛋白、hnRPNs等蛋白質。這些因素通過影響U1 snRNP,U2輔助因子和U2 snRNP剪接體的裝配來調控mRNA的序列,從而控制生物體的生長發(fā)育[45]。

    遞歸性剪接還有待進一步深入研究,目前尚待解決的問題有:(1)內含子遞歸剪接的機制,遞歸剪接位點間的片段是按前后順序依次切除,還是另有其他調控方式;(2)遞歸剪接事件的發(fā)生是否與生物發(fā)育的不同階段有關;(3)遞歸剪接位點突變與疾病間的關系。以上問題的研究必將有助于進一步理解RNA分子剪接機制,為疾病相關研究提供有價值的科研線索。

    4 結語與展望

    隨著RNA測序技術的發(fā)展與研究的深入,起初認為內含子是單步剪接,現(xiàn)在已發(fā)現(xiàn)部分長內含子存在多步遞歸剪接。目前,研究者已在人類、果蠅等多個真核生物基因組中發(fā)現(xiàn)了遞歸剪接現(xiàn)象,找到了多個遞歸剪接位點。雖然近年來對遞歸剪接的研究發(fā)展迅速,但是遞歸剪接位點突變或缺失造成的影響尚未可知,仍需要人們更進一步的探索。隨著遞歸剪接位點的不斷發(fā)現(xiàn),遞歸剪接位點突變或缺失在人類疾病發(fā)生發(fā)展中的功能作用將越來越清楚,為疾病的預防和治療提供新的科學線索。

    三代測序技術的迅猛發(fā)展對真核生物遞歸剪接的研究具有重要的推動作用。與二代測序技術相比,三代測序技術具有測序速度更快、測序長度激增、可直接測RNA序列等特點[46]。三代測序技術的測序長度由二代測序的上百堿基到現(xiàn)在的幾千堿基,因此其對基因長內含子的遞歸剪接事件的識別分析提供了更好的技術手段,相信未來會有更多有價值的遞歸剪接的生物學機制被發(fā)現(xiàn)。

    [1] Zhang J, Sun X, Qian Y, Maquat LE. Intron function in the nonsense-mediated decay of beta-globin mRNA: indications that pre-mRNA splicing in the nucleus can influence mRNA translation in the cytoplasm.,, 1998, 4(7): 801–815.

    [2] Hoffman MM, Birney E. Estimating the neutral rate of nucleotide substitution using introns., 2007, 24(2): 522–531.

    [3] Chevalier BS, Stoddard BL. Homing endonucleases: structural and functional insight into the catalysts of intron/intein mobility., 2001, 29(18): 3757–3774.

    [4] Zhang KH. The functions and applications of introns., 2012, 39(7): 80–83.張開慧. 內含子的功能及應用. 中國畜牧獸醫(yī), 2012, 39(7): 80–83.

    [5] Zhao JX, Wang F, Xu ZR, Fan YM. The epigenetic effect on pre-mRNA alternative splicing., 2014, 36(3): 248–255.趙金璇, 王芳, 徐崢嶸, 范怡梅. 表觀遺傳調控pre- mRNA的選擇性剪接. 遺傳, 2014, 36(3): 248–255.

    [6] Newman A. RNA splicing., 1998, 8(25): R903.

    [7] Yeo G, Burge CB. Maximum entropy modeling of short sequence motifs with applications to RNA splicing signals., 2004, 11(2–3): 377–394.

    [8] Li YI,van de Geijn B,Raj A,Knowles DA, Petti AA,Golan D, Gilad Y, Pritchard JK. RNA splicing is a primary link between genetic variation and disease., 2016, 352(6285): 600–604.

    [9] Hang J, Wan R, Yan C, Shi Y. Structural basis of pre-mRNA splicing., 2015, 349(6253): 1191– 1198.

    [10] Hatton AR, Subramaniam V, Lopez AJ. Generation of alternative Ultrabithorax isoforms and stepwise removal of a large intron by resplicing at exon-exon junctions., 1998, 2(6): 787–796.

    [11] Marashi SA, Goodarzi H, Sadeghi M, Eslahchi C, Pezeshk H. Importance of RNA secondary structure information for yeast donor and acceptor splice site predictions by neural networks., 2006, 30(1): 50–57.

    [12] Rautmann G, Matthes HW, Gait MJ, Breathnach R. Synthetic donor and acceptor splice sites function in an RNA polymerase B (II) transcription unit., 1984, 3(9): 2021–2028.

    [13] Burset M, Seledtsov IA, Solovyev VV. Analysis of canonical and non-canonical splice sites in mammalian genomes., 2000, 28(21): 4364–4375.

    [14] Schreier P, Aichinger C. Splicing as target for identifying new active substances. 2003.

    [15] Duff MO, Olson S, Wei X, Garrett SC, Osman A, Bolisetty M, Plocik A, Celniker SE, Graveley BR. Genome-wide identification of zero nucleotide recursive splicing in drosophila., 2015, 521(7552): 376–379.

    [16] Lopez PJ, Séraphin B. Uncoupling yeast intron recognition from transcription with recursive splicing., 2000, 1(4): 334–339.

    [17] Suzuki H, Kameyama T, Ohe K, Tsukahara T, Mayeda A. Nested introns in an intron: evidence of multi-step splicing in a large intron of the human dystrophin pre-mRNA., 2013, 587(6): 555–561.

    [18] Kelly S, Georgomanolis T, Zirkel A, Diermeier S, O’Reilly D, Murphy S, L?ngst G, Cook PR, Papantonis A. Splicing of many human genes involves sites embedded within introns., 2015, 43(9): 4721– 4732.

    [19] Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq., 2009, 25(9): 1105–1111.

    [20] Bradnam KR, Korf I. Longer first introns are a general property of eukaryotic gene structure., 2008, 3(8): e3093.

    [21] Shepard S, Mccreary M, Fedorov A. The peculiarities of large intron splicing in animals., 2009, 4(11): e7853.

    [22] Sibley CR, Emmett W, Blazquez L, Faro A, Haberman N, Briese M, Trabzuni D, Ryten M, Weale ME, Hardy J, Modic M, Curk T, Wilson SW, Plagnol V, Ule J. Recursive splicing in long vertebrate genes., 2015, 521(7552): 371.

    [23] Pulyakhina I, Gazzoli I, ’tHoen PA, Verwey N, den Dunnen JT, Aartsma-Rus A, Laros JF. SplicePie: a novel analytical approach for the detection of alternative, non-sequential and recursive splicing., 2015, 43(22): 11068.

    [24] Burnette JM, Miyamoto-Sato E, Schaub MA, Conklin J, Lopez AJ. Subdivision of large introns inby recursive splicing at nonexonic elements., 2005, 170(2): 661–674.

    [25] Conklin JF, Goldman A, Lopez AJ. Stabilization and analysis of intron lariats., 2005, 37(4): 368–375.

    [26] Pai AA, Paggi J, Adelman K, Burge CB. Numerous recursive sites contribute to accuracy of splicing of long introns in flies., 2018, 14(8): e1007588.

    [27] Taggart AJ, Desimone AM, Shih JS, Filloux ME, Fairbrother WG. Large-scale mapping of branchpoints in human pre-mRNA transcripts., 2012, 19(7): 719–721.

    [28] Ran FA, Hsu PD, Wright J, Agarwala V, Scott DA, Zhang F. Genome engineering using the CRISPR-Cas9 system., 2013, 8(11): 2281–2308.

    [29] Emmett WA. Computational identification of regulatory features affecting splicing in the human brain. 2016.

    [30] Hayashi T, Ozaki H, Sasagawa Y, Umeda M, Danno H, Nikaido I. Single-cell full-length total RNA sequencing uncovers dynamics of recursive splicing and enhancer RNAs., 2018, 9(1): 619.

    [31] Pai AA, Henriques Telmo, Paggi Joseph, Burkholder Adam, Adelman Karen, Burge Christopher B. Intron length and recursive sites are major determinants of splicing efficiency in flies., 2017, February 13.

    [32] Srndic, Ismet. Experimental verification and functional exploration of intrasplicing., 2016, urn:nbn:at:at-ubw: 1-10109.80624.509568-4.

    [33] Irimia M, Blencowe BJ. Alternative splicing: decoding an expansive regulatory layer., 2012, 24(3): 323–332.

    [34] Carrillo OF, Bieberstein N, Neugebauer KM. Pause locally, splice globally., 2011, 21(6): 328–335.

    [35] Maurano MT, Humbert R, Rynes E, Thurman RE, Haugen E, Wang H, Reynolds AP, Sandstrom R, Qu H, Brody J, Shafer A, Neri F, Lee K, Kutyavin T, Stehling-Sun S, Johnson AK, Canfield TK, Giste E, Diegel M, Bates D, Hansen RS, Neph S, Sabo PJ, Heimfeld S, Raubitschek A, Ziegler S, Cotsapas C, Sotoodehnia N, Glass I, Sunyaev SR, Kaul R, Stamatoyannopoulos JA. Systematic localization of common disease-associated variation in regulatory DNA., 2012, 337(6099): 1190–1195.

    [36] Luo Q, Gao C, Wang HL, Zhou JH, Gao TZ. Effect of a novel splicing mutation(IVS2-2A→C) of SEDL gene on RNA processing., 2005, 27(4): 544– 548.羅強, 高超, 王懷立, 周建華, 高鐵錚. X-連鎖遲發(fā)性脊椎骨骺發(fā)育不良基因剪接受體突變對mRNA加工的影響. 遺傳, 2005, 27(4): 544–548.

    [37] Savisaar R, Hurst LD. Estimating the prevalence of functional exonic splice regulatory information., 2017, 136(9): 1059–1078.

    [38] Editor T. Novel mechanism, called recursive splicing reveals how cells process large genes.,2005.

    [39] Zou YX, Gong YQ. Aberrant RNA splicing as the molecular basis of some pathogenic variants., 2017, 39(3): 200–207.鄒永新, 龔瑤琴. 影響RNA剪接的基因變異. 遺傳, 2017, 39(3): 200–207.

    [40] Burke JE, Longhurst AD, Merkurjev D, Sales-Lee J, Rao B, Moresco J, Yates JR 3rd, Li JJ, Madhani HD. Spliceosome profiling visualizes operations of a dynamic RNP at nucleotide resolution., 2018, 173(4): 1014– 1030.e17.

    [41] Donovan PD, Schr?der MS, Higgins DG, Butler G. Identification of non-coding RNAs in the candida parapsilosis species group., 2016, 11(9): e0163235.

    [42] Chabot B, Shkreta L. Defective control of pre-messenger RNA splicing in human disease., 2016, 212(1): 13–27.

    [43] Georgomanolis T, Sofiadis K, Papantonis A. Cutting a long intron short: recursive splicing and its implications., 2016, 7: 598.

    [44] Chen YM. Regulation and function of non-exonic recursive splicing. 2011.

    [45] Hollander D, Naftelberg S, Lev-Maor G, Kornblihtt AR, Ast G. How are short exons flanked by long introns defined and committed to splicing?, 2016, 32(10): 596–606.

    [46] Liu YH, Wang L, Yu L. The principle and application of the single-molecule real-time sequencing technology., 2015, 37(3): 259–268.柳延虎, 王璐, 于黎. 單分子實時測序技術的原理與應用. 遺傳, 2015, 37(3): 259–268.

    Molecular mechanisms of recursive splicing events in long introns of eukaryotes

    Jinchuan Wei1, Tianyi Xu1, Jing Wu1,2, Xiaofeng Song1

    Recursive splicing refers to the biological process that long introns are removed in multiple steps during pre-mRNA splicing. In comparison to large introns (>10 kb), most introns in higher eukaryotic genomes are removed in one step during transcription. Previous studies have revealed that recursive splicing events play important roles in many biological processes, including the pathogenesis and development of diseases. In recent years, more researchers have focused on recursive splicing events and found that recursive splicing occurs inand many other vertebrates. Multiple recursive splicing sites have been predicted by different bioinformatics methods and verified by experiments. Current researches focus on the process of recursive splicing, recursive splicing site recognition and its influence on biological processes. In this review, we summarize the molecular mechanism of recursive splicing events in eukaryotic genomes and the present development in this field, aiming tolay the basis for further understanding of the mechanisms of RNA splicing.

    recursive splicing; lariat structure; RS-sites; Total RNA-seq

    2018-10-30;

    2018-12-24

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

    魏金川,碩士研究生,專業(yè)方向:生物信息學。E-mail: weijc@nuaa.edu.cn

    宋曉峰,博士,教授,博士生導師,研究方向:生物信息學。E-mail: xfsong@nuaa.edu.cn

    吳靜,碩士,副教授,研究方向:生物信息學。E-mail: wujing@njmu.edu.cn

    10.16288/j.yczz.18-182

    2019/1/16 9:14:00

    URI: http://kns.cnki.net/kcms/detail/11.1913.R.20190116.0914.001.html

    (責任編委: 張根發(fā))

    猜你喜歡
    套索內含子堿基
    線粒體核糖體蛋白基因中內含子序列間匹配特性分析
    應用思維進階構建模型 例談培養(yǎng)學生創(chuàng)造性思維
    套索肽替代抗生素應用的研究進展
    飼料博覽(2020年10期)2020-11-20 09:52:48
    中國科學家創(chuàng)建出新型糖基化酶堿基編輯器
    套索肽
    ——一種可供多肽類藥物設計和改造的多功能載體
    藥學進展(2019年10期)2019-12-23 04:03:16
    不同方向內含子對重組CHO細胞中神經生長因子表達的影響
    數(shù)字化圖像的簡單合成教學分析
    新一代(2019年7期)2019-06-25 10:36:04
    生命“字母表”迎來4名新成員
    科學24小時(2019年5期)2019-06-11 08:39:38
    生命“字母表”迎來4名新成員
    更 正
    生物學通報(2019年1期)2019-02-15 16:33:43
    亚洲美女黄片视频| 日本撒尿小便嘘嘘汇集6| 午夜福利在线在线| 国产又黄又爽又无遮挡在线| 我要搜黄色片| 精华霜和精华液先用哪个| 99久久久亚洲精品蜜臀av| 欧美+日韩+精品| 91麻豆av在线| 亚洲成人久久爱视频| 91在线精品国自产拍蜜月 | 99国产综合亚洲精品| 深夜精品福利| 欧美一区二区国产精品久久精品| 亚洲美女视频黄频| 国产精品一及| 最新中文字幕久久久久| 日韩成人在线观看一区二区三区| 日韩精品青青久久久久久| 免费在线观看影片大全网站| 国产一区二区三区视频了| 白带黄色成豆腐渣| 国产免费男女视频| 日韩欧美在线乱码| 欧美乱码精品一区二区三区| 日韩亚洲欧美综合| 丁香欧美五月| 国产视频一区二区在线看| 国产激情欧美一区二区| 最后的刺客免费高清国语| 亚洲激情在线av| 日韩精品青青久久久久久| 国产精品一区二区三区四区免费观看 | 亚洲精品影视一区二区三区av| 亚洲五月天丁香| 欧美激情在线99| 色播亚洲综合网| 精品久久久久久久末码| 日本a在线网址| 女生性感内裤真人,穿戴方法视频| 男人舔女人下体高潮全视频| 欧美三级亚洲精品| 久久久久免费精品人妻一区二区| 色精品久久人妻99蜜桃| 亚洲av免费在线观看| 精品电影一区二区在线| 欧美绝顶高潮抽搐喷水| 中文在线观看免费www的网站| 国产97色在线日韩免费| 亚洲最大成人手机在线| а√天堂www在线а√下载| 噜噜噜噜噜久久久久久91| av欧美777| 大型黄色视频在线免费观看| 欧美日韩综合久久久久久 | 国产乱人伦免费视频| 国产在视频线在精品| 亚洲男人的天堂狠狠| 成人亚洲精品av一区二区| 国产成人啪精品午夜网站| 日本三级黄在线观看| 国产成人av激情在线播放| 国产亚洲精品久久久com| 久久亚洲真实| 亚洲av电影在线进入| 亚洲av美国av| 男女下面进入的视频免费午夜| 亚洲自拍偷在线| 国产精品av视频在线免费观看| 国产91精品成人一区二区三区| 露出奶头的视频| 变态另类丝袜制服| 免费高清视频大片| 亚洲精品美女久久久久99蜜臀| 国产高清视频在线观看网站| 亚洲精品国产精品久久久不卡| 国产精品自产拍在线观看55亚洲| 亚洲欧美一区二区三区黑人| 两个人视频免费观看高清| 99国产精品一区二区三区| 特级一级黄色大片| 岛国在线观看网站| 变态另类成人亚洲欧美熟女| 91麻豆av在线| 欧美成人a在线观看| 99在线人妻在线中文字幕| 亚洲精品日韩av片在线观看 | 久久久久国产精品人妻aⅴ院| 久久久久精品国产欧美久久久| 日本黄色片子视频| 首页视频小说图片口味搜索| svipshipincom国产片| 精品久久久久久久久久免费视频| 色综合婷婷激情| 欧美日韩瑟瑟在线播放| 日韩中文字幕欧美一区二区| 国产精品综合久久久久久久免费| 99久久久亚洲精品蜜臀av| 亚洲精品在线观看二区| 国产真实伦视频高清在线观看 | 深爱激情五月婷婷| 老汉色av国产亚洲站长工具| 久久中文看片网| 综合色av麻豆| 一本一本综合久久| 日本撒尿小便嘘嘘汇集6| 啦啦啦韩国在线观看视频| 国产av麻豆久久久久久久| 国产成人av激情在线播放| 久久草成人影院| 国产成人影院久久av| 日日干狠狠操夜夜爽| 嫩草影院精品99| 国产黄a三级三级三级人| 精品乱码久久久久久99久播| 久久精品国产亚洲av香蕉五月| 少妇丰满av| 中文字幕人成人乱码亚洲影| 色在线成人网| 好男人电影高清在线观看| 精品国产超薄肉色丝袜足j| 亚洲中文字幕日韩| 国产一区二区三区视频了| 国产乱人视频| 国产高清视频在线播放一区| 亚洲成人久久性| 午夜福利高清视频| 日韩欧美国产一区二区入口| 亚洲美女视频黄频| 亚洲精品国产精品久久久不卡| а√天堂www在线а√下载| 一边摸一边抽搐一进一小说| 国产中年淑女户外野战色| 老司机福利观看| 亚洲欧美日韩高清在线视频| 国产精品 欧美亚洲| xxxwww97欧美| 亚洲成a人片在线一区二区| svipshipincom国产片| 国产精品野战在线观看| 网址你懂的国产日韩在线| 一区二区三区激情视频| 亚洲 欧美 日韩 在线 免费| 久久精品人妻少妇| 在线免费观看的www视频| 狠狠狠狠99中文字幕| 亚洲五月婷婷丁香| 国产v大片淫在线免费观看| av专区在线播放| 两个人的视频大全免费| 日韩成人在线观看一区二区三区| 精品乱码久久久久久99久播| 91久久精品国产一区二区成人 | 免费在线观看亚洲国产| 国产aⅴ精品一区二区三区波| 少妇的逼水好多| 国产爱豆传媒在线观看| 色综合亚洲欧美另类图片| 亚洲avbb在线观看| 亚洲欧美日韩高清在线视频| 成年女人看的毛片在线观看| 一进一出好大好爽视频| 天天躁日日操中文字幕| tocl精华| 国产中年淑女户外野战色| 最近最新中文字幕大全电影3| 国产成+人综合+亚洲专区| 一进一出抽搐gif免费好疼| 男女床上黄色一级片免费看| 国产伦人伦偷精品视频| 成人av在线播放网站| 国产日本99.免费观看| 色播亚洲综合网| 欧美性猛交╳xxx乱大交人| 国产主播在线观看一区二区| 日本免费一区二区三区高清不卡| 一个人看的www免费观看视频| 一个人看的www免费观看视频| 人人妻人人看人人澡| 精品国产亚洲在线| 国产三级在线视频| 他把我摸到了高潮在线观看| 亚洲av电影在线进入| 色精品久久人妻99蜜桃| 欧美成人免费av一区二区三区| av女优亚洲男人天堂| 99久久精品国产亚洲精品| 午夜精品一区二区三区免费看| 小蜜桃在线观看免费完整版高清| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 午夜精品一区二区三区免费看| 午夜免费男女啪啪视频观看 | 亚洲精品乱码久久久v下载方式 | 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av中文字字幕乱码综合| 麻豆一二三区av精品| 日本五十路高清| 欧美一区二区精品小视频在线| 此物有八面人人有两片| 亚洲欧美日韩无卡精品| 精品久久久久久久毛片微露脸| 亚洲五月天丁香| 亚洲avbb在线观看| 国产精品一区二区三区四区免费观看 | 欧美在线黄色| 久久久国产精品麻豆| 色播亚洲综合网| 在线看三级毛片| 18禁黄网站禁片午夜丰满| 91麻豆av在线| 亚洲国产欧美网| 久久香蕉精品热| 国产真实伦视频高清在线观看 | 成人一区二区视频在线观看| 三级男女做爰猛烈吃奶摸视频| 国产成人系列免费观看| 黄色视频,在线免费观看| 最近在线观看免费完整版| av黄色大香蕉| 亚洲精品一卡2卡三卡4卡5卡| 国内少妇人妻偷人精品xxx网站| 国产一区在线观看成人免费| 色综合亚洲欧美另类图片| 长腿黑丝高跟| 日日干狠狠操夜夜爽| 午夜激情欧美在线| 成人一区二区视频在线观看| 九九在线视频观看精品| 欧美绝顶高潮抽搐喷水| 看免费av毛片| 成人精品一区二区免费| 国产伦一二天堂av在线观看| 成人欧美大片| 好看av亚洲va欧美ⅴa在| 成人高潮视频无遮挡免费网站| 波野结衣二区三区在线 | 好看av亚洲va欧美ⅴa在| 精品久久久久久久久久久久久| 尤物成人国产欧美一区二区三区| 亚洲avbb在线观看| 亚洲av电影不卡..在线观看| 久久草成人影院| 非洲黑人性xxxx精品又粗又长| 俺也久久电影网| 亚洲五月天丁香| 日本 欧美在线| 日本免费a在线| 叶爱在线成人免费视频播放| 极品教师在线免费播放| 国产精品精品国产色婷婷| 乱人视频在线观看| 在线观看免费午夜福利视频| 日日干狠狠操夜夜爽| 免费搜索国产男女视频| 69av精品久久久久久| 天美传媒精品一区二区| 精品一区二区三区视频在线 | 成人高潮视频无遮挡免费网站| 搡老熟女国产l中国老女人| 手机成人av网站| 亚洲乱码一区二区免费版| 一区福利在线观看| 国产精品爽爽va在线观看网站| 国产美女午夜福利| 啦啦啦观看免费观看视频高清| 999久久久精品免费观看国产| 丁香欧美五月| 国产精品女同一区二区软件 | 国产极品精品免费视频能看的| 久久婷婷人人爽人人干人人爱| 我要搜黄色片| 黄色片一级片一级黄色片| 免费看日本二区| 99久久精品一区二区三区| 日本免费一区二区三区高清不卡| 美女黄网站色视频| 日本免费a在线| 别揉我奶头~嗯~啊~动态视频| 女人被狂操c到高潮| 夜夜躁狠狠躁天天躁| 在线观看舔阴道视频| 亚洲专区国产一区二区| 日韩精品青青久久久久久| 国产蜜桃级精品一区二区三区| 亚洲无线在线观看| 性色avwww在线观看| 欧美黄色淫秽网站| 免费观看精品视频网站| 免费无遮挡裸体视频| 久久精品人妻少妇| 欧美黑人欧美精品刺激| 成人午夜高清在线视频| 麻豆成人午夜福利视频| 国产精品永久免费网站| 看免费av毛片| 欧美三级亚洲精品| 亚洲精品色激情综合| 99国产精品一区二区蜜桃av| 国产中年淑女户外野战色| 在线播放国产精品三级| 黄色日韩在线| 国产私拍福利视频在线观看| 99riav亚洲国产免费| 观看免费一级毛片| 色哟哟哟哟哟哟| 99久久精品热视频| 亚洲狠狠婷婷综合久久图片| 一进一出抽搐动态| 一区二区三区国产精品乱码| 国产乱人视频| 国产一区二区三区在线臀色熟女| 热99re8久久精品国产| 精品久久久久久久人妻蜜臀av| 亚洲一区二区三区色噜噜| 国产成人系列免费观看| 长腿黑丝高跟| 色播亚洲综合网| 成熟少妇高潮喷水视频| 老熟妇乱子伦视频在线观看| 色噜噜av男人的天堂激情| 18禁美女被吸乳视频| 亚洲人成伊人成综合网2020| 久久人妻av系列| 精品一区二区三区视频在线观看免费| www.999成人在线观看| 久久性视频一级片| 99精品在免费线老司机午夜| 一进一出抽搐gif免费好疼| 国产一区二区在线av高清观看| 18禁黄网站禁片午夜丰满| 欧美日本视频| 好男人电影高清在线观看| 国产精品三级大全| 狂野欧美白嫩少妇大欣赏| 看黄色毛片网站| 久久久久亚洲av毛片大全| 欧美性感艳星| 久久久国产成人精品二区| 国产熟女xx| 最新中文字幕久久久久| 国产精品精品国产色婷婷| 夜夜爽天天搞| 亚洲欧美精品综合久久99| 国产高清三级在线| 小说图片视频综合网站| 日韩免费av在线播放| 免费高清视频大片| 国产黄色小视频在线观看| 内地一区二区视频在线| 国产激情欧美一区二区| 全区人妻精品视频| 久久香蕉精品热| 国产一区在线观看成人免费| 老鸭窝网址在线观看| 最近最新中文字幕大全电影3| 国产一级毛片七仙女欲春2| 精品乱码久久久久久99久播| 国产淫片久久久久久久久 | 久久久久久大精品| 热99re8久久精品国产| 亚洲成av人片在线播放无| 91av网一区二区| 天天躁日日操中文字幕| 精品久久久久久久久久久久久| 亚洲欧美日韩东京热| 亚洲国产精品sss在线观看| 69人妻影院| 成人亚洲精品av一区二区| 在线十欧美十亚洲十日本专区| 无遮挡黄片免费观看| 蜜桃久久精品国产亚洲av| 欧美绝顶高潮抽搐喷水| 日韩免费av在线播放| 99精品欧美一区二区三区四区| www.www免费av| 国产综合懂色| 亚洲人成网站在线播放欧美日韩| 久久99热这里只有精品18| 欧美成人a在线观看| 成人18禁在线播放| 精品不卡国产一区二区三区| 国产精品一及| 美女cb高潮喷水在线观看| 亚洲精品乱码久久久v下载方式 | 国产高清视频在线观看网站| 高清毛片免费观看视频网站| 99久久成人亚洲精品观看| 人人妻人人看人人澡| 在线播放无遮挡| 波野结衣二区三区在线 | 国产精品亚洲一级av第二区| 亚洲欧美日韩无卡精品| 国产黄色小视频在线观看| 成年免费大片在线观看| 搡老熟女国产l中国老女人| 2021天堂中文幕一二区在线观| 国产av不卡久久| 成人性生交大片免费视频hd| 国产精品电影一区二区三区| 日本黄大片高清| 国产高清激情床上av| a级毛片a级免费在线| 亚洲av免费高清在线观看| 欧美一级毛片孕妇| 亚洲中文日韩欧美视频| 精品一区二区三区人妻视频| 精品国产三级普通话版| 日本a在线网址| 无人区码免费观看不卡| 国产综合懂色| 在线观看66精品国产| 88av欧美| 97超级碰碰碰精品色视频在线观看| 亚洲无线观看免费| 草草在线视频免费看| 成人无遮挡网站| 丁香六月欧美| 久久亚洲精品不卡| 无人区码免费观看不卡| 在线观看免费视频日本深夜| 19禁男女啪啪无遮挡网站| 首页视频小说图片口味搜索| 免费无遮挡裸体视频| 亚洲熟妇中文字幕五十中出| 亚洲人成伊人成综合网2020| 麻豆久久精品国产亚洲av| 亚洲av熟女| 真人做人爱边吃奶动态| 免费av毛片视频| 桃色一区二区三区在线观看| 久久香蕉精品热| 九九热线精品视视频播放| 精品久久久久久成人av| av在线蜜桃| 亚洲精华国产精华精| 在线播放无遮挡| 99热只有精品国产| 99久久综合精品五月天人人| 身体一侧抽搐| 日韩欧美免费精品| 最近最新中文字幕大全电影3| 免费看a级黄色片| 亚洲五月婷婷丁香| 亚洲一区二区三区不卡视频| 日韩欧美三级三区| 老司机午夜福利在线观看视频| 日韩中文字幕欧美一区二区| 99国产综合亚洲精品| 脱女人内裤的视频| 成人av一区二区三区在线看| 亚洲精品国产精品久久久不卡| 三级毛片av免费| 亚洲成av人片在线播放无| 18禁黄网站禁片免费观看直播| 欧美国产日韩亚洲一区| 五月玫瑰六月丁香| 国产亚洲精品久久久com| 精品久久久久久成人av| 亚洲国产精品合色在线| 美女大奶头视频| 男女之事视频高清在线观看| 亚洲av不卡在线观看| 九九在线视频观看精品| 有码 亚洲区| 夜夜夜夜夜久久久久| 日日摸夜夜添夜夜添小说| 最好的美女福利视频网| 中亚洲国语对白在线视频| 久久人妻av系列| 国产视频一区二区在线看| 男女午夜视频在线观看| 免费电影在线观看免费观看| 青草久久国产| 国产色爽女视频免费观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲av电影在线进入| www国产在线视频色| 一区二区三区激情视频| 天堂√8在线中文| 久久草成人影院| 夜夜看夜夜爽夜夜摸| 琪琪午夜伦伦电影理论片6080| 女人被狂操c到高潮| 99国产精品一区二区三区| 国产真实伦视频高清在线观看 | 欧美一级a爱片免费观看看| 丁香欧美五月| 91字幕亚洲| 麻豆国产av国片精品| 又紧又爽又黄一区二区| 欧美zozozo另类| 一级毛片高清免费大全| 五月玫瑰六月丁香| 黑人欧美特级aaaaaa片| 最好的美女福利视频网| 黄片大片在线免费观看| 久久精品国产亚洲av涩爱 | 一卡2卡三卡四卡精品乱码亚洲| 久久久久久人人人人人| 色精品久久人妻99蜜桃| 香蕉久久夜色| 欧美黄色片欧美黄色片| 久久草成人影院| 欧美三级亚洲精品| 欧美成人一区二区免费高清观看| 成人永久免费在线观看视频| 床上黄色一级片| 日本在线视频免费播放| 欧美日本亚洲视频在线播放| 黄片大片在线免费观看| 久久精品国产亚洲av涩爱 | 男女那种视频在线观看| 搡老岳熟女国产| 欧美色欧美亚洲另类二区| 91久久精品电影网| 午夜老司机福利剧场| 免费看光身美女| 叶爱在线成人免费视频播放| 天堂av国产一区二区熟女人妻| 日韩欧美三级三区| 亚洲国产精品久久男人天堂| 亚洲欧美日韩无卡精品| 99久久精品国产亚洲精品| 国产精品综合久久久久久久免费| 亚洲精华国产精华精| 国产免费男女视频| 亚洲内射少妇av| 亚洲第一电影网av| 最好的美女福利视频网| 亚洲人与动物交配视频| 三级国产精品欧美在线观看| 亚洲精品久久国产高清桃花| 国产精品乱码一区二三区的特点| 国产一区二区三区在线臀色熟女| 日本三级黄在线观看| 欧美av亚洲av综合av国产av| 国产日本99.免费观看| 在线观看av片永久免费下载| 青草久久国产| 九九在线视频观看精品| 亚洲五月天丁香| 亚洲av五月六月丁香网| 免费看美女性在线毛片视频| 午夜免费男女啪啪视频观看 | 精品人妻偷拍中文字幕| 精品欧美国产一区二区三| 午夜两性在线视频| 母亲3免费完整高清在线观看| 91九色精品人成在线观看| 国产亚洲欧美在线一区二区| 夜夜躁狠狠躁天天躁| 亚洲av免费在线观看| 久久人妻av系列| 99久久精品热视频| 欧洲精品卡2卡3卡4卡5卡区| 91av网一区二区| 老司机午夜福利在线观看视频| 免费看a级黄色片| 国产麻豆成人av免费视频| 午夜福利高清视频| 国产欧美日韩一区二区精品| 亚洲一区二区三区色噜噜| 两个人视频免费观看高清| 午夜视频国产福利| 国产一区二区亚洲精品在线观看| www.熟女人妻精品国产| 日本与韩国留学比较| 日本撒尿小便嘘嘘汇集6| 午夜福利在线观看免费完整高清在 | 757午夜福利合集在线观看| 国产亚洲欧美在线一区二区| 国产激情欧美一区二区| АⅤ资源中文在线天堂| 国内精品久久久久精免费| 久久精品国产自在天天线| 特大巨黑吊av在线直播| 在线视频色国产色| 久久精品国产99精品国产亚洲性色| 久久人人精品亚洲av| 日本黄色片子视频| 亚洲美女黄片视频| 夜夜爽天天搞| av福利片在线观看| 欧美成人一区二区免费高清观看| 亚洲专区国产一区二区| 欧美xxxx黑人xx丫x性爽| 又紧又爽又黄一区二区| 18美女黄网站色大片免费观看| 欧美日韩精品网址| 嫩草影院入口| 黄色成人免费大全| 中国美女看黄片| 日本黄大片高清| 操出白浆在线播放| 好男人电影高清在线观看| 99久久精品国产亚洲精品| 亚洲国产精品999在线| 国产久久久一区二区三区| 久久中文看片网| 真实男女啪啪啪动态图| 亚洲精品色激情综合| 脱女人内裤的视频| 九色成人免费人妻av| 99精品在免费线老司机午夜| 国产精品免费一区二区三区在线| 黄色丝袜av网址大全| 久久中文看片网| 桃色一区二区三区在线观看| 12—13女人毛片做爰片一| 亚洲国产色片| 嫩草影视91久久| 狂野欧美白嫩少妇大欣赏| 男人舔奶头视频| 国产69精品久久久久777片| 午夜福利18|