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

    楸樹葉綠體基因組密碼子偏性分析

    2021-02-09 00:32:08李鳳辛靜辛雅萱肖遙屈亞亞王軍輝麻文俊辛培堯
    關(guān)鍵詞:楸樹

    李鳳 辛靜 辛雅萱 肖遙 屈亞亞 王軍輝 麻文俊 辛培堯

    摘要:【目的】分析楸樹葉綠體基因組密碼子使用偏性,為楸樹乃至梓屬植物葉綠體基因組學(xué)、物種進(jìn)化特征及遺傳多樣性分析等研究提供理論依據(jù)?!痉椒ā坷肅odon W 1.4.2及EMBOSS expiorer的cusp和chips軟件對(duì)51個(gè)楸樹葉綠體基因的編碼區(qū)(CDS)序列進(jìn)行密碼子相關(guān)參數(shù)、中性繪圖、ENC-plot和PR2-plot等分析,并通過(guò)構(gòu)建高、低表達(dá)偏性基因庫(kù),以ENC值和RSCU值為參考標(biāo)準(zhǔn),最終篩選最優(yōu)密碼子。【結(jié)果】51個(gè)基因密碼子的平均GCall含量為38.48%,密碼子不同位置GC含量為:GC1(47.03%)>GC2(39.93%)>GC3(28.09%),說(shuō)明第3位的GC并非均勻分布;ENC為34.93~56.50,平均為46.87。GC2與GC1呈極顯著正相關(guān)(P<0.01),而GC3與GC1和GC2無(wú)顯著正相關(guān)(P>0.05);ENC與GCall、GC1、GC3和N均呈顯著正相關(guān)(P<0.05,下同);密碼子數(shù)(N)與GC3和ENC呈顯著正相關(guān)。RSCU>1的密碼子共有30個(gè),其中,以A或U結(jié)尾的密碼子有29個(gè)。中性繪圖分析結(jié)果顯示,GC3與GC12相關(guān)系數(shù)為0.0188,二者無(wú)顯著相關(guān)性,回歸系數(shù)為0.143。ENC-plot分析結(jié)果顯示,大部分基因偏離了ENC預(yù)期標(biāo)準(zhǔn)曲線,少數(shù)分布在預(yù)期曲線附近;有7個(gè)密碼子分布在GC3含量-0.05~0.15區(qū)間,其余區(qū)間共有44個(gè)密碼子,而這44個(gè)基因離ENC標(biāo)準(zhǔn)曲線較遠(yuǎn)。PR2-plot分析結(jié)果顯示,落在中線上的基因較少,平面圖右下方的基因分布較多,表明在各堿基使用頻率方面存在偏性,即T>A,G>C。各堿基使用頻率存在偏性(T>A,G>C),說(shuō)明基因組密碼子的偏性更多受選擇效應(yīng)的影響。最終篩選出7個(gè)為最優(yōu)密碼子UUU、CUU、CCU、AAA、GAU、CGA和GGU?!窘Y(jié)論】楸樹葉綠體基因組密碼子偏性較弱,偏好使用A或U結(jié)尾的密碼子,主要受選擇壓力的影響,受突變等其他因素影響較小。

    關(guān)鍵詞: 楸樹;葉綠體基因組;密碼子偏好性;最優(yōu)密碼子

    中圖分類號(hào):S792.99? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼: A 文章編號(hào):2095-1191(2021)10-2735-09

    Abstract:【Objective】Analyzed the codon usage bias of Catalpa bungei C. A. Mey. chloroplast genome to provide theoretical basis for the research of C. bungei and even Catalpa chloroplast genomics,species evolution characteristics and genetic diversity analysis. 【Method】Codon correlation parameters,neutral plots,ENC-plot and PR2-plot were analyzed for 51 chloroplast gene coding regions (CDS) sequences of C. bungei by Codon W 1.4.2 and cusp and chips of EMBOSS expiorer. And through the constitutive high and low expression-biased gene library, ENC value and RSCU value were used as reference criteria to screen the optimal codon. 【Result】The average GCall content of 51 gene codons was 38.48%,and GC content at different positions of codons was as follows: GC1(47.03%)>GC2(39.93%)>GC3(28.09%),indicating that GC at the third position was not evenly distributed. ENC ranged from 34.93 to 56.50,with an average of 46.87. GC2 was extremely significantly positively correlated with GC1 (P<0.01),while GC3 was not significantly positively correlated with GC1 and GC2 (P>0.05).? ENC was significantly positively correlated with GCall,GC1,GC3 and N (P<0.05,the same below).? The codon number (N) was positively correlated with GC3 and ENC. There were 30 codons with RSCU>1,among which 29 codons ended in A or U. Neutral plot analysis results showed that the correlation coefficient between GC3 and GC12 was 0.0188,and there was no significant correlation between them,and the regression coefficient was 0.143. ENC-plot analysis showed that most of the genes deviated from the ENC expected standard curve,and a few were distri-buted near the expected curve. There were 7 codons distributed in the range of GC3 content from -0.05 to 0.15,and there were 44 codons in the other regions,which were far from the ENC standard curve. PR2-plot analysis results showed that fewer genes fell on the middle line,and more genes were distributed in the lower right part of the graph,indicating that there was a bias in the use frequency of each base,namely T>A,G>C.There was bias (T>A,G>C) in the use frequency of each base,indicating that the bias of the genome codon was more influenced by the selection effect. Finally,7 optimal codons were selected as UUU,CUU,CCU,AAA,GAU,CGA and GGU. 【Conclusion】The chloroplast genomes of C. bungei show weak codon bias,preferring to use codons ending in A or U,which is mainly affected by selection pressure and less affected by mutation and other factors.

    Key words:Catalpa bungei C. A. Mey.; chloroplast genome; codon bias; optimal codons

    Foundation item: Fundamental Research Funds for the Central Public Welfare Research Institutes(CAFYBB2017 ZA001-8)

    0 引言

    【研究意義】楸樹(Catalpa bungei C. A. Mey.)是紫葳科(Bignoniaceae)梓屬(Catalpa Scop)落葉喬木,原產(chǎn)我國(guó),分布范圍較廣,主產(chǎn)地為黃河流域和長(zhǎng)江流域,素有“木王”之美稱,其材質(zhì)不翹不裂,用途廣泛,樹體高大挺拔,樹姿優(yōu)美,集用材、觀賞、環(huán)保和醫(yī)用價(jià)值等優(yōu)點(diǎn)于一身,是我國(guó)特有的優(yōu)質(zhì)、藥用和園林觀賞樹種(吳麗華等,2010)。葉綠體基因組不同于細(xì)胞核中的基因組DNA,通常是雙鏈環(huán)狀結(jié)構(gòu),不與蛋白質(zhì)結(jié)合,擁有一套完整的復(fù)制、轉(zhuǎn)錄和翻譯系統(tǒng)(辛雅萱等,2020),其結(jié)構(gòu)保守,并包含了長(zhǎng)、短單拷貝區(qū)(LSC和SSC)和2個(gè)反向重復(fù)區(qū)(IR)共4個(gè)大片段,這些特性有利于分析葉綠體基因組。目前葉綠體基因組已廣泛用于分子進(jìn)化、生物多樣性分析、DNA條形碼、系統(tǒng)發(fā)育和遺傳表達(dá)等研究領(lǐng)域(倪梁紅等,2015;范李強(qiáng)等,2018)。研究表明,同義密碼子在無(wú)任何自然選擇和基因突變等壓力影響的情況下以同樣的概率出現(xiàn),但如果受到影響則會(huì)導(dǎo)致同義密碼子的偏性使用,被稱為密碼子使用偏性(Codon usage bias,CUB),該現(xiàn)象廣泛存在于生物中(吳憲明等,2007;吳妙麗等,2019;王艷等,2021)。研究葉綠體密碼子偏性可為葉綠體基因工程相關(guān)研究打下基礎(chǔ)(季凱凱等,2020)。因此,分析楸樹葉綠體基因組密碼子使用偏性,確定密碼子的使用模式及其形成因素,對(duì)楸樹乃至梓屬植物葉綠體基因組學(xué)、物種進(jìn)化特征及遺傳多樣性分析等研究具有重要意義。【前人研究進(jìn)展】大量研究證實(shí),不同物種的葉綠體基因組密碼子偏性因物種而不同(原曉龍等,2021),但相同物種或者親緣關(guān)系較近的物種之間的密碼子使用偏性相似(王英等,2020)。葉綠體基因組密碼子使用偏性影響的因素除了受選擇壓力和自然壓力外,tRNA豐度(Moriyama and Powell,1997)、mRNA二級(jí)結(jié)構(gòu)(Gu et al.,2004)、外源基因堿基突變(Hershberg and Petrov,2008)及在寄主中的表達(dá)水平(武立偉等,2020)、蛋白表達(dá)效率和準(zhǔn)確性(Lee et al.,2009;Brandls and Hughes,2016)等均會(huì)影響密碼子使用偏性的產(chǎn)生。目前分析密碼子的使用偏性、篩選最優(yōu)密碼子及通過(guò)優(yōu)化密碼子來(lái)提高外源基因在宿主中的表達(dá)水平方面的研究較多(Stephen et al.,2013;Deng et al.,2020),如余桂容等(2015)通過(guò)密碼子優(yōu)化提高蘇云金芽胞桿菌(Bt)cry1Ah基因在轉(zhuǎn)基因玉米中的表達(dá)水平;王超等(2020)對(duì)朝倉(cāng)花椒(Zanthoxylum piperi-tum DC.var. inerme Makino)基因的密碼子使用偏性進(jìn)行分析,確定花椒基因使用偏性較高的密碼子,通過(guò)優(yōu)化密碼子的方式提高花椒基因在物種內(nèi)的表達(dá)水平;柳燕杰等(2020)通過(guò)分析美國(guó)紅梣(Fraxinus pennsylvanica Marsh)葉綠體基因組密碼子使用模式及影響密碼子使用偏性的因素,確定其密碼子偏性和最優(yōu)密碼子,從而預(yù)測(cè)基因的表達(dá)水平,確定未知基因的位置和改良外源基因。隨著葉綠體的基因工程在優(yōu)良品種繁育和種質(zhì)資源保護(hù)中的作用日益凸顯,且大量物種的葉綠體基因組序列相繼公布,對(duì)不同物種的葉綠體基因組的研究日益增多(梁璇等,2021)。近年來(lái),對(duì)葉綠體基因組密碼子偏性已在酸棗(Ziziphus jujuba)(胡曉艷等,2019)、巴山木竹(Arundinaria fargesii)(吳妙麗等,2019)、印度簕竹(Bambusa bambos)(吳妙麗等,2019)、杜梨(Pyrus betulifolia)(辛雅萱等,2020)等多種植物中進(jìn)行研究?!颈狙芯壳腥朦c(diǎn)】我國(guó)楸樹在優(yōu)良新品系的選擇(麻文俊等,2012)、抗旱脅迫的評(píng)價(jià)(董蕾等,2014)及規(guī)模化繁育(孟路等,2020)等方面的研究已較成熟。雖然楸樹葉綠體基因組的測(cè)序現(xiàn)已完成,但未見(jiàn)有關(guān)楸樹密碼子偏性的研究報(bào)道。【擬解決的關(guān)鍵問(wèn)題】利用Codon W 1.4.2及EMBOSS expiorer的cusp和chips軟件對(duì)楸樹葉綠體基因組密碼子的使用偏性進(jìn)行相關(guān)分析,從而確定最優(yōu)密碼子及密碼子偏性的主要影響因素,旨在為楸樹的遺傳變異及基因組學(xué)的研究提供強(qiáng)有力的分子理論依據(jù)。

    1 材料和方法

    1. 1 序列的獲取和篩選

    從NCBI數(shù)據(jù)庫(kù)(http//www.ncbi.nlm.gov/)中下載楸樹葉綠體基因組序列(GenBank登錄號(hào)為MT610898,序列全長(zhǎng)158238 bp),經(jīng)篩選共獲得84個(gè)基因的編碼區(qū)(CDS)序列。為了降低因序列較短而造成的誤差,從84個(gè)基因序列中去除重復(fù)序列、長(zhǎng)度小于300 bp的序列及內(nèi)部含有終止密碼子和不完整的序列,最終獲得起始密碼子為ATG,終止密碼子為TAA、TGA和TAG的51個(gè)基因序列,用于后續(xù)序列分析。

    1. 2 密碼子相關(guān)參數(shù)計(jì)算

    通過(guò)Codon W 1.4.2分析51個(gè)基因序列,并獲取其相對(duì)同義密碼子使用度(RSCU)(即某一特定密碼子所編碼氨基酸時(shí)的所有密碼子相對(duì)使用次數(shù)),把RSCU作為判斷密碼子是否存在偏性的標(biāo)準(zhǔn):當(dāng)RSCU=1時(shí),證明該密碼子使用上無(wú)偏性;RSCU>1時(shí),則表明其使用頻率相對(duì)較高,存在偏性(牛元等,2018;胡振民等,2020)。利用EMBOSS expiorer中的cusp軟件分析不同基因序列的GC1、GC2、GC3(分別代表密碼子第1位、第2位和第3位的GC含量)和GCall(密碼子總GC含量)含量;并用EMBOSS expiorer中的chips軟件分析有效密碼子數(shù)(ENC)(即密碼子非均衡使用的偏好程度),ENC的取值范圍為20~61(馮瑞云等,2019)。最后,使用SPSS 21.0作相關(guān)分析。

    1. 3 中性繪圖分析

    根據(jù)不同基因密碼子的GC3含量和GC12含量[即GC12=(GC1+GC2)/2]進(jìn)行中性繪圖,從而分析判斷影響密碼子偏性形成的因素(王鵬良等,2018)。若GC3與GC12相關(guān)性強(qiáng),則說(shuō)明突變壓力對(duì)密碼子3個(gè)位置的堿基的使用影響較大;若GC1、GC2與GC3相關(guān)性弱,則表明選擇壓力對(duì)密碼子使用影響較大(秦政等,2018;惠小涵等,2020)。

    1. 4 ENC-plot繪圖分析

    根據(jù)不同基因密碼子的GC含量和ENC,利用ENC-plot制作散點(diǎn)圖,并計(jì)算ENC的標(biāo)準(zhǔn)曲線,公式為:ENC=[29GC32+(1?GC3)2]+GC3+2,從而分析突變壓力在密碼子使用過(guò)程中的作用(尹為治等,2020)。如果各基因分布在ENC預(yù)期曲線的周圍,說(shuō)明葉綠體基因組密碼子的使用偏性受突變壓力的影響較大,而如果各基因分布在ENC預(yù)期曲線的下方或偏離預(yù)期曲線,則表明葉綠體基因組密碼子的偏性受到選擇壓力的影響較大(王鵬良等,2018)。

    1. 5 PR2-plot繪圖分析

    為了探究不同基因密碼子第3位堿基的嘌呤與嘧啶之間的突變壓力是否平衡,即A、T與G、C之間均衡分布的情況,對(duì)由同義密碼子編碼的氨基酸第3位上A、T、C、G 4種堿基組成情況進(jìn)行PR2-plot繪圖分析(金剛,2019),以確定楸樹葉綠體基因組密碼子偏性的影響因素。其中,以G3/(G3+C3)、A3/(A3+T3)分別為橫、縱坐標(biāo),PR2偏倚圖的中心(0.5,0.5)表示A=T且C=G(宋輝等,2015)。其余點(diǎn)由中心點(diǎn)向該點(diǎn)發(fā)出的矢量表明基因的偏倚程度及方向(陸奇豐等,2020)。

    1. 6 最優(yōu)密碼子分析

    為了選出楸樹葉綠體基因組中的最優(yōu)密碼子,利用Codon W 1.4.2計(jì)算每個(gè)密碼子的RSCU,把RSCU>1的密碼子定義為高頻密碼子。以5個(gè)基因的ENC作為偏性標(biāo)準(zhǔn),對(duì)其進(jìn)行從小到大排序,取前后兩端10%基因(前后各5個(gè)基因)構(gòu)建高、低表達(dá)偏性基因庫(kù)。兩組間的差值ΔRSCU≥0.08(ΔRSCU=RSCU高表達(dá)- RSCU低表達(dá))的密碼子作為高表達(dá)優(yōu)越密碼子,最終將同時(shí)滿足RSCU>1和ΔRSCU≥0.08的密碼子作為楸樹葉綠體基因組的最優(yōu)密碼子(胡莎莎等,2016)。

    2 結(jié)果與分析

    2. 1 密碼子堿基組成及偏好性分析

    由表1可知,51個(gè)基因密碼子的平均GCall含量為38.48%,平均GC1、GC2和GC3含量分別為47.03%、39.93%和28.09%,即GC3含量明顯低于GC1和GC2含量,說(shuō)明第3位的GC并非均勻分布,分布趨勢(shì)為GC1>GC2>GC3??梢?jiàn),密碼子第3位以A或U(T)結(jié)尾的頻率高于以G或C結(jié)尾的頻率。51個(gè)基因的ENC為34.93~56.50,平均為46.87,說(shuō)明楸樹葉綠體基因組的密碼子偏性較弱。

    對(duì)楸樹葉綠體基因的密碼子數(shù)量(N)、ENC和各密碼子的GC含量進(jìn)行相關(guān)分析,結(jié)果如表2所示。GC2與GC1呈極顯著正相關(guān)(P<0.01,下同),而GC3與GC1和GC2無(wú)顯著正相關(guān)(P>0.05,下同),說(shuō)明密碼子的第1位與第2位的堿基相似度較高,但二者與第3位的堿基相似度較低;GCall與GC1和GC3均呈極顯著正相關(guān),與GC2呈顯著正相關(guān)(P<0.05,下同);ENC除了與GC2無(wú)顯著負(fù)相關(guān)外,與GCall、GC1、GC3和N均呈顯著正相關(guān),說(shuō)明密碼子3個(gè)位置的堿基組成對(duì)ENC的影響有所不同,其中,第1位和第3位堿基組成對(duì)密碼子的偏性有影響。密碼子數(shù)(N)除了與GC3和ENC呈顯著正相關(guān)外,與其他數(shù)值均無(wú)顯著相關(guān)。

    由表3可知,RSCU>1的密碼子共有30個(gè),其中,以U(16個(gè))或A(13個(gè))結(jié)尾的密碼子數(shù)量明顯多于以G(1個(gè))結(jié)尾的密碼子,說(shuō)明以A或U結(jié)尾的密碼子是楸樹葉綠體基因組偏好的密碼子。

    2. 2 中性繪圖分析結(jié)果

    由圖1可知,GC12含量為0.20~0.37;GC3含量為0.35~0.55,二者的相關(guān)系數(shù)和回歸系數(shù)分別為0.0188和0.143,表明在密碼子的第1、2位堿基與第3位堿基組成相關(guān)性較弱,楸樹葉綠體基因組編碼基因保守性較高,密碼子偏好性受選擇壓力的影響較大。

    2. 3 ENC-plot繪圖分析結(jié)果

    由圖2可知,在51個(gè)基因中,大部分基因偏離了ENC預(yù)期標(biāo)準(zhǔn)曲線,少數(shù)分布在預(yù)期曲線附近,說(shuō)明楸樹葉綠體基因組的偏性受選擇壓力的影響較大,受基因突變壓力的影響較小;位于標(biāo)準(zhǔn)曲線上方的基因密碼子偏性較弱,而位于標(biāo)準(zhǔn)曲線下方的基因密碼子偏性較強(qiáng)。由不同GC3含量區(qū)間分布的密碼子數(shù)量(表4)可知,有7個(gè)密碼子分布在-0.05~0.15區(qū)間,其余區(qū)間共有44個(gè)密碼子,而這44個(gè)基因距離ENC標(biāo)準(zhǔn)曲線較遠(yuǎn)(圖2),表明楸樹葉綠體基因組密碼子的使用偏性主要受選擇壓力的影響較大,而受其他因素(如堿基突變等)的影響相對(duì)較小。

    2. 4 PR2-plot繪圖分析結(jié)果

    由圖3可知,51個(gè)基因在PR2平面圖的4個(gè)部分并非均勻分布,且落在中線上的基因較少,平面圖右下方的基因分布較多,表明各堿基使用頻率存在偏性,即T>A,G>C,說(shuō)明選擇壓力不是影響楸樹葉綠體基因組密碼子偏性的唯一因素,突變等其他因素也對(duì)楸樹葉綠體基因組的密碼子偏性產(chǎn)生一定的影響。

    2. 5 最優(yōu)密碼子的篩選

    從51個(gè)楸樹葉綠體基因選取中5個(gè)高表達(dá)基因(petA、ndhC、rps2、atpE和rps4)和5個(gè)低表達(dá)基因(cemA、atpF、psbA、rps18和rps14),然后利用Codon W 1.4.2計(jì)算出兩組基因同義密碼子的RSCU,再計(jì)算出ΔRSCU,結(jié)果(表5)表明,高頻密碼子有24個(gè)(RSCU≥1),高表達(dá)優(yōu)越密碼子有12個(gè)(ΔRSCU≥0.08,用*表示),ΔRSCU≥0.30(用**表示)的密碼子有10個(gè),ΔRSCU≥0.50(用***表示)密碼子有4個(gè)。同時(shí)滿足高頻(ΔRSCU≥0.08)和高表達(dá)(RSCU≥1)的密碼子有7個(gè),確定為最優(yōu)密碼子,分別是UUU、CUU、CCU、AAA、GAU、CGA和GGU,其中5個(gè)以U結(jié)尾,2個(gè)以A結(jié)尾。

    3 討論

    GC3是密碼子偏性的一個(gè)重要判斷指標(biāo)。本研究對(duì)楸樹葉綠體基因密碼子3個(gè)位置的GC含量進(jìn)行分析,結(jié)果表明,GC1和GC2含量顯著高于GC3含量,說(shuō)明楸樹葉綠體基因組偏好使用U或A結(jié)尾的密碼子,RSCU分析結(jié)果也證實(shí)該結(jié)論。在糜子(Panicum miliaceum L.)(劉慧等,2017)、普通油茶(Camellia oleifera Abel)(王鵬良等,2018)和西南樺(Betula alnoides Hamilton)(原曉龍等,2020b)等植物的葉綠體基因組偏性的研究也發(fā)現(xiàn),這些植物偏好使用A或U結(jié)尾的密碼子。本研究最終篩選出的7個(gè)最優(yōu)密碼子,均以U或A結(jié)尾,與杜梨(Pyrus betulifolia Bge.)(辛雅萱等,2020)的最優(yōu)密碼子研究結(jié)果一致。

    選擇壓力是引起密碼子使用偏性的主要因素,突變等其他方面的因素也發(fā)揮一定作用(趙森等,2020)。本研究發(fā)現(xiàn),楸樹葉綠體基因的ENC為34.93~56.50,平均為46.87,說(shuō)明楸樹葉綠體基因密碼子偏性較弱,與大花香水月季(Rosa odorata var. gigantean Rehd. & Wils)(牛元等,2018)、酸棗[Zizi-phus jujube var. spinosa(Bunge) Hu ex H. F. Chow](胡曉艷等,2019)、灰毛漿果楝[Cipadessa cinerascens(Pellegr.) H. -M.](原曉龍等,2020c)等植物的葉綠體密碼子使用偏性研究結(jié)果一致,由此推測(cè)葉綠體基因組結(jié)構(gòu)的保守性是導(dǎo)致密碼子使用偏性的主要因素之一。結(jié)合相關(guān)分析和中性繪圖分析結(jié)果發(fā)現(xiàn),堿基突變也是影響楸樹葉綠體基因組偏性的主要因素之一,與羅火林等(2016)研究結(jié)論一致。本研究的PR2-plot和ENC-plot分析結(jié)果均表明,楸樹葉綠體基因組密碼子的使用偏性主要受選擇壓力的影響較大,而受其他因素(如堿基突變等)的影響相對(duì)較小,與云南藍(lán)果樹(Nyssa yunnanensis W. C. Yin)(原曉龍等,2020a)、灰毛漿果楝(C. cinerascens (Pellegr.) H. -M.)(原曉龍等,2020b)和西南樺(Betula alnoides Hamilton)(原曉龍等,2020c)等研究結(jié)果一致,但與水稻(Oryza sativa Linn.)(劉慶坡和薛慶中,2004)研究結(jié)果不一致。結(jié)論存在差異可能因?yàn)榛虮磉_(dá)程度和密碼子堿基組成不同所引起。由此可知,物種的葉綠體基因組密碼子使用偏性的影響因素眾多,且不同物種的影響因素各不相同(傅建敏等,2017)。今后應(yīng)基于本研究結(jié)果對(duì)楸樹葉綠體功能基因的進(jìn)化產(chǎn)生是否受選擇壓力影響進(jìn)行深入研究。

    4 結(jié)論

    楸樹葉綠體基因組密碼子偏性較弱,偏好使用A或U結(jié)尾的密碼子,主要受選擇壓力的影響,受突變等其他因素影響較弱。

    參考文獻(xiàn):

    董蕾,陳博,李吉躍,王軍輝,謝堃. 2014. 5個(gè)楸樹無(wú)性系抗旱性的綜合評(píng)價(jià)[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào),34(4):30-35. [Dong L,Chen B,Li J Y,Wang J H,Xie K. 2014. Drought resistance evaluation of 5 Catalpa bungei clones[J]. Journal of Central South University of Forestry & Technology,34(4):30-35.] doi:10.3969/j.issn.1673-923X. 2014.04.007.

    范李強(qiáng),胡歡,鄭洪蕾,王田靜,王亞玲,馬濤,毛康珊. 2018. 響葉楊(楊屬)葉綠體基因組測(cè)序與比較分析[J]. 四川大學(xué)學(xué)報(bào)(自然科學(xué)版),55(1):165-171. [Fan L Q,Hu H,Zheng H L,Wang T J,Wang Y L,Ma T,Mao K S. 2018. Complete sequence and comparative analysis of the chloroplast genome of the Chinese aspen(Populus adenopoda,Salicaceae)[J]. Journal of Sichuan University(Naiural Seienee Edtion),55(1):165-171.] doi:10.3969/j.issn.0490-6756.2018.01.027.

    馮瑞云,梅超,王慧杰,雷夢(mèng)林,田翔,侯雅靜. 2019. 籽粒莧葉綠體基因組密碼子偏好性分析[J]. 中國(guó)草地學(xué)報(bào),41(4):8-15. [Feng R Y,Mei C,Wang H J,Lei M L,Tian X,Hou Y J. 2019. Analysis of codon usage in the chloroplast genome of grain amaranth(Amaranthus hypochondriacus L.)[J]. Chinese Journal of Grassland,41(4):8-15.] doi:10.16742/j.zgcdxb.20180227.

    傅建敏,索玉靜,劉慧敏,譚曉風(fēng). 2017. 柿屬植物葉綠體蛋白質(zhì)編碼基因密碼子用法[J]. 經(jīng)濟(jì)林研究,35(2):38-44. [Fu J M,Suo Y J,Liu H M,Tan X F. 2017. Analysis on codon usage in the chloroplast protein-coding genes of Diospyros spp.[J]. Nonwood Forest Research,35(2):38-44.] doi:10.14067/j.cnki. 1003-8981.2017.02.006.

    胡莎莎,羅洪,吳琦,姚慧鵬. 2016. 苦蕎葉綠體基因組密碼子偏愛(ài)性分析[J]. 分子植物育種,14(2):309-317. [Hu S S,Luo H,Wu Q,Yao H P. 2016. Analysis of codon bias of chloroplast genome of tartary buckwheat[J]. Molecular Plant Breeding,14(2):309-317.] doi:10.13271/j.mpb. 014.000309.

    胡曉艷,許艷秋,韓有志,杜淑輝. 2019. 酸棗葉綠體基因組密碼子使用偏性性分析[J]. 森林與環(huán)境學(xué)報(bào),19(6):621-628. [Hu X Y,Xu Y Q,Han Y Z,Du S H. 2019. Codon usage bias analysis of the chloroplast genome of Ziziphus jujuba var. spinosa[J]. Journal of Forest and Environment,19(6):621-628.] doi: 10.13324/j.cnki.jfcf.2019.06.009.

    胡振民,萬(wàn)青,李歡,李榮林,王楓,楊亦揚(yáng). 2019. 茶樹CsNRT1.1基因密碼子使用特性分析[J]. 江蘇農(nóng)業(yè)學(xué)報(bào),35(4):896-903. [Hu Z M,Wan Q,Li H,Li R L,Wang F,Yang Y Y. 2019. Analysis of codon usage features of CsNRT1.1 gene in Camellia sinensis[J]. Jiangsu Journal of Agricultural Sciences,35(4):896-903.] doi:10.3969/j.issn. 1000-4440.2019.04.021.

    惠小涵,程婷婷,柯衛(wèi)東,郭宏波. 2020. 蓮藕PPO基因密碼子偏好性特征分析[J]. 江蘇農(nóng)業(yè)學(xué)報(bào),36(2):438-446. [Hui X H,Cheng T T,Ke W D,Guo H B. 2020. Analysis on codon preference of PPO gene in lotus root[J]. Jiangsu Journal of Agricultural Sciences,36(2):438-446.] doi:10.3969 /j.issn.1000-4440.2020.02.026.

    季凱凱,宋希強(qiáng),陳春國(guó),李革,謝尚淺. 2020. 木蘭科葉綠體基因組的密碼子使用特征分析[J]. 中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào),22(11):52-62. [Ji K K,Song X Q,Chen C G,Li G,Xie S Q. 2020. Codon usage profiling of chloroplast genome in Magnoliaceae[J]. Journal of Agricultural Science and Technology,22(11):52-62.] doi:10.13304/j.nykjdb.2019. 0937.

    金剛,王麗萍,龍凌云,吳鳳,唐玉娟,覃劍鋒,危丹妮,黃秋偉,蘇文潘. 2019. 普通野生稻線粒體蛋白質(zhì)編碼基因密碼子使用偏好性的分析[J]. 植物科學(xué)學(xué)報(bào),37(2):188-197. [Jin G,Wang L P,Long L Y,Wu F,Tang Y J,Qin J F,Wei D N,Huang Q W,Su W P. 2019. Analysis of codon usage bias in the mitochondrial protein-coding genes of Oryza rufipogon[J]. Plant Science Journal,37(2):188-197.] doi:10.11913/PSJ.2095-0837.2019.20188.

    梁璇,李鵬,許冬梅,賈小云,王文斌. 2021. 紫蘇葉綠體全基因組特征分析[J]. 山西農(nóng)業(yè)科學(xué),49(3):265-272. [Liang X,Li P,Xu D M,Jia X Y,Wang W B. 2021. Characteristics analysis of whole chloroplast genome in Perilla frustescens[J]. Journal of Shanxi Agricultural Sciences,49(3):265-272.] doi:10.3969/j.issn.1002-2481. 2021.03.02.

    劉慧,王夢(mèng)醒,岳文杰,邢光偉,葛靈巧,聶小軍,宋衛(wèi)寧. 2017. 糜子葉綠體基因組密碼子使用偏性的分析[J]. 植物科學(xué)學(xué)報(bào),35(3):362-367. [Liu H,Wang M X,Yue W J,Xing G W,Ge L Q,Nie X J,Song W N. 2017. Analysis of codon usage in the chloroplast genome of broomcorn millet(Panicum miliaceum L.)[J]. Plant Scien-ce Journal,35(3):362-367.] doi:10.11913/PSJ.2095-0837.2017.30362.

    劉慶坡,薛慶中. 2004. 粳稻葉綠體基因組的密碼子用法[J]. 作物學(xué)報(bào),30(12):1220-1224. [Liu Q P,Xue Q Z. 2004. Codon usage in the chloroplast genome of rice(Oryza sativa L ssp japonica)[J]. Acta Agronomica Sinica,30(12):1220-1224.] doi:10.3321/j.issn:0496-3490. 2004.12.008.

    柳燕杰,田旭平,李倩. 2020. 美國(guó)紅梣葉綠體基因組密碼子偏好性分析[J]. 江蘇農(nóng)業(yè)科學(xué),48(15):83-88. [Liu Y J,Tian X P,Li Q. 2020. Codon usage bias analysis of the chloroplast genome of Fraxinus pennsylvanica[J]. Jiangsu Agricultural Sciences,48(15):83-88.] doi:10. 15889/ j. issn. 1002 - 1302. 2020. 15.014.

    陸奇豐,駱文華,黃至歡. 2020. 兩種梧桐葉綠體基因組密碼子使用偏性分析[J]. 廣西植物,40(2):173-183. [Lu Q F,Luo W H,Huang Z H. 2020. Codon usage bias of chloroplast genome from two species of Firmiana Marsili[J]. Guihaia,40(2):173-183.] doi:10.11931/guihaia.gxzw 201811035.

    羅火林,張亞楠,羅麗萍,熊冬金,楊柏云. 2016. 蘭科植物FT基因的密碼子偏性分析[J]. 分子植物育種,14(1):51-58. [Luo H L,Zhang Y N,Luo L P,Xiong D J,Yang B Y. 2016. The codon usage bias of FLOWERING LOCUS T gene(FT) in Orchidaceae[J]. Molecular Plant Bree-ding,14(1):51-58.] doi:10.13271/j.mpb.014.000051.

    麻文俊,張守攻,王軍輝,董菊蘭. 2012. 楸樹無(wú)性系苗期N素利用差異和高產(chǎn)無(wú)性系選擇[J]. 林業(yè)科學(xué),48(10):157-162. [Ma W J,Zhang S G,Wang J H,Dong J L. 2012. Variation of nitrogen utilization among Catalpa bungei clones at nursery stage and high-yield clones selection[J]. Scientia Silvae Sinicae,48(10):157-162.] doi:10.11707/ j.1001-7488.20121025.

    孟路,劉勇,王小平,李進(jìn)宇,賀國(guó)鑫,薛敦孟,邢麗霞,李世安. 2020. 4個(gè)楸樹品種組培快繁能力的比較研究[J]. 中南林業(yè)科技大學(xué)學(xué)報(bào),40(1):82-88. [Meng L,Liu Y,Wang X P,Li J Y,He G X,Xue D M,Xing L X,Li S A. 2020. Comparison study on rapid reproductive capacity of four cultivars of Catalpa bungei in tissue culture[J]. Journal of Central South University of Forestry & Technology,40(1):82-88.] doi:10.14067/j.cnki.1673-923x.2020. 01.010.

    倪梁紅,趙志禮,米瑪. 2015. 藥用植物葉綠體基因組研究進(jìn)展[J].中藥材,38(9):1990-1994. [Ni L H,Zhao Z L,Mi M.2015. Research progress on chloroplast genomes of medicinal plants[J]. Journal of Chinese Medicinal Materials,38(9):1990-1994.] doi:10.13863/j.issn1001-4454. 2015.09.046.

    牛元,徐瓊,王崳德,戴立蘭,莊健,趙亞蘭. 2018. 大花香水月季葉綠體基因組密碼子使用偏性分析[J]. 西北林學(xué)院學(xué)報(bào),33(3):123-130. [Niu Y,Xu Q,Wang Y D,Dai L L,Zhuang J,Zhao Y L. 2018. An analysis on codon usage bias of chloroplast genome of Rosa odorata var. gigantean[J]. Journal of Northwest Forestry University,33(3):123- 130.] doi:10.3969/j.issn.1001-7461.2018.03.19.

    秦政,鄭永杰,桂麗靜,謝谷艾,伍艷芳. 2018. 樟樹葉綠體基因組密碼子偏好性分析[J]. 廣西植物,38(10):1346-1355. [Qin Z,Zheng Y J,Gui L J,Xie G A,Wu Y F. 2018. Codon usage bias analysis of chloroplast genome of camphora tree(Cinnamomum camphora)[J]. Guihaia,38(10):1346-1355.] doi:10.11931/guihaia.gxzw201805 023.

    宋輝,王鵬飛,馬登超,夏晗,趙傳志,張燁,趙術(shù)珍. 2015. 蒺藜苜蓿WRKY轉(zhuǎn)錄因子密碼子使用偏好性分析[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào),23(2):203-212.[Song H,Wang P F,Ma D C,Xia H,Zhao C Z,Zhang Y,Zhao S Z. 2015. Analysis of codon usage bias of WRKY transcription factors in Medicago truncatula[J]. Journal of Agricultural Biotechnology,23(2):203-212.] doi:10.3969/j.issn.1674-7968.2015.02.007.

    王超,趙懿琛,趙德剛. 2020. 朝倉(cāng)花椒轉(zhuǎn)錄組密碼子偏好性分析[J]. 分子植物育種,18(21):6989-6996. [Wang C,Zhao Y C,Zhao D G. 2020. Preference analysis of transcriptome codon in Zanthoxylum piperitum DC. var. inerme Makino[J]. Molecular Plant Breeding,18(21):6989-6996.] doi:10.1016/S0309.1740(99)00142-4.

    王鵬良,楊利平,吳紅英,農(nóng)有良,吳雙成,肖玉菲,覃子海,王華宇,劉海龍. 2018. 普通油茶葉綠體基因組密碼子偏好性分析[J]. 廣西植物,38(2):135-144. [Wang P L,Yang L P,Wu H Y,Nong Y L,Wu S C,Xiao Y F,Qin Z H,Wang H Y,Liu H L. 2018. Condon preference of chloroplast genome in Camellia oleifera[J]. Guihaia,38(2):135-144.] doi:10.11931/guihaia.gxzw201708001.

    王艷,趙懿琛,趙德剛. 2021. 杜仲基因密碼子使用模式分析[J]. 廣西植物,41(2):274-282. [Wang Y,Zhao Y C,Zhao D G. 2021. Analysis of codon usage pattern of Eucommia ulmoides[J]. Guihaia,41(2):274-282.] doi:10. 11931/guihaia.gxzw201910008.

    王英,張?zhí)瑥埫蟼?,趙玉潔,袁偉程,劉翠玉,苑兆和. 2020. 石榴等陸生植物Actin1基因密碼子偏向性與進(jìn)化分析[J]. 分子植物育種,18(6):1799-1807. [Wang Y,Zhang T Q,Zhang M W,Zhao Y J,Yuan W C,Liu C Y,Yuan Z H. 2020. Codon usage bias and evolution analyses of Actin1 genes from Punica granatum and other land plants[J]. Molecular Plant Breecling,18(6):1799-1807.] doi:10. 13271/j.mpb.018.001799.

    吳麗華,王軍輝,林娟. 2010. 楸樹植物資源的研究概況[J]. 上海交通大學(xué)學(xué)報(bào)(農(nóng)業(yè)科學(xué)版),28(1):91-96. [Wu L H,Wang J H,Lin J. 2010. A survey of the studies on the resources of Catalpa bungei[J]. Journal of Shanghai Jiaotong University(Agricultural Science),28(1):91-96.] doi:10.3969/j.issn.1671-9964.2010.01.0017.

    吳妙麗,陳世品,陳輝. 2019. 竹亞科葉綠體基因組的密碼子使用偏性分析[J]. 森林與環(huán)境學(xué)報(bào),39(1):9-14. [Wu M L,Chen S P,Chen H. 2019. Condon preference of chloroplast genome of Bambusoideae[J]. Journal of Forest and Environment,39(1):9-14.] doi:10.13324/j.cnki.jfcf.2019.01.002.

    吳憲明,吳松鋒,任大明,朱云平,賀福初. 2007. 密碼子偏性的分析方法及相關(guān)研究進(jìn)展[J]. 遺傳,(4):420-426. [Wu X M,Wu S F,Ren D M,Zhu Y P,He F C. 2007. The analysis method and progress in the study of codon bias[J]. Hereditas(Beijing),(4):420-426.] doi:10.3321/j.issn:0253-9772.2007.04.006.

    武立偉,崔英賢,聶麗萍,徐志超,王瑀,宋經(jīng)元,焦連魁,姚輝. 2020. 細(xì)莖石斛葉綠體全基因組序列特征及系統(tǒng)發(fā)育分析[J]. 藥學(xué)學(xué)報(bào),55(5):302-312. [Wu L W,Cui Y X,Nie L P,Xu Z C,Wang Y,Song J Y,Jiao L K,Yao H. 2020. The characteristics of complete chloroplast genome sequence and phylogenetic analysis of Dendrobium moniliforme[J]. Acta Pharmaceutica Sinica,55(5):302-312.] doi:10.16438/j.0513-4870.2019-0940.

    辛雅萱,董章宏,瞿紹宏,劉成,葉鵬,辛培堯. 2020. 杜梨葉綠體基因組密碼子偏性分析[J]. 河北農(nóng)業(yè)大學(xué)學(xué)報(bào),43(6):51-59. [Xin Y X,Dong Z H,Qu S H,Liu C,Ye P,Xin P Y. 2020. Analysis on codon usage bias of chloroplast genome in Pyrus betulifolia Bge[J]. Journal of Hebei Agricultural University,43(6):51-59.] doi:10.13320/j.cnki.jauh.2020.0112.

    尹為治,方正,黃良鴻,龍文興,李佳靈. 2020. 海南2種龍腦香科植物葉綠體基因組密碼子偏好性分析[J].林業(yè)調(diào)查規(guī)劃,45(6):19-26. [Yin W Z,F(xiàn)ang Z,Huang L H,Long W X,Li J L. 2020. Codon usage bias in chloroplast genomes of two species of Dipterocarpaceae in Hainan[J]. Forest Inventory and Planning,45(6):19-26.] doi:10. 3969/j.issn.1671-3168.2020.06.005.

    余桂容,曹穎,杜文平,宋軍,陳謙,徐利遠(yuǎn). 2015. 抗蟲耐草甘膦雙價(jià)基因表達(dá)載體的構(gòu)建及玉米轉(zhuǎn)化[J]. 西南農(nóng)業(yè)學(xué)報(bào),28(2):475-479. [Yu G R,Cao Y,Du W P,Song J,Chen Q,Xu L Y. 2015. Construction of binary vector obtaining insect & herbicide resistance genes and transformation of maize[J]. Southwest China Journal of Agricultural Sciences,28(2):475-479.] doi:10.16213/j.cnki.scjas.2015.02.006.

    原曉龍,康洪梅,王毅. 2020a. 云南藍(lán)果樹葉綠體基因組密碼子偏好性分析[J]. 西北林學(xué)院學(xué)報(bào),35(4):26-31. [Yuan X L,Kang H M,Wang Y. 2020a. Codon usage bias analysis of chloroplast genome in Nyssa yunnanensis[J]. Journal of Northwest Forestry University,35(4):26-31.] doi:10.3969/j.issn.1001-7461.2020.04.05.

    原曉龍,李云琴,王毅,張勁峰. 2020b. 西南樺葉綠體基因組密碼子偏好性分析[J]. 基因組學(xué)與應(yīng)用生物學(xué),39(12):5758-5764. [Yuan X L,Li Y Q,Wang Y,Zhang J F. 2020b. Analysis of codon usage in the chloroplast genome of Betula alnoides[J]. Genomics and Applied Bio-logy,39(12):5758-5764.] doi:10.13417/j.gab.039.005758.

    原曉龍,王毅,張勁峰. 2020c. 灰毛漿果楝葉綠體基因組密碼子使用特征分析[J]. 福建林學(xué)院學(xué)報(bào),40(2):195-202. [Yuan X L,Wang Y,Zhang J F. 2020c. Characterization of codon usage in Cipadessa cinerascens chloroplast genome[J]. Journal of Fujian College of Forestry,40(2):195-202.] doi:10.13324/j.cnki.jfcf.2020.02.013.

    原曉龍,劉音,康洪梅,陳中華,李云琴,王毅. 2021. 蒜頭果葉綠體基因組密碼子偏好性分析[J]. 西南林業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)),41(3):15-22. [Yuan X L,Liu Y,Kang H M,Chen Z H,Li Y Q,Wang Y. 2021. Analysis of codon usage bias in chloroplast genome of Malania oleifera[J]. Journal of Southwest Forestry University,41(3):15-22.] doi:10.11929/j.swfu.202001047.

    趙森,鄧力華,陳芬. 2020. 秋茄葉綠體基因組密碼子使用偏好性分析[J]. 森林與環(huán)境學(xué)報(bào),40(5):534-541. [Zhao S,Deng L H,Chen F. 2020. Codon usage bias of chloroplast genome in Kandelia obovata[J]. Journal of Forest and Environment,40(5):534-541.] doi:10.13324/j.cnki.jfcf.2020.05.011.

    Brandls G,Hughes D. 2016. The selective advantage of synonymous codon usage bias in Salmonella[J]. PLoS Genetics,12(3):e1005926. doi:10.1371/journal.pgen.10059 26.

    Deng Y,Hedayioglu F L,Kalfon J,Chu D,Haar T. 2020. Hidden patterns of codon usage bias across kingdoms[J]. Journal of the Royal Society Interface,17(163):20190819.] doi:10.1098/rsif.2019.0819.

    Gu W J,Zhou T,Ma J M,Sun X, Lu Z H. 2004. The relationship between synonymous codon usage and protein structure in Escherichia coli and Homo sapiens[J]. Biosystems,73(2):89-97. doi:10.1016/j.biosystems.2003. 10.001.

    Hershberg R,Petrov D A. 2008. Selection on codon bias[J]. Annual Review of Genetics,42(1):287-299. doi:10.1021/ cm0010349.

    Lee S F,Li Y J,Halperin S A. 2009. Overcoming codon-usage bias in heterologous protein expression in Streptococcus gordonii[J]. Microbiology,155(11):3581-3588. doi:10.1099/mic.0.030064-0.

    Moriyama E N,Powell J R. 1997. Codon usage bias and tRNA abundance in Drosophila[J]. Journal of Molecular Evolution,45(5):514-523. doi:10.1007/PL00006256.

    Stephen J S,Jeremy D V,Robert M,Li F L,John M,Laszlo Z. 2013. Expression of chicken parvovirus VP2 in chicken embryo fibroblasts requires codon optimization for production of naked DNA and vectored Meleagrid herpesvirus type 1 vaccines[J]. Virus Genes,47(2):259-267. doi:10.1007/s11262-013-0944-9.

    (責(zé)任編輯 陳 燕)

    猜你喜歡
    楸樹
    勤耕不輟數(shù)十載 一片丹心育“木王”
    ——記楸樹育種專家趙鯤
    我和你的心仍如峭壁
    海燕(2021年10期)2021-11-12 23:02:04
    幸福的一日
    ——致秋天的花楸樹
    北方人(2021年19期)2021-10-29 00:46:00
    楸樹的應(yīng)用及育苗技術(shù)
    楸樹組織培養(yǎng)研究進(jìn)展
    我們?nèi)ツ吧拇迩f看楸樹
    楸樹
    草原(2018年10期)2018-12-21 12:28:38
    幸福一日
    幸福一日
    視野(2014年21期)2014-10-21 13:23:37
    致秋天的花楸樹
    日本av免费视频播放| 欧美精品一区二区免费开放| 国产成人aa在线观看| 久久综合国产亚洲精品| 亚洲美女视频黄频| 成人亚洲精品一区在线观看| 97在线人人人人妻| 国产精品一区www在线观看| 欧美精品人与动牲交sv欧美| 熟妇人妻不卡中文字幕| 久久久精品94久久精品| 成人手机av| 赤兔流量卡办理| 超碰97精品在线观看| 日本欧美国产在线视频| av福利片在线| 十八禁网站网址无遮挡| 久久精品夜色国产| 欧美性感艳星| 日本午夜av视频| 狠狠婷婷综合久久久久久88av| 亚洲av不卡在线观看| 午夜福利视频精品| 亚洲内射少妇av| av福利片在线| 美女大奶头黄色视频| 人人妻人人澡人人爽人人夜夜| 国产69精品久久久久777片| 国产精品一区二区在线不卡| 成人手机av| av有码第一页| 国产成人午夜福利电影在线观看| 晚上一个人看的免费电影| 91精品伊人久久大香线蕉| 亚洲美女视频黄频| 22中文网久久字幕| 赤兔流量卡办理| 久久久久国产精品人妻一区二区| 国国产精品蜜臀av免费| www.色视频.com| 婷婷色综合大香蕉| .国产精品久久| 久久狼人影院| 久久精品国产亚洲网站| 少妇猛男粗大的猛烈进出视频| 国产欧美另类精品又又久久亚洲欧美| 午夜av观看不卡| 国产欧美亚洲国产| 亚洲国产毛片av蜜桃av| 国产精品嫩草影院av在线观看| 欧美国产精品一级二级三级| av免费在线看不卡| 亚洲欧美成人综合另类久久久| 大香蕉久久网| 亚洲国产毛片av蜜桃av| 日韩 亚洲 欧美在线| 国产综合精华液| 中国三级夫妇交换| 国产一区二区三区av在线| 黄色一级大片看看| 丝袜喷水一区| 午夜精品国产一区二区电影| 一本大道久久a久久精品| 美女xxoo啪啪120秒动态图| 观看av在线不卡| 九草在线视频观看| 男人添女人高潮全过程视频| 国产精品蜜桃在线观看| 久久久精品区二区三区| 免费日韩欧美在线观看| 午夜福利视频精品| 色视频在线一区二区三区| 亚洲精品视频女| 国产精品99久久久久久久久| 国产日韩欧美在线精品| 亚洲国产欧美日韩在线播放| 少妇熟女欧美另类| 精品久久久精品久久久| 精品少妇久久久久久888优播| 黑人猛操日本美女一级片| 观看av在线不卡| 18在线观看网站| 人体艺术视频欧美日本| 99国产精品免费福利视频| 亚洲欧美成人精品一区二区| 三级国产精品片| 丰满迷人的少妇在线观看| 少妇人妻久久综合中文| 在线观看免费日韩欧美大片 | 在线播放无遮挡| 青春草视频在线免费观看| 日韩免费高清中文字幕av| 一级片'在线观看视频| 蜜臀久久99精品久久宅男| .国产精品久久| 久久av网站| 国产色婷婷99| 久久久久久久久久久丰满| 免费日韩欧美在线观看| 精品卡一卡二卡四卡免费| 天堂8中文在线网| 免费播放大片免费观看视频在线观看| av在线观看视频网站免费| 久久人妻熟女aⅴ| 午夜福利视频在线观看免费| 国产成人免费无遮挡视频| av天堂久久9| 九九久久精品国产亚洲av麻豆| 人妻一区二区av| 我要看黄色一级片免费的| 卡戴珊不雅视频在线播放| 五月玫瑰六月丁香| 一本—道久久a久久精品蜜桃钙片| 国产黄频视频在线观看| 久久久亚洲精品成人影院| 免费观看在线日韩| 日韩av在线免费看完整版不卡| av.在线天堂| 久久精品国产亚洲网站| a级毛片免费高清观看在线播放| 国产日韩一区二区三区精品不卡 | 高清黄色对白视频在线免费看| 在线播放无遮挡| 精品亚洲成国产av| 内地一区二区视频在线| 久久精品久久久久久久性| 精品人妻熟女av久视频| 欧美日韩av久久| 国产精品一区www在线观看| 色视频在线一区二区三区| 免费av中文字幕在线| 日韩av免费高清视频| 又大又黄又爽视频免费| 最新中文字幕久久久久| 制服诱惑二区| 又大又黄又爽视频免费| 一本色道久久久久久精品综合| 免费播放大片免费观看视频在线观看| 欧美精品一区二区免费开放| 菩萨蛮人人尽说江南好唐韦庄| av天堂久久9| 国模一区二区三区四区视频| 高清欧美精品videossex| 久久国内精品自在自线图片| 一区二区三区免费毛片| 一区二区三区免费毛片| 五月天丁香电影| 性色avwww在线观看| av电影中文网址| 在线观看免费高清a一片| 成人无遮挡网站| 国产精品久久久久成人av| 91久久精品国产一区二区成人| 五月玫瑰六月丁香| 免费不卡的大黄色大毛片视频在线观看| av免费在线看不卡| av天堂久久9| 国产高清国产精品国产三级| 国产欧美另类精品又又久久亚洲欧美| 爱豆传媒免费全集在线观看| 国产亚洲欧美精品永久| 国产片内射在线| 亚洲婷婷狠狠爱综合网| 中文乱码字字幕精品一区二区三区| 国产精品熟女久久久久浪| 久久国内精品自在自线图片| 免费黄色在线免费观看| 亚洲av综合色区一区| 婷婷色av中文字幕| 久久久久久人妻| 乱码一卡2卡4卡精品| 成人综合一区亚洲| 国产成人一区二区在线| 免费看光身美女| 国产成人精品在线电影| 大话2 男鬼变身卡| 国产熟女欧美一区二区| 22中文网久久字幕| av网站免费在线观看视频| 亚洲激情五月婷婷啪啪| 男女免费视频国产| 亚洲一区二区三区欧美精品| 日韩中字成人| 亚洲色图 男人天堂 中文字幕 | 我要看黄色一级片免费的| 最近的中文字幕免费完整| 成年人午夜在线观看视频| 国产精品久久久久久精品电影小说| 欧美xxxx性猛交bbbb| 九色亚洲精品在线播放| 国产高清有码在线观看视频| 插阴视频在线观看视频| 日韩一区二区三区影片| 黑人巨大精品欧美一区二区蜜桃 | 亚洲精品国产av成人精品| 51国产日韩欧美| 亚洲不卡免费看| 狠狠精品人妻久久久久久综合| 夜夜爽夜夜爽视频| 国产精品99久久久久久久久| 在线天堂最新版资源| 午夜福利网站1000一区二区三区| 最新的欧美精品一区二区| 极品人妻少妇av视频| 国产精品欧美亚洲77777| 激情五月婷婷亚洲| 人人澡人人妻人| 成年人午夜在线观看视频| 久久精品国产亚洲网站| 国产乱来视频区| 日本91视频免费播放| 日本黄色片子视频| 丝袜在线中文字幕| 亚洲图色成人| 免费黄网站久久成人精品| 秋霞伦理黄片| 18禁在线播放成人免费| 午夜久久久在线观看| 春色校园在线视频观看| 亚洲av成人精品一区久久| 久久亚洲国产成人精品v| 九九在线视频观看精品| 男女边摸边吃奶| 九九在线视频观看精品| 高清黄色对白视频在线免费看| 成人手机av| 国产伦精品一区二区三区视频9| 亚洲欧洲日产国产| 亚洲色图综合在线观看| av电影中文网址| 啦啦啦啦在线视频资源| 亚洲av成人精品一区久久| 午夜福利,免费看| 男女免费视频国产| 观看av在线不卡| 自拍欧美九色日韩亚洲蝌蚪91| 天美传媒精品一区二区| 亚洲欧美清纯卡通| 日本欧美国产在线视频| 精品亚洲成国产av| 国产精品一区二区三区四区免费观看| 高清午夜精品一区二区三区| 亚洲av中文av极速乱| 3wmmmm亚洲av在线观看| 国产成人免费观看mmmm| 久久久久久久精品精品| 热re99久久国产66热| 国产乱来视频区| 国产高清三级在线| 在线精品无人区一区二区三| 全区人妻精品视频| 亚洲av中文av极速乱| 国产精品秋霞免费鲁丝片| 成人18禁高潮啪啪吃奶动态图 | 9色porny在线观看| 一二三四中文在线观看免费高清| 亚洲欧美日韩另类电影网站| 丝袜喷水一区| 99精国产麻豆久久婷婷| 成人手机av| 一级片'在线观看视频| 国产免费一区二区三区四区乱码| 2021少妇久久久久久久久久久| 亚洲精品中文字幕在线视频| 少妇被粗大的猛进出69影院 | 精品久久国产蜜桃| 欧美少妇被猛烈插入视频| 午夜老司机福利剧场| 中文乱码字字幕精品一区二区三区| 国产深夜福利视频在线观看| kizo精华| 亚洲av日韩在线播放| 亚洲国产av新网站| 久久人人爽av亚洲精品天堂| av国产久精品久网站免费入址| 国产日韩欧美视频二区| 男女无遮挡免费网站观看| 久久久国产精品麻豆| 五月伊人婷婷丁香| 日本av手机在线免费观看| 最后的刺客免费高清国语| 午夜激情福利司机影院| 黄色欧美视频在线观看| 大陆偷拍与自拍| 国产永久视频网站| 国产精品人妻久久久影院| 精品少妇内射三级| 曰老女人黄片| 久久久久久久久久久丰满| 看免费成人av毛片| 最新的欧美精品一区二区| 久久久国产一区二区| 国产在视频线精品| 蜜桃久久精品国产亚洲av| 在线观看免费视频网站a站| 97超碰精品成人国产| 成年人免费黄色播放视频| 国产伦理片在线播放av一区| 亚洲国产日韩一区二区| 午夜日本视频在线| 免费久久久久久久精品成人欧美视频 | 下体分泌物呈黄色| 免费看不卡的av| 我的老师免费观看完整版| 亚洲综合色网址| 热re99久久精品国产66热6| kizo精华| 2022亚洲国产成人精品| 青青草视频在线视频观看| 另类精品久久| 亚洲,欧美,日韩| 日韩大片免费观看网站| 亚洲av综合色区一区| 国产高清不卡午夜福利| 久久久久久伊人网av| 精品卡一卡二卡四卡免费| 99久久人妻综合| 亚洲av二区三区四区| 久久精品国产自在天天线| 日日爽夜夜爽网站| 91精品国产九色| 国产探花极品一区二区| 日本vs欧美在线观看视频| 狂野欧美激情性bbbbbb| 免费大片黄手机在线观看| 男人操女人黄网站| 黄片播放在线免费| 亚洲av不卡在线观看| 国产精品.久久久| 搡女人真爽免费视频火全软件| 美女内射精品一级片tv| av在线老鸭窝| 欧美日韩视频高清一区二区三区二| 欧美亚洲日本最大视频资源| 国产av码专区亚洲av| 桃花免费在线播放| 五月天丁香电影| av视频免费观看在线观看| 成人手机av| 91久久精品国产一区二区成人| 国产日韩欧美亚洲二区| 成人毛片a级毛片在线播放| 免费av中文字幕在线| 日韩 亚洲 欧美在线| 成人午夜精彩视频在线观看| 人妻夜夜爽99麻豆av| 成人国语在线视频| 国产免费视频播放在线视频| 国产一级毛片在线| 国产一区二区三区综合在线观看 | 人妻制服诱惑在线中文字幕| 久久久午夜欧美精品| 国产av国产精品国产| 国产日韩欧美在线精品| 国产精品无大码| 免费观看无遮挡的男女| 久久av网站| 国产视频首页在线观看| 亚洲无线观看免费| 午夜免费鲁丝| 久久婷婷青草| 国产伦精品一区二区三区视频9| 亚洲国产精品成人久久小说| 国产探花极品一区二区| 精品久久久久久久久亚洲| 久久 成人 亚洲| 亚洲精品久久成人aⅴ小说 | 一本大道久久a久久精品| 国产欧美亚洲国产| 99九九线精品视频在线观看视频| av黄色大香蕉| 22中文网久久字幕| 亚洲内射少妇av| 亚洲精品国产av蜜桃| 国产高清不卡午夜福利| 色婷婷av一区二区三区视频| a级毛片黄视频| 久久久久国产网址| 日韩三级伦理在线观看| 亚洲av不卡在线观看| 18在线观看网站| 国内精品宾馆在线| 日韩欧美一区视频在线观看| 九九久久精品国产亚洲av麻豆| 狂野欧美激情性bbbbbb| 国产av精品麻豆| 久久久久人妻精品一区果冻| 国产精品麻豆人妻色哟哟久久| 亚洲av国产av综合av卡| 欧美日韩国产mv在线观看视频| 久久久欧美国产精品| 日韩精品有码人妻一区| 午夜福利视频精品| 另类亚洲欧美激情| 夫妻午夜视频| 97在线人人人人妻| 一级,二级,三级黄色视频| 91久久精品国产一区二区成人| 欧美成人午夜免费资源| 日本与韩国留学比较| 最近手机中文字幕大全| 热99久久久久精品小说推荐| 大码成人一级视频| 高清在线视频一区二区三区| 两个人免费观看高清视频| 美女国产视频在线观看| 欧美亚洲日本最大视频资源| 肉色欧美久久久久久久蜜桃| 一级片'在线观看视频| 精品亚洲乱码少妇综合久久| 国产一区二区三区av在线| 麻豆乱淫一区二区| 日韩一本色道免费dvd| 国产老妇伦熟女老妇高清| 最近中文字幕2019免费版| 99九九线精品视频在线观看视频| 国产高清三级在线| 人人妻人人澡人人看| 啦啦啦啦在线视频资源| 高清午夜精品一区二区三区| 成人影院久久| 亚洲欧美一区二区三区国产| 久久久欧美国产精品| 日本vs欧美在线观看视频| 美女视频免费永久观看网站| 狠狠婷婷综合久久久久久88av| 春色校园在线视频观看| 亚洲成色77777| xxxhd国产人妻xxx| 国产成人一区二区在线| 99九九线精品视频在线观看视频| kizo精华| 国产日韩一区二区三区精品不卡 | 不卡视频在线观看欧美| 精品久久国产蜜桃| 日韩一区二区视频免费看| 一边亲一边摸免费视频| 午夜福利视频在线观看免费| 3wmmmm亚洲av在线观看| 简卡轻食公司| 99精国产麻豆久久婷婷| 久久综合国产亚洲精品| 亚洲丝袜综合中文字幕| 国产精品国产三级国产av玫瑰| 91精品国产九色| 欧美日韩在线观看h| 简卡轻食公司| 日本免费在线观看一区| 男男h啪啪无遮挡| 91精品一卡2卡3卡4卡| 国产男人的电影天堂91| 精品亚洲乱码少妇综合久久| 中文字幕人妻丝袜制服| 亚洲国产精品一区三区| 午夜视频国产福利| www.色视频.com| 色婷婷久久久亚洲欧美| 国产毛片在线视频| 午夜激情久久久久久久| 99久久人妻综合| 国产精品一国产av| 男人爽女人下面视频在线观看| 99国产综合亚洲精品| 国产极品粉嫩免费观看在线 | 亚洲欧美日韩卡通动漫| 一个人免费看片子| 国产 一区精品| 中文乱码字字幕精品一区二区三区| 久久亚洲国产成人精品v| 久久精品国产自在天天线| 欧美xxxx性猛交bbbb| 最近中文字幕2019免费版| 国产成人av激情在线播放 | 久久久欧美国产精品| 国产成人91sexporn| 水蜜桃什么品种好| 波野结衣二区三区在线| 欧美日韩视频高清一区二区三区二| 热99久久久久精品小说推荐| 久久影院123| 欧美老熟妇乱子伦牲交| 久久久精品94久久精品| 80岁老熟妇乱子伦牲交| 一级a做视频免费观看| 精品视频人人做人人爽| xxx大片免费视频| 蜜桃久久精品国产亚洲av| 国产精品无大码| 午夜福利影视在线免费观看| 亚洲国产日韩一区二区| 国产精品女同一区二区软件| 国产伦理片在线播放av一区| 日本猛色少妇xxxxx猛交久久| av在线老鸭窝| 一级毛片 在线播放| 涩涩av久久男人的天堂| 国产伦精品一区二区三区视频9| 婷婷成人精品国产| 91国产中文字幕| 亚洲国产最新在线播放| 欧美老熟妇乱子伦牲交| 国产成人精品婷婷| av线在线观看网站| 国产黄频视频在线观看| 99热6这里只有精品| 三级国产精品欧美在线观看| 一区二区三区免费毛片| a级片在线免费高清观看视频| 精品一区二区免费观看| 欧美精品亚洲一区二区| 插逼视频在线观看| 国产成人a∨麻豆精品| 夫妻性生交免费视频一级片| 一区二区三区乱码不卡18| 全区人妻精品视频| 国内精品宾馆在线| 国产永久视频网站| 少妇 在线观看| 久久久久久久久久久久大奶| 日本欧美国产在线视频| 精品国产一区二区三区久久久樱花| 久久99热这里只频精品6学生| 国产成人午夜福利电影在线观看| 久久久久久久久大av| 欧美 日韩 精品 国产| 国产女主播在线喷水免费视频网站| 国产成人a∨麻豆精品| 亚洲,欧美,日韩| 丰满迷人的少妇在线观看| h视频一区二区三区| 精品久久蜜臀av无| 乱人伦中国视频| 欧美bdsm另类| 18禁观看日本| 免费观看a级毛片全部| 一边摸一边做爽爽视频免费| 中文天堂在线官网| 我要看黄色一级片免费的| videosex国产| 亚洲,一卡二卡三卡| 一个人看视频在线观看www免费| 久久久精品94久久精品| 免费人成在线观看视频色| 成人国语在线视频| 飞空精品影院首页| 国产 一区精品| 婷婷色av中文字幕| 一区二区三区精品91| 欧美日韩av久久| 国产成人精品无人区| 黑人欧美特级aaaaaa片| 99久久精品一区二区三区| 各种免费的搞黄视频| 亚洲av成人精品一区久久| av免费观看日本| 丝袜在线中文字幕| 欧美亚洲 丝袜 人妻 在线| 99热全是精品| 亚洲精品视频女| 男人爽女人下面视频在线观看| 日本-黄色视频高清免费观看| 高清欧美精品videossex| 黑丝袜美女国产一区| 国产精品久久久久久av不卡| 少妇高潮的动态图| 青青草视频在线视频观看| 91精品一卡2卡3卡4卡| 欧美日韩视频高清一区二区三区二| 亚洲精品一二三| 久久人人爽人人爽人人片va| 精品久久久噜噜| 18禁裸乳无遮挡动漫免费视频| 久久精品国产鲁丝片午夜精品| 制服诱惑二区| 亚洲五月色婷婷综合| 久久免费观看电影| 秋霞在线观看毛片| 免费高清在线观看视频在线观看| 久久人人爽人人片av| 亚洲激情五月婷婷啪啪| 亚洲精品成人av观看孕妇| 久久99精品国语久久久| 亚洲成人一二三区av| 国产精品一国产av| 欧美亚洲日本最大视频资源| 男人操女人黄网站| 欧美 日韩 精品 国产| 日韩精品有码人妻一区| 国产精品久久久久久久久免| 日本与韩国留学比较| 如何舔出高潮| 人人妻人人澡人人看| av线在线观看网站| 十分钟在线观看高清视频www| tube8黄色片| 国产男女内射视频| 亚洲三级黄色毛片| 大码成人一级视频| 国产色婷婷99| 国产男人的电影天堂91| 精品视频人人做人人爽| 国产亚洲最大av| 精品人妻熟女毛片av久久网站| 国产成人a∨麻豆精品| 高清午夜精品一区二区三区| 在现免费观看毛片| 午夜福利在线观看免费完整高清在| 精品少妇久久久久久888优播| 赤兔流量卡办理| 久久影院123| 看免费成人av毛片| 久久青草综合色| 色网站视频免费| 久久韩国三级中文字幕| 中国美白少妇内射xxxbb|