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

    鼠尾草葉綠體基因組序列分析

    2017-11-11 06:24:20何懿菡韓立敏劉玉萍王喆之
    植物研究 2017年4期
    關鍵詞:鼠尾草內含子葉綠體

    何懿菡 韓立敏 劉玉萍 田 娜 蘇 旭 王喆之*

    (1.陜西師范大學生命科學學院,藥用資源與天然藥物化學教育部重點實驗室,西北瀕危藥材資源開發(fā)國家工程實驗室,西安 710119; 2.青海師范大學生命與地理科學學院,西寧 810008; 3.陜西學前師范學院生命科學與食品工程學院,西安 710061)

    鼠尾草葉綠體基因組序列分析

    何懿菡1,2韓立敏3劉玉萍2田 娜1蘇 旭2王喆之1*

    (1.陜西師范大學生命科學學院,藥用資源與天然藥物化學教育部重點實驗室,西北瀕危藥材資源開發(fā)國家工程實驗室,西安 710119;2.青海師范大學生命與地理科學學院,西寧 810008;3.陜西學前師范學院生命科學與食品工程學院,西安 710061)

    鼠尾草(Salviajaponica)是唇形科(Labiatae)鼠尾草屬(Salvia)的一種多年生草本植物,具有十分重要的藥用和經(jīng)濟價值。本文采用第二代測序技術Illumina Hiseq平臺對鼠尾草的葉綠體基因組進行測序,同時以鼠尾草近緣物種丹參葉綠體基因組作為參考,組裝得到完整葉綠體基因組序列。結果表明,鼠尾草葉綠體基因組序列全長153 995 bp,呈典型的四段式結構,其中LSC區(qū)長84 573 bp,SSC區(qū)長19 874 bp,兩個IR區(qū)分別長24 774 bp;鼠尾草葉綠體基因組成功注釋13組葉綠體基因,基因的種類、數(shù)目及GC含量等與唇形科中其它物種較為類似。這些研究結果豐富了鼠尾草屬的葉綠體基因組數(shù)據(jù),為今后鼠尾草屬植物系統(tǒng)發(fā)育關系重建積累了基礎性數(shù)據(jù)。

    鼠尾草屬;鼠尾草;葉綠體基因組;基因結構;系統(tǒng)進化

    鼠尾草(Salviajaponica)是唇形科(Labiatae)、鼠尾草屬(Salvia)中的一種多年生草本植物,在我國主要分布于橫斷山區(qū),占鼠尾草屬物種總數(shù)的52.8%,尤其是特有種數(shù)目較多,高達23種[1]。鼠尾草藥材基原為全草,具有清熱利濕、活血調經(jīng)解毒消腫之功效[2],同時含有三萜、類固醇、酚類化合物等多種代謝產物[3],其根和葉含有咖啡酸、迷迭香酸、丹酚酸等8種水溶性成分,此外根中還含有丹參酮Ⅰ和丹參酮ⅡA等脂溶性成分[4]。

    鼠尾草屬植物因雄蕊形成具活動關節(jié)的杠桿狀結構,與昆蟲或蜂鳥的傳粉行為高度適應,因而它們在唇形科有較為重要的系統(tǒng)地位[5],并且Bentham[6]還曾利用這一特征將其劃分為四個亞屬,至今這種分類系統(tǒng)仍被廣泛采用。隨著植物分子生物學測序技術和分析手段的迅速發(fā)展,葉綠體基因的psbA-trnH、matK和rbcL等多個片段被廣泛用于鼠尾草屬植物的DNA分子鑒定[7]。然而,先前鼠尾草屬植物系統(tǒng)發(fā)育的分析均是基于葉綠體基因和核基因序列信息,DNA片段長度通常介于300~800 bp,信息位點較少,片段信息常常存在測定序列的分子進化速率與物種水平上的進化速率存在差異,以檢測序列位點的多樣性來代替物種的多樣性,未能真實反映物種間的系統(tǒng)發(fā)育關系[8];加之,鼠尾草屬種類繁多,種間變異復雜多樣,從而導致對該屬絕大部分物種的界定存在較多不足乃至張冠李戴。當前,隨著高通量測序技術的發(fā)展,葉綠體基因組序列呈現(xiàn)出具有更強分辨率且適用于不同分類階元的系統(tǒng)發(fā)育分析[9~10],可為藥用植物物種鑒定和混偽品鑒別提供較為準確的方法。鼠尾草屬細胞器基因組數(shù)據(jù)匱乏,進而限制了葉綠體基因組在物種鑒定和進化關系研究中的應用。藥用植物丹參(Salviamiltiorrhiza)是鼠尾草屬第一個獲得完整葉綠體基因組的物種,通過與唇形目(Lamiales)已知葉綠體基因組數(shù)據(jù)的綜合分析,發(fā)現(xiàn)它與芝麻(Sesamumindicum)親緣關系最近[11]。鼠尾草和丹參均隸屬于鼠尾草屬,是該屬的一個廣布種,并且具有復雜的遺傳多樣性[12],加之它的物種起源及系統(tǒng)進化等諸多科學問題至今未見報道,尤其是目前鼠尾草葉綠體基因組研究尚屬空白,據(jù)此本文以鼠尾草為研究對象,以丹參葉綠體基因組作為參考,通過對其葉綠體基因組的組裝和分析,探討鼠尾草葉綠體基因組的結構信息,旨在為鼠尾草種質資源的開發(fā)和利用提供科學依據(jù),同時也可為今后鼠尾草屬植物系統(tǒng)進化地位的確立及種質資源的鑒定提供基礎性數(shù)據(jù)。

    1 材料與方法

    1.1 測序材料

    鼠尾草種子撒播于室溫大棚,覆膜待出芽培養(yǎng)3個月后,取幼嫩葉片作為提取DNA的材料。DNA的提取參照Biospin植物基因組DNA提取試劑盒(Bioer,杭州)說明書進行。

    1.2 測序

    將檢測合格后的DNA,隨機打斷成1~2 kB的片段,富集大小一致的片段,連接至建庫載體中構建DNA文庫(DNA文庫制備試劑盒,利用Illumina公司產品),然后利用Illumina HiSeq 2500基因組測序平臺測序。

    1.3 葉綠體基因組的組裝與注釋

    測序獲得的raw reads數(shù)據(jù)利用CLC Genomics Workbench 9.0(CLC Bio,Aarhus,Denmark)軟件中Trim過濾去除低質量序列,過濾條件為Ambiguous limit=2、Quality limit=0.05;然后,以先前提交至NCBI的丹參葉綠體基因組(HF586694)作為參照序列,將Trim過濾后的8,877,028 reads用MITObim 1.7重建葉綠體基因組;最終得到由無gap的717 586條葉綠體序列,平均覆蓋702個堿基,長度為157 839 bp的葉綠體基因組,應用Geneious 8.0.3首先進行丹參與鼠尾草序列的CLUSTAL比對,然后進行堿基修正,修改蛋白質和tRNA等注釋,完成鼠尾草葉綠體基因組的分區(qū),利用OGDraw 1.2(http://ogdraw.mpimp-golm.mpg.de/index.shtml)軟件對鼠尾草葉綠體基因組序列進行繪圖。

    2 結果與分析

    2.1 鼠尾草葉綠體基因組基本特征

    使用Illumina Hiseq平臺對鼠尾草葉綠體DNA提取樣本進行測序,測序數(shù)據(jù)達2.4 G,目前基因組序列已提交至GenBank數(shù)據(jù)庫(登錄號KY646163)。鼠尾草葉綠體基因組序列全長153,995 bp,呈典型的四段式結構(圖1),其中LSC區(qū)(large single copy region,大的單拷貝區(qū))長84 573 bp,SSC區(qū)(small single copy region,小的單拷貝區(qū))長19 874 bp,兩個IR區(qū)(inverted repeat,反向重復區(qū))分別長24 774 bp(表1)??傮w看來A+T含量為62%,有明顯的AT偏向性。IR區(qū)的GC含量(43.4%)明顯高于LSC區(qū)(36.2%)和SSC區(qū)(32.0%),這主要由IR區(qū)包含的四個高GC含量的rRNA基因。

    2.2 鼠尾草葉綠體基因組各分區(qū)基因組成

    整條葉綠體基因組序列共編碼135個基因,去除重復后為118個,包括89(去除重復83)個蛋白編碼基因、37(去除重復30)個tRNA、8(去除重復4)個rRNA(表2)。6個蛋白編碼基因、7個tRNA基因及4個rRNA基因在兩個IR區(qū)各有一個拷貝。LSC區(qū)包含63個蛋白編碼基因、22個tRNA基因,而SSC區(qū)只含有12個蛋白編碼基因和1個tRNA基因(圖1)。

    表1鼠尾草葉綠體基因組堿基組成

    Table1BasecompositionofchloroplastgenomeinS.japonica

    T(U)(%)C(%)A(%)G(%)長Length(bp)GC(%)LSC32.618.531.217.78457336.2SSC33.916.734.115.31987432.0Ira28.222.628.320.82477443.4Irb28.320.828.322.62477443.4合計Total31.419.330.618.715399538.0CDS31.619.330.418.67828437.9

    注:CDS.蛋白編碼區(qū) 下同。

    Note:CDS.Protein-coding regions The same as below.

    表2鼠尾草葉綠體基因組各分區(qū)的基因組成

    Table2GenecompositionofLSC,SSC,IRintheS.japonicachloroplastgenome

    基因GeneCDStRNArRNA假基因PseudogeneLSC8563220SSC131210IRa197741psbaIRb19874合計Total135893781psba

    2.3 鼠尾草葉綠體基因組中的基因歸類

    鼠尾草葉綠體基因組全部基因功能分類如表3所示,它共含有18個基因內含子,其中3個基因(clpP、rpsl2、ycf3)含有2個內含子。rpsl2是一個反式剪接基因,5′端位于LSC區(qū),3′端在兩個IR區(qū)各有一個拷貝。

    圖1 鼠尾草葉綠體基因組環(huán)形基因圖 內側基因順時針;外側基因逆時針轉錄;不同功能的基因以不同顏色表示。Fig.1 Gene map of chloroplast genome in S.japonica Genes drawn inside the circle are transcribed clockwise,and those outsides are counter clockwise.Genes belonging to different functional groups are color-coded.

    基因分組Groupofgene基因名稱Genename1光系統(tǒng)ⅠPhotosysteinⅠpsaA,psaB,psaC,pasalI,psalJ2光系統(tǒng)ⅡPhotosystemⅡpsbA,psbB,psbC,psbD,psbE,psbF,psbH,psbl,psbJ,psbK,psbL,psbM,psbN,psbT,psbZ3細胞色素b/f復合體Cytochromeb/fcomplexpetA,petB?,petD?,petG,petL,petN4ATP合酶ATPsynthaseatpA,atpB,atpE,atpF?,atpH,atpI5NADH脫氫酶NADHdehydrogenasendhA?,ndhB(x2),ndhC,ndhD,ndhE,ndhF,ndhG,ndhH,ndhI,ndhJ,ndhK6二磷酸核酮糖羧化酶大亞基RubisCOlargesubunitrbcL7核糖核酸聚合酶RNApolymeraserpoA,rpoB,rpoC1?,rpoC28核糖體蛋白小亞基Ribosomalproteins(SSU)rps2,rps3,rps4,rps7(x2),rps8,rps11,rpsl2??(x2),rpsl4,rps15,rpsl16?,rpsl8,rpsl99核糖體蛋白小亞基Ribosomalproteins(LSU)rpl2?(x2),rpll4,rpl16?,rpl20,rpl22,rpl23(x2),rpl32,rpl33,rpl3610其他基因OthergenesclpP??,matK,accD,ccsA,infA,cemA11未知蛋白Proteinsofunknownfunctionycfl,ycf2(x2),ycf3??,ycf4,ycf15(x2)12轉運RNATransferRNAs37tRNAs(6containanintron,7intheIRs)13核糖體RNARibosomalRNAsrrn4.5(x2),rrn5(x2),rrn16(x2),rrn23(x2)

    注:*基因中有1或2個內含子

    Note:*gene contains one or two introns.

    表4鼠尾草葉綠體基因組內含子基因及外顯子、內含子長度

    Table4CharacterizationanalysisofgeneswithintronsinchloroplastgenomeofS.japonica

    基因Gene位置Location外顯子Exon(bp)內含子Intron(bp)ⅠⅡⅢⅠⅡatpFLSC411144704clpPLSC22829271612689ndhASSC5395531005ndhBIR756777675petBLSC6642730petDLSC8475708rpU6LSC19542863rpl2IR434391663rpoClLSC1620456755rps12?LSC11426232—539rpsl6LSC3399922trnA?UGCIR3835805trnG?UCCLSC1948688trnl?GAUIR3735948trnK?UUULSC35374390tmL?UAALSC3750483tmV?UACLSC3736593ycf3LSC153228129749699

    注:*rps12是一個反式剪接基因;5′端位于LSC區(qū);3′端位于兩個IR區(qū)

    Note:*therps12 is a trans-spliced gene with the 5′ end located in the LSC region and the duplicated 3′ end in the IR regions.

    2.4 鼠尾草葉綠體基因組內含子特征分析

    52.3%、1.8%、6.0%的基因組序列分別編碼蛋白、tRNA和rRNA基因,其余序列則歸為非編碼序列,包括內含子、基因間區(qū)及假基因。葉綠體基因組中內含子的基因特性如見表4所示,其中trnK-UUU內含子最長,多達4 390 bp,tmL-UAA內含子最短,為483 bp。

    2.5 鼠尾草葉綠體基因組密碼子使用情況

    分析整個葉綠體基因組的密碼子使用情況后發(fā)現(xiàn)共有2,797(10.6%)個密碼子編碼亮氨酸Leu,是編碼率最髙的氨基酸;而只有292(1.1%)個密碼子編碼半胱氨酸Cys,是編碼率最低的氨基酸(表5)。在所有密碼子中,相對同義密碼子使用頻率(relative synonymous codon usage,RSCU)最高的為UGG,其次為UAU、GAU,最低的為ATA、CTG和GTG。

    3 討論

    自地錢(Marchantiapolymorpha)完整葉綠體基因組發(fā)表以來[13],植物葉綠體基因組的研究得到國外內諸多學者的普遍關注,現(xiàn)成為葉綠體轉化技術和植物親緣關系鑒定的一個重要手段[14~15]。從先前的研究報道看,葉綠體基因組片段和核基因組序列已經(jīng)普遍用于鼠尾草屬植物的系統(tǒng)發(fā)育關系分析,并取得斐然成績。

    表5 鼠尾草葉綠體基因組密碼子使用情況

    注:相對同義密碼子使用頻率

    Note:Rrelative synonymous codon usage,RSCU

    譬如,Will等[16]通過ITS和ETS序列及葉綠體rpl32-trnL片段對非洲鼠尾草屬41個類群系統(tǒng)發(fā)育的分析,認為非洲鼠尾草屬物種并非單系起源(non-monophyletic);同樣,Walker等[17]利用葉綠體rbcL和trnL-F片段對美洲鼠尾草屬植物的研究,將其劃分為3個獨立的遺傳譜系分支。在中國,李虔全[18]通過ITS序列構建系統(tǒng)發(fā)育樹,認為中國鼠尾草屬植物與日本鼠尾草屬植物具有共同的祖先,新疆鼠尾草(S.deserta)與歐洲引進藥鼠尾草(S.officinalis)具有較近的親緣關系;然而,胡國雄[19]分子系統(tǒng)學的研究結果則表明鼠尾草屬是一個包含了4個獨立分支的多系類群,與Will等[20]的研究結果一致。盡管已有的分子系統(tǒng)學研究結果表明東亞鼠尾草屬為單系類群,但不支持國產鼠尾草屬的分類系統(tǒng),國產鼠尾草屬的分類系統(tǒng)不能通過合理地調整而與東亞鼠尾草的分子系統(tǒng)學研究結果相吻合,這可能與選取的研究樣本和所選片段有關[19]。鼠尾草是鼠尾草屬的一個廣布種和具有重要經(jīng)濟及藥用價值的物種,同時具有豐富的遺傳多樣性[12],因而加快對鼠尾草葉綠體基因組的分析,勢必可為今后利用基因遺傳標記、構建葉綠體超級條形碼,揭示亞洲地區(qū)乃至世界范圍內鼠尾草屬的進化歷史提供重要依據(jù)。目前,所有鼠尾草屬植物的系統(tǒng)發(fā)育分析都是基于核基因和葉綠體基因序列片段信息,然而因片段信息常常存在測定序列的分子進化速率與物種水平上的進化速率存在差異,以檢測序列位點的多樣性來代替物種的多樣性,這在一定程度上造成結果的不可靠性,基于基因組水平多態(tài)位點的檢測結果更能真實反映物種系統(tǒng)發(fā)育關系[8]。我們以丹參葉綠體基因組序列作為參考,對鼠尾草葉綠體基因組的研究表明,它的全長為153 995 bp,同多數(shù)已公布的被子植物葉綠體基因組大小范圍一樣[21],由一對反向重復區(qū)(Inverted repeat,IR)分隔一個大的單拷貝區(qū)(Large single copy region,LSC)和一個小的單拷貝區(qū)(Small single copy region,SSC),呈典型的四段式結構[6],鼠尾草葉綠體的GC含量為38.0%,這與先前報道的丹參的研究結果一致[11]。同時,鼠尾草葉綠體基因組大小同丹參葉綠體基因組相似,隸屬于是唇形目中較小的一族,葉綠體基因組大小的差異主要是因為LSC區(qū)的長度變異較大而導致的[11]。葉綠體基因組在基因組成及排列順序上具有高度保守性,根據(jù)功能可將葉綠體基因分為三大類[22],第一類與轉錄、轉譯有關;第二類與光合作用系統(tǒng)有關;第三類是與氨基酸、脂肪酸、色素等其它生物合成相關,鼠尾草基因組成與之類似。通常,葉綠體基因組中往往具有假基因現(xiàn)象,本研究所測的鼠尾草葉綠體基因組也不例外,同樣存在一些未知功能的基因,如ycf1基因,先前研究認為IRb/SSC邊界往往可以擴張到y(tǒng)cf1基因中,從而導致ycf1假基因的發(fā)生[23~25]。此外,先前研究還發(fā)現(xiàn),香豌豆屬(Lathyrus)[26]和棉屬(Gossypium)[27]植物的IR區(qū)都存在結構的變異,在進化過程中IR區(qū)部分的基因存在著丟失或萎縮現(xiàn)象,這些均有可能導致基因序列的變異,并極有可能致使物種進化事件的發(fā)生,因而我們推測它們很可能是鼠尾草葉綠體基因組中假基因產生的另一個重要原因之一。總之,本研究利用二代測序技術對鼠尾草的葉綠體基因組進行測序,并組裝得到完整的葉綠體基因組序列,是繼丹參后第二個鼠尾草屬植物葉綠體基因組的建成,不僅可為今后鼠尾草種質資源的綜合評價奠定了基礎,而且也豐富了鼠尾草屬的葉綠體基因組數(shù)據(jù)庫,為種質資源的合理開發(fā)利用提供基礎性數(shù)據(jù)。

    1.魏宇昆,王琦,黃艷波.唇形科鼠尾草屬的物種多樣性與分布[J].生物多樣性,2015,23(1):3-10.

    Wei Y Q,Wang Q,Huang Y B.Species diversity and distribution ofSalvia(Lamiaceae)[J].Biodiversity Science,2015,23(1):3-10.

    2.國家中醫(yī)藥管理局《中華本草》編委會.中華本草[M].上海:上??茖W技術出版社,1999.

    State Administration of Traditional Chinese Medicine.Chinese herbal medicine[M].Shanghai:Shanghai Science and Technology Publisher,1999.

    3.Isom T,Noda Y.Studies on the chemical constituents ofSalviajaponicaandS.glabrescens[J].Nippon Kagaku Zassi,1991(3):244-246.

    4.馮珂,趙鑫,鄭麗玲,等.3種鼠尾草屬植物有效成分分析[J].西北林學院學報,2010,25(6):140-143.

    Feng K,Zhao X,Zheng L L,et al.Analysis of active ingredients in threeSalviaplants[J].Journal of Northwest Forestry University,2010,25(6):140-143.

    5.Wester P,Claβen-bockhoff R.Pollination syndromes of new worldSalviaspecies with special reference to bird pollination[J].Annals of the Missouri Botanical Garden,2011,98(1):101-155.

    6.Bentham G.Labiatarum,genera et species[M].London:Ridgway,1832.

    7.陳士林,宋經(jīng)元,姚輝,等.藥用植物DNA條形碼鑒定策略及關鍵技術分析[J].中國天然藥物,2009(5):322-327.

    Chen S L,Song J Y,YAO H,et al.Strategy and key technique of identification of Chinese herbal medicine using DNA barcoding[J].Chinese Journal of Natural Medicines,2009(5):322-327.

    8.鄧科君,張勇,鄢丹,等.鼠尾草屬部分物種AFLP指紋圖譜構建及遺傳多樣性分析[J].植物遺傳資源學報,2010,11(4):483-486.

    Deng K J,Zhang Y,Yan D,et al.Construction of AFLP fingerprint and genetic diversity analysis of the genusSalvia[J].Journal of Plant Genetic Resources,2010,11(4):483-486.

    9.Clegg M T,Gaut B S,Learn G H J,et al.Rates and patterns of chloroplast DNA evolution[J].Proceedings of the National Academy of Sciences of the United States of America,1994,91(15):6795-6801.

    10.林小涵,劉志華,李卿,等.藥用植物葉綠體基因組研究[J].世界科學技術—中醫(yī)藥現(xiàn)代化,2010,12(3):442-446.

    Lin X H,Liu Z H,Li Q,et al.Chloroplast genome of medicinal plants[J].World Science & Technology,2010,12(3):442-446.

    11.Qian J,Song J Y,Gao H H,et al.The complete chloroplast genome sequence of the medicinal plantSalviamiltiorrhiza[J].PLoS One,2013,8(2):e57607.

    12.Sudarmono S,Okada H.Genetic variation ofSalviajaponicathunb.(Lamiaceae)and its related species[J].Floribunda Journal,2011,4(2):33-39.

    13.Umesono K,Inokuchi H,Ohyama K,et al.Nucleotide sequence ofMarchantiapolymorpha chloroplast DNA:a region possibly encoding three tRNAs and three proteins including a homologue ofE.coliribosomal protein S14[J].Nucleic Acids Research,1984,12(24):9551-9565.

    14.Li X W,Gao H H,Wang Y T,et al.Complete chloroplast genome sequence ofMagnoliagrandifloraand comparative analysis with related species[J].Science China Life Sciences,2013,56(2):189-198.

    15.邢少辰,Liu C J.葉綠體基因組研究進展[J].生物化學與生物物理進展,2008,35(1):21-28.

    Xin S C,Liu C J.Progress in chloroplast genome analysis[J].Progress in Biochemistry and Biophysics,2008,35(1):21-28.

    16.Will M,Claβen-bockhoff R.Why Africa matters:evolution of Old WorldSalvia(Lamiaceae)in Africa[J].Annals of Botany,2014,114(1):61-83.

    17.Walker J,Sytsma K J,Wink M.Salvia(Lamiaceae) is not monophyletic:implications for the systematics,radiation,and ecological specializations ofSalviaand tribe mentheae[J].American Journal of Botany,2004,91(7):1115-1125.

    18.李虔全.中國鼠尾草屬植物的分子系統(tǒng)發(fā)育及其亞屬間關系的研究[D].包頭:內蒙古科技大學包頭醫(yī)學院,2013.

    Li Q Q.Phylogenetic relationships ofSalvia(Lamiaceae)in China:evidence from DNA sequence datasets[D].Baotou:Baotou Medical College,2013.

    19.胡國雄.唇形科鼠尾草屬東亞分支的分子系統(tǒng)學研究[D].北京:中國科學院大學,2015.

    HU G X.Phylogeny ofSalvia(Lamiaceae) in East Asia[D].Beijing:Chinese Academy of Sciences University,2015.

    20.Will M,Claβen-bockhoff R.Time to splitSalvias.l.(Lamiaceae)-New insights from Old WorldSalviaphylogeny[J].Molecular Phylogenetics and Evolution,2017,109:33-58.

    21.Yang M,Zhang X W,Liu G M,et al.The complete chloroplast genome sequence of date palm(PhoenixdactyliferaL.)[J].PLoS One,2010,5(9):e12762.

    22.Wakasugi T,Tsudzuki T,Sugiura M.The genomics of land plant chloroplasts:Gene content and alteration of genomic information by RNA editing[J].Photosynthesis Research,2001,70(1):107-118.

    23.Raubeson L A,Peery R,Chumley T W,et al.Comparative chloroplast genomics:analyses including new sequences from the angiospermsNupharadvenaandRanunculusmacranthus[J].BMC Genomics,2007,8:174.

    24.Goulding S E,Wolfe K H,Olmstead R G,et al.Ebb and flow of the chloroplast inverted repeat[J].Molecular and General Genetics,1996,252(1-2):195-206.

    25.Wang R J,Cheng C L,Chang C C,et al.Dynamics and evolution of the inverted repeat-large single copy junctions in the chloroplast genomes of monocots[J].BMC Evolutionary Biology,2008,8:36.

    26.Magee A M,Aspinall S,Rice D W,et al.Localized hypermutation and associated gene losses in legume chloroplast genomes[J].Genome Research,2010,20(12):1700-1710.

    27.Talat F.三個D組棉花葉綠體全基因組測序及其進化研究[D].北京:中國農業(yè)科學院,2013.

    Talat F.Complete chloroplast genome sequences of Three D genome cotton species and their evolutionary implications[D].Beijing:Chinese Academy of Agricultural Sciences,2013.

    Under the auspices of the National Natural Science Foundation of China(32170338);The Key Project of Natural Science Foundation of Shaanxi Preschool Normal College(2015ZDKJ005)

    introduction:HE Yi-Han(1987—),female,Ph.D Candidates,engaged in research of medicinal plant secondary metabolic and molecular regulation.

    date:2017-03-06

    CompleteSequenceAnalysisofChloroplastGenomeofSalviajaponica

    HE Yi-Han1,2HAN Li-Min3LIU Yu-Ping3TIAN Na3SU Xu2WANG Zhe-Zhi1*

    (1.Key Laboratory of the Ministry of Education for Medicinal Resources and Natural Pharmaceutical Chemistry,National Engineering Laboratory for Resource Development of Endangered Crude Drugs in Northwest of China,College of Life Sciences,Shaanxi Normal University,Xi’an 710119;2.College of Geography and Life Science,Qinghai Normal University,Xining 810008;3.College of Life Sciences and Food Engineering,Shaanxi Xueqian Normal University,Xi’an 710061)

    Salviajaponicais a perennial herb of the genusSalvia(Lamiaceae), which is of great medicinal and economic value. The chloroplast genome ofS.japonicawas sequenced using Illumina Hiseq platform. The complete chloroplast genome sequence was obtained bySalviamiltiorrhizachloroplast genome as a reference. The full-length chloroplast genome ofS.japonicawas 153 995 bp. The LSC was 84 573 bp and 19 874 bp in SSC, and the length of IR region was 24 774 bp. The chloroplast genome ofS.japonicasuccessfully annotated 13 groups of genes. The gene types, number and GC content were similar with other species of Lamiaceae. These results enrich the chloroplast data and also provide some valuable molecular data for reconstructing evolutionary history ofSalvia.

    Salvia;Salviajaponica;Chloroplast genome;gene structure;phylogenetics and evolution

    國家自然科學基金項目(32170338)資助;陜西學前師范學院自然科學基金重點項目(2015ZDKJ005)資助

    何懿菡(1987—),女,博士研究生,主要研究方向為藥用植物次生代謝及分子調控。

    * 通信作者:E-mail:zhezhiwang@outlook.com

    2017-03-06

    * Corresponding author:E-mail:zhezhiwang@outlook.com

    Q949.777.6

    A

    10.7525/j.issn.1673-5102.2017.04.013

    猜你喜歡
    鼠尾草內含子葉綠體
    線粒體核糖體蛋白基因中內含子序列間匹配特性分析
    不同方向內含子對重組CHO細胞中神經(jīng)生長因子表達的影響
    更 正
    生物學通報(2019年1期)2019-02-15 16:33:43
    絨毛栗色鼠尾草根化學成分的研究
    中成藥(2018年10期)2018-10-26 03:41:06
    內含子的特異性識別與選擇性剪切*
    鼠尾草酸對β-淀粉樣蛋白損傷的神經(jīng)干細胞增殖的影響
    中成藥(2017年7期)2017-11-22 07:33:50
    亞臨界水萃取康定鼠尾草中丹參酮類成分工藝的優(yōu)化
    中成藥(2017年10期)2017-11-16 00:50:37
    南方紅豆杉葉綠體非編碼序列PCR體系優(yōu)化及引物篩選
    鼠尾草酸對自由基誘導蛋白質氧化損傷的保護作用
    食品科學(2013年15期)2013-03-11 18:25:50
    茶樹葉綠體DNA的PCR-RFLP反應體系優(yōu)化
    食品科學(2013年6期)2013-03-11 18:20:13
    成人av一区二区三区在线看| 国产区一区二久久| 日韩成人在线观看一区二区三区| 一区二区三区高清视频在线| 亚洲欧美激情在线| www.精华液| 黄片大片在线免费观看| 亚洲一码二码三码区别大吗| a级毛片在线看网站| 中文字幕另类日韩欧美亚洲嫩草| 女生性感内裤真人,穿戴方法视频| 亚洲精品国产精品久久久不卡| 国产精品精品国产色婷婷| 欧美日韩乱码在线| 少妇裸体淫交视频免费看高清 | 无人区码免费观看不卡| 美女 人体艺术 gogo| 又黄又粗又硬又大视频| 日韩高清综合在线| 最近最新免费中文字幕在线| 黄片大片在线免费观看| 成人手机av| 国产欧美日韩综合在线一区二区| 久久午夜亚洲精品久久| 中文字幕人妻熟女乱码| 女警被强在线播放| 免费在线观看黄色视频的| 亚洲人成伊人成综合网2020| 国产区一区二久久| 国产精品 欧美亚洲| 国产成人啪精品午夜网站| 淫妇啪啪啪对白视频| 精品国内亚洲2022精品成人| 黄色毛片三级朝国网站| 国产亚洲精品av在线| 美国免费a级毛片| 亚洲av五月六月丁香网| 国产乱人伦免费视频| av免费在线观看网站| 亚洲成人免费电影在线观看| 18禁观看日本| www.自偷自拍.com| 波多野结衣高清无吗| 日本欧美视频一区| 日本 av在线| 亚洲精品久久成人aⅴ小说| 男女床上黄色一级片免费看| 日韩三级视频一区二区三区| 国产精品二区激情视频| 99在线人妻在线中文字幕| 久久中文字幕人妻熟女| 久久久久久久久久久久大奶| www.www免费av| 亚洲人成伊人成综合网2020| 国产亚洲欧美精品永久| 高清在线国产一区| 久久 成人 亚洲| 亚洲国产欧美网| av电影中文网址| 一级黄色大片毛片| 日韩中文字幕欧美一区二区| 午夜精品国产一区二区电影| 国产男靠女视频免费网站| 此物有八面人人有两片| 一级a爱片免费观看的视频| 国产又色又爽无遮挡免费看| 午夜福利一区二区在线看| 99精品欧美一区二区三区四区| 欧美激情极品国产一区二区三区| 久久久久久久午夜电影| 极品教师在线免费播放| 电影成人av| 给我免费播放毛片高清在线观看| 中文字幕人妻熟女乱码| 一个人观看的视频www高清免费观看 | 亚洲国产精品sss在线观看| 99re在线观看精品视频| 欧美老熟妇乱子伦牲交| 在线天堂中文资源库| 一个人免费在线观看的高清视频| 91成年电影在线观看| 国产精品爽爽va在线观看网站 | 美国免费a级毛片| 国产成+人综合+亚洲专区| 在线观看66精品国产| 女人爽到高潮嗷嗷叫在线视频| 久久人人爽av亚洲精品天堂| 精品久久久久久久人妻蜜臀av | 12—13女人毛片做爰片一| 黄色 视频免费看| 波多野结衣av一区二区av| 乱人伦中国视频| 久久人人精品亚洲av| 亚洲伊人色综图| 两个人免费观看高清视频| 国产精品二区激情视频| 成在线人永久免费视频| 中国美女看黄片| 国产麻豆69| 日韩欧美国产在线观看| 久久久久久久久中文| 少妇熟女aⅴ在线视频| 久久香蕉精品热| 亚洲全国av大片| 亚洲片人在线观看| 成人免费观看视频高清| 亚洲精品中文字幕一二三四区| 大香蕉久久成人网| 久久久精品国产亚洲av高清涩受| 午夜福利视频1000在线观看 | 悠悠久久av| 国产精品影院久久| 亚洲精品一区av在线观看| 在线永久观看黄色视频| 亚洲精品国产精品久久久不卡| 色精品久久人妻99蜜桃| 日本 av在线| 国产黄a三级三级三级人| 一夜夜www| 99riav亚洲国产免费| 国内精品久久久久久久电影| 精品第一国产精品| 韩国av一区二区三区四区| 777久久人妻少妇嫩草av网站| 91字幕亚洲| 精品一品国产午夜福利视频| 国产精品二区激情视频| 日本免费a在线| 亚洲成av人片免费观看| 中文字幕色久视频| 91大片在线观看| 成人欧美大片| 精品国产国语对白av| 精品少妇一区二区三区视频日本电影| 精品电影一区二区在线| 波多野结衣av一区二区av| 欧美激情 高清一区二区三区| 国产精品爽爽va在线观看网站 | 丁香欧美五月| 深夜精品福利| а√天堂www在线а√下载| 一区在线观看完整版| 日本三级黄在线观看| 又黄又粗又硬又大视频| 久久久久久免费高清国产稀缺| 国产精品香港三级国产av潘金莲| 亚洲精品粉嫩美女一区| 亚洲国产欧美一区二区综合| 非洲黑人性xxxx精品又粗又长| 国产在线精品亚洲第一网站| 一边摸一边抽搐一进一出视频| 久久婷婷成人综合色麻豆| 亚洲第一av免费看| 国产亚洲欧美在线一区二区| 亚洲国产欧美一区二区综合| 国产麻豆成人av免费视频| 国产高清videossex| 美国免费a级毛片| 91字幕亚洲| 国产欧美日韩一区二区三| 亚洲七黄色美女视频| 久久这里只有精品19| 亚洲情色 制服丝袜| 变态另类成人亚洲欧美熟女 | 亚洲人成电影免费在线| 亚洲第一电影网av| 美国免费a级毛片| 老司机靠b影院| 国产蜜桃级精品一区二区三区| 黄频高清免费视频| 99久久国产精品久久久| 多毛熟女@视频| 丝袜美足系列| 妹子高潮喷水视频| 国产av一区二区精品久久| 亚洲九九香蕉| 真人一进一出gif抽搐免费| 久久影院123| 亚洲欧美精品综合一区二区三区| 最好的美女福利视频网| 国产亚洲欧美在线一区二区| 日韩大码丰满熟妇| 午夜福利成人在线免费观看| 午夜精品久久久久久毛片777| 免费不卡黄色视频| 日韩大尺度精品在线看网址 | 久久亚洲真实| 成人国语在线视频| netflix在线观看网站| 日韩中文字幕欧美一区二区| 亚洲精品久久成人aⅴ小说| 最近最新中文字幕大全免费视频| 日本一区二区免费在线视频| 久久影院123| 色老头精品视频在线观看| 国产av在哪里看| 90打野战视频偷拍视频| 青草久久国产| 国产片内射在线| 国产精品99久久99久久久不卡| 欧美+亚洲+日韩+国产| 国产精品香港三级国产av潘金莲| 免费高清视频大片| 欧美日韩乱码在线| 欧美日韩福利视频一区二区| 亚洲av电影不卡..在线观看| 麻豆一二三区av精品| 亚洲精品在线美女| 亚洲国产欧美日韩在线播放| 亚洲第一青青草原| 久久久久久久久久久久大奶| 十分钟在线观看高清视频www| 日日夜夜操网爽| 精品一区二区三区视频在线观看免费| 97人妻精品一区二区三区麻豆 | 丝袜在线中文字幕| 国产精品久久久人人做人人爽| 免费看美女性在线毛片视频| av福利片在线| 国产精品一区二区免费欧美| 757午夜福利合集在线观看| 97碰自拍视频| 久久精品国产亚洲av香蕉五月| 欧美日韩乱码在线| 国产熟女午夜一区二区三区| 亚洲第一电影网av| 国产成+人综合+亚洲专区| 免费搜索国产男女视频| 丝袜美腿诱惑在线| 一级作爱视频免费观看| 性欧美人与动物交配| 无限看片的www在线观看| 99国产极品粉嫩在线观看| 国产精品日韩av在线免费观看 | 一级黄色大片毛片| 一区福利在线观看| 高清黄色对白视频在线免费看| 欧美色欧美亚洲另类二区 | 国产一卡二卡三卡精品| 久久中文字幕一级| 日韩高清综合在线| 校园春色视频在线观看| 亚洲午夜精品一区,二区,三区| 国产亚洲精品第一综合不卡| www.999成人在线观看| 禁无遮挡网站| 在线播放国产精品三级| 99国产精品一区二区三区| 一进一出抽搐gif免费好疼| 久久狼人影院| 中文字幕另类日韩欧美亚洲嫩草| tocl精华| 精品一区二区三区av网在线观看| 又黄又爽又免费观看的视频| 久久狼人影院| 欧美在线一区亚洲| www国产在线视频色| 中文字幕av电影在线播放| 欧美日韩精品网址| 最好的美女福利视频网| 国产精品日韩av在线免费观看 | 精品少妇一区二区三区视频日本电影| 精品人妻1区二区| 天天添夜夜摸| 18禁裸乳无遮挡免费网站照片 | 国产成人精品在线电影| 99riav亚洲国产免费| 亚洲天堂国产精品一区在线| 97人妻天天添夜夜摸| 搡老岳熟女国产| 久久精品国产清高在天天线| 久久午夜亚洲精品久久| 中文字幕久久专区| 天天一区二区日本电影三级 | 一区在线观看完整版| 国产精品一区二区免费欧美| 亚洲精品中文字幕在线视频| 日本 av在线| 啪啪无遮挡十八禁网站| 亚洲一区高清亚洲精品| 午夜日韩欧美国产| 午夜视频精品福利| АⅤ资源中文在线天堂| 亚洲色图 男人天堂 中文字幕| 天堂影院成人在线观看| 欧美精品啪啪一区二区三区| 久久久久久人人人人人| 一区二区日韩欧美中文字幕| 国产一区二区在线av高清观看| 99国产精品一区二区蜜桃av| 国产亚洲欧美精品永久| 日韩国内少妇激情av| 婷婷六月久久综合丁香| 日本 欧美在线| 老司机午夜福利在线观看视频| 亚洲国产精品久久男人天堂| 午夜福利影视在线免费观看| 999久久久国产精品视频| 变态另类丝袜制服| 日韩av在线大香蕉| 色尼玛亚洲综合影院| 欧美成人午夜精品| av视频免费观看在线观看| 国产精品av久久久久免费| 两人在一起打扑克的视频| 欧美日本视频| 好男人电影高清在线观看| 亚洲七黄色美女视频| 国产精品久久电影中文字幕| 一级黄色大片毛片| 1024香蕉在线观看| 亚洲色图 男人天堂 中文字幕| 午夜成年电影在线免费观看| 久久久久国产精品人妻aⅴ院| 国产高清激情床上av| 在线观看午夜福利视频| 国产av一区在线观看免费| 精品国产超薄肉色丝袜足j| 欧美日本中文国产一区发布| 亚洲欧美激情综合另类| 免费观看精品视频网站| 色婷婷久久久亚洲欧美| 亚洲国产精品999在线| av视频在线观看入口| 熟女少妇亚洲综合色aaa.| 久久久久国产精品人妻aⅴ院| 久久亚洲精品不卡| 中国美女看黄片| 两个人看的免费小视频| 大型黄色视频在线免费观看| av在线天堂中文字幕| www.精华液| 日韩中文字幕欧美一区二区| 欧美激情 高清一区二区三区| 欧美乱妇无乱码| 国产精品野战在线观看| 国产精品久久久久久亚洲av鲁大| 中文字幕另类日韩欧美亚洲嫩草| 两人在一起打扑克的视频| 不卡av一区二区三区| 亚洲久久久国产精品| bbb黄色大片| 久久精品国产综合久久久| 亚洲美女黄片视频| 在线观看午夜福利视频| 精品久久久久久久人妻蜜臀av | 长腿黑丝高跟| 日日干狠狠操夜夜爽| 欧美国产日韩亚洲一区| 亚洲欧美精品综合一区二区三区| 天堂动漫精品| 国产一区在线观看成人免费| 母亲3免费完整高清在线观看| 久久久久精品国产欧美久久久| 一边摸一边做爽爽视频免费| 十八禁网站免费在线| aaaaa片日本免费| 老司机在亚洲福利影院| 午夜a级毛片| 久久久久国内视频| 中文字幕久久专区| 男女午夜视频在线观看| 亚洲国产欧美一区二区综合| 免费观看精品视频网站| 女人高潮潮喷娇喘18禁视频| 久久久久九九精品影院| 久久国产亚洲av麻豆专区| 亚洲美女黄片视频| av超薄肉色丝袜交足视频| 国产1区2区3区精品| 精品一区二区三区视频在线观看免费| 亚洲一区中文字幕在线| 久久国产精品男人的天堂亚洲| 18美女黄网站色大片免费观看| 亚洲 欧美一区二区三区| 在线视频色国产色| 欧美乱妇无乱码| xxx96com| 久久久久久人人人人人| 久久久久久大精品| 成人三级做爰电影| 久久久久久免费高清国产稀缺| 51午夜福利影视在线观看| 国产精品一区二区免费欧美| 变态另类成人亚洲欧美熟女 | 欧美色视频一区免费| 色老头精品视频在线观看| 88av欧美| 无人区码免费观看不卡| 免费无遮挡裸体视频| 久久精品影院6| 久久草成人影院| 岛国在线观看网站| АⅤ资源中文在线天堂| 最新在线观看一区二区三区| 久久国产精品影院| 美女高潮到喷水免费观看| 久久这里只有精品19| 在线观看66精品国产| 色综合站精品国产| 狂野欧美激情性xxxx| 亚洲中文字幕日韩| 国产一区二区三区在线臀色熟女| 校园春色视频在线观看| 激情在线观看视频在线高清| 国产欧美日韩精品亚洲av| 嫁个100分男人电影在线观看| 久久午夜综合久久蜜桃| 欧美一区二区精品小视频在线| 99国产精品免费福利视频| bbb黄色大片| 亚洲午夜精品一区,二区,三区| 天堂√8在线中文| 老司机福利观看| 成人av一区二区三区在线看| 免费人成视频x8x8入口观看| 50天的宝宝边吃奶边哭怎么回事| 涩涩av久久男人的天堂| 18美女黄网站色大片免费观看| 色播在线永久视频| av欧美777| 99精品久久久久人妻精品| 一级a爱视频在线免费观看| 天天添夜夜摸| 制服人妻中文乱码| 村上凉子中文字幕在线| 午夜福利视频1000在线观看 | 国产成人精品在线电影| 91麻豆av在线| 好男人电影高清在线观看| 国产成人影院久久av| 久久人人爽av亚洲精品天堂| 亚洲国产精品sss在线观看| 久久午夜亚洲精品久久| 中文字幕精品免费在线观看视频| 精品乱码久久久久久99久播| 免费不卡黄色视频| 欧美日韩黄片免| 国产成人精品久久二区二区91| 一区二区三区国产精品乱码| 国产精品免费一区二区三区在线| 国产精品亚洲av一区麻豆| 精品一区二区三区视频在线观看免费| 电影成人av| 91九色精品人成在线观看| 国产一卡二卡三卡精品| 国产区一区二久久| 夜夜看夜夜爽夜夜摸| 少妇裸体淫交视频免费看高清 | 国产午夜精品久久久久久| 亚洲av第一区精品v没综合| 国产人伦9x9x在线观看| 婷婷精品国产亚洲av在线| www.www免费av| 色综合婷婷激情| 国产成人系列免费观看| 大陆偷拍与自拍| 国产精品 国内视频| 在线观看一区二区三区| 国产精品1区2区在线观看.| 亚洲av片天天在线观看| 久久精品亚洲熟妇少妇任你| 国产精品98久久久久久宅男小说| 久久青草综合色| 一区二区三区高清视频在线| 色精品久久人妻99蜜桃| 少妇熟女aⅴ在线视频| 亚洲人成77777在线视频| 精品国内亚洲2022精品成人| 欧美黑人精品巨大| 大码成人一级视频| 精品一区二区三区视频在线观看免费| 男女下面插进去视频免费观看| 亚洲国产欧美网| 黄色视频不卡| 欧美最黄视频在线播放免费| 一区二区日韩欧美中文字幕| 久久久久国产一级毛片高清牌| 午夜福利高清视频| a级毛片在线看网站| 午夜激情av网站| 国产av精品麻豆| 中文字幕av电影在线播放| 午夜福利一区二区在线看| 老司机福利观看| 亚洲欧美日韩无卡精品| 亚洲人成伊人成综合网2020| 免费在线观看影片大全网站| 这个男人来自地球电影免费观看| 中文字幕精品免费在线观看视频| 欧美黑人欧美精品刺激| 两人在一起打扑克的视频| 亚洲精品国产区一区二| 国产成人精品在线电影| 精品熟女少妇八av免费久了| 日本一区二区免费在线视频| 亚洲av电影不卡..在线观看| 国产免费av片在线观看野外av| 日本一区二区免费在线视频| 人人妻人人澡欧美一区二区 | bbb黄色大片| 亚洲成人免费电影在线观看| 男人的好看免费观看在线视频 | 亚洲av美国av| 欧美黑人欧美精品刺激| 1024视频免费在线观看| 18禁黄网站禁片午夜丰满| 中文字幕高清在线视频| 欧美日本亚洲视频在线播放| 乱人伦中国视频| 久久人妻熟女aⅴ| 一级,二级,三级黄色视频| x7x7x7水蜜桃| 波多野结衣巨乳人妻| 欧美日韩黄片免| 88av欧美| 欧美色欧美亚洲另类二区 | 久久香蕉国产精品| 99精品在免费线老司机午夜| 老汉色av国产亚洲站长工具| 亚洲午夜理论影院| 一个人免费在线观看的高清视频| 9191精品国产免费久久| 色播亚洲综合网| 欧美不卡视频在线免费观看 | 国产97色在线日韩免费| 日日夜夜操网爽| 久久久久国产精品人妻aⅴ院| 久久人妻福利社区极品人妻图片| 真人做人爱边吃奶动态| 男女床上黄色一级片免费看| 免费不卡黄色视频| 国产高清视频在线播放一区| 如日韩欧美国产精品一区二区三区| 欧美成人性av电影在线观看| 可以在线观看毛片的网站| 深夜精品福利| 午夜久久久在线观看| 亚洲专区中文字幕在线| 国产色视频综合| 黄色成人免费大全| 色综合站精品国产| 午夜久久久久精精品| 男女午夜视频在线观看| 日本一区二区免费在线视频| 12—13女人毛片做爰片一| 久久九九热精品免费| 成人国产综合亚洲| 久久久精品欧美日韩精品| 国产激情欧美一区二区| 男女床上黄色一级片免费看| 亚洲欧洲精品一区二区精品久久久| 精品国产美女av久久久久小说| 日本在线视频免费播放| 日韩欧美一区视频在线观看| 精品午夜福利视频在线观看一区| 亚洲第一av免费看| 熟女少妇亚洲综合色aaa.| 久久天堂一区二区三区四区| 国产精品永久免费网站| 亚洲国产精品999在线| 国产成人一区二区三区免费视频网站| 一进一出抽搐gif免费好疼| 国产欧美日韩一区二区精品| 免费在线观看影片大全网站| tocl精华| 免费在线观看影片大全网站| 国产极品粉嫩免费观看在线| 成人av一区二区三区在线看| 一级作爱视频免费观看| 日韩中文字幕欧美一区二区| 亚洲色图 男人天堂 中文字幕| 亚洲一区二区三区色噜噜| 真人一进一出gif抽搐免费| 精品不卡国产一区二区三区| 十八禁人妻一区二区| 亚洲熟女毛片儿| 午夜免费激情av| 搞女人的毛片| 亚洲无线在线观看| 欧美久久黑人一区二区| 欧美黄色片欧美黄色片| xxx96com| 一二三四社区在线视频社区8| 黄片播放在线免费| 99在线人妻在线中文字幕| 怎么达到女性高潮| 亚洲熟妇中文字幕五十中出| 成人18禁高潮啪啪吃奶动态图| 在线永久观看黄色视频| 精品午夜福利视频在线观看一区| 美女大奶头视频| 亚洲一区高清亚洲精品| 波多野结衣高清无吗| 久久久精品国产亚洲av高清涩受| 亚洲电影在线观看av| 真人一进一出gif抽搐免费| 亚洲 国产 在线| 久久性视频一级片| 精品国产超薄肉色丝袜足j| 午夜免费激情av| 亚洲色图av天堂| 这个男人来自地球电影免费观看| 日韩精品中文字幕看吧| 免费观看人在逋| 国产亚洲精品av在线| 国产精品自产拍在线观看55亚洲| 色播在线永久视频| 国产精品爽爽va在线观看网站 | 日韩欧美三级三区| 亚洲国产精品久久男人天堂| 亚洲片人在线观看| 亚洲第一欧美日韩一区二区三区| 国产精品久久久av美女十八| 一卡2卡三卡四卡精品乱码亚洲|