姜汶君 郭夢(mèng)月 龐曉慧

摘要?葉綠體是細(xì)胞內(nèi)具有自主遺傳信息的重要細(xì)胞器，近年來(lái)葉綠體基因組已發(fā)展成為近緣物種分子標(biāo)記、揭示植物進(jìn)化和系統(tǒng)發(fā)育關(guān)系的有力工具。藥用植物是人類賴以生存和發(fā)展的自然資源，葉綠體基因組在藥用植物研究中的應(yīng)用價(jià)值受到廣泛關(guān)注。在總結(jié)葉綠體基因組結(jié)構(gòu)研究的基礎(chǔ)上，現(xiàn)就葉綠體基因組在藥用植物物種鑒定、系統(tǒng)進(jìn)化等方面的應(yīng)用進(jìn)行綜述，以期為藥用植物種質(zhì)資源的開發(fā)和利用提供參考。

關(guān)鍵詞?藥用植物;葉綠體基因組;物種鑒定;DNA條形碼;高變異區(qū);系統(tǒng)進(jìn)化;分歧時(shí)間;基因工程

Application of Chloroplast Genome in Identification and Phylogenetic Analysis of Medicinal Plants

JIANG Wenjun， GUO Mengyue， PANG Xiaohui

（Key Lab of Chinese Medicine Resources Conservation， State Administration of Traditional Chinese Medicine of the People′s Republic of China， Institute of Medicinal Plant Development， Chinese Academy of Medical Sciences and Peking Union Medical College， Beijing 100193， China; Engineering Research Center of Chinese Medicine Resources， Ministry of Education， Beijing 100193， China）

Abstract?Chloroplast is an important organelle with autonomous genetic information in the cell.In recent years， the chloroplast genome has been used as a powerful tool for molecular markers of related species and uncovering plant evolution and phylogeny.Medicinal plants are natural resources that human beings depend on for survival and development， and the application value of chloroplast genomes in the research of medicinal plants has been widely concerned.The paper reviewed the application of the chloroplast genome in species identification and phylogenetic evolution of medicinal plants based on summarizing the research of chloroplast genome structure， providing a reference for the development and utilization of germplasm resources of medicinal plants.

Keywords?Medicinal plants; Chloroplast genome; Highly divergent regions; Species identification; DNA barcode; Phylogenetic analysis; Divergence times; Genetic engineering

中圖分類號(hào)：R284文獻(xiàn)標(biāo)識(shí)碼：Adoi：10.3969/j.issn.1673-7202.2020.05.008

葉綠體（Chloroplast）被普遍認(rèn)為起源于內(nèi)共生的藍(lán)藻[1]，不僅是細(xì)胞內(nèi)具有自主遺傳信息的半自主細(xì)胞器，還是植物進(jìn)行光合作用和生物合成的關(guān)鍵場(chǎng)所。近年來(lái)隨著新一代測(cè)序技術(shù)的迅速發(fā)展和測(cè)序成本的下降，已發(fā)表的葉綠體基因組數(shù)量急劇增加。截至2019年12月，GenBank共收錄超過(guò)3 300條葉綠體基因組序列，距離煙草和地錢葉綠體基因組（Chloroplast Genome）的首次成功測(cè)序僅有短短30余年[2-3]。葉綠體基因組包含豐富的遺傳信息，以結(jié)構(gòu)、大小和基因種類均較保守而著稱[4]，被廣泛應(yīng)用于物種鑒定及系統(tǒng)進(jìn)化研究，對(duì)于篩選藥用植物密切相關(guān)物種的分子標(biāo)記以及破譯親緣關(guān)系較近植物類群間的系統(tǒng)發(fā)育關(guān)系具有重要貢獻(xiàn)。本文對(duì)近年來(lái)葉綠體基因組的結(jié)構(gòu)研究以及其在藥用植物鑒定及系統(tǒng)進(jìn)化研究中的應(yīng)用進(jìn)行了綜述，分析了研究中存在的問(wèn)題，同時(shí)對(duì)葉綠體基因組研究的發(fā)展方向提出展望。

1?葉綠體基因組結(jié)構(gòu)研究

得益于高通量測(cè)序技術(shù)和生物信息技術(shù)的發(fā)展，越來(lái)越多物種的葉綠體基因組被成功測(cè)序和組裝，研究表明大多數(shù)陸生植物的葉綠體基因組基本上是被2個(gè)反向重復(fù)序列（Inverted Repeat Sequence，IR）分成大單拷貝區(qū)（Large Single Copy，LSC）和小單拷貝區(qū)（Small Single Copy，SSC）的環(huán)狀四分體結(jié)構(gòu)[4]。伴隨著IR區(qū)的擴(kuò)展和收縮，葉綠體基因組大小通常為120～160 kb[4-5]。葉綠體基因組共編碼120～130個(gè)基因，主要參與光合作用的光反應(yīng)、轉(zhuǎn)錄和翻譯等相關(guān)功能[5-6]。盡管短片段倒位（Small Inversion）、插入/缺失（Insertion/Deletion or Indel）、基因組結(jié)構(gòu)重排、基因轉(zhuǎn)移與丟失等葉綠體結(jié)構(gòu)變異事件已被報(bào)道[7]，但葉綠體基因組編碼基因種類、排列順序和含量通常被認(rèn)為是保守的。

蛋白編碼基因和tRNA的數(shù)量非常相似，基因組大小的差異主要來(lái)源于IR區(qū)基因數(shù)量的不同[5]。Guo等發(fā)現(xiàn)與其他被子植物相比，五味子ycf15基因丟失，IR區(qū)明顯收縮10 kb，但I(xiàn)R區(qū)收縮機(jī)制尚不明確，收縮區(qū)的GC含量明顯低于IR區(qū)，推測(cè)IR區(qū)收縮可能與GC含量相關(guān)[8];與2種其他傘形科藥用植物相比，濱海前胡的IR區(qū)邊界明顯擴(kuò)增，并且具有一個(gè)625 bp的含有trnD-GUC、trnY-GUA、trnE-UUC基因的短片段倒位[9];此外也有部分物種丟失了IR區(qū)，例如豆科植物黃耆只有一個(gè)反向重復(fù)序列[10]，屬于蝶形花亞科中的反向重復(fù)序列缺失分支（Inverted Repeat-Lacking Clade）。大多數(shù)寄生植物的葉綠體基因組中存在基因丟失現(xiàn)象[11]，這與其光合作用能力的部分或全部退化、從寄主上吸收養(yǎng)分和水分的生活方式有關(guān)。例如廣寄生和桑寄生，與煙草相比其葉綠體基因組丟失了ndh基因、infA基因、3個(gè)rpl基因（rpl32、rps15和rps16）、7個(gè)tRNA基因（trnA-UGC、trnG-UCC、trnH-GUG、trnL-GAU、trnK-UUU、trnL-UAA和trnV-UAC）和4個(gè)ycf基因（ycf1、ycf5、ycf9和ycf10）[12]，研究表明ndh基因在從自養(yǎng)型向異養(yǎng)型轉(zhuǎn)化過(guò)程中最先丟失[13]。其他自養(yǎng)型物種在進(jìn)化過(guò)程中也存在基因丟失、短片段倒位等結(jié)構(gòu)變異現(xiàn)象，其中infA基因的丟失在被子植物中最為普遍[14]。黃耆葉綠體基因組中存在3個(gè)基因（infA、rps16、rpl22）丟失，其在rps16基因和rbcL基因之間有一個(gè)50 kb的倒位，這是豆科植物中廣泛存在的現(xiàn)象[15]，此外還發(fā)現(xiàn)其ndhF基因和ycf1基因之間有一個(gè)20 kb的倒位[10]。黃花蒿具有一個(gè)完全倒置的SSC區(qū)（18 267 bp），這一現(xiàn)象在其他菊科蒿屬植物中也曾被報(bào)道，此外rps19基因在LSC區(qū)而非IR區(qū)，rps19假基因在IR區(qū)缺失[16]。菊科橐吾屬藥用植物葉綠體基因組的LSC區(qū)具有一個(gè)約3.4 kb的小倒位嵌套在23 kb的大倒位中，這也是菊科植物的一個(gè)獨(dú)特特征[17]。He等報(bào)道忍冬葉綠體基因組在trnI-CAU和trnN-GUU之間有一個(gè)獨(dú)特的重排現(xiàn)象，ycf1假基因丟失，rps19、rpl2和rpl23基因從IR區(qū)移動(dòng)到LSC區(qū)，并首次發(fā)現(xiàn)ycf2和rps18基因中含有內(nèi)含子[18]。

2?藥用植物葉綠體基因組的應(yīng)用

2.1?物種鑒定

近年來(lái)，DNA條形碼（DNA Barcoding）已發(fā)展成為物種鑒定的強(qiáng)有力工具。在植物物種鑒定中，常用的DNA片段包括葉綠體基因組序列rbcL、matK、psbA-trnH等與核基因組的ITS和ITS2[19]。其中，鑒于ITS2序列較強(qiáng)的鑒定能力，已被建議作為藥用植物的通用DNA條形碼[20]。然而，對(duì)于一些特殊藥用植物類群，普通DNA條形碼不具備足夠的變異信息以實(shí)現(xiàn)物種鑒定，例如淫羊藿屬和貝母屬[21-22]。與普通的DNA條形碼短片段相比，葉綠體全基因組包含更豐富的變異位點(diǎn)，鑒定效率更高。葉綠體全基因組一方面被直接用作物種鑒定的超級(jí)條形碼（Super Barcode），另一方面則基于其篩選出高變區(qū)作為鑒定物種的潛在分子標(biāo)記[23-25]。

2.1.1?超級(jí)條形碼?Cui等比對(duì)分析了砂仁、縮砂仁、海南砂仁等32種豆蔻屬葉綠體基因組，結(jié)果表明葉綠體全基因組可準(zhǔn)確鑒定豆蔻屬物種[26]。以往研究表明ITS、matK、psbA-trnH和rbcL對(duì)橐吾屬的物種鑒定效率不理想[27]，Chen等報(bào)道了6種橐吾屬植物的葉綠體基因組并提出以全基因組作為超級(jí)條形碼可有效鑒定橐吾屬植物[17]。鐘志敏發(fā)現(xiàn)基于葉綠體全基因組不僅可以有效鑒定石斛屬41個(gè)物種，還可區(qū)分來(lái)自3個(gè)不同產(chǎn)地的鐵皮石斛[28]。Xia等的研究表明了基于葉綠體全基因組可將野菊和近緣物種準(zhǔn)確區(qū)分開，可作為鑒定這些物種的超級(jí)條形碼[29]。

2.1.2?高變區(qū)篩選?李雪佩等比較分析了黃草烏和展毛黃草烏的葉綠體全基因組，發(fā)現(xiàn)trnR-atpA和atpF-atpH基因間隔區(qū)上的缺失可作為分子標(biāo)記鑒定2個(gè)物種[30]。葉綠體InDel分子標(biāo)記也被證明是鑒定白首烏[31]、6種忍冬屬植物[32]和人參屬[33]的有效工具。馬雙姣基于測(cè)定的薯蕷和叉蕊薯蕷以及從GenBank下載的4條薯蕷屬植物葉綠體基因組比對(duì)分析，篩選出了5個(gè)蛋白編碼基因和5個(gè)基因間隔區(qū)作為潛在分子標(biāo)記[34]。Chen等通過(guò)對(duì)川貝母及其近緣物種的8條葉綠體全基因組進(jìn)行比對(duì)分析，篩選出6個(gè)高變區(qū)作為鑒定貝母屬植物的潛在DNA條形碼，并提出葉綠體全基因組作為超級(jí)條形碼可能更適合該屬的物種鑒定和系統(tǒng)發(fā)育分析[35]。楊俏俏等通過(guò)比對(duì)分析6種蔥屬植物的葉綠體基因組，共篩選到7個(gè)變異較大的候選DNA條形碼序列用于鑒定蔥屬植物[36]。李曉娟等發(fā)現(xiàn)華重樓葉綠體基因組中cemA基因是假基因，該基因編碼序列中的終止密碼子位置差異可作為區(qū)分華重樓和北重樓的依據(jù)[37]。Do等基于黃瓜假還陽(yáng)參的葉綠體基因組開發(fā)了新的來(lái)自于rbcL和matK基因的SNP分子標(biāo)記，可將其與相關(guān)物種進(jìn)行有效區(qū)分[38]。Chen等基于3個(gè)麻黃屬藥用物種的葉綠體基因組篩選出11個(gè)高變區(qū)，可作為鑒定麻黃屬植物的潛在DNA條形碼[39]。Zhou等基于掌葉大黃、藥用大黃、唐古特大黃的葉綠體基因組開發(fā)了包含鑒定位點(diǎn)的潛在DNA條形碼并成功驗(yàn)證，為超級(jí)DNA條形碼在大黃中的應(yīng)用奠定了基礎(chǔ)[40]。

2.2?系統(tǒng)進(jìn)化

葉綠體基因組在大多數(shù)被子植物中為母系遺傳，重組率低、核苷酸置換率適中[41]，具有良好的系統(tǒng)發(fā)育重建潛力。葉綠體基因組的遺傳信息在解決系統(tǒng)發(fā)育位置和不同物種間的親緣關(guān)系方面做出了重要貢獻(xiàn)，同時(shí)也被廣泛用于揭示物種起源和估計(jì)各譜系間的分歧時(shí)間。

2.2.1?系統(tǒng)發(fā)育位置?魏似婕等基于藏藥“巴萵色保”的基原植物獐牙菜屬輪葉獐牙菜、川西獐牙菜及6種龍膽屬植物葉綠體全基因組構(gòu)建ML樹，結(jié)果支持經(jīng)典分類學(xué)觀點(diǎn)對(duì)2個(gè)屬物種的分類[42]。沈立群基于實(shí)驗(yàn)獲得的薄荷、白毛夏枯草和活血丹3種唇形科藥用植物的葉綠體基因組及從NCBI下載的唇形科序列構(gòu)建了系統(tǒng)發(fā)育樹，3種植物在進(jìn)化樹中的位置均被支持且與傳統(tǒng)分類學(xué)相符[43]。張慧等基于ML系統(tǒng)進(jìn)化分析確定了益母草在野芝麻亞科中的進(jìn)化位置，發(fā)現(xiàn)益母草屬和水蘇屬的親緣關(guān)系較近[44]。在Chen等的研究中，基于對(duì)川貝母及其密切相關(guān)物種的系統(tǒng)發(fā)育分析顯示，8個(gè)貝母屬植物聚成3支且均有較強(qiáng)的支持率，川貝母與甘肅貝母、中華貝母密切相關(guān)[35]。Chen等報(bào)道基于葉綠體基因組構(gòu)建的ML和MP系統(tǒng)發(fā)育樹可以鑒定3種麻黃屬植物，同時(shí)發(fā)現(xiàn)麻黃屬植物與買麻藤屬植物、百歲蘭具有更近的系統(tǒng)發(fā)育關(guān)系[39]。Cui等基于已發(fā)表姜目物種的葉綠體基因組序列，探討了姜在姜目中的系統(tǒng)發(fā)育位置，結(jié)果表明姜屬為山柰屬的姊妹分支[45]。Gao等分析了益智的葉綠體基因組并構(gòu)建了姜科的系統(tǒng)發(fā)育樹，結(jié)果強(qiáng)烈支持姜科植物喜陰和喜陽(yáng)的特征分類[46]。Gichira等對(duì)非洲一種隸屬于薔薇科的藥用植物Hagenia abyssinica的葉綠體基因組進(jìn)行了測(cè)定并構(gòu)建了系統(tǒng)發(fā)育樹，結(jié)果支持將該屬分在薔薇亞科下[47]。Guo等首次報(bào)道了五味子的葉綠體基因組序列，通過(guò)構(gòu)建ML和BI系統(tǒng)發(fā)育樹確定了五味子屬與八角屬形成姊妹群[8]。Gu等測(cè)定了藥用植物黃薇的葉綠體基因組，基于12種桃金娘目物種共有的68個(gè)CDS編碼基因構(gòu)建了系統(tǒng)發(fā)育樹，結(jié)果表明其與紫薇屬物種的親緣關(guān)系更近[48]。Li等基于廣寄生、桑寄生等10個(gè)物種的54個(gè)蛋白編碼基因和matK基因構(gòu)建的NJ和ML系統(tǒng)發(fā)育樹拓?fù)浣Y(jié)構(gòu)基本一致，結(jié)果有力支持桑寄生科和槲寄生科為獨(dú)立的2個(gè)單系類群[12]。

2.2.2?分歧時(shí)間估計(jì)?Lee等以Daucus carata和人參為標(biāo)定點(diǎn)，基于14個(gè)傘形科物種的76個(gè)葉綠體基因組共有蛋白編碼基因構(gòu)建貝葉斯進(jìn)化樹，結(jié)果表明防風(fēng)和濱海前胡的分歧時(shí)間大約在2.3～3.4百萬(wàn)年前（Million Years Ago，Mya），而珊瑚菜與紫花前胡、大葉山芹的分歧時(shí)間大約在1.6～2.8 Mya[9]。He等基于8種唇形科植物及外類群的葉綠體基因組構(gòu)建系統(tǒng)發(fā)育樹，采用分子鐘模型估算廣藿香分歧時(shí)間，結(jié)果表明其大約在29.45 Mya從黃芩亞科分支分化[49]。采用最大似自然法對(duì)21個(gè)毛茛目物種和2個(gè)木蘭科物種的葉綠體基因組進(jìn)行了系統(tǒng)發(fā)育推斷分析，He等發(fā)現(xiàn)黃連與其他毛茛科物種的分歧時(shí)間約為81 Mya，而毛茛科和小檗科的分歧時(shí)間約為111 Mya，在侏羅紀(jì)時(shí)期（大約153 Mya）兩者有一個(gè)共同的祖先[50]。Kim等基于葉綠體基因組比較分析了五加科的系統(tǒng)發(fā)育關(guān)系，樹參屬與五加屬的分歧時(shí)間大約在4.48～5.60 Mya，人參屬與楤木屬的分歧時(shí)間大約在2.58～3.20 Mya，人參和西洋參的分歧時(shí)間大約在0.72～0.87 Mya[51]。Fan等基于葉綠體基因組的蛋白編碼序列計(jì)算推斷了14個(gè)川續(xù)斷目物種的分歧時(shí)間，五?；婆c忍冬科的分歧時(shí)間大約為81.14 Mya，莢蒾屬與華福花屬、五福花屬、四福花屬3個(gè)屬的分歧時(shí)間大約為47.16 Mya，華福花屬、五福花屬、四福花屬的分歧時(shí)間大約在15～16 Mya[52]。

3?小結(jié)與展望

新一代測(cè)序技術(shù)的出現(xiàn)和進(jìn)步是葉綠體基因組學(xué)和遺傳學(xué)領(lǐng)域快速發(fā)展的重要因素之一，眾多植物葉綠體全基因組序列的研究和報(bào)道，不僅增強(qiáng)了我們對(duì)于葉綠體結(jié)構(gòu)功能、系統(tǒng)進(jìn)化的認(rèn)知，也為葉綠體轉(zhuǎn)化提供了重要信息。葉綠體基因組編碼植物參與光合作用和其他必要代謝過(guò)程的關(guān)鍵蛋白，以應(yīng)對(duì)干旱、鹽、熱、光等環(huán)境壓力及防御入侵病原體[53]。葉綠體轉(zhuǎn)化具有蛋白表達(dá)水平高、從多順?lè)醋觤RNA中表達(dá)多種蛋白的可行性、避免花粉傳播造成的基因污染等優(yōu)點(diǎn)[54]。葉綠體轉(zhuǎn)基因工程在增強(qiáng)藥用植物的優(yōu)良農(nóng)藝性狀或生產(chǎn)高價(jià)值的生物醫(yī)學(xué)產(chǎn)品等方面具有巨大潛力。

由于葉綠體自身的復(fù)制機(jī)制和相對(duì)獨(dú)立的進(jìn)化，現(xiàn)已廣泛應(yīng)用于藥用植物的物種鑒定和系統(tǒng)進(jìn)化研究中。見表1。物種鑒定是生物多樣性研究的重要前提和基礎(chǔ)，DNA條形碼不依賴于傳統(tǒng)的形態(tài)學(xué)鑒定經(jīng)驗(yàn)，此概念一經(jīng)提出便受到國(guó)內(nèi)外學(xué)者的廣泛關(guān)注和逐步認(rèn)可，但目前通用的DNA條形碼不足以準(zhǔn)確鑒別所有的藥用植物物種。隨著新一代測(cè)序技術(shù)和生物信息技術(shù)的進(jìn)步，葉綠體全基因組作為超級(jí)條形碼或基于葉綠體測(cè)序分析篩選候選DNA條形碼逐步發(fā)展起來(lái)，植物葉綠體基因組的數(shù)量也在逐年增加。就藥用植物而言，目前已發(fā)表的葉綠體基因組序列的數(shù)量還存在不足。葉綠體基因組目前已顯示出了巨大的物種識(shí)別潛力，尤其是在親緣關(guān)系較近的類群之間，因此有必要對(duì)更多的藥用植物進(jìn)行葉綠體基因組測(cè)序，深入研究其鑒定與系統(tǒng)進(jìn)化問(wèn)題。

參考文獻(xiàn)

[1]Dyall S D，Brown M T，Johnson P J.Ancient invasions：from endosymbionts to organelles[J].Science，2004，304（5668）：253-257.

[2]Shinozaki K，Ohme M，Tanaka M，et al.The complete nucleotide sequence of the tobacco chloroplast genome：its gene organization and expression[J].Embo J，1986，5（9）：2043-2049.

[3]Ohyama K，F(xiàn)ukuzawa H，Kohchi T，et al.Chloroplast gene organization deduced from complete sequence of liverwort Marchantia polymorpha chloroplast DNA[J].Nature，1986，322（6079）：572-574.

[4]Wicke S，Schneeweiss G M，Depamphilis C W，et al.The evolution of the plastid chromosome in land plants：gene content，gene order，gene function[J].Plant Mol Biol，2011，76（3-5）：273-297.

[5]Green B R.Chloroplast genomes of photosynthetic eukaryotes[J].Plant J，2011，66（1）：34-44.

[6]Tonti-Filippini J，Nevill P G，Dixon K，et al.What can we do with 1000 plastid genomes？[J].Plant J，2017，90（4）：808-818.

[7]王玲，董文攀，周世良.被子植物葉綠體基因組的結(jié)構(gòu)變異研究進(jìn)展[J].西北植物學(xué)報(bào)，2012，32（6）：1282-1288.

[8]Guo H，Liu J，Luo L，et al.Complete chloroplast genome sequences of Schisandra chinensis：genome structure，comparative analysis，and phylogenetic relationship of basal angiosperms[J].Sci China Life Sci，2017，60：1286-1290.

[9]Lee H O，Joh H J，Kim K，et al.Dynamic Chloroplast Genome Rearrangement and DNA Barcoding for Three Apiaceae Species Known as the Medicinal Herb “Bang-Poong”[J].Int J Mol Sci，2019，20（9）：2196.

[10]Lei W，Ni D，Wang Y，et al.Intraspecific and heteroplasmic variations，gene losses and inversions in the chloroplast genome of Astragalus membranaceus[J].Sci Rep，2016，6：21669.

[11]Wicke S，Müller K F，de Pamphilis C W，et al.Mechanisms of functional and physical genome reduction in photosynthetic and nonphotosynthetic parasitic plants of the broomrape family[J].Plant Cell，2013，25（3711）：113373.

[12]Li Y，Zhou J，Chen X，et al.Gene losses and partial deletion of small single-copy regions of the chloroplast genomes of two hemiparasitic Taxillus species[J].Sci Rep，2017，7（1）：12834.

[13]Martín M，Sabater B.Plastid ndh genes in plant evolution[J].Plant Physiol Bioch，2010，48（8）：636-645.

[14]Millen R S，Olmstead R G，Adams K L，et al.Many parallel losses of infA from chloroplast DNA during angiosperm evolution with multiple independent transfers to the nucleus[J].Plant Cell，2001，13（3）：645-658.

[15]Doyle J J，Doyle J L，Ballenger J A，et al.The distribution and phylogenetic significance of a 50-kb chloroplast DNA inversion in the flowering plant family Leguminosae[J].Mol Phylogenet Evol，1996，5（2）：429-438.

[16]Shen X，Wu M，Liao B，et al.Complete chloroplast genome sequence and phylogenetic analysis of the medicinal plant Artemisia annua[J].Molecules，2017，22（8）：1330.

[17]Chen X，Zhou J，Cui Y，et al.Identification of Ligularia Herbs using the complete chloroplast genome as a super-barcode[J].Front Pharmacol，2018，9：695.

[18]He L，Qian J，Li X，et al.Complete chloroplast genome of medicinal plant Lonicera japonica：Genome rearrangement，intron gain and loss，and implications for phylogenetic studies[J].Molecules，2017，22（2）：249.

[19]Hollingsworth P M，Li D Z，van der Bank M，et al.Telling plant species apart with DNA：from barcodes to genomes[J].Philos T R Soc B，2016，371（1702）：20150338.

[20]Chen S，Yao H，Han J，et al.Validation of the ITS2 region as a novel DNA barcode for identifying medicinal plant species[J].PloS one，2010，5（1）：e8613.

[21]Guo M，Xu Y，Ren L，et al.A Systematic Study on DNA Barcoding of Medicinally Important Genus Epimedium L.（Berberidaceae）[J].Genes，2018，9（12）：637.

[22]羅焜，馬培，姚輝，等.基于ITS2序列鑒定川貝母及其混偽品基原植物[J].世界科學(xué)技術(shù)-中醫(yī)藥現(xiàn)代化，2012，14（1）：1153-1158.

[23]Nock C J，Waters D L E，Edwards M A，et al.Chloroplast genome sequences from total DNA for plant identification[J].Plant Biotechnol J，2011，9（3）：328-333.

[24]Kane N，Sveinsson S，Dempewolf H，et al.Ultra-barcoding in cacao（Theobroma spp.;Malvaceae）using whole chloroplast genomes and nuclear ribosomal DNA[J].Am J Bot，2012，99（2）：320-329.

[25]Li X，Yang Y，Henry R J，et al.Plant DNA barcoding：from gene to genome[J].Biol Rev，2015，90（1）：157-166.

[26]Cui Y，Chen X，Nie L，et al.Comparison and Phylogenetic Analysis of Chloroplast Genomes of Three Medicinal and Edible Amomum Species[J].Int J Mol Sci，2019，20（16）：4040.

[27]He W Y，Pan Y Z.Study on the DNA barcoding of genus Ligularia Cass.Asteraceae）[J].Plant Divers，2015，37：693-703.

[28]鐘志敏.石斛DNA條形碼鑒定及系統(tǒng)分類研究[D].廣州：廣州中醫(yī)藥大學(xué)，2018.

[29]Xia Y，Hu Z，Li X，et al.The complete chloroplast genome sequence of Chrysanthemum indicum[J].Mitochondrial DNA A，2016，27（6）：4668-4669.

[30]李雪佩，孟靜，張琳娜，等.黃草烏與展毛黃草烏葉綠體全基因組結(jié)構(gòu)的比較分析[J].中藥材，2018，41（8）：1812-1820.

[31]Kim Y，Choi H，Shin J，et al.Molecular Discrimination of Cynanchum wilfordii and Cynanchum auriculatum by InDel Markers of Chloroplast DNA[J].Molecules，2018，23（6）：1337.

[32]Kang K B，Kang S J，Kim M S，et al.Chemical and genomic diversity of six Lonicera species occurring in Korea[J].Phytochemistry，2018，155：126-135.

[33]Nguyen V B，Park H S，Lee S C，et al.Authentication markers for five major Panax species developed via comparative analysis of complete chloroplast genome sequences[J].J Agr Food Chem，2017，65（30）：6298-6306.

[34]馬雙姣.草藥分子鑒定：從條形碼到基因組[D].北京：北京協(xié)和醫(yī)學(xué)院，2017.

[35]Chen Q，Wu X，Zhang D.Phylogenetic analysis of Fritillaria cirrhosa D.Don and its closely related species based on complete chloroplast genomes[J].PeerJ，2019，7：e7480.

[36]楊俏俏，姜梅，王立強(qiáng)，等.藥食兩用藠頭葉綠體基因組解析、比較基因組學(xué)及系統(tǒng)發(fā)育研究[J].藥學(xué)學(xué)報(bào)，2019，54（1）：173-181.

[37]李曉娟，楊振艷，黃玉玲，等.藥用植物華重樓（黑藥花科）葉綠體全基因組研究[J].熱帶亞熱帶植物學(xué)報(bào)，2015，23（6）：601-613.

[38]Do H D K，Jung J，Hyun J Y，et al.The newly developed single nucleotide polymorphism（SNP）markers for a potentially medicinal plant，Crepidiastrum denticulatum（Asteraceae），inferred from complete chloroplast genome data[J].Mol Biol Rep，2019，46（3）：3287-3297.

[39]Chen X，Cui Y，Nie L，et al.Identification and Phylogenetic Analysis of the Complete Chloroplast Genomes of Three Ephedra Herbs Containing Ephedrine[J].Biomed Res Int，2019，2019.

[40]Zhou Y，Nie J，Xiao L，et al.Comparative Chloroplast Genome Analysis of Rhubarb Botanical Origins and the Development of Specific Identification Markers[J].Molecules，2018，23（11）：2811.

[41]Drouin G，Daoud H，Xia J.Relative rates of synonymous substitutions in the mitochondrial，chloroplast and nuclear genomes of seed plants[J].Mol Phylogenetics Evol，2008，49（3）：827-831.

[42]魏似婕，趙志禮，倪梁紅，等.藏藥“巴萵色?！被b定及其葉綠體全基因組序列分析[J].藥學(xué)學(xué)報(bào)，2018，53（6）：1009-1015.

[43]沈立群.唇形科三種藥用植物葉綠體全基因組及科內(nèi)的比較與進(jìn)化分析[D].杭州：浙江大學(xué)，2018.

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

[45]Cui Y，Nie L，Sun W，et al.Comparative and Phylogenetic Analyses of Ginger（Zingiber officinale）in the Family Zingiberaceae Based on the Complete Chloroplast Genome[J].Plants，2019，8（8）：283.

[46]Gao B，Yuan L，Tang T，et al.The complete chloroplast genome sequence of Alpinia oxyphylla Miq.and comparison analysis within the Zingiberaceae family[J].PloS one，2019，14（6）：e0218817.

[47]Gichira A W，Li Z，Saina J K，et al.The complete chloroplast genome sequence of an endemic monotypic genus Hagenia（Rosaceae）：Structural comparative analysis，gene content and microsatellite detection[J].PeerJ，2017，5：e2846.

[48]Gu C，Dong B，Xu L，et al.The complete chloroplast genome of Heimia myrtifolia and comparative analysis within Myrtales[J].Molecules，2018，23（4）：846.

[49]He Y，Xiao H，Deng C，et al.The complete chloroplast genome sequences of the medicinal plant Pogostemon cablin[J].Int J Mol Sci，2016，17（6）：820.

[50]He Y，Xiao H，Deng C，et al.Complete chloroplast genome sequence of Coptis chinensis Franch.and its evolutionary history[J].Biomed Res Int，2017，2017.

[51]Kim K，Nguyen V B，Dong J，et al.Evolution of the Araliaceae family inferred from complete chloroplast genomes and 45S nrDNAs of 10 Panax-related species[J].Sci Rep，2017，7（1）：4917.

[52]Fan W B，Wu Y，Yang J，et al.Comparative chloroplast genomics of Dipsacales species：insights into sequence variation，adaptive evolution，and phylogenetic relationships[J].Front Plant Sci，2018，9：689.

[53]Liu X，Li Y，Yang H，et al.Chloroplast genome of the folk medicine and vegetable plant Talinum paniculatum（Jacq.）Gaertn.：Gene organization，comparative and phylogenetic analysis[J].Molecules，2018，23（4）：857.

[54]Wang H H，Yin W B，Hu Z M.Advances in chloroplast engineering[J].J Genet Genomics，2009，36（7）：387-398.

[55]Gao X，Zhang X，Meng H，et al.Comparative chloroplast genomes of Paris Sect.Marmorata：insights into repeat regions and evolutionary implications[J].BMC genomics，2018，19（10）：878.

[56]Guo Z Y，Zhang H R，Shrestha N，et al.Complete chloroplast genome of a valuable medicinal plant，Huperzia serrata（Lycopodiaceae），and comparison with its congener[J].Appl Plant Sci，2016，4（11）：1600071.

[57]Huang J，Chen R，Li X.Comparative analysis of the complete chloroplast genome of four known Ziziphus species[J].Genes，2017，8（12）：340.

[58]Kang S H，Lee J H，Lee H O，et al.Complete chloroplast genome and 45S nrDNA sequences of the medicinal plant species Glycyrrhiza glabra and Glycyrrhiza uralensis[J].Genes Genet Syst，2017：17-00002.

[59]Kong H，Liu W，Yao G，et al.A comparison of chloroplast genome sequences in Aconitum（Ranunculaceae）：A traditional herbal medicinal genus[J].PeerJ，2017，5：e4018.

[60]Li D M，Zhao C Y，Liu X F.Complete Chloroplast Genome Sequences of Kaempferia Galanga and Kaempferia Elegans：Molecular Structures and Comparative Analysis[J].Molecules，2019，24（3）：474.

[61]Li X，Zuo Y，Zhu X，et al.Complete chloroplast genomes and comparative analysis of sequences evolution among seven Aristolochia（Aristolochiaceae）medicinal species[J].Int J Mol Sci，2019，20（5）：1045.

[62]Li B，Lin F，Huang P，et al.Complete chloroplast genome sequence of Decaisnea insignis：Genome organization，genomic resources and comparative analysis[J].Sci Rep，2017，7（1）：10073.

[63]Li P，Lu R S，Xu W Q，et al.Comparative genomics and phylogenomics of East Asian tulips（Amana，Liliaceae）[J].Front Plant Sci，2017，8：451.

[64]Li Y，Zhang Z，Yang J，et al.Complete chloroplast genome of seven Fritillaria species，variable DNA markers identification and phylogenetic relationships within the genus[J].PloS one，2018，13（3）：e0194613.

[64]Meng D，Xiaomei Z，Wenzhen K，et al.Detecting useful genetic markers and reconstructing the phylogeny of an important medicinal resource plant，Artemisia selengensis，based on chloroplast genomics[J].PloS one，2019，14（2）：e0211340.

[65]Meng J，Li X，Li H，et al.Comparative analysis of the complete chloroplast genomes of four Aconitum medicinal species[J].Molecules，2018，23（5）：1015.

[66]Shahzadi I，Mehmood F，Ali Z，et al.Chloroplast genome sequences of Artemisia maritima and Artemisia absinthium：Comparative analyses，mutational hotspots in genus Artemisia and phylogeny in family Asteraceae[J].Genomics，2019，112（2）：1454-1463.

[67]Zhao Z，Wang X，Yu Y，et al.Complete chloroplast genome sequences of Dioscorea：Characterization，genomic resources，and phylogenetic analyses[J].PeerJ，2018，6：e6032.

[68]Wang X，Zhou T，Bai G，et al.Complete chloroplast genome sequence of Fagopyrum dibotrys：genome features，comparative analysis and phylogenetic relationships[J].Sci Rep，2018，8（1）：12379.

[69]Park I，Yang S，Choi G，et al.The complete chloroplast genome sequences of Aconitum pseudolaeve and Aconitum longecassidatum，and development of molecular markers for distinguishing species in the Aconitum subgenus Lycoctonum[J].Molecules，2017，22（11）：2012.

[70]Song Y，Wang S，Ding Y，et al.Chloroplast genomic resource of Paris for species discrimination[J].Sci Rep，2017，7（1）：3427.

（2020-02-10收稿?責(zé)任編輯：徐穎）