• <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
    别揉我奶头~嗯~啊~动态视频| 女生性感内裤真人,穿戴方法视频| 国产亚洲av高清不卡| www.熟女人妻精品国产| 别揉我奶头~嗯~啊~动态视频| 日日干狠狠操夜夜爽| 五月玫瑰六月丁香| 精品久久久久久,| 麻豆久久精品国产亚洲av| 亚洲av成人一区二区三| 精品久久久久久久久久久久久| 午夜激情欧美在线| 我要搜黄色片| 99精品欧美一区二区三区四区| 亚洲电影在线观看av| 成人永久免费在线观看视频| 精品一区二区三区av网在线观看| 久久草成人影院| 亚洲人与动物交配视频| 亚洲欧美日韩卡通动漫| 波多野结衣高清作品| 亚洲国产看品久久| 三级国产精品欧美在线观看 | 天天一区二区日本电影三级| xxx96com| 久久这里只有精品19| 久久天躁狠狠躁夜夜2o2o| 国产精品自产拍在线观看55亚洲| 国产高潮美女av| 夜夜躁狠狠躁天天躁| 99久久精品国产亚洲精品| 嫁个100分男人电影在线观看| 国产黄色小视频在线观看| 成年女人永久免费观看视频| 婷婷六月久久综合丁香| 女人高潮潮喷娇喘18禁视频| 麻豆久久精品国产亚洲av| 国产精品98久久久久久宅男小说| 亚洲国产精品合色在线| 色综合婷婷激情| 免费观看精品视频网站| av黄色大香蕉| 99热这里只有精品一区 | www.熟女人妻精品国产| 国产成人一区二区三区免费视频网站| 99riav亚洲国产免费| 亚洲欧美日韩卡通动漫| 一个人免费在线观看的高清视频| 在线永久观看黄色视频| 久久欧美精品欧美久久欧美| 日韩中文字幕欧美一区二区| 国产精品一区二区三区四区免费观看 | 青草久久国产| 欧美成人性av电影在线观看| 综合色av麻豆| 亚洲国产欧美人成| 亚洲精品美女久久久久99蜜臀| 午夜a级毛片| 精品久久久久久,| 波多野结衣高清无吗| 黑人巨大精品欧美一区二区mp4| 日韩精品中文字幕看吧| 午夜日韩欧美国产| 又大又爽又粗| 视频区欧美日本亚洲| 午夜免费观看网址| 久久久成人免费电影| 怎么达到女性高潮| 国产成人一区二区三区免费视频网站| 亚洲国产色片| 香蕉国产在线看| 国产高清有码在线观看视频| 一区福利在线观看| 国产精品久久久久久亚洲av鲁大| 色老头精品视频在线观看| av黄色大香蕉| 啪啪无遮挡十八禁网站| 亚洲九九香蕉| 亚洲狠狠婷婷综合久久图片| 麻豆国产97在线/欧美| 99热只有精品国产| 国产欧美日韩一区二区三| 欧美日本亚洲视频在线播放| 人人妻,人人澡人人爽秒播| 国产三级中文精品| 老熟妇仑乱视频hdxx| 国产一区二区三区视频了| 久久久久精品国产欧美久久久| 中出人妻视频一区二区| 桃色一区二区三区在线观看| 亚洲欧美日韩高清专用| 俺也久久电影网| 中亚洲国语对白在线视频| 国产精品 国内视频| 免费搜索国产男女视频| 亚洲av熟女| 91av网站免费观看| 黄频高清免费视频| 一夜夜www| 日本免费a在线| 狠狠狠狠99中文字幕| 国产精品av久久久久免费| 久久精品综合一区二区三区| 国产av不卡久久| 天堂动漫精品| АⅤ资源中文在线天堂| 天堂动漫精品| 成年女人永久免费观看视频| 一个人看的www免费观看视频| 人人妻人人看人人澡| 麻豆一二三区av精品| 日日夜夜操网爽| 99热精品在线国产| 制服人妻中文乱码| 国产毛片a区久久久久| 国产亚洲精品av在线| 精品久久久久久久人妻蜜臀av| 国产激情欧美一区二区| 国产探花在线观看一区二区| 美女被艹到高潮喷水动态| 亚洲 欧美一区二区三区| 国产激情偷乱视频一区二区| 一区二区三区国产精品乱码| 色综合婷婷激情| 久久精品国产综合久久久| 欧美乱色亚洲激情| 在线观看舔阴道视频| 亚洲午夜理论影院| 少妇的逼水好多| 色综合婷婷激情| 国产免费男女视频| 搞女人的毛片| 亚洲国产精品合色在线| 免费无遮挡裸体视频| 少妇熟女aⅴ在线视频| 白带黄色成豆腐渣| 国产精品99久久99久久久不卡| 三级男女做爰猛烈吃奶摸视频| 女同久久另类99精品国产91| 一二三四在线观看免费中文在| 国产精品98久久久久久宅男小说| 国产精品98久久久久久宅男小说| 日韩中文字幕欧美一区二区| 麻豆久久精品国产亚洲av| 欧美+亚洲+日韩+国产| 琪琪午夜伦伦电影理论片6080| 久久这里只有精品中国| 久久这里只有精品中国| 深夜精品福利| 成人特级av手机在线观看| av在线天堂中文字幕| 精品国产乱码久久久久久男人| 18禁美女被吸乳视频| 亚洲国产精品成人综合色| 丝袜人妻中文字幕| 毛片女人毛片| 精品一区二区三区视频在线 | 免费大片18禁| 免费人成视频x8x8入口观看| 国产成人精品久久二区二区91| 丝袜人妻中文字幕| 国产精品99久久99久久久不卡| 757午夜福利合集在线观看| 手机成人av网站| 欧美日韩综合久久久久久 | 免费在线观看亚洲国产| 老司机午夜十八禁免费视频| 国产精品亚洲av一区麻豆| 老司机午夜福利在线观看视频| 日韩欧美国产在线观看| 最近最新中文字幕大全电影3| 国产1区2区3区精品| 在线观看舔阴道视频| 亚洲va日本ⅴa欧美va伊人久久| 国产精品国产高清国产av| 在线观看免费视频日本深夜| 亚洲国产高清在线一区二区三| 欧美一区二区精品小视频在线| 久久国产乱子伦精品免费另类| 国产精品香港三级国产av潘金莲| 久久天堂一区二区三区四区| 99精品欧美一区二区三区四区| 国产成人福利小说| 99国产综合亚洲精品| 51午夜福利影视在线观看| 俺也久久电影网| 国产真实乱freesex| 久久这里只有精品19| 成在线人永久免费视频| 91在线精品国自产拍蜜月 | 欧美性猛交╳xxx乱大交人| 精品免费久久久久久久清纯| 一边摸一边抽搐一进一小说| 1024手机看黄色片| 免费在线观看亚洲国产| 欧美一级毛片孕妇| 国产aⅴ精品一区二区三区波| 免费观看精品视频网站| 国产伦一二天堂av在线观看| 小说图片视频综合网站| 黄色 视频免费看| 嫩草影视91久久| 亚洲成人中文字幕在线播放| 97超视频在线观看视频| 亚洲av日韩精品久久久久久密| 国产黄色小视频在线观看| 国产三级中文精品| 亚洲美女黄片视频| 久久天堂一区二区三区四区| 国产精品爽爽va在线观看网站| 日本a在线网址| 欧美成狂野欧美在线观看| 国产亚洲欧美98| 法律面前人人平等表现在哪些方面| 夜夜躁狠狠躁天天躁| 成人永久免费在线观看视频| 亚洲电影在线观看av| 成人三级黄色视频| 在线永久观看黄色视频| 亚洲精品美女久久av网站| 亚洲欧美日韩卡通动漫| 亚洲国产色片| 欧洲精品卡2卡3卡4卡5卡区| 国产精品电影一区二区三区| 白带黄色成豆腐渣| 大型黄色视频在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 国产成年人精品一区二区| 制服丝袜大香蕉在线| 好男人在线观看高清免费视频| 99久久综合精品五月天人人| 一级毛片女人18水好多| 国产精品一区二区三区四区免费观看 | 精品一区二区三区视频在线观看免费| 日韩欧美精品v在线| 在线免费观看的www视频| 一个人免费在线观看电影 | 成年版毛片免费区| 我的老师免费观看完整版| 国产不卡一卡二| 亚洲无线在线观看| 一边摸一边抽搐一进一小说| 狂野欧美白嫩少妇大欣赏| 99久久综合精品五月天人人| 国产欧美日韩精品一区二区| 免费看光身美女| av在线蜜桃| 国产主播在线观看一区二区| 嫩草影院精品99| 久久午夜亚洲精品久久| 成熟少妇高潮喷水视频| 精品久久久久久久人妻蜜臀av| 久久国产精品影院| 国产野战对白在线观看| 久久久久久人人人人人| 婷婷精品国产亚洲av在线| 黄色女人牲交| 观看美女的网站| 成人av在线播放网站| 蜜桃久久精品国产亚洲av| 无限看片的www在线观看| 免费在线观看影片大全网站| av国产免费在线观看| 亚洲精品美女久久久久99蜜臀| 国产精品 国内视频| 两人在一起打扑克的视频| 法律面前人人平等表现在哪些方面| 亚洲国产欧美网| 久久精品91蜜桃| 男女下面进入的视频免费午夜| 色综合站精品国产| 一本久久中文字幕| 精品99又大又爽又粗少妇毛片 | 日本免费一区二区三区高清不卡| 亚洲欧美日韩高清在线视频| 国产精品国产高清国产av| 国产精品影院久久| 美女高潮喷水抽搐中文字幕| cao死你这个sao货| 99久久精品国产亚洲精品| 久久久精品大字幕| 亚洲人成网站在线播放欧美日韩| 午夜免费成人在线视频| 一区二区三区高清视频在线| 国产精品久久久人人做人人爽| 叶爱在线成人免费视频播放| 后天国语完整版免费观看| 精品午夜福利视频在线观看一区| 看黄色毛片网站| 两个人的视频大全免费| 午夜免费观看网址| 岛国在线免费视频观看| 亚洲成人久久爱视频| 亚洲欧美日韩高清在线视频| 亚洲专区国产一区二区| 久久人人精品亚洲av| 成年女人永久免费观看视频| 国产亚洲av嫩草精品影院| 成人三级做爰电影| 国产成+人综合+亚洲专区| 亚洲一区二区三区色噜噜| 久久婷婷人人爽人人干人人爱| or卡值多少钱| 男女床上黄色一级片免费看| 欧美zozozo另类| 狠狠狠狠99中文字幕| 亚洲成av人片免费观看| 亚洲欧洲精品一区二区精品久久久| 天堂av国产一区二区熟女人妻| 国产一区二区三区视频了| 久久久成人免费电影| 久久精品91无色码中文字幕| 桃色一区二区三区在线观看| 免费一级毛片在线播放高清视频| 亚洲电影在线观看av| 97碰自拍视频| 久久精品夜夜夜夜夜久久蜜豆| 1024香蕉在线观看| 天天躁日日操中文字幕| 免费看日本二区| 搡老熟女国产l中国老女人| av在线蜜桃| 99热只有精品国产| 嫁个100分男人电影在线观看| 亚洲男人的天堂狠狠| 成人国产综合亚洲| 亚洲国产精品999在线| 亚洲中文av在线| 欧美日韩一级在线毛片| 亚洲在线自拍视频| 亚洲 欧美 日韩 在线 免费| 精品国产美女av久久久久小说| 啪啪无遮挡十八禁网站| 91麻豆av在线| 亚洲一区二区三区不卡视频| 国产精品影院久久| 久久久久久大精品| 欧美色欧美亚洲另类二区| 啦啦啦韩国在线观看视频| 动漫黄色视频在线观看| АⅤ资源中文在线天堂| 精品欧美国产一区二区三| 在线播放国产精品三级| 国产亚洲av高清不卡| 在线免费观看的www视频| 又大又爽又粗| 国产高清videossex| h日本视频在线播放| avwww免费| 啪啪无遮挡十八禁网站| 精品99又大又爽又粗少妇毛片 | 熟女人妻精品中文字幕| 亚洲av成人精品一区久久| 中文亚洲av片在线观看爽| 国产蜜桃级精品一区二区三区| 中国美女看黄片| 久久久久精品国产欧美久久久| 欧美一级a爱片免费观看看| 欧美一级a爱片免费观看看| 日韩欧美在线乱码| 亚洲av成人不卡在线观看播放网| 久久久久久国产a免费观看| 熟女人妻精品中文字幕| 国产欧美日韩精品一区二区| 99久久久亚洲精品蜜臀av| АⅤ资源中文在线天堂| bbb黄色大片| 精品国内亚洲2022精品成人| 国内精品美女久久久久久| 亚洲熟妇中文字幕五十中出| 午夜福利视频1000在线观看| 啦啦啦观看免费观看视频高清| 亚洲精品在线美女| 国产成年人精品一区二区| 一本一本综合久久| 亚洲在线观看片| 国内精品美女久久久久久| 十八禁网站免费在线| 熟女人妻精品中文字幕| 99在线人妻在线中文字幕| 男女之事视频高清在线观看| 国产视频内射| 欧美激情久久久久久爽电影| 午夜a级毛片| 成人无遮挡网站| 制服人妻中文乱码| 色哟哟哟哟哟哟| 琪琪午夜伦伦电影理论片6080| 亚洲人成电影免费在线| 99久久精品一区二区三区| 人人妻,人人澡人人爽秒播| 99国产综合亚洲精品| 欧美中文日本在线观看视频| 欧美不卡视频在线免费观看| 欧美午夜高清在线| 久久性视频一级片| 99精品在免费线老司机午夜| 男人舔奶头视频| 一夜夜www| 美女扒开内裤让男人捅视频| 成年版毛片免费区| 久久这里只有精品中国| 俺也久久电影网| 国产探花在线观看一区二区| 国产高潮美女av| 中文字幕熟女人妻在线| 国产精品亚洲av一区麻豆| 757午夜福利合集在线观看| 国产精品久久久人人做人人爽| 亚洲欧美日韩高清专用| 丁香六月欧美| 99国产极品粉嫩在线观看| 我要搜黄色片| 国内精品久久久久久久电影| 一本综合久久免费| 91av网一区二区| 亚洲黑人精品在线| 久久人妻av系列| 亚洲专区中文字幕在线| 精品国产乱码久久久久久男人| 色综合欧美亚洲国产小说| 国产精品1区2区在线观看.| 国产又色又爽无遮挡免费看| 老鸭窝网址在线观看| 两个人看的免费小视频| 亚洲成人免费电影在线观看| 亚洲av免费在线观看| 一进一出抽搐动态| 香蕉久久夜色| 草草在线视频免费看| www.www免费av| 欧美最黄视频在线播放免费| 国内少妇人妻偷人精品xxx网站 | 亚洲国产中文字幕在线视频| 99视频精品全部免费 在线 | 法律面前人人平等表现在哪些方面| 国产视频一区二区在线看| 免费高清视频大片| 国产成人av教育| 欧美乱色亚洲激情| 国产免费av片在线观看野外av| 欧美成狂野欧美在线观看| 亚洲av中文字字幕乱码综合| 人人妻人人澡欧美一区二区| 少妇裸体淫交视频免费看高清| 日韩欧美在线二视频| 动漫黄色视频在线观看| 怎么达到女性高潮| 成年版毛片免费区| 国产精品,欧美在线| 三级毛片av免费| 国产精品乱码一区二三区的特点| 亚洲国产精品合色在线| 久久久成人免费电影| 精品一区二区三区视频在线 | 国产人伦9x9x在线观看| 久久久久久人人人人人| 日韩欧美免费精品| 熟女少妇亚洲综合色aaa.| 国产亚洲精品久久久久久毛片| 18禁国产床啪视频网站| 亚洲人成电影免费在线| 99国产精品一区二区三区| 18禁国产床啪视频网站| 午夜激情欧美在线| 免费无遮挡裸体视频| 宅男免费午夜| 亚洲欧美激情综合另类| 99视频精品全部免费 在线 | 午夜成年电影在线免费观看| 欧美日韩黄片免| 中亚洲国语对白在线视频| 99精品在免费线老司机午夜| 一级a爱片免费观看的视频| 色在线成人网| a在线观看视频网站| 国产激情偷乱视频一区二区| 久久精品夜夜夜夜夜久久蜜豆| 久久久色成人| 757午夜福利合集在线观看| 波多野结衣高清作品| 又爽又黄无遮挡网站| 亚洲中文字幕一区二区三区有码在线看 | 欧美乱妇无乱码| 久久精品亚洲精品国产色婷小说| 免费电影在线观看免费观看| 国产亚洲欧美98| 精品国产亚洲在线| 狠狠狠狠99中文字幕| 夜夜躁狠狠躁天天躁| 一级a爱片免费观看的视频| 午夜精品在线福利| 搞女人的毛片| 亚洲精品美女久久久久99蜜臀| 美女高潮的动态| 人妻夜夜爽99麻豆av| 国产成人精品无人区| 99精品在免费线老司机午夜| 91av网一区二区| 1024手机看黄色片| 国产精品乱码一区二三区的特点| 亚洲国产欧洲综合997久久,| 成年女人毛片免费观看观看9| 成年女人看的毛片在线观看| 亚洲av熟女| 中文字幕高清在线视频| 他把我摸到了高潮在线观看| bbb黄色大片| 日韩中文字幕欧美一区二区| 久久久精品欧美日韩精品| 日本一本二区三区精品| 亚洲精品一区av在线观看| 国产亚洲精品久久久com| 欧美乱色亚洲激情| 国产精品九九99| 黄色丝袜av网址大全| 午夜福利18| 青草久久国产| 在线十欧美十亚洲十日本专区| 国产一区二区三区在线臀色熟女| 亚洲欧洲精品一区二区精品久久久| 久久国产乱子伦精品免费另类| 久久久色成人| 国产精品一区二区精品视频观看| 久久久精品欧美日韩精品| 国产aⅴ精品一区二区三区波| 成人三级黄色视频| 国产麻豆成人av免费视频| 九九在线视频观看精品| 黄频高清免费视频| 99久久精品一区二区三区| 国产精品久久久人人做人人爽| 成年女人毛片免费观看观看9| 91字幕亚洲| 日韩高清综合在线| 又黄又爽又免费观看的视频| 天堂动漫精品| 亚洲国产中文字幕在线视频| 国产av在哪里看| 给我免费播放毛片高清在线观看| 成人永久免费在线观看视频| 日韩三级视频一区二区三区| 黄色丝袜av网址大全| 国产精品香港三级国产av潘金莲| 在线免费观看不下载黄p国产 | 国产亚洲精品av在线| 亚洲熟妇熟女久久| 精品久久久久久成人av| 亚洲精品一区av在线观看| 国产欧美日韩一区二区精品| 国产精品99久久99久久久不卡| 狠狠狠狠99中文字幕| 欧美高清成人免费视频www| 久久天堂一区二区三区四区| 老司机福利观看| 日韩欧美精品v在线| 国产成人一区二区三区免费视频网站| av天堂中文字幕网| 日韩欧美免费精品| 免费一级毛片在线播放高清视频| 亚洲av美国av| 精品久久久久久久久久久久久| 少妇人妻一区二区三区视频| 国产欧美日韩一区二区精品| 亚洲成人免费电影在线观看| 嫩草影院入口| 超碰成人久久| 国产午夜精品久久久久久| 真实男女啪啪啪动态图| 精品国产三级普通话版| 国产成人欧美在线观看| 亚洲真实伦在线观看| 可以在线观看的亚洲视频| 午夜两性在线视频| 亚洲精华国产精华精| 色老头精品视频在线观看| 香蕉丝袜av| 欧美乱码精品一区二区三区| 麻豆国产97在线/欧美| 91老司机精品| 88av欧美| 国产又黄又爽又无遮挡在线| 国产精品av视频在线免费观看| 精品久久蜜臀av无| 神马国产精品三级电影在线观看| 成人性生交大片免费视频hd| 亚洲国产高清在线一区二区三| 欧美xxxx黑人xx丫x性爽| h日本视频在线播放| 精品免费久久久久久久清纯| 国产伦一二天堂av在线观看| 视频区欧美日本亚洲| 国产激情欧美一区二区| 免费高清视频大片| 精品久久久久久久久久免费视频| 他把我摸到了高潮在线观看| www.熟女人妻精品国产| 国产蜜桃级精品一区二区三区| 欧美国产日韩亚洲一区| 国产高清激情床上av| 国产激情久久老熟女| 日本 av在线| 欧美一级毛片孕妇| 亚洲自拍偷在线| 夜夜夜夜夜久久久久| 久久人人精品亚洲av| 亚洲一区高清亚洲精品| 日韩成人在线观看一区二区三区| 桃色一区二区三区在线观看| 久久国产乱子伦精品免费另类| 国产69精品久久久久777片 | 久久久久国内视频| 久久国产精品影院| 一边摸一边抽搐一进一小说| 婷婷精品国产亚洲av|