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

    膽道閉鎖潛在分子亞型及其核心基因識別

    2022-06-20 01:30:40陳治宏羅鳳媛馬麗媛李紅東
    贛南醫(yī)學院學報 2022年4期
    關(guān)鍵詞:秩次差異基因亞型

    陳治宏,羅鳳媛,馬麗媛,李紅東

    (1. 贛南醫(yī)學院2019級碩士研究生;2. 贛南醫(yī)學院2020級碩士研究生;3. 贛南醫(yī)學院2021級碩士研究生;4. 贛南醫(yī)學院信息工程學院,江西 贛州 341000)

    膽道閉鎖(Biliary atresia,BA)是新生兒常見的一種進行性阻塞性膽道疾病,臨床表現(xiàn)為梗阻性黃疸,如果BA 患兒得不到及時治療,易發(fā)展為膽汁淤積性肝硬化、門靜脈高壓,最終導(dǎo)致肝衰竭而死亡[1]。膽道閉鎖具有地區(qū)和種族差異,一般認為亞洲人發(fā)病率較高,尤以日本和中國的發(fā)病率最高[2]。目前對BA 患兒應(yīng)用最廣泛且療效較好的治療手段是Kasai 手術(shù),但由于進行性肝纖維化、膽汁淤積性肝硬化和門脈高壓的出現(xiàn),多數(shù)患兒最終須選擇肝臟移植救治[3]。BA 因其早期診斷困難、病理改變特殊以及治療效果不理想,已成為威脅新生兒生命健康的重要危險因素。

    目前,研究者提出了許多理論來解釋BA的病因,包括遺傳變異、病毒感染和免疫介導(dǎo)等。CHENG G等[4]對89 例BA 患者開展全基因組拷貝數(shù)突變分析發(fā)現(xiàn)了29 個BA 相關(guān)的拷貝數(shù)變異,其中包括富集在炎癥相關(guān)通路的基因,提示了炎癥在BA 中的重要作用。SHEN C 等[5]研究表明,60%的BA 患兒在診斷時檢測到巨細胞病毒(CMV)的DNA,同時感染CMV 的BA 患兒患有膽管炎和較高的膽管纖維化程度。SQUIRES J E 等[6]通過對輪狀病毒(RV)誘導(dǎo)的BA小鼠模型研究發(fā)現(xiàn),BA小鼠的NK細胞數(shù)量增加會促進慢性膽管炎癥。SAXENA V 等[7]發(fā)現(xiàn)RV 誘導(dǎo)的BA 小鼠模型中的漿細胞樣樹突狀細胞(pDCs)顯著增加,同時pDCs產(chǎn)生白細胞介素15(IL-15)激活CD80,加速膽管纖維化進程。BA 的復(fù)雜致病因素,提示BA中可能存在不同的分子亞型。

    目前,基于高通量數(shù)據(jù)對疾病進行分子亞型的研究主要是根據(jù)樣本基因表達水平的相似性。通過無監(jiān)督聚類的方法,例如非負矩陣分解(NMF)、層次聚類和潛在因子分析(LF)等挖掘潛在的分子亞型[8-9]。但這類方法往往基于較為復(fù)雜的數(shù)量統(tǒng)計模型,難以理解和進行生物學解釋,且不易進行臨床實踐和驗證。有研究表明,基因間相對表達秩次(大小)關(guān)系在正常組織樣本內(nèi)存在廣泛的穩(wěn)定性,在疾病組織樣本中受到了廣泛的擾動?;谶@一現(xiàn)象,本文提出了一種基于樣本內(nèi)基因間相對表達秩次關(guān)系的方法識別潛在疾病分子的亞型方法。這種方法不僅易于理解和解釋,還可整合來自不同實驗室的樣本進行生物信息學分析,對基因表達水平的個體間生物變異具有魯棒性[10-11]。本研究將該方法應(yīng)用于挖掘BA 潛在的分子亞型,進而識別亞型相關(guān)的核心基因。

    1 材料與方法

    1.1 數(shù)據(jù)來源及其預(yù)處理本文所分析的基因表達譜數(shù)據(jù)集,GSE46960 和GSE15235 均來自基因表達綜合數(shù)據(jù)庫(GEO,https://www.ncbi.nlm.nih.gov/geo/)。GSE46960數(shù)據(jù)集檢測平臺為GPL6244,共檢測了64 例BA 患兒的肝組織、7例正常嬰兒的肝組織。GSE15235 數(shù)據(jù)集檢測平臺為GPL570,共檢測了43例BA患兒的肝組織。

    各數(shù)據(jù)集均從GEO 下載原始表達譜數(shù)據(jù)(. CEL 文件),利用R 軟件Affy 軟件包提供的RMA(Robust Multichip Analysis)算法對其進行背景校正[12],并利用相應(yīng)的平臺注釋文件將探針I(yè)D 對應(yīng)到基因ID 上。如果多個探針映射到同一基因,則取探針的平均表達值作為該基因的表達水平。最終,數(shù)據(jù)集GSE46960 共檢測18 859 個基因,數(shù)據(jù)集GSE15235 共檢測20 684 個基因。兩套數(shù)據(jù)集共同檢測的17 179個基因作為背景基因用作后續(xù)分析。

    1.2 基于基因間相對表達秩次關(guān)系識別BA 分子亞型假設(shè)樣本內(nèi)基因間相對表達秩次關(guān)系與疾病的分子亞型相關(guān),那么應(yīng)存在以下幾種現(xiàn)象:⑴在不同分子亞型樣本內(nèi),存在基因間相對表達秩次關(guān)系不同的現(xiàn)象。⑵具有潛在鑒別疾病分子亞型能力的基因?qū)Γ雌湎鄬Ρ磉_秩次關(guān)系劃分的兩組樣本間存在差異表達的基因,同時涉及生物學功能的改變。⑶識別具有鑒別疾病分子亞型能力的潛在基因?qū)Υ嬖诰奂F(xiàn)象。

    基于上述假設(shè),本文識別疾病分子亞型的具體步驟如下。

    1.2.1 篩選正常組織樣本中基因間相對表達秩次關(guān)系穩(wěn)定的基因?qū)⒈尘盎騼蓛山M合,如果一對基因(基因i 和基因j)在正常樣本內(nèi)滿足P(Gi>Gj)≥0.9[公式⑴],那么該基因?qū)Χx為穩(wěn)定基因?qū)Α?/p>

    其中n表示所有的正常樣本數(shù),t代表第t個樣本,Gi>Gj表示基因i和基因j的表達秩次關(guān)系。

    1.2.2 基于穩(wěn)定基因?qū)Y選疾病組織樣本中的顯著逆轉(zhuǎn)基因?qū)谡颖局械玫降姆€(wěn)定基因?qū)?,評估其在疾病中相對表達秩次關(guān)系改變情況。對于任一穩(wěn)定對(Gi,Gj),假設(shè)在正常樣本中觀察到Gi>Gj的樣本有n1個、Gi≤Gj的樣本有n2個,在疾病樣本中滿足Gi>Gj的樣本有m1個、Gi≤Gj的樣本有m2個。利用Fisher 精確檢驗,可檢驗Gi、Gj間相對表達秩次關(guān)系在正常和疾病樣本間分布是否存在差異。經(jīng)多重假設(shè)檢驗校正后,滿足假陽性發(fā)現(xiàn)率(FDR)<5%的基因?qū)Χx為逆轉(zhuǎn)基因?qū)Α?/p>

    1.2.3 篩選具有潛在疾病分子亞型鑒別能力的基因?qū)τ谠诩膊≈泻Y選得到的逆轉(zhuǎn)基因?qū)?,按如下步驟進一步篩選具有潛在疾病分子亞型鑒別能力的基因?qū)Γ孩拍孓D(zhuǎn)對的秩次關(guān)系的改變只在20%~80%的疾病患者之中發(fā)生。這是為了確?;?qū)哂需b別疾病分子亞型的能力。⑵對于任意逆轉(zhuǎn)對(Gi,Gj),根據(jù)其相對表達秩次關(guān)系,可將疾病樣本分為兩組,Gi>Gj為1 組,Gi≤Gj為1 組。然后利用t檢驗,控制FDR 為5%,識別兩組間的差異基因,并基于KEGG(京都基因與基因組百科全書,https://www. kegg. jp/)通路數(shù)據(jù)庫對兩組間識別的差異基因進行功能富集分析。若一個生物學通路顯著富集了差異基因,則認為該通路發(fā)生了生物學功能的擾動。對于一個逆轉(zhuǎn)基因?qū)?,可按其在疾病樣本中的相對表達秩次關(guān)系將疾病樣本分成兩組。如果這兩組樣本間存在差異表達基因,同時存在生物學功能的擾動,則認為該基因?qū)哂需b別潛在疾病分子亞型的能力,稱其為分子亞型候選基因?qū)Α?/p>

    1.2.4 分子亞型識別首先,按照分子亞型候選基因?qū)_動通路的數(shù)目,將基因?qū)M行降序排列。然后,構(gòu)建一個分子亞型候選基因?qū)Α良膊颖镜木仃?,其中元素rtij為矩陣中第t個樣本中基因?qū)Γ╥,j)的相對秩次關(guān)系,其值為1 或0(如果Gi>Gj,則為1,如果Gi≤Gj,則為0)。以1 000 為梯度,分別取前k個基因?qū)Γ╧=1 000*n,n=1,2,3…)對應(yīng)的矩陣,利用歐式距離進行聚類分析,識別疾病分子亞型。為了保證結(jié)果的可靠性,對k個基因?qū)垲惤Y(jié)果,利用隨機實驗進行驗證:隨機從所有分子亞型候選基因?qū)χ刑暨xk 個基因?qū)M行聚類分析,重復(fù)100 次。如果隨機分類結(jié)果與前k個基因?qū)Φ姆诸惤Y(jié)果一致率<95%,則重新選擇標記。反之,將其作為分子亞型預(yù)測標記,用作后續(xù)分析。

    1.3 基于加權(quán)基因共表達網(wǎng)絡(luò)識別疾病分子亞型相關(guān)基因模塊基于加權(quán)基因共表達網(wǎng)絡(luò)(WGCNA)分析提取與疾病亞型顯著相關(guān)的模塊(用R 包WGCNA 實現(xiàn))[13]。首先,對背景基因表達矩陣相關(guān)系數(shù)進行加權(quán),使基因間的相互作用關(guān)系符合無標度分布。然后對基因進行分類,并將具有相似表達模式的基因分成一個模塊,最小模塊大?。╩in-ModuleSize)設(shè)置為100,其他參數(shù)設(shè)置為默認值。同一模塊的基因往往表現(xiàn)出相似的表達模式和功能[14]。

    1.4 模塊基因功能富集分析基于京都基因和基因組百科全書(KEGG)數(shù)據(jù)庫對模塊基因進行通路富集分析,用R包clusterProfile實現(xiàn)[15]。

    1.5 PPI 網(wǎng)絡(luò)構(gòu)建和核心基因篩選基于蛋白質(zhì)互作在線數(shù)據(jù)庫STRING(www. string-db. org)對候選模塊內(nèi)的基因進行蛋白與蛋白之間的互作分析(PPI),利用網(wǎng)絡(luò)最大集團度(MCC)評分算法[公式⑵][16],識別PPI網(wǎng)絡(luò)中的關(guān)鍵基因。

    其中,集合S表示集合中元素的數(shù)量,S(ν)是包含ν 的最大集團的集合,(|C|-1)!是所有小于|C|的正整數(shù)的乘積。

    2 結(jié) 果

    2.1 具有潛在疾病分子亞型鑒別能力的基因?qū)Π凑辗椒▽W篩選具有潛在疾病分子亞型鑒別能力的基因?qū)α鞒??;跀?shù)據(jù)集GSE46960,在每個樣本中分別對基因的表達水平進行兩兩比較,形成基因?qū)Α翗颖镜?-1 矩陣X。Xij表示第i 個基因?qū)Γ℅a,Gb)在第j 個樣本中的相對表達秩次關(guān)系,其值為1或0(1代表Ga>Gb,0代表Ga≤Gb)。在7例正常嬰兒樣本中,有144 822 185個基因?qū)υ?0%以上的正常樣本組織中具有Ga>Gb的表達模式,其中有5 105 085個基因?qū)Φ南鄬Ρ磉_秩次關(guān)系在BA 樣本中發(fā)生顯著逆轉(zhuǎn)Ga≤Gb(Fisher 精確檢驗,F(xiàn)DR<0.05),暗示了這些基因?qū)Φ南鄬Ρ磉_秩次關(guān)系可能與BA 的發(fā)生發(fā)展相關(guān)。對每一個逆轉(zhuǎn)基因?qū)?,進一步根據(jù)其在疾病樣本中的相對表達秩次關(guān)系,將疾病樣本分為兩組(Ga>Gb為一組,Ga≤Gb為一組),然后識別兩組間的差異基因并基于KEGG數(shù)據(jù)庫進行功能富集分析,發(fā)現(xiàn)有4 836 202 個基因?qū)M足條件,用作后續(xù)分析。

    2.2 識別BA 的分子亞型基于亞型候選基因?qū)?,首先?gòu)建亞型候選基因?qū)Α良膊颖镜?-1 矩陣X。矩陣中xij表示第i 個基因?qū)υ诘趈 個樣本中的相對表達秩次關(guān)系,其值為1 或0(1 代表>,0 代表≤)。按照亞型候選基因?qū)Ω患降纳飳W功能通路個數(shù)由大到小進行排序,然后分別選取前1 000 對、2 000 對、3 000 對、4 000 對和5 000 對基因?qū)Γ瑢?4例BA 樣本分別進行聚類分析。結(jié)果顯示前3 000 2A),采用動態(tài)剪切法劃分模塊,合并相似度>75%的模塊,最終構(gòu)建了17 個共表達模塊(圖2B)。進一步分析了BA 不同亞型與各種基因模塊的相關(guān)性,發(fā)現(xiàn)洋紅色(Magenta)模塊與BA 亞型最相關(guān)(|r|=0.58,P=6×10-7)(圖2C)。

    圖2 BA亞型基因共表達網(wǎng)絡(luò)模塊的構(gòu)建

    基于KEGG數(shù)據(jù)庫,對Magenta模塊中涉及的基因進行富集分析,發(fā)現(xiàn)這些基因主要涉及PI3K-Akt信號通路、黏附斑激酶通路(FAK)和ECM 受體相互作用通路等(圖3A)。對、4 000 對、5 000 對的聚類結(jié)果與隨機情況一致。因此,本文將前3 000 對基因?qū)Φ木垲惤Y(jié)果作為BA分子亞型識別依據(jù)。聚類結(jié)果顯示64例BA 樣本明顯聚集為兩類:一類樣本個數(shù)為54 個,另一類樣本個數(shù)為10個(圖1A)。

    為了進一步驗證BA 中存在的不同分子亞型,我們把得到的前3 000 個亞型候選基因?qū)?yīng)用于GSE15235 數(shù) 據(jù)集中,GSE15235 數(shù)據(jù)集中43 例BA樣本中依然明顯聚集為兩大類:一類樣本個數(shù)15個,另一類樣本個數(shù)為28 個(圖1B)。通過數(shù)據(jù)集GSE46960、GSE15235 間兩類樣本間的差異基因比較發(fā)現(xiàn),在GSE46960 兩類間識別了436 個差異基因,GSE15235 兩類間識別了584 個差異基因,共交疊了75個差異基因(超幾何分布P=4.8×10-32),并且這75 個差異基因有69 個基因在兩組間失調(diào)方向一致(二項分布P=5.77×10-15)。

    圖1 BA樣本分子亞型分類

    2.3 BA 疾病亞型共表達模塊構(gòu)建及核心基因識別為了進一步分析BA 中不同亞型的基因表達模式,利用WGCNA 方法對數(shù)據(jù)集GSE46960構(gòu)建加權(quán)共表達網(wǎng)絡(luò)。使用無標度拓撲標準選擇β=8(圖

    將所有Magenta 模塊內(nèi)的基因?qū)隨TRING 數(shù)據(jù)庫以構(gòu)建PPI 網(wǎng)絡(luò),通過MCC 算法發(fā)現(xiàn)了10 個BA 亞型相關(guān)的核心基因(LUM、COL6A3、FBN1、SPARC、DCN、LAMA4、FAP、ANTXR1、LAMA2 和COL1A2)。LAMA2、LAMA4 和LUM 基因(圖中紅色基因)與其他蛋白質(zhì)的相互作用更頻繁(圖3B)。

    圖3 Magenta模塊中基因的KEGG通路富集分析和PPI網(wǎng)絡(luò)構(gòu)建

    對10個核心基因在GSE46960數(shù)據(jù)集的正常組織、BA亞型Ⅰ和亞型Ⅱ三組間的表達情況進行分析發(fā)現(xiàn),相對于正常組織,10 個基因在亞型Ⅰ中顯著高表達,而在亞型Ⅱ中顯著低表達(圖4A)。在GSE15235 數(shù)據(jù)集中,我們還發(fā)現(xiàn)這10 個核心基因的表達水平在纖維化和炎癥兩類樣本間有顯著差異,且纖維化組的表達水平顯著高于炎癥組(圖4B),與BA亞型Ⅰ和亞型Ⅱ一致。

    圖4 10個核心基因在不同亞型間的表達分布情況

    3 討 論

    本研究基于樣本內(nèi)基因間相對表達秩次關(guān)系,開發(fā)了疾病分子亞型識別算法,確定了在BA 疾病中存在兩組不同的疾病亞型,并利用WGCNA 分析篩選出了與亞型相關(guān)共表達基因模塊。通路富集分析顯示,Magenta 模塊內(nèi)基因主要涉及PI3K-Akt信號通路、黏附斑激酶通路(FAK)和ECM 受體相互作用通路等。有研究表明,PI3K-Akt 信號通路可能抑制NF-κB 和NLRP3 炎癥通路的蛋白表達,從而減少炎癥因子的分泌[17]。FAK 和ECM 受體相互作用通路抑制TGF-β1信號傳導(dǎo),降低成纖維細胞分化[18]。最后,本文還識別了與BA 分子亞型相關(guān)的10 個核心基因。其中,LAMA2 和LAMA4 可促進轉(zhuǎn)化生長因子(TGF-β)信號傳導(dǎo),而TGF-β 信號會阻止慢性損傷期間促纖維化FAP 的細胞凋亡,從而促進纖維化環(huán)境[19-20]。KRISHNAN A 等[21]研究表明,LUM 可在體外誘導(dǎo)生長因子β1(TGF-β1)轉(zhuǎn)化為促纖維化細胞因子。COL1A2和COL6A3編碼膠原蛋白,而膠原蛋白的過度表達是肝纖維化的基本特征。當肝星形細胞(HSCs)被激活時,基質(zhì)金屬蛋白(MMPs)的表達增加,同時MMPs 的抑制劑基質(zhì)金屬蛋白酶(TIMPs)的表達增加。如果TIMPs 增加過快,MMPs/TIMPs 的比值會發(fā)生變化,使細胞外基質(zhì)(ECM)的合成和降解失衡,促進纖維化的發(fā)展[22-23]。其余5個關(guān)鍵基因(FBN1、SPARC、DCN、FAP 和ANTXR1)雖然還未有其與BA 相關(guān)報道,但其參與了蛋白質(zhì)消化吸收、ECM-受體相互作用和PI3K-Akt 信號通路。其中ECM 的變化由多種介質(zhì)和生長因子誘導(dǎo),從而調(diào)節(jié)各種效應(yīng),如刺激血管生成、炎癥反應(yīng)和促進基質(zhì)侵襲,可導(dǎo)致成纖維細胞的遷移或分化,最終使細胞纖維化[24-25]。提示這10 個核心基因與BA 進行性肝炎、肝纖維化有關(guān)。

    在本研究中,正常樣本相對較少,在穩(wěn)定基因?qū)Y選過程中,可能引入一定的假陽性基因?qū)?,后續(xù)需要更多的正常樣本來驗證。同時,本研究篩選的亞型候選基因?qū)?shù)量較多,不利于后續(xù)聚類鑒別疾病亞型。因此,如何篩選出更少、更具有生物學意義的亞型候選基因?qū)κ俏磥淼难芯恐攸c。后續(xù)研究我們將考慮如基因在候選對中出現(xiàn)的頻率、基因間的表達相關(guān)性等因素進一步改進算法。

    猜你喜歡
    秩次差異基因亞型
    ICR鼠肝和腎毒性損傷生物標志物的篩選
    基于RNA 測序研究人參二醇對大鼠心血管內(nèi)皮細胞基因表達的影響 (正文見第26 頁)
    非整秩次秩和比法綜合評價醫(yī)院工作效率
    Ikaros的3種亞型對人卵巢癌SKOV3細胞增殖的影響
    山東省農(nóng)業(yè)可持續(xù)發(fā)展能力評價研究——基于非整秩次WRSR
    ABO亞型Bel06的分子生物學鑒定
    SSH技術(shù)在絲狀真菌功能基因篩選中的應(yīng)用
    HeLa細胞中Zwint-1選擇剪接亞型v7的表達鑒定
    抑郁的新亞型:閾下抑郁
    2009—2010賽季CBA各隊攻防能力的分析
    久久久久精品久久久久真实原创| 一级毛片久久久久久久久女| 久久6这里有精品| 久久99热这里只频精品6学生| 青春草亚洲视频在线观看| 精品久久久久久久久av| 天天躁日日操中文字幕| 五月天丁香电影| 久热久热在线精品观看| 男女无遮挡免费网站观看| 高清视频免费观看一区二区| 午夜福利网站1000一区二区三区| 亚洲欧美精品专区久久| 十分钟在线观看高清视频www | 纵有疾风起免费观看全集完整版| 3wmmmm亚洲av在线观看| 一本—道久久a久久精品蜜桃钙片| 成人影院久久| 春色校园在线视频观看| 久久久久久人妻| 99热这里只有精品一区| 精品一品国产午夜福利视频| 成人亚洲欧美一区二区av| 日韩一本色道免费dvd| 久久久成人免费电影| 久久久精品94久久精品| 高清欧美精品videossex| 麻豆精品久久久久久蜜桃| 亚洲国产精品国产精品| 人人妻人人爽人人添夜夜欢视频 | 亚洲第一av免费看| 日日摸夜夜添夜夜添av毛片| 三级国产精品欧美在线观看| 久久久久久久久大av| 国产精品爽爽va在线观看网站| 亚洲婷婷狠狠爱综合网| 男人狂女人下面高潮的视频| 少妇高潮的动态图| 少妇人妻久久综合中文| 夜夜骑夜夜射夜夜干| videos熟女内射| 欧美最新免费一区二区三区| 我的女老师完整版在线观看| 2021少妇久久久久久久久久久| 久久99热6这里只有精品| 中文字幕人妻熟人妻熟丝袜美| 国产精品麻豆人妻色哟哟久久| 日日摸夜夜添夜夜添av毛片| 这个男人来自地球电影免费观看 | 国产美女午夜福利| 中文乱码字字幕精品一区二区三区| 99热这里只有是精品在线观看| 乱码一卡2卡4卡精品| 人人妻人人爽人人添夜夜欢视频 | 青春草视频在线免费观看| 99国产精品免费福利视频| 天堂中文最新版在线下载| 国产真实伦视频高清在线观看| 18禁裸乳无遮挡动漫免费视频| 综合色丁香网| 男人爽女人下面视频在线观看| 久久久久久久大尺度免费视频| 日韩电影二区| 亚洲第一区二区三区不卡| 国产高潮美女av| 26uuu在线亚洲综合色| 国国产精品蜜臀av免费| 简卡轻食公司| 久久久久久久久久久丰满| 国产精品99久久久久久久久| 免费高清在线观看视频在线观看| 国产精品一区二区在线不卡| 欧美精品亚洲一区二区| 精品一区二区三卡| 肉色欧美久久久久久久蜜桃| 日韩精品有码人妻一区| 99九九线精品视频在线观看视频| 国产亚洲精品久久久com| 狠狠精品人妻久久久久久综合| 联通29元200g的流量卡| 欧美区成人在线视频| 亚洲欧美一区二区三区黑人 | 婷婷色av中文字幕| 蜜桃在线观看..| 全区人妻精品视频| 日本vs欧美在线观看视频 | 一级毛片久久久久久久久女| 久久精品久久久久久久性| 日韩电影二区| 国产免费一级a男人的天堂| 日本wwww免费看| 一级毛片久久久久久久久女| 国产无遮挡羞羞视频在线观看| 精品一区二区三卡| 国产免费一级a男人的天堂| 精品久久久久久电影网| 深爱激情五月婷婷| 日韩一区二区三区影片| 午夜免费观看性视频| h日本视频在线播放| 免费人妻精品一区二区三区视频| 伦理电影免费视频| 91久久精品国产一区二区成人| 欧美xxxx性猛交bbbb| 观看免费一级毛片| 成年美女黄网站色视频大全免费 | 国产欧美日韩一区二区三区在线 | 99热网站在线观看| 亚洲国产av新网站| 在线看a的网站| 国内少妇人妻偷人精品xxx网站| 一级毛片黄色毛片免费观看视频| 91狼人影院| av专区在线播放| 国产亚洲av片在线观看秒播厂| 性色av一级| 国产69精品久久久久777片| 三级国产精品片| 综合色丁香网| 视频中文字幕在线观看| 亚洲精品久久午夜乱码| 午夜激情久久久久久久| 国产男女内射视频| av国产久精品久网站免费入址| 韩国av在线不卡| 免费少妇av软件| 人体艺术视频欧美日本| 日韩制服骚丝袜av| 国产成人免费无遮挡视频| 国产精品久久久久久精品古装| 成人午夜精彩视频在线观看| 色网站视频免费| av播播在线观看一区| 少妇人妻久久综合中文| 午夜福利影视在线免费观看| 日日啪夜夜爽| 在线看a的网站| 我要看日韩黄色一级片| 亚洲国产精品国产精品| 免费高清在线观看视频在线观看| 亚洲综合精品二区| 国产精品99久久久久久久久| 亚洲av福利一区| 97精品久久久久久久久久精品| 人人妻人人爽人人添夜夜欢视频 | 久久久久久伊人网av| 热99国产精品久久久久久7| 亚洲av中文字字幕乱码综合| 成人国产麻豆网| 国产精品熟女久久久久浪| 亚洲人成网站在线观看播放| 国产精品av视频在线免费观看| 春色校园在线视频观看| 国产探花极品一区二区| 在线观看三级黄色| 一级av片app| 大陆偷拍与自拍| xxx大片免费视频| 午夜激情福利司机影院| av福利片在线观看| 成人影院久久| 一个人免费看片子| 日日啪夜夜撸| 2021少妇久久久久久久久久久| 国产精品人妻久久久久久| 日本色播在线视频| 国产淫片久久久久久久久| 在线播放无遮挡| 日本wwww免费看| 国产精品一区www在线观看| 伦精品一区二区三区| 最新中文字幕久久久久| 国产精品熟女久久久久浪| 91狼人影院| 日日撸夜夜添| 国产精品久久久久久精品古装| 人体艺术视频欧美日本| 精品少妇久久久久久888优播| 老师上课跳d突然被开到最大视频| 国产爱豆传媒在线观看| 日本爱情动作片www.在线观看| 国产在线一区二区三区精| 久久热精品热| 亚洲欧美清纯卡通| 国产成人freesex在线| 中国三级夫妇交换| 午夜精品国产一区二区电影| 久久久亚洲精品成人影院| 精品少妇久久久久久888优播| 又黄又爽又刺激的免费视频.| 亚洲一区二区三区欧美精品| 亚洲高清免费不卡视频| 伦精品一区二区三区| 最近最新中文字幕免费大全7| 免费黄色在线免费观看| 高清午夜精品一区二区三区| 一级片'在线观看视频| 一级毛片我不卡| 建设人人有责人人尽责人人享有的 | 国产精品熟女久久久久浪| 欧美日韩亚洲高清精品| 精品久久久噜噜| 黄色日韩在线| 日本午夜av视频| 人人妻人人添人人爽欧美一区卜 | 免费av中文字幕在线| 精品一区二区三区视频在线| 国产日韩欧美亚洲二区| 黄色日韩在线| 精品国产乱码久久久久久小说| 国产黄片视频在线免费观看| 国产v大片淫在线免费观看| 熟女av电影| 男女边摸边吃奶| 国产精品三级大全| 国产无遮挡羞羞视频在线观看| 亚洲国产高清在线一区二区三| 丰满乱子伦码专区| videos熟女内射| 久久99蜜桃精品久久| 午夜老司机福利剧场| 国产精品人妻久久久久久| 搡女人真爽免费视频火全软件| 我要看黄色一级片免费的| 国产伦理片在线播放av一区| 国产精品三级大全| 亚洲人成网站高清观看| 视频区图区小说| 精品人妻熟女av久视频| 亚洲成人av在线免费| 欧美另类一区| 在线观看一区二区三区激情| 男女下面进入的视频免费午夜| 国产黄色视频一区二区在线观看| 日本免费在线观看一区| 国产免费视频播放在线视频| 久久亚洲国产成人精品v| 国产成人免费观看mmmm| 欧美日韩亚洲高清精品| 街头女战士在线观看网站| 婷婷色综合www| 日韩,欧美,国产一区二区三区| 啦啦啦在线观看免费高清www| 午夜精品国产一区二区电影| 成人亚洲欧美一区二区av| 国内揄拍国产精品人妻在线| 狂野欧美激情性xxxx在线观看| 热99国产精品久久久久久7| 久久久久久人妻| 女人久久www免费人成看片| 有码 亚洲区| 亚洲美女搞黄在线观看| 自拍偷自拍亚洲精品老妇| 一本一本综合久久| 国产精品久久久久久久电影| 少妇猛男粗大的猛烈进出视频| 亚洲人成网站在线观看播放| 国产男女内射视频| 99热这里只有是精品50| 日日撸夜夜添| 韩国高清视频一区二区三区| 精品人妻偷拍中文字幕| 国内精品宾馆在线| 大陆偷拍与自拍| 久久精品久久精品一区二区三区| 伦精品一区二区三区| 久久人人爽人人爽人人片va| 日日啪夜夜爽| 九九久久精品国产亚洲av麻豆| 99热这里只有是精品在线观看| 亚洲av日韩在线播放| 天天躁日日操中文字幕| 日韩国内少妇激情av| 高清日韩中文字幕在线| 国产91av在线免费观看| 亚洲欧美日韩无卡精品| 国产69精品久久久久777片| 国产免费一级a男人的天堂| 美女内射精品一级片tv| 欧美激情极品国产一区二区三区 | 看免费成人av毛片| 18禁在线播放成人免费| 少妇的逼水好多| 国产精品国产三级国产av玫瑰| 欧美成人精品欧美一级黄| 亚洲av国产av综合av卡| 国产精品一区www在线观看| 全区人妻精品视频| 亚洲精品国产av成人精品| 女性生殖器流出的白浆| 欧美zozozo另类| 日韩在线高清观看一区二区三区| 国产精品三级大全| av专区在线播放| 久久久久视频综合| 国产 一区 欧美 日韩| 日韩伦理黄色片| 99热这里只有是精品在线观看| 亚洲av欧美aⅴ国产| 性色avwww在线观看| 久久鲁丝午夜福利片| 我要看日韩黄色一级片| 国产av一区二区精品久久 | 尤物成人国产欧美一区二区三区| 少妇人妻 视频| 啦啦啦啦在线视频资源| 欧美精品一区二区免费开放| 欧美3d第一页| 日本av免费视频播放| 51国产日韩欧美| 国产成人精品久久久久久| 六月丁香七月| 国产白丝娇喘喷水9色精品| 日韩成人av中文字幕在线观看| 成人午夜精彩视频在线观看| 欧美另类一区| 国产成人午夜福利电影在线观看| 日韩强制内射视频| 少妇的逼水好多| 高清av免费在线| 国产深夜福利视频在线观看| 久久久久精品久久久久真实原创| 国产91av在线免费观看| 国产精品av视频在线免费观看| 女的被弄到高潮叫床怎么办| 亚洲欧美精品专区久久| 国产精品久久久久久精品古装| 免费高清在线观看视频在线观看| 99热这里只有是精品在线观看| 亚洲精品自拍成人| 最近的中文字幕免费完整| 日本黄大片高清| 春色校园在线视频观看| 永久网站在线| 亚洲av免费高清在线观看| 最近中文字幕2019免费版| 卡戴珊不雅视频在线播放| 成人高潮视频无遮挡免费网站| 久久影院123| 国产成人午夜福利电影在线观看| 青春草亚洲视频在线观看| 午夜免费男女啪啪视频观看| 丰满乱子伦码专区| 91精品一卡2卡3卡4卡| 国产精品爽爽va在线观看网站| 久久这里有精品视频免费| 少妇的逼好多水| 中文字幕久久专区| xxx大片免费视频| 亚洲国产毛片av蜜桃av| 日韩中文字幕视频在线看片 | 国产精品爽爽va在线观看网站| 成人午夜精彩视频在线观看| 我要看日韩黄色一级片| 欧美激情国产日韩精品一区| 校园人妻丝袜中文字幕| 大香蕉久久网| 男男h啪啪无遮挡| www.av在线官网国产| 观看美女的网站| 日韩在线高清观看一区二区三区| 国产综合精华液| 国产精品99久久久久久久久| 亚洲,一卡二卡三卡| 一个人看视频在线观看www免费| 亚洲国产精品一区三区| 国产在线视频一区二区| 国产精品久久久久久av不卡| 欧美日韩综合久久久久久| 91精品国产国语对白视频| 久久鲁丝午夜福利片| 亚洲,一卡二卡三卡| www.av在线官网国产| 久久 成人 亚洲| 欧美激情极品国产一区二区三区 | 日韩欧美 国产精品| 国产精品久久久久久精品电影小说 | 欧美高清性xxxxhd video| 日韩中字成人| 九九爱精品视频在线观看| 在线观看一区二区三区| 亚洲国产精品国产精品| av又黄又爽大尺度在线免费看| 夜夜爽夜夜爽视频| 色哟哟·www| 日韩一区二区视频免费看| 一边亲一边摸免费视频| 日韩在线高清观看一区二区三区| 天美传媒精品一区二区| 看免费成人av毛片| 亚洲国产最新在线播放| 身体一侧抽搐| 最近的中文字幕免费完整| 亚洲精品久久久久久婷婷小说| 免费在线观看成人毛片| 日本av免费视频播放| 欧美 日韩 精品 国产| 亚洲欧美日韩无卡精品| 午夜精品国产一区二区电影| 七月丁香在线播放| 免费人成在线观看视频色| 精品久久久久久久久av| 国产成人午夜福利电影在线观看| 狂野欧美激情性xxxx在线观看| 尾随美女入室| 久久国产亚洲av麻豆专区| 性高湖久久久久久久久免费观看| 亚洲欧美日韩卡通动漫| 欧美最新免费一区二区三区| 热99国产精品久久久久久7| 亚洲精品aⅴ在线观看| 最近的中文字幕免费完整| 自拍偷自拍亚洲精品老妇| 91精品国产九色| 亚洲欧洲国产日韩| 精品亚洲成国产av| 大码成人一级视频| kizo精华| 亚洲aⅴ乱码一区二区在线播放| 夫妻性生交免费视频一级片| 国语对白做爰xxxⅹ性视频网站| 国产成人精品福利久久| 日本一二三区视频观看| 国产高潮美女av| 99热全是精品| 在线观看一区二区三区| 国产一区有黄有色的免费视频| xxx大片免费视频| 人妻 亚洲 视频| 色哟哟·www| 久久国产精品男人的天堂亚洲 | 国产成人精品婷婷| 亚洲av在线观看美女高潮| 国产黄频视频在线观看| 国产人妻一区二区三区在| 在线天堂最新版资源| av国产精品久久久久影院| 91精品一卡2卡3卡4卡| 日本猛色少妇xxxxx猛交久久| 97在线人人人人妻| 日韩中字成人| 国产爱豆传媒在线观看| 免费高清在线观看视频在线观看| 男人爽女人下面视频在线观看| 全区人妻精品视频| 午夜福利视频精品| 国产精品蜜桃在线观看| 亚洲真实伦在线观看| 国产成人aa在线观看| 男人舔奶头视频| 女性被躁到高潮视频| 美女脱内裤让男人舔精品视频| 六月丁香七月| 国产成人一区二区在线| 亚洲精品久久午夜乱码| 欧美高清成人免费视频www| 国产亚洲午夜精品一区二区久久| 又黄又爽又刺激的免费视频.| 菩萨蛮人人尽说江南好唐韦庄| 中文字幕亚洲精品专区| 欧美日韩在线观看h| 国产无遮挡羞羞视频在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美一区二区三区国产| 男女免费视频国产| 我要看黄色一级片免费的| 美女脱内裤让男人舔精品视频| videossex国产| 国产男女超爽视频在线观看| 六月丁香七月| 女性生殖器流出的白浆| 国精品久久久久久国模美| 久久久久国产网址| 久久99热这里只频精品6学生| 欧美+日韩+精品| 色视频在线一区二区三区| 久久av网站| 亚洲国产精品一区三区| 亚洲精华国产精华液的使用体验| 国产精品爽爽va在线观看网站| 涩涩av久久男人的天堂| 亚洲精品亚洲一区二区| 一本一本综合久久| 久久国产精品大桥未久av | 国产一区二区在线观看日韩| 如何舔出高潮| 中文欧美无线码| 国产乱来视频区| 亚洲精品自拍成人| 简卡轻食公司| 欧美精品人与动牲交sv欧美| 老司机影院成人| 人妻夜夜爽99麻豆av| 一本一本综合久久| 亚洲电影在线观看av| 久久午夜福利片| 国产高清不卡午夜福利| 亚洲欧美日韩卡通动漫| 伦精品一区二区三区| 久久青草综合色| 九草在线视频观看| 久久6这里有精品| 五月开心婷婷网| 久久 成人 亚洲| 精品酒店卫生间| 欧美zozozo另类| 亚洲中文av在线| 高清不卡的av网站| 国产在线视频一区二区| 亚洲,一卡二卡三卡| 波野结衣二区三区在线| 一区二区三区免费毛片| 国产在线免费精品| 91精品伊人久久大香线蕉| 亚洲经典国产精华液单| 18禁在线无遮挡免费观看视频| 夫妻性生交免费视频一级片| 精品久久久精品久久久| 国产精品久久久久久久电影| 国产精品久久久久成人av| 欧美最新免费一区二区三区| 精品久久久精品久久久| 夜夜看夜夜爽夜夜摸| 日韩欧美精品免费久久| 色网站视频免费| 国产精品久久久久成人av| 国产成人精品久久久久久| 一本一本综合久久| 亚洲av日韩在线播放| 麻豆成人午夜福利视频| 国产精品一区二区性色av| 简卡轻食公司| 亚洲精品日韩在线中文字幕| 性色av一级| 一区在线观看完整版| 建设人人有责人人尽责人人享有的 | 男人舔奶头视频| 国产乱来视频区| 亚洲精品日韩av片在线观看| 香蕉精品网在线| 亚洲欧美日韩无卡精品| 99热这里只有精品一区| 一区二区av电影网| 极品少妇高潮喷水抽搐| 中国三级夫妇交换| 精品一区二区三卡| 欧美少妇被猛烈插入视频| 人妻 亚洲 视频| 久久人人爽av亚洲精品天堂 | 天美传媒精品一区二区| 五月开心婷婷网| 亚洲综合精品二区| 狂野欧美激情性bbbbbb| 大话2 男鬼变身卡| 免费黄色在线免费观看| 国产精品成人在线| 精品99又大又爽又粗少妇毛片| 搡女人真爽免费视频火全软件| 又爽又黄a免费视频| 午夜激情久久久久久久| kizo精华| 免费观看性生交大片5| 黄色配什么色好看| 一级毛片我不卡| 永久免费av网站大全| 国产成人aa在线观看| 只有这里有精品99| 香蕉精品网在线| 日日撸夜夜添| 美女中出高潮动态图| 亚洲精品国产成人久久av| 国产成人一区二区在线| 在线观看免费高清a一片| 国产av码专区亚洲av| 国产淫片久久久久久久久| 国产精品国产三级专区第一集| 亚洲av综合色区一区| 日本wwww免费看| 黄色欧美视频在线观看| 欧美国产精品一级二级三级 | 一个人看视频在线观看www免费| 免费大片18禁| 国产在线免费精品| 日韩免费高清中文字幕av| 欧美三级亚洲精品| 中文乱码字字幕精品一区二区三区| 中文天堂在线官网| 欧美丝袜亚洲另类| 哪个播放器可以免费观看大片| 直男gayav资源| 自拍欧美九色日韩亚洲蝌蚪91 | 成人免费观看视频高清| 亚洲精品国产av成人精品| 成人毛片a级毛片在线播放| 亚洲精品成人av观看孕妇| 深夜a级毛片| 亚洲图色成人| 国产国拍精品亚洲av在线观看| 91精品一卡2卡3卡4卡| 毛片一级片免费看久久久久| av.在线天堂| 久久久久性生活片| 3wmmmm亚洲av在线观看| 最近最新中文字幕大全电影3| 亚洲欧洲日产国产| 人妻一区二区av| 一边亲一边摸免费视频| 一级毛片我不卡| 欧美三级亚洲精品| 51国产日韩欧美| 亚洲精品日韩av片在线观看| 久久久久久久久久久丰满| 波野结衣二区三区在线| 欧美一区二区亚洲| 欧美区成人在线视频|