蔣 明，王軍峰，朱 晏，黃雯馨，應(yīng)夢(mèng)豪，馬佳瑩，戴晨宇

射干葉綠體基因組結(jié)構(gòu)、序列特征與系統(tǒng)發(fā)育分析

蔣 明1，王軍峰2，朱 晏1，黃雯馨1，應(yīng)夢(mèng)豪1，馬佳瑩1，戴晨宇1

1. 臺(tái)州學(xué)院生命科學(xué)學(xué)院，浙江 臺(tái)州 318000 2. 華東藥用植物園科研管理中心，浙江 麗水 323000

以射干為材料，在測(cè)序、組裝獲得葉綠體基因組的基礎(chǔ)上，明確其結(jié)構(gòu)、序列特征及系統(tǒng)發(fā)育關(guān)系。利用PE150雙末端策略進(jìn)行建庫(kù)測(cè)序，用NOVOPlasty組裝完整的葉綠體基因組，經(jīng)PCR驗(yàn)證邊界，借助生物信息學(xué)工具進(jìn)行序列分析和系統(tǒng)發(fā)育研究。射干的葉綠體基因組全長(zhǎng)為153 816 bp，大單拷貝區(qū)、反向重復(fù)區(qū)和小單拷貝區(qū)的長(zhǎng)度分別為83 143、26 214、18 245 bp。射干葉綠體基因組共有133個(gè)基因，編碼基因、tRNA和rRNA的數(shù)量分別為92、38和8；有2個(gè)拷貝，其中一個(gè)為假基因。系統(tǒng)發(fā)育分析結(jié)果表明，7種植物葉綠體基因組在發(fā)育樹(shù)上可分為4組，射干與同為鳶尾科的溪蓀聚為一組，支持率達(dá)100%。射干葉綠體基因組的組裝、序列分析和系統(tǒng)發(fā)育分析，為該藥用植物的遺傳結(jié)構(gòu)和遺傳多樣性研究奠定了基礎(chǔ)。

射干；葉綠體基因組；結(jié)構(gòu)；序列分析；鳶尾科

射干(L.) Redouté為鳶尾科（Iridaceae）射干屬多年生草本植物，具結(jié)節(jié)狀的根狀莖，葉片劍形，花橙紅色，花瓣上散生紫褐色斑點(diǎn)，分布于我國(guó)的吉林、山東、安徽、江蘇、浙江和福建等20多個(gè)省份[1]。射干分布廣泛，是一種十分常見(jiàn)的藥用植物，有著悠久的栽培和藥用歷史；射干的主要藥用部位為其根狀莖，它富含酚類(lèi)化合物，尤其是類(lèi)黃酮和異黃酮類(lèi)物質(zhì)，如射干苷、白射干素、鳶尾黃素和野鳶尾苷等[2]。射干具有清熱解毒、散結(jié)消炎和利咽消腫等功效，可用于治療扁桃腺炎、喉痹咽痛和腰痛等[3-5]。射干葉形優(yōu)美、花色亮麗、花期較長(zhǎng)，具有較好的觀賞價(jià)值，常用于盆栽、庭院美化或藥草園片植等，是一種觀賞和藥用兼用植物[6]。近年來(lái)，有關(guān)射干的研究主要集中在栽培技術(shù)、組培快繁、病蟲(chóng)害防治、化學(xué)成分、藥理作用和轉(zhuǎn)錄組等方面[7-11]。

葉綠體是綠色植物進(jìn)行光合作用和能量轉(zhuǎn)換的細(xì)胞器，它將太陽(yáng)能轉(zhuǎn)換為化學(xué)能，用于物質(zhì)積累和植物的生長(zhǎng)發(fā)育[12]。葉綠體在植物逆境防御中也起著十分重要的作用，在不良環(huán)境下，產(chǎn)生大量的活性氧（reactive oxygen species，ROS），ROS作為一種調(diào)控基因表達(dá)水平的信號(hào)分子，促使植物產(chǎn)生新的適應(yīng)[13-14]。葉綠體基因組以雙鏈環(huán)狀形式存在于葉綠體中，被子植物的葉綠體基因組大小為120～170 ?kb，測(cè)序和組裝比核基因組容易[15]。隨著測(cè)序技術(shù)的快速發(fā)展及測(cè)序費(fèi)用的降低，越來(lái)越多的葉綠體基因組得以注釋?zhuān)壳耙延?000多種植物完成葉綠體基因組的測(cè)序，它們?cè)谶z傳多樣性、基因組進(jìn)化、基因水平轉(zhuǎn)移、系統(tǒng)發(fā)育分析和群體遺傳學(xué)等方面得到了應(yīng)用[16-18]。目前，有關(guān)射干葉綠體基因組的測(cè)序、組裝和注釋等研究未見(jiàn)報(bào)道。本研究以射干葉片DNA為材料，在利用高通量測(cè)序的基礎(chǔ)上拼接葉綠體基因組，明確其序列特征、基因組結(jié)構(gòu)及與近緣特種間的系統(tǒng)發(fā)育關(guān)系，為后續(xù)開(kāi)展遺傳多樣性和群體遺傳學(xué)研究奠定基礎(chǔ)。

1 材料與儀器

1.1 材料

射干葉片采自浙江舟山東福山島，伴生植物有山菅(L.) DC.、小茄Thunb.、蘚狀景天Hemsl.和濱海前胡Thunb.等。采集健康葉片，裝入樣品袋后帶回實(shí)驗(yàn)室備用。

1.2 儀器

艾本德Eppendorf移液槍?zhuān)怀瑑艄ぷ髋_(tái)（蘇州安泰空氣技術(shù)有限公司）；ThinkPad P52移動(dòng)工作站；伯樂(lè)C1000型PCR儀（Bio-Rad公司，美國(guó)）；Covaris超聲波DNA破碎儀（Chromatin Shearing，美國(guó)）；Illumina HiSeq X Ten測(cè)序儀；北京六一DYY-12型電泳儀及電泳槽（北京市六一儀器廠）；伯樂(lè)Gel Doc XR+凝膠成像系統(tǒng)（Bio-Rad公司，美國(guó)）。

2 方法

2.1 DNA的提取和文庫(kù)構(gòu)建

在無(wú)菌研砵中加入適量液氮，用研棒將葉片磨成細(xì)粉末，再利用十六烷基三甲基溴化銨（Cetyltrimethylammonium bromide，CTAB）方法提取基因組DNA，經(jīng)電泳檢測(cè)后用于構(gòu)建文庫(kù)。用超聲波破碎儀將基因組DNA片段化，經(jīng)末端修復(fù)、添加A尾、兩端加測(cè)序接頭、產(chǎn)物純化及PCR擴(kuò)增等過(guò)程，完成文庫(kù)的構(gòu)建。

2.2 高通量測(cè)序

采用雙端（paired-end，PE）策略，用Illumina HiSeq X Ten高通量測(cè)序儀進(jìn)行測(cè)序，讀長(zhǎng)為2×150 bp，共獲得3.55 G原始數(shù)據(jù)。利用NGS QC Toolkit v2.3.3對(duì)原始數(shù)據(jù)進(jìn)行過(guò)濾，去除接頭和低質(zhì)量的序列[19]。最終共得到11 814 172條clean reads，20值達(dá)97.53%，30為92.73%，序列質(zhì)量較高，可用于后續(xù)的拼接和注釋。

2.3 葉綠體基因組的拼接和注釋

葉綠體基因組的拼接在ThinkPad P52移動(dòng)工作站上進(jìn)行，拼接采用NOVOPlasty程序[20]。利用在線工具DOGMA（dual organellar genoMe annotator）對(duì)序列進(jìn)行基因注釋?zhuān)W(wǎng)址為http:// dogma.ccbb.utexas.edu/，起始和終止密碼子通過(guò)手工方式進(jìn)行調(diào)整[21]。tRNA用tRNAscan-SE（http:// www.lowelab.ucsc.edu/tRNAscan-SE/）和ARAGORN預(yù)測(cè)[22-23]。OGDRAW（organellar genome draw）用于生成葉綠體基因組圈圖，網(wǎng)址為http://www. ogdraw. mpimp-golm.mpg.de[24]。

2.4 邊界序列的PCR克隆和測(cè)序驗(yàn)證

葉綠體基因組的4個(gè)邊界采用PCR方法進(jìn)行鑒定，根據(jù)拼接的草圖，共設(shè)計(jì)了4對(duì)引物，分別是YGUP1：5’-GGGCGAACCAAAAAGAATGAT- G-3’、YGDN1：5’-CTTTTGTAGCCAATCATTTAT- CGGG-3’；YGUP2：5’-GGTTATGGAAGAAGG- AACCGAGAA-3’、YGDN2：5’-GCTATTTCCTC-TGCTTGTATTGGT-3’；YGUP3：5’-CTATTTTA- CGTCTTTGCGCGC-3’、YGDN3：5’-CCGAGCT- CGGGTTATGGAAG-3’；YGUP4：5’-CTGTAGA- CCCACGGAAAAATGT-3’、YGDN4：5’-GGTA- GAGCCGGATCGAAGT-3’。

Eppendorf管中，依次加入15 μL ddH2O、2.0 μL 10×緩沖液、0.4 μL dNTPs（10 mmol/L）、0.3 μL的上游引物（20 μmol/L）、0.3 μL下游引物（20 μmol/L）、1.0 μL DNA模板（50 ng/μL）和0.5 μLDNA聚合酶（2 U/μL）。PCR反應(yīng)在伯樂(lè)C1000型PCR儀上進(jìn)行，程序?yàn)椋?4 ℃變性5 min；32個(gè)循環(huán)中，94 ℃、30 s，54.3 ℃、45 s，72 ℃、100 s，共32個(gè)循環(huán)。PCR產(chǎn)物經(jīng)電泳、割膠、回收和純化后，將其與p-GEM T-easy載體（Promega）連接，置于4 ℃過(guò)夜。把連接產(chǎn)物導(dǎo)入大腸桿菌DH5α感受態(tài)細(xì)胞，經(jīng)PCR檢測(cè)后，各取3份菌液測(cè)序。

2.5 系統(tǒng)發(fā)育分析

從NCBI數(shù)據(jù)庫(kù)中下載6條葉綠體基因組的序列，它們分別來(lái)自溪蓀Donn ex Horn.（KT626943）、龍須菜Kunth（KX790361）、蔻第百子蓮F. M. Leight（KX790363）、惠普爾絲蘭(Torr.) Trel.、鹿蔥Maxim.（NC_040164）和馬褂木(Hemsl.) Sargent.（KU170538）。其中的溪蓀為鳶尾科植物；蔻第百子蓮和鹿蔥屬石蒜科（Amaryllidaceae）；龍須菜和惠普爾絲蘭為百合科（Liliaceae）植物；而外類(lèi)群馬褂木隸屬木蘭科（Magnoliaceae）。

利用Geneious 11.1.5內(nèi)置的Find repeats工具預(yù)測(cè)葉綠體基因組的反向重復(fù)（inverted repeat，IR）序列、大單拷貝區(qū)（large single copy，LSC）和小單拷貝區(qū)（small single copy，SSC），并用EXCEL計(jì)算各自的GC值。Geneious 11.1.5軟件中的MAFFT 7.388程序用于葉綠體基因組的多重比對(duì)，導(dǎo)出的比對(duì)序列用PhyML-SMS（Smart Model Selection in PhyML）在線工具獲得最佳替代模型，最后生成最大似然（maximum likelihood，ML）樹(shù)[25-26]。利用Mega X軟件構(gòu)建最大簡(jiǎn)約（maximum parsimony）樹(shù)，自舉檢測(cè)值為1000。

3 結(jié)果與分析

3.1 葉綠體基因組的邊界驗(yàn)證

利用NOVOPlasty對(duì)射干葉片DNA的Clean reads進(jìn)行拼接，并通過(guò)PCR手段對(duì)4個(gè)邊界進(jìn)行克隆和測(cè)序驗(yàn)證。電泳結(jié)果表明，4對(duì)引物均能擴(kuò)增出單1條帶，PCR產(chǎn)物分別位于800、1300、1800、1300 bp處，大小與預(yù)期一致（圖1）。序列測(cè)定結(jié)果表明，4條序列的長(zhǎng)度分別為812、1328、1859、1325 bp，堿基排列與組裝完成的葉綠體基因組邊界序列完全一致。

M-Marker 1～4-LSC/IRb、IRB/SSC、SSC/IRa和IRa/LSC邊界序列

3.2 葉綠體基因組的特點(diǎn)

射干葉綠體基因組全長(zhǎng)為153 816 bp，具有一個(gè)典型的四分體結(jié)構(gòu)，即分別由LSC、SSC、IRa和IRb 4個(gè)部分組成（圖2）。LSC與SSC的大小為83 143 bp和18 245 bp，IRa和IRb長(zhǎng)度均為26 214 bp（表1）。利用Excel計(jì)算GC值，結(jié)果表明，反向重復(fù)區(qū)的GC值最大，達(dá)43.0%，LSC次之，為36.0%，而SSC的GC值最小，僅31.4%；射干整個(gè)葉綠體基因組的GC值為37.8%。

3.3 葉綠體基因的組成和特點(diǎn)

射干葉綠體基因組上共有133個(gè)基因，包括38個(gè)轉(zhuǎn)運(yùn)RNA（tRNA）、8個(gè)核糖體RNA（rRNA）、14個(gè)核糖體蛋白小亞基基因、11個(gè)核糖體蛋白大亞基基因、4個(gè)RNA聚合酶基因、12個(gè)NADH脫氫酶亞基基因、20個(gè)光系統(tǒng)I/光系統(tǒng)II亞基基因、6個(gè)細(xì)胞色素b/f復(fù)合物亞基基因和6個(gè)ATP合成酶亞基基因等。另有7個(gè)未知功能基因，它們是、、和，除外，其它各有2份拷貝。

rRNA基因中，、、和各有2份拷貝，分別位于2個(gè)反向重復(fù)區(qū)域，4種rRNA的長(zhǎng)度分別為103、121、1491、2810 bp。tRNA中，、、、、、、和各有2份拷貝，其余tRNA均只有一個(gè)。具有2份拷貝的基因還有、、、、和基因等，也有2份拷貝，但其中1個(gè)為假基因（表2）。射干葉綠體基因組中，大部分基因沒(méi)有內(nèi)含子，少部分具1～2個(gè)內(nèi)含子。、和基因有2個(gè)內(nèi)含子，、、、、、、、、、、、、和具1個(gè)內(nèi)含子（表2）。

內(nèi)圈深色部分為GC含量

表1 射干葉綠體基因組的堿基組成

3.4 基因組特征比較分析

從NCBI下載了溪蓀、龍須菜、蔻第百子蓮、惠普爾絲蘭、馬褂木和鹿蔥6種植物的葉綠體基因組，它們的全長(zhǎng)分別為152 408、156 875、157 055、157 832、159 429和158 459 bp。馬褂木葉綠體基因組的GC值最大，為39.2%，溪蓀次之，為38.0%，蔻第百子蓮的GC值最小，僅37.5%（表3）。馬褂木的LSC和SSC最長(zhǎng)，分別為87 766、18 997 bp，但它的IR較短，為26 333 bp。蔻第百子蓮的IR最長(zhǎng)，長(zhǎng)度為26 869 bp，鹿蔥次之，為26 764 bp，而溪蓀的IR最短，僅26 026 bp。

射干與其他6種植物的LSC/IRb、IRb/SSC、SSC/IRa和IRa/LSC邊界及基因分布如圖3所示，基因在邊界的排列情況基本相同。除馬褂木外，LSC/IRb邊界兩側(cè)分布有和基因，而馬褂木LSC/IRb兩側(cè)為和基因。和基因位于IRb/SSC的邊界，7種植物在該位置的基因全長(zhǎng)為899～1112 bp，顯著短于正常的基因，它們均為假基因。另一個(gè)位于SSC/IRa交界處，長(zhǎng)度為5276～5489 bp，均為正?；?。馬褂木IRa/LSC邊界兩側(cè)分布有和基因，而其他6種植物在該處為和基因。

表2 射干葉綠體基因組的基因

×2-2份拷貝 Ψ-假基因*-一個(gè)內(nèi)含子**-2個(gè)內(nèi)含子

×2-2 copies Ψ-pseudogene*-one intron**-2 introns

表3 7種植物葉綠體基因組的特征

3.5 系統(tǒng)發(fā)育分析

利用在線工具PhyML-SMS獲得最佳分子進(jìn)化模型，結(jié)果表明，7條葉綠體基因組的最佳模型為GTR＋G＋I(xiàn)，該模型的赤池信息標(biāo)準(zhǔn)（Akaike information criterion，AIC）與貝葉斯信息標(biāo)準(zhǔn)（Bayesian information criterion，BIC）分別為820 620.543 72和820 831.677 20，與其他模型相比，它們的數(shù)值最低。以馬褂木葉綠體基因組為外類(lèi)群，構(gòu)建系統(tǒng)發(fā)育樹(shù)（圖4）。7種植物葉綠體基因組在系統(tǒng)發(fā)育樹(shù)上可分為4組，石蒜科的鹿蔥和蔻第百子蓮聚于組I，支持率為100%；百合科的龍須菜和惠普爾絲蘭聚于組II，支持率達(dá)100%；鳶尾科的溪蓀和射干聚于組III，支持率也為100%；而外類(lèi)群馬褂木單獨(dú)處于分支IV。同時(shí)，利用Mega X構(gòu)建了最大簡(jiǎn)約樹(shù)，結(jié)果與最大似然樹(shù)基本一致，除龍須菜和惠普爾絲蘭的支持率為99%外，其余均達(dá)100%（圖4）。

圖3 射干葉綠體基因組的反向重復(fù)區(qū)域/單拷貝區(qū)域邊界

圖4 基于葉綠體全基因組構(gòu)建的系統(tǒng)發(fā)育樹(shù)

4 討論

被子植物葉綠體基因組的結(jié)構(gòu)通常十分保守，為雙鏈環(huán)狀，具2個(gè)反向重復(fù)序列IRa和IRb，它們將LSC和SSC隔開(kāi)，最終形成四分體結(jié)構(gòu)[27]。葉綠體基因組的長(zhǎng)度也十分保守，大部分陸生植物的葉綠體基因組大小在135～160 kb[28]。禾本科（Gramineae）植物剪股穎L.的葉綠體基因組大小為136 584 bp[29]；百合科洋蔥L.的葉綠體基因大小為153 538 bp[30]。有些植物的葉綠體基因組較小，如藥用植物木賊麻黃Bge.僅109 518 bp[31]；而天竺葵Bailey的葉綠體基因組大小達(dá)217～942 bp，在目前已完成測(cè)序的陸生植物中最大[32]。本研究以射干為材料，在高通量測(cè)序的基礎(chǔ)上，對(duì)其葉綠體基因組進(jìn)行了組裝，經(jīng)PCR克隆和測(cè)序驗(yàn)證，得到射干葉綠體基因組的全長(zhǎng)序列，它的長(zhǎng)度為153 816 bp。

葉綠體基因組的基因組成十分保守，通常包含130個(gè)左右的基因，這些基因的功能涉及光合作用、轉(zhuǎn)錄和翻譯等[33]。葉綠體基因組通常有114種基因，包括4種rRNA基因、30種tRNA基因和80種蛋白質(zhì)編碼基因[34]。在一些寄生或半寄生植物中，基因丟失現(xiàn)象十分普遍。山毛櫸寄生(L.) W. P. C. Barton的葉綠體基因組大小僅70 kb，大部分基因丟失，僅剩下42個(gè)，與光合作用和葉綠體呼吸相關(guān)的基因全部缺失[35]。在半寄生植物廣寄生(DC) Danser和桑寄生(Lecomte) Danser葉綠體基因組中，僅注釋到106個(gè)基因，包括66個(gè)蛋白編碼基因、28個(gè)tRNA、8個(gè)rRNA和4個(gè)假基因[36]。在射干葉綠體基因組中，共有4種rRNA基因，它們各有2份拷貝，分布在不同的IR區(qū)域，另有30種tRNA基因，其中的、和等各有2份拷貝。

植物在進(jìn)化過(guò)程中，葉綠體基因組LSC/IRb、IRb/SSC、SSC/IRa和IRa/LSC邊界常發(fā)生擴(kuò)張或收縮事件，使葉綠體基因組大小出現(xiàn)一定的差異，甚至發(fā)生假基因化。比較魯桑Perr.、蒙桑Schneid.、印度桑L.和川桑Schneid.的葉綠體基因組發(fā)現(xiàn)，它們的SSC/IRb交界處均存在1個(gè)假基因；在鼠尾草Thunb. SSC/IRb處的也是一個(gè)假基因。本研究中，也發(fā)現(xiàn)類(lèi)似的現(xiàn)象，在射干葉綠體基因組中共有2份拷貝，位于SSC/IRb邊界的為假基因，而另一份正常。葉綠體基因組可用于植物的進(jìn)化分析和親緣關(guān)系鑒定，目前已應(yīng)用于不同分類(lèi)級(jí)別的系統(tǒng)進(jìn)化分析。張慧等[37]利用15個(gè)野芝麻亞科（Lamioideae）65個(gè)共有葉綠體蛋白序列構(gòu)建最大似然法樹(shù)（maximum likelihood，ML），發(fā)現(xiàn)益母草(Laur.) S. Y. Hu和水蘇屬Linn.的親緣關(guān)系較近，大部分節(jié)點(diǎn)的支持率達(dá)100%。射干與其他6種植物葉綠體基因組的系統(tǒng)發(fā)育分析結(jié)果表明，射干與溪蓀的關(guān)系最為接近，支持率達(dá)100%。

射干種質(zhì)資源十分豐富，但不同種源藥材的質(zhì)量參差不齊。葉綠體基因組的組裝和序列分析，為開(kāi)發(fā)高分辨率的遺傳標(biāo)記提供了依據(jù)，也為后續(xù)開(kāi)展該植物的群體遺傳學(xué)和遺傳多樣性研究奠定了基礎(chǔ)。

利益沖突 所有作者均聲明不存在利益沖突

[1] 中國(guó)科學(xué)院中國(guó)植物志編輯委員會(huì). 中國(guó)植物志 [第16(1)卷] [M]. 北京: 科學(xué)出版社, 1985: 26.

[2] 張婧涵, 張曉瑞, 李國(guó)信, 等. 線性回歸色譜峰定位法在射干藥材多組分同時(shí)測(cè)定中的應(yīng)用 [J]. 藥物分析雜志, 2014, 34(7): 1149-1155.

[3] 廉偉偉, 熊維政, 賈玉梅, 等. 射干及其方劑臨床應(yīng)用探討 [J]. 中醫(yī)學(xué)報(bào), 2018, 33(11): 2184-2190.

[4] 王迪, 辛旭陽(yáng), 尤獻(xiàn)民, 等. 射干效用演變探析 [J]. 遼寧中醫(yī)藥大學(xué)學(xué)報(bào), 2015, 17(9): 67-69.

[5] 展銳, 焦正花, 王紅麗, 等. 射干的藥理作用研究概況 [J]. 甘肅中醫(yī), 2011, 24(1): 78-80.

[6] 王曉榮, 朱序弼, 劉金. 耐寒耐旱射干西部開(kāi)發(fā)佳卉 [J]. 中國(guó)花卉盆景, 2002(9): 33.

[7] 王繼東, 盧杉. 延慶區(qū)射干栽培技術(shù) [J]. 農(nóng)業(yè)科技通訊, 2019(1): 192-193.

[8] 燕晨宇, 張翔, 秦民堅(jiān). 射干真葉、花蕾誘導(dǎo)愈傷組織的研究 [J]. 中國(guó)野生植物資源, 2018, 37(5): 16-19.

[9] 莫雪梅, 焦明姚, 文永剛, 等. 射干葉枯病防治藥劑篩選試驗(yàn)初探 [J]. 耕作與栽培, 2017(4): 36-37.

[10] Li J Y, Ni G, Li L,. New iridal-type triterpenoid derivatives with cytotoxic activities from[J]., 2019, 83: 20-28.

[11] Tian M, Zhang X, Zhu Y,. Global transcriptome analyses reveal differentially expressed genes of six organs and putative genes involved in (iso) flavonoid biosynthesis in[J]., 2018, 9: 1160.

[12] Suzuki N, Koussevitzky S, Mittler R,. ROS and redox signalling in the response of plants to abiotic stress [J]., 2012, 35(2): 259-270.

[13] Liu X M, Zhou Y L, Xiao J W,. Effects of chilling on the structure, function and development of chloroplasts [J]., 2018, 9: 1715.

[14] Miller G, Shulaev V, Mittler R. Reactive oxygen signaling and abiotic stress [J]., 2008, 133(3): 481-489.

[15] Du Y P, Bi Y, Yang F P,. Complete chloroplast genome sequences of Lilium: Insights into evolutionary dynamics and phylogenetic analyses [J]., 2017, 7(1): 5751.

[16] Xiong A S, Peng R H, Zhuang J,. Gene duplication, transfer, and evolution in the chloroplast genome [J]., 2009, 27(4): 340-347.

[17] Daniell H, Lin C S, Yu M,. Chloroplast genomes: Diversity, evolution, and applications in genetic engineering [J]., 2016, 17(1): 134.

[18] Twyford A D, Ness R W. Strategies for complete plastid genome sequencing [J]., 2017, 17(5): 858-868.

[19] Patel R K, Jain M. NGS QC Toolkit: A toolkit for quality control of next generation sequencing data [J]., 2012, 7(2): e30619.

[20] Dierckxsens N, Mardulyn P, Smits G. NOVOPlasty: de novo assembly of organelle genomes from whole genome data [J]., 2017, 45(4): e18.

[21] Wyman S K, Jansen R K, Boore J L. Automatic annotation of organellar genomes with DOGMA [J]., 2004, 20(17): 3252-3255.

[22] Lowe T M, Eddy S R. tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence [J]., 1997, 25(5): 955-964.

[23] Laslett D, Canback B. ARAGORN, a program to detect tRNA genes and tmRNA genes in nucleotide sequences [J]., 2004, 32(1): 11-16.

[24] Lohse M, Drechsel O, Kahlau S,OrganellarGenomeDRAW: A suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets [J]., 2013, 41: W575-W581.

[25] Katoh K, Standley D M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability [J]., 2013, 30(4): 772-780.

[26] Guindon S, Dufayard J F, Lefort V,. New algorithms and methods to estimate maximum-likelihood phylogenies: Assessing the performance of PhyML 3.0 [J]., 2010, 59(3): 307-321.

[27] Bendich A J. Circular chloroplast chromosomes: The grand illusion [J]., 2004, 16(7): 1661-1666.

[28] Raubeson L A, Peery R, Chumley T W,. Comparative chloroplast genomics: Analyses including new sequences from the angiospermsand[J]., 2007, 8: 174.

[29] Saski C, Lee S B, Fjellheim S,. Complete chloroplast genome sequences of Hordeum vulgare,and, and comparative analyses with other grass genomes [J]., 2007, 115(4): 571-590.

[30] von Kohn C, Kie?kowska A, Havey M J. Sequencing and annotation of the chloroplast DNAs and identification of polymorphisms distinguishing normal male-fertile and male-sterile cytoplasms of onion [J]., 2013, 56(12): 737-742.

[31] Wu C S, Lai Y T, Lin C P,. Evolution of reduced and compact chloroplast genomes (cpDNAs) in gnetophytes: Selection toward a lower-cost strategy [J]., 2009, 52(1): 115-124.

[32] Chumley T W, Palmer J D, Mower J P,. The complete chloroplast genome sequence of: Organization and evolution of the largest and most highly rearranged chloroplast genome of land plants [J]., 2006, 23(11): 2175-2190.

[33] Daniell H, Lee S B, Grevich J,. Complete chloroplast genome sequences of,and comparative analyses with other[J]., 2006, 112(8): 1503-1518.

[34] 宋菊, 龍?jiān)录t, 林麗梅, 等. 五加科植物葉綠體基因組結(jié)構(gòu)與進(jìn)化分析 [J]. 中草藥, 2017, 48(24): 5070-5075.

[35] 李巧麗, 延娜, 宋瓊, 等. 魯桑葉綠體基因組序列及特征分析 [J]. 植物學(xué)報(bào), 2018, 53(1): 94-103.

[36] 何懿菡, 韓立敏, 劉玉萍, 等. 鼠尾草葉綠體基因組序列分析 [J]. 植物研究, 2017, 37(4): 572-578.

[37] 張慧, 何帥兵, 孔繁德, 等. 益母草葉綠體基因組序列與系統(tǒng)進(jìn)化位置分析 [J]. 中醫(yī)藥信息, 2018, 35(4): 21-27.

Structure, sequence characteristics, and phylogenetic evolution analysis ofchloroplast genome

JIANG Ming1, WANG Jun-feng2, ZHU Yan1, HUANG Wen-xin1, YING Meng-hao1, MA Jia-ying1, DAI Chen-yu1

1. College of Life Science, Taizhou University, Taizhou 318000, China 2. Scientific Research Management Center, East China Medicinal Botanical Garden, Lishui 323000, China

To confirm the genome structure, sequence characteristics, and phylogenetic relationship by sequencing and assembling the chloroplast genome of a medicinal plant.A PE150 strategy was applied to construct library. The complete chloroplast genome was generated using NOVOPlasty, followed by PCR confirmation of the borders, and sequence analysis, as well as phylogenetic study was conducted by bioinformatic tools.The full-length chloroplast genome was 153 816 bp in length, with a large single copy of 83 143 bp, an inverted repeat of 26 214 bp, and a small single copy of 18 245 bp. Thechloroplast genome consisted of 133 genes, including 92 protein-coding genes, 38 tRNAs, and 8 rRNAs, respectively. There were twogenes, one of which was a pseudogene. Phylogenetic evolution analysis results indicated that the seven chloroplast genomes can be divided into four groups,andfrom Iridaceae were found to cluster in the same clade, with a support rate of 100%.Assembly, sequence analysis and phylogenetic evolution ofchloroplast genome provides an insight into studies on both genetic structure and genetic diversity.

(L.) Redouté; chloroplast genome; structure; sequence analysis; Iridaceae

R282.12

A

0253 - 2670(2021)13 - 4039 - 08

10.7501/j.issn.0253-2670.2021.13.027

2020-12-03

臺(tái)州市211人才工程經(jīng)費(fèi)資助（2012年度）

蔣 明（1973—），男，浙江嵊州人，博士，教授，碩士生導(dǎo)師，研究方向?yàn)橹参锘蚪M學(xué)、植物逆境生物學(xué)及其分子調(diào)控。E-mail: jiangming1973@139.com

[責(zé)任編輯 時(shí)圣明]