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

    2種羌活藥用植物的葉綠體基因組密碼子使用偏性分析

    2022-03-04 00:36:22王琳璇王高烽谷威趙志敏馬永貴
    安徽農(nóng)業(yè)科學(xué) 2022年3期
    關(guān)鍵詞:羌活

    王琳璇 王高烽 谷威 趙志敏 馬永貴

    摘要 對(duì)羌活(Notopterygium incisum)與寬葉羌活(N.franchetii)的葉綠體編碼基因的密碼子使用偏性(codon usage bias)及其影響因素進(jìn)行了分析。結(jié)果表明:羌活和寬葉羌活的葉綠體基因組總ENC值分別為49.79和49.70,偏性較弱,密碼子第3位堿基傾向于使用A/T。PR2-plot分析表明,選擇壓力和突變都是影響2種羌活葉綠體密碼子使用偏性的主要因素;而中性繪圖分析和ENC-plot分析表明選擇壓力是影響2種羌活葉綠體密碼子使用偏性的主要因素。相對(duì)同義密碼子使用度(RSCU)分析確定了2種羌活中使用頻率較高的30個(gè)高頻密碼子,幾乎都以A或U結(jié)尾,只有UUG以G結(jié)尾。羌活的最優(yōu)密碼子9個(gè),寬葉羌活的最優(yōu)密碼子7個(gè)。2種羌活有5個(gè)最優(yōu)密碼子相同,分別是UUG、CAU、CGU、AGU、GGA。該研究結(jié)果為后續(xù)藥用羌活植物的基因表達(dá)、蛋白質(zhì)翻譯和遺傳育種研究提供了研究依據(jù)和基礎(chǔ)數(shù)據(jù)。

    關(guān)鍵詞 羌活;寬葉羌活;葉綠體基因組;密碼子使用偏性;最優(yōu)密碼子

    中圖分類號(hào) Q 943.2? 文獻(xiàn)標(biāo)識(shí)碼 A? 文章編號(hào) 0517-6611(2022)03-0095-07

    doi:10.3969/j.issn.0517-6611.2022.03.025

    Codon Usage Bias in Chloroplast Genomes of Two Medicinal Species of Notopterygium H. de Boissieu

    WANG Lin-xuan, WANG Gao-feng, GU Wei et al

    (College of Life Sciences, Qinghai Normal University, Xining, Qinghai 810000)

    Abstract This study analyzed the codon usage bias and affecting factors of the protein coding genes (CDS) of the chloroplast genomes of two medicinal plants, Notopterygium incisum and N. franchetii. The results showed that the chloroplast genomes of N. incisum and N. franchetii had a total ENC value of 49.79 and 49.70, respectively, indicating weak codon usage bias. The third base of the codon of the two chloroplast genomes preferred to use A/T. The analysis of PR2-plot showed that selection pressure and nutation were the main factors affecting the bias of codon usage in N. incisum and N. franchetii CDS, while neutrality plot and ENC-plot analyses both illustrated selection pressure was the main factors affecting codon usage bias in the two species. Relative Synonymous Codon Usage (RSCU) analysis determined 30 high-frequency codons used in N. incisum and N. franchetii, almost all ending with A or U, with only UUG ends with G. There were nine optimal codons identified in N. incisum, and seven optimal codons found in N. franchetii. The five optimal codons shared by N. incisum and N. franchetii were namely UUG, CAU, CGU, AGU, GGA. The results provide important basal data for the future study of gene expression, protein translation and molecular breeding of medicinal Notopterygium plants.

    Key words Notopterygium incisum;Notopterygium franchetii;Chloroplast genome;Codon usage bias;Optimal codons

    基金項(xiàng)目 青海省基礎(chǔ)研究計(jì)劃項(xiàng)目(2018-ZJ-741)。

    作者簡(jiǎn)介 王琳璇(1994—),女,河北廊坊人,碩士研究生,研究方向:生物化學(xué)與分子生物學(xué)。

    通信作者,教授,博士,從事動(dòng)植物分子進(jìn)化與適應(yīng)研究。

    收稿日期 2021-04-03

    密碼子在蛋白質(zhì)合成過程中起著傳遞遺傳信息的關(guān)鍵作用,在大多數(shù)情況下,編碼同一氨基酸的同一密碼子的使用頻率是不同的,這一現(xiàn)象稱為密碼子使用偏性(codon usage bias)。某一物種或某一基因通常傾向于使用一種或幾種特定的同義密碼子,這些密碼子被稱為最優(yōu)密碼子(optimal codon)。影響密碼子使用偏性的因素很多,最主要的因素是定向突變和中性選擇[1-3],還包括GC含量、同義替換率[4]、tRNA豐度[5]、有效翻譯起始位點(diǎn)的選擇[6]、密碼子親水性和DNA復(fù)制起始位點(diǎn)[7]、基因長(zhǎng)度[8]和表達(dá)水平[9]等。在不同物種中,突變和選擇的影響強(qiáng)度不同。例如在藍(lán)藻中,影響密碼子使用偏性較大的是GC組成與選擇壓力[10],但在果蠅中則是突變對(duì)的密碼子使用偏性影響較大[3]。

    葉綠體是植物進(jìn)行光合作用的器官,也參與到植物的氨基酸合成等通路中。作為植物重要的細(xì)胞器,葉綠體有自己一套獨(dú)特的遺傳信息[11]。針對(duì)葉綠體基因組的密碼子使用偏性研究,可以幫助了解葉綠體蛋白質(zhì)編碼基因的使用偏好,為后續(xù)遺傳改造和基因表達(dá)研究提供基礎(chǔ)數(shù)據(jù)。研究發(fā)現(xiàn),乳油木葉綠體基因組密碼子使用偏性更多受突變影響[12],選擇壓力對(duì)其使用偏性影響較小。杜梨葉綠體密碼子使用偏性則主要受選擇的影響[13]。胡曉艷等[14]在對(duì)酸棗葉綠體基因組的研究中獲得了7個(gè)最優(yōu)密碼子。

    羌活屬(Notopterygium Boissieu)植物隸屬于傘形科(Apiaceae)美味芹族,是我國(guó)特有的多年生草本藥用植物,主要分布在我國(guó)西南部高山地區(qū)[15]。據(jù)《中國(guó)植物志(英文版)》記載,該屬包括寬葉羌活、羌活、羽苞羌活、瀾滄羌活、卵葉羌活和細(xì)葉羌活6個(gè)種[16]?!吨袊?guó)藥典》規(guī)定的羌活基源植物為羌活與寬葉羌活,入藥部位為干燥根及根莖,具有解表散寒、祛風(fēng)除濕和止痛的功效。羌活主要產(chǎn)于四川、云南,次產(chǎn)于青海及甘肅;寬葉羌活主產(chǎn)于甘肅、青海,次產(chǎn)于四川、山西和內(nèi)蒙古等地。近年來,由于過度采挖及生境受到嚴(yán)重破壞,野生羌活的自然群體已經(jīng)面臨瀕危,被載入《中國(guó)生物多樣性紅色名錄(高等植物卷)》[17]。羌活作為一味傳統(tǒng)的中藥,富含異歐芹素、羌活醇、阿魏酸等活性物質(zhì),具有抗炎、抗氧化、抗心律失常、抗菌、抗癌細(xì)胞增殖等藥效作用[18]。目前有關(guān)羌活屬葉綠體基因組密碼子使用偏性研究鮮見報(bào)道。該研究擬對(duì)2種可以入藥的羌活屬植物,即羌活和寬葉羌活進(jìn)行葉綠體基因組的密碼子使用模式分析,研究結(jié)果不但可以闡明羌活光合作用相關(guān)基因的使用偏好,也可為羌活基因表達(dá)研究及遺傳育種研究提供理論基礎(chǔ)。

    1 材料與方法

    1.1 數(shù)據(jù)獲取

    羌活樣品采集于青海省果洛藏族自治州達(dá)日縣,99.46° E,33.77° N,海拔3 943.24 m。寬葉羌活樣品采自青海省瑪可河麥浪溝,100.96° E,32.84° N,海拔3 790.78 m。保存于干冰中,送于北京百邁客生物科技有限公司進(jìn)行葉綠體基因組測(cè)序。利用Illumina HiSeq X Ten 平臺(tái)測(cè)序,使用MIRA序列組裝軟件的MIRA v3.4.1.1模塊進(jìn)行,最終得到羌活的葉綠體全基因組序列長(zhǎng)度為 157 050 bp。寬葉羌活的葉綠體全基因組序列長(zhǎng)度為 142 184 bp。測(cè)序結(jié)果提交至NCBI(national center for biotechnology information),登錄號(hào):羌活MW820163、寬葉羌活MW820162。使用Geneious11.1.5 提取編碼基因(coding sequence,CDS)并篩選,剔除小于300 bp的基因,并選擇起始密碼子為ATG的編碼基因進(jìn)行后續(xù)分析,最終羌活獲得51條編碼基因,寬葉羌活獲得50條編碼基因。

    1.2 分析方法

    1.2.1 密碼子相關(guān)參數(shù)計(jì)算。

    使用軟件codonW1.4.2分析獲得的編碼基因,計(jì)算有效密碼子數(shù)(effective number of codon,ENC)。ENC范圍一般在20~61,越接近20,表示密碼子使用偏性越大。使用在線分析軟件CUSP(http://embossgui.sourceforge.net/demo/)對(duì)葉綠體基因組密碼子的GC含量等堿基組成信息進(jìn)行研究,并對(duì)相關(guān)參數(shù)進(jìn)行相關(guān)性分析[19]。

    1.2.2 中性繪圖分析。

    分別以CGC1、CGC2、CGC3表示密碼子第1、2、3位的GC含量,CGC12表示CGC1與CGC2的平均數(shù),以CGC3和CGC12為橫、縱坐標(biāo)繪制散點(diǎn)圖,分析葉綠體基因組密碼子第3位和前2位的相關(guān)性,從而分析影響密碼子偏性的因素[20]。以CGC3、CGC12為橫、縱坐標(biāo)繪制散點(diǎn)圖,如果CGC3和CGC12呈顯著相關(guān),說明密碼子第3位的堿基組成與前2位相似,密碼子偏性主要受突變的影響;如果CGC3和CGC12無顯著相關(guān)關(guān)系,趨勢(shì)線斜率接近0,則說明密碼子第3位與前2位堿基組成差異較大,密碼子堿基相對(duì)保守,這種情況密碼子偏性更多受選擇的影響[21]。

    1.2.3 PR2-plot繪圖分析。

    為進(jìn)一步分析影響2種羌活葉綠體基因組密碼子偏性的因素,以CA3、CT3、CC3、CG3分別表示密碼子第3位A、T、C、G 4種堿基的含量,以CG3/(CG3+CC3)為橫坐標(biāo),CA3/(CA3+CT3)為縱坐標(biāo)繪制散點(diǎn)圖。如果大部分基因分布在圖像上半部,說明密碼子第3位堿基含量A>T,反之則T>A,同理如果大部分基因分布在圖像右半部,說明密碼子第3位堿基含量G>C,反之則C>G。若點(diǎn)在圖中均勻分布,則密碼子偏性完全受突變影響[20-21]。

    1.2.4 ENC-plot繪圖分析。

    以GC3為橫坐標(biāo),以ENC值為縱坐標(biāo)繪制散點(diǎn)圖,并以ENC期望繪制標(biāo)準(zhǔn)曲線。

    ENC期望=2+CGC32+29CGC32+(1-CGC3)2

    標(biāo)準(zhǔn)曲線表示在無選擇壓力的情況下ENC與CGC3的關(guān)系,如果點(diǎn)主要分布于曲線附近,說明密碼子偏性更多受突變的影響,如果點(diǎn)分布于距標(biāo)準(zhǔn)曲線較遠(yuǎn)的位置,則說明密碼子偏性更多受選擇的影響。根據(jù)公式(ENC期望-ENC實(shí)際)/ENC期望制作ENC分布頻數(shù)表,統(tǒng)計(jì)ENC實(shí)際與ENC期望的差距[22]。

    1.2.5 RSCU分析。

    使用CodonW計(jì)算2種羌活葉綠體基因組的RSCU(relative synonymous codon usage)。RSCU是實(shí)際觀察到的某密碼子使用頻率與預(yù)期(該密碼子編碼的氨基酸的同義密碼子使用頻率相等時(shí))使用頻率的比值。理論值與預(yù)期值相等時(shí),密碼子不存在偏性。所以當(dāng)RSCU=1.00,表明密碼子使用無偏性;RSCU>1.00時(shí),表明該密碼子實(shí)際使用頻率較高;RSCU<1.00時(shí),表明實(shí)際使用頻率較低[23]。

    1.2.6 最優(yōu)密碼子分析。

    選出2種羌活葉綠體編碼基因密碼子ENC值最大的10%和最小的10%,分別建立高低表達(dá)庫(kù),選出庫(kù)中ΔRSCU>0.08且整體RSCU>1.00的密碼子[24]。

    2 結(jié)果與分析

    2.1 密碼子組成

    2種羌活葉綠體基因組密碼子第3位的堿基組成結(jié)果顯示,羌活的葉綠體基因組密碼子第3位的CA3、CT3分別為42.75%與46.97%(表1),明顯高于CC3、CG3,說明羌活葉綠體基因組密碼子多以AT結(jié)尾;寬葉羌活的密碼子第3位堿基組成與羌活類似,AT含量較高。 羌活與寬葉羌活的葉綠體密碼子第1位GC含量從第1位到第3位逐漸下降,第1位分別為46.16%、45.99%,第2位38.17%、37.91%,第3位29.97%、29.67%。2種羌活的ENC均大于49.00,說明2種羌活葉綠體基因組密碼子偏性較弱。

    使用CUSP統(tǒng)計(jì)2種羌活葉綠體基因組編碼基因的CGC、CGC1、CGC2、CGC3、ENC,結(jié)果表明,羌活和寬葉羌活的葉綠體基因組編碼基因的CGC平均值為38.64%和38.44%(表2),差異不大。羌活和寬葉羌活中CGC最高的基因皆為rps11,CGC最低的基因?yàn)閥cf1。羌活和寬葉羌活編碼基因的ENC值范圍分別是36.33~55.03、37.97~57.09,而且羌活約有83%編碼基因的ENC值大于45.00,寬葉羌活有79%編碼基因的ENC值大于45.00,說明多數(shù)編碼基因密碼子使用偏性較弱。2種羌活的rps8的ENC值最小,分別為36.33和37.97,表明rps8在2種羌活葉綠體編碼基因中偏性最強(qiáng)。

    對(duì)2種羌活的葉綠體編碼基因的CGC1、CGC2、CGC3、GC、ENC、葉綠體編碼基因序列長(zhǎng)度(sequence length,SL)等參數(shù)進(jìn)行相關(guān)性分析,結(jié)果表明2種羌活的堿基組成相關(guān)性較接近,CGC1與CGC2皆為極顯著相關(guān)(表3),且CGC1、CGC2都與CGC3無顯著相關(guān)性;2種羌活的GC含量與第1、2、3位組成均極顯著相關(guān)。羌活與寬葉羌活的SL與CGC1、CGC2、CGC3、GC、ENC均無顯著相關(guān)性。2種羌活的相關(guān)性分析結(jié)果也存在一定差異,羌活的ENC與CGC3極顯著相關(guān),而寬葉羌活的ENC除了與CGC3極顯著相關(guān)外,還與GC顯著相關(guān)。

    2.2 中性繪圖分析

    羌活的CGC12取值范圍為 0.324~0.550(圖1),CGC3與CGC12的回歸系數(shù)為0.218 7,皮爾森系數(shù)為0.147,無顯著相關(guān)性,說明羌活的葉綠體密碼子第3位堿基組成與第1、2位差異較大,密碼子偏性主要受選擇的影響。寬葉羌活的CGC12的取值范圍0.321~0.547(圖1),CGC3與CGC12的回歸系數(shù)為0.431 7,皮爾森系數(shù)為0.316,相關(guān)性較低,說明寬葉羌活的葉綠體密碼子第3位堿基組成與第1、2位差異較大,密碼子偏性受選擇影響較大。

    2.3 PR2-plot分析

    PR2-plot分析以CG3/(CG3+CC3)為橫坐標(biāo),CA3/(CA3+CT3)為縱坐標(biāo)繪制散點(diǎn)圖,結(jié)果表明,2種羌活大部分基因分布在右下方位置(圖2),說明羌活和寬葉羌活的葉綠體基因組密碼子不僅受到突變的影響,還受到選擇的影響,并且密碼子第3位的堿基頻率T>A,G>C。

    2.4 ENC-plot分析

    對(duì)2種羌活進(jìn)行ENC-plot繪圖分析,結(jié)果見圖3。標(biāo)準(zhǔn)曲線表示在無選擇壓力的情況下,ENC與CGC3的關(guān)系。如果點(diǎn)主要分布于曲線附近,說明密碼子偏

    性更多受突變的影響,如果點(diǎn)分布于距標(biāo)準(zhǔn)曲線較遠(yuǎn)的位置,說明密碼子偏性更多受選擇的影響。2種羌活的編碼基因在標(biāo)準(zhǔn)曲線兩側(cè)均有基因分布,大部分基因位于標(biāo)準(zhǔn)曲線下側(cè),部分基因距標(biāo)準(zhǔn)曲線較遠(yuǎn)。

    為了進(jìn)一步分析ENC實(shí)際值與標(biāo)準(zhǔn)曲線的偏離情況,制作ENC分布頻數(shù)表,結(jié)果見表4,羌活有24個(gè)基因的比值在-0.05~<0.05,表明有24個(gè)基因的ENC實(shí)際值與期望值接近(ENC比值-0.05~<0.05),其余27個(gè)基因的ENC實(shí)際值與期望值相差較大(ENC比值>0.05或<-0.05),其中psbA、rps8、rp116、petD、psbC、psbD、petB的ENC實(shí)際值距標(biāo)準(zhǔn)曲線較遠(yuǎn)(ENC比值>0.15)。寬葉羌活有23個(gè)基因的比值在-0.05~<0.05,ENC實(shí)際值與期望值接近,27個(gè)基因ENC實(shí)際值與期望值差距較大,其中rps8、petB、petD、psbC、ndhE、ndhI的ENC實(shí)際值距標(biāo)準(zhǔn)曲線較遠(yuǎn)。說明2種羌活的葉綠體基因組密碼子偏性受突變的影響較小,更多受選擇的影響。

    2.5 RSCU分析

    由表5可知,羌活和寬葉羌活各有30個(gè)密碼子RSCU>1.00,說明各有30個(gè)密碼子的使用有偏好性。實(shí)際使用頻率偏高(RSCU>1.00)的密碼子中結(jié)尾為U的有16個(gè),結(jié)尾為A的有13個(gè),只有UUG以G結(jié)尾。說明2種羌活葉綠體基因組密碼子偏好以A、U結(jié)尾,其中U的偏好性略強(qiáng)于A。

    2.6 最優(yōu)密碼子確定

    2種羌活的最優(yōu)密碼子均符合總RSCU>1.00,ΔRSCU≥0.08。2種羌活的最優(yōu)密碼子都多以A/U結(jié)尾,只有UUG以G結(jié)尾。確認(rèn)了羌活的最優(yōu)密碼子9個(gè)(UUG、CUU、ACU、UAU、CAU、AAU、CGU、AGU、GGA),寬葉羌活的最優(yōu)密碼子7個(gè)(UUG、CCA、CAU、CGU、AGU、AGA、GGA)。2種羌活有5個(gè)最優(yōu)密碼子相同,分別是UUG、CAU、CGU、AGU、GGA(表6)。

    3 討論

    該研究對(duì)羌活和寬葉羌活基因組密碼子使用偏好性及其主要影響因素進(jìn)行了探討和比較。結(jié)果表明,ENC小于45的基因較少,多數(shù)基因ENC較高,說明2種羌活基因組密碼子使用偏性均較弱。相關(guān)性分析結(jié)果表明,2種羌活的葉綠體基因組ENC都與CGC3顯著相關(guān),說明密碼子第3位堿基組成對(duì)密碼子偏性的影響較大,這與原曉龍等[25]對(duì)云南藍(lán)果樹的研究結(jié)論一致。而尹為治等[26]研究結(jié)果表明,坡壘的葉綠體基因組密碼子第1位堿基組成對(duì)密碼子偏性的影響更大。密碼子組成分析顯示,2種羌活的總GC和CGC3較低,且CGC1>CGC2>CGC3,密碼子的第3位傾向于使用A/T,第3位堿基使用頻率T>A,G>C。

    突變和選擇壓力是2個(gè)影響密碼子偏性的主要因素。PR2-plot分析結(jié)果發(fā)現(xiàn),基因分布在圖中發(fā)生了偏移,可知2種羌活的葉綠體密碼子使用偏性不僅受到了突變的影響,還受到了選擇壓力的影響。中性繪圖結(jié)果顯示,2種羌活CGC12與CGC3相關(guān)性不顯著,說明密碼子第3位的堿基組成與前兩位差異大,密碼子相對(duì)保守,表明選擇壓力對(duì)密碼子的影響較大。ENC-plot結(jié)果顯示,2種羌活的多數(shù)葉綠體基因偏離標(biāo)準(zhǔn)曲線進(jìn)一步證實(shí)選擇壓力是影響2種羌活葉綠體密碼子使用偏性的主要因素。綜上所述,2種羌活的葉綠體基因組密碼子使用偏性受到突變和選擇壓力的影響,最主要的影響因素是選擇壓力。影響竹亞科的巴山木竹、印度簕竹、Olyra latifolia L.以及苦蕎的葉綠體基因組密碼子使用偏性的主要因素也是選擇壓力[27-28]。而突變對(duì)槲蕨屬的川滇槲蕨、櫟葉槲蕨的葉綠體基因組密碼子使用偏性的影響更大[29]。

    羌活和寬葉羌活的葉綠體基因組密碼子使用模式較為相似。2種羌活的總GC含量與密碼子各位的GC含量差異不大;葉綠體基因組密碼子偏性都較弱;密碼子第3位都傾向于使用A/T,且堿基含量T>A,G>C;影響密碼子使用偏性的主要因素都是選擇壓力。

    通過RSCU分析,確定了2種羌活中使用頻率較高的30個(gè)高頻密碼子(RSCU>1.00),幾乎都以A或U結(jié)尾,只有UUG以G結(jié)尾。確認(rèn)了羌活的最優(yōu)密碼子9個(gè)(UUG、CUU、ACU、UAU、CAU、AAU、CGU、AGU、GGA),寬葉羌活的最優(yōu)密碼子7個(gè)(UUG、CCA、CAU、CGU、AGU、AGA、GGA)。2種羌活有5個(gè)最優(yōu)密碼子相同,分別是UUG、CAU、CGU、AGU、GGA。

    該研究通過對(duì)2種羌活藥用植物的葉綠體基因組編碼基因的使用偏性進(jìn)行比較分析,得出主要結(jié)論:①2種藥用羌活葉綠體基因組的密碼子使用偏性整體較為相似,且偏性較弱;②密碼子的第3位傾向于使用A、T;③突變和選擇壓力在2種羌活的密碼子使用偏性中起主要作用,其中選擇壓力對(duì)密碼子使用偏性影響更大。該研究結(jié)果為后續(xù)藥用羌活植物的基因表達(dá)、蛋白質(zhì)翻譯和遺傳育種研究提供了研究依據(jù)和基礎(chǔ)數(shù)據(jù),有潛在的經(jīng)濟(jì)價(jià)值和應(yīng)用價(jià)值。

    參考文獻(xiàn)

    [1] MARAIS G,MOUCHIROUD D,DURET L.Neutral effect of recombination on base composition in Drosophila[J].Genetical research,2003,81(2):79-87.

    [2] PRABHA R,SINGH D P,SINHA S,et al.Genome-wide comparative analysis of codon usage bias and codon context patterns among cyanobacterial genomes[J].Marine genomics,2017,32:31-39.

    [3] VICARIO S,MORIYAMA E N,POWELL J R.Codon usage in twelve species of Drosophila[J].BMC evolutionary biology,2007,7:1-17.

    [4] SHARP P M,LI W H.The rate of synonymous substitution in enterobacterial genes is inversely related to codon usage bias[J].Molecular biology and evolution,1987,4(3):222-230.

    [5] OLEJNICZAK M,UHLENBECK O C.tRNA residues that have coevolved with their anticodon to ensure uniform and accurate codon recognition[J].Biochimie,2006,88(8):943-950.

    [6] ZALUCKI Y M,POWER P M,JENNINGS M P.Selection for efficient translation initiation biases codon usage at second amino acid position in secretory proteins[J].Nucleic acids research,2007,35(17):5748-5754.

    [7] HUANG Y,KOONIN E V,LIPMAN D J,et al.Selection for minimization of translational frameshifting errors as a factor in the evolution of codon usage[J].Nucleic acids research,2009,37(20):6799-6810.

    [8] 孫錚,馬亮,MURPHY R,等.Wolbachia pipientis wMel 基因組水平上的密碼子使用分析[J].中國(guó)科學(xué)(C輯:生命科學(xué)),2009,39(10):948-953.

    [9] HIRAOKA Y,KAWAMATA K,HARAGUCHI T,et al.Codon usage bias is correlated with gene expression levels in the fission yeast Schizosaccharomyces pombe[J].Genes to cells,2009,14(4):499-509.

    [10] SUBRAMANIAN S.Nearly neutrality and the evolution of codon usage bias in eukaryotic genomes[J].Genetics,2008,178(4):2429-2432.

    [11] 黃瓊林.高良姜葉綠體基因組測(cè)序與特征分析[J].熱帶作物學(xué)報(bào),2021,42(1):1-6.

    [12] 原曉龍,李云琴,張勁峰,等.乳油木葉綠體基因組密碼子偏好性分析[J].分子植物育種,2020,18(17):5658-5664.

    [13] 辛雅萱,董章宏,瞿紹宏,等.杜梨葉綠體基因組密碼子偏好性分析[J].河北農(nóng)業(yè)大學(xué)學(xué)報(bào),2020,43(6):51-59.

    [14] 胡曉艷,許艷秋,韓有志,等.酸棗葉綠體基因組密碼子使用偏性分析[J].森林與環(huán)境學(xué)報(bào),2019,39(6):621-628.

    [15] 劉文武,蔣曉文,張帥,等.羌活中香豆素類化學(xué)成分及其抗氧化活性研究[J].中草藥,2019,50(6):1310-1315.

    [16] 楊嬌.羌活屬葉綠體基因組進(jìn)化及物種界定研究[D].西安:西北大學(xué),2018:5-6.

    [17] 閻夢(mèng)穎.珍稀瀕危羌活屬藥用植物的譜系地理學(xué)研究[D].西安:西北大學(xué),2016:5-6.

    [18] 郭培,郎擁軍,張國(guó)桃.羌活化學(xué)成分及藥理活性研究進(jìn)展[J].中成藥,2019,41(10):2445-2459.

    [19] 趙森,鄧力華,陳芬.秋茄葉綠體基因組密碼子使用偏好性分析[J].森林與環(huán)境學(xué)報(bào),2020,40(5):534-541.

    [20] 吳憲明,吳松鋒,任大明,等.密碼子偏性的分析方法及相關(guān)研究進(jìn)展[J].遺傳,2007,29(4):420-426.

    [21] 尚明照,劉方,華金平,等.陸地棉葉綠體基因組密碼子使用偏性的分析[J].中國(guó)農(nóng)業(yè)科學(xué),2011,44(2):245-253.

    [22] WRIGHT F.The‘effective number of codons’used in a gene[J].Gene,1990,87(1):23-29.

    [23] SHARP P M,LI W H.The codon Adaptation Index—A measure of directional synonymous codon usage bias,and its potential applications[J].Nucleic Acids Res,1987,15(3):1281-1295.

    [24] 陸奇豐,駱文華,黃至歡.兩種梧桐葉綠體基因組密碼子使用偏性分析[J].廣西植物,2020,40(2):173-183.

    [25] 原曉龍,康洪梅,王毅.云南藍(lán)果樹葉綠體基因組密碼子偏好性分析[J].西北林學(xué)院學(xué)報(bào),2020,35(4):26-31,124.

    [26] 尹為治,方正,黃良鴻,等.海南2種龍腦香科植物葉綠體基因組密碼子偏好性分析[J].林業(yè)調(diào)查規(guī)劃,2020,45(6):19-26.

    [27] 吳妙麗,陳世品,陳輝.竹亞科葉綠體基因組的密碼子使用偏性分析[J].森林與環(huán)境學(xué)報(bào),2019,39(1):9-14.

    [28] 胡莎莎,羅洪,吳琦,等.苦蕎葉綠體基因組密碼子偏愛性分析[J].分子植物育種,2016,14(2):309-317.

    [29] 沈宗芳,陸添權(quán),張志榮,等.槲蕨屬葉綠體基因組密碼子偏好性分析[J].廣西植物,2021,41(2):266-273.

    猜你喜歡
    羌活
    羌活本草考證及炮制歷史沿革分析
    觀察羌活勝濕湯加減治療腰痛臨床療效
    羌活治風(fēng)濕的故事
    保健與生活(2021年8期)2021-04-19 00:08:04
    云南寬葉羌活揮發(fā)油成分的GC-MS分析及應(yīng)用前景探討
    不同基源羌活含量測(cè)定和指紋圖譜結(jié)合化學(xué)模式識(shí)別研究
    羌活的臨床應(yīng)用及其用量探究
    羌活
    柱層析用硅膠對(duì)羌活中主要成分的影響
    中成藥(2017年10期)2017-11-16 00:50:15
    羌活屬植物在不同海拔引種的生態(tài)適宜性研究
    野生藥材羌活馴化高效育苗探析
    種子科技(2017年4期)2017-01-13 07:03:37
    欧美黑人精品巨大| 天天躁夜夜躁狠狠久久av| 亚洲精品日韩在线中文字幕| 国产精品偷伦视频观看了| 国产男人的电影天堂91| 国产成人免费无遮挡视频| 久久久久久人妻| 国产亚洲午夜精品一区二区久久| 成年女人毛片免费观看观看9 | 99热国产这里只有精品6| 国产精品免费大片| 午夜福利影视在线免费观看| 熟妇人妻不卡中文字幕| 一本一本久久a久久精品综合妖精| 欧美精品一区二区大全| 爱豆传媒免费全集在线观看| 亚洲欧美一区二区三区久久| 中文字幕亚洲精品专区| 国产又色又爽无遮挡免| 久久99一区二区三区| 色婷婷av一区二区三区视频| 伊人亚洲综合成人网| 成人18禁高潮啪啪吃奶动态图| 男女下面插进去视频免费观看| 久久久欧美国产精品| 久久久久久免费高清国产稀缺| 丝袜脚勾引网站| 一本久久精品| 婷婷色av中文字幕| 赤兔流量卡办理| 免费观看av网站的网址| 极品人妻少妇av视频| 9色porny在线观看| 少妇人妻精品综合一区二区| 日韩电影二区| 国产免费现黄频在线看| 99精品久久久久人妻精品| 成人影院久久| 亚洲精品国产av成人精品| 久久久久久免费高清国产稀缺| 亚洲第一青青草原| 男女边摸边吃奶| 十分钟在线观看高清视频www| av卡一久久| 国产欧美日韩一区二区三区在线| 在线精品无人区一区二区三| 青草久久国产| 亚洲av电影在线观看一区二区三区| 亚洲欧美一区二区三区国产| 青春草视频在线免费观看| 欧美xxⅹ黑人| 精品一区二区免费观看| 国产99久久九九免费精品| 国产精品久久久久久久久免| 亚洲在久久综合| 亚洲美女黄色视频免费看| 制服人妻中文乱码| 一级黄片播放器| 少妇被粗大的猛进出69影院| 国产精品 国内视频| 日本猛色少妇xxxxx猛交久久| 老司机在亚洲福利影院| 国产又爽黄色视频| 在线天堂最新版资源| 国产成人精品久久二区二区91 | 男女高潮啪啪啪动态图| 咕卡用的链子| 新久久久久国产一级毛片| 午夜福利在线免费观看网站| 久久这里只有精品19| 国产1区2区3区精品| 飞空精品影院首页| 大香蕉久久网| 成人国产麻豆网| 亚洲成人手机| 亚洲av成人不卡在线观看播放网 | 久久久精品免费免费高清| 日韩一本色道免费dvd| 伦理电影免费视频| 久久久精品区二区三区| 国产在线视频一区二区| 美女主播在线视频| 一二三四中文在线观看免费高清| 午夜精品国产一区二区电影| 国产精品.久久久| 日韩精品有码人妻一区| 啦啦啦中文免费视频观看日本| 晚上一个人看的免费电影| 国产一区二区三区综合在线观看| 免费观看a级毛片全部| 国产成人av激情在线播放| 另类精品久久| 十八禁网站网址无遮挡| 欧美日韩一级在线毛片| 国产免费视频播放在线视频| 国产亚洲精品第一综合不卡| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕色久视频| 777久久人妻少妇嫩草av网站| 国产成人午夜福利电影在线观看| 丝袜人妻中文字幕| 国产一区二区 视频在线| 另类精品久久| 久久久国产精品麻豆| 69精品国产乱码久久久| 国产精品国产三级专区第一集| 国产不卡av网站在线观看| 欧美精品av麻豆av| 丝袜美腿诱惑在线| 99re6热这里在线精品视频| 亚洲精品视频女| 国产人伦9x9x在线观看| 啦啦啦啦在线视频资源| 啦啦啦 在线观看视频| 久久精品久久久久久噜噜老黄| 大片免费播放器 马上看| 国产亚洲av片在线观看秒播厂| 狠狠婷婷综合久久久久久88av| 欧美日韩亚洲国产一区二区在线观看 | 久久av网站| 纵有疾风起免费观看全集完整版| 激情五月婷婷亚洲| 欧美日韩亚洲高清精品| 成人免费观看视频高清| 少妇人妻久久综合中文| 老司机深夜福利视频在线观看 | 精品亚洲乱码少妇综合久久| 丝袜美腿诱惑在线| 精品一区二区免费观看| 国产国语露脸激情在线看| 日日啪夜夜爽| 黑人巨大精品欧美一区二区蜜桃| 飞空精品影院首页| 纯流量卡能插随身wifi吗| 久久久久精品人妻al黑| 精品国产一区二区久久| 国产亚洲最大av| 亚洲色图综合在线观看| 久热这里只有精品99| 如何舔出高潮| 亚洲欧美一区二区三区久久| 亚洲综合精品二区| 天堂俺去俺来也www色官网| 精品亚洲成国产av| 日韩av在线免费看完整版不卡| 精品一区在线观看国产| 999精品在线视频| 90打野战视频偷拍视频| 国产精品人妻久久久影院| 秋霞在线观看毛片| 日本黄色日本黄色录像| 中文字幕色久视频| 天堂俺去俺来也www色官网| 国产精品.久久久| av又黄又爽大尺度在线免费看| 美女国产高潮福利片在线看| 高清视频免费观看一区二区| bbb黄色大片| 老司机深夜福利视频在线观看 | 超碰97精品在线观看| 久久久国产精品麻豆| 亚洲综合色网址| 亚洲精品第二区| 捣出白浆h1v1| 成人国语在线视频| 国产精品国产三级国产专区5o| 看免费成人av毛片| 亚洲成国产人片在线观看| 欧美精品亚洲一区二区| 十八禁网站网址无遮挡| 少妇的丰满在线观看| 最近2019中文字幕mv第一页| 亚洲精品美女久久av网站| 一本—道久久a久久精品蜜桃钙片| 欧美日韩精品网址| 不卡视频在线观看欧美| 最新在线观看一区二区三区 | 日韩制服丝袜自拍偷拍| 国产成人精品久久二区二区91 | 日韩一区二区视频免费看| 青草久久国产| 婷婷色综合www| 熟女av电影| 人人妻人人澡人人看| 曰老女人黄片| 国产在线免费精品| 国产男人的电影天堂91| 飞空精品影院首页| 亚洲一码二码三码区别大吗| 制服人妻中文乱码| 91成人精品电影| 国产精品国产三级专区第一集| 99久久综合免费| 男女无遮挡免费网站观看| 欧美精品一区二区免费开放| 一本久久精品| 国产成人一区二区在线| 18禁裸乳无遮挡动漫免费视频| 亚洲熟女毛片儿| 最近中文字幕2019免费版| 一级片'在线观看视频| 另类精品久久| 高清黄色对白视频在线免费看| 黑人欧美特级aaaaaa片| 在线观看人妻少妇| 国产av精品麻豆| 成人毛片60女人毛片免费| 一二三四在线观看免费中文在| 欧美成人精品欧美一级黄| 不卡视频在线观看欧美| 街头女战士在线观看网站| 亚洲 欧美一区二区三区| 国产精品蜜桃在线观看| 1024香蕉在线观看| 丝袜美腿诱惑在线| 一级毛片 在线播放| 国产精品嫩草影院av在线观看| 日韩电影二区| 日韩,欧美,国产一区二区三区| 国产精品 国内视频| 久久精品久久久久久噜噜老黄| 男男h啪啪无遮挡| 婷婷色综合www| 精品人妻在线不人妻| 日韩,欧美,国产一区二区三区| 人体艺术视频欧美日本| 色精品久久人妻99蜜桃| 七月丁香在线播放| 欧美乱码精品一区二区三区| 婷婷成人精品国产| 晚上一个人看的免费电影| 午夜福利免费观看在线| 亚洲国产欧美在线一区| 麻豆乱淫一区二区| 99re6热这里在线精品视频| 午夜av观看不卡| 午夜91福利影院| 欧美在线黄色| 国产精品成人在线| 午夜日韩欧美国产| 精品人妻一区二区三区麻豆| 制服丝袜香蕉在线| 最近最新中文字幕免费大全7| 一边摸一边抽搐一进一出视频| 亚洲精品一二三| 国产精品一国产av| 最近中文字幕高清免费大全6| 在现免费观看毛片| 久热爱精品视频在线9| 高清在线视频一区二区三区| 亚洲精华国产精华液的使用体验| 日韩av不卡免费在线播放| 国产高清国产精品国产三级| 最近的中文字幕免费完整| 久久精品aⅴ一区二区三区四区| 香蕉国产在线看| 欧美乱码精品一区二区三区| 少妇精品久久久久久久| 精品国产一区二区三区久久久樱花| 免费观看a级毛片全部| 高清不卡的av网站| 日本爱情动作片www.在线观看| 国产乱人偷精品视频| 国产精品 国内视频| 黄色怎么调成土黄色| 侵犯人妻中文字幕一二三四区| 日本欧美视频一区| 日本黄色日本黄色录像| 国产成人精品久久二区二区91 | 国产精品欧美亚洲77777| 黄色毛片三级朝国网站| 激情视频va一区二区三区| 99re6热这里在线精品视频| 亚洲天堂av无毛| 国产精品偷伦视频观看了| 亚洲人成网站在线观看播放| 色94色欧美一区二区| 日韩av在线免费看完整版不卡| 2018国产大陆天天弄谢| 最近2019中文字幕mv第一页| 人妻 亚洲 视频| 国产精品一国产av| 王馨瑶露胸无遮挡在线观看| 99re6热这里在线精品视频| 男女边吃奶边做爰视频| 久久国产精品男人的天堂亚洲| 免费在线观看黄色视频的| 一边摸一边做爽爽视频免费| 欧美xxⅹ黑人| 女人爽到高潮嗷嗷叫在线视频| 国产一区二区三区综合在线观看| 青春草亚洲视频在线观看| 在线观看人妻少妇| 日本wwww免费看| 激情视频va一区二区三区| 一级片免费观看大全| 欧美乱码精品一区二区三区| 亚洲av日韩精品久久久久久密 | 男女高潮啪啪啪动态图| 美女中出高潮动态图| 成人黄色视频免费在线看| 在线 av 中文字幕| 欧美最新免费一区二区三区| 国产1区2区3区精品| 亚洲,欧美,日韩| 国产亚洲欧美精品永久| 精品亚洲成a人片在线观看| 哪个播放器可以免费观看大片| 国产精品免费视频内射| 午夜福利网站1000一区二区三区| 美国免费a级毛片| 亚洲国产日韩一区二区| 亚洲国产精品一区三区| 亚洲视频免费观看视频| 日韩av免费高清视频| 韩国精品一区二区三区| 免费在线观看黄色视频的| 婷婷色av中文字幕| 国产精品偷伦视频观看了| 久久久久国产精品人妻一区二区| 国产成人精品久久二区二区91 | 久久精品人人爽人人爽视色| 亚洲精品av麻豆狂野| 国产精品二区激情视频| 欧美精品av麻豆av| 国产精品久久久人人做人人爽| 国产精品人妻久久久影院| 新久久久久国产一级毛片| 视频区图区小说| 国产乱来视频区| 美女大奶头黄色视频| 老司机深夜福利视频在线观看 | 国产一区二区激情短视频 | 国产人伦9x9x在线观看| 久久毛片免费看一区二区三区| 啦啦啦在线观看免费高清www| 亚洲国产精品国产精品| 亚洲伊人久久精品综合| 十八禁高潮呻吟视频| 久久久精品国产亚洲av高清涩受| 自线自在国产av| 精品国产乱码久久久久久小说| 美女视频免费永久观看网站| 哪个播放器可以免费观看大片| 免费看av在线观看网站| 亚洲色图 男人天堂 中文字幕| 青春草国产在线视频| 电影成人av| 亚洲国产精品成人久久小说| 亚洲七黄色美女视频| 国语对白做爰xxxⅹ性视频网站| 伦理电影大哥的女人| 国产av国产精品国产| 亚洲国产日韩一区二区| 久久 成人 亚洲| 国产探花极品一区二区| 综合色丁香网| 国产精品二区激情视频| 久久99热这里只频精品6学生| 人人妻人人爽人人添夜夜欢视频| 亚洲欧洲精品一区二区精品久久久 | 欧美最新免费一区二区三区| 男男h啪啪无遮挡| 亚洲av日韩精品久久久久久密 | 午夜激情av网站| av网站免费在线观看视频| 色婷婷久久久亚洲欧美| 最近最新中文字幕大全免费视频 | 久久国产精品男人的天堂亚洲| 少妇的丰满在线观看| 国产精品麻豆人妻色哟哟久久| 亚洲av男天堂| 亚洲欧美色中文字幕在线| 日韩中文字幕视频在线看片| 日韩 欧美 亚洲 中文字幕| 久久久国产一区二区| 国产精品香港三级国产av潘金莲 | 一二三四在线观看免费中文在| 日韩人妻精品一区2区三区| 亚洲国产精品成人久久小说| 色精品久久人妻99蜜桃| 亚洲熟女毛片儿| 精品福利永久在线观看| 国产精品亚洲av一区麻豆 | 一边亲一边摸免费视频| 国产成人免费无遮挡视频| 免费黄网站久久成人精品| 夫妻性生交免费视频一级片| 一边亲一边摸免费视频| 国产乱来视频区| 99国产综合亚洲精品| 可以免费在线观看a视频的电影网站 | 日本91视频免费播放| 一级片免费观看大全| 国精品久久久久久国模美| 日日啪夜夜爽| 亚洲国产精品一区二区三区在线| 国产精品无大码| 啦啦啦中文免费视频观看日本| 黄色视频不卡| 久久精品国产综合久久久| 黑丝袜美女国产一区| 国产精品女同一区二区软件| 十八禁高潮呻吟视频| 成人免费观看视频高清| 亚洲成国产人片在线观看| 丝袜脚勾引网站| 亚洲人成77777在线视频| 亚洲精品乱久久久久久| 亚洲四区av| 午夜福利一区二区在线看| 伊人亚洲综合成人网| 中文字幕色久视频| 日韩视频在线欧美| 欧美国产精品一级二级三级| 另类亚洲欧美激情| 两个人看的免费小视频| 国产在线视频一区二区| 精品卡一卡二卡四卡免费| 色网站视频免费| 老司机亚洲免费影院| 日韩av不卡免费在线播放| 青春草国产在线视频| 欧美人与性动交α欧美精品济南到| 黄色 视频免费看| 国产精品 国内视频| 免费在线观看完整版高清| 亚洲av日韩精品久久久久久密 | 久久鲁丝午夜福利片| 亚洲欧美一区二区三区国产| 日韩av不卡免费在线播放| 天堂俺去俺来也www色官网| 69精品国产乱码久久久| 精品亚洲成国产av| 免费av中文字幕在线| 国产毛片在线视频| 午夜日韩欧美国产| 永久免费av网站大全| 一二三四在线观看免费中文在| 国产精品偷伦视频观看了| 久久久久久久大尺度免费视频| 亚洲av成人精品一二三区| 女人高潮潮喷娇喘18禁视频| 日韩av免费高清视频| 久久亚洲国产成人精品v| 校园人妻丝袜中文字幕| 叶爱在线成人免费视频播放| 亚洲国产精品成人久久小说| 91aial.com中文字幕在线观看| 女的被弄到高潮叫床怎么办| 亚洲一级一片aⅴ在线观看| 亚洲精品久久成人aⅴ小说| 丰满迷人的少妇在线观看| 美女高潮到喷水免费观看| 日韩免费高清中文字幕av| 久久久久久久久久久久大奶| 丝袜在线中文字幕| 日韩大片免费观看网站| 亚洲视频免费观看视频| 久久国产亚洲av麻豆专区| 免费女性裸体啪啪无遮挡网站| 国产黄色视频一区二区在线观看| 美女主播在线视频| 一级爰片在线观看| av天堂久久9| 欧美日本中文国产一区发布| 中文字幕精品免费在线观看视频| 又大又爽又粗| 午夜精品国产一区二区电影| 在线观看三级黄色| 高清视频免费观看一区二区| 丰满少妇做爰视频| 日韩av在线免费看完整版不卡| 飞空精品影院首页| 少妇的丰满在线观看| 免费看不卡的av| 久久99热这里只频精品6学生| 人体艺术视频欧美日本| 侵犯人妻中文字幕一二三四区| 别揉我奶头~嗯~啊~动态视频 | 久久热在线av| 久久av网站| 黑人欧美特级aaaaaa片| 大香蕉久久网| a级毛片黄视频| 黄片无遮挡物在线观看| 国产亚洲av片在线观看秒播厂| 久久久精品94久久精品| 久久鲁丝午夜福利片| 我要看黄色一级片免费的| 十八禁人妻一区二区| 伊人久久大香线蕉亚洲五| 成年女人毛片免费观看观看9 | 精品国产一区二区三区久久久樱花| 一本大道久久a久久精品| 欧美精品人与动牲交sv欧美| 美女视频免费永久观看网站| 日韩中文字幕视频在线看片| 欧美黑人欧美精品刺激| 婷婷色综合www| 国产极品天堂在线| 十八禁高潮呻吟视频| 免费av中文字幕在线| 精品酒店卫生间| 亚洲精品日本国产第一区| 久久99精品国语久久久| videos熟女内射| 欧美少妇被猛烈插入视频| 久久久久精品人妻al黑| 热99国产精品久久久久久7| 91国产中文字幕| 99久久人妻综合| 欧美国产精品一级二级三级| 看十八女毛片水多多多| 九色亚洲精品在线播放| 日韩不卡一区二区三区视频在线| 人妻 亚洲 视频| 在线观看免费高清a一片| 亚洲七黄色美女视频| 亚洲av综合色区一区| 精品一区二区免费观看| 国产不卡av网站在线观看| netflix在线观看网站| 亚洲av欧美aⅴ国产| 美女视频免费永久观看网站| 在线观看三级黄色| 亚洲国产日韩一区二区| 曰老女人黄片| 成人国产麻豆网| 久久久精品免费免费高清| 午夜福利视频精品| 亚洲精品视频女| 精品亚洲乱码少妇综合久久| 国产精品一区二区在线不卡| 叶爱在线成人免费视频播放| 性色av一级| 亚洲免费av在线视频| 国产伦理片在线播放av一区| 丰满少妇做爰视频| 久久性视频一级片| 午夜福利视频在线观看免费| 久久久久人妻精品一区果冻| 久久久久久久久久久免费av| av女优亚洲男人天堂| 日韩伦理黄色片| 久久久精品免费免费高清| 色播在线永久视频| 亚洲人成电影观看| 波多野结衣av一区二区av| 欧美日韩视频高清一区二区三区二| 成人亚洲欧美一区二区av| 欧美精品亚洲一区二区| 精品久久久久久电影网| 久久亚洲国产成人精品v| 夜夜骑夜夜射夜夜干| 一区二区三区精品91| 亚洲欧美色中文字幕在线| xxxhd国产人妻xxx| 男的添女的下面高潮视频| 久久人人97超碰香蕉20202| 国产男人的电影天堂91| 伦理电影大哥的女人| 少妇人妻精品综合一区二区| 桃花免费在线播放| 精品一品国产午夜福利视频| 色播在线永久视频| 国产黄频视频在线观看| 新久久久久国产一级毛片| 黄片播放在线免费| 欧美日韩综合久久久久久| 亚洲精品在线美女| 亚洲精品乱久久久久久| 黄色视频不卡| 侵犯人妻中文字幕一二三四区| 亚洲av男天堂| 精品少妇黑人巨大在线播放| 亚洲四区av| 日本一区二区免费在线视频| 国产精品国产av在线观看| 一本久久精品| 观看av在线不卡| av视频免费观看在线观看| av电影中文网址| 哪个播放器可以免费观看大片| h视频一区二区三区| 9191精品国产免费久久| 只有这里有精品99| 精品一区二区三卡| 国精品久久久久久国模美| 久久天躁狠狠躁夜夜2o2o | 成人亚洲精品一区在线观看| 97精品久久久久久久久久精品| 最近最新中文字幕免费大全7| 日本黄色日本黄色录像| 狂野欧美激情性bbbbbb| 成人国产av品久久久| 成人18禁高潮啪啪吃奶动态图| 免费女性裸体啪啪无遮挡网站| 久久久久久人人人人人| 午夜久久久在线观看| 免费高清在线观看日韩| 婷婷色av中文字幕| svipshipincom国产片| 亚洲成人av在线免费| 黄色 视频免费看| 丰满乱子伦码专区| 国产精品av久久久久免费| 欧美成人精品欧美一级黄| 国产亚洲最大av| 九色亚洲精品在线播放| 一本—道久久a久久精品蜜桃钙片| 黄色视频不卡| 亚洲熟女毛片儿| 美女扒开内裤让男人捅视频| 美女脱内裤让男人舔精品视频|