• <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
    99久久久亚洲精品蜜臀av| 久久久国产精品麻豆| 免费看美女性在线毛片视频| 久久久国产成人免费| 99久国产av精品| 夜夜躁狠狠躁天天躁| 日本熟妇午夜| 老司机深夜福利视频在线观看| 精品国内亚洲2022精品成人| 欧美另类亚洲清纯唯美| 国产激情久久老熟女| 高清毛片免费观看视频网站| 精品久久久久久久末码| 美女cb高潮喷水在线观看 | 真人一进一出gif抽搐免费| 亚洲精品国产精品久久久不卡| 欧美日韩综合久久久久久 | 成人无遮挡网站| 偷拍熟女少妇极品色| 免费电影在线观看免费观看| 美女免费视频网站| 欧美一区二区精品小视频在线| 欧美最黄视频在线播放免费| 天天躁日日操中文字幕| 国产三级在线视频| 神马国产精品三级电影在线观看| 亚洲 国产 在线| 精品国产超薄肉色丝袜足j| 手机成人av网站| 老司机深夜福利视频在线观看| netflix在线观看网站| 日日摸夜夜添夜夜添小说| 久久久久久久久免费视频了| 啦啦啦韩国在线观看视频| 亚洲第一电影网av| 男女那种视频在线观看| 国产激情偷乱视频一区二区| 人人妻人人看人人澡| 国产探花在线观看一区二区| 国产1区2区3区精品| 两个人看的免费小视频| 一级毛片精品| 中文字幕高清在线视频| 亚洲专区字幕在线| 成年版毛片免费区| 五月伊人婷婷丁香| 国产高清有码在线观看视频| 成年女人毛片免费观看观看9| 99热精品在线国产| 亚洲成人久久性| 免费av毛片视频| 51午夜福利影视在线观看| 午夜激情欧美在线| 99久久精品一区二区三区| 亚洲av五月六月丁香网| 97人妻精品一区二区三区麻豆| 亚洲av成人av| 久9热在线精品视频| www.熟女人妻精品国产| 国产精品98久久久久久宅男小说| 久久精品夜夜夜夜夜久久蜜豆| 久久中文字幕人妻熟女| 无人区码免费观看不卡| 9191精品国产免费久久| 91字幕亚洲| 国产成年人精品一区二区| 国产精品 国内视频| e午夜精品久久久久久久| 国产精品电影一区二区三区| 国产精品综合久久久久久久免费| АⅤ资源中文在线天堂| 一级作爱视频免费观看| 真人做人爱边吃奶动态| 亚洲av成人不卡在线观看播放网| 舔av片在线| 无人区码免费观看不卡| 欧美高清成人免费视频www| 久久久久久人人人人人| 19禁男女啪啪无遮挡网站| 午夜精品在线福利| 国产又黄又爽又无遮挡在线| 黄色成人免费大全| 久久精品国产清高在天天线| 国模一区二区三区四区视频 | 免费av毛片视频| 免费人成视频x8x8入口观看| 欧美黄色淫秽网站| 欧美乱码精品一区二区三区| 成人国产综合亚洲| 精品乱码久久久久久99久播| 五月玫瑰六月丁香| 香蕉av资源在线| 国产黄a三级三级三级人| 天天添夜夜摸| 18禁黄网站禁片午夜丰满| 免费观看精品视频网站| 黄色 视频免费看| 岛国视频午夜一区免费看| 久久精品91无色码中文字幕| 日本在线视频免费播放| 国产亚洲av高清不卡| 国产不卡一卡二| 黄色视频,在线免费观看| 可以在线观看的亚洲视频| 精品国产乱子伦一区二区三区| 日本成人三级电影网站| 亚洲精品色激情综合| 又大又爽又粗| 男女那种视频在线观看| 日韩免费av在线播放| 长腿黑丝高跟| 在线免费观看不下载黄p国产 | 精品国产亚洲在线| 男女午夜视频在线观看| 美女高潮喷水抽搐中文字幕| 欧美一区二区国产精品久久精品| 99国产精品一区二区蜜桃av| 最近在线观看免费完整版| 老鸭窝网址在线观看| 人人妻人人澡欧美一区二区| 变态另类丝袜制服| 久久久久国产精品人妻aⅴ院| 狂野欧美激情性xxxx| 久久久久国内视频| 午夜免费成人在线视频| 国产 一区 欧美 日韩| 欧美成人一区二区免费高清观看 | 亚洲avbb在线观看| 亚洲天堂国产精品一区在线| 欧美高清成人免费视频www| 在线十欧美十亚洲十日本专区| 亚洲精品在线美女| 99riav亚洲国产免费| 青草久久国产| 成年人黄色毛片网站| 一边摸一边抽搐一进一小说| 亚洲欧美日韩高清专用| 香蕉国产在线看| 国产成人一区二区三区免费视频网站| av片东京热男人的天堂| 国产v大片淫在线免费观看| 亚洲天堂国产精品一区在线| 91在线观看av| 亚洲人成网站在线播放欧美日韩| 久久久国产欧美日韩av| 91麻豆精品激情在线观看国产| 欧美日韩中文字幕国产精品一区二区三区| 成人特级av手机在线观看| 久久久国产成人免费| 午夜成年电影在线免费观看| 欧美乱色亚洲激情| 久久久久亚洲av毛片大全| 国产三级黄色录像| 亚洲成人久久爱视频| 啦啦啦韩国在线观看视频| 中文亚洲av片在线观看爽| 搡老岳熟女国产| 波多野结衣高清作品| 90打野战视频偷拍视频| 狂野欧美白嫩少妇大欣赏| 久久午夜综合久久蜜桃| 岛国在线观看网站| 此物有八面人人有两片| 九色成人免费人妻av| av女优亚洲男人天堂 | 在线a可以看的网站| 男女下面进入的视频免费午夜| 性色av乱码一区二区三区2| 美女cb高潮喷水在线观看 | 亚洲美女视频黄频| 午夜福利在线观看吧| 久久久久国产一级毛片高清牌| 高清在线国产一区| 三级男女做爰猛烈吃奶摸视频| 国产av在哪里看| 亚洲,欧美精品.| 国产亚洲av嫩草精品影院| 久久精品国产99精品国产亚洲性色| 精品久久久久久久久久免费视频| 中文亚洲av片在线观看爽| 午夜影院日韩av| 亚洲av第一区精品v没综合| 婷婷丁香在线五月| 成人国产综合亚洲| 999久久久国产精品视频| 亚洲av美国av| 亚洲国产色片| 国产人伦9x9x在线观看| 在线a可以看的网站| 一a级毛片在线观看| a级毛片a级免费在线| 美女 人体艺术 gogo| 久久久久久久久中文| 两个人的视频大全免费| 免费av毛片视频| 日本与韩国留学比较| 亚洲中文av在线| 久久久久性生活片| 国产亚洲精品综合一区在线观看| 国产精品亚洲一级av第二区| 中文字幕人成人乱码亚洲影| 成人三级黄色视频| 国产蜜桃级精品一区二区三区| 观看美女的网站| 国产成人一区二区三区免费视频网站| 中文字幕久久专区| 白带黄色成豆腐渣| 精品福利观看| 看免费av毛片| 日韩欧美在线乱码| 亚洲欧美精品综合一区二区三区| 成人国产一区最新在线观看| 国产精品久久久久久久电影 | 国内精品一区二区在线观看| 草草在线视频免费看| 久久天躁狠狠躁夜夜2o2o| 窝窝影院91人妻| 婷婷精品国产亚洲av在线| 视频区欧美日本亚洲| 亚洲成a人片在线一区二区| 午夜久久久久精精品| 变态另类成人亚洲欧美熟女| 88av欧美| 国产成人aa在线观看| 91在线观看av| 国产极品精品免费视频能看的| 免费看美女性在线毛片视频| 99热这里只有精品一区 | 久久精品国产清高在天天线| 国产精品久久久久久人妻精品电影| 高清在线国产一区| 91字幕亚洲| 久久精品综合一区二区三区| 99热这里只有是精品50| 久久久久久人人人人人| 我要搜黄色片| 久久久久精品国产欧美久久久| 国产男靠女视频免费网站| 日本免费a在线| 午夜成年电影在线免费观看| 亚洲专区国产一区二区| 国产av不卡久久| 无人区码免费观看不卡| 熟妇人妻久久中文字幕3abv| 日韩成人在线观看一区二区三区| 一个人观看的视频www高清免费观看 | 欧美黄色淫秽网站| 欧美日韩瑟瑟在线播放| 久久中文看片网| 国产精品日韩av在线免费观看| 美女 人体艺术 gogo| 99热这里只有精品一区 | 久久性视频一级片| 99国产精品99久久久久| 国产三级在线视频| 一本久久中文字幕| 欧美绝顶高潮抽搐喷水| 一本精品99久久精品77| 悠悠久久av| 精华霜和精华液先用哪个| 久久精品国产亚洲av香蕉五月| 久久久久国内视频| 欧美黄色淫秽网站| h日本视频在线播放| 久久天堂一区二区三区四区| 欧美中文综合在线视频| 亚洲电影在线观看av| 国产成人精品久久二区二区91| 国产又黄又爽又无遮挡在线| 国产精品av视频在线免费观看| 一本一本综合久久| 亚洲人成网站在线播放欧美日韩| 国产一区二区在线av高清观看| 日本三级黄在线观看| 久久午夜综合久久蜜桃| 亚洲专区中文字幕在线| 日韩欧美免费精品| 黄色丝袜av网址大全| 精品久久久久久久人妻蜜臀av| av天堂在线播放| 国产精品女同一区二区软件 | 免费在线观看亚洲国产| 亚洲色图av天堂| 91麻豆精品激情在线观看国产| 久久久久国产一级毛片高清牌| 精品福利观看| 久久国产精品影院| 亚洲欧美一区二区三区黑人| 国产精品永久免费网站| 少妇熟女aⅴ在线视频| 亚洲激情在线av| 真人一进一出gif抽搐免费| 制服人妻中文乱码| 波多野结衣高清无吗| 免费观看的影片在线观看| 搡老熟女国产l中国老女人| 十八禁网站免费在线| xxxwww97欧美| 久久国产精品影院| or卡值多少钱| xxx96com| 99久久精品热视频| 亚洲成人久久爱视频| 美女扒开内裤让男人捅视频| 色哟哟哟哟哟哟| 亚洲精品国产精品久久久不卡| 婷婷丁香在线五月| 色哟哟哟哟哟哟| 久久精品91无色码中文字幕| 啦啦啦韩国在线观看视频| 国产精品免费一区二区三区在线| 国内精品久久久久久久电影| 黑人操中国人逼视频| www.999成人在线观看| 性欧美人与动物交配| 亚洲一区二区三区色噜噜| 国产伦人伦偷精品视频| 久久精品国产99精品国产亚洲性色| 婷婷丁香在线五月| 悠悠久久av| 日韩成人在线观看一区二区三区| 男女床上黄色一级片免费看| or卡值多少钱| 精品久久久久久成人av| 国产精品久久电影中文字幕| 精品一区二区三区av网在线观看| 国产激情欧美一区二区| 免费无遮挡裸体视频| av中文乱码字幕在线| 淫秽高清视频在线观看| 国产成人av教育| 夜夜躁狠狠躁天天躁| 婷婷精品国产亚洲av| 色噜噜av男人的天堂激情| 国产精品精品国产色婷婷| 熟女人妻精品中文字幕| 国产aⅴ精品一区二区三区波| 不卡av一区二区三区| 三级国产精品欧美在线观看 | 岛国视频午夜一区免费看| 亚洲黑人精品在线| 免费在线观看日本一区| 中文字幕人成人乱码亚洲影| 精品免费久久久久久久清纯| 国产av麻豆久久久久久久| 亚洲美女黄片视频| 欧美色视频一区免费| 亚洲成av人片在线播放无| 精品国产三级普通话版| 日韩欧美在线二视频| 日韩欧美 国产精品| 国产一区二区在线av高清观看| 99国产综合亚洲精品| 人人妻人人看人人澡| 国产精品九九99| 久久天堂一区二区三区四区| av视频在线观看入口| 中文字幕精品亚洲无线码一区| 91麻豆av在线| 亚洲国产欧美网| 欧美一区二区精品小视频在线| 免费av毛片视频| 国产高清有码在线观看视频| 18禁黄网站禁片午夜丰满| 亚洲av成人精品一区久久| 国产视频一区二区在线看| 亚洲激情在线av| 久久精品91蜜桃| 欧美性猛交黑人性爽| 黄色片一级片一级黄色片| 欧美日韩精品网址| 岛国视频午夜一区免费看| 欧美中文综合在线视频| 国产野战对白在线观看| 日本一二三区视频观看| 免费看光身美女| 99久久精品热视频| 日本撒尿小便嘘嘘汇集6| 国产综合懂色| 欧美又色又爽又黄视频| 又大又爽又粗| 一个人看视频在线观看www免费 | 久久久久久人人人人人| 美女cb高潮喷水在线观看 | 成在线人永久免费视频| a在线观看视频网站| 国产免费男女视频| 亚洲精品国产精品久久久不卡| 校园春色视频在线观看| 又黄又爽又免费观看的视频| 夜夜爽天天搞| 搡老熟女国产l中国老女人| 国产私拍福利视频在线观看| 国产高清三级在线| 亚洲国产看品久久| 麻豆国产av国片精品| 亚洲va日本ⅴa欧美va伊人久久| 国产精品,欧美在线| 欧美3d第一页| 此物有八面人人有两片| a在线观看视频网站| 99国产精品一区二区三区| 国产精品,欧美在线| 久久久久国产一级毛片高清牌| 午夜福利在线观看免费完整高清在 | 亚洲人成网站高清观看| 身体一侧抽搐| 色播亚洲综合网| 日日摸夜夜添夜夜添小说| 老司机在亚洲福利影院| www.www免费av| 天堂av国产一区二区熟女人妻| 国产伦精品一区二区三区视频9 | 在线播放国产精品三级| 精品人妻1区二区| 丝袜人妻中文字幕| 亚洲国产欧美网| 久久中文字幕人妻熟女| 国产在线精品亚洲第一网站| 亚洲欧美日韩无卡精品| 在线免费观看不下载黄p国产 | 一进一出抽搐gif免费好疼| 欧美日韩福利视频一区二区| 国产视频一区二区在线看| 极品教师在线免费播放| 欧美一级a爱片免费观看看| 国产高清视频在线播放一区| 久久久久性生活片| 999久久久国产精品视频| 搡老妇女老女人老熟妇| 熟女少妇亚洲综合色aaa.| 国产精品久久久久久亚洲av鲁大| 国产精品 国内视频| 好看av亚洲va欧美ⅴa在| 天堂动漫精品| 精品99又大又爽又粗少妇毛片 | 嫩草影院精品99| 久久人妻av系列| 欧美在线黄色| 国产美女午夜福利| 一级黄色大片毛片| 亚洲 欧美一区二区三区| 国产男靠女视频免费网站| 午夜免费观看网址| 99国产极品粉嫩在线观看| 亚洲,欧美精品.| 麻豆成人午夜福利视频| 国产久久久一区二区三区| 亚洲色图av天堂| 毛片女人毛片| 日本一二三区视频观看| 窝窝影院91人妻| 老汉色∧v一级毛片| 精品久久久久久久毛片微露脸| 淫妇啪啪啪对白视频| 2021天堂中文幕一二区在线观| 男人和女人高潮做爰伦理| 精品福利观看| cao死你这个sao货| 老司机深夜福利视频在线观看| 国产精品日韩av在线免费观看| 蜜桃久久精品国产亚洲av| 18禁国产床啪视频网站| 伦理电影免费视频| 啦啦啦免费观看视频1| 亚洲人成伊人成综合网2020| 欧美一区二区精品小视频在线| 久久久久性生活片| 国产乱人视频| 天天添夜夜摸| 夜夜躁狠狠躁天天躁| 91麻豆精品激情在线观看国产| svipshipincom国产片| 亚洲专区字幕在线| 噜噜噜噜噜久久久久久91| 观看免费一级毛片| 久久中文看片网| 一级作爱视频免费观看| 久久久精品大字幕| 一区福利在线观看| 欧美黄色淫秽网站| 国产伦在线观看视频一区| 一个人观看的视频www高清免费观看 | 欧美一区二区国产精品久久精品| 国产精品一区二区免费欧美| 国产精品av视频在线免费观看| 免费看日本二区| 亚洲国产精品999在线| 久久伊人香网站| 国产成人系列免费观看| 日本黄色片子视频| 久久久久亚洲av毛片大全| 后天国语完整版免费观看| 日本五十路高清| 亚洲人成电影免费在线| 三级国产精品欧美在线观看 | 99久国产av精品| 亚洲国产欧美网| 99国产精品一区二区蜜桃av| 日韩精品中文字幕看吧| 在线观看美女被高潮喷水网站 | 中文字幕久久专区| 久久久精品欧美日韩精品| 啦啦啦韩国在线观看视频| 十八禁网站免费在线| 亚洲专区字幕在线| 久久久久久大精品| 国产激情偷乱视频一区二区| 日日干狠狠操夜夜爽| 成人高潮视频无遮挡免费网站| 久久这里只有精品中国| 一卡2卡三卡四卡精品乱码亚洲| 五月伊人婷婷丁香| 国产成人一区二区三区免费视频网站| 亚洲真实伦在线观看| 91av网站免费观看| 男人的好看免费观看在线视频| 午夜两性在线视频| 欧美xxxx黑人xx丫x性爽| 午夜福利成人在线免费观看| 国产成人福利小说| 天堂网av新在线| 日本免费a在线| 黑人操中国人逼视频| 免费观看人在逋| 成年女人看的毛片在线观看| 国产精华一区二区三区| 日本免费一区二区三区高清不卡| 99精品久久久久人妻精品| 九九在线视频观看精品| 黄色视频,在线免费观看| 欧美高清成人免费视频www| 香蕉丝袜av| 国产伦在线观看视频一区| 两个人视频免费观看高清| 免费看光身美女| 亚洲欧美日韩无卡精品| 精品无人区乱码1区二区| 国语自产精品视频在线第100页| 极品教师在线免费播放| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲美女视频黄频| 在线国产一区二区在线| 九色成人免费人妻av| 脱女人内裤的视频| 一级黄色大片毛片| 亚洲国产精品成人综合色| 欧美精品啪啪一区二区三区| 最近视频中文字幕2019在线8| 熟女少妇亚洲综合色aaa.| av国产免费在线观看| 亚洲成a人片在线一区二区| 黄色成人免费大全| 亚洲aⅴ乱码一区二区在线播放| 日韩国内少妇激情av| 国产成人精品无人区| 久久国产精品人妻蜜桃| 欧美黄色片欧美黄色片| 熟妇人妻久久中文字幕3abv| 久久久久久久精品吃奶| 精品久久久久久久久久免费视频| а√天堂www在线а√下载| www日本在线高清视频| 宅男免费午夜| 可以在线观看的亚洲视频| 国产免费av片在线观看野外av| 特大巨黑吊av在线直播| 欧美日韩黄片免| 操出白浆在线播放| 999精品在线视频| 午夜精品在线福利| 我的老师免费观看完整版| 九九热线精品视视频播放| 99热这里只有精品一区 | 看免费av毛片| 亚洲av第一区精品v没综合| 香蕉久久夜色| 久久中文字幕一级| 麻豆成人av在线观看| 18美女黄网站色大片免费观看| 亚洲av五月六月丁香网| 亚洲黑人精品在线| 最近最新中文字幕大全电影3| 亚洲av成人精品一区久久| 欧美黄色片欧美黄色片| 久久午夜综合久久蜜桃| 国产一区二区三区在线臀色熟女| 亚洲国产欧美网| 听说在线观看完整版免费高清| 亚洲国产中文字幕在线视频| 色吧在线观看| 成人欧美大片| 亚洲,欧美精品.| 午夜久久久久精精品| 成人国产一区最新在线观看| 噜噜噜噜噜久久久久久91| 国产伦一二天堂av在线观看| 亚洲七黄色美女视频| 精品一区二区三区av网在线观看| 国产成人福利小说| 可以在线观看毛片的网站| 亚洲中文字幕日韩| 香蕉国产在线看| 搞女人的毛片| 亚洲欧美日韩无卡精品| 一本久久中文字幕| 国产乱人伦免费视频| 国产伦人伦偷精品视频| 国产 一区 欧美 日韩| 亚洲片人在线观看| 特大巨黑吊av在线直播| 又紧又爽又黄一区二区| 五月玫瑰六月丁香| 亚洲精品乱码久久久v下载方式 | 久久午夜亚洲精品久久|