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

    橄欖葉綠體基因組密碼子偏好性特征

    2022-11-01 04:20:26賴瑞聯(lián)覃振師陳義挺沈朝貴田奇琳吳如健
    關鍵詞:差異影響分析

    賴瑞聯(lián), 陳 瑾, 馮 新, 覃振師, 陳義挺, 沈朝貴, 田奇琳, 吳如健

    (1.福建省農業(yè)科學院果樹研究所,福建 福州 350013;2.廣西農業(yè)科學院廣西南亞熱帶農業(yè)科學研究所,廣西 崇左532415;3.閩南師范大學生物科學與技術學院,福建 漳州 363000)

    密碼子是生物體遺傳信息的基本單元,自然界中的20種氨基酸除酪氨酸和甲硫氨酸外,其余氨基酸均對應多個密碼子,稱為同義密碼子.同義密碼子在生物體中的使用并不是隨機和均等的,這種使用頻率的差異便形成了密碼子偏好性.在生物體中,同義密碼子的選擇并不會改變蛋白質的種類和特征,但會影響蛋白質的翻譯效率,在長期的物種進化過程中,形成的最優(yōu)密碼子往往集中在高表達基因或多態(tài)性位點中,從而維持選擇的有效性[1-2].在密碼子偏好性特征的形成機理研究中,比較公認的是中性理論和選擇—突變—漂移假說,即密碼子偏好性是對非同義密碼子的突變—漂移以及對最優(yōu)密碼子選擇壓力綜合作用的結果.因此,目前廣泛認為,突變壓力和自然選擇是密碼子偏好性形成的最主要的影響因素[3].不同生物體密碼子偏好性模式不同,通常近緣物種的偏好模式較為一致,開展密碼子偏好性研究有利于進一步揭示物種的起源和進化規(guī)律[4-6].

    葉綠體是高等植物或藻類特有的能量轉換系統(tǒng),葉綠體基因組是存在于葉綠體中,呈雙鏈環(huán)狀,能夠進行基因轉錄和翻譯的相對獨立的遺傳信息載體.由于植物中的葉綠體基因組往往具有高保守性和低進化速率等特征,其常常被用于物種進化規(guī)律和遺傳鑒定等方面的研究.其中,葉綠體基因組密碼子編碼規(guī)律的研究可為揭示植物的分子特征提供重要的參考依據[7].

    橄欖(Canariumalbum)是我國熱帶和亞熱帶地區(qū)的特色果樹,其果實兼具藥食兩用,具有很好的開發(fā)利用前景.通常認為,中國是橄欖種質資源起源和遺傳多樣性分布中心,然而,目前橄欖的相關研究工作主要集中在種質資源、栽培技術、藥理活性、保鮮加工和果實品質等方面[8].橄欖遺傳背景的研究較少,其密碼子偏好性的研究僅在基因個體或轉錄本中有少數(shù)報道[9-10].鑒于此,本研究在前期組裝的橄欖葉綠體基因組的基礎上,進一步對其密碼子偏好性及其影響因素進行研究,同時確定橄欖葉綠體基因組最優(yōu)密碼子,旨在為橄欖葉綠體基因組的進化遺傳學和葉綠體基因工程研究提供參考.

    1 材料與方法

    1.1 基因序列的獲取

    在前期工作中,本項目組采用BGISEQ-500測序平臺完成了橄欖葉綠體基因組的測序,組裝獲得了163 140 bp基因組序列(附件圖Ⅰ,掃OSID碼可見),包含83個蛋白編碼基因,GenBank登錄號為MN217684[11].在基因篩選時,首先剔除重復序列,隨后選擇以ATG作為起始密碼子,TAA、TAG和TGA為終止密碼子,且編碼區(qū)總長度超過300 bp的葉綠體基因序列用于密碼子偏好性分析,最終獲得54個符合條件的蛋白編碼基因.篩選結果如表1所示.

    1.2 橄欖葉綠體基因組密碼子偏好性分析

    1.2.1 密碼子偏好性參數(shù)分析 采用CodonW程序和EMBOSS在線軟件(http://bioinformatics.nl/emboss-explorer/)分析橄欖葉綠體基因組密碼子的偏好性參數(shù),包括有效密碼子數(shù)(effective number of codons, ENc)、密碼子第1、2、3位上的GC含量(GC1s、GC2s、GC3s)、總GC含量(GC)、同義密碼子相對使用度(relative synonymous codon usage, RSCU)、密碼子出現(xiàn)次數(shù)(codon number, CN)、密碼子適應指數(shù)(codon adaptation index, CAI)、密碼子偏好指數(shù)(codon bias index, CBI)和最優(yōu)密碼子使用頻率(frequency of optimal codons, FOP).同時,采用統(tǒng)計學軟件SPSS 19.0進行各參數(shù)間的相關性和顯著性分析.

    1.2.2 中性繪圖分析 密碼子第1、2位上的堿基通常會發(fā)生非同義突變,改變編碼蛋白,影響基因的功能,而第3位上的堿基通常會發(fā)生同義突變.中性繪圖分析通過比較GC12(密碼子第1、2位上的平均GC含量)與GC3s的相關性,衡量密碼子偏好性形成的影響因素.通常認為,GC12與GC3s之間存在顯著相關性,且回歸系數(shù)近似1.0時,密碼子不同位點上的堿基組成沒有差異,密碼子偏好性的形成主要受突變壓力的影響;反之,密碼子不同位點上的堿基組成存在差異,密碼子偏好性的形成主要受自然選擇的影響.

    1.2.3 ENc對應分析 ENc是衡量密碼子偏好性程度的重要指標,取值20~61,值越小說明密碼子偏好性越強,反之則偏好性越弱.根據Wright[12]提出的理論,通過GC3s可以計算獲得ENc的理論值.ENc對應分析的是ENc與GC3s之間的關系,以GC3s為橫坐標,若ENc分布在標準曲線附近,則說明密碼子偏好性受突變壓力的影響;若ENc分布在標準曲線下方較遠的位置,則說明密碼子偏好性主要受自然選擇的影響.

    1.2.4 偏倚分析 偏倚用于衡量密碼子第3位上的AT與GC之間的突變是否平衡.以A3/(A3+T3)為縱坐標,G3/(G3+C3)為橫坐標,由中心點(A=T、C=G)發(fā)出的矢量反映該堿基的偏移程度和方向.若基因集中分布在中心位置附近,則說明4種堿基使用的頻率相對一致,密碼子偏好性主要受突變壓力的影響;若基因偏離中心,位置較遠,說明可能還存在其他影響因素.

    1.2.5 RSCU對應性分析 采用CodonW程序計算各基因的RSCU,基于各葉綠體基因的功能分類將RSCU分布到一個59維的向量空間,進行各指標間的對應性分析.

    1.2.6 最優(yōu)密碼子分析 以CAI作為參考標準,參照Liu[13]的方法將54個基因的CAI由高到低排序,篩選前10%和后10%的基因分別作為高、低表達的基因,分析不同基因組間各密碼子RSCU的差異,同時采用統(tǒng)計學軟件SPSS 19.0分析組間的差異顯著性.若高表達基因的RSCU大于低表達基因,且存在差異顯著性,則將該密碼子確定為基因組的最優(yōu)密碼子.

    2 結果與分析

    2.1 橄欖葉綠體基因組密碼子的偏好性

    橄欖葉綠體基因組密碼子偏好性參數(shù)如表2所示.首先,以54個基因作為一個整體進行分析時,基因組整體的ENc、CAI和CBI分別為50.360、0.186和0.167.不同基因之間的ENc差異較大,取值36.797~58.462,ENc均大于35.0,且CAI和CBI遠小于1.0,說明無論是基因組還是基因個體,密碼子偏好性均較弱,但所有基因的ENc均小于61.0,說明所有基因均存在一定程度的密碼子偏好性.從堿基組成上看,GC、GC1s、GC2s和GC3s分別為0.385、0.466、0.383和0.305,其中,GC1s、GC2s、GC3s和GC在不同基因個體間的取值分別為0.344~0.585、0.283~0.612、0.219~0.480和0.320~0.488,說明GC在密碼子不同位置上的含量并不均勻.從整體上看,基因組更偏向使用A和T,并以A或T結尾的密碼子.

    表2 橄欖葉綠體基因密碼子的偏好性參數(shù)Table 2 Parameters of codon usage preference in C.album chloroplast genome

    橄欖葉綠體基因組密碼子偏好性參數(shù)間的相關性如表3所示.從堿基組成上看,GC與GC1s、GC2s、GC3s之間存在極顯著正相關,GC1s與GC2s之間也存在顯著相關,GC3s與GC1s、GC2s之間的相關性則不顯著,說明密碼子第1位與第2位的堿基組成較相似,但與第3位之間存在較大差異,但任意位置上的堿基組成均會影響密碼子GC含量;ENc只與GC3s存在極顯著相關,表明密碼子第3位上的堿基組成會影響ENc;CN與所有指標之間的相關性都不顯著,說明密碼子出現(xiàn)次數(shù)對整體偏好性沒有影響;GC1s、GC均與CAI、CBI、FOP之間存在顯著或極顯著相關,GC3s與FOP之間也存在顯著相關,說明GC1s和GC會影響CAI和CBI,而GC1s、GC3s和GC均會影響FOP.此外,CAI、CBI、FOP兩兩之間的相關性均達到極顯著水平,說明3個指標間均會互相影響.

    表3 橄欖葉綠體基因密碼子偏好性參數(shù)間的相關系數(shù)1)Table 3 Correlation analysis on parameters of codon usage preference for C.album chloroplast genome

    2.2 橄欖葉綠體基因組密碼子中性繪圖分析

    橄欖葉綠體基因組密碼子中性繪圖分析結果(圖1)顯示,GC3s和GC12分別為0.219~0.480和0.344~0.568,所有基因在平面圖上的分布范圍較為分散,許多基因分布在回歸曲線[y=0.268x+0.356(R2=0.054)]的兩側,此外,也有部分基因的分布位置距離回歸曲線較遠.表明基因組密碼子偏好性在受到突變壓力影響的同時,可能還受到其他作用力的共同影響.

    2.3 橄欖葉綠體基因組密碼子ENc對應分析

    橄欖葉綠體基因組密碼子ENc對應分析結果(圖2)顯示:大部分基因均勻分布在ENc理論值曲線附近,表明這些基因ENc的實際值與理論值較為相近;此外,也有部分基因分布在ENc理論值下方較遠的位置,說明ENc的實際值與理論值之間仍然存在一定的差異.進一步分析發(fā)現(xiàn),ENc比值頻率[(ENc理論值-ENc實際值)/ENc理論值)]為-0.085~0.320.其中,多數(shù)基因分布在-0.1~0.1(表4),占所有基因的74.07%,說明大部分基因ENc的實際值與理論值相差較小,即基因組密碼子偏好性與第3位上堿基組成的差異相關.由此可見,基因組密碼子偏好性受突變壓力的影響,同時可能還受自然選擇的影響.

    2.4 橄欖葉綠體基因組密碼子偏倚分析

    橄欖葉綠體基因組密碼子偏倚分析結果(圖3)顯示,所有基因G3/(G3+C3)、A3/(A3+T3)平均取值的分布位置為(0.479、0.514),說明在總體上,基因組密碼子中A與T的比例接近,G與C的比例接近.然而,從圖4也可以看出,很多基因并非聚集在平面圖中心,而是分布在距離中心點較遠的位置.其中,處于平面圖上半部分基因的A的使用頻率高于T,處于左半圖基因的C的使用頻率高于G,反之亦然.研究認為,當密碼子使用偏好只受突變壓力的影響時,不同堿基的使用頻率相似[14].由此可見,基因組密碼子的偏好性除了受突變壓力的影響,還受自然選擇等其他作用力的影響.

    2.5 橄欖葉綠體基因組密碼子RSCU對應性分析

    以RSCU大于1.0為標準,獲得橄欖葉綠體基因組中偏好性較強的密碼子有30個,然而這些密碼子的RSCU均小于2.0,說明基因組中不存在極強偏好性的密碼子,可能也是造成基因組整體密碼子偏好性較弱的原因.進一步基于基因功能,將54個基因分為4類(光合系統(tǒng)基因、遺傳系統(tǒng)基因、其他基因、未知功能基因),基于RSCU將各基因分布到一個59維的向量空間,第1向量軸(第1軸)、第2向量軸(第2軸)的分布結果(圖4)顯示,4個類型的基因在平面分布上相對分散.其中,第1向量軸顯示有9.71%的差異,第2、3、4向量軸分別存在9.09%、7.90%和6.98%的差異,前4個向量軸累計差異為33.68%.關聯(lián)分析發(fā)現(xiàn),除GC3s(R=0.269,P<0.05)外,各基因在第1軸上的坐標值與ENc(R=0.266)、GC(R=0.046)、CAI(R=-0.325)、CBI(R=-0.216)、CN(R=0.058)之間的相關性并不顯著(P>0.05),進一步說明基因組密碼子的偏好性不僅受單一因素的影響,可能還受突變壓力、自然選擇和其他因素的共同影響.

    2.6 橄欖葉綠體基因組最優(yōu)密碼子

    以CAI作為橄欖葉綠體高、低表達基因的衡量標準,分析高、低表達基因之間的RSCU差異.結果(表5)顯示,高、低表達基因之間的RSCU普遍相似,可能與基因組密碼子偏好性普遍較弱有關.在進行差異顯著性分析后發(fā)現(xiàn),高表達基因CGT、CGC、TTG、CCT、GGT、CAC、TAC中的RSCU顯著高于低表達基因,其中,CGT、CCT、GGT在0.01水平上存在極顯著差異,且RSCU均大于1.0,可作為最優(yōu)密碼子.值得注意的是,這些密碼子均以A或T結尾,與基因組整體密碼子偏好性一致.此外,低表達基因中AGA、AGG、GGA、ATA、TAT的RSCU顯著高于高表達基因,可能是造成這些基因低表達的原因之一.

    表5 橄欖葉綠體基因組高、低表達基因的RSCU比較1)Table 5 Comparison of RSCU in C.album chloroplast genome with high and low levels of expressions

    2.7 橄欖葉綠體基因密碼子偏好性聚類分析

    基于RSCU對橄欖葉綠體基因進行歸類.結果(圖5)顯示,54個基因整體上可以分為3個大類,分別包含26、13和15個基因.然而基于密碼子偏好性的聚類并非按基因功能或類別進行劃分,不同家族成員之間的聚類結果相互交叉,說明基因個體間的密碼子偏好性存在特殊的進化方式,與基因本身功能或所屬類別無明顯相關.

    同一基因在不同密碼子中的RSCU值轉換為Z-score值后以熱圖展示,顏色越紅,基因RSCU值越高;顏色越綠,基因RSCU值越低.圖5 橄欖葉綠體基因密碼子偏好性聚類結果Fig.5 Cluster analysis on preferred codons of C.album chloroplast genome

    3 討論

    密碼子使用偏好性會影響生物體mRNA穩(wěn)定性、mRNA轉錄、蛋白質翻譯準確性以及蛋白折疊等,從而精細調控基因表達[15],其中,ENc、GC和GC3s是密碼子偏好性研究中最常用的指標.葉綠體基因組作為植物中相對獨立的遺傳系統(tǒng),往往具有高度保守性,開展葉綠體基因組密碼子偏好性研究對于揭示植物的分子特征具有重要意義.本研究中,橄欖葉綠體基因組的ENc為50.36,遠遠大于35.0,說明其密碼子偏好性較弱;而從密碼子組成上看,GC和GC3s分別為0.385和0.305,說明基因組偏好使用A和T,并以A或T結尾的密碼子,符合雙子葉植物密碼子偏好性規(guī)律[16].在劍麻(Agavehybrid)[17]、秋茄(Kandeliaobovata)[18]、灰毛漿果楝(Cipadessacinerascens)[19]和翠雀(Delphiniumgrandiflorum)[20]等其他雙子葉植物的研究中也發(fā)現(xiàn),其葉綠體基因組密碼子偏好以A或T結尾,與本研究結果一致.

    RSCU指在編碼同一氨基酸時,某一特定密碼子在所有同義密碼子中出現(xiàn)的相對概率,能夠直觀反映密碼子偏好性.基于對RSCU的分析發(fā)現(xiàn),橄欖葉綠體基因組不存在偏好性極強的密碼子,這可能也是導致整體密碼子偏好性較弱的直接因素.在橄欖葉綠體基因組編碼氨基酸的過程中,對所有類型的密碼子均有使用,與橄欖單基因密碼子編碼特征不同[10].在最優(yōu)密碼子選擇時,橄欖葉綠體基因組偏好使用以A或T結尾的密碼子,與基因組整體堿基選擇偏好一致,而與核基因最優(yōu)密碼子的偏好規(guī)律存在一定的差異[9].研究表明,與核基因組相比,葉綠體基因組較小,結果相對穩(wěn)定,序列高度保守,遺傳重組率低[21],可能是造成橄欖葉綠體基因組和核基因組密碼子偏好性差異的重要原因.從橄欖葉綠體基因個體上看,這種進化方式與基因本身功能或歸屬類別無明顯相關性,在枯草芽孢桿菌(Bacillussubtilis)中也曾發(fā)現(xiàn)不同基因間的密碼子偏好性不完全一致[22].表明橄欖葉綠體基因組與核基因組密碼子偏好性規(guī)律存在相似性和差異性,而基因個體之間的密碼子偏好性進化方式也相對獨立.進一步分析發(fā)現(xiàn),橄欖葉綠體基因組中4個NCG型密碼子的RSCU均遠小于1.0.有研究認為,密碼子CG中的C在甲基化過程中容易脫去氨基變成T,G也容易發(fā)生變異,物種中高NCG容易發(fā)生DNA甲基化導致突變[23].因此,橄欖葉綠體基因組中的NCG型密碼子偏好性維持在較低水平,可能與該物種本身甲基化程度有關.

    密碼子偏好性是研究植物基因組進化規(guī)律的一項重要的參考依據.目前廣泛認為,突變壓力和自然選擇是導致物種密碼子偏好性差異形成的最重要的作用力,然而不同物種中密碼子偏好性形成的主導影響因素不一樣.已有的葉綠體基因組密碼子偏好性的相關研究認為:刺榆(Hemipteleadavidii)和大戟科(Euphorbiaceae)植物主要受自然選擇壓力的影響[24-25];涼粉草(Mesonachinensis)同時受自然選擇和突變壓力的共同作用[26];紫菜(Porphyraumbilicalis)和沙棗(Elaeagnusangustifolia)除了受自然選擇、突變壓力的作用,可能還受堿基組成等其他因素的影響[27-28].本研究結果表明,橄欖葉綠體基因組密碼子偏好性不是單一因素影響的結果,可能受突變壓力、自然選擇和其他因素的共同作用.綜合前人[15,29-30]及本試驗的研究結果認為,在突變壓力、自然選擇和其他作用力的共同影響下,橄欖葉綠體基因組密碼子的部分堿基逐步由GC向AT方向突變,從而形成特定的偏好性特征.隨著生物技術的不斷發(fā)展,越來越多的果樹基因組被組裝完成,在基因組測序的基礎上,系統(tǒng)結合細胞核基因組和細胞器基因組深入研究密碼子使用規(guī)律,進而揭示橄欖的起源和進化也是后續(xù)研究的重要方向.

    猜你喜歡
    差異影響分析
    相似與差異
    音樂探索(2022年2期)2022-05-30 21:01:37
    是什么影響了滑動摩擦力的大小
    哪些顧慮影響擔當?
    當代陜西(2021年2期)2021-03-29 07:41:24
    隱蔽失效適航要求符合性驗證分析
    找句子差異
    電力系統(tǒng)不平衡分析
    電子制作(2018年18期)2018-11-14 01:48:24
    生物為什么會有差異?
    電力系統(tǒng)及其自動化發(fā)展趨勢分析
    擴鏈劑聯(lián)用對PETG擴鏈反應與流變性能的影響
    中國塑料(2016年3期)2016-06-15 20:30:00
    M1型、M2型巨噬細胞及腫瘤相關巨噬細胞中miR-146a表達的差異
    久久精品人妻少妇| 少妇熟女aⅴ在线视频| 精品国内亚洲2022精品成人| 网址你懂的国产日韩在线| 国产午夜精品久久久久久一区二区三区 | 久久久久久九九精品二区国产| av福利片在线观看| 久久热精品热| 中文字幕人成人乱码亚洲影| 日韩欧美精品v在线| 久久久久久久久大av| 美女高潮喷水抽搐中文字幕| 一边摸一边抽搐一进一小说| 国产aⅴ精品一区二区三区波| 99久久九九国产精品国产免费| 大型黄色视频在线免费观看| 国产成人福利小说| 日韩亚洲欧美综合| 美女被艹到高潮喷水动态| 国产毛片a区久久久久| 村上凉子中文字幕在线| 国产伦人伦偷精品视频| 桃色一区二区三区在线观看| 国产成+人综合+亚洲专区| 国产精品久久久久久久电影| 国产成年人精品一区二区| 久久天躁狠狠躁夜夜2o2o| 色综合欧美亚洲国产小说| 搡老熟女国产l中国老女人| 亚洲精品亚洲一区二区| 黄色视频,在线免费观看| 亚洲成人久久爱视频| 亚洲三级黄色毛片| 亚洲 国产 在线| 一二三四社区在线视频社区8| 十八禁人妻一区二区| 国内久久婷婷六月综合欲色啪| 国产精品国产高清国产av| 少妇人妻一区二区三区视频| 欧美在线一区亚洲| 亚洲av电影不卡..在线观看| 婷婷精品国产亚洲av| 国产不卡一卡二| 国产黄色小视频在线观看| 三级毛片av免费| 免费人成在线观看视频色| 熟妇人妻久久中文字幕3abv| 亚洲,欧美,日韩| 国产精品98久久久久久宅男小说| 亚洲国产精品sss在线观看| 久久久久久大精品| 一级av片app| 国产亚洲精品av在线| 男人舔奶头视频| 欧美一区二区国产精品久久精品| 亚洲人成电影免费在线| 午夜福利成人在线免费观看| 中文字幕人妻熟人妻熟丝袜美| 97超级碰碰碰精品色视频在线观看| 美女cb高潮喷水在线观看| 一个人观看的视频www高清免费观看| 日韩欧美精品免费久久 | 国产 一区 欧美 日韩| 欧美潮喷喷水| 日韩欧美在线乱码| 精品久久久久久久人妻蜜臀av| 国产伦精品一区二区三区四那| 久久久久久久亚洲中文字幕 | 变态另类成人亚洲欧美熟女| 女人被狂操c到高潮| 久久午夜亚洲精品久久| 成人av一区二区三区在线看| 国产精品国产高清国产av| h日本视频在线播放| 麻豆一二三区av精品| 天堂网av新在线| 最近最新中文字幕大全电影3| 9191精品国产免费久久| 亚洲av二区三区四区| 好男人电影高清在线观看| 精品一区二区免费观看| 久久久久久久午夜电影| 亚洲精品日韩av片在线观看| 深夜精品福利| 极品教师在线免费播放| 男女之事视频高清在线观看| 舔av片在线| 成人特级av手机在线观看| 国产精品三级大全| 亚洲熟妇熟女久久| 99国产综合亚洲精品| 国产成人影院久久av| 国产视频内射| 成人美女网站在线观看视频| 乱人视频在线观看| 有码 亚洲区| 美女黄网站色视频| 国产真实伦视频高清在线观看 | 高潮久久久久久久久久久不卡| 91狼人影院| 欧美乱色亚洲激情| 看黄色毛片网站| 色5月婷婷丁香| 脱女人内裤的视频| 搡老熟女国产l中国老女人| 欧美区成人在线视频| 国产精品综合久久久久久久免费| 久久久久久大精品| 国产精品久久电影中文字幕| 精品国产亚洲在线| 欧美一区二区国产精品久久精品| 亚洲一区高清亚洲精品| 成年版毛片免费区| 欧美xxxx性猛交bbbb| 国产aⅴ精品一区二区三区波| 国产一区二区三区视频了| 成年女人永久免费观看视频| 国产精品亚洲av一区麻豆| 少妇熟女aⅴ在线视频| 国产精品av视频在线免费观看| 欧美xxxx黑人xx丫x性爽| 午夜久久久久精精品| 午夜福利在线观看免费完整高清在 | 国产成人aa在线观看| 天堂av国产一区二区熟女人妻| 国产精品久久久久久亚洲av鲁大| 日本a在线网址| 国产久久久一区二区三区| 久久香蕉精品热| 88av欧美| 久久久久性生活片| 日本熟妇午夜| 97超视频在线观看视频| 99热这里只有是精品在线观看 | 欧美成狂野欧美在线观看| 精品久久久久久久人妻蜜臀av| 亚洲成人精品中文字幕电影| 天天躁日日操中文字幕| 欧美潮喷喷水| 一区福利在线观看| 毛片女人毛片| 一本综合久久免费| 久久精品国产清高在天天线| 99久久精品国产亚洲精品| 日本熟妇午夜| 久久人人爽人人爽人人片va | 动漫黄色视频在线观看| 我要看日韩黄色一级片| 欧美日韩亚洲国产一区二区在线观看| 成年女人看的毛片在线观看| 国产一区二区激情短视频| 桃色一区二区三区在线观看| 九九在线视频观看精品| 天堂动漫精品| 午夜精品在线福利| 欧美绝顶高潮抽搐喷水| 国产成年人精品一区二区| 99热6这里只有精品| 国产精品亚洲美女久久久| 99久久精品热视频| 在线看三级毛片| 国产毛片a区久久久久| 五月玫瑰六月丁香| 国产单亲对白刺激| 国产精品一区二区性色av| 听说在线观看完整版免费高清| 如何舔出高潮| 欧美三级亚洲精品| 国产三级黄色录像| 在现免费观看毛片| 亚洲人成网站在线播| 啦啦啦韩国在线观看视频| 欧美日韩黄片免| 久久精品国产自在天天线| 亚洲内射少妇av| 免费人成视频x8x8入口观看| 日本一二三区视频观看| 亚洲av成人精品一区久久| 一进一出抽搐gif免费好疼| 欧美精品国产亚洲| 国产人妻一区二区三区在| 精品一区二区三区av网在线观看| 日韩精品中文字幕看吧| 亚洲色图av天堂| 日本黄大片高清| eeuss影院久久| 国产黄a三级三级三级人| 俺也久久电影网| 中文字幕人妻熟人妻熟丝袜美| 成人国产综合亚洲| av天堂在线播放| 免费搜索国产男女视频| 91在线观看av| 精品乱码久久久久久99久播| xxxwww97欧美| 亚洲一区二区三区色噜噜| 99国产综合亚洲精品| 精品久久久久久久久久久久久| 精品人妻偷拍中文字幕| 一个人观看的视频www高清免费观看| 免费在线观看日本一区| 99久久精品热视频| 免费看日本二区| 国产一级毛片七仙女欲春2| 亚洲av一区综合| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美成人a在线观看| 中文字幕av在线有码专区| 色哟哟哟哟哟哟| 亚洲七黄色美女视频| 国产色婷婷99| 国产高清激情床上av| 成人av在线播放网站| 成人高潮视频无遮挡免费网站| 国产成年人精品一区二区| 欧美中文日本在线观看视频| 91午夜精品亚洲一区二区三区 | www.熟女人妻精品国产| 日韩成人在线观看一区二区三区| 亚洲人与动物交配视频| 99视频精品全部免费 在线| 亚洲最大成人av| 亚洲天堂国产精品一区在线| 精华霜和精华液先用哪个| 1024手机看黄色片| 国产视频一区二区在线看| 听说在线观看完整版免费高清| av黄色大香蕉| av国产免费在线观看| 真人一进一出gif抽搐免费| 欧美日韩瑟瑟在线播放| 中文字幕熟女人妻在线| 九色国产91popny在线| 国产真实伦视频高清在线观看 | 国产欧美日韩精品一区二区| 最后的刺客免费高清国语| 91九色精品人成在线观看| 亚洲av免费在线观看| 久久草成人影院| 日韩欧美 国产精品| 久久久久久久精品吃奶| 免费一级毛片在线播放高清视频| 久久中文看片网| 色视频www国产| 日本黄色视频三级网站网址| 啪啪无遮挡十八禁网站| 熟妇人妻久久中文字幕3abv| 噜噜噜噜噜久久久久久91| 九色国产91popny在线| 丁香欧美五月| 很黄的视频免费| 成年女人看的毛片在线观看| 欧美成人一区二区免费高清观看| 9191精品国产免费久久| 亚洲国产欧洲综合997久久,| 午夜久久久久精精品| 亚洲成人精品中文字幕电影| 国内精品久久久久久久电影| 99久久久亚洲精品蜜臀av| 亚洲五月天丁香| 看片在线看免费视频| 一a级毛片在线观看| 国产成人aa在线观看| 欧美区成人在线视频| 国产伦人伦偷精品视频| 免费一级毛片在线播放高清视频| www.999成人在线观看| 国产69精品久久久久777片| 国内精品一区二区在线观看| 国内揄拍国产精品人妻在线| www.www免费av| 最近最新中文字幕大全电影3| 国产毛片a区久久久久| avwww免费| 久久久久九九精品影院| 久久国产精品人妻蜜桃| 国产一区二区三区视频了| 欧美一区二区国产精品久久精品| 日本黄色片子视频| 激情在线观看视频在线高清| 久久99热这里只有精品18| 夜夜夜夜夜久久久久| 两个人的视频大全免费| 高清毛片免费观看视频网站| 日韩国内少妇激情av| АⅤ资源中文在线天堂| 淫妇啪啪啪对白视频| 亚洲成人精品中文字幕电影| 桃色一区二区三区在线观看| 在线观看av片永久免费下载| 午夜免费激情av| 人妻制服诱惑在线中文字幕| 我要看日韩黄色一级片| 欧美高清成人免费视频www| 久久人人爽人人爽人人片va | 麻豆久久精品国产亚洲av| 亚洲五月天丁香| 欧美高清成人免费视频www| 欧美黑人欧美精品刺激| 精品久久久久久成人av| 深夜a级毛片| 美女 人体艺术 gogo| 三级国产精品欧美在线观看| 给我免费播放毛片高清在线观看| 欧美区成人在线视频| 日韩av在线大香蕉| 欧美绝顶高潮抽搐喷水| 噜噜噜噜噜久久久久久91| 亚洲在线观看片| 久久婷婷人人爽人人干人人爱| 国产精品伦人一区二区| 日韩成人在线观看一区二区三区| 午夜激情欧美在线| 97碰自拍视频| 精品无人区乱码1区二区| 12—13女人毛片做爰片一| 日韩欧美精品免费久久 | 一区二区三区免费毛片| 中文字幕av在线有码专区| 五月伊人婷婷丁香| 成人午夜高清在线视频| 国产69精品久久久久777片| 9191精品国产免费久久| 婷婷精品国产亚洲av在线| 日本a在线网址| 国产成人啪精品午夜网站| 亚洲 欧美 日韩 在线 免费| a级毛片a级免费在线| 色综合欧美亚洲国产小说| 日韩中文字幕欧美一区二区| 久久人人爽人人爽人人片va | 国产av一区在线观看免费| 一区二区三区高清视频在线| 90打野战视频偷拍视频| 蜜桃亚洲精品一区二区三区| 伦理电影大哥的女人| 男人舔女人下体高潮全视频| 国产亚洲精品久久久com| 国产精品不卡视频一区二区 | 国产成人aa在线观看| 日日摸夜夜添夜夜添av毛片 | 在线观看一区二区三区| 国产高清有码在线观看视频| 亚洲美女搞黄在线观看 | 亚洲成av人片免费观看| 波多野结衣巨乳人妻| 国产野战对白在线观看| 天美传媒精品一区二区| 欧美另类亚洲清纯唯美| 欧美日韩综合久久久久久 | 欧美日韩国产亚洲二区| 欧美在线黄色| 日本黄大片高清| 亚洲性夜色夜夜综合| 亚洲第一欧美日韩一区二区三区| 婷婷精品国产亚洲av| 最近视频中文字幕2019在线8| 黄色丝袜av网址大全| 嫩草影视91久久| 久久久久久久精品吃奶| 欧美最黄视频在线播放免费| 日韩欧美一区二区三区在线观看| 色视频www国产| 日本免费一区二区三区高清不卡| 国产精品久久视频播放| 内射极品少妇av片p| 色综合站精品国产| 国产在视频线在精品| 在线观看av片永久免费下载| 狠狠狠狠99中文字幕| 欧美高清成人免费视频www| 亚洲综合色惰| 在线免费观看不下载黄p国产 | 成人高潮视频无遮挡免费网站| 久久久久久国产a免费观看| 最近最新中文字幕大全电影3| 亚洲美女黄片视频| 老鸭窝网址在线观看| 搡老熟女国产l中国老女人| 亚洲av.av天堂| 日本撒尿小便嘘嘘汇集6| 欧美日韩综合久久久久久 | 3wmmmm亚洲av在线观看| 色播亚洲综合网| 男女那种视频在线观看| 99国产极品粉嫩在线观看| 两个人视频免费观看高清| 中文字幕av在线有码专区| 香蕉av资源在线| 国产精品美女特级片免费视频播放器| 欧美日本视频| 变态另类丝袜制服| 国产在线精品亚洲第一网站| 国产午夜精品久久久久久一区二区三区 | 中文字幕熟女人妻在线| 欧美日韩亚洲国产一区二区在线观看| 日本免费一区二区三区高清不卡| www.www免费av| 亚洲va日本ⅴa欧美va伊人久久| 国产三级黄色录像| av在线老鸭窝| 中文字幕人成人乱码亚洲影| 最新在线观看一区二区三区| 国产爱豆传媒在线观看| 性色avwww在线观看| 麻豆一二三区av精品| 99热只有精品国产| 久久精品国产亚洲av天美| 最好的美女福利视频网| 一进一出抽搐动态| 午夜福利在线观看免费完整高清在 | 日韩av在线大香蕉| 欧美日本视频| 国产91精品成人一区二区三区| 色av中文字幕| 欧美又色又爽又黄视频| 亚洲av第一区精品v没综合| 国产69精品久久久久777片| 黄色一级大片看看| 久久久久久国产a免费观看| 中文在线观看免费www的网站| 91字幕亚洲| 极品教师在线免费播放| 国产久久久一区二区三区| 2021天堂中文幕一二区在线观| 免费人成视频x8x8入口观看| 老女人水多毛片| 精华霜和精华液先用哪个| 能在线免费观看的黄片| 精品一区二区三区视频在线观看免费| 性插视频无遮挡在线免费观看| 一个人看视频在线观看www免费| or卡值多少钱| 久久久精品大字幕| 他把我摸到了高潮在线观看| 欧美另类亚洲清纯唯美| 国产真实伦视频高清在线观看 | 国产成人啪精品午夜网站| 国产真实伦视频高清在线观看 | 直男gayav资源| 一区福利在线观看| 国产单亲对白刺激| av中文乱码字幕在线| 无人区码免费观看不卡| 午夜两性在线视频| 亚洲av熟女| 久久久久久久久久成人| 我要看日韩黄色一级片| 啦啦啦韩国在线观看视频| 色哟哟哟哟哟哟| 噜噜噜噜噜久久久久久91| 免费观看精品视频网站| 国产高清视频在线观看网站| 丰满人妻熟妇乱又伦精品不卡| 精品久久久久久久久亚洲 | 国产 一区 欧美 日韩| 国产免费av片在线观看野外av| 99热这里只有精品一区| 婷婷精品国产亚洲av在线| 51国产日韩欧美| 免费在线观看影片大全网站| 亚洲欧美清纯卡通| 久久国产精品影院| 免费搜索国产男女视频| 欧美日韩瑟瑟在线播放| 免费大片18禁| 在线观看舔阴道视频| 欧美性猛交╳xxx乱大交人| 美女被艹到高潮喷水动态| 国产av不卡久久| 国内精品久久久久精免费| 欧美黑人欧美精品刺激| 乱码一卡2卡4卡精品| 久久精品夜夜夜夜夜久久蜜豆| 性欧美人与动物交配| 人人妻,人人澡人人爽秒播| 欧美成狂野欧美在线观看| 欧美日韩黄片免| 国产欧美日韩一区二区精品| 国内精品美女久久久久久| 成年人黄色毛片网站| 日本黄大片高清| 国产欧美日韩一区二区精品| 成年女人看的毛片在线观看| 国产亚洲精品久久久com| 久久精品91蜜桃| 国产精品嫩草影院av在线观看 | 少妇熟女aⅴ在线视频| 99国产精品一区二区蜜桃av| 特级一级黄色大片| 国内毛片毛片毛片毛片毛片| 中文资源天堂在线| 色哟哟·www| 中文字幕人妻熟人妻熟丝袜美| 日韩欧美国产一区二区入口| 日韩欧美精品免费久久 | 色综合亚洲欧美另类图片| 非洲黑人性xxxx精品又粗又长| 欧美日韩瑟瑟在线播放| АⅤ资源中文在线天堂| 亚洲男人的天堂狠狠| 91午夜精品亚洲一区二区三区 | 极品教师在线免费播放| 久久国产精品人妻蜜桃| 日韩欧美精品v在线| 男人的好看免费观看在线视频| 91在线观看av| 亚州av有码| 成年女人毛片免费观看观看9| 亚洲 欧美 日韩 在线 免费| 舔av片在线| 国产精品一区二区免费欧美| 久久国产乱子伦精品免费另类| 久久九九热精品免费| 黄色一级大片看看| 欧美一区二区国产精品久久精品| 99热这里只有精品一区| 日韩人妻高清精品专区| 最后的刺客免费高清国语| 亚洲成av人片在线播放无| 日韩欧美三级三区| 亚洲欧美精品综合久久99| 动漫黄色视频在线观看| 亚洲七黄色美女视频| 久久久精品大字幕| avwww免费| 一级黄片播放器| 熟女电影av网| 免费观看人在逋| 国产av麻豆久久久久久久| 精品国产三级普通话版| 国产午夜福利久久久久久| 一区二区三区四区激情视频 | 少妇人妻一区二区三区视频| 简卡轻食公司| 好男人电影高清在线观看| 国产欧美日韩一区二区精品| 色哟哟哟哟哟哟| 久9热在线精品视频| 精品一区二区三区视频在线观看免费| 一边摸一边抽搐一进一小说| 一级黄色大片毛片| 亚洲精品影视一区二区三区av| 成年版毛片免费区| 色播亚洲综合网| 国产精品98久久久久久宅男小说| 性插视频无遮挡在线免费观看| 变态另类成人亚洲欧美熟女| 性插视频无遮挡在线免费观看| 国产精品1区2区在线观看.| 一个人观看的视频www高清免费观看| 国产毛片a区久久久久| 性插视频无遮挡在线免费观看| 婷婷亚洲欧美| 69av精品久久久久久| 亚洲色图av天堂| 精品久久久久久久末码| 久久香蕉精品热| 午夜精品在线福利| 亚洲电影在线观看av| 国产精华一区二区三区| 亚洲一区二区三区色噜噜| 可以在线观看的亚洲视频| 国产免费一级a男人的天堂| 好男人在线观看高清免费视频| 女同久久另类99精品国产91| 国模一区二区三区四区视频| 国产老妇女一区| 一本一本综合久久| 老熟妇仑乱视频hdxx| 美女 人体艺术 gogo| 最好的美女福利视频网| 如何舔出高潮| 久久99热这里只有精品18| 在线播放国产精品三级| 国产成年人精品一区二区| 国产精品野战在线观看| 欧美性猛交黑人性爽| 欧美+日韩+精品| 欧美在线一区亚洲| 国产精品美女特级片免费视频播放器| 动漫黄色视频在线观看| 午夜影院日韩av| 欧美高清性xxxxhd video| 免费高清视频大片| 人人妻人人看人人澡| 搡老妇女老女人老熟妇| 亚洲国产精品合色在线| 亚洲黑人精品在线| 欧美三级亚洲精品| 亚洲人成伊人成综合网2020| 亚洲成人精品中文字幕电影| 色av中文字幕| 亚洲,欧美精品.| 免费人成在线观看视频色| 精品国内亚洲2022精品成人| 男女做爰动态图高潮gif福利片| 成人亚洲精品av一区二区| 色噜噜av男人的天堂激情| 99热这里只有是精品在线观看 | 国产精品女同一区二区软件 | 久久久久久久久久成人| 少妇裸体淫交视频免费看高清| 国产免费男女视频| 国产高清有码在线观看视频| 一个人看视频在线观看www免费| 最近视频中文字幕2019在线8| 亚洲av中文字字幕乱码综合| 变态另类丝袜制服| 91九色精品人成在线观看| 啪啪无遮挡十八禁网站| 国产精品98久久久久久宅男小说| 男插女下体视频免费在线播放| 国产午夜福利久久久久久| 十八禁人妻一区二区|