亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        人類(lèi)微小RNA-208a的靶基因預(yù)測(cè)及生物學(xué)過(guò)程和信號(hào)通路的生物信息學(xué)分析

        2016-02-17 05:26:34黎珊珊鐘國(guó)強(qiáng)蔣智淵溫麗娜
        廣西醫(yī)學(xué) 2016年12期
        關(guān)鍵詞:信息學(xué)心肌病生物學(xué)

        黎珊珊 鐘國(guó)強(qiáng) 蔣智淵 溫麗娜

        (廣西醫(yī)科大學(xué)第一附屬醫(yī)院心血管內(nèi)科,南寧市 530021,E-mail:xiaobaishanshan33@126.com)

        論著·基礎(chǔ)研究

        人類(lèi)微小RNA-208a的靶基因預(yù)測(cè)及生物學(xué)過(guò)程和信號(hào)通路的生物信息學(xué)分析

        黎珊珊 鐘國(guó)強(qiáng) 蔣智淵 溫麗娜

        (廣西醫(yī)科大學(xué)第一附屬醫(yī)院心血管內(nèi)科,南寧市 530021,E-mail:xiaobaishanshan33@126.com)

        目的 采用生物信息學(xué)技術(shù)預(yù)測(cè)人類(lèi)微小RNA-208a(hsa-miRNA-208a)靶基因及分析其可能參與的生物學(xué)過(guò)程及信號(hào)通路。方法 應(yīng)用miRbase數(shù)據(jù)庫(kù)和UCSC Genome Browser分析工具獲取hsa-miRNA-208a的染色體定位、堿基序列和物種保守性等基本信息,利用miRanda、miRDB、TargetScan和miRwalk進(jìn)行靶基因預(yù)測(cè),采用miRwalk網(wǎng)站對(duì)靶基因取交集,合并有文獻(xiàn)支持和經(jīng)實(shí)驗(yàn)證實(shí)的靶基因作為基因集,對(duì)基因集進(jìn)行功能富集分析(GO分析)和信號(hào)通路分析。結(jié)果 hsa-miRNA-208a序列在各物種間高度保守。對(duì)TargetScan、miRDB、miRanda和miRwalk預(yù)測(cè)的靶基因進(jìn)行交集后共得到16個(gè)靶基因(CPEB2、CSNK2A2、DLD、MTF2、FBXO28、LEP、SMAD4、NLK、GPR88、VPS13D、ZNF215、CHD9、SLC7A5、C19orf44、SLC45A3、MAP4K4),miRwalk、DIANA Lab TarBase檢索到已驗(yàn)證的靶基因分別為12個(gè)(CASP3、MRGPRX3、FPRL1、MED13、IL10、GATA4、HSH2D、MYH6、MYH7、TNNI3、CDKN1A、PAK3)和3個(gè)(SOX6、MTM1、TAB3)。GO分析結(jié)果顯示靶基因富集于心室肌組織發(fā)育、心室形態(tài)發(fā)育、心腔發(fā)育、心肌發(fā)育、心臟發(fā)育以及細(xì)胞發(fā)育等生物學(xué)過(guò)程(P<0.05);信號(hào)通路分析結(jié)果顯示靶基因富集于黏附連接、緊密連接、心肌收縮、Wnt信號(hào)通路以及病毒性心肌炎、肥厚型心肌病和擴(kuò)張型心肌病的相關(guān)信號(hào)通路中(P<0.05)。結(jié)論 hsa-miRNA-208a參與心臟生長(zhǎng)發(fā)育的生物學(xué)過(guò)程,與心肌疾病的發(fā)生密切相關(guān)。

        微小RNA-208a;人類(lèi);生物信息學(xué);靶基因;心臟生長(zhǎng)發(fā)育;心肌疾病

        微小RNA(microRNA,miRNA)是一類(lèi)長(zhǎng)度約20~24個(gè)核苷酸的內(nèi)源性、非編碼RNA,它通過(guò)影響mRNA的穩(wěn)定性和翻譯,參與多細(xì)胞生物表達(dá)基因的轉(zhuǎn)錄后調(diào)控。成熟miRNA整合入RNA誘導(dǎo)的沉默復(fù)合體(RNA-induced silencing complex,RISC),miRNA 5′端第2~8個(gè)核苷酸(種子序列)通過(guò)與靶基因mRNA 3′非編碼區(qū)序列完全或不完全互補(bǔ)聯(lián)系于RISC,導(dǎo)致靶基因mRNA降解或抑制其翻譯,從而在轉(zhuǎn)錄后水平負(fù)性調(diào)控靶基因的表達(dá)。近年來(lái),對(duì)miRNA的研究已成為生命科學(xué)領(lǐng)域中的一個(gè)重要方向。miRNA分布廣泛,參與的生物學(xué)過(guò)程復(fù)雜,調(diào)控的靶基因眾多,因此研究miRNA的功能及作用機(jī)制有重要意義。目前確定miRNA靶基因的方法主要包括生物信息學(xué)軟件預(yù)測(cè)和生物學(xué)實(shí)驗(yàn)方法[1]。生物信息學(xué)是一門(mén)交叉學(xué)科,它包含了生物信息的獲取、處理、儲(chǔ)存、分析和解釋等內(nèi)容,綜合運(yùn)用數(shù)學(xué)、計(jì)算機(jī)科學(xué)和生物學(xué)的各種工具來(lái)闡明和理解大量數(shù)據(jù)所包含的生物學(xué)意義[2]。預(yù)先進(jìn)行生物信息學(xué)研究能為研究者提供較大的幫助,為后續(xù)實(shí)驗(yàn)提供思路,提高實(shí)驗(yàn)效率[3]。已有研究證實(shí)miRNA-208a與心房顫動(dòng)、心肌肥厚、急性心肌梗死、查加斯病心肌病變、中東呼吸綜合征及肌強(qiáng)直性營(yíng)養(yǎng)不良等有關(guān)[4-9],但其生物功能仍有待進(jìn)一步闡明。本研究擬利用生物信息學(xué)方法對(duì)人類(lèi)miRNA(homo sapiens-miRNA)-208a進(jìn)行分析,為進(jìn)一步研究其生物功能提供線索。

        1 資料與方法

        1.1 資料來(lái)源 使用PubMed數(shù)據(jù)庫(kù)(http://www.ncbi.nlm.nih.gov/pubmed/)hsa-miRNA-208a相關(guān)文獻(xiàn)進(jìn)行檢索,使用miRbase(http://www.mirbase.org/)數(shù)據(jù)庫(kù)和UCSC Genome Browser(http://genome.ucsc.edu/)分析工具獲取hsa-miRNA-208a的染色體定位、堿基序列和物種保守性等基本信息。

        1.2 研究方法 綜合miRNAbase和UCSC Genome Browser 獲取hsa-miRNA-208a的基本信息。利用miRanda (http://www.microrna.org/microrna/home.do)、miRDB (http://www.mirdb.org/miRDB/)、TargetScan (http://www.targetscan.org/)、miRNAwalk(http://www.umm.uni-heidelberg.de/apps/zmf/mirwalk/)4種方法預(yù)測(cè)hsa-miRNA-208a靶基因后,使用miRNAwalk(http://www.umm.uni-heidelberg.de/apps/zmf/mirwalk/)對(duì)靶基因取交集。然后合并有文獻(xiàn)支持以及由miRwalk和DIANA Lab TarBase檢索到的經(jīng)實(shí)驗(yàn)證實(shí)的靶基因共同作為進(jìn)一步生物信息學(xué)分析的基因集合。使用基因本體(gene ontology,GO)分析網(wǎng)站(http://geneontology.org/)和DAVID(https://david.ncifcrf.gov/)數(shù)據(jù)庫(kù)對(duì)此基因集進(jìn)行功能富集分析,并利用京都基因與基因組百科全書(shū)(Kyoto Encyclopedia of Genes and Genomes,KEGG)數(shù)據(jù)庫(kù)(http://www.genome.jp/kegg/)進(jìn)行信號(hào)通路富集分析,即得到基因集在GO類(lèi)別上的高頻率注釋及有統(tǒng)計(jì)學(xué)意義的信號(hào)通路。以P<0.05為差異有統(tǒng)計(jì)學(xué)意義。

        2 結(jié) 果

        2.1 has-miRNA-208a染色體定位及物種保守性分析 hsa-miRNA-208a定位于染色體14q11.2。使用miRbase檢索到包括人類(lèi)在內(nèi)共10個(gè)不同物種的miRNA-208a成熟miRNA序列并進(jìn)行分析,miRNA-208a的“AUAAGACGAGCAAAAAGC”18個(gè)堿基序列在這10個(gè)物種中高度保守,見(jiàn)表1。

        表1 不同物種miRNA-208a成熟序列

        注:mmu:mus musculus(小鼠);rno:rattus norvegicus(大鼠);mdo:monodelphis domestica(負(fù)鼠);mml:macaca mulatta(獼猴);ptr:pan troglodytes(黑猩猩);bta:bos taurus(牛);cfa:canis familiaris(犬);eca:equus caballus(馬):ppy:pongo pygmaeus(猩猩)。

        2.2 hsa-miRNA-208a靶基因預(yù)測(cè)TargetScan、miRDB、miRanda和miRwalk軟件預(yù)測(cè) hsa-miRNA-208a靶基因分別為2 799、548、4 757、634個(gè),然后使用miRwalk對(duì)hsa-miRNA-208a靶基因進(jìn)行交集(基于miRanda、miRDB、miRwalk、TargetScan),共得到16個(gè)基因,分別為CPEB2、CSNK2A2、DLD、MTF2、FBXO28、LEP、SMAD4、NLK、GPR88、VPS13D、ZNF215、CHD9、SLC7A5、C19orf44、SLC45A3、MAP4K4。miRwalk中查到經(jīng)驗(yàn)證的靶基因有12個(gè),分別為CASP3、MRGPRX3、FPRL1、MED13、IL10、GATA4、HSH2D、MYH6、MYH7、TNNI3、CDKN1A、PAK3。DIANA Lab TarBase數(shù)據(jù)庫(kù)中查到已收錄的實(shí)驗(yàn)數(shù)據(jù)支持的靶基因有SOX6、MTM1、TAB3。2.3 hsa-miRNA-208a靶基因GO分析和信號(hào)通路分析結(jié)果 將預(yù)測(cè)的靶基因和已證實(shí)的靶基因進(jìn)行合并后進(jìn)行GO功能富集分析,結(jié)果顯示這些靶基因富集在心室肌組織發(fā)育、心室形態(tài)發(fā)育、心腔發(fā)育、心肌發(fā)育、心臟發(fā)育、細(xì)胞發(fā)育過(guò)程以及細(xì)胞發(fā)育等生物學(xué)過(guò)程和功能上(P<0.05),見(jiàn)表2。通過(guò)DAVID數(shù)據(jù)庫(kù)對(duì)hsa-miRNA-208a靶基因集進(jìn)行基于KEGG的信號(hào)通路富集分析,結(jié)果顯示hsa-miRNA-208a靶基因富集在黏附連接、緊密連接、Wnt信號(hào)通路、心肌收縮通路等通路上,以及病毒性心肌炎、肥厚型心肌病和擴(kuò)張型心肌病的疾病通路中(P<0.05),見(jiàn)表3。

        表2 hsa-miRNA-208a靶基因的GO富集分析結(jié)果

        表3 hsa-miRNA-208a靶基因的Pathway分析結(jié)果

        3 討 論

        自miRNAs的生物學(xué)活動(dòng)在哺乳動(dòng)物細(xì)胞中被描述以來(lái),現(xiàn)已發(fā)現(xiàn)的miRNA超過(guò)2 500個(gè),其中仍有1 000余個(gè)miRNAs的調(diào)控功能未被定義[10]。miRNA調(diào)控作用遍及生命體的各種活動(dòng)中,對(duì)miRNA作用機(jī)制研究主要是針對(duì)其與靶基因的相互作用。盡管比較基因組學(xué)和高通量實(shí)驗(yàn)研究能預(yù)測(cè)一個(gè)miRNA在轉(zhuǎn)錄組的數(shù)百個(gè)結(jié)合位點(diǎn),但樣本穩(wěn)定性和可重復(fù)性差、實(shí)驗(yàn)價(jià)格昂貴,且miRNA的表達(dá)受干擾后能導(dǎo)致多個(gè)預(yù)測(cè)靶基因和下游蛋白的微小變化,靶基因參與的調(diào)控機(jī)制也比較復(fù)雜,所以通過(guò)生物信息學(xué)對(duì)miRNA進(jìn)行分析,預(yù)測(cè)其靶基因,并對(duì)其靶基因進(jìn)行生物功能富集分析和通路富集分析,對(duì)確定miRNA的研究方向并指導(dǎo)生物實(shí)驗(yàn)尤為重要。

        miRNA-208a是肌球蛋白基因編碼的miRNAs家族中的一員,于2003年首次被發(fā)現(xiàn),當(dāng)時(shí)命名其為miRNA-208[11]。在2009年,發(fā)現(xiàn)miRNA-208b后miRNA-208重新命名為miRNA-208a。miRNA-208a在心臟細(xì)胞特異性高表達(dá),由Myh6基因的29號(hào)內(nèi)含子編碼。其生理功能是協(xié)同miRNA-208b、miRNA-499參與控制心肌肌球蛋白含量、肌纖維特性和肌性能表現(xiàn)[12]。

        本研究通過(guò)對(duì)多個(gè)物種間miRNA-208a成熟序列進(jìn)行比對(duì),發(fā)現(xiàn)多個(gè)種屬間該序列高度保守,表明它可完成許多生命形式所必須的基本功能。同時(shí),我們采用miRanda、miRDB、TargetScan和miRwalk軟件進(jìn)行靶基因預(yù)測(cè),發(fā)現(xiàn)靶基因數(shù)目龐大,因此應(yīng)用miRwalk軟件篩選各預(yù)測(cè)軟件的交集部分,并合并DIANA Lab TarBase數(shù)據(jù)庫(kù)中已驗(yàn)證的靶基因,進(jìn)行GO富集分析和基于KEGG的信號(hào)通路分析。GO富集分析顯示hsa-miRNA-208a靶基因富集于心室肌組織發(fā)育、心室形態(tài)發(fā)育、心腔發(fā)育、心肌發(fā)育、心臟發(fā)育、細(xì)胞發(fā)育過(guò)程以及細(xì)胞發(fā)育等生物學(xué)過(guò)程和功能上。信號(hào)通路分析結(jié)果表明hsa-miRNA-208a存在黏附連接、緊密連接、心肌收縮、Wnt信號(hào)通路等通路上,以及病毒性心肌炎、肥厚型心肌病和擴(kuò)張型心肌病的疾病通路中,hsa-miRNA-208a示其與心肌疾病的發(fā)生密切相關(guān)。目前其他學(xué)者也發(fā)現(xiàn)miRNA-208a表達(dá)異常主要是與各種心肌疾病的發(fā)生有關(guān)。其中,van Rooij等[13]研究發(fā)現(xiàn)miRNA-208a與心肌肥大、纖維化有關(guān),且miRNA-208a通過(guò)調(diào)節(jié)靶蛋白—甲狀腺激素受體相關(guān)蛋白參與壓力負(fù)荷下心肌肥厚與纖維化。Callis等[14]通過(guò)研究再次證明miRNA-208a與心肌肥厚相關(guān),且發(fā)現(xiàn)miRNA-208a過(guò)表達(dá)的轉(zhuǎn)基因鼠出現(xiàn)房室傳導(dǎo)阻滯,而敲除miRNA-208a基因的小鼠出現(xiàn)心房顫動(dòng),提示miRNA-208a是心臟正常傳導(dǎo)所必需的,并與其調(diào)控唯同源異型結(jié)構(gòu)域蛋白、鋅指蛋白轉(zhuǎn)錄因子4和縫隙連接蛋白40有關(guān)。有學(xué)者通過(guò)尸檢發(fā)現(xiàn),與正常成人和胎兒相比,心肌梗死者心肌組織miRNA-208a表達(dá)上調(diào),在心肌梗死者和胎兒心臟中miRNA的表達(dá)模式相似,提示在心臟重塑中心臟遺傳因子發(fā)生了重編程[15]。另一項(xiàng)關(guān)于擴(kuò)張型心肌病的研究結(jié)果顯示,人心肌組織中miRNA-208a表達(dá)量與β-MHC的mRNA表達(dá)量以及心肌膠原含量有關(guān),miRNA-208a表達(dá)升高提示預(yù)后不佳[16]。

        綜上所述,hsa-miRNA-208a參與心臟生長(zhǎng)發(fā)育的生物學(xué)過(guò)程,與心肌疾病的發(fā)生密切相關(guān)。我們的研究結(jié)果提示,hsa-miRNA-208a除了與心肌病與缺血性心肌病變發(fā)生有關(guān),可能在先天性心臟病和各種病理因素導(dǎo)致的心肌細(xì)胞凋亡中也起著重要的作用。此外,我們的研究還為進(jìn)一步探索心肌病發(fā)生的下游信號(hào)通路(如MAPK/NKL 信號(hào)通路、TGF-β/Smad信號(hào)通路、Ca/Myosin信號(hào)通路)指出了潛在的方向。這對(duì)于我們進(jìn)一步研究miRNA-208a在心血管疾病發(fā)生中的作用具有一定的指導(dǎo)意義。

        [1] 夏 偉,曹?chē)?guó)軍,邵寧生.MicroRNA靶基因的尋找及鑒定方法研究進(jìn)展[J].中國(guó)科學(xué)C輯:生命科學(xué),2009,39(1):121-128.

        [2] 張春霆.生物信息學(xué)的現(xiàn)狀與展望[J].世界科技研究與發(fā)展,2000,22(6):17-20.

        [3] 方德福.microRNA的研究方法與應(yīng)用[M].北京:中國(guó)協(xié)和醫(yī)科大學(xué)出版社,2008:34.

        [4] Shyu KG,Wang BW,Cheng WP,et al.MicroRNA-208a increases myocardial endoglin expression and myocardial fibrosis in acute myocardial infarction[J].Can J Cardiol,2015,31(5):679-690.

        [5] Slagsvold KH,Johnsen AB,Rognmo O,et al.Mitochondrial respiration and microRNA expression in right and left atrium of patients with atrial fibrillation[J].Physiol Genomics,2014,46(14):505-511.

        [6] Ooi JY,Bernardo BC,Mcmullen JR.The therapeutic potential of miRNAs regulated in settings of physiological cardiac hypertrophy[J].Future Med Chem,2014,6(2):205-222.

        [7] Hasan MM,Akter R,Ullah MS,et al.A computational approach for predicting role of human microRNAs in MERS-CoV genome[J].Adv Bioinformatics,2014,2014:967 946.

        [8] Ferreira LR,Frade AF,Santos RH,et al.MicroRNAs miR-1,miR-133a,miR-133b,miR-208a and miR-208b are dysregulated in Chronic Chagas disease Cardiomyopathy[J].Int J Cardiol,2014,175(3):409-417.

        [9] Greco S,Perfetti A,Fasanaro P,et al.Deregulated microRNAs in myotonic dystrophy type 2[J].PLoS One,2012,7(6):e39732.

        [10]Min PK,Chan SY.The biology of circulating microRNAs in cardiovascular disease[J].Eur J Clin Invest,2015,45(8):860-874.

        [11]Lagos-Quintana M,Rauhut R,Meyer J,et al.New microRNAs from mouse and human[J].RNA,2003,9(2):175-179.

        [12]Van Rooij E,Quiat D,Johnson BA,et al.A family of microRNAs encoded by myosin genes governs myosin expression and muscle performance[J].Dev Cell,2009,17(5):662-673.

        [13]van Rooij E,Sutherland LB,Qi X,et al.Control of stress-dependent cardiac growth and gene expression by a microRNA[J].Science,2007,316(5824):575-579.

        [14]Callis TE,Pandya K,Seok HY,et al.MicroRNA-208a is a regulator of cardiac hypertrophy and conduction in mice[J].J Clin Invest,2009,119(9):2 772-2 786.

        [15]Bostjancic E,Zidar N,Stajer D,et al.MicroRNAs miR-1,miR-133a,miR-133b and miR-208 are dysregulated in human myocardial infarction[J].Cardiology,2010,115(3):163-169.

        [16]Satoh M,Minami Y,Takahashi Y,et al.Expression of microRNA-208 is associated with adverse clinical outcomes in human dilated cardiomyopathy[J].J Card Fail,2010,16(5):404-410.

        Prediction of target gene of homo sapiens-microRNA-208a and bioinformatics analysis of its biological processes and signaling pathways

        LIShan-shan,ZHONGGuo-qiang,JIANGZhi-yuan,WENLi-na

        (DepartmentofCardiovascularMedicine,theFirstAffiliatedHospitalofGuangxiMedicalUniversity,Nanning530021,China)

        Objective To predict the target genes of homo sapiens-microRNA-208a(hsa-miRNA-208a),and analyze their possible biological processes and signaling pathways using bioinformatics technology.Methods Basic data of hsa-miRNA-208a including chromosomal location,base sequence and conservation of species were obtained using miRbase database and UCSC Genome Browser tool.The target genes were predicted using miRanda,miRDB,TargetScan and miRwalk.An intersection of target genes was acquired using miRwalk,and then was combined with the target genes supported by literature and confirmed by experiments as a target gene set,which were analyzed with Gene Ontology(GO) analysis and signal pathway analysis for target gene set.Results The sequence of has-miRNA-208a was highly conserved among various species.Sixteen target genes(including CPEB2,CSNK2A2,DLD,MTF2,FBXO28,LEP,SMAD4,NLK,GPR88,VPS13D,ZNF215,CHD9,SLC7A5,C19orf44,SLC45A3 and MAP4K4) were obtained from the intersection of target genes predicted by TargetScan,miRDB,miRanda and miRwalk.Twelve(CASP3,MRGPRX3,FPRL1,MED13,IL10,GATA4,HSH2D,MYH6,MYH7,TNNI3,CDKN1A and PAK3) and three(SOX6,MTM1 and TAB3) confirmed target genes were retrieved from miRwalk and DIANA Lab TarBase,respectively.GO analysis showed that the target genes were enriched in biological process including ventricular cardiac muscle tissue development,cardiac ventricle morphogenesis,cardiac chamber development,cardiac muscle tissue development,heart development,cell development,ect(P<0.05).The signal pathway analysis showed that the target genes were enriched in the pathways including adherent junction,tight junction,cardiac muscle contraction,Wnt signaling pathway,viral myocarditis,hypertrophic cardiomyopathy,dilated cardiomyopathy,ect(P<0.05).Conclusion Has-miRNA-208a participates in the biological process of cardiac growth and development,and is closely related to the occurrence of myocardial disease. 【Key words】 MicroRNA-208a,Homo sapiens,Bioinformatics,Target gene,Heart,Growth,Myocardial disease

        黎珊珊(1990~),女,在讀碩士研究生,研究方向:心電生理學(xué)、心律失常及冠心病的基礎(chǔ)研究及介入治療。通信作者:鐘國(guó)強(qiáng)(1958~),男,博士,教授,研究方向:心電生理學(xué)、心律失常及冠心病的基礎(chǔ)研究及介入治療,E-mail:gq_zhong@126.com。

        R 342.7

        A

        0253-4304(2016)12-1634-04

        10.11675/j.issn.0253-4304.2016.12.02

        2016-05-28

        2016-08-08)

        猜你喜歡
        信息學(xué)心肌病生物學(xué)
        谷稗的生物學(xué)特性和栽培技術(shù)
        雞NRF1基因啟動(dòng)子區(qū)生物信息學(xué)分析
        初中生物學(xué)糾錯(cuò)本的建立與使用
        初中生物學(xué)糾錯(cuò)本的建立與使用
        伴有心肌MRI延遲強(qiáng)化的應(yīng)激性心肌病1例
        初論博物館信息學(xué)的形成
        擴(kuò)張型心肌病中醫(yī)辨治體會(huì)
        PEDF抗腫瘤的生物學(xué)作用
        TAKO-TSUBO心肌病研究進(jìn)展
        miRNA-148a在膀胱癌組織中的表達(dá)及生物信息學(xué)分析
        国产精品日韩av一区二区三区| av一区无码不卡毛片| 官网A级毛片| 国产精品视频一区国模私拍| 亚洲成人av一区二区三区| 亚洲专区一区二区三区四区五区| 欧美国产激情18| 成人区人妻精品一熟女 | 中文字幕亚洲好看有码| 久久精品国产亚洲av夜夜| 精品亚洲成av人在线观看| 无套内射无矿码免费看黄| 亚洲VA不卡一区| 亚洲av中文字字幕乱码软件| 久久不见久久见www日本网| 国产精品无码a∨精品影院| 久久99热精品免费观看欧美| 国产亚洲av一线观看| 很黄很色很污18禁免费| 国产自偷亚洲精品页65页| 9丨精品国产高清自在线看| 精品精品国产三级av在线| 欧美人与禽2o2o性论交| 亚洲av无码av在线播放| 农村国产毛片一区二区三区女| 国产精品一区二区av麻豆日韩| 老师粉嫩小泬喷水视频90| 国产男女插插一级| 国产精女同一区二区三区久| 99麻豆久久久国产精品免费| 久久人人97超碰超国产| 亚洲最大成av人网站| 国内激情一区二区视频| 亚洲乱码无人区卡1卡2卡3| 蜜臀久久99精品久久久久久小说| 欧美精品高清在线xxxx| 日韩美女av一区二区| 亚洲avav天堂av在线网毛片| 小12箩利洗澡无码视频网站| 亚洲粉嫩视频在线观看| 日本大乳高潮视频在线观看|