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

    lncRNA與疾病關(guān)聯(lián)關(guān)系預(yù)測研究進展

    2019-07-01 02:35:55徐志鑫王亞東
    智能計算機與應(yīng)用 2019年3期
    關(guān)鍵詞:計算方法疾病

    徐志鑫 王亞東

    摘 要:lncRNAs(long noncoding RNAs)是一類長度超過200個核苷酸的非編碼RNAs,越來越多的證據(jù)表明lncRNAs在許多生物過程中起關(guān)鍵作用。lncRNAs也與多種疾病的發(fā)生和發(fā)展有密切的聯(lián)系。近年來,很多研究人員關(guān)注于預(yù)測lncRNA與疾病的關(guān)聯(lián)預(yù)測的計算模型,這可以有效地減少生物實驗的時間和成本。本文總結(jié)了近年來有關(guān)lncRNA與疾病關(guān)聯(lián)預(yù)測的一些數(shù)據(jù)資源和有代表性的計算方法,并對當前存在的問題進行了分析,討論了未來的挑戰(zhàn)和發(fā)展方向。

    關(guān)鍵詞: lncRNAs;疾病;lncRNA-疾病關(guān)聯(lián);計算方法

    文章編號: 2095-2163(2019)03-0211-05 中圖分類號: TP311.13 文獻標志碼: A

    0 引 言

    隨著人類基因組計劃的完成,蛋白質(zhì)編碼基因僅占整個基因組的約1.5%,這意味著超過98%的人類基因組還不是編碼蛋白質(zhì)序列,這些基因組區(qū)域通常轉(zhuǎn)錄為非編碼RNA(ncRNA)[1]。在很長一段時間里,非編碼RNA在生物學(xué)上被認定是轉(zhuǎn)錄噪聲[2]。然而,越來越多的研究表明,非編碼RNA在許多生物過程中發(fā)揮著關(guān)鍵作用。其中,lncRNAs(long noncoding RNAs)就是一類長度超過200個核苷酸的非編碼RNAs。近年來,lncRNAs引起了研究人員的高度關(guān)注,因為已有證據(jù)指出lncRNA在多種生物機制中扮演重要角色,如表觀遺傳調(diào)控、染色質(zhì)重構(gòu)、基因轉(zhuǎn)錄、蛋白質(zhì)轉(zhuǎn)運、細胞分化、細胞運輸?shù)萚3]。此外,lncRNAs的突變和失調(diào)與肺癌、乳腺癌、卵巢癌、腎病、心血管疾病、自身免疫疾病等多種疾病的發(fā)生和發(fā)展有密切的關(guān)系[4]。盡管已經(jīng)取得了一定的研究成果,但對于生物學(xué)家和遺傳學(xué)家來說,確定那些未知的關(guān)聯(lián)關(guān)系仍可堪稱是一個重大的挑戰(zhàn)。傳統(tǒng)的生物實驗和臨床方法會耗費大量的時間和精力,因此時下的眾多研究者正致力于開發(fā)有效的計算方法來解決這一問題。近年來,研究人員提出了一些計算模型和方法來預(yù)測lncRNA與疾病的關(guān)聯(lián)關(guān)系。

    本文綜述了近年來涌現(xiàn)的一些預(yù)測疾病相關(guān)的候選lncRNA的算法和工具,闡述了lncRNA-疾病關(guān)聯(lián)關(guān)系預(yù)測的研究進展。預(yù)測方法主要分為2類:基于生物網(wǎng)絡(luò)的方法和基于機器學(xué)習(xí)的方法,本文重點探討了其中一些方法。最后,則總結(jié)了當前的關(guān)鍵問題和未來的研究方向。

    1 相關(guān)數(shù)據(jù)資源

    隨著生物數(shù)據(jù)的快速增長,研究人員建立了一些特定的數(shù)據(jù)庫來存儲和管理數(shù)據(jù)。本節(jié)擬集中推介一些和lncRNA、疾病相關(guān)的公共數(shù)據(jù)資源,詳情可見如下。

    1.1 lncRNA數(shù)據(jù)資源

    NONCODE(http://www.bioinfo.org/noncode/)[5]是一個包含幾乎所有非編碼RNA的綜合知識庫。其中涵蓋16種(人類、小鼠、牛、大鼠、黑猩猩、大猩猩、猩猩、酵母、果蠅、擬南芥等)和527 336個lncRNA轉(zhuǎn)錄本。

    LNCipedia(http://www.lncipedia.org/)[6]是一個人類lncRNA的綜合數(shù)據(jù)庫,通過整合不同的數(shù)據(jù),當前版本包含了120 353個人類lncRNA轉(zhuǎn)錄物。而且,也研發(fā)提供了預(yù)測蛋白質(zhì)編碼能力的實用工具。

    lncRNAdb(http://www.lncrnadb.org)[7]是一種知識資源,旨在提供真核生物非編碼RNA的全面注釋。在該數(shù)據(jù)庫中,每個條目包含lncRNA的序列、結(jié)構(gòu)、基因組語境、基因表達、亞細胞定位、保守性和功能信息。用戶可以搜索和下載lncRNA信息并提交新條目。

    LncRBase(http://bicresources.jcbose.ac.in/zhumur/lncrbase)[8]是人類和小鼠lncRNA的存儲庫,總共包含216 562個lncRNA轉(zhuǎn)錄物條目。在該數(shù)據(jù)庫中,已廣泛收集了基本的lncRNA轉(zhuǎn)錄特征、基因組位置、重疊小非編碼RNA、相關(guān)重復(fù)元件、相關(guān)印記基因和lncRNA啟動子信息。

    1.2 疾病數(shù)據(jù)資源

    Disease Ontology(http://disease-ontology.org)[9]是馬里蘭大學(xué)醫(yī)學(xué)院基因組科學(xué)研究所主辦的與人類疾病相關(guān)的生物醫(yī)學(xué)數(shù)據(jù)整合的開放數(shù)據(jù)庫。設(shè)計上包含了8 043個遺傳,發(fā)育和獲得性人類疾病,旨在通過疾病概念連接不同的數(shù)據(jù)集。

    Medical Subject Headings(https://www.ncbi.nlm.nih.gov/mesh)是用于索引生命科學(xué)文章的綜合受控詞匯表詞庫。有27 883個描述符和87 000個按照層次結(jié)構(gòu)樹排列的條目。

    OMIM(http://www.omim.org/)[10]是由約翰霍普金思大學(xué)醫(yī)學(xué)院McKusick-Nathans遺傳醫(yī)學(xué)研究所維護的人類基因和遺傳疾病在線知識庫。截至2017年5月26日,共有24 065個條目。

    1.3 lncRNA-疾病關(guān)聯(lián)關(guān)系相關(guān)數(shù)據(jù)資源

    LncRNADisease(http://www.cuilab.cn/lncrnadisease)[11]是手動收集的經(jīng)過實驗驗證的lncRNA-疾病關(guān)聯(lián)關(guān)系數(shù)據(jù)庫。當前版本包含大約3000個lncRNA-疾病對和475個lncRNA相互作用對,其中涉及大約2 000個出版文獻中的914個lncRNA和329種疾病。

    Lnc2Cancer(http://www.bio-bigdata.net/lnc2cancer)[12]是一個癌癥相關(guān)lncRNAs資源,搜集包含了666個人類lncRNA和97個人類癌癥之間的1 488條關(guān)聯(lián)條目。與此同時,還提供了lncRNA表達模式、實驗技術(shù)、簡要功能描述、原始參考和附加注釋信息。此外,也允許用戶搜索、瀏覽和下載數(shù)據(jù)以及向數(shù)據(jù)庫提交新數(shù)據(jù)等操作。

    MNDR(http://www.rna-society.org/mndr)[13]是哺乳動物非編碼RNA-疾病關(guān)系在線知識庫,當前版本包含了807個lncRNA相關(guān)、229個miRNA相關(guān)、13個piRNA相關(guān)和100個snoRNA相關(guān)的數(shù)據(jù)條目。

    2 交叉驗證

    為了評估算法的效果,對于lncRNA-疾病關(guān)聯(lián)關(guān)系預(yù)測,通常對研究提出的模型進行留一交叉驗證(LOOCV)或5-折交叉驗證。其中,LOOCV在已知的、經(jīng)過實驗驗證的lncRNA-疾病關(guān)聯(lián)關(guān)系上設(shè)計發(fā)生,關(guān)聯(lián)關(guān)系中的每一對作為測試樣本,而剩余的用于訓(xùn)練模型。如果測試樣本的排名高于特定閾值,就可以認為正確預(yù)測了這個lncRNA-疾病關(guān)聯(lián)關(guān)系對。有2種類型的LOOCV,分別為:全局LOOCV和局部LOOCV。對于全局LOOCV,測試樣本的排名被列入所有調(diào)查疾病的全部未知候選微生物-疾病關(guān)聯(lián)對中。對于局部LOOCV,測試樣本僅在給定疾病的其它未知的疾病相關(guān)lncRNA中進行排名。即2種LOOCV的主要區(qū)別在于是否考慮了所有被調(diào)查的疾病。而-5折交叉驗證將所有已知的、經(jīng)過驗證的lncRNA-疾病關(guān)聯(lián)對隨機分為5組。不是選擇一個lncRNA-疾病關(guān)聯(lián)對作為測試樣本,而是將這5組中的每一組輪流作為測試集,而其它4組作為訓(xùn)練集。為了直觀地評估性能,研究者一般選擇接受者操作特征曲線(ROC),這是評估二元分類模型的常用方法。研究中,敏感度(sensitivity)和特異度(specificity)是ROC曲線中使用的2個關(guān)鍵指標。對于lncRNA-疾病關(guān)聯(lián)關(guān)系的預(yù)測,敏感度表示出一個測試被正確識別的比例,而特異度表示出一個測試被錯誤地認為有關(guān)聯(lián)關(guān)系的比例。利用這種方式,通過不斷改變閾值,使用真陽性率(敏感度)與假陽性率(1-特異度)繪制ROC曲線來做出比較。ROC曲線下的面積也通常用于測試性能。一般地,AUC=0.5表示效果是隨機的,AUC=1表示效果完美。

    3 計算方法

    越來越多的證據(jù)表明,lncRNA的突變和失調(diào)與多種人類疾病密切相關(guān),lncRNA功能及其與人類疾病的關(guān)系引起了更多研究者的興趣與重視。計算模型可能是識別潛在lncRNA功能和lncRNA-疾病關(guān)聯(lián)關(guān)系的有效方法。本節(jié)論述了一些用于識別lncRNA-疾病關(guān)聯(lián)的計算方法,詳情可見如下。

    3.1 基于生物網(wǎng)絡(luò)的方法

    具有相似功能lncRNA可能與表型相似的疾病相關(guān),根據(jù)這個假設(shè),一些研究人員提出了數(shù)種基于生物網(wǎng)絡(luò)的方法來預(yù)測疾病相關(guān)的lncRNA,對此可做闡釋解析如下。

    Sun等人[14]基于lncRNA功能相似網(wǎng)絡(luò)提出了一種基于全局網(wǎng)絡(luò)的計算方法RWRlncD。研究中,在相繼構(gòu)建了lncRNA-疾病關(guān)聯(lián)網(wǎng)絡(luò)、疾病相似性網(wǎng)絡(luò)和lncRNA功能相似網(wǎng)絡(luò)后,RWRlncD通過在lncRNA功能相似網(wǎng)絡(luò)上進行隨機游走重啟(RWR)來預(yù)測潛在的lncRNA-疾病關(guān)系。基于已驗證的lncRNA-疾病關(guān)聯(lián)關(guān)系,RWRlncD在LOOCV下獲得0.822的AUC。然而此方法不能應(yīng)用于沒有任何已知相關(guān)lncRNA的疾病。當未來有更多l(xiāng)ncRNA-疾病關(guān)聯(lián)和更準確的lncRNA功能相似性度量時,RWRlncD的預(yù)測性能將會改善。

    此后,研究立足于有更多共同的miRNA相互作用的lncRNA傾向于與相似的疾病相關(guān)聯(lián)的假設(shè),Zhou等人[15]又提出了RWRHLD模型預(yù)測潛在的lncRNA-疾病關(guān)聯(lián)關(guān)系。RWRHLD將3個網(wǎng)絡(luò)(miRNA相關(guān)聯(lián)的lncRNA-lncRNA串擾網(wǎng)絡(luò)、疾病相似性網(wǎng)絡(luò)和已知的lncRNA疾病關(guān)聯(lián)網(wǎng)絡(luò))整合為一個網(wǎng)絡(luò),并在其上進行隨機游走。基于已知實驗驗證的lncRNA-疾病關(guān)聯(lián),RWRHLD在LOOCV下獲得0.871的AUC值。然而,RWRHLD僅適用于具有已知的lncRNA-miRNA相互作用的lncRNA,此外lncRNA串擾網(wǎng)絡(luò)和lncRNA-疾病關(guān)聯(lián)網(wǎng)絡(luò)的不完全覆蓋可能會降低模型預(yù)測性。

    Yang等人[16]基于已知疾病基因和lncRNA-疾病關(guān)聯(lián)構(gòu)建了編碼-非編碼基因-疾病二分網(wǎng)絡(luò),并進一步實現(xiàn)了該二分網(wǎng)絡(luò)上的傳播算法,以推測潛在的lncRNA-疾病關(guān)聯(lián)。該方法在LOOCV下獲得了0.788 1的AUC。然而,非編碼基因、蛋白質(zhì)編碼基因和lncRNA功能注釋之間相互作用的缺失卻影響了方法的性能。

    除了單層的網(wǎng)絡(luò),一些研究人員嘗試構(gòu)建了lncRNA-疾病的多級網(wǎng)絡(luò),并據(jù)此基于多級網(wǎng)絡(luò)來研究識別新的疾病相關(guān)lncRNA。

    Yao等人[17]提出了基于多層復(fù)合網(wǎng)絡(luò)來預(yù)測疾病相關(guān)lncRNA的算法LncPriCNet。研究通過將表型-表型相互作用、lncRNA-lncRNA相互作用以及基因-基因相互作用與疾病-ncRNA關(guān)系相結(jié)合來構(gòu)建復(fù)合網(wǎng)絡(luò),而后再使用隨機游走重啟算法(RWR)來預(yù)測候選的疾病相關(guān)lncRNA。當已知疾病的lncRNAs信息不足時,LncPriCNet依然表現(xiàn)良好。原因可能是多層復(fù)合網(wǎng)絡(luò)能夠支持更多的信息交互。

    Zhang等人[18]基于多種生物信息提出了LncRDNetFlow算法。該算法整合了多個網(wǎng)絡(luò),包括lncRNA相似性網(wǎng)絡(luò)、蛋白質(zhì)相互作用網(wǎng)絡(luò)、疾病相似性網(wǎng)絡(luò)以及相互之間的關(guān)聯(lián)網(wǎng)絡(luò),并在異構(gòu)網(wǎng)絡(luò)上利用流傳播算法來預(yù)測lncRNA-疾病關(guān)聯(lián),最終在LOOCV下達到0.841的AUC。此方法能夠在沒有已知關(guān)聯(lián)的情況下預(yù)測新的關(guān)聯(lián)關(guān)系。究其原因就在于其整合了蛋白質(zhì)的信息及其與lncRNAs和疾病的關(guān)聯(lián)。

    3.2 基于機器學(xué)習(xí)的方法

    機器學(xué)習(xí)對于預(yù)測疾病候選lncRNA有著重要作用,可根據(jù)已知疾病相關(guān)的lncRNA和未知的lncRNA特征來訓(xùn)練分類器,同時根據(jù)各項集之間的生物學(xué)特征差異來對候選lncRNA進行排名。

    Chen等人[19]開發(fā)了LRLSLDA(Laplacian Regularized Least Squares for LncRNA–Disease Association)計算模型,該模型是基于半監(jiān)督學(xué)習(xí)框架的預(yù)測潛在疾病相關(guān)lncRNA模型。模型基于一個假設(shè)——相似的疾病傾向于與功能相似的lncRNA相關(guān)聯(lián)。LRLSLDA結(jié)合了已知的疾病-lncRNA關(guān)聯(lián)關(guān)系和lncRNA表達譜,在留一交叉驗證(LOOCV)下獲得0.776的AUC,同時也不需要負樣本的信息,而負樣本往往難以獲得。但LRLSLDA仍有一些限制,例如,模型中出現(xiàn)很多參數(shù),如何選擇參數(shù)并未得到根本性的解決。此外,對同一個lncRNA-疾病關(guān)聯(lián)對會分別從lncRNA和疾病空間得到2個不同的分數(shù),如何高效地結(jié)合2個分數(shù)也已成為時下的研究課題。

    Zhao等人[20]使用已知的癌癥相關(guān)lncRNA,基于多元數(shù)據(jù)、基因組、調(diào)節(jié)物組和轉(zhuǎn)錄組的整合,開發(fā)了一種基于樸素貝葉斯分類器的模型,以識別新的癌癥相關(guān)lncRNA。該模型基于可公開獲得的多個癌癥類型的外顯子數(shù)據(jù)和小鼠的直系同源lncRNA knockdown數(shù)據(jù),采用10折交叉驗證進行評估。測試后表明該模型顯示出良好的效果,并成功鑒定707種潛在的癌癥相關(guān)lncRNA。文獻中使用的監(jiān)督分類器,如支持向量機(SVM)和樸素貝葉斯分類器的關(guān)鍵限制是需要負樣本的信息,而這在當前的研究中是無法獲得的。因此,總是隨機選擇未被標記的lncRNA-疾病關(guān)聯(lián)對作為負樣本,這就嚴重影響了預(yù)測性能。

    基于功能相似的lncRNA總是與相似的疾病相關(guān)聯(lián),Chen等人[1]通過計算與lncRNA相關(guān)的疾病組之間的語義相似性開發(fā)了2個新的模型LNCSIM1和LNCSIM2。這2個模型之間的差別就在基于疾病有向無環(huán)圖的疾病語義相似度的計算上,而這對如何有效表示不同疾病之間的關(guān)系是至關(guān)重要的。當疾病語義相似性和lncRNA功能相似性(由LNCSIM計算)與之前的LRLSLDA相結(jié)合時,得到新的lncRNA-疾病關(guān)聯(lián)預(yù)測模型LRLSLDA-LNCSIM,可以在更大程度上提升LRLSLDA的預(yù)測性能。

    Biswas等人[21]提出了基于矩陣分解的lncRNA-疾病關(guān)聯(lián)預(yù)測計算模型。模型整合了lncRNA-疾病關(guān)聯(lián)關(guān)系,經(jīng)過實驗驗證的基因-疾病關(guān)聯(lián)關(guān)系,基因-基因相互作用數(shù)據(jù)以及l(fā)ncRNA和基因的表達譜。使用非負矩陣分解方法預(yù)測lncRNA和疾病的相互作用,并使用雙聚類來識別lncRNA模塊。

    對于有監(jiān)督的機器學(xué)習(xí)算法,需要假設(shè)與疾病相關(guān)的lncRNA和不相關(guān)的lncRNA是分離的,然而已被證明與疾病相關(guān)的lncRNA數(shù)量仍然較少,而且?guī)缀鯖]有實驗?zāi)軌蜃C明那些關(guān)系不存在的,故而研究人員設(shè)計開發(fā)了一些基于半監(jiān)督學(xué)習(xí)的模型。機器學(xué)習(xí)方法面臨的挑戰(zhàn)是如何選擇有用的生物特征來訓(xùn)練分類器。因此,整合多個數(shù)據(jù)資源是提高性能的有效方法。但是,一些冗余或不相關(guān)的生物信息可能是無用的,甚至?xí)档托阅?。同時,不同的分類算法也可能僅適合不同的數(shù)據(jù)資源。因此,利用多種學(xué)習(xí)算法可以獲得更好的預(yù)測性能。

    3.3 其它方法

    在前述研究中分析可知,所有的計算模型均需要與已知的lncRNA-疾病關(guān)聯(lián)來實現(xiàn)預(yù)測。然而,即使是現(xiàn)在,已知的經(jīng)實驗驗證的lncRNA-疾病關(guān)聯(lián)關(guān)系仍然非常有限。因此,研究人員已轉(zhuǎn)而開始探索基于已知的疾病相關(guān)基因、miRNA和lncRNA與基因或miRNA之間的關(guān)系來預(yù)測lncRNA-疾病關(guān)聯(lián)。

    Liu等人[22]開發(fā)了一種不需要依賴已知lncRNA-疾病關(guān)系,而通過整合已知人類疾病基因和人類lncRNA與基因表達譜來預(yù)測潛在的人類lncRNA-疾病關(guān)聯(lián)的計算方法。在該方法中,將lncRNA分為2部分:組織特異性和非組織特異性lncRNA。首先根據(jù)不同組織中所有l(wèi)ncRNA的表達水平計算組織特異性評分。然后,對于組織特異性lncRNA,該計算框架推測這些lncRNA與各種人類組織相關(guān)的疾病間可能存在的關(guān)聯(lián)。此外,還可以基于疾病-基因關(guān)聯(lián)和基因-lncRNA共表達關(guān)系獲得非組織特異性lncRNA的相關(guān)疾病。該模型在LOOCV下獲得0.764 5的AUC,非組織特異性lncRNA的預(yù)測精度為0.89。然而,該方法仍然不能預(yù)測無相關(guān)基因記錄的疾病相關(guān)lncRNA。

    Li等人[23]提出了一種基于基因組位置的計算方法,運行后能夠更為全面地預(yù)測與血管疾病相關(guān)的人類lncRNA。10個被預(yù)測為與血管平滑肌細胞相關(guān)的lncRNA將被選擇用于下一步的實驗驗證以測試該方法的準確性。結(jié)果,10個中共有8個得到了證實。實驗結(jié)果表明該方法的可靠預(yù)測性能及鑒定新型lncRNA用于診斷和治療血管疾病的潛在價值。然而,該方法的應(yīng)用范圍尚屬有限,因為并非所有的lncRNA都具有相鄰基因,即使該lncRNA具有相鄰基因,也可能與其鄰近基因無功能相關(guān)性。

    另外,統(tǒng)計學(xué)方法也已應(yīng)用在lncRNA-疾病關(guān)聯(lián)的預(yù)測中。Chen[24]基于超幾何分布提出了HGLDA模型。通過整合疾病語義相似性,miRNA-疾病關(guān)聯(lián)和miRNA-lncRNA相互作用來計算lncRNA功能相似性。過程中,將測試lncRNA和疾病是否顯著共享可以與這兩者相互作用的常見miRNA,模型對每個lncRNA-疾病對進行超幾何分布測試。HGLDA在LOOCV下獲得0.762 1的AUC,然而,HGLDA不能應(yīng)用于那些沒有任何已知miRNA相互作用的lncRNA。

    4 結(jié)束語

    越來越多的證據(jù)表明非編碼RNA(lncRNA)在許多人類疾病中發(fā)揮重要作用。預(yù)測新的lncRNA-疾病關(guān)聯(lián)將有助于生物學(xué)家理解疾病的致病機理。此外,這也有助于人類疾病的診斷、預(yù)防和治療。最近,很多研究人員構(gòu)建計算模型來預(yù)測新的lncRNA-疾病關(guān)聯(lián),使得顯著減少生物實驗的時間和成本已成為現(xiàn)實可能。本文中,擷選了一些重要的公開可用的數(shù)據(jù)庫,然后,詳述了近年來一些預(yù)測疾病相關(guān)lncRNA的計算方法,包括基于生物網(wǎng)絡(luò)的方法、基于機器學(xué)習(xí)的方法以及其它類型的方法,這些方法取得了一些成果,但還需要后續(xù)更大力度投入來改進現(xiàn)有研究。很多方法整合了其它類型的生物數(shù)據(jù)來提高性能,然而,將適合與lncRNA-疾病關(guān)聯(lián)預(yù)測的多種數(shù)據(jù)整合起來仍然頗顯困難。因此,需要開發(fā)更具針對性的方法來有效整合這些數(shù)據(jù)。對于一些機器學(xué)習(xí)模型,負樣本往往難以獲得,而且在參數(shù)選擇、分類器組合等方面也存在一定問題與不足。在很多模型中,相似性計算對關(guān)聯(lián)預(yù)測的準確性有重要影響,因此,如何設(shè)計科學(xué)合理的方法來計算相似性,同時處理整合來自不同生物信息的相似性分值也是亟待深入探討研究的關(guān)鍵問題。

    參考文獻

    [1]CHEN X, YAN C C, LUO C, et al. Constructing lncRNA functional similarity network based on lncRNA-disease associations and disease semantic similarity[J]. Scientific Reports, 2015, 5:11338.

    [2] LOURO R, SMIRNOVA A S, VERJOVSKI-ALMEIDA S. Long intronic noncoding RNA transcription: Expression noise or expression choice?[J]. Genomics, 2009, 93(4):291-298.

    [3] GEISLER S, COLLER J. RNA in unexpected places: Long non-coding RNA functions in diverse cellular contexts[J]. Nature Reviews Molecular Cell Biology, 2013, 14(11):699-712.

    [4] CHEN X, YAN C C, ZHANG X, et al. Long non-coding RNAs and complex diseases: From experimental results to computational models[J]. Briefings in Bioinformatics, 2016, 18(4):558-576.

    [5] ZHAO Yi, LI Hui, FANG Shuangsang, et al. NONCODE 2016: An informative and valuable data source of long non-coding RNAs[J]. Nucleic Acids Research, 2016, 44(Database issue):D203-D208.

    [6] VOLDERS P J, HELSENS K, WANG X, et al. LNCipedia: A database for annotated human lncRNA transcript sequences and structures[J]. Nucleic Acids Research, 2013, 41(Database issue):D246-D251.

    [7] QUEK X C, THOMSON D W, MAAG J L, et al. lncRNAdb v2.0: Expanding the reference database for functional long noncoding RNAs[J]. Nucleic Acids Research, 2015, 43(Database issue):168-173.

    [8] CHAKRABORTY S, DEB A, MAJI R K, et al. LncRBase: An enriched resource for lncRNA information[J]. PloS One, 2014, 9(9):e108010.

    [9] LYNN S, ARZE C, NADENDLA S, et al. Disease ontology: A backbone for disease semantic integration[J]. Nucleic Acids Research, 2011, 40(Database issue):940-946.

    [10]AMBERGER J S, BOCCHINI C A, SCHIETTECATTE F, et al. OMIM.org: Online mendelian inheritance in man (OMIM), an online catalog of human genes and genetic disorders[J]. Nucleic Acids Research, 2015, 43(Database issue):789-798.

    [11]CHEN G, WANG Z, WANG D, et al. LncRNADisease: A database for long-non-coding RNA-associated diseases[J]. Nucleic Acids Research, 2013, 41(Database issue):983-986.

    [12]NING S, ZHANG J, WANG P, et al. Lnc2Cancer: A manually curated database of experimentally supported lncRNAs associated with various human cancers[J]. Nucleic Acids Research, 2016, 44(Database issue):D980-D985.

    [13]WANG Y, CHEN L, CHEN B, et al. Mammalian ncRNA-disease repository: A global view of ncRNA-mediated disease network[J]. Cell Death & Disease, 2013, 4(8):e765.

    [14]SUN Jie, SHI Hongbo, WANG Zhenzhen, et al. Inferring novel lncRNA-disease associations based on a random walk model of a lncRNA functional similarity network[J]. Molecular Biosystems, 2014, 10(8):2074-2081.

    [15]ZHOU M, WANG X, LI J, et al. Prioritizing candidate disease-related long non-coding RNAs by walking on the heterogeneous lncRNA and disease network[J]. Molecular bioSystems, 2015, 11(3): 760-769.

    [16]YANG X, GAO L, GUO X, et al. A network based method for analysis of lncRNA-disease associations and prediction of lncRNAs implicated in diseases[J]. PLoS One, 2014, 9(1): e87797.

    [17]YAO Qianlan, WU Leilei, JIA Li, et al. Global prioritizing disease candidate lncRNAs via a multi-level composite network[J]. Scientific Reports, 2017, 7:39516.

    [18]ZHANG Jingpu, ZHANG Zuping, CHEN Zhigang, et al. Integrating multiple heterogeneous networks for novel LncRNA-disease association inference[J]. IEEE/ACM Transactions on Computational Biology & Bioinformatics, 2017, PP(99):1.

    [19]CHEN Xing, YAN Guiying. Novel human lncRNA-disease association inference based on lncRNA expression profiles[J]. Bioinformatics, 2013, 29(20): 2617-2624.

    [20]ZHAO T, XU J, LIU L, et al. Identification of cancer-related lncRNAs through integrating genome, regulome and transcriptome features[J]. Molecular BioSystems, 2014, 11(1): 126-136.

    [21]BISWAS A K, KANG M, KIM D C, et al. Inferring disease associations of the long non-coding RNAs through non-negative matrix factorization[J]. Network Modeling Analysis in Health Informatics & Bioinformatics, 2015, 4(1):1-17.

    [22]LIU Mingxi, CHEN Xing, CHEN Geng, et al. A computational framework to infer human disease-associated long noncoding RNAs[J]. PLoS One, 2014, 9(1):e84408.

    [23]LI Jianwei, CHENG Gao, WANG Yuchen, et al. A bioinformatics method for predicting long noncoding RNAs associated with vascular disease[J]. Science China Life Sciences, 2014, 57(8):852-857.

    [24]Chen Xing. Predicting lncRNA-disease associations and constructing lncRNA functional similarity network based on the information of miRNA[J]. Scientific Reports, 2015, 5:13186.

    猜你喜歡
    計算方法疾病
    浮力計算方法匯集
    進擊的疾病
    尿碘與甲狀腺疾病的相關(guān)性
    易與豬大腸桿菌病混淆的腹瀉類疾病鑒別診斷
    立秋后要注意調(diào)養(yǎng)預(yù)防呼吸疾病
    夏季養(yǎng)生之疾病篇
    隨機振動試驗包絡(luò)計算方法
    基于能量的球壓痕硬度的計算方法
    不同應(yīng)變率比值計算方法在甲狀腺惡性腫瘤診斷中的應(yīng)用
    疾病與人
    av在线天堂中文字幕| 免费看日本二区| 久久热精品热| 日韩欧美国产在线观看| av视频在线观看入口| 久久国产乱子免费精品| 国产毛片a区久久久久| 久久久a久久爽久久v久久| 在线看三级毛片| 舔av片在线| eeuss影院久久| 身体一侧抽搐| 国产精品不卡视频一区二区| 成人鲁丝片一二三区免费| 欧美xxxx黑人xx丫x性爽| 久久精品久久久久久噜噜老黄 | 老师上课跳d突然被开到最大视频| 成人av在线播放网站| 亚洲精品国产成人久久av| 日本色播在线视频| 在线观看免费视频日本深夜| 啦啦啦韩国在线观看视频| 中出人妻视频一区二区| 美女 人体艺术 gogo| 91在线精品国自产拍蜜月| 亚洲av.av天堂| 久久久午夜欧美精品| 亚洲av五月六月丁香网| 欧美最黄视频在线播放免费| 女人被狂操c到高潮| 麻豆国产av国片精品| 久久久久久大精品| 日韩精品有码人妻一区| 日韩av不卡免费在线播放| 国产伦一二天堂av在线观看| av专区在线播放| 亚洲天堂国产精品一区在线| 一级av片app| 欧美色视频一区免费| 18禁在线无遮挡免费观看视频 | 三级毛片av免费| 少妇猛男粗大的猛烈进出视频 | 欧美成人a在线观看| 亚洲av成人av| 婷婷六月久久综合丁香| 亚洲av中文av极速乱| 国产视频一区二区在线看| 亚洲国产日韩欧美精品在线观看| 午夜久久久久精精品| 99热这里只有是精品在线观看| 国产亚洲av嫩草精品影院| 免费人成视频x8x8入口观看| 亚洲国产精品sss在线观看| 亚洲av五月六月丁香网| 久久久国产成人精品二区| 久久综合国产亚洲精品| 久久精品久久久久久噜噜老黄 | 亚洲欧美成人综合另类久久久 | 亚洲av美国av| av福利片在线观看| 国产伦一二天堂av在线观看| 国产久久久一区二区三区| 久久午夜亚洲精品久久| avwww免费| 久久久久久久久久成人| 国产高清三级在线| 日韩欧美精品免费久久| 亚洲性夜色夜夜综合| 日本熟妇午夜| 亚洲中文日韩欧美视频| 欧美高清性xxxxhd video| 麻豆成人午夜福利视频| 99久久无色码亚洲精品果冻| 国产精品三级大全| 精品乱码久久久久久99久播| 又粗又爽又猛毛片免费看| 一个人看的www免费观看视频| 麻豆乱淫一区二区| 女生性感内裤真人,穿戴方法视频| 黄片wwwwww| 日本黄大片高清| 欧美绝顶高潮抽搐喷水| 成人美女网站在线观看视频| 欧美色欧美亚洲另类二区| 一边摸一边抽搐一进一小说| 久久久久久久午夜电影| 村上凉子中文字幕在线| 精品一区二区三区视频在线观看免费| 免费av观看视频| 麻豆成人午夜福利视频| 久久6这里有精品| 十八禁网站免费在线| 久久亚洲精品不卡| 国产精品无大码| av在线蜜桃| 久久亚洲国产成人精品v| 日本在线视频免费播放| 精品午夜福利视频在线观看一区| 成人av在线播放网站| 俺也久久电影网| 日韩av在线大香蕉| 听说在线观看完整版免费高清| 高清日韩中文字幕在线| 国产三级在线视频| 在线观看av片永久免费下载| 亚洲专区国产一区二区| 亚洲欧美成人精品一区二区| 欧美国产日韩亚洲一区| 一进一出抽搐动态| 欧美成人a在线观看| 日本a在线网址| АⅤ资源中文在线天堂| av免费在线看不卡| www.色视频.com| 少妇人妻一区二区三区视频| 国产乱人视频| 亚洲第一区二区三区不卡| 国产在线男女| 亚洲av成人精品一区久久| 成熟少妇高潮喷水视频| 一级黄片播放器| 毛片一级片免费看久久久久| 亚洲自偷自拍三级| 国产私拍福利视频在线观看| 久久精品国产自在天天线| 亚洲精品亚洲一区二区| 12—13女人毛片做爰片一| 国产午夜福利久久久久久| 亚洲人成网站在线播放欧美日韩| 熟妇人妻久久中文字幕3abv| 波多野结衣巨乳人妻| 美女内射精品一级片tv| 乱人视频在线观看| 十八禁网站免费在线| 黑人高潮一二区| 亚州av有码| 12—13女人毛片做爰片一| 成人美女网站在线观看视频| 99久国产av精品| 欧美极品一区二区三区四区| a级毛色黄片| 男人的好看免费观看在线视频| 亚洲aⅴ乱码一区二区在线播放| 一区二区三区四区激情视频 | 乱系列少妇在线播放| 成人漫画全彩无遮挡| 天堂√8在线中文| 午夜福利18| 无遮挡黄片免费观看| 97碰自拍视频| 国产私拍福利视频在线观看| 国产老妇女一区| 国产黄片美女视频| 精品熟女少妇av免费看| 日韩一本色道免费dvd| 亚洲精品一区av在线观看| 国产不卡一卡二| 美女黄网站色视频| 日本爱情动作片www.在线观看 | 精品久久久久久久久av| 中文字幕熟女人妻在线| 免费人成视频x8x8入口观看| 插阴视频在线观看视频| 国产 一区 欧美 日韩| 大又大粗又爽又黄少妇毛片口| 久久精品国产亚洲av涩爱 | 老熟妇乱子伦视频在线观看| 久久精品国产清高在天天线| 亚洲欧美日韩高清专用| 亚洲人成网站在线观看播放| 日本五十路高清| 草草在线视频免费看| 女人十人毛片免费观看3o分钟| 国产精品伦人一区二区| 国产精品野战在线观看| 99国产极品粉嫩在线观看| 亚洲成人中文字幕在线播放| 性插视频无遮挡在线免费观看| 成年女人毛片免费观看观看9| 久久久久久久久久久丰满| 免费无遮挡裸体视频| 久久99热这里只有精品18| 亚洲精品久久国产高清桃花| 精品久久久噜噜| 老熟妇仑乱视频hdxx| 亚洲国产欧美人成| 色视频www国产| 日韩国内少妇激情av| 亚洲真实伦在线观看| 俄罗斯特黄特色一大片| 成人亚洲精品av一区二区| 亚洲人成网站在线播放欧美日韩| 久久久久久久午夜电影| 啦啦啦啦在线视频资源| 看十八女毛片水多多多| 国产蜜桃级精品一区二区三区| 麻豆国产97在线/欧美| 国产在线男女| 性欧美人与动物交配| 免费看日本二区| 国内久久婷婷六月综合欲色啪| 欧美中文日本在线观看视频| 免费观看的影片在线观看| 久久精品国产亚洲网站| 午夜免费激情av| 欧美xxxx性猛交bbbb| 国产精品嫩草影院av在线观看| 99热这里只有是精品在线观看| 又爽又黄a免费视频| 日韩欧美在线乱码| 两个人的视频大全免费| 亚洲美女视频黄频| 春色校园在线视频观看| 波多野结衣巨乳人妻| av在线观看视频网站免费| 国产又黄又爽又无遮挡在线| 精品一区二区三区视频在线观看免费| 中文字幕熟女人妻在线| 女的被弄到高潮叫床怎么办| 九九在线视频观看精品| or卡值多少钱| 免费看日本二区| 免费av不卡在线播放| 久久久久久久久久黄片| 三级毛片av免费| av在线亚洲专区| 精华霜和精华液先用哪个| 亚洲七黄色美女视频| 少妇熟女欧美另类| 久久久久国产网址| 91麻豆精品激情在线观看国产| 毛片女人毛片| 亚洲人成网站在线观看播放| 国产高清激情床上av| 一边摸一边抽搐一进一小说| 亚洲人成网站在线观看播放| 日韩中字成人| 免费观看的影片在线观看| 麻豆一二三区av精品| 老熟妇仑乱视频hdxx| 美女内射精品一级片tv| 国产伦精品一区二区三区四那| 亚洲自拍偷在线| 97在线视频观看| 99久久成人亚洲精品观看| 国产高潮美女av| 一个人免费在线观看电影| 男人狂女人下面高潮的视频| 国产免费男女视频| 99热网站在线观看| 久久国内精品自在自线图片| 日韩一区二区视频免费看| 久久人人精品亚洲av| 丝袜美腿在线中文| 真人做人爱边吃奶动态| 久久精品久久久久久噜噜老黄 | 国产麻豆成人av免费视频| 亚洲av成人精品一区久久| 欧美+亚洲+日韩+国产| 欧美一区二区亚洲| 国产精品永久免费网站| 国产黄片美女视频| 床上黄色一级片| 午夜激情福利司机影院| 成年版毛片免费区| 波多野结衣高清作品| 桃色一区二区三区在线观看| 尾随美女入室| 91av网一区二区| 一个人观看的视频www高清免费观看| 老女人水多毛片| 美女免费视频网站| 最后的刺客免费高清国语| 亚洲在线观看片| 成人毛片a级毛片在线播放| 男女做爰动态图高潮gif福利片| 99久久精品国产国产毛片| 中文字幕av在线有码专区| 亚洲国产精品sss在线观看| 国模一区二区三区四区视频| 国产精品电影一区二区三区| 日韩欧美精品v在线| 亚洲五月天丁香| 亚洲在线观看片| 婷婷精品国产亚洲av在线| 女人十人毛片免费观看3o分钟| 国产色爽女视频免费观看| av天堂在线播放| 嫩草影院入口| 久久久久久久亚洲中文字幕| 无遮挡黄片免费观看| 午夜老司机福利剧场| 精品久久久久久久人妻蜜臀av| 97人妻精品一区二区三区麻豆| 97热精品久久久久久| 亚洲丝袜综合中文字幕| 免费看a级黄色片| videossex国产| 成人一区二区视频在线观看| 久久国产乱子免费精品| 久久久国产成人免费| 成人永久免费在线观看视频| 日韩三级伦理在线观看| 久久久国产成人精品二区| 在线观看免费视频日本深夜| 久久久精品欧美日韩精品| 亚洲不卡免费看| 亚洲精品日韩av片在线观看| 欧美丝袜亚洲另类| 不卡视频在线观看欧美| 男女做爰动态图高潮gif福利片| 久久午夜亚洲精品久久| 黄色一级大片看看| 俄罗斯特黄特色一大片| 在线看三级毛片| 日韩国内少妇激情av| 五月玫瑰六月丁香| 免费看av在线观看网站| 日本黄色视频三级网站网址| 噜噜噜噜噜久久久久久91| 久久精品国产亚洲av天美| ponron亚洲| 免费大片18禁| 天堂影院成人在线观看| 2021天堂中文幕一二区在线观| 99riav亚洲国产免费| 在线播放无遮挡| 色综合亚洲欧美另类图片| 97超级碰碰碰精品色视频在线观看| 搞女人的毛片| 亚洲精品日韩av片在线观看| 国产高清有码在线观看视频| 欧美高清成人免费视频www| 婷婷六月久久综合丁香| 成年免费大片在线观看| 亚洲人成网站在线播放欧美日韩| 欧美性猛交╳xxx乱大交人| 国产av一区在线观看免费| 国语自产精品视频在线第100页| 国产亚洲91精品色在线| 亚洲第一电影网av| 国产熟女欧美一区二区| 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 国产精品人妻久久久影院| 如何舔出高潮| 久久精品人妻少妇| 男人的好看免费观看在线视频| 看黄色毛片网站| 国产精品,欧美在线| 日本一二三区视频观看| 国产激情偷乱视频一区二区| 久久久久久久久大av| 欧美3d第一页| 国产午夜福利久久久久久| 久久6这里有精品| 少妇的逼水好多| 免费在线观看影片大全网站| 九九久久精品国产亚洲av麻豆| 一级毛片我不卡| 尾随美女入室| 免费高清视频大片| 久久久久精品国产欧美久久久| 日韩人妻高清精品专区| 俺也久久电影网| 国产色婷婷99| 美女免费视频网站| 午夜精品在线福利| 老司机福利观看| 大型黄色视频在线免费观看| 亚洲国产精品成人综合色| 观看免费一级毛片| 亚洲欧美日韩无卡精品| 亚洲国产精品成人久久小说 | 麻豆国产97在线/欧美| 国产不卡一卡二| 少妇被粗大猛烈的视频| 亚州av有码| 少妇熟女欧美另类| 人人妻人人看人人澡| 欧美潮喷喷水| 欧美高清成人免费视频www| 小说图片视频综合网站| 日韩精品有码人妻一区| 久久午夜福利片| 深爱激情五月婷婷| 国产av麻豆久久久久久久| 悠悠久久av| 观看美女的网站| 老师上课跳d突然被开到最大视频| 91久久精品电影网| 乱人视频在线观看| 亚洲精品456在线播放app| 免费av不卡在线播放| 国产精品福利在线免费观看| 六月丁香七月| 精品久久久久久久久久久久久| 不卡视频在线观看欧美| 国产蜜桃级精品一区二区三区| 亚洲va在线va天堂va国产| 亚洲三级黄色毛片| 欧美日韩综合久久久久久| 成人高潮视频无遮挡免费网站| 国产一级毛片七仙女欲春2| 亚洲一级一片aⅴ在线观看| 免费观看的影片在线观看| 国产成人a区在线观看| 国产一级毛片七仙女欲春2| 精品久久久久久久久av| 一进一出抽搐动态| 菩萨蛮人人尽说江南好唐韦庄 | 俄罗斯特黄特色一大片| 精品午夜福利视频在线观看一区| 99国产极品粉嫩在线观看| 一本一本综合久久| av福利片在线观看| 亚洲av美国av| 日韩av在线大香蕉| 最后的刺客免费高清国语| 看免费成人av毛片| 免费看美女性在线毛片视频| 国产在线男女| 亚洲一级一片aⅴ在线观看| 亚洲欧美精品综合久久99| 久久亚洲国产成人精品v| 男女边吃奶边做爰视频| 亚洲在线自拍视频| 国产精品免费一区二区三区在线| 精品人妻一区二区三区麻豆 | 大香蕉久久网| 成人漫画全彩无遮挡| 成人美女网站在线观看视频| 97超视频在线观看视频| 精品久久久久久成人av| 国产av在哪里看| 国产真实乱freesex| 欧美不卡视频在线免费观看| 日韩成人av中文字幕在线观看 | 欧美日韩综合久久久久久| 两个人视频免费观看高清| 午夜免费男女啪啪视频观看 | 国产 一区 欧美 日韩| 国产v大片淫在线免费观看| 特大巨黑吊av在线直播| 少妇的逼水好多| 久久人人精品亚洲av| 久久久色成人| 色尼玛亚洲综合影院| 亚洲av二区三区四区| 日日撸夜夜添| 欧美成人一区二区免费高清观看| 国产aⅴ精品一区二区三区波| 国产精品亚洲一级av第二区| 天堂动漫精品| 少妇人妻精品综合一区二区 | 日日干狠狠操夜夜爽| 国产精品电影一区二区三区| 蜜臀久久99精品久久宅男| 一个人观看的视频www高清免费观看| 丝袜喷水一区| av中文乱码字幕在线| 少妇猛男粗大的猛烈进出视频 | 国产伦一二天堂av在线观看| 国产精品亚洲一级av第二区| 99久久精品热视频| 欧美日韩在线观看h| 午夜福利在线在线| 男女视频在线观看网站免费| 国产精品伦人一区二区| 久久鲁丝午夜福利片| 五月玫瑰六月丁香| 一级av片app| 国产精品人妻久久久久久| 97超视频在线观看视频| 在线免费观看不下载黄p国产| 国产爱豆传媒在线观看| 日韩亚洲欧美综合| 天堂av国产一区二区熟女人妻| 99久久久亚洲精品蜜臀av| 日韩精品中文字幕看吧| 神马国产精品三级电影在线观看| 亚洲最大成人av| 日本在线视频免费播放| 亚洲一区二区三区色噜噜| 国产精品爽爽va在线观看网站| 欧美极品一区二区三区四区| 在线观看66精品国产| 一区二区三区高清视频在线| 偷拍熟女少妇极品色| 亚洲欧美中文字幕日韩二区| 亚洲av不卡在线观看| 97碰自拍视频| 午夜亚洲福利在线播放| 国产91av在线免费观看| 熟女电影av网| 欧美bdsm另类| 午夜老司机福利剧场| 免费av观看视频| 99久久精品热视频| 狂野欧美白嫩少妇大欣赏| 国产成人freesex在线 | 91久久精品国产一区二区成人| 亚洲最大成人中文| 三级经典国产精品| 精品一区二区三区视频在线观看免费| 国产成人91sexporn| 亚洲最大成人av| 国产日本99.免费观看| 婷婷亚洲欧美| 久久99热这里只有精品18| 国内精品宾馆在线| 免费av观看视频| 久久精品国产清高在天天线| 国产淫片久久久久久久久| 免费av观看视频| 少妇丰满av| 亚洲熟妇中文字幕五十中出| 亚洲最大成人手机在线| 久久久a久久爽久久v久久| 麻豆久久精品国产亚洲av| 亚洲精品国产成人久久av| 欧美色欧美亚洲另类二区| 国产精品人妻久久久久久| 人人妻人人澡欧美一区二区| 亚洲高清免费不卡视频| 天天一区二区日本电影三级| 最近视频中文字幕2019在线8| 啦啦啦韩国在线观看视频| av免费在线看不卡| 欧美色欧美亚洲另类二区| 一级黄色大片毛片| 国产精品免费一区二区三区在线| 插阴视频在线观看视频| 99久久精品一区二区三区| 久久国产乱子免费精品| 国模一区二区三区四区视频| 一卡2卡三卡四卡精品乱码亚洲| 中文亚洲av片在线观看爽| 亚洲18禁久久av| 成人二区视频| 久久久精品大字幕| 精品一区二区三区视频在线观看免费| 床上黄色一级片| 久久久色成人| 中文字幕av在线有码专区| 免费看a级黄色片| 美女内射精品一级片tv| 亚洲人成网站高清观看| 在线观看av片永久免费下载| 色播亚洲综合网| 天堂影院成人在线观看| 国产女主播在线喷水免费视频网站 | 欧美xxxx性猛交bbbb| 亚洲性久久影院| 国产精品,欧美在线| 99热这里只有是精品在线观看| 欧美日韩综合久久久久久| 亚洲国产欧美人成| 国产91av在线免费观看| 国产女主播在线喷水免费视频网站 | 美女免费视频网站| 亚洲va在线va天堂va国产| 看黄色毛片网站| 少妇人妻一区二区三区视频| 午夜福利视频1000在线观看| 精品国内亚洲2022精品成人| 美女免费视频网站| 校园人妻丝袜中文字幕| 精品久久久久久久久久免费视频| 午夜影院日韩av| 中国国产av一级| 高清毛片免费看| 久久热精品热| 欧美在线一区亚洲| 国内精品久久久久精免费| 欧美日本亚洲视频在线播放| 久久国内精品自在自线图片| 欧美日本视频| 99热网站在线观看| av黄色大香蕉| 亚州av有码| av专区在线播放| 亚洲内射少妇av| 一级毛片aaaaaa免费看小| 99热只有精品国产| 老师上课跳d突然被开到最大视频| 亚洲一级一片aⅴ在线观看| 国产aⅴ精品一区二区三区波| 丰满的人妻完整版| 亚洲国产精品国产精品| 久久99热这里只有精品18| 亚洲久久久久久中文字幕| www日本黄色视频网| 亚洲欧美精品综合久久99| 国产高潮美女av| videossex国产| 色5月婷婷丁香| 国产精品无大码| 国产精品久久视频播放| 啦啦啦韩国在线观看视频| 看非洲黑人一级黄片| 91在线精品国自产拍蜜月| 人妻久久中文字幕网| 99视频精品全部免费 在线| 女生性感内裤真人,穿戴方法视频| 欧美色欧美亚洲另类二区| 久久精品91蜜桃| 男女之事视频高清在线观看| 啦啦啦观看免费观看视频高清| 亚洲欧美日韩无卡精品| 精品一区二区三区视频在线| 免费看光身美女| 亚洲国产欧洲综合997久久,| 校园春色视频在线观看| 国产乱人视频|