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

    疾病相關miRNAs識別方法研究進展

    2019-05-07 01:33:264
    廣州大學學報(自然科學版) 2019年1期
    關鍵詞:課題組靶向編碼

    4

    (1.廣州大學 計算科技研究院, 廣東 廣州 510006; 2.溫州職業(yè)技術學院 信息技術系, 浙江 溫州 325035;3.廣州大學 人事處, 廣東廣州 510006; 4.廣東省數(shù)學教育軟件工程技術研究中心, 廣東 廣州 510006)

    0 引 言

    非編碼RNA是一類不能編碼蛋白質(zhì)的RNA分子總稱,約占人類基因組總長度的97%,在過去很長一段時間內(nèi)被稱為“垃圾序列”或者“暗物質(zhì)”.根據(jù)序列長度的大小可以將非編碼RNA分為3類:<50 nt, 包括microRNA,siRNA,piRNA等等;50~500 nt,包括lncRNA,rRNA,snRNA,snoRNA等等;大于500 nt,包括長的mRNA-like 的非編碼RNA,長的不帶polyA 尾巴的非編碼RNA等等.隨著研究的深入,人們發(fā)現(xiàn)雖然非編碼RNA不能直接編碼蛋白質(zhì),但是其中一些對于維持細胞內(nèi)部的平衡,以及與疾病的產(chǎn)生和惡化有著密切的關系.

    miRNAs(microRNAs)是一類重要的非編碼RNAs,是基因表達和蛋白翻譯過程中的調(diào)節(jié)因子.已有研究表明,miRNA廣泛參與細胞的增殖、凋亡及分化生物過程,在腫瘤的發(fā)生、形成過程中扮演著十分重要的角色.腫瘤的發(fā)生、發(fā)展是多條通路聯(lián)合破壞的結果,因此,腫瘤細胞中的“基因網(wǎng)絡特征”比“單分子特征”更能揭示疾病惡化和進展的機制[1-3].由于一個miRNA能夠調(diào)控多個基因的表達,影響多個信號通路的活性,因此,Chen[4]在著名的新英格蘭醫(yī)學期刊上預測認為,將miRNA當作腫瘤生物治療的靶分子或者藥物將比編碼基因更加有效.

    權威數(shù)據(jù)庫miRBase收錄的miRNAs數(shù)目從創(chuàng)辦之初的218條,到目前已經(jīng)增加到38 589條,數(shù)量增加了177倍(圖1),這一定程度上反映了miRNAs研究領域發(fā)展之迅速.此外,筆者通過Web of Science檢索關鍵詞“miRNA OR microRNA”梳理了從2001 年至今的發(fā)表的相關文獻,結果表明,miRNA 相關的研究文章從2001 年8篇到2018年增加至15 587篇,數(shù)量增長倍數(shù)高達1 900 多倍.最近幾年有關miRNA的研究文章,仍處于逐年上升的趨勢,這說明miRNA領域的研究依然是當前的熱點課題.

    目前miRNA研究領域主要集中在以下4個方面:①靶基因的預測.通過開發(fā)新的計算方法預測miRNA與基因之間的靶向關系;②miRNA生物學功能作用機制.通過生物學實驗的方式研究miRNA在生物體內(nèi)的功能作用;③疾病相關miRNAs的篩選方法.以miRNA表達譜數(shù)據(jù)為基礎,計算篩選與疾病密切相關的miRNAs;④miRNA小分子藥物的開發(fā)及應用.以miRNA為藥物或者藥物靶點,開發(fā)適用于治療特定疾病的新型藥物.本文擬從以上幾方面闡述疾病相關miRNAs的篩選方法研究進展.

    圖1 miRNAs相關研究及已發(fā)現(xiàn)的miRNAs數(shù)量統(tǒng)計數(shù)據(jù)Fig.1 miRNAs related studies and the statistic data for miRNAs 數(shù)據(jù)檢索時間2019年1月6日

    1 miRNA基本概念

    1.1 miRNA定義

    miRNAs是一類重要的內(nèi)源性非編碼微小RNAs分子,其長度大小約為18~25核苷酸(Nucleotide,nt).在細胞核內(nèi),miRNA基因的初始轉(zhuǎn)錄產(chǎn)物(pri-miRNA)在RNA聚合酶Ⅲ(RNase Ⅲ)Drosha的作用下切割成為前體miRNA即pre-miRNA[5-7].然后pre-miRNA在轉(zhuǎn)運蛋白(exportin-5)的作用下,從細胞核內(nèi)轉(zhuǎn)移到細胞質(zhì)中[8].最后,在另一種RNA聚合酶Ⅲ(RNase Ⅲ)Dicer的作用下,切割生成成熟的單鏈miRNAs[9].

    miRNA在生物體內(nèi)發(fā)揮的功能作用是相當復雜的,總結起來有如下特性:①結合位點差異性.成熟的miRNAs能夠通過堿基互補配對的方式結合到靶基因mRNA 3’ 端非轉(zhuǎn)錄區(qū)域,當完全匹配時候,會促進靶基因的降解,而當非完全匹配時候,則會抑制其靶基因翻譯生成對應的蛋白質(zhì)[10-11];②時空動態(tài)特性.當細胞周期阻滯的時候,miRNA能夠引導AGO等蛋白與AREs(AU-rich elemens)相結合促進基因的翻譯,而在增殖細胞中miRNA對靶基因又起到抑制的作用[12];③位置的影響.除了在細胞質(zhì)中發(fā)揮功能,Xiao等[13]發(fā)現(xiàn)miR-24-1在細胞核內(nèi)通過與基因的增強子區(qū)域結合,從而促進基因的表達.

    1.2 miRNA命名規(guī)則

    在研究miRNA之初,miRNA的名稱主要是根據(jù)其表型進行命名,如lin-4、let-7等等.隨著越來越多miRNA的出現(xiàn),為了方便注釋和研究,人們提出了一些關于miRNA的命名標準,目前最新版本的miRBase數(shù)據(jù)庫中關于miRNA的命名準則見圖2.

    圖2 miRNA命名規(guī)則簡要示意圖Fig.2 A brief sketch of the naming rules for miRNA.

    miRNA成熟體命名規(guī)則如下(以動物miRNA為例):

    (1)在統(tǒng)一命名規(guī)則制定之前發(fā)現(xiàn)的miRNA,則繼續(xù)保留原來名字,如hsa-lin-4;

    (2)miRNA的成熟體表示為miR,然后根據(jù)其物種類型,以及被發(fā)現(xiàn)的先后順序標注阿拉伯數(shù)字,如hsa-miR-21;

    (3)若多個miRNAs高度同源,則在數(shù)字后面加上英文小寫字母(a,b,c,…),如hsa-miR-34a,hsa-miR-34b等;

    (4)由不同染色體上的DNA序列轉(zhuǎn)錄加工而成的具有相同成熟體序列的miRNAs,就在后面添加上阿拉伯數(shù)字,如hsa-miR-199a-1和hsa-miR-199a-2;

    (5)通常一個miRNA前體長度大約為70~80nt,有時候2個臂分別會產(chǎn)生成熟體miRNAs.在之前人們的做法是:表達量比較高的miRNA后面不添加任何符號,而對于表達量比較低的miRNA后面則加上*號,如rno-miR-9*.而最新版本的miRBase數(shù)據(jù)庫中則以“-5p”和“-3p”對成熟miRNAs分別命名.如hsa-miR-21-5p和hsa-miR-21-3p,分別表示從前體has-mir-21的5’端臂和3’端臂產(chǎn)生而來.

    2 miRNA靶基因預測方法

    已有研究表明,人類相關基因中大約有三分之一的基因會受到miRNA的靶向調(diào)控.到目前為止,雖然陸續(xù)有生物學實驗證實了一些miRNAs與mRNAs之間的靶向調(diào)節(jié)關系,但是仍然還有許多靶基因沒有經(jīng)過生物學實驗驗證,這一定程度上阻礙了miRNAs生物學的功能被充分的發(fā)掘[10, 14-15].因此,關于miRNAs一項重要的研究內(nèi)容同時也是極具挑戰(zhàn)性的工作,就是如何快速、準確地預測以及驗證miRNAs的靶基因.該項研究對于探究miRNAs生物學功能作用、參與涉及的生物信號通路以及它在疾病的發(fā)展過程中扮演的角色具有重要的意義,與此同時,也能夠加快miRNAs作為一種治療靶點或者藥物走向臨床實踐.

    借助計算生物學的專業(yè)優(yōu)勢,將其用來替代生物學家大量的分子克隆實驗工作,基于這種研究模式獲得了非常好的效果[16-18].但是,預測動物的miRNAs靶基因是一項極具挑戰(zhàn)性的工作,其主要原因是由于許多動物的miRNAs與其靶基因mRNAs結合位點并非完全互補的.除此以外,由于已知的經(jīng)過實驗驗證的靶向關系比較少,所以在評價預測結果的時候,沒有足夠的金標準來作為客觀的參考對象.到目前為止,對于計算機方法預測、推斷出來的靶基因,還沒有一個快速、可靠的高通量驗證鑒定方案,這在一定程度上影響了算法準確性的評價,也阻礙了計算預測算法的更進一步的優(yōu)化及改進[19].從成為生命科學界的研究熱點到目前為止,miRNAs相關生物學功能的探索發(fā)現(xiàn)研究工作已經(jīng)歷了10多年的發(fā)展歷程,關于miRNAs靶基因的預測方法已經(jīng)發(fā)展出數(shù)十種之多.這些靶基因預測計算方法主要可以劃分為2大類:ab initio計算方法和機器學習預測方法.總體上來說,剛開始提出的計算預測方法均屬于ab initio計算方法,這些預測方法都是基于實驗獲得的結構特征來引導預測方法的開發(fā).機器學習預測方法主要是基于實驗訓練集的相似度來預測潛在的靶向目標.

    2.1 ab initio方法

    (1)miRanda[20-21]:MiRanda以及相應的網(wǎng)站microRNA.org是由美國著名的癌癥研究中心——斯隆-凱特琳研究所(Sloan-Kettering Institute,SKI)開發(fā)的在線查詢數(shù)據(jù)庫,該數(shù)據(jù)資源可以非常方便地免費在線獲取.miRanda的最新版本發(fā)布時間是在2010年8月,該最新版本也稱為mirSVR.該方法適用范圍廣,不受物種限制.MiRanda采用Smith-Whatman算法的思想[22],利用該方法來檢測互補匹配的序列(允許位于miRNA的5p端第2-7nt種子序列有一個G∶U匹配或者錯誤匹配的存在),除此以外,這一算法還增加了關于RNA二級結構預測的程序,從而可以從動力學的角度去考察miRNA與其候選靶基因之間的結合穩(wěn)定性;根據(jù)前面的準則得到了靶基因結合位點,然后考察這些結合位點在不同物種間的保守性,最后給出自由能(△G)以及序列匹配分數(shù)(Score).miRanda靶基因預測方法主要考慮了以下幾個方面因素:①靶向結合序列的互補性;②mRNA序列的二級結構對miRNA-mRNA結合穩(wěn)定性的影響;③靶基因結合位點在跨物種間的保守性.

    (2)TargetScan[23-24]:是由Lewis等于2003年開發(fā)出來的miRNA靶基因預測工具,該方法首次引入假陽性率來評估預測結果.TargetScan靶基因預測算法主要考慮如下2個方面特點:①考察miRNA靶向基因mRNA的3p端非轉(zhuǎn)錄區(qū)序列在跨基因跨物種間的保守性;②研究miRNA與靶基因mRNA形成的miRNA-mRNA雙鏈結構的動力學穩(wěn)定性[25].該方法要求miRNA的5p端的第2-8位堿基(種子序列)與靶基因mRNA的3p端必須滿足互補匹配,然后要求種子序列完全互補的情況下從種子序列向兩端擴展,直到遇到不能配對的堿基為止,這個過程中允許G-U配對.該方法通過計算靶基因結合位點產(chǎn)生的自由能以及信噪比來評價計算機預測結果的準確性.

    (3)PicTar[26]:該算法同樣要求miRNA 5P端1-7nt或者2-8nt種子序列在靶基因位點識別中的關鍵因素,主要強調(diào)miRNA與其靶基因形成的二聚體結合能在靶基因翻譯抑制過程中的影響.PicTar將種子序列分為2類:“完全匹配的種子序列”和“不完全匹配的種子序列”,前者需要種子序列與靶基因序列完全匹配,后者則可以在滿足miRNA靶基因二聚體結合能要求的情況下允許種子序列出現(xiàn)錯配(不允許G-U配對).

    (4)DIANA-microT:基于實驗和計算預測方法,Kiriakidou等于2004年開發(fā)出一種“DIANA-microT”miRNA靶基因預測程序[27].該方法結合了實驗與計算機計算方法來預測哺乳動物的靶基因,這個方法的一大特點是它著重考慮了單一結合位點(miRNA∶mRNA)的miRNA靶基因.在篩選miRNA靶基因的過程中,DIANA-microT不僅考慮了重要的種子序列區(qū)域,而且考慮了miRNA的3p端與靶基因mRNA的結合.在識別靶基因時候,不僅考慮保守位點,而且也考察了非保守結合位點.

    (5)PITA[28]:該算法由Kertesz等于2007年提出,考慮了特定二聚體互補匹配信息,并且引入了miRNA位點可接近性的概念.miRNA位點可接近性表示整個二聚體的最小自由能與互補匹配區(qū)域的原始能量之間的差值.用戶可以通過調(diào)整約束條件來選擇候選靶基因列表(種子序列最小長度,G∶U錯配與未配對個數(shù)).

    2.2 機器學習預測方法

    (1)TargetBoost[29]:該算法采用GPboost模型,基于miRNA與候選靶基因mRNA形成的miRNA-mRNA二聚體的序列互補匹配程度、熱動力學穩(wěn)定性、以及跨物種保守性等特征,來推測線蟲和果蠅的miRNA的靶基因.該算法使用的正樣本為36個實驗驗證的miRNA-mRNA靶向關系數(shù)據(jù),而用于訓練的負樣本數(shù)據(jù)集為300個隨機生成長度為30 nt的基因序列.

    (2)miTarget[30]:該算法使用支持向量機方法,以及徑向基函數(shù),預測miRNA靶基因.該方法考慮了miRNA-mRNA二聚體的結構特征、熱動力學性質(zhì)和堿基互補匹配程度等特征.miTarget算法沒有考慮靶基因在多種物種間的序列保守性.miTarget算法中用于訓練支持向量機的負樣本數(shù)據(jù)集,涉及83個經(jīng)過實驗驗證的miRNA-mRNA靶向關系,以及163個通過實驗數(shù)據(jù)推理得到的miRNA-mRNA靶向關系數(shù)據(jù),正樣本數(shù)據(jù)集包含152個miRNA-mRNA靶向關系數(shù)據(jù).

    (3)MiRTif[31]:該方法首先整合miRanda、PicTar以及TargetScan這3種miRNA靶基因預測方法得到的各種特征得分,然后使用支持向量機方法預測miRNA候選靶基因,核函數(shù)采用徑向基函數(shù).MiRTif方法使用到的正樣本數(shù)據(jù)集由195個實驗驗證的miRNA-mRNA靶向關系數(shù)據(jù)構成,負樣本數(shù)據(jù)集囊括21個生物學實驗驗證的和17個假定的miRNA-mRNA靶向關系數(shù)據(jù).

    (4)MTar[32]:該算法重點考慮了3類區(qū)域的靶基因結合位點:僅僅考慮5’端種子區(qū)域;以5’端種子區(qū)域為主;以及以3’端種子區(qū)域為主.首先計算相應區(qū)域miRNA-mRNA雙鏈的特征得分,然后采用人工神經(jīng)網(wǎng)絡方法計算預測miRNA靶向基因.Mtar算法使用了340個miRNA-mRNA靶基因數(shù)據(jù)作為正樣本數(shù)據(jù)集,隨機產(chǎn)生了400個miRNA-mRNA靶基因數(shù)據(jù)作為負樣本數(shù)據(jù)集.

    (5)miRTDL[33]:該算法主要考慮了miRNA-mRNA之間的堿基互補匹配性、可接近性以及序列的保守性等特征.本方法選擇了1 297個實驗驗證的miRNA靶基因數(shù)據(jù)作為正樣本,309個實驗驗證負例樣本.因為該方法中負樣本遠少于正樣本的數(shù)量,所以該方法首先利用約束松弛方法構建了均衡的正、負樣本數(shù)據(jù)集,然后再利用深度學習模型,推測miRNA的靶基因.

    miRNA從1993年被首次發(fā)現(xiàn),到成為研究熱點,持續(xù)至今已經(jīng)歷20多年的歷史,在這個過程中,圍繞miRNA,科研學者積累了大量的數(shù)據(jù)資源,并建立了數(shù)據(jù)庫供同行分享使用,以便人們能夠在miRNA研究道路上走的更好、更遠.關于miRNA的數(shù)據(jù)庫系統(tǒng)已有數(shù)十種甚至多達百種,本文總結了一些常用的數(shù)據(jù)庫(表1).這些數(shù)據(jù)庫中有些是關于miRNA功能注釋的(miRbase,miRDB),有些是有關miRNA靶基因信息預測的數(shù)據(jù)庫(miRanda,TargetScan,PITA等).

    通過回顧總結miRNA靶基因計算預測的工作,可以發(fā)現(xiàn)這些方法都有如下一些特點:miRNA與其靶向mRNA序列跨物種、跨基因間的保守性;miRNA的5p端在miRNA-mRNA互補匹配中的重要作用;miRNA與靶基因mRNA形成的雙鏈的二級結構以及能量特征.不同的預測方法有其一定的使用范圍和獨特優(yōu)點,相信隨著實驗驗證方法的不斷發(fā)展,各種計算預測算法也可以相互補充完善,從而提高靶基因預測的精確度.

    表1 miRNA靶基因預測數(shù)據(jù)庫

    3 miRNA調(diào)控機制

    Lee等[34]于1993年在做篩選線蟲時序相關基因的研究時,發(fā)現(xiàn)了一個長度僅有22nt的RNA小分子——lin-4,這個小分子不編碼蛋白質(zhì),而是以堿基互補配對的方式結合到lin-14基因的3’端非轉(zhuǎn)錄區(qū)域(3’UTR)的特定位點,抑制lin-14對應mRNA的翻譯,最終使得其蛋白質(zhì)合成的量明顯降低,因此,稱其為轉(zhuǎn)錄后抑制,這是首次發(fā)現(xiàn)miRNA具有調(diào)控基因表達功能的研究.當時人們并不能完全理解為何如此微小的RNA分子會具有調(diào)控基因表達的功能,隨后關于該領域的研究一直沉寂了7年之久才再次有新的發(fā)現(xiàn).2000年,Reinhart等[35]同樣是在研究線蟲C.elegan時序數(shù)據(jù)時,發(fā)現(xiàn)了另一個長度只有21 nt小分子RNA——let-7,經(jīng)過研究發(fā)現(xiàn),該小分子能夠靶向基因lin-41,從而導致該基因?qū)牡鞍踪|(zhì)表達量減少.后來,隨著不斷的深入研究,生物學家總結出miRNA的功能作用機制,即成熟的miRNA能夠通過堿基互補配對的方式結合到靶基因mRNA上,當種子序列完全互補匹配時會促進靶基因的降解,而當不完全互補匹配時則會抑制靶基因的翻譯生成對應的蛋白質(zhì).

    長期以來,人們一直認為miRNA對于靶基因的表達是起到負調(diào)控的作用,因此,在研究miRNA功能的時候會選擇性地忽略正相關的miRNA-Gene調(diào)控關系,僅僅考慮具有負相關的miRNA-Gene調(diào)控關系.然而也有一些打破傳統(tǒng)約束,揭示出令人驚喜的研究成果,Vasudevan實驗室2007年發(fā)表在Science雜志上的一篇文章,發(fā)現(xiàn)miRNA并非總是負調(diào)控基因的表達,在一些特定情況下miRNA能夠起到促進基因表達的作用.該課題組通過實驗發(fā)現(xiàn),miR369-3能夠引導AGO等蛋白與AREs(AU-rich elemens)相結合促進基因的翻譯,除此以外,該課題組還發(fā)現(xiàn)在靜態(tài)細胞中即G0期狀態(tài)下let-7能夠與人工合成的miRcxcr4協(xié)同促進靶基因的翻譯表達,而在增殖細胞中它們對靶基因又起到抑制作用[12].2017年復旦大學于文強課題組在RNA Biology上發(fā)表了一篇文章,表明miRNA在細胞核內(nèi)可以與增強子結合,從而改變增強子的染色質(zhì)狀態(tài),進一步達到激活基因的轉(zhuǎn)錄表達的結果[13].曾獲得諾貝爾生理學獎的Phillip Sharp同樣于2017年在Cell上提出了miRNA激活理論[36],文章指出miRNA在細胞核內(nèi)則會與增強子相互作用,這也一定程度上肯定了復旦大學于文強課題組的研究成果.

    4 疾病相關miRNAs篩選方法

    已有研究表明,包括腫瘤在內(nèi)的許多復雜疾病的發(fā)生、發(fā)展過程中往往伴隨著編碼基因以及非編碼基因尤其是miRNA的差異化表達,如何識別在疾病進展過程中異常表達的miRNAs,對于疾病的診斷以及新藥的開發(fā)具有十分重要的意義.生物學家通過實驗的方法發(fā)現(xiàn)了一些疾病相關的miRNAs,哈爾濱工業(yè)大學蔣慶華課題組[37],北京大學崔慶華課題組[38],以及Yang等[39]基于文獻搜索方式,人工整理了與疾病相關的miRNAs(表2).除了通過實驗的方法以外,近年來關于篩選疾病相關miRNAs,出現(xiàn)了多種基于計算的方法,本節(jié)主要分3個部分來敘述:基于miRNA表達濃度的分析方法;基于miRNA表達相對位置排序的分析方法;基于miRNA功能相似性網(wǎng)絡分析方法.

    表2 人類疾病相關的miRNAs數(shù)據(jù)庫

    (1)基于miRNA表達差異分析方法

    直接法:直接法的意思是直接基于miRNA表達水平做差異表達分析,該方法假設如果某個miRNA差異表達,那么它對應的下游靶基因就是差異表達的.分析miRNAs表達譜最簡單常用的方法是倍數(shù)法(FC,F(xiàn)old Change)以及T-test差異檢驗方法[40-41].倍數(shù)法是基于miRNAs表達濃度計算變化的倍數(shù),選取變化大的miRNAs作為候選的疾病相關miRNAs.倍數(shù)方法操作簡單方便,同時它也有顯而易見的缺陷,就是這種方法缺少差異顯著性估計值,使得倍數(shù)閾值的選擇比較隨意,會受到主觀意愿的影響(到底是1.5倍、2倍或者還是3倍?).T-test差異檢驗方法可以彌補倍數(shù)法的不足,該方法可以通過設置顯著性閾值來獲得表達差異比較顯著的miRNAs.T-test差異檢驗方法是基于miRNAs的表達濃度呈正態(tài)分布假設為前提的.不管是倍數(shù)法還是T-test差異檢驗方法,他們都會受到實驗批次效應(Batch effect)和數(shù)據(jù)標準化過程的影響[42].

    間接法:miRNA與靶基因mRNA之間是一種“多對多”靶向調(diào)控關系,因此,簡單地通過單個miRNA的表達變化情況來決定其靶基因是否差異變化存在明顯的不足之處,因此,一些綜合的分析方法被提出,即通過多個miRNAs的協(xié)同作用間接確定下游靶基因的差異變化.東南大學肖忠黨課題組于2016年提出了基于miRNA表達濃度以及miRNA-mRNA結合系數(shù),綜合評估m(xù)iRNA表達譜整體對下游靶基因的影響[43],為系統(tǒng)性篩選受miRNA表達譜全局性顯著調(diào)控的編碼基因提供了解決方案.Garcia-Garcia等[44]提出基于多個miRNAs表達變化疊加效應來確定下游靶基因是否差異表達,該模型中考慮了miRNA水平差異表達變化的方向(+/-)以及顯著性值P-value.

    (2)基于miRNA表達量的排序方法

    哈爾濱醫(yī)科大學顧云燕課題組等認為在個性化樣本內(nèi),一些miRNAs在正常樣本中是存在著相對穩(wěn)定排序的,而疾病狀態(tài)下有些miRNAs會發(fā)生排序的逆轉(zhuǎn),基于該設想,顧云燕課題組提出基于miRNA表達相對位置排序篩選差異變化的miRNAs.該課題組首先在正常樣本群體內(nèi)篩選具有穩(wěn)定表達順序的miRNAs基因?qū)?,然后檢驗個性化樣體中miRNA排序的差異變化來篩選差異表達的miRNAs[45].在編碼基因水平,已有研究指出,基于基因排序篩選差異表達基因的方法可以大大地降低對于實驗批次效應以及數(shù)據(jù)標準化過程帶來的負面影響[46-47].雖然該方法對于抵抗批次效應以及數(shù)據(jù)標準化過程帶來的負面影響,是個非常有效的方法,然而它卻是犧牲了表達譜中關于miRNA的精確的量化信息.

    (3)基于miRNA功能相似性網(wǎng)絡方法

    單個miRNA能夠靶向調(diào)控多個mRNAs的表達翻譯,這些mRNAs可能富集在某些GO功能模塊中.如果2個miRNAs能夠協(xié)同調(diào)控某一個公共的GO功能模塊,那么這2個miRNAs就被認為具有功能協(xié)同的作用.基于該思想,哈爾濱醫(yī)科大學李霞課題組等構建了miRNA-miRNA協(xié)同調(diào)控網(wǎng)絡,從miRNAs協(xié)同調(diào)控網(wǎng)絡的角度探究腫瘤疾病相關的miRNAs,為開發(fā)miRNAs藥物或者發(fā)現(xiàn)新的治療靶點提供了一個新的視角[48].Li等[49]采用多個miRNAs協(xié)同調(diào)控靶基因mRNAs的思想,基于miRNA、mRNA表達譜數(shù)據(jù)提出了一種構建miRNA-mRNA調(diào)控網(wǎng)絡的分析方法.功能相似的miRNAs 調(diào)控的疾病也應該比較相似,基于這樣的假設,郭茂祖等[50]提出了BNPDCMDA算法.該方法首先構建了miRNA-疾病雙層網(wǎng)絡;然后,通過miRNAs 的功能相似度對其進行基于密度的聚類;最終,將聚類后的miRNAs與疾病構成的miRNA-疾病雙層網(wǎng)絡,采用二分網(wǎng)絡投影方法預測miRNA-疾病關聯(lián).

    除了以上3類從miRNAs表達譜數(shù)據(jù)出發(fā)逐步篩選疾病相關的miRNAs以外,也有一些從編碼基因表達譜數(shù)據(jù)出發(fā)的逆向反推導的方法.趙興明課題組僅利用編碼基因表達譜數(shù)據(jù),提出了篩選疾病相關miRNAs的研究方法[51].本課題組在前期的工作中提出基于基因子通路的方法篩選疾病相關miRNAs,該方法用到的數(shù)據(jù)是編碼基因表達譜數(shù)據(jù),首先對基因數(shù)據(jù)集做差異表達分析,然后對差異表達的基因做子通路富集分析,最后根據(jù)miRNA-mRNA作用關系篩選疾病相關的miRNAs[52].

    5 miRNA與小分子藥物

    已有研究表明,生物活性小分子或者藥物能夠調(diào)控miRNA的表達,這預示通過小分子藥物靶向調(diào)控miRNA或許能夠為疾病帶來新的治療方案.2008年Gumireddy實驗室開發(fā)了首個靶向miRNA的小分子抑制劑,該抑制劑被設計主要用來靶向抑制miR-21的表達[53].自此以后,人們利用高通量篩選或者通過計算機序列比對方法設計出大量的針對非編碼miRNA的小分子抑制劑.2012年哈爾濱醫(yī)科大學李霞課題組基于miRNA表達數(shù)據(jù)以及藥物對轉(zhuǎn)錄組的影響結果,提出一種方法,將miRNA與小分子藥物聯(lián)系起來,構建了一個網(wǎng)絡SMirN[54],為發(fā)現(xiàn)新的治療靶點或者治療藥物提供了輔助性工作.2013年哈爾濱醫(yī)科大學李霞課題組建立了一個miRNA與小分子藥物關系數(shù)據(jù)庫SM2miR,專門收集實驗驗證的能夠靶向調(diào)控miRNA表達的小分子藥物[55],該數(shù)據(jù)庫包含了2 925個作用關系(涉及151個小分子藥物,747個miRNAs以及17個物種).2017年上海同濟大學趙興明課題組構建了一個miRNA影響藥物治療效果的關系數(shù)據(jù)庫[56],該數(shù)據(jù)庫不僅存儲了miRNA與小分子藥物的對應關系,而且還給出了miRNA對藥物治療效果的影響.

    隨著miRNA相關研究的不斷深入,一些針對包括腫瘤在內(nèi)的多種疾病的miRNAs藥物被開發(fā)出來,并且陸續(xù)有一些miRNA藥物從實驗室走向臨床試驗階段(表3).miRNA小分子藥物相關數(shù)據(jù)庫的構建必將有效地加速推進miRNA在臨床醫(yī)學中的具體應用.

    表3 進入臨床試驗中的miRNAs藥物[57]

    6 總 結

    隨著miRNA研究的不斷深入以及高通量測序技術的發(fā)展,人們在靶基因預測、miRNA生物學功能機制的探究以及疾病相關miRNAs篩選等方面,發(fā)展了大量的分析研究方法,并發(fā)現(xiàn)了許多對疾病具有重要影響的miRNAs,其中一些對疾病治療極具潛力的miRNAs藥物已經(jīng)進入臨床試驗階段.

    雖然關于miRNA的研究已有十幾年的歷史,并且取得了喜人的豐碩成果,但是還有一些方面有待繼續(xù)跟進:①靶基因預測方法普遍存在高假陽性的缺點,靶基因預測準確性對于miRNA功能的研究至關重要,因此,亟待發(fā)展具有突破性的計算方法來預測miRNA的靶向基因;②miRNA對靶基因調(diào)控的復雜性,miRNA與靶基因之間是一種多對多的靶向關系,且miRNA對靶基因具有抑制或者促進的雙重作用,因此,如何準確量化miRNA對其靶基因的影響是一個具有挑戰(zhàn)性的重要課題;③加速推進miRNA藥物走向臨床,篩選疾病相關的miRNAs,其重要目的在于將其推向臨床醫(yī)學實際應用,可喜的是已有些miRNAs藥物進入臨床試驗階段,遺憾的是目前還沒有一個合格的miRNA藥物被批準進入市場.因此,有待新的有效篩選方法被提出,從而加速miRNA藥物早日進入臨床治病救人.

    猜你喜歡
    課題組靶向編碼
    陽城縣“耕心微寫”課題組
    如何判斷靶向治療耐藥
    中國臨床醫(yī)學影像雜志(2021年6期)2021-08-14 02:21:56
    基于SAR-SIFT和快速稀疏編碼的合成孔徑雷達圖像配準
    毛必靜:靶向治療,你了解多少?
    肝博士(2020年5期)2021-01-18 02:50:18
    《全元詩》未編碼疑難字考辨十五則
    原科技大學新能源開發(fā)與應用課題組介紹
    子帶編碼在圖像壓縮編碼中的應用
    電子制作(2019年22期)2020-01-14 03:16:24
    Genome and healthcare
    課題組成員
    支點(2015年11期)2015-11-16 10:25:03
    国产精品福利在线免费观看| 给我免费播放毛片高清在线观看| 欧美黑人欧美精品刺激| 成人av一区二区三区在线看| 不卡视频在线观看欧美| 老女人水多毛片| 麻豆av噜噜一区二区三区| 欧美激情在线99| 亚洲熟妇中文字幕五十中出| av在线观看视频网站免费| 看十八女毛片水多多多| 亚洲在线观看片| 国产亚洲欧美98| 俺也久久电影网| 国产欧美日韩精品亚洲av| 美女大奶头视频| 亚洲中文字幕日韩| 久久久久久久久久成人| 99热这里只有是精品50| 91麻豆av在线| 熟女人妻精品中文字幕| 亚洲黑人精品在线| 成年女人永久免费观看视频| 午夜a级毛片| 亚洲精品色激情综合| 亚洲最大成人av| 亚洲精品日韩av片在线观看| 亚洲av二区三区四区| 国产色婷婷99| 国产伦在线观看视频一区| 天堂影院成人在线观看| bbb黄色大片| 99久久精品国产国产毛片| 舔av片在线| 精品一区二区免费观看| 国产一级毛片七仙女欲春2| 亚洲成av人片在线播放无| 人妻久久中文字幕网| 99久久精品国产国产毛片| 日本免费a在线| 国产精品,欧美在线| 人人妻人人澡欧美一区二区| 亚洲av二区三区四区| 看片在线看免费视频| 色av中文字幕| 国产真实乱freesex| 国产伦人伦偷精品视频| 99久久中文字幕三级久久日本| 国产爱豆传媒在线观看| 国产精品乱码一区二三区的特点| 观看免费一级毛片| 99热这里只有是精品在线观看| 午夜福利在线观看吧| 综合色av麻豆| 国产精品av视频在线免费观看| 免费在线观看影片大全网站| 国产精品乱码一区二三区的特点| 最近最新中文字幕大全电影3| 午夜激情欧美在线| 久久午夜福利片| АⅤ资源中文在线天堂| 国产精品无大码| 欧美三级亚洲精品| 亚洲,欧美,日韩| 亚洲欧美日韩高清在线视频| 日本黄色片子视频| 男人和女人高潮做爰伦理| 国产黄色小视频在线观看| 亚洲av美国av| 精品久久久久久久久亚洲 | 18禁在线播放成人免费| 亚洲美女视频黄频| 悠悠久久av| 亚洲av五月六月丁香网| 亚洲在线自拍视频| 亚洲av不卡在线观看| 好男人在线观看高清免费视频| 天美传媒精品一区二区| 高清日韩中文字幕在线| 亚洲乱码一区二区免费版| 午夜福利在线观看免费完整高清在 | 欧美精品国产亚洲| 成人国产综合亚洲| 麻豆一二三区av精品| 国产精品久久电影中文字幕| 琪琪午夜伦伦电影理论片6080| 日日啪夜夜撸| 欧美绝顶高潮抽搐喷水| 日韩大尺度精品在线看网址| 国产精品久久电影中文字幕| 最新在线观看一区二区三区| 精品人妻熟女av久视频| 中国美女看黄片| 少妇裸体淫交视频免费看高清| 美女高潮喷水抽搐中文字幕| 男人舔女人下体高潮全视频| 国产高清视频在线播放一区| 国产精品久久久久久久电影| 黄色配什么色好看| 国产精品99久久久久久久久| 色综合婷婷激情| 窝窝影院91人妻| 成人精品一区二区免费| 久久久久久久久久成人| 国产精品一区二区三区四区免费观看 | 欧美国产日韩亚洲一区| 91精品国产九色| 欧美日韩亚洲国产一区二区在线观看| 免费观看精品视频网站| 中文字幕高清在线视频| 国产伦一二天堂av在线观看| 一个人看的www免费观看视频| 成人欧美大片| 亚洲性久久影院| 欧美+日韩+精品| 亚洲经典国产精华液单| 欧美区成人在线视频| 国产免费一级a男人的天堂| 老熟妇仑乱视频hdxx| 国产高潮美女av| 日本-黄色视频高清免费观看| 一个人看的www免费观看视频| 搡老熟女国产l中国老女人| 99久久精品一区二区三区| 欧美高清成人免费视频www| 亚洲性久久影院| 国产午夜精品论理片| 亚洲av二区三区四区| 欧美又色又爽又黄视频| eeuss影院久久| 亚洲自拍偷在线| 香蕉av资源在线| 精品午夜福利在线看| a级毛片a级免费在线| 综合色av麻豆| 综合色av麻豆| 黄色视频,在线免费观看| 日韩欧美免费精品| 国产精品爽爽va在线观看网站| 国产精品,欧美在线| 亚洲无线在线观看| 国产黄a三级三级三级人| 欧美日韩中文字幕国产精品一区二区三区| 有码 亚洲区| 国产黄a三级三级三级人| 天美传媒精品一区二区| 乱人视频在线观看| 校园人妻丝袜中文字幕| 精品久久久久久久末码| 国产精品久久视频播放| а√天堂www在线а√下载| 亚洲狠狠婷婷综合久久图片| 精品免费久久久久久久清纯| 最后的刺客免费高清国语| 国内精品久久久久久久电影| 精品久久久久久成人av| 日韩欧美在线二视频| 久久草成人影院| 亚洲人与动物交配视频| 亚洲一区高清亚洲精品| 国产在视频线在精品| 在线国产一区二区在线| 久久久久久久久久成人| 色综合色国产| 网址你懂的国产日韩在线| 内射极品少妇av片p| .国产精品久久| 国产精品一区二区三区四区免费观看 | 亚洲黑人精品在线| 亚洲国产精品成人综合色| 色播亚洲综合网| 毛片一级片免费看久久久久 | 日韩在线高清观看一区二区三区 | 日韩中字成人| 97热精品久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 色综合站精品国产| 亚洲内射少妇av| 变态另类丝袜制服| 日韩欧美一区二区三区在线观看| 村上凉子中文字幕在线| 国产免费一级a男人的天堂| 夜夜看夜夜爽夜夜摸| 日本黄色视频三级网站网址| 小说图片视频综合网站| av天堂中文字幕网| 精品人妻视频免费看| 我要看日韩黄色一级片| 嫩草影院新地址| 男人狂女人下面高潮的视频| 国产亚洲精品av在线| av中文乱码字幕在线| 国产成人福利小说| 听说在线观看完整版免费高清| av女优亚洲男人天堂| 男女那种视频在线观看| 黄色日韩在线| 99久久精品热视频| 男女之事视频高清在线观看| 午夜激情福利司机影院| 精品久久久久久久久久免费视频| 日韩欧美在线乱码| 久久久久久久久中文| 伦理电影大哥的女人| 日韩欧美 国产精品| 2021天堂中文幕一二区在线观| 男女下面进入的视频免费午夜| 国产精品精品国产色婷婷| 欧美色视频一区免费| 热99re8久久精品国产| 国内精品久久久久久久电影| 啦啦啦韩国在线观看视频| 内射极品少妇av片p| 国产精品1区2区在线观看.| av在线蜜桃| 无遮挡黄片免费观看| 久久婷婷人人爽人人干人人爱| 国产av不卡久久| 嫩草影院新地址| 国产亚洲av嫩草精品影院| aaaaa片日本免费| 亚洲美女搞黄在线观看 | 欧美高清性xxxxhd video| 色尼玛亚洲综合影院| 欧美人与善性xxx| 小说图片视频综合网站| 麻豆久久精品国产亚洲av| 精品久久久久久,| 美女黄网站色视频| 亚洲精品国产成人久久av| 成人午夜高清在线视频| 亚洲欧美清纯卡通| av女优亚洲男人天堂| 成人鲁丝片一二三区免费| 欧美国产日韩亚洲一区| 国产色婷婷99| netflix在线观看网站| 亚洲av二区三区四区| 免费在线观看影片大全网站| 国产三级在线视频| 三级国产精品欧美在线观看| 日韩人妻高清精品专区| .国产精品久久| 91久久精品电影网| 国产不卡一卡二| 哪里可以看免费的av片| 狠狠狠狠99中文字幕| 真人一进一出gif抽搐免费| 蜜桃亚洲精品一区二区三区| 亚洲人成网站高清观看| 国产精品99久久久久久久久| 国产老妇女一区| 神马国产精品三级电影在线观看| 国内揄拍国产精品人妻在线| 亚洲国产欧美人成| 丝袜美腿在线中文| 乱系列少妇在线播放| 一卡2卡三卡四卡精品乱码亚洲| 99精品在免费线老司机午夜| 欧美激情久久久久久爽电影| 亚洲一区二区三区色噜噜| 亚洲精品日韩av片在线观看| 亚洲av成人av| 91在线精品国自产拍蜜月| 亚洲四区av| 婷婷丁香在线五月| 99热精品在线国产| 网址你懂的国产日韩在线| 此物有八面人人有两片| 久久精品国产清高在天天线| 俄罗斯特黄特色一大片| 日本欧美国产在线视频| 91麻豆av在线| 波多野结衣巨乳人妻| 日本成人三级电影网站| 精品久久久久久久人妻蜜臀av| 国内精品久久久久精免费| 亚洲中文日韩欧美视频| 亚洲专区国产一区二区| 亚洲真实伦在线观看| 亚洲精品色激情综合| 欧美激情在线99| 非洲黑人性xxxx精品又粗又长| 中文在线观看免费www的网站| 国产成人一区二区在线| 人妻夜夜爽99麻豆av| 久久这里只有精品中国| 亚洲最大成人av| 午夜精品在线福利| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美在线乱码| 免费观看在线日韩| 日本成人三级电影网站| 国产精品久久久久久亚洲av鲁大| 国产蜜桃级精品一区二区三区| 国产一区二区三区在线臀色熟女| 五月玫瑰六月丁香| 欧美成人a在线观看| 久久6这里有精品| 久久午夜亚洲精品久久| 精品99又大又爽又粗少妇毛片 | 美女高潮喷水抽搐中文字幕| 一区二区三区四区激情视频 | 大型黄色视频在线免费观看| 日韩欧美一区二区三区在线观看| 日本免费a在线| 麻豆国产av国片精品| 久久久久久久久久久丰满 | 神马国产精品三级电影在线观看| 免费电影在线观看免费观看| 99国产精品一区二区蜜桃av| 国产成人aa在线观看| 内地一区二区视频在线| 亚洲无线观看免费| 91久久精品电影网| 亚洲一级一片aⅴ在线观看| 丝袜美腿在线中文| 国产69精品久久久久777片| 天堂av国产一区二区熟女人妻| 女同久久另类99精品国产91| av在线观看视频网站免费| 精品久久国产蜜桃| 给我免费播放毛片高清在线观看| 国产精品永久免费网站| 国产成人精品福利久久| av在线播放精品| 久久国内精品自在自线图片| 伊人久久国产一区二区| 成人高潮视频无遮挡免费网站| 91午夜精品亚洲一区二区三区| 亚洲自偷自拍三级| 午夜老司机福利剧场| 欧美极品一区二区三区四区| 中文资源天堂在线| 久久久亚洲精品成人影院| 啦啦啦在线观看免费高清www| 一本久久精品| 伊人久久国产一区二区| 精品酒店卫生间| 国产黄频视频在线观看| 夫妻午夜视频| 日韩亚洲欧美综合| 又黄又爽又刺激的免费视频.| .国产精品久久| 亚洲欧洲国产日韩| 国产成人免费无遮挡视频| 男人狂女人下面高潮的视频| 99re6热这里在线精品视频| av播播在线观看一区| 舔av片在线| 韩国av在线不卡| 日韩一区二区视频免费看| 亚洲欧美精品自产自拍| 国产一区亚洲一区在线观看| 我要看黄色一级片免费的| 老女人水多毛片| 99精国产麻豆久久婷婷| 精品久久久久久电影网| 美女福利国产在线 | 免费看日本二区| 内射极品少妇av片p| 亚洲欧美日韩无卡精品| 激情 狠狠 欧美| 人妻系列 视频| 尾随美女入室| 亚洲精品,欧美精品| 久久综合国产亚洲精品| 尤物成人国产欧美一区二区三区| 亚洲av国产av综合av卡| 性色avwww在线观看| 国产伦精品一区二区三区视频9| 免费观看性生交大片5| 亚洲国产毛片av蜜桃av| 亚洲av国产av综合av卡| 国产视频首页在线观看| 99热全是精品| 人人妻人人爽人人添夜夜欢视频 | 欧美高清成人免费视频www| 日日啪夜夜撸| 成人18禁高潮啪啪吃奶动态图 | 国产午夜精品一二区理论片| 午夜福利网站1000一区二区三区| 久久久国产一区二区| 男人爽女人下面视频在线观看| 一级毛片aaaaaa免费看小| 久久久亚洲精品成人影院| 久久ye,这里只有精品| 99热这里只有是精品在线观看| 亚洲国产精品999| 哪个播放器可以免费观看大片| 免费观看无遮挡的男女| 国产乱人偷精品视频| 九九在线视频观看精品| 日韩中文字幕视频在线看片 | 欧美最新免费一区二区三区| 亚洲av.av天堂| 国产精品嫩草影院av在线观看| 精品一品国产午夜福利视频| 亚洲激情五月婷婷啪啪| 国产精品精品国产色婷婷| 嘟嘟电影网在线观看| h日本视频在线播放| av免费在线看不卡| 久久久久网色| 99久久精品一区二区三区| av网站免费在线观看视频| 久久久精品94久久精品| 性色av一级| 亚洲中文av在线| 菩萨蛮人人尽说江南好唐韦庄| 精品久久久精品久久久| 男人添女人高潮全过程视频| 中文字幕av成人在线电影| 亚洲精品国产成人久久av| 99热这里只有是精品在线观看| 久久久久久久久久成人| 看非洲黑人一级黄片| 欧美xxxx黑人xx丫x性爽| 精品人妻一区二区三区麻豆| 成年美女黄网站色视频大全免费 | 久久婷婷青草| 免费人妻精品一区二区三区视频| av国产精品久久久久影院| 又爽又黄a免费视频| 永久免费av网站大全| 亚洲av电影在线观看一区二区三区| 国产精品伦人一区二区| 夫妻性生交免费视频一级片| 亚洲国产欧美在线一区| 亚洲精品第二区| 蜜桃在线观看..| 欧美日韩精品成人综合77777| 97热精品久久久久久| 日韩电影二区| 免费人妻精品一区二区三区视频| 嫩草影院新地址| 直男gayav资源| 欧美精品一区二区大全| 日本爱情动作片www.在线观看| 精品99又大又爽又粗少妇毛片| 国产av码专区亚洲av| 成人特级av手机在线观看| 如何舔出高潮| 夜夜看夜夜爽夜夜摸| 成人影院久久| 国产久久久一区二区三区| 国产高潮美女av| 亚洲精品久久午夜乱码| 尤物成人国产欧美一区二区三区| 欧美精品一区二区大全| 观看美女的网站| 老女人水多毛片| 婷婷色麻豆天堂久久| 国产免费视频播放在线视频| 黄色欧美视频在线观看| 欧美精品一区二区免费开放| 18禁裸乳无遮挡免费网站照片| 午夜视频国产福利| 亚洲成人一二三区av| 网址你懂的国产日韩在线| 免费久久久久久久精品成人欧美视频 | 大话2 男鬼变身卡| 国产av国产精品国产| 一级毛片aaaaaa免费看小| 一级片'在线观看视频| 涩涩av久久男人的天堂| 亚洲怡红院男人天堂| 韩国av在线不卡| 黑丝袜美女国产一区| 久久久国产一区二区| 特大巨黑吊av在线直播| 人妻夜夜爽99麻豆av| 精品久久久精品久久久| 午夜免费男女啪啪视频观看| 日本vs欧美在线观看视频 | 天堂俺去俺来也www色官网| 欧美日韩一区二区视频在线观看视频在线| 激情 狠狠 欧美| av在线老鸭窝| 成人黄色视频免费在线看| 最近手机中文字幕大全| 欧美最新免费一区二区三区| 在线观看国产h片| 久久这里有精品视频免费| 一边亲一边摸免费视频| 成年av动漫网址| 观看av在线不卡| 国产视频内射| 久久久午夜欧美精品| 观看美女的网站| 中文字幕久久专区| 国产欧美另类精品又又久久亚洲欧美| 久久久久久久久久久丰满| 久久精品久久久久久噜噜老黄| 国产精品伦人一区二区| 免费高清在线观看视频在线观看| 国产成人aa在线观看| 男女免费视频国产| 国产免费一区二区三区四区乱码| 日韩,欧美,国产一区二区三区| 欧美日韩视频高清一区二区三区二| 日韩电影二区| 亚洲国产欧美在线一区| 一个人免费看片子| 国产精品国产三级国产av玫瑰| 日日啪夜夜爽| 久久久久精品久久久久真实原创| 免费看日本二区| 一区二区三区精品91| 少妇人妻精品综合一区二区| 日韩人妻高清精品专区| 人妻少妇偷人精品九色| 91久久精品电影网| 国内精品宾馆在线| 你懂的网址亚洲精品在线观看| 丰满少妇做爰视频| av福利片在线观看| 亚洲欧美日韩卡通动漫| 精品少妇黑人巨大在线播放| 在线观看国产h片| 亚洲av中文字字幕乱码综合| 午夜激情久久久久久久| 男人添女人高潮全过程视频| 久久久久视频综合| 国产日韩欧美亚洲二区| 精品久久久久久久久亚洲| 狂野欧美激情性bbbbbb| 久久久久久久久久人人人人人人| 丰满迷人的少妇在线观看| 久久久久国产精品人妻一区二区| 久久久精品免费免费高清| 少妇被粗大猛烈的视频| 我要看黄色一级片免费的| 成人美女网站在线观看视频| h日本视频在线播放| 国产国拍精品亚洲av在线观看| 国产在视频线精品| 国产精品麻豆人妻色哟哟久久| 亚洲av国产av综合av卡| 亚洲人与动物交配视频| 国产黄片美女视频| 黑人猛操日本美女一级片| 亚洲人成网站在线播| 中文字幕av成人在线电影| 嫩草影院入口| kizo精华| 亚洲精品乱久久久久久| 国产精品久久久久久精品古装| 永久网站在线| 最近中文字幕高清免费大全6| 国产成人精品一,二区| 观看av在线不卡| 精品一区二区三区视频在线| 女人久久www免费人成看片| 激情五月婷婷亚洲| 亚洲人成网站高清观看| 精品国产乱码久久久久久小说| 三级国产精品片| 久久久亚洲精品成人影院| 国产精品熟女久久久久浪| 欧美最新免费一区二区三区| 免费观看无遮挡的男女| 亚洲在久久综合| 美女xxoo啪啪120秒动态图| 午夜视频国产福利| 国产男人的电影天堂91| 老师上课跳d突然被开到最大视频| 国产精品无大码| 三级国产精品片| freevideosex欧美| 国产黄色视频一区二区在线观看| av在线蜜桃| 日韩亚洲欧美综合| 午夜福利在线观看免费完整高清在| 日韩欧美精品免费久久| 波野结衣二区三区在线| 免费观看无遮挡的男女| 亚洲国产精品999| 久久人人爽av亚洲精品天堂 | 色哟哟·www| 成年女人在线观看亚洲视频| 一边亲一边摸免费视频| 日韩一区二区视频免费看| 色哟哟·www| 欧美变态另类bdsm刘玥| 丰满乱子伦码专区| kizo精华| 久久综合国产亚洲精品| 国产探花极品一区二区| 啦啦啦中文免费视频观看日本| 亚洲国产日韩一区二区| 卡戴珊不雅视频在线播放| 人人妻人人澡人人爽人人夜夜| videossex国产| 精品国产露脸久久av麻豆| 最近手机中文字幕大全| 成年美女黄网站色视频大全免费 | 久久热精品热| 少妇熟女欧美另类| 街头女战士在线观看网站| 如何舔出高潮| 免费播放大片免费观看视频在线观看| 免费黄网站久久成人精品| 99久久精品热视频| av在线老鸭窝| 2018国产大陆天天弄谢| av视频免费观看在线观看| 成人18禁高潮啪啪吃奶动态图 | 99热国产这里只有精品6| 蜜桃久久精品国产亚洲av| 精品久久久久久电影网| 国产精品麻豆人妻色哟哟久久| 国产亚洲91精品色在线| 日韩一区二区视频免费看| 国产高清国产精品国产三级 |