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

    芍藥花色調(diào)控基因的密碼子使用模式及其影響因素分析

    2016-07-18 09:32:30吳彥慶趙大球
    關(guān)鍵詞:芍藥影響因素

    吳彥慶,趙大球,陶 俊

    (揚(yáng)州大學(xué)園藝與植物保護(hù)學(xué)院/江蘇省作物遺傳生理重點(diǎn)實(shí)驗(yàn)室,江蘇揚(yáng)州 22509)

    ?

    芍藥花色調(diào)控基因的密碼子使用模式及其影響因素分析

    吳彥慶,趙大球,陶俊

    (揚(yáng)州大學(xué)園藝與植物保護(hù)學(xué)院/江蘇省作物遺傳生理重點(diǎn)實(shí)驗(yàn)室,江蘇揚(yáng)州 22509)

    摘要:【目的】芍藥花色的優(yōu)劣影響其觀賞價(jià)值和商業(yè)價(jià)值,研究芍藥花色調(diào)控基因的密碼子使用偏好性和密碼子使用模式的影響因素,為芍藥花色調(diào)控基因在mRNA翻譯、轉(zhuǎn)基因設(shè)計(jì)、新基因表達(dá)與功能預(yù)測(cè)以及分子生物進(jìn)化研究提供參考?!痉椒ā扛鶕?jù)前期芍藥花色嵌合體品種‘金輝’轉(zhuǎn)錄組測(cè)序篩選的6 345個(gè)芍藥花色調(diào)控基因,并根據(jù)CDS序列特征和大于300 bp原則進(jìn)行過濾后最終獲得的2 234個(gè)基因序列作為研究對(duì)象,利用Mobyle軟件計(jì)算GC含量、第1與2位密碼子的平均GC含量(GC12)、第3位密碼子的GC含量(GC3s)、有效密碼子數(shù)ENC、密碼子適應(yīng)指數(shù)CAI、相對(duì)同義密碼子使用度RSCU等密碼子偏性指標(biāo),其次進(jìn)行中性繪圖(GC12 vs. GC3)、ENC-GC3s繪圖以及PR2(Parity Rule 2)繪圖分析,并運(yùn)用多元統(tǒng)計(jì)分析方法探討突變壓力和選擇作用對(duì)密碼子使用模式的影響程度,最后以5%CAI值作為高、低表達(dá)樣本組,計(jì)算這兩個(gè)樣本組的同義密碼子相對(duì)使用度,利用卡方檢驗(yàn)Chi-square test分析兩組之間的顯著性差異來確定最優(yōu)密碼子?!窘Y(jié)果】芍藥花色相關(guān)基因的密碼子GC3s含量為46.37%,大部分基因GC含量主要分布在30%—55%;中性繪圖分析表明GC3s與GC12呈極顯著的正相關(guān)(R2=0.202,P<0.01);ENC-GC3s繪圖表明大部分基因分布在標(biāo)準(zhǔn)曲線周圍,也有一部分基因分布在標(biāo)準(zhǔn)曲線下方較遠(yuǎn)的位置,同時(shí)大部分基因(ENCexp-ENCobs)/ENCexp比值集中分布在0.0—0.4;PR2繪圖分析顯示密碼子第三位T的使用頻率高于A,C使用頻率高于G,表明嘌呤(A和G)與嘧啶(T和C)的使用頻率并不均衡;對(duì)應(yīng)性分析COA(Correspondence Analysis)表明,第一軸上顯示了38.09%的差異,其他3個(gè)軸分別為18.42%、15.09%、14.59%,表明芍藥花色調(diào)控基因的密碼子使用模式評(píng)價(jià)以第一軸(Axis 1)為主;突變壓力和選擇作用分析發(fā)現(xiàn),第一主軸與GC3s、CAI的相關(guān)系數(shù)均達(dá)到極顯著正相關(guān)(R2=0.736,P<0.01;R2=0.286,P<0.01);利用△RSCU和卡方顯著性檢驗(yàn)的方法,確定了21個(gè)為芍藥花色調(diào)控相關(guān)基因的最優(yōu)密碼子,其中18個(gè)以G或C結(jié)尾,僅CGU、GGU等2個(gè)密碼子以U結(jié)尾。【結(jié)論】芍藥花色調(diào)控基因的最優(yōu)密碼子多數(shù)以G/C結(jié)尾,并且密碼子使用模式主要受到堿基差異(R2=0.736)和基因表達(dá)水平(R2=0.286)共同作用的影響,其中堿基差異占主導(dǎo)因素。本研究了解了芍藥花色調(diào)控基因的密碼子使用模式情況,為通過密碼子改造開展芍藥花色遺傳改良以及分子進(jìn)化研究提供了一定的理論依據(jù)。

    關(guān)鍵詞:芍藥;花色調(diào)控基因;密碼子使用模式;影響因素

    聯(lián)系方式:吳彥慶,E-mail:yqwu19880928@126.com。通信作者陶俊,E-mail:taojun@yzu.edu.cn

    0 引言

    【研究意義】遺傳密碼是連接DNA和蛋白質(zhì)的重要橋梁,每種氨基酸至少對(duì)應(yīng)一個(gè)遺傳密碼子(一般不超過6個(gè)),編碼同一種氨基酸的密碼子稱為同義密碼子。同義密碼子在使用頻率中存在一定的差異,這種現(xiàn)象稱為密碼子使用偏好性(codon usage bias,CUB),并且在某一特定物種或者基因中比較傾向使用的同義密碼子稱為最優(yōu)密碼子(optimal codons)[1]。不同物種間基因組的 CUB現(xiàn)象可能不同,表明這些物種在進(jìn)化過程中可能受到的突變壓力和選擇作用的影響不同[2],因此,密碼子偏好性的量化,能夠幫助有效理解物種機(jī)體的進(jìn)化發(fā)展[2],同時(shí)密碼子偏好性的分析還能夠有效的幫助同義密碼子使用偏性相關(guān)機(jī)制的理解[3]。其次,密碼子偏性的分析有助于了解轉(zhuǎn)錄和翻譯進(jìn)程中的調(diào)控機(jī)制,在預(yù)測(cè)外源基因的最適宿主和通過改良外源基因以提高其表達(dá)水平等方面也具有重要生物學(xué)意義[4]?!厩叭搜芯窟M(jìn)展】影響密碼子偏好性的因素較多,目前主要集中在突變壓力(如GC含量[5]、基因堿基組成[6])和選擇作用(如基因表達(dá)水平[7]、tRNA豐度[8]、蛋白結(jié)構(gòu)與長(zhǎng)度[9]、翻譯的起始信號(hào)[10])兩個(gè)方面。目前密碼子偏好性的研究在不斷發(fā)展,從研究方向來看,對(duì)于特定基因密碼子偏好性,以前主要研究原核生物以及低等真核生物中密碼子偏好性對(duì)基因表達(dá)的影響[11],隨后開始關(guān)注高等動(dòng)植物的密碼子偏好性[4,12],而對(duì)于基因組密碼子偏好性來說,到目前為止主要集中在單細(xì)胞和模式生物中,包括大腸桿菌(Escherichia coli)[13]、酵母(Yeast)[14]、衣原體(Chlamydiae)[15]、新桿狀線蟲(Caenorhabditis)、果蠅(Drosophila)和擬南芥(Arabidopsis)[16]等。從研究理論來看,密碼子偏好性現(xiàn)象存在2種理論,即中性理論和選擇-突變-漂變學(xué)說[17]?!颈狙芯壳腥朦c(diǎn)】芍藥(Paeonia lactiflora Pall.)為芍藥科(Paeoniaceae)芍藥屬植物,是原產(chǎn)中國(guó)的傳統(tǒng)名花,與花中之王牡丹并稱為“花中二絕”,有“花相”之美譽(yù),目前關(guān)于芍藥基因密碼子偏好性分析的報(bào)道相對(duì)較少[18]。隨著新一代測(cè)序技術(shù)的不斷成熟,許多物種的基因組信息被揭示,然而芍藥基因組信息目前尚未公布,課題組前期利用轉(zhuǎn)錄組測(cè)序(RNA-Seq)在芍藥花色嵌合體品種‘金輝’中篩選出6 345個(gè)花色調(diào)控相關(guān)的差異基因[19],作為本次分析密碼子使用模式的研究對(duì)象。芍藥花色的優(yōu)劣不僅影響到觀賞植物的觀賞價(jià)值,而且直接關(guān)系到其商業(yè)開發(fā)價(jià)值,因此,分析芍藥花色調(diào)控基因的密碼子使用模式對(duì)芍藥花色育種具有重要的現(xiàn)實(shí)意義?!緮M解決的關(guān)鍵問題】利用多元統(tǒng)計(jì)分析(multivariate statistical analysis)方法探討芍藥花色調(diào)控基因的密碼子使用模式以及其影響因素,不僅對(duì)提高芍藥花色基因的外源表達(dá)水平具有指導(dǎo)意義,同時(shí)也為今后深入研究芍藥花色調(diào)控的分子機(jī)制和改善芍藥花色育種奠定理論基礎(chǔ)。

    1 材料與方法

    1.1 序列數(shù)據(jù)來源

    芍藥花色嵌合體品種‘金輝’花瓣材料,于2013 年5月采自揚(yáng)州大學(xué)園藝與植物保護(hù)學(xué)院芍藥種質(zhì)資源圃,課題組前期對(duì)其進(jìn)行了轉(zhuǎn)錄組測(cè)序,2014年完成測(cè)序工作和數(shù)據(jù)分析。6 345個(gè)芍藥花色調(diào)控基因序列來自芍藥花色嵌合體品種‘金輝’的轉(zhuǎn)錄組測(cè)序(RNA-Seq)[19],為了降低取樣誤差,編寫C程序?qū)π蛄羞M(jìn)行篩選處理,提取滿足下列條件的CDS序列共2 234條:以ATG為起始密碼子;并以TAA、TAG 或TGA為終止密碼子;編碼序列長(zhǎng)度大于300 bp[20]。

    1.2 堿基組成分析及中性繪圖

    利用 Mobyle軟件[21](http://www.molbiol.ox.ac. uk/cu,version 1.4.2)計(jì)算每個(gè)基因序列的堿基組成,并統(tǒng)計(jì)以下指標(biāo):(1)第3位密碼子各堿基含量(A3s、U3s、C3s、G3s);(2)第1密碼子GC含量(GC1)與第2位密碼子GC含量(GC2);(3)密碼子整體GC含量(GC);第1與2位密碼子的平均GC含量(GC12);(4)第三位密碼子的GC含量(GC3s)。中性繪圖(neutrality plot)即GC12與GC3s的相關(guān)性分析,是衡量選擇與突變對(duì)密碼子使用模式影響程度的一種分析方法,繪圖以GC12位縱坐標(biāo),GC3s為橫坐標(biāo)。如果GC12與GC3s之間呈顯著相關(guān),則說明密碼子3個(gè)位置上堿基組成無(wú)差異,密碼子的使用主要受到突變壓力的影響;如果GC12與GC3s相關(guān)性不顯著,說明密碼子第1、2位和第3位堿基組成不同,基因組GC含量高度保守,密碼子的使用更多地是受選擇作用的影響[22]。

    1.3 同義密碼子使用偏好性分析

    有效密碼子數(shù)ENC(effective number of codon),是評(píng)估基因整體密碼子偏好性的一個(gè)有效指標(biāo),其數(shù)值范圍為20(每個(gè)氨基酸只使用一個(gè)同義密碼子的極端偏好情況)到61(每個(gè)同義密碼子被平均使用的無(wú)偏好情況),ENC值越小表明密碼子偏好性越強(qiáng)[23]。目前研究發(fā)現(xiàn)當(dāng)ENC值小于或者等于35時(shí),基因密碼子具有顯著的使用偏好性[24]。密碼子適應(yīng)指數(shù)CAI (codon adaptation index)是評(píng)估密碼子偏好性程度從而揭示基因表達(dá)水平的一個(gè)重要指標(biāo)[25],CAI數(shù)值范圍在0—1,其數(shù)值越大表明密碼子偏好性越強(qiáng)。相對(duì)同義密碼子使用度RSCU(relative synonymous codon usage)是指某一特定密碼子在使用頻率與其無(wú)偏好性使用時(shí)預(yù)期頻率之間的比值,也是衡量密碼子偏好性程度的有效指標(biāo)[26],RSCU比值等于1說明該密碼子無(wú)使用偏好性,若RSCU比值大于1說明該密碼子的使用頻率較高。本研究中2 234條基因序列的ENC、CAI以及RSCU值均由Mobyle軟件中CodonW程序進(jìn)行計(jì)算。

    1.4 ENC繪圖分析

    ENC繪圖(ENC-plot)以ENC值為縱坐標(biāo),GC3s為橫坐標(biāo),是分析各基因密碼子使用特征,并探究基因堿基組成和密碼子偏好性之間關(guān)系的一個(gè)有效手段[23]。ENC-plot被廣泛用于分析影響密碼子使用的主要因素:如果密碼子偏好性只受突變壓力的影響,基因沿標(biāo)準(zhǔn)曲線分布或落在標(biāo)準(zhǔn)曲線附近;如果密碼子偏好性只受到選擇作用的影響,基因應(yīng)該落在標(biāo)準(zhǔn)曲線下方較遠(yuǎn)的位置。

    1.5 對(duì)應(yīng)性分析

    對(duì)應(yīng)性分析COA(correspondence analysis)是利用多元統(tǒng)計(jì)方法探究樣本各變量之間的關(guān)系,從而揭示影響芍藥花色調(diào)控基因密碼子使用模式的主要因素,COA分析基于RSCU值將樣本中所有基因按密碼子的使用頻率分布在一個(gè)59維(64個(gè)密碼子去除3個(gè)終止密碼子以及甲硫氨酸和色氨酸的密碼子)的向量空間[27]。

    1.6 PR2繪圖分析

    在計(jì)算第3位密碼子核苷酸含量(A3s、U3s、C3s、G3s)基礎(chǔ)上,分別計(jì)算A3/(A3+U3)與G3/(G3+C3)值,PR2(Parity Rule 2),繪圖是以A3/(A3+U3)為縱坐標(biāo),G3/(G3+C3)為橫坐標(biāo)[28]。

    1.7 最優(yōu)密碼子測(cè)定

    本研究中最優(yōu)密碼子的確定參照YANG[29]的分析方法,最優(yōu)密碼子的選擇以CAI值為偏性標(biāo)準(zhǔn),確定高低表達(dá)樣本。具體方法如下:通過密碼子適應(yīng)指數(shù)(CAI)對(duì)所篩選出的所有序列進(jìn)行排序,從排列好樣本總數(shù)的前后端各取5%作為高、低表達(dá)樣本組,計(jì)算這兩個(gè)樣本組的同義密碼子相對(duì)使用度(RSCU),利用卡方檢驗(yàn)Chi-square test分析兩組之間的顯著性差異,將高表達(dá)組中RSCU值極顯著高于低表達(dá)組(P<0.01)的密碼子定義為最優(yōu)密碼子。

    1.8 統(tǒng)計(jì)分析

    利用Mobyle軟件中CodonW程序計(jì)算所有基因序列堿基組成和密碼子偏好性指標(biāo),利用 SPSS 18.0 (http://www.spss.com/)軟件基于Spearman方法進(jìn)行相關(guān)性和卡方檢驗(yàn)Chi-square test分析。

    2 結(jié)果

    2.1 GC含量分布以及中性繪圖

    密碼子偏好性 CUB很大程度上受到各基因整體堿基含量的影響,整體GC含量往往反映了方向性突變的強(qiáng)弱,尤其是同義密碼子的主要差別體現(xiàn)在第3位堿基上(GC3s)[30]。在2 234個(gè)芍藥花色調(diào)控基因序列中,GC含量變化范圍為24.6%—73.6%(標(biāo)準(zhǔn)差SD=6.56),其中大部分基因GC含量主要分布在30%—55%(圖1)。整體來看,平均GC和AU含量分別為46.26%和53.74%。GC1為37.05%, GC2為29.90%,GC3s為46.37%,GC3s含量與整體基因編碼區(qū)一致,與GC1和GC2存在一定的差異。中性繪圖分析(Neutrality analysis)以GC3s為橫坐標(biāo),GC12為縱坐標(biāo)(圖 2),相關(guān)性分析顯示 GC3s與GC12呈極顯著的正相關(guān)(R2=0.202,P<0.01),表明芍藥花色調(diào)控基因密碼子的使用主要受到突變壓力的影響。

    圖1 芍藥花色調(diào)控基因的GC含量分布Fig. 1 Distribution of the GC contents of regulating color genes in P. lactiflora

    圖2 中性繪圖(GC12 vs. GC3)Fig. 2 Neutrality plots (GC12 vs. GC3)

    2.2 ENC與GC3s的關(guān)聯(lián)分析

    ENC-GC3s繪圖以各基因 ENC值為縱坐標(biāo),以GC3s值為橫坐標(biāo)(圖3),基因ENC值分布在14—61,GC3s值分布在0.085—0.916。由圖3可知,大部分基因分布在標(biāo)準(zhǔn)曲線周圍,也有一部分基因分布在標(biāo)準(zhǔn)曲線下方較遠(yuǎn)的位置,由此表明芍藥花色調(diào)控基因的密碼子偏好性不僅受到突變壓力的影響,也一定程度上受到選擇作用的影響。為了準(zhǔn)確評(píng)估ENC觀察值(ENCobs)與預(yù)期值(ENCexp)的差異,本研究計(jì)算了(ENCexp-ENCobs)/ENCexp的比值,結(jié)果見圖4,(ENCexp-ENCobs)/ENCexp比值最高峰分布在 0.1—0.2,大部分基因其比值集中分布在 0.0—0.4,表明大部分基因的ENCexp與ENCobs值相差很小,芍藥花色調(diào)控基因密碼子偏好性差異與GC3s的差異有關(guān),表明密碼子偏好性主要受突變壓力的影響。

    2.3 PR2-plot分析

    利用PR2-plot(Parity Rule 2-plot)的方法分析了各基因密碼子中4個(gè)氨基酸家族嘌呤(A和G)與嘧啶(T和C)之間的關(guān)系,由圖5表明,密碼子第3 位T的使用頻率高于A,C使用頻率高于G,4個(gè)堿基的不均衡使用表明芍藥花色調(diào)控基因密碼子偏好性的影響因素除突變壓力外,還受到其他方面的影響(如選擇作用)。

    圖3 ENC-GC3s繪圖Fig. 3 ENC-GC3s plot

    圖4 ENC比值頻率分布Fig. 4 Frequency distribution of effective number of codons (ENC) ratio

    紅色圓圈表示每個(gè)點(diǎn)平均分布位置,值為:x=0.4983072±0.12896476,y=0.4224539±0.1327275Red open circle indicates the average position for each plot, calculated as follows: x=0.4983072±0.12896476, y=0.4224539±0.1327275圖5 PR2繪圖[(A3/(A3+T3) vs. G3/(G3+C3)]Fig. 5 PR2-bias plot [A3/(A3+T3) against G3/(G3+C3)]

    2.4 對(duì)應(yīng)性分析COA

    基于RSCU值進(jìn)行了對(duì)應(yīng)性分析COA,結(jié)果發(fā)現(xiàn)第一軸顯示了 38.09%差異,其他 3個(gè)軸分別為18.42%、15.09%、14.59%。芍藥花色調(diào)控基因密碼子使用模式評(píng)價(jià)分析以第一軸(Axis 1)為主,本研究以各基因第一軸(Axis 1)與第二軸(Axis 2)進(jìn)行GC含量分布的對(duì)應(yīng)性分析(圖6),大部分基因集中在Axis 1(-0.5—1.5)和Axis 2(-0.5—1.0)范圍內(nèi),GC含量小于45%的基因主要分布在軸左下方區(qū)域,GC含量位于 45%—60%的基因主要分布在軸中間區(qū)域。此外,還分析了兩軸中密碼子G/C和A/U結(jié)尾情況(圖 7),結(jié)果發(fā)現(xiàn)第一軸中各基因密碼子偏好性差異可能由于以G/C和A/U結(jié)尾的使用頻率不同導(dǎo)致。

    表1 第一主軸與堿基組成的相關(guān)性分析Table 1 Summary of correlation between the first major principle axes and nucleotide composition

    以第一軸和第二軸為主的基因分布;紅色、藍(lán)色和綠色圓點(diǎn)分別代表GC≥60%、45%≤GC<60%、GC<45%The distribution of genes is shown along the first and second axes. Red, blue and green dots indicate genes with G+C content ≥60%, ≥45% but <60%, and <45%, respectively圖6 密碼子使用模式的對(duì)應(yīng)性分析Fig. 6 Correspondence analysis of codon usage patterns

    2.5 突變壓力和選擇作用對(duì)芍藥花色調(diào)控基因的密碼子使用模式的影響

    目前突變壓力和選擇作用是影響密碼子使用模式的主要因素[17]。為了評(píng)估突變壓力對(duì)芍藥花色調(diào)控基因密碼子使用模式的影響,本研究根據(jù)密碼子使用模式評(píng)價(jià)主要軸(Axis 1)為對(duì)象,探討堿基組成對(duì)其影響程度。結(jié)果發(fā)現(xiàn)(表1),第一軸Axis 1與GC3s呈極顯著的正相關(guān)(R2=0.736,P<0.01),表明突變壓力中堿基組成(R2=0.736)是影響芍藥花色調(diào)控基因密碼子使用模式的主要因素。此外,為了探討選擇作用中基因表達(dá)水平、蛋白長(zhǎng)度等影響因素,同時(shí)目前CAI值是被用于評(píng)估基因表達(dá)水平的重要指標(biāo)[31-32]。因此,本研究首先分析了CAI值、蛋白長(zhǎng)度(Protein length)對(duì)密碼子偏好性(ENC值)的影響(圖8、圖9),發(fā)現(xiàn)CAI值與ENC值均表現(xiàn)為極顯著的負(fù)相關(guān)(R2=-0.278,P<0.01),而蛋白長(zhǎng)度與ENC值呈負(fù)相關(guān),但不顯著(R2=-0.005,P>0.05),表明基因表達(dá)水平對(duì)密碼子偏好性具有顯著的影響。因此,本研究以基因表達(dá)水平作為選擇作用因素,進(jìn)一步分析CAI值對(duì)主要軸Axis 1的影響程度,發(fā)現(xiàn)Axis 1與CAI值表現(xiàn)為極顯著的正相關(guān)(R2=0.286,P<0.01),表明選擇作用中基因表達(dá)水平(R2=0.286)對(duì)芍藥花色調(diào)控基因的密碼子使用模式具有一定的影響。

    以第一軸和第二軸為主的密碼子分布;紅色圓點(diǎn)代表以A/U結(jié)尾的密碼子,藍(lán)色圓點(diǎn)代表C/G結(jié)尾的密碼子The distribution of codons is shown along the first and second axes. Note: codons ending with A and U are shown in red dots; codons ending with C and G are shown in green dots圖7 密碼子使用模式的對(duì)應(yīng)性分析Fig. 7 Correspondence analysis of codon usage patterns

    圖8 ENC與蛋白長(zhǎng)度的關(guān)系繪圖Fig. 8 Plot of ENC versus encoded protein length

    圖9 ENC與基因表達(dá)水平的關(guān)系繪圖Fig. 9 Plot of ENC versus gene expression level

    2.6 最優(yōu)密碼子(Optimal codons)分析

    本研究依據(jù)CAI值大小對(duì)芍藥花色調(diào)控基因進(jìn)行高、低表達(dá)樣本庫(kù)比較(表 2),結(jié)果表明,共發(fā)現(xiàn)21個(gè)密碼子在高、低表達(dá)樣本庫(kù)中具有極顯著的差異(P<0.01),這些密碼子為最優(yōu)密碼子,其中19個(gè)最優(yōu)密碼子以G或C結(jié)尾,僅CGU和GGU 2個(gè)密碼子以U結(jié)尾。

    表2 芍藥花色調(diào)控基因中高低表達(dá)水平的密碼子使用頻率比較Table 2 Comparison of codon usage frequencies between regulating color gene in P. lactiflora with high and low levels of expression

    3 討論

    每個(gè)生物在長(zhǎng)期進(jìn)化過程中都會(huì)形成一種特定的密碼子使用模式,其中GC含量是生物基因組中堿基組成的一個(gè)重要指標(biāo),在基因組的演變中具有重要意義。GC含量往往反映了方向性突變的強(qiáng)弱,尤其是同義密碼子的主要差別體現(xiàn)在第3位堿基上(GC3s),由于密碼子第3位上堿基受到的突變壓力較小,因此,GC3s通常被作為分析密碼子使用模式的一個(gè)重要參數(shù)[33]。本研究分析了2 234個(gè)芍藥花色調(diào)控基因的GC含量分布情況,發(fā)現(xiàn)平均GC含量和GC3s比較接近(均略小于 50%),表明所有花色調(diào)控基因中整體AU含量略高于GC,且密碼子稍微偏向以A/U結(jié)尾。目前研究發(fā)現(xiàn),在小麥(Triticum aestivum)、大麥(Hordeum vulgare)、水稻(Oryza sativa)等植物基因組中均表現(xiàn)較高的GC含量和偏向以G/C結(jié)尾[34],然而在真菌、一些單細(xì)胞微生物如鐮狀瘧原蟲(Plasmodium falciparum)和支原體(Mycoplasma capricolum)以及植物線粒體微生物基因組中AU含量明顯高于GC含量[35-37]。此外,密碼子使用模式在形成過程中往往受到很多因素的影響,其中主要包括突變和選擇[17]。在本研究中,中性繪圖顯示GC12與GC3s之間具有極顯著的正相關(guān);ENC-plot分析發(fā)現(xiàn)大部分基因分布在標(biāo)準(zhǔn)曲線周圍,也有一部分基因分布在標(biāo)準(zhǔn)曲線下方較遠(yuǎn)的位置;PR2-plot分析發(fā)現(xiàn)A、C、T 和G 4個(gè)堿基使用不均衡,結(jié)合中性繪圖、ENC-plot 和PR2-plot綜合分析,表明芍藥花色調(diào)控基因的密碼子使用模式可能受到突變和選擇等多重因素的影響。

    本研究發(fā)現(xiàn)芍藥花色調(diào)控基因序列堿基相關(guān)的因素(A3s、U3s、C3s、G3s、GC、GC3s)與第一主軸(Axis 1)間的相關(guān)系數(shù)均呈現(xiàn)極顯著相關(guān)(R2=-0.417,P<0.01;R2=-0.246,P<0.01;R2=0.549,P<0.01;R2=-0.054,P<0.01;R2=0.525,P<0.01;R2=0.736,P<0.01),表明密碼子使用特點(diǎn)受堿基組成差異影響較大。此外,在已分化的多細(xì)胞真核生物不同組織以及不同發(fā)育階段中,確定某個(gè)基因的表達(dá)水平是非常困難的。在未知的芍藥基因組中,很難通過EST序列來預(yù)測(cè)單個(gè)基因的表達(dá)水平,密碼子適應(yīng)指數(shù)CAI一直被廣泛用于評(píng)估基因表達(dá)水平的重要指標(biāo)[31-32]?;虮磉_(dá)水平作為選擇作用的主要因素,本研究分析CAI值與第一向量主軸間的相關(guān)系數(shù)呈現(xiàn)極顯著的正相關(guān)(R2=0.286,P<0.01),以上兩者分析得到芍藥花色調(diào)控基因的密碼子使用模式主要受到堿基差異(R2=0.736)和基因表達(dá)水平(R2=0.286)共同作用的影響,其中堿基差異占了主導(dǎo)因素。目前已完成全基因組測(cè)序的植物并不多,續(xù)晨等[38]分析葡萄基因組偏好性主要受到堿基差異(R2=0.925)和選擇作用(R2=0.193)共同作用的影響,突變壓力占了主導(dǎo)因素;劉慶坡等[39]在水稻基因組研究中,發(fā)現(xiàn)基因表達(dá)水平(R2=-0.608)的影響程度明顯大于堿基組成(R2=0.344);續(xù)晨等[40]在蝴蝶蘭葉綠體基因組中發(fā)現(xiàn),其密碼子偏好性受堿基差異(R2=0.614)和選擇作用(R2=0.508)共同影響,兩者影響程度相差不大。分析基因組密碼子使用模式及其影響因素是一個(gè)復(fù)雜的過程,相對(duì)基因組序列來說,由于本研究分析的芍藥花色調(diào)控基因數(shù)量較小,結(jié)果可能與芍藥基因組密碼子使用特點(diǎn)有所區(qū)別,并且目前芍藥基因組序列尚未公布,因此,本研究初步揭示了芍藥花色調(diào)控相關(guān)基因的密碼子使用特點(diǎn),為今后進(jìn)一步在分子水平上研究毛茛科植物的系統(tǒng)進(jìn)化以及開展芍藥花色相關(guān)基因的功能研究具有指導(dǎo)意義。

    關(guān)于最優(yōu)密碼子的確定方法,目前報(bào)道不完全一致,許多學(xué)者比較了高偏性與低偏性基因庫(kù)RSCU值,利用△RSCU大小等級(jí)劃分來確定最優(yōu)密碼子[39-42]。此外,YANG等[29]在此基礎(chǔ)上,進(jìn)一步利用卡方顯著性檢驗(yàn)來尋找高、低表達(dá)組具有極顯著差異的密碼子作為最優(yōu)密碼子,本研究參考這種方法,最終將CGC、GCC、UGC和GGU等21個(gè)密碼子確定為芍藥花色調(diào)控基因主要偏愛的密碼子。在最終確定的21個(gè)密碼子中除CGU和GGU外,其余密碼子均以G或C結(jié)尾,表明芍藥花色調(diào)控基因密碼子偏好性可能與第 3 位GC含量呈正相關(guān)。在對(duì)水稻[43]、小麥[44]、玉米[20]等高等植物基因組密碼子使用的研究中發(fā)現(xiàn),最優(yōu)密碼子較多的是以 C/G結(jié)尾,且通常表現(xiàn)出對(duì)嘧啶堿基C的偏好強(qiáng)于嘌呤堿基G[45],這與本研究結(jié)果相一致。因此,本研究確定了芍藥花色調(diào)控基因的最優(yōu)密碼子,今后可以通過密碼子優(yōu)化設(shè)計(jì)對(duì)外源基因進(jìn)行密碼子改造,從而提高外源基因在芍藥中的表達(dá)水平,為今后從遺傳本質(zhì)上進(jìn)行芍藥花色育種改良提供理論參考。

    4 結(jié)論

    利用 Mobyle軟件并結(jié)合多元統(tǒng)計(jì)方法分析了芍藥花色調(diào)控基因的密碼子使用模式及其主要影響因素,發(fā)現(xiàn)GC3s與GC含量比較接近(約為46.3%)。中性繪圖、ENC-plot以及PR2-plot綜合分析表明芍藥花色相關(guān)基因密碼子使用模式主要受突變壓力和選擇作用的影響,進(jìn)一步通過對(duì)應(yīng)性和主軸相關(guān)性分析發(fā)現(xiàn),第一主軸與堿基突變指標(biāo)(GC3s)、基因表達(dá)水平指標(biāo)(CAI)的相關(guān)系數(shù)均達(dá)到極顯著正相關(guān)(R2=0.736,P<0.01;R2=0.286,P<0.01),表明密碼子偏好性主要受到堿基差異(R2=0.736)和基因表達(dá)水平(R2=0.286)共同作用的影響,其中堿基差異占了主導(dǎo)因素。此外,確定了21個(gè)為芍藥花色調(diào)控相關(guān)基因的最優(yōu)密碼子,且大多數(shù)以G/C結(jié)尾,對(duì)基因工程中外源基因的密碼子改良及提高其表達(dá)水平研究提供了參考。

    References

    [1] IKEMURA T. Codon usage and tRNA content in unicellular and multicellular organisms. Molecular Biology and Evolution, 1985, 2:13-34.

    [2] ANGELLOTTI M C, BHUIYAN S B, CHEN G, WAN X F. CodonO:Codon usage bias analysis within and across genomes. Nucleic Acids Research, 2007, 35: 132-136.

    [3] LU H, ZHAO W M, ZHENG Y, WANG H, QI M, YU X P. Analysis of synonymous codon usage bias in Chlamydia. Acta Biochimica Biophysica Sinica (Shanghai), 2005, 37(1): 1-10.

    [4] 郭秀麗, 王玉, 楊路成, 丁兆堂. 茶樹 CBF1基因密碼子使用特性分析. 遺傳, 2012, 34(12): 1614-1623. GUO X L, WANG Y, YANG L C, DING Z T. Analysis of codon use features of CBF gene in Camellia sinensis. Hereditas, 2012, 34(12):1614-1623. (in Chinese)

    [5] 石秀凡, 黃京飛, 柳樹群, 劉次全. 人類基因同義密碼子偏好的特征以及與基因 GC含量的關(guān)系. 生物化學(xué)與生物物理進(jìn)展, 2002,29(3): 411-414. SHI X F, HUANG J F, LIU S Q, LIU C Q. The features of synonymous codon bias and GC-content relationship in human genes. Progress in Biochemistry and Biophysics, 2002, 29(3): 411-414. (in Chinese)

    [6] FEDOROV A, SAXONOV S, GILBERT W. Regular ties of contextdependent codon bias in eukaryotic genes. Nucleic Acids Research,2002, 30: 1192-1197.

    [7] HIRAOKA Y, KAWAMATA K, HARAGUCHI T, CHIKASHIGE Y. Codon usage bias is correlated with gene expression levels in the fission yeast Schizosaccharomyces pombe. Genes to Cells, 2009, 14:499-509.

    [8] MORIYAMA E, POWELL J. Codon usage bias and tRNA abundance in Drosophila. Journal of Molecular Evolution, 1997, 45(5): 514-523.

    [9] ORESIC M, DEHN M, KORENBLUM D, SHALLOWAY D. Tracing specific synonymous codon-secondary structure correlations through evolution. Journal of Molecular Evolution, 2003, 56: 473-484.

    [10] SAKAI H, WASHIO T, SAITO R, SHINAGAWA A, ITOH M,SHIBATA K, CARNINCI P, KONNO H, KAWAI J, HAYASHIZAKI Y, TOMITA M. Correlation between sequENCe conservation of the 5′untranslated region and codon usage bias in Mus musculus genes. Gene, 2001, 276(1/2): 101-105.

    [11] GUSTAFSSON C, GOVINDARAJAN S, MINSHULL J. Codon bias and heterologous protein expression. Trends in Biotechnology, 2004,22(7): 346-353.

    [12] 吳正常, 王靖, 趙喬輝, 朱世平, 訾臣, 吳圣龍, 包文斌. 豬脂多糖結(jié)合蛋白基因(LBP)的密碼子偏好性分析. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào),2013, 21(10): 1135-1144. WU Z C, WANG J, ZHAO Q H, ZHU S P, ZI C, WU S L, BAO W B. Analysis of codon bias of lipopolysaccharide-binding protein gene (LBP) in pigs (Sus scrofa). Journal of Agricultural Biotechnology,2013, 21(10): 1135-1144. (in Chinese)

    [13] Sharp P M, Bailes E, Grocock R J, Peden J F, Sockett R E. Variation in the strength of selected codon usage bias among bacteria. Nucleic Acids Research, 2005, 33: 1141-1153.

    [14] SHARP P M, TUOHY T M, MOSURSKI K R. Codon usage in yeast:Cluster analysis clearly differentiates highly and lowly expressed genes. Nucleic Acids Research, 1986, 14: 5125-5143.

    [15] ZHOU T, SUN X, LU Z. Synonymous codon usage in environmental chlamydia UWE25 reflects an evolutional divergence from pathogenic chlamydiae. Gene, 2006, 368: 117-125.

    [16] DURET L, MOUCHIROUD D. Expression pattern and, surprisingly,gene length shape codon usage in Caenorhabditis, Drosophila, and Arabidopsis. Proceedings of the National Academy of Sciences of the United States of America, 1999, 96: 4482-4487.

    [17] BULMER M. The selection-mutation-drift theory of synonymous codon usage. Genetics, 1991, 129(3): 897-907.

    [18] 吳彥慶, 葛金濤, 陶俊. 芍藥 AP1(APETALA1)基因密碼子使用的偏好性分析. 湖南農(nóng)業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版), 2015, 41(6):610-615. WU Y Q, GE J T, TAO J. Analysis the codon bias of AP1 gene in Paeonia lactiflora. Journal of Hunan Agricultural University (Natural Sciences Edition), 2015, 41(6): 610-615. (in Chinese)

    [19] ZHAO D Q, JIANG Y, NING C L, MENG J S, LIN S, DING W, TAO J. Transcriptome sequencing of a chimaera reveals coordinated expression of anthocyanin biosynthetic genes mediating yellow formation in herbaceous peony (Paeonia lactiflora Pall.). BMC Genomics, 2014, 15(1): 689.

    [20] LIU H M, HE R, ZHANG H Y, HUANG Y B, TIAN M L, ZHANG J J. Analysis of synonymous codon usage in Zea mays. Molecular Biology Reports, 2010, 37(2): 677-684.

    [21] PEDEN J F. CodonW. University of Nottingham, 1999.

    [22] SUEOKA N. Directional mutation pressure and neutral molecular evolution. Proceedings of the National Academy of Sciences, 1988,85(8): 2653-2657.

    [23] WRIGHT F. The ‘effective number of codons' used in a gene. Gene,1990, 87: 23-29.

    [24] COMERON J M, AGUADE M. An evaluation of measures of synonymous codon usage bias. Journal of Molecular Evolution, 1998,47: 268-274.

    [25] SHARP P M, LI W H. The codon adaptation index-a measure of directional synonymous codon usage bias, and its potential applications. Nucleic Acids Research, 1987, 15: 1281-1295.

    [26] SHARP P M, LI W H. An evolutionary perspective on synonymous codon usage in unicellular organisms. Journal of Molecular Evolution,1986, 24(1/2): 28-38.

    [27] GREENACRE M J. Theory and Applications of Correspondence Analysis. London: Academic Press 1984.

    [28] NOBORU S. Near homogeneity of PR2-bias fingerrints in the human genome and their implications in phylogenetic analyses. Journal of Molecular Evolution, 2001, 53: 469-476.

    [29] YANG X, LUO X N, CAI X P. Analysis of codon usage pattern in Taenia saginata based on a transcriptome dataset. Parasites & Vectors,2014, 7: 527.

    [30] ZHOU J H, DING Y Z, HE Y, CHU Y F, ZHAO P, ZHAO P, MA L Y,WANG X J, LI X R, LIU Y S. The effect of multiple evolutionary selections on synonymous codon usage of genes in the Mycoplasma bovis Genome. Plos One, 2014, 9: e108949.

    [31] NAYA H, ROMERO H, CARELS N, ZAVALA A, MUSTO H. Translational selection shapes codon usage in the GC-rich genome of Chlamydomonas reinhardtii. Febs Letters, 2001, 501: 127-130.

    [32] GUPTA S, BHATTACHARYYA T, GHOSH T C. Synonymous codon usage in Lactococcus lactis: mutational bias versus translational selection. Journal of Biomolecular Structure & Dynamics, 2004, 21:527-535.

    [33] BELLGARD M, SCHIBECI D, TRIFONOV E, GOJOBORI T. Early detection of G+C differences in bacterial species inferred from the comparative analysis of the two completely sequenced Helicobacter pylori strains. Journal of Molecular Evolution, 2001, 53(4/5):465-468.

    [34] KAWABE A, MIYASHITA N T. Patterns of codon usage bias in three dicot and four monocot plant species. Genes & Genetic Systems, 2003,78: 343-352.

    [35] SAUL A, BATTISTUTTA D. Codon usage in Plasmodium falciparum. Molecular and Biochemical Parasitology, 1988, 27: 35-42.

    [36] Muto A, Yamao F, Osawa S. The genome of Mycoplasma capricolum. Progress in Nucleic Acid Research and Molecular Biology, 1987, 34:29-58.

    [37] DYBVIG K, VOELKER L L. Molecular biology of mycoplasmas. Annual Reviews in Microbiology, 1996, 50(1): 25-57.

    [38] 續(xù)晨, 蔡小寧, 錢保俐, 賁愛玲. 葡萄基因組密碼子使用偏好模式研究. 西北植物學(xué)報(bào), 2012, 32(2): 409-415. XU C, CAI X N, QIAN B L, BEN A L. Codon usage bias in Vitis vinifera. Acta Botanica Boreali-Occidentalia Sinica, 2012, 32(2):409-415. (in Chinese)

    [39] 劉慶坡, 薛慶中. 粳稻葉綠體基因組的密碼子用法. 作物學(xué)報(bào),2004, 30(12): 1220-1224. LIU Q P, XUE Q Z. Codon Usage in the chloroplast genome of rice (Oryza sativa L.ssp. japonica). Acta Agronomica Sinica, 2004, 30(12):1220-1224. (in Chinese)

    [40] 續(xù)晨, 賁愛玲, 蔡曉寧. 蝴蝶蘭葉綠體基因組密碼子使用的相關(guān)分析. 分子植物育種, 2010, 8(5): 945-950. XU C, BEN A L, CAI X L. Analysis of synonymous codon usage in chloroplast genome of Phalaenopsis aphrodite subsp. Formosana. Molecular Plant Breeding, 2010, 8(5): 945-950. (in Chinese)

    [41] 尚明照, 劉方, 華金平, 王坤波. 陸地棉葉綠體基因組密碼子使用偏性的分析. 中國(guó)農(nóng)業(yè)科學(xué), 2011, 44(2): 245-253. SHANG M Z, LIU F, HUA J P, WANG K B. Analysis on codon usage of chloroplast genome of Gossypium hirsutum. Scientia Agricultura Sinica, 2011, 44(2): 245-253. (in Chinese)

    [42] 李秀璋, 宋輝, 李春杰. 茄腐鐮孢(Fusarium solani)線粒體基因組密碼子偏好性分析. 基因組學(xué)與應(yīng)用生物學(xué), 2015, 34(11):2465-2472. LI X Z, SONG H, LI C J. Analysis of codon usage bias in mitochondria genome of Fusarium solani. Genomics and Applied Biology, 2015, 34(11): 2465-2472. (in Chinese)

    [43] LIU Q P, FENG Y, ZHAO X A, DONG H, XUE Q Z. Synonymous codon usage bias in Oryza sativa. Plant Science, 2004, 167(1): 101-105.

    [44] ZHANG W J, ZHOU J, LI Z F, WANG L, GU X, ZHONG Y. Comparative analysis of codon usage patterns among mitochondrion,chloroplast and nuclear genes in Triticum aestivum L. Journal of Integrative Plant Biology, 2007, 49(2): 246-254.

    [45] WANG L J, ROOSSINCK M J. Comparative analysis of expressed sequences reveals a conserved pattern of optimal codon usage in plants. Plant Molecular Biology, 2006, 61(4): 699-710.

    (責(zé)任編輯 趙伶俐)

    Analysis of Codon Usage Pattern of Paeonia lactiflora Genes Regulating Flower Color and Its Influence Factors

    WU Yan-qing, ZHAO Da-qiu, TAO Jun
    (Jiangsu Key Laboratory of Crop Genetics and Physiology/College of Horticulture and Plant Protection, Yangzhou University,Yangzhou 225009, Jiangsu)

    Abstract:【Objective】 The quality of Paeonia lactiflora flower color affects its ornamental value and the commercial value of ornamental plants. This study aims to understand the codon usage pattern of genes regulating flower color and probe into the mainfactors affecting the formation of codon bias, which has important biological significance for mRNA translation, design of transgenes,the prediction of expression level and functions of new genes, and studies of molecular biology and evolution, etc. 【Method】In a previous study, 6,345 differential genes were screened out by transcriptome sequencing of a flower color chimaera cultivar “Jinhui” with a consistent genetic background red outer-petal and yellow inner-petal, followed by a further filtering analysis according to the principle of CDS sequence characteristics and greater than 300 bp. We finally obtained 2,234 genes as our research object. Mobyle software was used to calculate different parameters for the codon usage, such as GC content, average GC content of the first and second positions (GC12), GC content of the third position (GC3s), effective number of codon (ENC), codon adaptation index (CAI), and relative synonymous codon usage (RSCU). Further analysis of a neutrality plot (GC12 vs. GC3), an ENC-GC3s plot, and a Parity Rule 2 (PR2) plot were performed. Additionally, we probed into the influence of mutational pressure and translational selection by a multivariate statistical analysis. Finally, we took 5% CAI value as high-expression and low-expression sample groups, then calculated the RSCU value, and analyzed the significant difference to determine the optimal codons by a chi-square test.【Result】The results showed that the GC content at the third position of codons was 46.37%. The GC content of most genes was mainly distributed between 30% and 55%. Neutrality analysis showed that there was a significant positive correlation (R2=0.202, P<0.01) between GC3s and GC12 value. The ENC-plot showed most of the genes on or close to the expected curve, but also some points with low-ENC values were below it. The (ENCexp-ENCobs)/ENCexp ratio of most genes ranged from -0.05 to 0.05. The Parity Rule 2-plot showed that the frequency of T nucleotide at the third position was higher than A, and C was higher than G,suggesting that the use frequencies of four nucleotide were not balanced. Correspondence analysis showed that the first axis showed a 38.09% variation, while the other three axes showed 18.42%, 15.09%, and 14.59%, respectively, suggesting that the first axis was the main index evaluating the codon usage bias of Paeonia lactiflora genes regulating flower color. Mutation pressure and selection analysis showed there were significant negative correlations (R2=0.736, P<0.01. R2=0.286, P<0.01) between the first axis and GC3s, CAI value, respectively. Using the delta RSCU and significant chi-square test methods, we defined 21 codons as the major preference codons in the Paeonia lactiflora genes regulating flower color, of 18 codons ending with G or C, only CGU and GGU ending with U.【Conclusion】In conclusion, most of optimal codons ended with G or C. Meanwhile, the codon usage pattern of Paeonia lactiflora genes regulating flower color is formed under the effect of mutational pressure (R2=0.736) and translational selection (R2=0.286), but a mutational bias was the major influence on codon usage. This study not only preliminarily reveals the codon usage pattern of Paeonia lactiflora genes, but also provides a certain theoretical basis for further carrying out genetic improvement of Paeonia lactiflora flower color by codon reconstruction and analyzing the molecular evolution.

    Key words:Paeony; flower color controlling gene; codon usage pattern; influence factor

    收稿日期:2015-12-31;接受日期:2016-03-24

    基金項(xiàng)目:國(guó)家自然科學(xué)基金(31372097,31400592)、江蘇省高校自然科學(xué)研究重大項(xiàng)目(13KJA210005)

    猜你喜歡
    芍藥影響因素
    芍藥鮮切花 美景變“錢”景
    陸抑非《芍藥》
    美麗芍藥化學(xué)成分的研究
    中成藥(2019年12期)2020-01-04 02:02:44
    水驅(qū)油效率影響因素研究進(jìn)展
    突發(fā)事件下應(yīng)急物資保障能力影響因素研究
    環(huán)衛(wèi)工人生存狀況的調(diào)查分析
    農(nóng)業(yè)生產(chǎn)性服務(wù)業(yè)需求影響因素分析
    商(2016年27期)2016-10-17 07:09:07
    村級(jí)發(fā)展互助資金組織的運(yùn)行效率研究
    商(2016年27期)2016-10-17 04:40:12
    基于系統(tǒng)論的煤層瓦斯壓力測(cè)定影響因素分析
    科技視界(2016年20期)2016-09-29 13:45:52
    芍藥為誰(shuí)生
    火花(2015年3期)2015-02-27 07:40:49
    久久久a久久爽久久v久久| 日韩成人av中文字幕在线观看 | 一本一本综合久久| 中文字幕免费在线视频6| 中文资源天堂在线| 国产高潮美女av| 身体一侧抽搐| 欧美3d第一页| 国产亚洲91精品色在线| 午夜福利在线观看免费完整高清在 | 国产欧美日韩精品一区二区| 老司机福利观看| 我要搜黄色片| 成人永久免费在线观看视频| 色尼玛亚洲综合影院| 插阴视频在线观看视频| 深爱激情五月婷婷| 我要看日韩黄色一级片| 男女下面进入的视频免费午夜| 日本三级黄在线观看| av在线天堂中文字幕| 久久精品影院6| 日产精品乱码卡一卡2卡三| 免费无遮挡裸体视频| 91久久精品电影网| 国产精品一区二区三区四区久久| 亚洲av成人av| 一级a爱片免费观看的视频| 免费无遮挡裸体视频| 天堂影院成人在线观看| 免费观看的影片在线观看| 麻豆国产97在线/欧美| 热99在线观看视频| 国产不卡一卡二| 欧美bdsm另类| 久久午夜福利片| 欧美3d第一页| 最近最新中文字幕大全电影3| 女的被弄到高潮叫床怎么办| 少妇猛男粗大的猛烈进出视频 | 午夜爱爱视频在线播放| 亚洲av免费高清在线观看| 最新在线观看一区二区三区| 亚洲一区高清亚洲精品| 少妇猛男粗大的猛烈进出视频 | 国产精品嫩草影院av在线观看| 国产欧美日韩精品亚洲av| 嫩草影院入口| 日韩欧美国产在线观看| 精品一区二区三区视频在线观看免费| 日本色播在线视频| 嫩草影视91久久| 大又大粗又爽又黄少妇毛片口| 成人性生交大片免费视频hd| 日本a在线网址| 俄罗斯特黄特色一大片| 男女之事视频高清在线观看| 一边摸一边抽搐一进一小说| 成人午夜高清在线视频| 亚洲精品粉嫩美女一区| 亚洲精品国产成人久久av| 亚洲成人久久爱视频| 久久久精品大字幕| 亚洲中文日韩欧美视频| 久久久精品大字幕| videossex国产| 最近中文字幕高清免费大全6| 天堂√8在线中文| 99热网站在线观看| 99久久精品热视频| 十八禁网站免费在线| 十八禁网站免费在线| 日产精品乱码卡一卡2卡三| 亚洲色图av天堂| 网址你懂的国产日韩在线| 波野结衣二区三区在线| 亚洲色图av天堂| 亚洲性久久影院| .国产精品久久| 日本 av在线| or卡值多少钱| 亚洲av熟女| 成人漫画全彩无遮挡| 日韩精品有码人妻一区| 一a级毛片在线观看| 日产精品乱码卡一卡2卡三| 国产探花极品一区二区| a级一级毛片免费在线观看| 色综合色国产| 欧美不卡视频在线免费观看| 精品人妻一区二区三区麻豆 | 国产精品av视频在线免费观看| 寂寞人妻少妇视频99o| 少妇熟女aⅴ在线视频| 全区人妻精品视频| 免费不卡的大黄色大毛片视频在线观看 | 1024手机看黄色片| or卡值多少钱| avwww免费| 亚洲精品成人久久久久久| 丰满人妻一区二区三区视频av| av免费在线看不卡| 国产精品av视频在线免费观看| 日韩国内少妇激情av| 级片在线观看| 中国美白少妇内射xxxbb| 精品久久久久久久人妻蜜臀av| 国产伦在线观看视频一区| 久久人人精品亚洲av| 波多野结衣高清作品| 成年女人永久免费观看视频| avwww免费| av在线蜜桃| 在线a可以看的网站| 99久久九九国产精品国产免费| 看免费成人av毛片| 亚洲一区二区三区色噜噜| 1000部很黄的大片| 身体一侧抽搐| 欧美一区二区精品小视频在线| 亚洲第一区二区三区不卡| 黄色日韩在线| 你懂的网址亚洲精品在线观看 | 久久久久久大精品| 69av精品久久久久久| 女的被弄到高潮叫床怎么办| 身体一侧抽搐| 嫩草影院新地址| 国产精品永久免费网站| 中文字幕av成人在线电影| 舔av片在线| 国产精品伦人一区二区| 网址你懂的国产日韩在线| 嫩草影视91久久| 最好的美女福利视频网| 婷婷精品国产亚洲av| 日韩欧美在线乱码| 老师上课跳d突然被开到最大视频| 成熟少妇高潮喷水视频| 亚洲美女搞黄在线观看 | 国产精品久久久久久亚洲av鲁大| 99国产极品粉嫩在线观看| 午夜福利在线在线| 国产av在哪里看| 一a级毛片在线观看| 国产精品一及| 91麻豆精品激情在线观看国产| 色综合站精品国产| 人妻久久中文字幕网| 在线观看66精品国产| 欧美精品国产亚洲| 免费看av在线观看网站| 精品熟女少妇av免费看| 日本 av在线| 国产av不卡久久| 婷婷精品国产亚洲av在线| 亚洲精品在线观看二区| av国产免费在线观看| 亚洲国产欧洲综合997久久,| 精品免费久久久久久久清纯| 女人十人毛片免费观看3o分钟| 中文字幕精品亚洲无线码一区| 伦精品一区二区三区| 村上凉子中文字幕在线| 国产亚洲av嫩草精品影院| 一级毛片aaaaaa免费看小| 免费看a级黄色片| 亚洲av一区综合| 国产日本99.免费观看| 人妻久久中文字幕网| 性欧美人与动物交配| 成熟少妇高潮喷水视频| 国产色婷婷99| 国产成人一区二区在线| 色播亚洲综合网| 性色avwww在线观看| 成熟少妇高潮喷水视频| 久久鲁丝午夜福利片| 国产精品福利在线免费观看| 美女大奶头视频| 丰满人妻一区二区三区视频av| 国产黄a三级三级三级人| 一个人观看的视频www高清免费观看| 成人性生交大片免费视频hd| 在线播放国产精品三级| 免费一级毛片在线播放高清视频| 欧美区成人在线视频| 亚洲精品色激情综合| 简卡轻食公司| 精品国内亚洲2022精品成人| 在线a可以看的网站| 欧美区成人在线视频| 中文亚洲av片在线观看爽| 99久久久亚洲精品蜜臀av| 日本成人三级电影网站| 悠悠久久av| 国产成人91sexporn| 夜夜爽天天搞| 老司机福利观看| 亚洲性久久影院| 身体一侧抽搐| avwww免费| 嫩草影院精品99| 性色avwww在线观看| 校园春色视频在线观看| 国产精品野战在线观看| 亚洲专区国产一区二区| 九九久久精品国产亚洲av麻豆| 精品免费久久久久久久清纯| 国产真实乱freesex| 亚洲人成网站在线播| 日韩成人伦理影院| 高清午夜精品一区二区三区 | АⅤ资源中文在线天堂| 校园春色视频在线观看| 自拍偷自拍亚洲精品老妇| 日韩高清综合在线| 特级一级黄色大片| 十八禁网站免费在线| 国产v大片淫在线免费观看| 国产成人精品久久久久久| 欧美一区二区亚洲| 男女那种视频在线观看| 人人妻,人人澡人人爽秒播| 美女内射精品一级片tv| 熟妇人妻久久中文字幕3abv| 日韩av不卡免费在线播放| 久久精品综合一区二区三区| 岛国在线免费视频观看| 51国产日韩欧美| 成人二区视频| 一级毛片aaaaaa免费看小| 亚洲专区国产一区二区| 在线观看一区二区三区| 天美传媒精品一区二区| 国产视频内射| 欧美人与善性xxx| 日日啪夜夜撸| 乱人视频在线观看| 麻豆精品久久久久久蜜桃| 久久久欧美国产精品| 日韩av不卡免费在线播放| 色哟哟·www| 亚洲av免费高清在线观看| 亚洲av成人精品一区久久| 日韩精品中文字幕看吧| 亚洲精品亚洲一区二区| 亚洲国产欧美人成| 国产av一区在线观看免费| 国产爱豆传媒在线观看| 亚洲性夜色夜夜综合| 国产69精品久久久久777片| a级一级毛片免费在线观看| 大香蕉久久网| 精品一区二区三区人妻视频| 国产伦精品一区二区三区视频9| 搡老妇女老女人老熟妇| 精品无人区乱码1区二区| 午夜免费男女啪啪视频观看 | 淫秽高清视频在线观看| 夜夜夜夜夜久久久久| 人人妻人人澡人人爽人人夜夜 | 久久人妻av系列| av中文乱码字幕在线| 精品久久久久久久末码| 国产伦在线观看视频一区| 欧美精品国产亚洲| 丝袜喷水一区| 男女下面进入的视频免费午夜| 国产真实乱freesex| 国产精品乱码一区二三区的特点| 日韩强制内射视频| 欧美国产日韩亚洲一区| 欧美日韩综合久久久久久| 免费av毛片视频| 俄罗斯特黄特色一大片| 97热精品久久久久久| 97超级碰碰碰精品色视频在线观看| 淫妇啪啪啪对白视频| 亚洲成人av在线免费| 国产精品无大码| 69人妻影院| 国产精品一区二区免费欧美| 久久精品国产自在天天线| 亚洲七黄色美女视频| 男人狂女人下面高潮的视频| 淫妇啪啪啪对白视频| 国产毛片a区久久久久| 亚洲精品在线观看二区| 国产亚洲91精品色在线| 欧美色视频一区免费| 欧美潮喷喷水| 神马国产精品三级电影在线观看| 国产精品女同一区二区软件| 精品午夜福利视频在线观看一区| 欧美日韩国产亚洲二区| 美女大奶头视频| 18禁在线播放成人免费| 99热这里只有是精品在线观看| 亚洲,欧美,日韩| 国产精品一区二区三区四区久久| 少妇熟女aⅴ在线视频| 长腿黑丝高跟| 国产aⅴ精品一区二区三区波| 一级a爱片免费观看的视频| 国产精品美女特级片免费视频播放器| 亚洲第一区二区三区不卡| 一区二区三区四区激情视频 | 国产亚洲精品综合一区在线观看| 一区二区三区高清视频在线| 精品乱码久久久久久99久播| 卡戴珊不雅视频在线播放| 麻豆精品久久久久久蜜桃| 亚洲精华国产精华液的使用体验 | 国产亚洲精品综合一区在线观看| 国产三级在线视频| 在线观看美女被高潮喷水网站| 国产大屁股一区二区在线视频| 少妇人妻一区二区三区视频| 国内精品宾馆在线| 免费一级毛片在线播放高清视频| 99久久久亚洲精品蜜臀av| 哪里可以看免费的av片| 久久6这里有精品| 国产v大片淫在线免费观看| 天天一区二区日本电影三级| 国产 一区精品| 欧美性猛交╳xxx乱大交人| 成人二区视频| av天堂中文字幕网| 日日干狠狠操夜夜爽| 亚洲人成网站在线观看播放| 亚洲国产欧洲综合997久久,| 国产又黄又爽又无遮挡在线| 99久久中文字幕三级久久日本| 午夜激情欧美在线| 亚洲,欧美,日韩| 国产真实伦视频高清在线观看| 91久久精品国产一区二区三区| 国产成人影院久久av| av在线老鸭窝| 在线观看66精品国产| 亚洲av成人av| 久久久国产成人免费| 欧美潮喷喷水| 成人国产麻豆网| 中文字幕精品亚洲无线码一区| 精品一区二区三区视频在线观看免费| 国产精品久久视频播放| 熟妇人妻久久中文字幕3abv| 国产一区二区激情短视频| 如何舔出高潮| 99视频精品全部免费 在线| АⅤ资源中文在线天堂| 国产日本99.免费观看| 可以在线观看毛片的网站| 国产精品综合久久久久久久免费| 国产国拍精品亚洲av在线观看| 久99久视频精品免费| 搞女人的毛片| 色噜噜av男人的天堂激情| 欧美xxxx黑人xx丫x性爽| 国产伦精品一区二区三区四那| 免费电影在线观看免费观看| 精品久久久久久久久久久久久| 3wmmmm亚洲av在线观看| 一级黄片播放器| 成人av一区二区三区在线看| 亚洲国产色片| 又黄又爽又免费观看的视频| 美女黄网站色视频| 97超级碰碰碰精品色视频在线观看| 亚洲色图av天堂| 十八禁网站免费在线| 日本欧美国产在线视频| 啦啦啦韩国在线观看视频| 国产午夜福利久久久久久| 色综合站精品国产| 亚洲人成网站在线播| 在线天堂最新版资源| 天堂动漫精品| 色视频www国产| 女同久久另类99精品国产91| 亚洲七黄色美女视频| 久久久久久久久中文| 99久久精品热视频| 国产片特级美女逼逼视频| 精品久久久久久久人妻蜜臀av| 日韩精品青青久久久久久| 深爱激情五月婷婷| 性欧美人与动物交配| 国产高清不卡午夜福利| 九九热线精品视视频播放| 最近最新中文字幕大全电影3| 免费看a级黄色片| 六月丁香七月| 国产成人a区在线观看| 国产蜜桃级精品一区二区三区| 亚洲精品乱码久久久v下载方式| www日本黄色视频网| 精品久久久久久成人av| avwww免费| 国产精品永久免费网站| 日本精品一区二区三区蜜桃| 最近最新中文字幕大全电影3| 精品人妻视频免费看| 一本久久中文字幕| 欧美潮喷喷水| 午夜影院日韩av| 亚洲av熟女| 精品久久久噜噜| 国产午夜精品久久久久久一区二区三区 | 亚洲精品成人久久久久久| 国产亚洲精品久久久com| 99久久九九国产精品国产免费| 国产一区二区亚洲精品在线观看| 成人亚洲精品av一区二区| 99视频精品全部免费 在线| 有码 亚洲区| 性欧美人与动物交配| 国产极品精品免费视频能看的| a级毛片a级免费在线| 欧美成人精品欧美一级黄| 高清午夜精品一区二区三区 | 嫩草影院入口| 欧美日本视频| 97超视频在线观看视频| 国产单亲对白刺激| 在线a可以看的网站| 女生性感内裤真人,穿戴方法视频| 亚洲一区高清亚洲精品| 国产白丝娇喘喷水9色精品| 久久久久久久久久黄片| 桃色一区二区三区在线观看| 亚洲av二区三区四区| 一级毛片电影观看 | 日韩欧美三级三区| 在线看三级毛片| 国产毛片a区久久久久| 久久中文看片网| 三级毛片av免费| 少妇高潮的动态图| 国产男人的电影天堂91| 热99re8久久精品国产| 天天一区二区日本电影三级| 欧美精品国产亚洲| 国模一区二区三区四区视频| 欧美潮喷喷水| av国产免费在线观看| 老司机影院成人| 内射极品少妇av片p| 我的女老师完整版在线观看| 久99久视频精品免费| 国产成人一区二区在线| 亚洲综合色惰| 国产日本99.免费观看| 国产极品精品免费视频能看的| 国产人妻一区二区三区在| 国产单亲对白刺激| 亚洲第一区二区三区不卡| 天堂动漫精品| 亚洲一区高清亚洲精品| 人人妻人人澡人人爽人人夜夜 | 99在线人妻在线中文字幕| 日韩一区二区视频免费看| 激情 狠狠 欧美| 免费人成在线观看视频色| 香蕉av资源在线| 99热6这里只有精品| 亚洲经典国产精华液单| 国产精品一区二区三区四区久久| 久久久国产成人免费| 老熟妇乱子伦视频在线观看| 五月玫瑰六月丁香| 女的被弄到高潮叫床怎么办| 男人舔奶头视频| 一本久久中文字幕| 欧美高清成人免费视频www| 午夜a级毛片| 欧洲精品卡2卡3卡4卡5卡区| 高清日韩中文字幕在线| 免费观看精品视频网站| 美女xxoo啪啪120秒动态图| h日本视频在线播放| 俺也久久电影网| 身体一侧抽搐| 日本黄色片子视频| 亚洲激情五月婷婷啪啪| 麻豆成人午夜福利视频| 精品不卡国产一区二区三区| 日韩强制内射视频| 精品日产1卡2卡| 白带黄色成豆腐渣| 长腿黑丝高跟| 久久久国产成人精品二区| 日韩欧美 国产精品| 国产精品女同一区二区软件| 丝袜喷水一区| 97人妻精品一区二区三区麻豆| 久久精品国产99精品国产亚洲性色| 婷婷色综合大香蕉| 国产高清三级在线| 热99在线观看视频| 国产精品伦人一区二区| 天天躁夜夜躁狠狠久久av| 亚洲国产精品成人综合色| 联通29元200g的流量卡| 搡女人真爽免费视频火全软件 | 韩国av在线不卡| 精品一区二区三区人妻视频| 免费人成视频x8x8入口观看| 国产精品人妻久久久久久| 村上凉子中文字幕在线| 一本一本综合久久| 婷婷亚洲欧美| 久久久精品欧美日韩精品| 青春草视频在线免费观看| 亚洲av.av天堂| 国产精品久久电影中文字幕| 亚洲自偷自拍三级| 国产亚洲精品综合一区在线观看| 久久6这里有精品| 午夜福利在线观看免费完整高清在 | 好男人在线观看高清免费视频| 亚洲精品在线观看二区| 日本在线视频免费播放| www日本黄色视频网| 亚洲精品456在线播放app| 在线观看av片永久免费下载| 嫩草影院入口| 欧美+日韩+精品| 有码 亚洲区| 日本三级黄在线观看| 22中文网久久字幕| 久久国内精品自在自线图片| 黄色一级大片看看| 一夜夜www| 成人一区二区视频在线观看| 丰满的人妻完整版| 欧美另类亚洲清纯唯美| 黄片wwwwww| 最近最新中文字幕大全电影3| 色哟哟·www| 白带黄色成豆腐渣| 日韩 亚洲 欧美在线| avwww免费| 成人欧美大片| 99久久无色码亚洲精品果冻| 中文字幕av成人在线电影| 免费黄网站久久成人精品| 亚洲av成人av| 午夜免费男女啪啪视频观看 | 亚洲人成网站高清观看| 久久久久久久久大av| 欧美极品一区二区三区四区| 亚洲激情五月婷婷啪啪| 精品久久久久久久末码| 国产精品乱码一区二三区的特点| 国产精品精品国产色婷婷| 亚洲欧美日韩高清在线视频| 热99在线观看视频| 中出人妻视频一区二区| 亚洲真实伦在线观看| 亚洲人与动物交配视频| 亚洲欧美精品自产自拍| 性色avwww在线观看| 亚洲人成网站在线播放欧美日韩| 午夜免费激情av| 色哟哟哟哟哟哟| 国产淫片久久久久久久久| 我的老师免费观看完整版| 国产精品嫩草影院av在线观看| 午夜福利在线观看免费完整高清在 | 国产成人精品久久久久久| 卡戴珊不雅视频在线播放| 亚洲国产日韩欧美精品在线观看| 少妇的逼好多水| 精品一区二区三区av网在线观看| www.色视频.com| 在线看三级毛片| 亚洲av中文av极速乱| 美女cb高潮喷水在线观看| 日本成人三级电影网站| 美女内射精品一级片tv| 日韩人妻高清精品专区| 国产精品伦人一区二区| 亚洲美女黄片视频| 亚洲精品456在线播放app| 亚洲国产精品sss在线观看| 成年女人毛片免费观看观看9| 欧美日韩在线观看h| 99热6这里只有精品| 国产精品综合久久久久久久免费| 日韩亚洲欧美综合| 悠悠久久av| 麻豆国产av国片精品| 悠悠久久av| 熟女人妻精品中文字幕| ponron亚洲| 久久久久性生活片| 人人妻,人人澡人人爽秒播| 久久精品国产亚洲网站| 日韩制服骚丝袜av| 国产伦精品一区二区三区四那| 少妇猛男粗大的猛烈进出视频 | 国产乱人视频| 少妇熟女欧美另类| 男人舔奶头视频| 夜夜夜夜夜久久久久| 国产真实伦视频高清在线观看| 少妇熟女欧美另类| av在线观看视频网站免费| 麻豆乱淫一区二区| 深夜a级毛片| 色av中文字幕|