摘要：為了解貴州海桐（Pittosporum kweichowense）葉綠體基因組基本特征及系統(tǒng)發(fā)育關(guān)系，本研究基于Illumina高通量測(cè)序技術(shù)，采用生物信息學(xué)方法對(duì)序列進(jìn)行組裝，獲得貴州海桐完整葉綠體基因組，并對(duì)其功能特征、密碼子偏好性、簡(jiǎn)單重復(fù)序列、基因組比對(duì)信息及系統(tǒng)發(fā)育關(guān)系等進(jìn)行分析。結(jié)果表明，（1）貴州海桐葉綠體基因組呈典型的四分體結(jié)構(gòu)，總長(zhǎng)度為153 582 bp，其中大單拷貝區(qū)（LSC）為84 944 bp，小單拷貝區(qū)（SSC）為18 740 bp，此外還有2個(gè)長(zhǎng)度為24 949 bp的反向重復(fù)序列區(qū)（IR），總GC含量為38.94%。（2）注釋到130個(gè)基因，其中包括85個(gè)蛋白質(zhì)編碼基因、37個(gè)tRNA基因和8個(gè)rRNA基因。（3）共檢測(cè)到25 910個(gè)密碼子，其中編碼亮氨酸（Leu）的密碼子數(shù)量最多，密碼子第3位堿基有較高的A/U偏好性。（4）通過(guò)微衛(wèi)星分析鑒定到44個(gè)簡(jiǎn)單重復(fù)序列（SSR）位點(diǎn)，且明顯偏好使用A/T堿基。（5）海桐花屬植物的葉綠體基因組IR/SC邊界相對(duì)保守，變異主要存在于LSC。（6）系統(tǒng)發(fā)育分析結(jié)果顯示，貴州海桐、昆明海桐（P. kunmingense）的親緣關(guān)系最近。

關(guān)鍵詞：貴州海桐；葉綠體基因組；結(jié)構(gòu)特征；密碼子偏好性；系統(tǒng)發(fā)育分析

中圖分類(lèi)號(hào)：S718.46" 文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1002-1302（2025）03-0035-09

滕" 堯，楊加文，何選澤，等. 貴州海桐葉綠體基因組特征及系統(tǒng)發(fā)育分析［J］. 江蘇農(nóng)業(yè)科學(xué)，2025，53（3）：35-44.

doi：10.15889/j.issn.1002-1302.2025.03.005

收稿日期：2024-08-28

基金項(xiàng)目：貴州省科技計(jì)劃（編號(hào)：黔科合基礎(chǔ)-ZK［2024］一般625）；貴州省林業(yè)科研項(xiàng)目（編號(hào)：黔林科合［2022］01號(hào)）。

作者簡(jiǎn)介：滕" 堯（1991—），男，貴州羅甸人，碩士，助理研究員，主要研究方向?yàn)橹参镔Y源栽培利用。E-mail：574177089@qq.com。

通信作者：彭" 熙，碩士，副研究員，主要研究方向?yàn)橥寥琅c植物互作。E-mail：pengxi7979@163.com。

貴州海桐（Pittosporum kweichowense）是海桐花科（Pittosporaceae）海桐花屬（Pittosporum）常綠灌木，分布于貴州、四川和湖南部分地區(qū)，該科在全世界共有9屬約360種，我國(guó)有1屬約50種，主要分布于陜西、湖南、四川、貴州、云南等省份［1］。海桐花屬植物的根、樹(shù)皮、葉和種子均可入藥，種子還可用來(lái)榨油，是工業(yè)用油脂原料［1-3］?；诤Ｍ┗▽僦参镙^高的醫(yī)藥、工業(yè)價(jià)值，國(guó)內(nèi)外對(duì)該屬的研究方向主要集中在藥理作用（如抗結(jié)核、抗菌、抗癌、抗炎癥測(cè)試［4-5］，對(duì)新冠病毒蛋白酶的抑制效果［6］）和次生代謝物種類(lèi)及生物活性［2，7］等方面，而對(duì)該科（屬）植物系統(tǒng)進(jìn)化及資源分類(lèi)鑒定的報(bào)道較少。因此，開(kāi)展更多海桐花屬植物葉綠體基因組圖譜構(gòu)建和比較分析，既能豐富該屬植物細(xì)胞器基因組信息庫(kù)，又能為該屬的種間關(guān)系、種質(zhì)資源利用提供參考依據(jù)。

葉綠體是植物進(jìn)行光合生理所需的細(xì)胞器［8］，它除了具有典型的母系遺傳特征［9］外，還具有結(jié)構(gòu)穩(wěn)定、序列保守、變異位點(diǎn)信息較豐富等特點(diǎn)［10］。大多數(shù)被子植物葉綠體基因組都是雙鏈環(huán)狀四分體結(jié)構(gòu)，包含1個(gè)大單拷貝區(qū)（LSC）、1個(gè)小單拷貝區(qū)（SSC）和2個(gè)反向重復(fù)區(qū)（IR）［11］，加上葉綠體基因組的測(cè)序成本不高，因此常被用于資源鑒定、分類(lèi)和系統(tǒng)進(jìn)化等研究［12］。隨著研究的發(fā)展，近年來(lái)發(fā)布的海桐花科植物葉綠體基因組有近30種，但是開(kāi)展葉綠體基因組特征分析的僅有10種［13-15］，這在一定程度上限制了海桐花科植物的分類(lèi)學(xué)研究。

本研究擬基于Illumina二代測(cè)序平臺(tái)對(duì)貴州海桐的葉綠體基因組進(jìn)行測(cè)序，采用生物信息學(xué)方法對(duì)貴州海桐的葉綠體基因組進(jìn)行組裝，并初步對(duì)其基因組結(jié)構(gòu)特征、密碼子使用模式、重復(fù)序列特征、種間差異、核苷酸多態(tài)性和系統(tǒng)進(jìn)化關(guān)系等展開(kāi)分析，以期揭示貴州海桐的遺傳特征，為今后海桐花屬植物的葉綠體基因組學(xué)研究、分子標(biāo)記開(kāi)發(fā)、種質(zhì)資源分類(lèi)及利用等提供參考。

1" 材料與方法

1.1" 試驗(yàn)材料

本試驗(yàn)所用材料為貴州海桐（P. kweichowense），樣品于2024年5月10日采自貴州省望謨縣樂(lè)旺鎮(zhèn)（106°25′02.07″E，25°20′01.43″N，海拔1 076 m），取樣后將葉片保存于液氮中并立即帶回貴州科學(xué)院貴州省山地資源研究所實(shí)驗(yàn)室進(jìn)行預(yù)處理，將標(biāo)本保存于貴州科學(xué)院貴州省植物園標(biāo)本館內(nèi)，標(biāo)本號(hào)Y240505。貴州海桐的葉綠體基因組相關(guān)數(shù)據(jù)已上傳至CNGBdb（https：//db.cngb.org/），項(xiàng)目身份標(biāo)志號(hào)（ID）為CNP0005941，樣本ID為CNS1153709，測(cè)序ID為CNX1085786。

在本研究中，14個(gè)用于校正和比對(duì)分析的海桐花科植物的完整葉綠體基因組數(shù)據(jù)下載自美國(guó)國(guó)家生物技術(shù)信息中心（NCBI），包括羊脆木（P.kerrii，登錄號(hào)：NC_046847）、滇越海桐（P.merrillianum，登錄號(hào)：NC_068010）、縫線(xiàn)海桐（P.perryanum，登錄號(hào)：NC_068014）、昆明海桐（P.kunmingense，登錄號(hào)：NC_068009）、短萼海桐（P.brevicalyx，登錄號(hào)：NC_058267）、皺葉海桐（P.crispulum，登錄號(hào)：NC_068008）、海金子（P.illicioides，登錄號(hào)：NC_088055）、小果海桐（P.parvicapsulare，登錄號(hào)：NC_068013）、棱果海桐（P.trigonocarpum，登錄號(hào)：NC_068015）、海桐（P.tobira，登錄號(hào)：NC_057165）、滇藏海桐（P.napaulense，登錄號(hào)：NC_068011）、圓錐海桐（P.paniculiferum，登錄號(hào)：NC_068012）、崖花子（P.truncatum，登錄號(hào)：NC_068016）和P.serpentinum（登錄號(hào)：NC_084331）。

1.2" DNA提取、建庫(kù)及測(cè)序

用十六烷基三甲基溴化銨（CTAB）法［16］提取貴州海桐葉片的全基因組DNA，用DNA提取試劑盒（諾貝萊）獲得高質(zhì)量基因組DNA，用凝膠電泳初步檢測(cè)DNA的完整性，用NanoDrop微量分光光度計(jì)檢測(cè)DNA的純度，用安捷倫2100檢測(cè)DNA的完整性，RIN（RNA完整值）gt;7.0表明樣品完整性較好，可以用來(lái)進(jìn)行文庫(kù)構(gòu)建。隨后用 NEBNext UltraTMⅡDNA Library Prep Kit經(jīng)隨機(jī)擴(kuò)增、片段分選，獲得長(zhǎng)度為350～400 bp的隨機(jī)片段用于上機(jī)測(cè)序。在NGS重測(cè)序中，采用華大基因的DNBSEQ-T7平臺(tái)進(jìn)行雙端測(cè)序，下機(jī)數(shù)據(jù)經(jīng)過(guò)質(zhì)控、過(guò)濾，最終獲得150 bp的雙端測(cè)序讀段。

1.3" 葉綠體基因組組裝及注釋

用MAFFT軟件［17］對(duì)數(shù)據(jù)進(jìn)行校正后，用GetOrganelle軟件［18］組裝貴州海桐的葉綠體基因組，用在線(xiàn)工具CPGAVAS2（http：//47.96.249.172：16019/analyzer/annotate）［19］進(jìn)行注釋分析，用在線(xiàn)工具Chloroplot（https：//irscope.shinyapps.io/Chloroplot/）［20］繪制貴州海桐的葉綠體基因組圈圖。

1.4" 密碼子偏好性分析

用CodonW 1.4.2軟件［21］分析貴州海桐葉綠體基因組同義密碼子相對(duì)使用度（relative synonymous codon usage，RSCU）。

1.5" 簡(jiǎn)單重復(fù)序列（SSR）分析

通過(guò)在線(xiàn)工具M(jìn)ISA（https：//webblast.ipkgatersleben.de/misa/）［22］分析貴州海桐葉綠體基因組的SSR，單核苷酸重復(fù)數(shù)量≥10個(gè)，二核苷酸重復(fù)數(shù)量≥6個(gè)，三核苷酸、四核苷酸、五核苷酸重復(fù)數(shù)量≥5個(gè)，2個(gè)SSR的最小距離設(shè)為100 bp。

1.6" 基因組比較分析

基于MAFFT數(shù)據(jù)比對(duì)后的系統(tǒng)發(fā)育分析結(jié)果，通過(guò)在線(xiàn)工具IRSCOPE（https：//irscope.shinyapps.io/irapp/）［23］進(jìn)行葉綠體基因組的比較分析。

1.7" 屬內(nèi)核苷酸的多態(tài)性分析

用DnaSP軟件［24］分析海桐花屬植物核苷酸的多態(tài)性（π），根據(jù)π值篩選出變異較高的位點(diǎn)。

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

將MAFFT軟件的比對(duì)結(jié)果導(dǎo)入RAxML-ng-v1.2.0軟件［25］，基于GTRCAT模型構(gòu)建最大似然（ML）系統(tǒng)發(fā)育樹(shù)，Bootstrap重復(fù)1 000次。

2" 結(jié)果與分析

2.1" 葉綠體基因組結(jié)構(gòu)的基本特征

貴州海桐葉綠體基因組為典型的環(huán)狀四分體結(jié)構(gòu)，長(zhǎng)度為153 582 bp，包含長(zhǎng)度為84 944 bp的大單拷貝區(qū)（LSC）、長(zhǎng)度為18 740 bp的小單拷貝區(qū)（SSC）及長(zhǎng)度為24 949 bp的2個(gè)反向重復(fù)序列區(qū)（IR）（圖1）。GC總含量為38.94%，IR的GC含量高于LSC、SSC（表1）。

注釋結(jié)果顯示，貴州海桐葉綠體基因組共編碼130個(gè)基因，包括85個(gè)蛋白質(zhì)編碼基因、8個(gè)rRNA基因和37個(gè)tRNA基因。按功能可將貴州海桐葉綠體基因組分為4類(lèi)，其中參與光合作用的基因45個(gè)，參與自我復(fù)制的基因74個(gè)，其他基因6個(gè)及功能未知的基因5個(gè)，其中有18個(gè)基因具有內(nèi)含子，ndhA、ndhB、petB、petD、atpF、rpl16、rpl2、rps16、rpoC1、trnA-UGC、trnG-UCC、trnI-GAU、trnK-UUU、trnL-UAA、trnV-UAC各含1個(gè)內(nèi)含子，rps12、clpP、ycf3有2個(gè)內(nèi)含子（表2）。

2.2" 密碼子偏好性分析結(jié)果

測(cè)序結(jié)果顯示，從貴州海桐葉綠體基因組中共檢測(cè)出25 910個(gè)密碼子， 其中編碼亮氨酸（Leu）的密碼子數(shù)量最多，為2 722個(gè)，占總密碼子數(shù)的10.51%，終止密碼子（TER）數(shù)量最少，為85個(gè)，占總密碼子數(shù)的0.33%，編碼亮氨酸（Leu）、絲氨酸（Ser）和精氨酸（Arg）的密碼子種類(lèi)最多，均為6種（表3、圖2）。密碼子具有較明顯的偏好性，RSCU=1的氨基酸有2個(gè)，分別為甲硫氨酸（Met）和色氨酸（Trp），RSCU＞1的高頻密碼子有30個(gè)，其中沒(méi)有密碼子第3位堿基以C結(jié)尾，只有1個(gè)以G結(jié)尾的密碼子（Leu的UUG），剩余29個(gè)密碼子均以A/U結(jié)尾。

2.3" 簡(jiǎn)單重復(fù)序列分析結(jié)果

微衛(wèi)星位點(diǎn)分析結(jié)果（表4）顯示，從貴州海桐葉綠體基因組中共識(shí)別出44個(gè)SSR位點(diǎn)，其中單核苷酸43個(gè)，占位點(diǎn)總數(shù)的97.73%，雙核苷酸1個(gè)，三核苷酸及以上的SSR位點(diǎn)未被識(shí)別出，SSR位點(diǎn)偏好使用A/T堿基。SSR位點(diǎn)主要分布在LSC，占比達(dá)74%，其余分布在SSC、IR（圖3）。

2.4" 葉綠體基因組比較分析結(jié)果

對(duì)貴州海桐（P.kweichowense）、羊脆木（P.kerrii）、短萼海桐（P.brevicalyx）、海金子（P.illicioides）、海桐（P.tobira）及P. serpentinum等6種海桐花屬植物葉綠體基因組IR/SC邊界分析的結(jié)果（圖4）顯示，6種海桐花屬葉綠體基因組長(zhǎng)度無(wú)明顯差異， 即無(wú)明顯的邊界收縮或擴(kuò)張，但邊界基因的位置和類(lèi)型有差別。短萼海桐邊在JLB邊界上的基因?yàn)閞pl22、rpl2，其余5種均為rps19、rpl2，且rps19基因均跨越JLB邊界，P.serpentinum的rps19基因位置與其余4種略有不同；trnN、ndhF基因分布在JSB邊界；ycf1、trnN基因分布在JSA邊界，大小相同且ycf1均橫跨該區(qū)域；rpl2、trnH基因分布在JLA邊界，大小相同且無(wú)基因跨越該區(qū)域；P. serpentinum 在JSB、JSA及JLA邊界的基因位置與其余5種稍有差異。

以貴州海桐葉綠體基因組為參照，對(duì)海桐花屬其他5種植物葉綠體基因組進(jìn)行比較。mVISTA比對(duì)結(jié)果（圖5）表明，海桐花屬植物葉綠體基因組的基因結(jié)構(gòu)整體上十分接近，P.serpentinum在51～54 kb（LSC）處存在變異，變異基因有trnV-UAC、trnM-CAU，海金子在94～99 kb（IRA）處也存在明顯變異，變異基因包括ndhB、rps7和rps12，其余位置則無(wú)明顯差異。

2.5" 屬內(nèi)核苷酸多態(tài)性分析結(jié)果

核苷酸多態(tài)性π峰值圖（圖6）顯示，海桐花屬植物葉綠體基因組變異較高的位點(diǎn)主要存在于LSC，其次是SSC，IR則相對(duì)保守，其中變異最高的4個(gè)位點(diǎn)（π＞0.008）均位于LSC，分別為rps4-trnT（47 446～48 084 bp）、trnT-trnL-exon1（47 885～48 485 bp）、ndhC-trnV-exon2（52 318～52 917 bp）和rpl33-rps18（69 042～69 645 bp）。

2.6" 系統(tǒng)進(jìn)化分析結(jié)果

ML系統(tǒng)發(fā)育樹(shù)結(jié)果（圖7）表明，羊脆木（P.kerrii）和滇越海桐（P.merrillianum）聚為一支，自展值為72，二者又與縫線(xiàn)海桐（P.perryanum）聚為一支，自展值為65，再與昆明海桐（P.kunmingense）聚為一支，自展值為53，該支最終與貴州海桐（P.kweichowense）聚為一支，自展值為100。上述結(jié)果說(shuō)明，貴州海桐與羊脆木、滇越海桐、縫線(xiàn)海桐和昆明海桐的親緣關(guān)系較近，其中與昆明海桐的親緣關(guān)系最近，而與小果海桐（P.parvicapsulare）、海桐（P.tobira）和崖花子（P.truncatum）等同屬植物的親緣關(guān)系較遠(yuǎn)。

3" 討論與結(jié)論

貴州海桐葉綠體基因組總長(zhǎng)度為153 582 bp，包含1個(gè)LSC（84 944 bp）、1個(gè)SSC（18 740 bp）和2個(gè)IR（24 949 bp）， 為典型的四分體結(jié)構(gòu)， 與已報(bào)道的大多數(shù)被子植物葉綠體基因組結(jié)構(gòu)特征一致［26］。G+C堿基對(duì)含量在植物細(xì)胞中比較穩(wěn)定，不容易受外界因素影響，常用于物種、屬的分類(lèi)［27-28］，貴州海桐葉綠體基因組總GC含量為38.94%，與海桐花屬已完成葉綠體基因組測(cè)序的物種間差異很小［13-15］。有趣的是，貴州海桐葉綠體基因組中共有85個(gè)蛋白質(zhì)編碼基因，與Zhang等已完成測(cè)序和注釋得出的海桐花近緣屬種均為86個(gè)蛋白質(zhì)編碼基因的結(jié)果不同，而tRNA和rRNA的基因數(shù)量則沒(méi)有差異［15］?？紤]到Zhang的試驗(yàn)材料中縫線(xiàn)海桐（P. perryanum）和小果海桐（P. parvicapsulare）樣品也采自貴州，蛋白質(zhì)編碼基因存在微小差異的原因應(yīng)該不是貴州喀斯特氣候環(huán)境，該原因需要進(jìn)一步研究。從貴州海桐葉綠體基因組中共鑒定到5個(gè)未知功能基因，均為ycf基因家族， 該基因家族在不同植物中的種類(lèi)、數(shù)量和功能均有差異［29］，關(guān)于海桐屬植物ycf基因家族具體功能的研究目前尚未見(jiàn)報(bào)道。

密碼子偏好性隨生物的突變、自然選擇和基因漂變而進(jìn)化，密碼子使用偏好是決定基因表達(dá)、細(xì)胞功能的關(guān)鍵因素，它通過(guò)影響RNA加工、蛋白質(zhì)翻譯和蛋白質(zhì)折疊等多種過(guò)程來(lái)實(shí)現(xiàn)， 密碼子偏好性分析主要應(yīng)用在轉(zhuǎn)基因設(shè)計(jì)中，通過(guò)優(yōu)化密碼子來(lái)增加基因的表達(dá)水平，以開(kāi)發(fā)轉(zhuǎn)基因作物［30］。貴州海桐葉綠體基因組中RSCU＞1的密碼子絕大部分以A/U堿基結(jié)尾，這與大多數(shù)雙子葉植物密碼子偏好使用A/T（U）堿基結(jié)尾的結(jié)論［31］一致。研究貴州海桐密碼子偏好性可以了解其遺傳信息傳遞規(guī)律，在后續(xù)生物技術(shù)的開(kāi)發(fā)和應(yīng)用中可幫助提高外源基因在宿主細(xì)胞中的表達(dá)效率。

SSR分子標(biāo)記具有隨機(jī)分布、多態(tài)性強(qiáng)、顯性遺傳等特性，對(duì)SSR分子標(biāo)記的比較研究有助于發(fā)掘關(guān)鍵性狀相關(guān)QTL，進(jìn)而確定候選基因［32］，是植物功能基因組學(xué)中最具信息量和用途較多的遺傳標(biāo)記之一，在分子生物學(xué)中有廣泛應(yīng)用［33］。本研究從貴州海桐葉綠體基因組中共檢測(cè)出44個(gè)SSR位點(diǎn)，包含單核苷酸重復(fù)、二核苷酸重復(fù)，其中單核苷酸重復(fù)序列A/T的重復(fù)次數(shù)最多，這與目前已報(bào)到的多種植物葉綠體基因組SSR特征［10，34］基本一致。本研究首次完成了貴州海桐SSR的位點(diǎn)研究，可為海桐花屬植物鑒定、遺傳圖譜構(gòu)建、QTL定位及分子標(biāo)記開(kāi)發(fā)等提供參考信息。

植物在長(zhǎng)期適應(yīng)環(huán)境的過(guò)程中，其葉綠體基因組邊界的收縮/擴(kuò)張時(shí)有發(fā)生，這會(huì)導(dǎo)致序列改變、基因減少或增加等，是植物葉綠體基因組進(jìn)化的主要機(jī)制［35-36］。序列比對(duì)結(jié)果顯示，貴州海桐與5個(gè)近緣物種的葉綠體基因組LSC/IR和SSC/IR邊界分布的基因類(lèi)型、結(jié)構(gòu)及排列差異性較小，表明本研究對(duì)貴州海桐葉綠體基因組的組裝完成度較高，邊界的收縮或擴(kuò)張僅有幾十至幾百bp不等，證明海桐花屬植物的邊界保守性較高。但Zhang等研究發(fā)現(xiàn)，海桐花科植物葉綠體基因組邊界區(qū)擴(kuò)張或收縮差異顯著［15］，其一原因可能是筆者根據(jù)現(xiàn)有的海桐花屬植物葉綠體基因組數(shù)據(jù)庫(kù)對(duì)注釋結(jié)果進(jìn)行了矯正以降低分析誤差；其二原因可能是Zhang的材料及比對(duì)數(shù)據(jù)來(lái)源與本研究的差異較大， 可能由于

不同地域的海桐花科植物為了適應(yīng)氣候環(huán)境產(chǎn)生不同程度的變異，導(dǎo)致邊界擴(kuò)張或收縮，因而在比對(duì)分析時(shí)差異顯著，這說(shuō)明海桐花科植物葉綠體基因組中基因及其排序保守性和邊界擴(kuò)張或收縮多樣性并存。

核苷酸多態(tài)性（π值）用于量化種群內(nèi)部或不同種群間的遺傳多樣性，π值越大，說(shuō)明核苷酸的多樣性越高，也可以用來(lái)推演進(jìn)化關(guān)系［37］。核苷酸多態(tài)性分析結(jié)果顯示，海桐花屬植物葉綠體基因組變異較高的位點(diǎn)主要存在于LSC、SSC，其中LSC變異位點(diǎn)數(shù)量最多且π值更大，IR的變異程度最低，植物葉綠體基因組序列的變異程度可能與A/T密碼子使用偏好性呈正相關(guān)，即GC含量低的LSC、SSC的變異程度較高，GC含量高的IR的變異程度較低。變異度較高的區(qū)域和位點(diǎn)可為種群遺傳學(xué)提供潛在的分子標(biāo)記，用于構(gòu)建DNA指紋圖譜、種間關(guān)系研究及基因克隆等［35］。

了解物種之間的親緣關(guān)系是生物學(xué)研究的基礎(chǔ)，一個(gè)精確的系統(tǒng)進(jìn)化樹(shù)有助于理解物種進(jìn)化過(guò)程中的重大轉(zhuǎn)變，是推測(cè)新基因起源、檢測(cè)分子適應(yīng)性、理解形態(tài)特征進(jìn)化和重建最近分化物種等的關(guān)鍵［38］。本研究中的系統(tǒng)進(jìn)化分析結(jié)果與Zhang等對(duì)海桐屬系統(tǒng)發(fā)育研究結(jié)果相似，如短萼海桐和皺葉海桐聚為一支，棱果海桐和滇藏海桐聚為一支等，但個(gè)別物種的親緣關(guān)系遠(yuǎn)近有一定差異［15］。貴州海桐所在分支的其他物種自展值不高，這可能是因?yàn)楸狙芯克玫奈锓N數(shù)量較小，且部分物種的系統(tǒng)發(fā)育位置仍存在爭(zhēng)議，未來(lái)需完善物種取樣覆蓋度，并結(jié)合更多的物種信息來(lái)探索海桐科植物的進(jìn)化和發(fā)展，明確屬內(nèi)植物分類(lèi)地位，為后續(xù)系統(tǒng)發(fā)育、資源鑒定和進(jìn)化分析提供依據(jù)。

參考文獻(xiàn)：

［1］中國(guó)科學(xué)院中國(guó)植物志編輯委員會(huì). 中國(guó)植物志：第二十七卷［M］. 北京：科學(xué)出版社，1979：1979-1980.

［2］樊鈺虎，劉" 江，王澤秀. 海桐花屬植物化學(xué)成分及藥理活性研究進(jìn)展［J］. 中草藥，2013，9（42）： 1842-1851.

［3］國(guó)家中醫(yī)藥管理局《中華本草》編委會(huì). 中華本草［M］. 上海：上海科學(xué)技術(shù)出版社，1999：2522-2523.

［4］吳虢東，錢(qián)中清，戴" 軍，等. 短萼海桐抗結(jié)核實(shí)驗(yàn)研究［J］. 中國(guó)防癆雜志，2007，29（1）：41-43.

［5］Madikizela B，McGaw L J. Pittosporum viridiflorum Sims （Pittosporaceae）：a review on a useful medicinal plant native to South Africa and tropical Africa［J］. Journal of Ethnopharmacology，2017，205：217-230.

［6］Thodi R C，Ibrahim J M，Surendran V A，et al. Rutaretin1′-（6″-sinapoylglucoside）：promising inhibitor of COVID 19 mpro catalytic dyad from the leaves of Pittosporum dasycaulon Miq （Pittosporaceae）［J］. Journal of Biomolecular Structure amp; Dynamics，2022，40（23）：12557-12573.

［7］夏建開(kāi). 尖萼海桐次生代謝產(chǎn)物及其生物活性研究［D］. 楊凌：西北農(nóng)林科技大學(xué)，2022：2-13.

［8］Hu Y Y，Sun Y Q，Zhu Q H，et al. Poaceae chloroplast genome sequencing：great leap forward in recent ten years［J］. Current Genomics，2023，23（6）：369-384.

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

［10］龔意輝，謝雪陽(yáng)，魏媛媛，等. 金冠8-18油桃葉綠體全基因組序列特征及其系統(tǒng)發(fā)育［J］. 江蘇農(nóng)業(yè)科學(xué)，2023，51（19）：30-36.

［11］Daniell H，Lin C S，Yu M，et al. Chloroplast genomes：diversity，evolution，and applications in genetic engineering［J］. Genome Biology，2016，17（1）：134.

［12］宋" 蕓，賈孟君，曹亞萍，等. 連翹葉綠體基因組特征分析［J］. 園藝學(xué)報(bào)，2022，49（1）：187-199.

［13］Wang Y，Zhao C，Li Y Q. The complete chloroplast genome sequence of Pittosporum kerrii：the first Pittosporaceae plastome［J］. Mitochondrial DNA Part B（Resources），2019，4（2）：3959-3960.

［14］Xiao M K，Xiong X K，Shen S B，et al. The complete chloroplast genome of Pittosporum brevicalyx［J］. Mitochondrial DNA Part B（Resources），2021，6（7）：1919-1920.

［15］Zhang S D，Ling L Z，Zhang Q H. Comparative and phylogenetic analysis of the complete chloroplast genomes of ten Pittosporum species from East Asia［J］. Functional amp; Integrative Genomics，2024，24（2）：64.

［16］陳昆松，李" 方，徐昌杰，等. 改良CTAB法用于多年生植物組織基因組DNA的大量提?。跩］. 遺傳，2004，26（4）：529-531.

［17］Katoh K，Standley D M. MAFFT multiple sequence alignment software version 7：improvements in performance and usability［J］. Molecular Biology and Evolution，2013，30（4）：772-780.

［18］Jin J J，Yu W B，Yang J B，et al. GetOrganelle：a fast and versatile toolkit for accurate de novo assembly of organelle genomes［J］. Genome Biology，2020，21（1）：241.

［19］Shi L C，Chen H M，Jiang M，et al. CPGAVAS2，an integrated plastome sequence annotator and analyzer［J］. Nucleic Acids Research，2019，47（W1）：W65-W73.

［20］Zheng S Y，Poczai P，Hyvnen J，et al. Chloroplot：an online program for the versatile plotting of organelle genomes［J］. Frontiers in Genetics，2020，11：576124.

［21］王子豪，郭佳樂(lè)，范" 琪，等. 紅邊龍血樹(shù)葉綠體基因組特征及其系統(tǒng)發(fā)育分析［J］. 生物工程學(xué)報(bào)，2023，39（7）：2926-2938.

［22］Habibi N，Uddin S，Behebehani M，et al. Data on SSR markers and SNPs filtered from transcriptome of Parvocalanus crassirostris［J］. Data in Brief，2023，50：109449.

［23］Amiryousefi A，Hyvnen J，Poczai P. IRscope：an online program to visualize the junction sites of chloroplast genomes［J］. Bioinformatics，2018，34（17）：3030-3031.

［24］張浩蓉，禿玉翔，李啟少，等. 桑屬植物葉綠體基因組特征比較及系統(tǒng)發(fā)育分析［J］. 江蘇農(nóng)業(yè)科學(xué)，2024，52（7）：34-40.

［25］Hhler D，Pfeiffer W，Ioannidis V，et al. RAxML Grove：an empirical phylogenetic tree database［J］. Bioinformatics，2022，38（6）：1741-1742.

［26］Maheswari P，Kunhikannan C，Yasodha R. Chloroplast genome analysis of angiosperms and phylogenetic relationships among Lamiaceae members with particular reference to teak （Tectona grandis L.f）［J］. Journal of Biosciences，2021，46：43.

［27］Mann S，Chen Y P P. Bacterial genomic G+C composition-eliciting environmental adaptation［J］. Genomics，2010，95（1）：7-15.

［28］Meier-Kolthoff J P，Klenk H P，Gker M. Taxonomic use of DNA G+C content and DNA-DNA hybridization in the genomic age［J］. International Journal of Systematic and Evolutionary Microbiology，2014，64（Pt2）：352-356.

［29］高娜娜，趙志禮，倪梁紅. 植物葉綠體ycf15基因應(yīng)用于藥用植物鑒定的前景展望［J］. 中草藥，2017，48（15）：3210-3217.

［30］Parvathy S T，Udayasuriyan V，Bhadana V. Codon usage bias［J］. Molecular Biology Reports，2022，49（1）：539-565.

［31］龔意輝，唐詩(shī)眙，周桂花，等. 黃荊坪竹根椒葉綠體基因組密碼子使用偏好性及影響因素分析［J］. 江蘇農(nóng)業(yè)科學(xué)，2023，51（20）：28-34.

［32］Kumar M，Choi J Y，Kumari N，et al. Molecular breeding in Brassica for salt tolerance：importance of microsatellite （SSR） markers for molecular breeding in Brassica［J］. Frontiers in Plant Science，2015，6：688.

［33］Taheri S，Abdullah T L，Yusop M R，et al. Mining and development of novel SSR markers using next generation sequencing （NGS） data in plants［J］. Molecules，2018，23（2）：399.

［34］趙鵬宇，趙" 威，侯智揚(yáng)，等. 黃背草與2種菅屬植物葉綠體基因組特征比較及系統(tǒng)發(fā)育分析［J］. 中草藥，2023，54（10）：3261-

3272.

［35］馬錄花，寧佳奇，王永杰，等. 桃兒七葉綠體比較基因組學(xué)分析［J］. 生物工程學(xué)報(bào)，2022，38（10）：3695-3712.

［36］Danecek P，Auton A，Abecasis G，et al. The variant call format and VCFtools［J］. Bioinformatics，2011，27（15）：2156-2158.

［37］Yang T G，Wu Z H，Tie J，et al. A comprehensive analysis of chloroplast genome provides new insights into the evolution of the genus Chrysosplenium［J］. International Journal of Molecular Sciences，2023，24（19）：14735.

［38］Kapli P，Yang Z H，Telford M J. Phylogenetic tree building in the genomic age［J］. Nature Reviews Genetics，2020，21（7）：428-444.