• <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
    致秋天的花楸樹
    三级男女做爰猛烈吃奶摸视频| 可以在线观看的亚洲视频| 99久国产av精品国产电影| 非洲黑人性xxxx精品又粗又长| 蜜桃久久精品国产亚洲av| 午夜免费男女啪啪视频观看| 国内精品美女久久久久久| 久久婷婷人人爽人人干人人爱| 国产黄色小视频在线观看| 精品人妻熟女av久视频| 看片在线看免费视频| 免费观看在线日韩| 色吧在线观看| 好男人视频免费观看在线| 亚洲一区高清亚洲精品| 国产精品电影一区二区三区| 悠悠久久av| 99久久中文字幕三级久久日本| 69av精品久久久久久| 99热全是精品| 久久久久免费精品人妻一区二区| 国内少妇人妻偷人精品xxx网站| 免费电影在线观看免费观看| 免费搜索国产男女视频| 成人av在线播放网站| 春色校园在线视频观看| 熟妇人妻久久中文字幕3abv| 赤兔流量卡办理| 一个人看的www免费观看视频| 精品人妻偷拍中文字幕| 国产乱人视频| 蜜臀久久99精品久久宅男| 我的女老师完整版在线观看| 亚洲欧洲日产国产| 精品免费久久久久久久清纯| 一区二区三区免费毛片| 少妇人妻精品综合一区二区 | 99热全是精品| 亚洲欧美日韩高清在线视频| 成人永久免费在线观看视频| 寂寞人妻少妇视频99o| 欧美一区二区精品小视频在线| 又粗又爽又猛毛片免费看| 色噜噜av男人的天堂激情| 久久久久久国产a免费观看| 日本免费a在线| 丝袜美腿在线中文| 三级毛片av免费| 中文字幕免费在线视频6| 欧美最黄视频在线播放免费| 国产成人午夜福利电影在线观看| 亚洲成人av在线免费| 在线观看66精品国产| 久久精品国产亚洲av涩爱 | 亚洲欧美日韩卡通动漫| 能在线免费看毛片的网站| 久久九九热精品免费| 午夜免费激情av| 搡女人真爽免费视频火全软件| 午夜a级毛片| 国产私拍福利视频在线观看| 少妇的逼好多水| 欧美一区二区亚洲| 99久国产av精品国产电影| 日韩 亚洲 欧美在线| 91精品国产九色| 欧美激情国产日韩精品一区| 男女下面进入的视频免费午夜| 日韩成人av中文字幕在线观看| 亚洲一区高清亚洲精品| .国产精品久久| 久久亚洲精品不卡| 一级毛片我不卡| 又爽又黄a免费视频| 国产白丝娇喘喷水9色精品| 97人妻精品一区二区三区麻豆| 久久综合国产亚洲精品| 老司机福利观看| 天天一区二区日本电影三级| 成人国产麻豆网| av卡一久久| 在线播放国产精品三级| 身体一侧抽搐| 狂野欧美激情性xxxx在线观看| 欧美zozozo另类| 免费av毛片视频| 成年免费大片在线观看| 可以在线观看的亚洲视频| 中出人妻视频一区二区| 级片在线观看| 啦啦啦观看免费观看视频高清| 亚洲欧美成人精品一区二区| 欧美激情国产日韩精品一区| av福利片在线观看| 你懂的网址亚洲精品在线观看 | 三级经典国产精品| 男人舔奶头视频| 26uuu在线亚洲综合色| 狂野欧美激情性xxxx在线观看| 国产淫片久久久久久久久| 日韩欧美国产在线观看| 欧美最新免费一区二区三区| 精品久久久久久久久亚洲| 国产一区二区在线观看日韩| 一级二级三级毛片免费看| 热99re8久久精品国产| 国产午夜精品论理片| 亚洲国产精品sss在线观看| 岛国在线免费视频观看| 国产精品蜜桃在线观看 | 九九热线精品视视频播放| 少妇人妻精品综合一区二区 | 精品久久国产蜜桃| 欧美色欧美亚洲另类二区| 免费看光身美女| 搡女人真爽免费视频火全软件| 国产成人影院久久av| 免费观看精品视频网站| 欧美一区二区精品小视频在线| 日韩强制内射视频| 国产91av在线免费观看| 久久精品91蜜桃| 99热网站在线观看| 在现免费观看毛片| 精品久久久久久久人妻蜜臀av| 成年版毛片免费区| 亚洲四区av| 美女国产视频在线观看| 亚洲最大成人av| 1024手机看黄色片| 国产在线精品亚洲第一网站| 尤物成人国产欧美一区二区三区| 久久久欧美国产精品| 人妻夜夜爽99麻豆av| 亚洲综合色惰| 能在线免费看毛片的网站| 亚洲天堂国产精品一区在线| 欧美日本视频| 国内精品美女久久久久久| 天堂√8在线中文| 禁无遮挡网站| 中文字幕免费在线视频6| 亚洲国产精品合色在线| 最好的美女福利视频网| 成人漫画全彩无遮挡| 国产精品一区www在线观看| 99热这里只有是精品50| .国产精品久久| 久久99精品国语久久久| 亚洲av免费在线观看| 一边亲一边摸免费视频| 69人妻影院| 国产成人a∨麻豆精品| 久久久国产成人精品二区| 国内精品久久久久精免费| 欧美不卡视频在线免费观看| 亚洲欧美精品专区久久| or卡值多少钱| 国产淫片久久久久久久久| 国产在线精品亚洲第一网站| 亚洲经典国产精华液单| eeuss影院久久| 中文字幕熟女人妻在线| 午夜精品在线福利| 乱码一卡2卡4卡精品| 国产黄片美女视频| 麻豆国产97在线/欧美| 青春草国产在线视频 | 国产精品av视频在线免费观看| 亚洲av免费高清在线观看| 午夜视频国产福利| 国产高清激情床上av| 精品一区二区三区视频在线| 少妇熟女aⅴ在线视频| 色吧在线观看| 久久草成人影院| 小说图片视频综合网站| 免费看美女性在线毛片视频| 18禁在线播放成人免费| 嘟嘟电影网在线观看| 国产精品一区二区在线观看99 | 日本黄大片高清| 狂野欧美激情性xxxx在线观看| 精品久久久久久成人av| 国产精品综合久久久久久久免费| 久久韩国三级中文字幕| 久久亚洲国产成人精品v| 久久精品夜夜夜夜夜久久蜜豆| 日韩亚洲欧美综合| 伦精品一区二区三区| 国产蜜桃级精品一区二区三区| 亚洲成人久久爱视频| 岛国毛片在线播放| 欧美一区二区国产精品久久精品| 国内久久婷婷六月综合欲色啪| 99久久久亚洲精品蜜臀av| а√天堂www在线а√下载| 国产精品永久免费网站| 国产av在哪里看| 一级毛片我不卡| 亚洲成av人片在线播放无| 日韩欧美一区二区三区在线观看| 女人被狂操c到高潮| 婷婷亚洲欧美| 狂野欧美激情性xxxx在线观看| 亚洲图色成人| 日本与韩国留学比较| 能在线免费看毛片的网站| 黄色欧美视频在线观看| 黄色配什么色好看| 日本黄大片高清| 欧美在线一区亚洲| av卡一久久| 美女 人体艺术 gogo| 欧美日本视频| 日产精品乱码卡一卡2卡三| 青春草亚洲视频在线观看| 最近的中文字幕免费完整| a级一级毛片免费在线观看| 22中文网久久字幕| 久久久精品94久久精品| 精品无人区乱码1区二区| 人人妻人人看人人澡| 91午夜精品亚洲一区二区三区| 久久精品综合一区二区三区| 国产精品一区www在线观看| 小蜜桃在线观看免费完整版高清| 亚洲精品久久久久久婷婷小说 | 国产高潮美女av| 天美传媒精品一区二区| 亚洲丝袜综合中文字幕| 亚洲美女视频黄频| 蜜桃亚洲精品一区二区三区| 91午夜精品亚洲一区二区三区| 好男人视频免费观看在线| 欧美日本视频| 欧美精品国产亚洲| 午夜久久久久精精品| 亚洲精品久久久久久婷婷小说 | 久久精品国产亚洲网站| 国内精品美女久久久久久| 欧美激情在线99| 可以在线观看毛片的网站| 国产精品av视频在线免费观看| 久久人人爽人人片av| 小说图片视频综合网站| 亚洲国产精品成人综合色| 国产午夜精品一二区理论片| АⅤ资源中文在线天堂| 亚洲天堂国产精品一区在线| 午夜a级毛片| 韩国av在线不卡| 欧美日韩综合久久久久久| 欧美一区二区亚洲| 人人妻人人看人人澡| 国产精品国产三级国产av玫瑰| 麻豆av噜噜一区二区三区| 身体一侧抽搐| 国产精华一区二区三区| 成人亚洲欧美一区二区av| 久久久久久大精品| 99久久成人亚洲精品观看| 欧美日韩在线观看h| 国内精品宾馆在线| 在线a可以看的网站| 欧美性猛交╳xxx乱大交人| 国产色爽女视频免费观看| 看非洲黑人一级黄片| 熟女人妻精品中文字幕| 国产午夜福利久久久久久| 一夜夜www| 少妇的逼水好多| 美女黄网站色视频| 精品日产1卡2卡| 天天躁夜夜躁狠狠久久av| 天天躁日日操中文字幕| 国内精品一区二区在线观看| 偷拍熟女少妇极品色| 免费人成在线观看视频色| 色综合亚洲欧美另类图片| 性欧美人与动物交配| 天堂中文最新版在线下载 | 黄色欧美视频在线观看| 99热只有精品国产| 内地一区二区视频在线| av天堂在线播放| 亚洲国产色片| 99久久九九国产精品国产免费| 综合色av麻豆| 亚洲av成人精品一区久久| 搡女人真爽免费视频火全软件| 夜夜夜夜夜久久久久| 国产成人精品久久久久久| 欧美成人a在线观看| 国产免费男女视频| 综合色丁香网| 亚洲成av人片在线播放无| 亚洲av免费在线观看| 日本三级黄在线观看| 人人妻人人澡欧美一区二区| 乱系列少妇在线播放| 欧美最新免费一区二区三区| 嫩草影院入口| 免费av不卡在线播放| 校园人妻丝袜中文字幕| 亚洲欧美成人精品一区二区| 狂野欧美激情性xxxx在线观看| 青春草视频在线免费观看| 午夜a级毛片| 亚洲精品456在线播放app| 日本五十路高清| 好男人在线观看高清免费视频| 一个人看的www免费观看视频| 一个人免费在线观看电影| a级一级毛片免费在线观看| 国产高清视频在线观看网站| 直男gayav资源| 成人三级黄色视频| 51国产日韩欧美| 国产午夜福利久久久久久| 中国美女看黄片| www.av在线官网国产| 国产视频首页在线观看| 青春草视频在线免费观看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 日韩欧美 国产精品| 一本久久精品| 免费观看在线日韩| 久久久色成人| 一区二区三区高清视频在线| 女人十人毛片免费观看3o分钟| 欧美日本视频| 国产精品女同一区二区软件| 久久久久免费精品人妻一区二区| 欧美一区二区亚洲| 一个人看视频在线观看www免费| 中国国产av一级| 少妇的逼水好多| 亚洲欧洲日产国产| 久久久久久久久中文| 国产一区二区激情短视频| 国产成人a区在线观看| 色吧在线观看| 亚洲在线观看片| 国产精品久久久久久精品电影| 免费黄网站久久成人精品| 国产真实伦视频高清在线观看| 国产视频内射| 深夜a级毛片| 亚洲一区高清亚洲精品| 97人妻精品一区二区三区麻豆| 国产69精品久久久久777片| 久久99精品国语久久久| 欧美丝袜亚洲另类| 国产极品天堂在线| 啦啦啦啦在线视频资源| 18禁黄网站禁片免费观看直播| 少妇熟女欧美另类| 蜜臀久久99精品久久宅男| 99热6这里只有精品| 在线观看av片永久免费下载| 99热全是精品| 少妇人妻一区二区三区视频| 国产精品蜜桃在线观看 | 热99re8久久精品国产| 国产乱人视频| 青春草国产在线视频 | 亚洲欧美中文字幕日韩二区| 国产精品综合久久久久久久免费| 一级毛片我不卡| 亚洲无线观看免费| 最新中文字幕久久久久| 超碰av人人做人人爽久久| 精品99又大又爽又粗少妇毛片| 成年女人看的毛片在线观看| 嫩草影院精品99| 国产高清三级在线| 97人妻精品一区二区三区麻豆| 久久久国产成人免费| 亚洲无线观看免费| 国产成年人精品一区二区| 嘟嘟电影网在线观看| 国产精品综合久久久久久久免费| 亚洲精品影视一区二区三区av| 亚洲欧美成人精品一区二区| 高清日韩中文字幕在线| 久久久久性生活片| 亚洲国产高清在线一区二区三| 日日摸夜夜添夜夜爱| 成人欧美大片| 菩萨蛮人人尽说江南好唐韦庄 | 日韩高清综合在线| 国产精品不卡视频一区二区| 日韩欧美三级三区| 一卡2卡三卡四卡精品乱码亚洲| 日韩精品青青久久久久久| 亚洲精品影视一区二区三区av| 色5月婷婷丁香| 国产精品人妻久久久影院| 尤物成人国产欧美一区二区三区| 一区福利在线观看| 成人二区视频| 日日啪夜夜撸| 乱系列少妇在线播放| 91午夜精品亚洲一区二区三区| 中文字幕精品亚洲无线码一区| 国产亚洲5aaaaa淫片| 麻豆av噜噜一区二区三区| 少妇人妻一区二区三区视频| 亚洲精品乱码久久久v下载方式| www日本黄色视频网| 欧美高清性xxxxhd video| 99热网站在线观看| 免费观看精品视频网站| 舔av片在线| 91狼人影院| 两个人视频免费观看高清| 国内精品一区二区在线观看| 一级毛片我不卡| 亚洲中文字幕一区二区三区有码在线看| 亚洲欧美日韩高清专用| 99国产极品粉嫩在线观看| 三级毛片av免费| 中文欧美无线码| 晚上一个人看的免费电影| av国产免费在线观看| 特大巨黑吊av在线直播| av国产免费在线观看| 日本在线视频免费播放| 99久久中文字幕三级久久日本| 黄色欧美视频在线观看| 久久久精品欧美日韩精品| 国产探花极品一区二区| 亚洲精品影视一区二区三区av| 在线天堂最新版资源| 国产片特级美女逼逼视频| 两个人的视频大全免费| 一级黄片播放器| 中文欧美无线码| 亚洲久久久久久中文字幕| 免费电影在线观看免费观看| a级一级毛片免费在线观看| 一级二级三级毛片免费看| 搞女人的毛片| 在线观看午夜福利视频| 亚洲一区二区三区色噜噜| 热99re8久久精品国产| 国产午夜精品久久久久久一区二区三区| 国产单亲对白刺激| 高清日韩中文字幕在线| 天堂av国产一区二区熟女人妻| 国产精品女同一区二区软件| 黄色日韩在线| 国内精品宾馆在线| 国产av在哪里看| 午夜精品国产一区二区电影 | 晚上一个人看的免费电影| 久久久成人免费电影| 亚洲人成网站高清观看| 春色校园在线视频观看| 色哟哟哟哟哟哟| 国产真实乱freesex| 国产av在哪里看| 精品99又大又爽又粗少妇毛片| 国产精品,欧美在线| 国产激情偷乱视频一区二区| 亚洲av第一区精品v没综合| 免费无遮挡裸体视频| 高清毛片免费观看视频网站| 国内少妇人妻偷人精品xxx网站| 色视频www国产| 久久精品人妻少妇| 女人被狂操c到高潮| 日韩欧美一区二区三区在线观看| av专区在线播放| 日韩成人av中文字幕在线观看| 欧美最新免费一区二区三区| 国产真实乱freesex| 91aial.com中文字幕在线观看| 99久久人妻综合| 亚洲av男天堂| 国产激情偷乱视频一区二区| av视频在线观看入口| 一级黄色大片毛片| 国产成人aa在线观看| 哪个播放器可以免费观看大片| 日本在线视频免费播放| 国产一区二区激情短视频| 不卡视频在线观看欧美| 亚洲精品久久国产高清桃花| 亚洲成人久久性| 免费在线观看成人毛片| 亚洲欧美成人综合另类久久久 | 熟女电影av网| 麻豆成人午夜福利视频| 亚洲成人中文字幕在线播放| 老女人水多毛片| 国产精品,欧美在线| 中文字幕av在线有码专区| 日日摸夜夜添夜夜爱| 身体一侧抽搐| 亚洲精品456在线播放app| 丰满人妻一区二区三区视频av| 久久这里有精品视频免费| 日日干狠狠操夜夜爽| 日本黄色视频三级网站网址| 非洲黑人性xxxx精品又粗又长| 午夜a级毛片| 欧美最黄视频在线播放免费| 国产高清有码在线观看视频| 国内精品美女久久久久久| 啦啦啦啦在线视频资源| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 99久久精品热视频| 少妇的逼水好多| 久久久久网色| 少妇人妻一区二区三区视频| 男女那种视频在线观看| 我要搜黄色片| 天堂中文最新版在线下载 | 欧美精品一区二区大全| 国产精品伦人一区二区| 免费电影在线观看免费观看| 久久久久免费精品人妻一区二区| 插阴视频在线观看视频| 2021天堂中文幕一二区在线观| 中文字幕人妻熟人妻熟丝袜美| 亚洲内射少妇av| 日韩国内少妇激情av| 一本久久精品| 国产午夜精品久久久久久一区二区三区| 中文字幕人妻熟人妻熟丝袜美| 日韩强制内射视频| 精品不卡国产一区二区三区| a级毛色黄片| 国产精品一及| 麻豆成人av视频| 日韩高清综合在线| 久久人人精品亚洲av| 又爽又黄a免费视频| 亚洲精品粉嫩美女一区| 大香蕉久久网| 亚洲人与动物交配视频| www.色视频.com| 在线播放无遮挡| 一本久久精品| 欧美日韩在线观看h| 99热这里只有是精品50| 69人妻影院| 边亲边吃奶的免费视频| 99热这里只有精品一区| 成人特级av手机在线观看| 免费av不卡在线播放| 精品无人区乱码1区二区| 美女国产视频在线观看| 久久久久免费精品人妻一区二区| 国产亚洲精品久久久com| 美女大奶头视频| 青春草视频在线免费观看| 日日摸夜夜添夜夜爱| 国产乱人视频| 尤物成人国产欧美一区二区三区| 99九九线精品视频在线观看视频| 欧美最新免费一区二区三区| 久99久视频精品免费| 国产精品无大码| 国产淫片久久久久久久久| 国产亚洲5aaaaa淫片| 久久久国产成人免费| 日韩av不卡免费在线播放| 在线免费观看不下载黄p国产| 久久久午夜欧美精品| 国产精品伦人一区二区| 少妇丰满av| 人妻久久中文字幕网| 色综合亚洲欧美另类图片| 精品久久久噜噜| 亚洲精品国产成人久久av| 国产精品99久久久久久久久| 欧美性感艳星| 日韩,欧美,国产一区二区三区 | 国产精品嫩草影院av在线观看| 日韩一本色道免费dvd| 国产精品伦人一区二区| 久久精品国产亚洲av香蕉五月| 亚洲av成人精品一区久久| 一级毛片电影观看 | 久久久久久久亚洲中文字幕| 97在线视频观看| 中国美白少妇内射xxxbb| 国产黄片视频在线免费观看| 欧美最新免费一区二区三区| 久久久久网色| 久久人人精品亚洲av| 久久久久久久久久久丰满| 黑人高潮一二区| kizo精华| 99国产精品一区二区蜜桃av| 成人毛片60女人毛片免费| 桃色一区二区三区在线观看| 哪个播放器可以免费观看大片| 亚洲三级黄色毛片| 日韩av在线大香蕉| 只有这里有精品99| 免费观看精品视频网站| 精品久久久久久久末码| 亚洲精品自拍成人| 亚洲av成人精品一区久久| 少妇猛男粗大的猛烈进出视频 | 蜜桃久久精品国产亚洲av| 我要看日韩黄色一级片| av在线蜜桃| 高清毛片免费看| 欧美日韩乱码在线| 日韩一本色道免费dvd|