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

    在人轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò)中識別組織特異性功能模塊

    2015-04-29 20:39:16胡躍兵李彭平
    安徽農(nóng)業(yè)科學(xué) 2015年16期
    關(guān)鍵詞:模塊

    胡躍兵 李彭平

    摘要從人體的7個組織調(diào)控網(wǎng)絡(luò)出發(fā),通過MCODE算法獲得各個組織的調(diào)控模塊,再通過計算模塊間共有的節(jié)點(diǎn)數(shù)目獲得組織間相似的模塊或組織特異性模塊。結(jié)果表明,腎特異性模塊與腎癌密切相關(guān);腦特異性模塊在室間隔的形成和心肌細(xì)胞功能的實(shí)現(xiàn)中起著重要作用。

    關(guān)鍵詞共調(diào)控;轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò);模塊

    中圖分類號S03文獻(xiàn)標(biāo)識碼A文章編號0517-6611(2015)16-039-04

    Identification of Tissuespecific Moduels within Human Tissue Regulatory Networks

    HU Yuebing1, LI Pengping2

    (1. Jinling Hospital, Medical School of Nanjing University, Nanjing, Jiangsu 210002; 2. School of Life Science, Nanjing University, Nanjing, Jiangsu 210023)

    AbstractStarting from seven tissue regulatory network, using MCODE method, the regulatory module of each tissue was obtained. Through calculating common node number among modules, kidneyspecific module is closely correlated with renal cancer. The results showed that heartspecifc module is important for formation of ventricular septum and it is crucial for functions of cardiac cells.

    Key words Coregulation; Tanscriptional regulatory network; Module

    作者簡介

    胡躍兵(1983-),男,四川攀枝花人,碩士,從事神經(jīng)外科研究。

    收稿日期20150430

    基因表達(dá)是調(diào)控子共同作用使特定的mRNA維持在穩(wěn)定水平的過程。其中一類調(diào)控子是轉(zhuǎn)錄因子,它通過結(jié)合在其調(diào)控基因的啟動子區(qū),促進(jìn)或抑制該基因的轉(zhuǎn)錄。它在調(diào)控細(xì)胞生長和凋亡中起到重要作用。很多轉(zhuǎn)錄因子還可作為疾病治療的靶標(biāo)[1]。而microRNA (miRNA)則是通過抑制mRNA表達(dá)或促進(jìn)mRNA降解,在轉(zhuǎn)錄后水平調(diào)解基因的表達(dá)[2]。由此,轉(zhuǎn)錄因子和miRNA分別通過在轉(zhuǎn)錄水平和轉(zhuǎn)錄后水平調(diào)控基因表達(dá),從而在影響基因功能過程中發(fā)揮重要作用。而轉(zhuǎn)錄因子本身的表達(dá)在其轉(zhuǎn)錄后水平也受到miRNA的調(diào)控,而miRNA自身的轉(zhuǎn)錄也要接受轉(zhuǎn)錄因子的調(diào)控。因此,細(xì)胞內(nèi)基因的表達(dá)系統(tǒng)實(shí)際上是由轉(zhuǎn)錄因子、miRNA以及它們靶基因之間組成的一個復(fù)雜的調(diào)控網(wǎng)絡(luò)。許多研究表明,特定生物學(xué)功能的實(shí)現(xiàn)是由一個包含多種相似功能分子的功能模塊完成的。模塊內(nèi)部分子功能的相似度高于模塊之間分子功能的相似度,并且分子在模塊內(nèi)的連接度遠(yuǎn)遠(yuǎn)高于模塊間分子之間的連接度[3]。轉(zhuǎn)錄調(diào)控模塊是由不同調(diào)控子調(diào)控特定靶基因集或同一類調(diào)控子調(diào)控不同靶基因集而形成的一個調(diào)控子網(wǎng),是多個基因及其產(chǎn)物之間相互作用形成轉(zhuǎn)錄模塊的結(jié)果[4]。研究轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò)模塊結(jié)構(gòu)有利于進(jìn)一步了解細(xì)胞的轉(zhuǎn)錄調(diào)控、細(xì)胞分化及其功能。

    組織特異性基因的表達(dá)在組織的發(fā)育、細(xì)胞類型、功能和轉(zhuǎn)錄調(diào)控中起到重要作用[5]。在已經(jīng)識別的miRNA中,有很多表現(xiàn)出了組織特異性或發(fā)育階段特異性的表達(dá)模式,并維持組織的特點(diǎn)和功能[6]。并且由轉(zhuǎn)錄因子形成的調(diào)控網(wǎng)絡(luò)會控制組織特異性基因的表達(dá)。組織特異性的轉(zhuǎn)錄因子和miRNA對于組織特異性的調(diào)控和功能極其重要。基因轉(zhuǎn)錄表達(dá)過程實(shí)際上是調(diào)控子和靶基因之間形成功能模塊作用的結(jié)果。由單個基因的組織特異性的表達(dá)調(diào)控的研究上升到模塊的組織特異性調(diào)控模式的挖掘,可以進(jìn)一步加深對基因組織特異性轉(zhuǎn)錄調(diào)控和組織功能的理解。為此,筆者從人體的7個組織調(diào)控網(wǎng)絡(luò)出發(fā),通過MCODE算法獲得各個組織的調(diào)控模塊,再通過計算模塊間共有的節(jié)點(diǎn)數(shù)目獲得組織間相似的模塊或組織特異性模塊,從而獲得組織特異性的調(diào)控模式。

    1 數(shù)據(jù)來源與方法

    1.1各組織表達(dá)譜數(shù)據(jù)來源

    從UniGene[7]和 CGAP[8]獲取人體的大腦、心臟、腎臟、肝臟、卵巢、脾和睪丸7個組織的蛋白質(zhì)編碼基因(包括TF和nonTF)的表達(dá)譜數(shù)據(jù),而miRNA 的表達(dá)譜數(shù)據(jù)則取自Landgraf等[6]的研究。

    1.2網(wǎng)絡(luò)構(gòu)建

    通過以下方式構(gòu)建包括人體大腦、心臟、腎臟、肝臟、卵巢、脾和睪丸7個組織的轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò) (TRN):首先在人基因組范圍內(nèi)預(yù)測TF、miRNA以及非蛋白質(zhì)編碼基因之間的調(diào)控關(guān)系,獲得人的參考網(wǎng)絡(luò),然后通過對應(yīng)組織的表達(dá)譜數(shù)據(jù)得到在7個組織中表達(dá)的基因,并從參考網(wǎng)絡(luò)中獲取這些基因之間的調(diào)控關(guān)系從而構(gòu)建各個組織網(wǎng)絡(luò)。在參考網(wǎng)絡(luò)中,TF與靶基因之間的調(diào)控關(guān)系是先通過TRANSFAC數(shù)據(jù)庫[9]找到靶基因上游1 kb的序列,之后在該段序列內(nèi)尋找TF的靶點(diǎn),該研究選取的是5個物種(人、狗、牛、小鼠和負(fù)鼠)保守的結(jié)果。而miRNA與靶基因之間的調(diào)控關(guān)系是通過取Targetscan[10]、Pictar[11]和Tarbase[12]3個預(yù)測結(jié)果的交集得到。

    1.3組織模塊挖掘

    生物網(wǎng)絡(luò)可以直觀地看成一個圖形,圖形節(jié)點(diǎn)和邊分別表示生物分子和分子之間的關(guān)系。對于轉(zhuǎn)錄調(diào)控,則可以理解為一個有向的信號傳遞流,用有向圖表示。一般從圖中尋找子簇或局部緊密連接區(qū)域主要是依賴于網(wǎng)絡(luò)流最小割理論或譜聚類。MCODE[13]是用節(jié)點(diǎn)聚類系數(shù)大小作為節(jié)點(diǎn)權(quán)重的衡量,從網(wǎng)絡(luò)中獲取分子復(fù)合體的模塊挖掘方法。該算法主要分為節(jié)點(diǎn)權(quán)重計算、復(fù)合體預(yù)測、通過一定規(guī)則對復(fù)合體增加或減少節(jié)點(diǎn)3個階段。MCODE在挖掘模塊時會給每個模塊計算一個得分并按照該得分對模塊進(jìn)行排序。模塊得分越高,則排位越靠前,該模塊內(nèi)部連接也更加緊密。

    1.4模塊相似性得分

    用公式(1)計算模塊相似性,式中N1、N2分別表示模塊1和模塊2的節(jié)點(diǎn)數(shù),C表示2個模塊共有節(jié)點(diǎn)個數(shù)。當(dāng)2個模塊節(jié)點(diǎn)大小和類型完全相同時,R等于1,表示2個模塊完全相同。當(dāng)2個模塊之間的節(jié)點(diǎn)都不同時,共有節(jié)點(diǎn)數(shù)目C為0,此時R等于0。當(dāng)R越接近1時,模塊越相似。

    R=C2/N1×N2

    (1)

    得到模塊相似得分后,以該得分作為模塊間的距離,對模塊進(jìn)行系統(tǒng)聚類。聚類到一起的模塊就是組織間相似的模塊,而單個成簇的模塊就是組織特異性模塊。

    2 結(jié)果與分析

    2.1組織模塊

    通過MCODE算法對人體7個組織的調(diào)控網(wǎng)絡(luò)進(jìn)行挖掘,得到大腦、心臟、腎臟、肝臟、卵巢、脾和睪丸的模塊數(shù)目分別是20、11、15、12、16、16和15(表1)。得到模塊后,統(tǒng)計模塊中調(diào)控子的組成情況,發(fā)現(xiàn)所有組織大部分模塊的調(diào)控子都包含轉(zhuǎn)錄因子和miRNA,只有少數(shù)模塊里面的調(diào)控子只是轉(zhuǎn)錄因子或miRNA(表2),說明模塊的作用模式主要是轉(zhuǎn)錄因子和miRNA協(xié)同調(diào)控,這與前人的研究結(jié)果“轉(zhuǎn)錄因子和miRNA在乳腺癌和前列腺癌中共同調(diào)控基因的表達(dá)”[14]是一致的。

    表1網(wǎng)絡(luò)邊的大小和模塊數(shù)目

    組織網(wǎng)絡(luò)邊數(shù)目網(wǎng)絡(luò)模塊數(shù)目

    大腦63 07020

    心臟17 65011

    腎臟23 66015

    肝臟16 54212

    卵巢31 35916

    脾28 57516

    睪丸60 36215

    表2模塊的作用模式

    組織調(diào)控子只是轉(zhuǎn)錄因子調(diào)控子只是miRNA調(diào)控子包含轉(zhuǎn)錄因子和miRNA

    大腦3116

    心臟1010

    腎臟1014

    肝臟0111

    卵巢0016

    脾0115

    睪丸0114

    計算每2個組織每2個模塊之間的相似性得分。例如,計算大腦和心臟模塊的相似性,則將大腦20個模塊分別和心臟的11個模塊一一進(jìn)行計算。之后將組織相似性得分作為距離,對模塊進(jìn)行系統(tǒng)聚類。105個模塊共聚為40類(圖1)。其中,卵巢的第14個模塊、腎臟的第14和第15個模塊、心臟的第10個模塊都各自聚類成一個簇,說明該4個模塊和其他模塊相似性很低,因此將其定義為組織特異性模塊。

    注:黑色箭頭所指處為組織特異性模塊,分別為:卵巢的第14個模塊(o14)、腎臟的第15個模塊(k15)、腎臟的第14個模塊(k14)、心臟的第10個模塊(h10)。

    圖1根據(jù)組織模塊相似性得分對模塊聚類

    2.2腎臟組織特異性模塊

    腎臟的第15個模塊(圖2)包含調(diào)控子MTF1、hsamiR651、hsamiR708以及靶基因RAP1B和MTMR3。在該模塊中,hsamiR651和 hsamiR708共同調(diào)控RAP1B和MTMR3的表達(dá),同時hsamiR651又調(diào)控轉(zhuǎn)錄因子MTF1的表達(dá)。RAP1B是RASlike小GTP結(jié)合蛋白超家族成員之一,其主要功能是調(diào)控細(xì)胞粘附、生長、分化、遷移增殖以及整合素介導(dǎo)的細(xì)胞信號轉(zhuǎn)導(dǎo)[15]。有研究表明由RAP1B參與的cAMP/Rap1B/BRaf通路可以調(diào)控結(jié)節(jié)性硬化癥腎細(xì)胞癌p27的表達(dá)和 p27cyclin D1的細(xì)胞質(zhì)的錯誤定位,從而影響腎細(xì)胞癌的發(fā)生和癌癥的預(yù)后[16]。KEGG通路分析顯示,RAP1家族與腎細(xì)胞癌通路密切相關(guān)[17]。RAP1家族主要通過細(xì)胞粘附、增殖和遷移而參與到腎細(xì)胞癌通路中(圖3)。HsamiR708被發(fā)現(xiàn)在腎細(xì)胞癌中能夠誘導(dǎo)細(xì)胞凋亡并抑制腫瘤生長,是腎細(xì)胞癌中非常重要的腫瘤抑制子[18]。MTF1是一種核質(zhì)穿梭作為轉(zhuǎn)錄因子,聚集在細(xì)胞核內(nèi)并結(jié)合到含有金屬反應(yīng)元件(MRE)啟動子上,主要是誘導(dǎo)金屬硫蛋白及其他維持金屬穩(wěn)態(tài)相關(guān)蛋白的表達(dá)。而金屬硫蛋白主要功能是抵抗重金屬、清除自由基,許多研究表明其表達(dá)異常與腫瘤密切相關(guān),已經(jīng)作為結(jié)腸直腸癌和前列腺癌診斷的分子標(biāo)記物[19]。綜上所述,RAP1B和hasmiR708都與癌癥密切相關(guān),根據(jù)模塊中分子功能相似性特點(diǎn)推測 hsamiR651可能是和hasmiR708協(xié)同調(diào)控RAP1B和MTMR3的表達(dá)。而hasmiR651又調(diào)控轉(zhuǎn)錄因子MTF1的表達(dá),該轉(zhuǎn)錄因子又進(jìn)一步作用于金屬硫蛋白等相關(guān)蛋白,這些轉(zhuǎn)錄因子和miRNA 形成的調(diào)控通路可能在腎細(xì)胞癌的形成發(fā)展中起重要作用。該模塊中的節(jié)點(diǎn)及其調(diào)控關(guān)系需要通過試驗(yàn)進(jìn)一步證實(shí)。

    注:紅色圓圈表示轉(zhuǎn)錄因子;綠色圓圈表示miRNA;黃色圓圈表示靶基因;箭頭表示調(diào)控子對靶基因的調(diào)控。

    圖2腎臟特異性模塊——腎臟的第15個模塊

    圖3Rap1家族參與到腎細(xì)胞癌通路中

    2.3心臟組織特異性模塊

    心臟的第10個模塊由GSK3B、CPEB3、hsamiR196b和轉(zhuǎn)錄因子HAND1組成(圖4),在系統(tǒng)聚類中獨(dú)自成簇,是心臟的一個組織特異性模塊。HAND1是心臟特異性轉(zhuǎn)錄因子,主要功能是參與血管生成、發(fā)育以及肌肉組織形成。并且在心室隔和心臟形態(tài)學(xué)形成中也起著重要作用。 HAND1和其另外一個家族成員以互補(bǔ)的方式工作,調(diào)節(jié)右心室和主動脈弓動脈的形成,提示HAND1和其家族成員可能會介導(dǎo)為先天性心臟病。Hatemi等研究表明HAND1序列的變異在心房異構(gòu)的發(fā)生和發(fā)展中可能起重要作用[20]。而在先天性心臟病中也發(fā)現(xiàn)HAND1等基因的拷貝數(shù)變異,這些變異可能與心臟功能缺陷緊密相關(guān)[21]。

    GSK3B是一種絲氨酸蘇氨酸激酶,屬于糖原合酶激酶亞家族,是為數(shù)不多的信號分子,調(diào)節(jié)關(guān)鍵的細(xì)胞內(nèi)信號轉(zhuǎn)導(dǎo)通路。在大鼠的胚胎發(fā)育過程中,GSK3B的缺失會導(dǎo)致胚胎干細(xì)胞中受損的心肌細(xì)胞的分化并促使心肌細(xì)胞顯著增生,從而使心室被心肌細(xì)胞填充[22]。研究表明,GSK3B是調(diào)節(jié)心肌纖維化的重要調(diào)控子,心臟成纖維細(xì)胞中GSK3β特異性缺失會導(dǎo)致心臟纖維化、左心室功能障礙和在缺血性心臟中產(chǎn)生過多的疤痕[23]。

    HsamiR196b主要功能是參與粒細(xì)胞生成、脂質(zhì)代謝和造血功能[24]。而粒細(xì)胞升高會導(dǎo)致心內(nèi)膜炎和心肌病。CPEB3、hsamiR196b可能也與心室隔形成、心肌細(xì)胞功能相關(guān)。具體作用機(jī)制可能是hsamiR196b作用于轉(zhuǎn)錄因子HAND1,HAND1再作用于GSK3B,該調(diào)控通路在心室隔的形成、維持心肌細(xì)胞正常功能中發(fā)揮重要作用,需要通過試驗(yàn)進(jìn)一步證實(shí)。

    43卷16期

    胡躍兵等在人轉(zhuǎn)錄調(diào)控網(wǎng)絡(luò)中識別組織特異性功能模塊

    注:紅色圓圈表示轉(zhuǎn)錄因子;綠色圓圈表示miRNA;黃色圓圈表示靶基因;箭頭表示調(diào)控子對靶基因的調(diào)控。

    圖4心臟特異性模塊

    3結(jié)論

    通過構(gòu)建人體7個組織的調(diào)控網(wǎng)絡(luò),并對各調(diào)控網(wǎng)絡(luò)進(jìn)行模塊挖掘與比較分析,得到了與腎癌密切相關(guān)的腎臟特異性功能模塊,以及和室間隔的形成有關(guān)并參與心肌細(xì)胞功能的心臟特異性模塊。這些模塊中的基因以及基因之間的調(diào)控關(guān)系對于相應(yīng)組織疾病的治療可能有重要作用。

    參考文獻(xiàn)

    [1]

    WESTON V.Large scale screening for DNA damageinduced transcription factors as potential targets for treatment of CLL with p53 apoptotic defect[J].Blood,2009,114(22):1336.

    [2] ERHARD F,HAAS J,LIEBER D,et al.Widespread context dependency of microRNAmediated regulation[J].Genome Res,2014,24(6):906-919.

    [3] TORNOW S,MEWES H W.Functional modules by relating protein interaction networks and gene expression[J].Nucleic Acids Res,2003,31(21):6283-6289.

    [4] SEGAL E,F(xiàn)RIEDMAN N,KOLLER D,et al.A module map showing conditional activity of expression modules in cancer[J].Nat Genet,2004,36(10):1090-1098.

    [5] NIEHRS C,POLLET N.Synexpression groups in eukaryotes[J].Nature,1999,402(6761):483-487.

    [6] LANDGRAF P,RUSU M,SHERIDAN R,et al.A mammalian microRNA expression atlas based on small RNA library sequencing[J].Cell,2007,129(7):1401-1414.

    [7] SAYERS E W,BARRETT T,BENSON D A,et al.Database resources of the national center for biotechnology information[J].Nucleic Acids Res,2012,36:13-21.

    [8] LASH A E,TOLSTOSHEV C M,WAGNER L,et al.SAGEmap:A public gene expression resource[J].Genome Res,2000,10(7):1051-1060.

    [9] WINGENDER E.The TRANSFAC project as an example of framework technology that supports the analysis of genomic regulation[J].Briefings in Bioinformatics,2008,9(4):326-332.

    [10] LEWIS B P,SHIH I H,JONESRHOADES M,et al.Prediction of mammalian microRNA targets[J].Cell,2003,115(7):787-798.

    [11] KREK A,GRUN D,POY M N.Combinatorial microRNA target predictions[J].Nature Genetics,2005,37(5):495-500.

    [12] SETHUPATHY P,CORDA B,HATZIGEORGIOU A G.TarBase:A comprehensive database of experimentally supported animal microRNA targets[J].RNA,2006,12(2):192-197.

    [13] BADER G D,HOGUE C W.An automated method for finding molecular complexes in large protein interaction networks[J].BMC Bioinformatics,2003,4:2.

    [14] WU J H,SUN Y J,HSIEH P H,et al.Inferring coregulation of transcription factors and microRNAs in breast cancer[J].Gene,2013,518(1):139-144.

    [15] DUBE N,KOOISTRA M R,PANNEKOEK W J,et al.The RapGEF PDZGEF2 is required for maturation of cellcell junctions[J].Cellular Signalling,2008,20(9):1608-1615.

    [16] COHEN J D,THAM K Y,MASTRANDREA N J,et al.cAMPdependent cytosolic mislocalization of p27(kip)cyclin D1 during quinolthioetherinduced tuberous sclerosis renal ceu carcinoma[J].Toxicol Sci,2011,122(2):361-371.

    [17] KANEHISA M,GOTO S,SATO Y,et al.Data,information,knowledge and principle:Back to metabolism in KEGG[J].Nucleic Acids Research,2014,42(D1):199-205.

    [18] BIRKHAUSER F D,KOYA R C,NEUFELD C,et al.Dendritic cellbased immunotherapy in prevention and treatment of renal cell carcinoma:Efficacy,safety,and activity of AdGM.CAIX in immunocompetent mouse models[J].J Immunother,2013,36(2):102-111.

    [19] ARRIAGA J M,BRAVO I A,BRUNO L,et al.Combined metallothioneins and p53 proteins expression as a prognostic marker in patients with Dukes stage B and C colorectal cancer[J].Human Pathology,2012,43(10):1695-1703.

    [20] HATEMI A C,GULEC C C ,CINE N,et al.Sequence variations of NKX25 and HAND1 genes in patients with atrial isomerism[J].Anadolu Kardiyoloji Dergisithe Anatolian Journal of Cardiology,2011,11(4):319-328.

    [21] GOLDMUNTZ E,PALURU P,GLESSNER J,et al.Microdeletions and microduplications in patients with congenital heart disease and multiple congenital anomalies[J].Congenital Heart Disease,2011,6(6):592-602.

    [22] KERKELA R,KOCKERITZ L,MACAULAY K,et al.Deletion of GSK3beta in mice leads to hypertrophic cardiomyopathy secondary to cardiomyoblast hyperproliferation[J].J Clin Invest,2008,118(11):3609-3618.

    [23] LAL H,AHMAD F,WOODGETT J,et al.The GSK3 family as therapeutic target for myocardial diseases[J].Circulation Research,2015,116(1):138-149.

    [24] LU M,ZHANG Q,DENG M,et al.An analysis of human microRNA and disease associations[J].PLoS One,2008,3(10):3420.

    猜你喜歡
    模塊
    例談有機(jī)化學(xué)知識在其他模塊試題中的滲透
    28通道收發(fā)處理模塊設(shè)計
    “選修3—3”模塊的復(fù)習(xí)備考
    Module 4 Music Born in America
    基于NX的數(shù)控加工模塊設(shè)計與開發(fā)
    高精度模塊式自清洗過濾系統(tǒng)
    必修模塊相關(guān)知識過關(guān)訓(xùn)練
    三個必修模塊 教學(xué)各有側(cè)重
    選修6 第三模塊 International Relationships
    膜蒸餾模塊
    石油化工(2015年11期)2015-08-15 00:43:05
    黄色成人免费大全| 国产成人精品在线电影| 国产主播在线观看一区二区| 人人妻,人人澡人人爽秒播| 999久久久精品免费观看国产| 国产无遮挡羞羞视频在线观看| 97人妻天天添夜夜摸| av电影中文网址| 亚洲精品久久午夜乱码| 校园春色视频在线观看| 一进一出抽搐动态| 一级作爱视频免费观看| 中国美女看黄片| 精品免费久久久久久久清纯 | 久久香蕉激情| 成年人黄色毛片网站| 亚洲精华国产精华精| 99国产精品一区二区蜜桃av | 亚洲精品久久午夜乱码| 两个人看的免费小视频| 免费av中文字幕在线| 国产激情久久老熟女| 一级片'在线观看视频| 国产欧美日韩综合在线一区二区| 欧美亚洲 丝袜 人妻 在线| 在线看a的网站| 露出奶头的视频| 天堂动漫精品| 岛国毛片在线播放| 十八禁网站免费在线| 少妇被粗大的猛进出69影院| 国产精品久久久人人做人人爽| 在线视频色国产色| 欧美日韩黄片免| 亚洲五月天丁香| 极品教师在线免费播放| 国产av精品麻豆| 精品人妻熟女毛片av久久网站| 亚洲精品美女久久久久99蜜臀| 国产精品国产av在线观看| 在线永久观看黄色视频| 在线永久观看黄色视频| 一a级毛片在线观看| 91麻豆精品激情在线观看国产 | 成人免费观看视频高清| 亚洲精品国产区一区二| 超碰成人久久| 国产亚洲精品久久久久久毛片 | 麻豆国产av国片精品| www.熟女人妻精品国产| 老司机亚洲免费影院| 久久婷婷成人综合色麻豆| 在线观看一区二区三区激情| 黄色视频不卡| 在线观看免费日韩欧美大片| svipshipincom国产片| 国产精品乱码一区二三区的特点 | 亚洲在线自拍视频| 久久精品亚洲熟妇少妇任你| 不卡一级毛片| 久久精品熟女亚洲av麻豆精品| 成人亚洲精品一区在线观看| av网站免费在线观看视频| 三级毛片av免费| 人人妻人人添人人爽欧美一区卜| 中文字幕高清在线视频| 国产一区二区三区视频了| 另类亚洲欧美激情| 免费观看人在逋| 在线观看免费午夜福利视频| 欧美久久黑人一区二区| 色播在线永久视频| 亚洲精品国产区一区二| 欧美av亚洲av综合av国产av| 久久香蕉国产精品| 最新的欧美精品一区二区| 久久香蕉激情| 亚洲国产精品一区二区三区在线| 亚洲午夜精品一区,二区,三区| 中文字幕精品免费在线观看视频| 国产亚洲精品一区二区www | 啦啦啦免费观看视频1| 国产淫语在线视频| av欧美777| 国产精品香港三级国产av潘金莲| 国产精品久久久人人做人人爽| 久久中文看片网| 午夜福利,免费看| 五月开心婷婷网| 五月开心婷婷网| 91av网站免费观看| 又大又爽又粗| 欧美乱妇无乱码| 香蕉国产在线看| 国产成人一区二区三区免费视频网站| 很黄的视频免费| 老司机福利观看| 中文字幕最新亚洲高清| 国产成人啪精品午夜网站| 国产高清videossex| xxxhd国产人妻xxx| 久久午夜综合久久蜜桃| 免费在线观看日本一区| 国产成人精品久久二区二区免费| 99国产精品99久久久久| 成人18禁高潮啪啪吃奶动态图| 黄色a级毛片大全视频| 国产99久久九九免费精品| 欧美乱码精品一区二区三区| a级毛片黄视频| 侵犯人妻中文字幕一二三四区| 国产精品1区2区在线观看. | 午夜福利,免费看| 天堂√8在线中文| 色94色欧美一区二区| 久久99一区二区三区| 久9热在线精品视频| 欧美色视频一区免费| 性少妇av在线| 午夜久久久在线观看| 黄片播放在线免费| 免费一级毛片在线播放高清视频 | 国产精品自产拍在线观看55亚洲 | 天天躁日日躁夜夜躁夜夜| 欧美激情久久久久久爽电影 | а√天堂www在线а√下载 | 国产精品98久久久久久宅男小说| 亚洲熟女毛片儿| www.自偷自拍.com| 免费av中文字幕在线| 丝袜美腿诱惑在线| 另类亚洲欧美激情| 岛国毛片在线播放| 俄罗斯特黄特色一大片| 亚洲五月天丁香| 别揉我奶头~嗯~啊~动态视频| 女人被狂操c到高潮| 久久国产精品人妻蜜桃| 热re99久久国产66热| 一边摸一边抽搐一进一出视频| 在线观看免费视频日本深夜| 丝瓜视频免费看黄片| 精品视频人人做人人爽| 欧美激情极品国产一区二区三区| 精品国产乱码久久久久久男人| 亚洲av欧美aⅴ国产| 国产av又大| 首页视频小说图片口味搜索| 丰满迷人的少妇在线观看| 99久久99久久久精品蜜桃| 亚洲精品成人av观看孕妇| ponron亚洲| 精品一品国产午夜福利视频| 黑人巨大精品欧美一区二区mp4| 一边摸一边抽搐一进一出视频| 中文亚洲av片在线观看爽 | 国产精品久久久久久人妻精品电影| 国产精品久久视频播放| 久久久精品国产亚洲av高清涩受| 欧美中文综合在线视频| 亚洲在线自拍视频| 免费在线观看视频国产中文字幕亚洲| 麻豆成人av在线观看| 日韩制服丝袜自拍偷拍| 精品视频人人做人人爽| 亚洲第一青青草原| 久久人妻福利社区极品人妻图片| 久久久久国内视频| 黄色片一级片一级黄色片| 亚洲avbb在线观看| 人妻丰满熟妇av一区二区三区 | 亚洲熟女毛片儿| 91成人精品电影| 天堂动漫精品| 亚洲av欧美aⅴ国产| 午夜精品久久久久久毛片777| 欧美一级毛片孕妇| 怎么达到女性高潮| 两性午夜刺激爽爽歪歪视频在线观看 | 黑人操中国人逼视频| 日本黄色视频三级网站网址 | 精品熟女少妇八av免费久了| 两人在一起打扑克的视频| 国产高清videossex| 国产在视频线精品| 亚洲国产欧美日韩在线播放| 色在线成人网| 男女床上黄色一级片免费看| 国内毛片毛片毛片毛片毛片| 青草久久国产| 精品午夜福利视频在线观看一区| 两个人看的免费小视频| 不卡av一区二区三区| 欧美成人午夜精品| 欧美精品高潮呻吟av久久| 亚洲av成人av| 中亚洲国语对白在线视频| 夜夜夜夜夜久久久久| 国产成人免费无遮挡视频| 黑人巨大精品欧美一区二区mp4| 午夜精品国产一区二区电影| 少妇被粗大的猛进出69影院| 午夜福利在线免费观看网站| 中国美女看黄片| 精品午夜福利视频在线观看一区| 怎么达到女性高潮| 久热爱精品视频在线9| 伊人久久大香线蕉亚洲五| 久久中文看片网| 久久性视频一级片| 国产单亲对白刺激| 国产熟女午夜一区二区三区| 亚洲 欧美一区二区三区| 91字幕亚洲| 黑丝袜美女国产一区| 国产精品亚洲av一区麻豆| 国产免费男女视频| 啦啦啦免费观看视频1| av超薄肉色丝袜交足视频| 国产人伦9x9x在线观看| 身体一侧抽搐| 人人妻人人澡人人爽人人夜夜| 精品久久久久久电影网| 一级毛片精品| 成人av一区二区三区在线看| 亚洲精品国产精品久久久不卡| 国产精品偷伦视频观看了| a级片在线免费高清观看视频| 伊人久久大香线蕉亚洲五| 欧美在线黄色| 亚洲熟妇熟女久久| 久久天躁狠狠躁夜夜2o2o| 老司机午夜十八禁免费视频| 黄色丝袜av网址大全| 9191精品国产免费久久| 亚洲熟女毛片儿| 999久久久精品免费观看国产| 国产精品香港三级国产av潘金莲| 99精品欧美一区二区三区四区| a级毛片在线看网站| 老司机福利观看| 男男h啪啪无遮挡| 亚洲中文字幕日韩| 日韩欧美一区视频在线观看| 极品教师在线免费播放| 亚洲人成电影观看| √禁漫天堂资源中文www| 两性午夜刺激爽爽歪歪视频在线观看 | 久久久国产欧美日韩av| 欧美黑人欧美精品刺激| 久久精品91无色码中文字幕| 中出人妻视频一区二区| 中文字幕高清在线视频| 欧美 日韩 精品 国产| 一进一出好大好爽视频| 男女免费视频国产| 国产区一区二久久| 美女福利国产在线| 天天躁夜夜躁狠狠躁躁| 99re6热这里在线精品视频| 色精品久久人妻99蜜桃| 久久久国产一区二区| 搡老熟女国产l中国老女人| 亚洲av成人不卡在线观看播放网| 精品人妻熟女毛片av久久网站| 黄片大片在线免费观看| 高清av免费在线| 女人爽到高潮嗷嗷叫在线视频| 日韩欧美国产一区二区入口| 19禁男女啪啪无遮挡网站| 国产av又大| 91老司机精品| 久久久精品区二区三区| 欧美老熟妇乱子伦牲交| 十八禁人妻一区二区| 久久人人爽av亚洲精品天堂| 精品一区二区三区视频在线观看免费 | 国产精品免费大片| 免费日韩欧美在线观看| 涩涩av久久男人的天堂| 一个人免费在线观看的高清视频| 国产精品偷伦视频观看了| 最新美女视频免费是黄的| 一二三四在线观看免费中文在| 国产真人三级小视频在线观看| 天天添夜夜摸| 岛国在线观看网站| 国产视频一区二区在线看| 一a级毛片在线观看| 日本a在线网址| 成年女人毛片免费观看观看9 | 超碰97精品在线观看| 国产精品国产av在线观看| 亚洲欧美一区二区三区黑人| 夜夜夜夜夜久久久久| 久久精品国产综合久久久| 亚洲精品在线观看二区| 黄色视频,在线免费观看| 日本撒尿小便嘘嘘汇集6| 黄色怎么调成土黄色| 在线观看66精品国产| 精品人妻在线不人妻| 夜夜爽天天搞| 婷婷丁香在线五月| 曰老女人黄片| 久久精品aⅴ一区二区三区四区| 久久人妻熟女aⅴ| 国产97色在线日韩免费| 99国产精品一区二区三区| av天堂在线播放| 99精国产麻豆久久婷婷| 国产亚洲欧美98| 女人爽到高潮嗷嗷叫在线视频| 波多野结衣av一区二区av| 亚洲aⅴ乱码一区二区在线播放 | 国产精品免费一区二区三区在线 | 99久久精品国产亚洲精品| 成人精品一区二区免费| 乱人伦中国视频| 免费观看人在逋| 精品欧美一区二区三区在线| 久久精品亚洲精品国产色婷小说| 人人妻,人人澡人人爽秒播| 少妇裸体淫交视频免费看高清 | 免费日韩欧美在线观看| av线在线观看网站| 国产av精品麻豆| 欧美成人免费av一区二区三区 | 亚洲国产精品sss在线观看 | 国产成人系列免费观看| 亚洲黑人精品在线| 国产精品 欧美亚洲| 麻豆乱淫一区二区| 一本一本久久a久久精品综合妖精| 男人的好看免费观看在线视频 | 国产男女超爽视频在线观看| 欧美日韩福利视频一区二区| 欧美乱码精品一区二区三区| 日韩成人在线观看一区二区三区| 少妇裸体淫交视频免费看高清 | 国产在线观看jvid| 国产免费现黄频在线看| 精品国产一区二区久久| 国产成人精品久久二区二区免费| 国产日韩欧美亚洲二区| 狠狠狠狠99中文字幕| av网站在线播放免费| 变态另类成人亚洲欧美熟女 | 丰满迷人的少妇在线观看| 一级黄色大片毛片| 国产高清视频在线播放一区| 18禁国产床啪视频网站| 国产男女超爽视频在线观看| 两性夫妻黄色片| 中文字幕av电影在线播放| 亚洲 国产 在线| 99国产综合亚洲精品| 熟女少妇亚洲综合色aaa.| 老司机在亚洲福利影院| 日本精品一区二区三区蜜桃| 国产成人一区二区三区免费视频网站| 欧美乱码精品一区二区三区| 欧美在线黄色| 色精品久久人妻99蜜桃| 久久国产精品人妻蜜桃| 天天躁夜夜躁狠狠躁躁| 天天添夜夜摸| 国产亚洲欧美精品永久| 亚洲中文字幕日韩| videosex国产| 熟女少妇亚洲综合色aaa.| 色婷婷久久久亚洲欧美| 12—13女人毛片做爰片一| 亚洲精品国产区一区二| 美女福利国产在线| 十分钟在线观看高清视频www| 欧美另类亚洲清纯唯美| 亚洲成a人片在线一区二区| 精品国产国语对白av| 久热这里只有精品99| 如日韩欧美国产精品一区二区三区| 9191精品国产免费久久| 国产精品国产av在线观看| 久久这里只有精品19| 97人妻天天添夜夜摸| 在线天堂中文资源库| 欧美日韩亚洲国产一区二区在线观看 | 亚洲午夜精品一区,二区,三区| 欧美黑人欧美精品刺激| 精品欧美一区二区三区在线| 国产亚洲一区二区精品| 99精国产麻豆久久婷婷| 亚洲欧洲精品一区二区精品久久久| 十八禁高潮呻吟视频| 在线观看免费日韩欧美大片| 女人爽到高潮嗷嗷叫在线视频| 久久精品aⅴ一区二区三区四区| 欧美日韩一级在线毛片| 王馨瑶露胸无遮挡在线观看| 18禁国产床啪视频网站| 精品一区二区三区四区五区乱码| 欧美日韩精品网址| 精品久久久久久久毛片微露脸| 免费av中文字幕在线| 12—13女人毛片做爰片一| 欧美精品人与动牲交sv欧美| 国产1区2区3区精品| 欧美黑人精品巨大| 国产欧美日韩综合在线一区二区| 又黄又爽又免费观看的视频| 久久久久精品人妻al黑| 亚洲精品国产一区二区精华液| 美女视频免费永久观看网站| 丝袜美腿诱惑在线| 制服人妻中文乱码| 三上悠亚av全集在线观看| 在线观看午夜福利视频| 亚洲欧美一区二区三区久久| 久久中文看片网| 女人精品久久久久毛片| 飞空精品影院首页| 在线观看午夜福利视频| 亚洲国产精品合色在线| 国产精品av久久久久免费| 免费一级毛片在线播放高清视频 | 精品少妇久久久久久888优播| 村上凉子中文字幕在线| 亚洲五月婷婷丁香| 一本大道久久a久久精品| 人妻一区二区av| 久久久国产成人精品二区 | www日本在线高清视频| 亚洲av成人av| 校园春色视频在线观看| 男人舔女人的私密视频| 这个男人来自地球电影免费观看| 免费观看a级毛片全部| 亚洲国产中文字幕在线视频| 亚洲专区中文字幕在线| 国产精品 国内视频| 亚洲av成人一区二区三| 亚洲第一青青草原| 大码成人一级视频| 中文字幕制服av| 夜夜爽天天搞| 精品久久久精品久久久| 国产精品二区激情视频| 久久久久国产一级毛片高清牌| 欧美亚洲日本最大视频资源| 99国产精品一区二区蜜桃av | 亚洲精品国产精品久久久不卡| 国产精品免费大片| 好男人电影高清在线观看| 一级片免费观看大全| 久久精品国产综合久久久| 欧美黄色片欧美黄色片| 国产精品 欧美亚洲| 国产在线观看jvid| 国产极品粉嫩免费观看在线| 亚洲精品成人av观看孕妇| 91av网站免费观看| 日日夜夜操网爽| 久久久水蜜桃国产精品网| 淫妇啪啪啪对白视频| 天堂中文最新版在线下载| 99久久99久久久精品蜜桃| 亚洲精品国产精品久久久不卡| 免费观看精品视频网站| a级片在线免费高清观看视频| 制服人妻中文乱码| 亚洲视频免费观看视频| 国产片内射在线| 黄色女人牲交| 午夜久久久在线观看| 亚洲第一av免费看| 亚洲成a人片在线一区二区| 欧美亚洲日本最大视频资源| 国产精品 国内视频| 欧美乱色亚洲激情| 露出奶头的视频| 亚洲色图综合在线观看| a级片在线免费高清观看视频| 高清av免费在线| 巨乳人妻的诱惑在线观看| 欧美丝袜亚洲另类 | 日日爽夜夜爽网站| 精品乱码久久久久久99久播| 亚洲第一欧美日韩一区二区三区| 日韩欧美国产一区二区入口| 后天国语完整版免费观看| 女同久久另类99精品国产91| 啪啪无遮挡十八禁网站| 国产精品久久久久成人av| 99精国产麻豆久久婷婷| 真人做人爱边吃奶动态| 国产欧美日韩精品亚洲av| 欧美久久黑人一区二区| 国产精品香港三级国产av潘金莲| 成年女人毛片免费观看观看9 | 国产精品一区二区在线观看99| 久久人人爽av亚洲精品天堂| 人妻 亚洲 视频| 成熟少妇高潮喷水视频| 亚洲成人免费av在线播放| 欧美日韩一级在线毛片| 不卡一级毛片| 久久精品亚洲av国产电影网| 在线观看66精品国产| 国产成人精品久久二区二区91| 黑人操中国人逼视频| av有码第一页| 黄色怎么调成土黄色| 欧美黄色片欧美黄色片| 亚洲欧洲精品一区二区精品久久久| 美女午夜性视频免费| 中出人妻视频一区二区| 日本五十路高清| 亚洲人成伊人成综合网2020| 伊人久久大香线蕉亚洲五| 人人妻人人澡人人看| 身体一侧抽搐| 亚洲色图 男人天堂 中文字幕| 亚洲三区欧美一区| 日韩人妻精品一区2区三区| 国产精品1区2区在线观看. | 老司机福利观看| 欧美精品高潮呻吟av久久| 国产精品亚洲av一区麻豆| 久久久国产欧美日韩av| 女人精品久久久久毛片| 精品亚洲成a人片在线观看| 操美女的视频在线观看| av线在线观看网站| 免费在线观看黄色视频的| 动漫黄色视频在线观看| 超色免费av| 黄网站色视频无遮挡免费观看| 18禁裸乳无遮挡免费网站照片 | av免费在线观看网站| 又大又爽又粗| 美女视频免费永久观看网站| 国产国语露脸激情在线看| 极品少妇高潮喷水抽搐| 可以免费在线观看a视频的电影网站| 亚洲精品在线观看二区| 亚洲五月婷婷丁香| 久久国产精品影院| 91国产中文字幕| 国产精品乱码一区二三区的特点 | 中文字幕高清在线视频| 老司机影院毛片| 亚洲av成人一区二区三| 欧美日韩中文字幕国产精品一区二区三区 | 在线国产一区二区在线| 色94色欧美一区二区| 热re99久久国产66热| 精品国产一区二区三区久久久樱花| 国产亚洲精品一区二区www | 免费一级毛片在线播放高清视频 | 丰满饥渴人妻一区二区三| 国产精华一区二区三区| 激情视频va一区二区三区| 国产黄色免费在线视频| 丰满饥渴人妻一区二区三| 成人黄色视频免费在线看| 亚洲专区中文字幕在线| 99久久精品国产亚洲精品| 免费看a级黄色片| 国产亚洲欧美精品永久| 久久久精品免费免费高清| 国产精品 欧美亚洲| 国产一卡二卡三卡精品| 欧美久久黑人一区二区| 欧美日本中文国产一区发布| 中文字幕高清在线视频| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲中文日韩欧美视频| 国内毛片毛片毛片毛片毛片| 老司机深夜福利视频在线观看| 99国产精品一区二区蜜桃av | 高清黄色对白视频在线免费看| 757午夜福利合集在线观看| 啪啪无遮挡十八禁网站| 国产精品一区二区精品视频观看| 人成视频在线观看免费观看| 精品少妇一区二区三区视频日本电影| 一区二区三区精品91| 免费在线观看影片大全网站| 欧美色视频一区免费| 亚洲 国产 在线| 亚洲一区高清亚洲精品| 亚洲国产看品久久| 中文字幕最新亚洲高清| 亚洲精品美女久久久久99蜜臀| xxxhd国产人妻xxx| 人人妻人人澡人人爽人人夜夜| 91九色精品人成在线观看| 久久久国产精品麻豆| 一本大道久久a久久精品| 亚洲国产精品一区二区三区在线| a级毛片黄视频| 99riav亚洲国产免费| 亚洲精品在线美女| 女警被强在线播放| 波多野结衣av一区二区av| 美女扒开内裤让男人捅视频| 欧美午夜高清在线| 久久青草综合色| 精品国产一区二区三区久久久樱花| 女性被躁到高潮视频| 1024视频免费在线观看| 亚洲成人免费电影在线观看| 免费观看a级毛片全部| 国产日韩欧美亚洲二区| 欧洲精品卡2卡3卡4卡5卡区|