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

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

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

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

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

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

鼠尾草屬；鼠尾草；葉綠體基因組；基因結(jié)構(gòu)；系統(tǒng)進(jìn)化

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

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

1 材料與方法

1.1 測序材料

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

1.2 測序

將檢測合格后的DNA，隨機(jī)打斷成1～2 kB的片段，富集大小一致的片段，連接至建庫載體中構(gòu)建DNA文庫(DNA文庫制備試劑盒，利用Illumina公司產(chǎn)品)，然后利用Illumina HiSeq 2500基因組測序平臺(tái)測序。

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

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

2 結(jié)果與分析

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

使用Illumina Hiseq平臺(tái)對(duì)鼠尾草葉綠體DNA提取樣本進(jìn)行測序，測序數(shù)據(jù)達(dá)2.4 G，目前基因組序列已提交至GenBank數(shù)據(jù)庫(登錄號(hào)KY646163)。鼠尾草葉綠體基因組序列全長153,995 bp，呈典型的四段式結(jié)構(gòu)(圖1)，其中LSC區(qū)(large single copy region，大的單拷貝區(qū))長84 573 bp，SSC區(qū)(small single copy region，小的單拷貝區(qū))長19 874 bp，兩個(gè)IR區(qū)(inverted repeat，反向重復(fù)區(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ū)包含的四個(gè)高GC含量的rRNA基因。

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

整條葉綠體基因組序列共編碼135個(gè)基因，去除重復(fù)后為118個(gè)，包括89(去除重復(fù)83)個(gè)蛋白編碼基因、37(去除重復(fù)30)個(gè)tRNA、8(去除重復(fù)4)個(gè)rRNA(表2)。6個(gè)蛋白編碼基因、7個(gè)tRNA基因及4個(gè)rRNA基因在兩個(gè)IR區(qū)各有一個(gè)拷貝。LSC區(qū)包含63個(gè)蛋白編碼基因、22個(gè)tRNA基因，而SSC區(qū)只含有12個(gè)蛋白編碼基因和1個(gè)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合計(jì)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合計(jì)Total135893781psba

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

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

圖1 鼠尾草葉綠體基因組環(huán)形基因圖 內(nèi)側(cè)基因順時(shí)針；外側(cè)基因逆時(shí)針轉(zhuǎ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細(xì)胞色素b/f復(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轉(zhuǎn)運(yùn)RNATransferRNAs37tRNAs(6containanintron，7intheIRs)13核糖體RNARibosomalRNAsrrn4．5(x2)，rrn5(x2)，rrn16(x2)，rrn23(x2)

注：*基因中有1或2個(gè)內(nèi)含子

Note:*gene contains one or two introns.

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

Table4CharacterizationanalysisofgeneswithintronsinchloroplastgenomeofS.japonica

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

注：*rps12是一個(gè)反式剪接基因；5′端位于LSC區(qū)；3′端位于兩個(gè)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 鼠尾草葉綠體基因組內(nèi)含子特征分析

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

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

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

3 討論

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

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

注：相對(duì)同義密碼子使用頻率

Note:Rrelative synonymous codon usage,RSCU

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

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

19.胡國雄.唇形科鼠尾草屬東亞分支的分子系統(tǒng)學(xué)研究[D].北京:中國科學(xué)院大學(xué),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.三個(gè)D組棉花葉綠體全基因組測序及其進(jìn)化研究[D].北京:中國農(nóng)業(yè)科學(xué)院,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

國家自然科學(xué)基金項(xiàng)目(32170338)資助；陜西學(xué)前師范學(xué)院自然科學(xué)基金重點(diǎn)項(xiàng)目(2015ZDKJ005)資助

何懿菡(1987—)，女，博士研究生，主要研究方向?yàn)樗幱弥参锎紊x及分子調(diào)控。

* 通信作者：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