叶绿体是植物细胞中一种半自主复制的细胞器，是植物进行光合作用的场所，能将光能转化为化学能进行储存，因此，叶绿体对植物乃至人类的生存至关重要。目前，多数观点认为叶绿体的起源来源于真核细胞早期捕获蓝藻或光合细菌共生形成(内共生假说^[1])。叶绿体基因组是叶绿体具有的独立遗传物质，最早完成叶绿体基因组测序的物种为烟草(Nicotiana tabacum)^[2]与地钱(Marchantia polymorpha)^[3]。在早期，叶绿体基因组序列需首先提取叶绿体DNA，并随机打断或酶切后才能通过克隆和桑格测序技术获得，获取过程费时费力^[4]，随后，通过设计大量保守引物及长PCR方法^[5]使叶绿体基因组学研究进入了快车道。

随着高通量测序技术的快速发展，其测序数据通量和准确率都得到了极大提升，而且测序费用较为低廉，使该技术已普遍应用于当前的生物学研究中。基因组浅层测序是一种二代全基因组测序方法，可高效便捷地获取基因组中高重复的基因或区域(如重复元件、核DNA的核糖体内转录间隔区等)^[5]。叶绿体基因组在光合组织中的含量较高，约为细胞内总DNA含量的10%~20%^[6]，因而基因组浅层测序技术可获得高丰度的叶绿体基因组信息。此外，不断发展的组装及注释技术^[7–9]也使近年来叶绿体基因组数据呈爆发式增长。

植物叶绿体基因组在园艺植物研究方面具有重要价值^[6]，叶绿体基因组可用于园艺植物的起源、演化历史及遗传多样性方面的研究。叶绿体基因组通常为单亲遗传(且绝大多数为母系遗传)，其序列结构及组成相对保守，核苷酸变异虽低于核基因但高于植物线粒体基因。叶绿体基因组的有效群体大小是雌雄异株植物的1/4或雌雄同株植物的1/2，因此叶绿体单倍型的固定时间及溯祖时间都较短，有利于我们发现其遗传变异并追溯其演化历史^[6]。因此，叶绿体基因组在植物起源、演化历史、遗传多样性、杂交方向等领域具有重要研究价值。叶绿体基因组有利于园艺植物的分子定向育种和杂交选育研究。叶绿素含量是植物叶片的重要生理指标, 但杂交育种或基因工程会导致核质不亲和或叶绿素合成不正常而形成花叶或白化，增加园艺植物的观赏价值^{[6, 10–11]}。叶绿体的合成或质体形成过程中的自发突变是导致一些园艺嵌合体植物[如天竺葵(Pelargonium×hybridum)、牵牛花(Ipomoea nil)、紫茉莉(Mirabilis jalapa)、大丽菊(Dahlia pinnata)等]花斑现象形成的原因^[12]。最近，有研究比较了20种嵌合体植物的叶组织(绿叶组织和白化组织)的叶绿体基因组，发现了8个叶绿体编码的RNA聚合酶或光系统有关的基因中存在14个与白化现象产生相关的点突变^[11]。

蔷薇科(Rosaceae)绣线菊属(Spiraea)是一类具有重要园艺观赏价值的植物类群。根据最新的蔷薇科分类系统，绣线菊属归属于桃亚科(Amygda- loideae)绣线菊族(Spiraeeae)^[13–14]，为北温带分布的大属，全球包含约90种(http://www.plantsoftheworldonline.org)，间断分布于欧、亚和北美洲地区。我国是绣线菊属的现代分布和分化中心，多达70种(其中47种为我国特有^[15])。因此，我国拥有丰富的绣线菊属的种质资源。绣线菊属植物多为伞形、伞形总状、复伞房或圆锥花序，多数种类具有美丽的花朵、细致的叶片，生性强健，具有耐寒、耐旱、繁殖简便等优点，可应用于花篱、花径、花坛的栽植及园林主景的构造^[16]，是具有重要经济价值的观赏类园艺植物。然而，我国绣线菊属植物的园林应用仅限于几种常见种类，与我国丰富的种质资源十分不符。基于菱叶绣线菊(S.×vanhouttei)培育出了多种彩叶园艺品种，其中‘粉霜’彩叶绣线菊(‘Pink Ice’)和‘黄金喷泉’菱叶绣线菊(‘Gold Fountain’)是两个性状优良的菱叶绣线菊品种，表现出较好的园艺价值。

近年来，绣线菊属植物的叶绿体基因组研究得到了一些发展。目前，已公开8种绣线菊属植物的叶绿体基因组数据^{[14, 17–23]}，绣线菊属叶绿体基因组具有典型的四分体结构，长度为156 226~155 981 kb，编码基因的数量变化较大，菱叶绣线菊有129个编码基因^[23]，而渐尖叶粉花绣线菊(S. japonica var. acumi- nata)有125个编码基因^[22]，这可能是注释不完整导致的，因此需进行更详尽的注释及比较基因组学研究。此外，‘粉霜’和‘黄金喷泉’的叶色形成机制及与绣线菊属其他植物的演化关系尚未开展研究。鉴于叶绿体基因组在叶片色素形成及植物亲缘关系及进化方面的重要价值，本研究基于二代测序的浅层测序技术，对‘粉霜’和‘黄金喷泉’的叶绿体基因组进行组装、注释、绘制其叶绿体基因组图谱；结合网上已有的绣线菊属植物的叶绿体全基因组开展比较基因组学研究，探讨‘粉霜’和‘黄金喷泉’的叶绿体基因组的组成及高变异位点；以及这两个品种与绣线菊属植物的亲缘关系，为深入解析菱叶绣线菊及栽培品种的遗传基础和开发绣线菊属的园艺资源奠定基础。

‘粉霜’彩叶绣线菊(凭证标本号：Yan Hai-Fei Y1940)和‘黄金喷泉’菱叶绣线菊(凭证标本号：Yan Hai-Fei Y1948)植株购自浙江桐乡四季花园家庭农场，并栽种于杭州植物园。凭证标本保存在中国科学院华南植物园标本馆。采集新鲜健康叶片保存在透气的纸袋中，并放入带有变色硅胶的封口袋中干燥保藏。在实验室采用CTAB法^[24]从干燥叶片中提取总DNA，并在深圳华大基因股份有限公司(武汉)的MGISEQ-2000基因组测序平台上开展基因组浅层测序，获得原始读长数据(约2 Gbp)。

原始读长数据在去除接头和低质量数据后，使用软件GetOrganelle v1.7.3.5^[8]进行叶绿体基因组的从头组装。将组装好的序列使用在线注释网站Geseq (https://chlorobox.mpimp-golm.mpg.de/geseq.html ^[25])进行基因注释，并在Geneious Prime 2019软件^[26]中进行适当的手工调整。通过线上工具OGDRAW v1.3.1 (https://chlorobox.mpimp-golm.mpg.de/OGDraw.html)绘制叶绿体基因组序列物理图谱^[27]。

将叶绿体基因组导入Geneious Prime 2019软件, 比较叶绿体基因组总长度及大单拷贝区(LSC)、小单拷贝区(SSC)和反向重复区(IR)的长度、基因组成、GC含量等叶绿体基因组的基本信息。利用mVISTA在线软件(https://genome.lbl.gov/vista/mvista/submit.shtml)的Shuffle-LAGAN模型检测叶绿体基因组的序列异质性。利用在线软件Reputer (https://bibiserv.cebitec.uni-bielefeld.de/repuer)鉴定叶绿体基因组的重复序列类型(如正向重复、倒置重复、互补重复、回文重复)^[28]；其参数设置为：最小重复长度设置为8，最小排列值为50。

使用简单重复序列(SSR)识别软件MISA^[29]鉴定叶绿体基因组的单核苷酸、二核甘酸、三核苷酸、四核苷酸、五核苷酸和六核苷酸SSR；每种类型最小重复参数设置依次为10、5、4、3、3、3；两个SSR之间的最小距离设置为100 bp；当两个SSR之间的距离小于100 bp时识别为复合微卫星(compound SSR)。核苷酸多态性分析采用DnaSP 6.0软件^[30]，窗口大小设置为500 bp，步移大小设置为200 bp。

从NCBI的GenBank中下载所有绣线菊属物种的叶绿体基因组序列数据，以中华绣线梅(Neillia sinensis) (MT683856)和华北珍珠梅(Sorbaria kirilowii) (NC_062461)为外类群(表 1)。本研究提取出所有叶绿体基因组中共有的蛋白质编码区(Protein-coding genes, PCGs)，并采用MAFFT v7.308软件^[31]比对后形成串联序列矩阵。利用RAxML v8.2.10软件^[32]序列矩阵进行最大似然树的构建，核苷酸替代模型设置为GTRGAMMA，执行1 000次重复自举抽样以获得系统树上每个分支的支持率。

表 1(Table 1)

表 1 绣线菊属叶绿体基因组的基本特征 Table 1 Chloroplast genome feature of genus Spiraea 序号No. 物种Species 总长Total length (bp) 大拷贝区Large single-copy region (bp) 重复区Inverted repeat region (bp) 小拷贝区Small single-copy region (bp) GC /% 总基因数Total gene number 蛋白质编码基因Protein- coding gene rRNA基因rRNA gene tRNA基因tRNA gene 1 绣球绣线菊S. blumei 155 957 84 384 26 343 18 887 36.8 130 85 8 37 2 三裂绣线菊S. trilobata 155 981 84 417 26 343 18 878 36.8 130 85 8 37 3 渐尖叶粉花绣线菊S. japonica var. acuminata 153 822 82 226 26 344 18 908 36.6 129 84 8 37 4 S. insularis 158 637 86 997 26 365 18 910 36.9 130 85 8 37 5 蒙古绣线菊S. mongolica 155 949 84 375 26 340 18 894 36.7 130 85 8 37 6 毛枝绣线菊S. martini 156 226 84 378 26 492 18 864 36.8 130 85 8 37 7 珍珠绣线菊S. thunbergii 155 922 84 360 26 342 18 879 36.8 130 85 8 37 8 菱叶绣线菊S.×vanhouttei 155 957 84 384 26 340 18 893 36.8 130 85 8 37 9 ‘粉霜’ ‘Pink Ice’ 155 953 84 384 26 343 18 883 36.8 130 85 8 37 10 ‘黄金喷泉’ ‘Gold Fountain’ 155 941 84 395 26 343 18 860 36.7 130 85 8 37

表 1 绣线菊属叶绿体基因组的基本特征 Table 1 Chloroplast genome feature of genus Spiraea

‘粉霜’和‘黄金喷泉’的叶绿体基因组长度分别为155 953和155 941 bp，均稍小于菱叶绣线菊(155 957 bp)，也小于绝大多数已知绣线菊属植物，除毛枝绣线菊(Spiraea martini)和蒙古绣线菊(S. mongolica)外(表 1), ‘粉霜’和‘黄金喷泉’的叶绿体基因组与其他绣线菊属物种的一致，均为典型的四分体结构，即含有1个LSC、1个SSC和2个IR (图 1)。‘粉霜’叶绿体基因组的LSC、SSC和IR的长度分别为84 384、18 883和26 342 bp，‘黄金喷泉’分别为84 395、18 860和26 343 bp。‘粉霜’与菱叶绣线菊叶绿体基因组的GC含量均为36.8%，而‘黄金喷泉’为36.7%，说明绣线菊属植物叶绿体基因组的GC含量相近。通过重新注释并仔细检查，除渐尖叶粉花绣线菊缺失psaA基因外，该属叶绿体基因组的基因组成高度保守。

图 1(Fig. 1)

图 1 ‘粉霜’和‘黄金喷泉’叶绿体基因组 Fig. 1 Chloroplast genomes of 'Pink Ice' and 'Gold Fountain'

基因注释结果表明, ‘粉霜’与‘黄金喷泉’叶绿体基因组各含有130个基因，包括85个蛋白质编码基因，37个转运RNA基因(tRNA gene)和8个核糖体RNA基因(rRNA gene)。这些基因可划分成光合作用相关基因、自身表达相关基因和其他基因3类(表 2)。‘粉霜’与‘黄金喷泉’叶绿体基因组中这3类基因的种类和数量均完全一致。光合作用相关的基因均有45个，其中5个基因(atpF、petB、petD、ndhA和ndhB)含1个内含子，ndhB分布在IR上, 因此具有2个拷贝。自身表达相关基因有61个基因, 其中rpoC1、rpl2、rps12、rps16、trnA-UGC、trnG- UCC、trnI-GAU、trnK-UUU、trnL-UAA、trnV- UAC均有1个内含子，而clpP含有2个内含子, 其中, 有15个基因分布于IR上(rpl2、rpl23、rps7、rps12、rrn4.5、rrn5、rrn16、rrn23、trnA-UGC、trnI-CAU、trnI-GAU、trnL-CAA、trnN-GUU、trnR-ACG、trnV- GAC)，因此均有2个拷贝。其他基因中的ycf1和ycf2分布于IR上，均有2个拷贝，ycf3含有2个内含子。

表 2(Table 2)

表 2 ‘粉霜’和‘黄金喷泉’叶绿体基因组的注释基因及其归类 Table 2 Annotated genes and their classification in the chloroplast genomes of 'Pink Ice' and 'Gold Fountain' 类别Category 基因功能Gene function 基因Gene 光合作用相关Photosynthesis 二磷酸核酮糖羧化酶大亚基Large subunit of Rubisco rbcL 光系统I Subunit of photosystem I psaA, psaB, psaC, psaI, psaJ 光系统II Subunit of photosystem II psbA, psbB, psbC, psbD, psbE, psbF, psbH, psbI, psbJ, psbK, psbL, psbM, psbN, psbT, psbZ ATP合酶Subunit of ATP synthase atpA, atpB, atpE, atpF*, atpH, atpI 细胞色素b/f复合体Subunit of Cytochrome petA, petB*, petD*, petG, petL, petN NADPH脱氢酶Subunit of NADH ndhA*, ndhB*(2), ndhC, ndhD, ndhE, ndhF, ndhG, ndhH, ndhI, ndhJ, ndhK 细胞色素c C-type cytochrome synthesis ccsA 自身复制相关Self-replication RNA聚合酶RNA polymerase rpoA, rpoB, rpoC1*, rpoC2 核糖体蛋白Large subunit of ribosome rpl2*(2), rpl14, rpl16, rpl20, rpl22, rpl23(2), rpl32, rpl33, rpl36, rps2, rps3, rps4, rps7(2), rps8, rps11, rps12*(2), rps14, rps15, rps16*, rps18, rps19 核糖体RNAs Ribosomal RNA rrn4.5(2), rrn5(2), rrn16(2), rrn23(2) 转运RNAs Transefer RNA trnA-UGC*(2), trnC-GCA, trnD-GUC, trnE-UUC, trnF-GAA, trnfM-CAU, trnG-GCC, trnG-UCC*, trnH-GUG, trnI-CAU (2), trnI-GAU*(2), trnK-UUU*, trnL-CAA(2), trnL-UAA*, trnL-UAG, trnM-CAU, trnN-GUU(2), trnP-UGG, trnQ-UUG, trnR-ACG (2), trnR-UCU, trnS-CGA, trnS-GCU, trnS-GGA, trnS- UGA, trnT-GGU, trnT-UGU, trnV-GAC(2), trnV-UAC*, trnW-CCA, trnY-GUA ATP蛋白酶ATP dependent protease clpP** 其他Other 核糖体Maturase matK 乙酰辅酶A羧化酶亚基Subunit of Acetyl-CoACarboxyase accD 囊膜蛋白基因Envelop membrane protein cemA 保守的开放阅读框Conserved open reading frames ycf1(2), ycf2(2)，ycf3**, ycf4 (2): 含有2个重复单元；*: 含有1个内含子；**: 含有2个内含子。 (2): Two repeat units; *: With one intron; **: With two introns.

表 2 ‘粉霜’和‘黄金喷泉’叶绿体基因组的注释基因及其归类 Table 2 Annotated genes and their classification in the chloroplast genomes of 'Pink Ice' and 'Gold Fountain'

在‘粉霜’叶绿体基因组中共鉴定出49个重复序列，其中正向重复和回文重复各19个，反向重复9个，互补重复2个; 在‘黄金喷泉’中鉴定出49个重复序列，其中正向重复和反向重复均与‘粉霜’一致，但回文重复和互补重复分别为20个和1个; 在菱叶绣线菊中鉴定到49个重复序列，其中2个互补重复、21个回文重复、19个正向重复、7个反向重复(表 3)。

表 3(Table 3)

表 3 绣线菊属植物叶绿体基因组的重复序列 Table 3 Repetitive sequence in chloroplast genomes of Spiraea 物种Species 互补Complement 回文Palindromic 正向Forward 反向Reverse 毛枝绣线菊S. martini 1 21 21 6 绣球绣线菊S. blumei 2 21 19 7 三裂绣线菊S. trilobata 2 21 18 8 渐尖叶粉花绣线菊S. japonica var. acuminata 0 21 15 14 S. insularis 1 17 27 4 蒙古绣线菊S. mongolica 0 21 18 10 珍珠绣线菊S. thunbergii 0 21 18 10 菱叶绣线菊S.×vanhouttei 2 21 19 7 ‘粉霜’ ‘Pink Ice’ 2 19 19 9 ‘黄金喷泉’ ‘Gold Fountain’ 1 20 19 9

表 3 绣线菊属植物叶绿体基因组的重复序列 Table 3 Repetitive sequence in chloroplast genomes of Spiraea

‘粉霜’、‘黄金喷泉’和菱叶绣线菊的叶绿体基因组分别含有67、69和65个简单重复序列(表 4)。其中，单核苷酸重复序列最多，4碱基重复序列次之，2碱基重复序列均为6个。‘粉霜’叶绿体基因组有1个5碱基重复序列，而其余2种菱叶绣线菊均不存在这类重复。本研究的3种菱叶绣线菊均不存在3碱基重复和6碱基重复，其叶绿体基因组的简单重复特征也基本反映了绣线菊属的情况。

表 4(Table 4)

表 4 绣线菊属叶绿体基因组的简单重复序列SSR类型 Table 4 SSR type in chloroplast genomes of Spiraea 物种Species 单碱基重复Mononucleotide 2碱基重复Dinucleotide 3碱基重复Trinucleotide 4碱基重复Tetranucleotide 5碱基重复Pentanucleotide 6碱基重复Hexanucleotide 总计Total 毛枝绣线菊S. martini 55 7 1 10 0 0 73 绣球绣线菊S. blumei 49 6 0 10 0 0 65 三裂绣线菊S. trilobata 51 6 0 10 0 0 67 渐尖叶粉花绣线菊S. japonica var. acuminata 51 8 0 9 1 0 69 S. insularis 55 7 1 9 1 1 74 蒙古绣线菊S. mongolica 50 6 0 9 0 0 65 珍珠绣线菊S. thunbergii 51 5 0 8 0 0 64 菱叶绣线菊S.×vanhouttei 49 6 0 10 0 0 65 ‘粉霜’ ‘Pink Ice’ 50 6 0 10 1 0 67 ‘黄金喷泉’ ‘Gold Fountain’ 52 6 0 11 0 0 69

表 4 绣线菊属叶绿体基因组的简单重复序列SSR类型 Table 4 SSR type in chloroplast genomes of Spiraea

以三裂绣线菊(S. trilobata)的叶绿体全基因组为参考序列，利用mVISTA在线软件基于Shuffle- LAGAN模型检测绣线菊属叶绿体基因组的序列异质性(图 2)。结果表明, 绣线菊属叶绿体基因组较为保守，特别是IR区的变异程度整体较低。变异多发生在大小拷贝区(LSC和SSC区)，特别是基因间隔区(Intergenic spacers: IGSs)的变异信号更多。

图 2(Fig. 2)

图 2 绣线菊属植物叶绿体基因组的比对分析 Fig. 2 Alignment of chloroplast genomes of Spiraea

使用DnaSP6软件的滑动窗口方法筛选‘粉霜’和‘黄金喷泉’叶绿体高变区域(图 3)，结果存在7个高变异区域，分别为trnH_GUG-psbA、trnK_UUU、trnR_ UCU-atpA、trnT_GCU-psbD、ndhC、rpl32和ycf1。

图 3(Fig. 3)

图 3 ‘粉霜’和‘黄金喷泉’叶绿体基因组核苷酸多态性(Pi)的滑动窗口分析 Fig. 3 Sliding window analysis of nucleotide diversity (Pi) in chloroplast genomes of 'Pink Ice' and 'Gold Fountain'

对10个有叶绿体基因组数据的绣线菊属物种进行了系统发育分析(图 4)。结果表明，除1个节点外，其余分支节点均有很高的支持率(> 95%)；菱叶绣线菊、‘粉霜’和‘黄金喷泉’虽有非常近缘的关系，但并未聚成单系；菱叶绣线菊与绣球绣线菊(S. blumei)形成姐妹群，该分支再与‘黄金喷泉’组成姐妹关系，但其支持率非常低(37%)，其关系并未解决。

图 4(Fig. 4)

图 4 基于叶绿体基因组中编码基因构建的绣线菊属植物最大似然树。分支上的数字为支持率 Fig. 4 ML tree of Spiraea based on coding genes of chloroplast genomes. Bootstrap values on each branch

当前，欧美国家已在原产于亚洲的绣线菊属物种中进行了大量的新品种培育，基于粉花绣线菊(S. japonica)培育出‘金焰’绣线菊(‘Gold Flame’)、‘金山’绣线菊(‘Gold Mound’)、‘魔毡’绣线菊(‘Walbuma’^PBR)等大量园艺品种。菱叶绣线菊是麻叶绣线菊(S. cantoniensis)和三裂绣线菊的人工杂交种，广泛栽培于我国的山东、江苏、广东、广西、四川^[33]。基于菱叶绣线菊培育出了很多的彩叶园艺品种，其中‘粉霜’彩叶绣线菊和‘黄金喷泉’菱叶绣线菊是两个性状优良的菱叶绣线菊品种，具有极高的园艺价值。我国是绣线菊属的现代分布和分化中心，种类多达70种(其中47种为我国特有种^[15])，拥有丰富的绣线菊植物资源，但其园艺资源亟待开发。绣线菊属叶绿体基因组的研究还未深入开展，据统计，仅有8个绣线菊属物种具有叶绿体基因组数据^{[14, 17–23]}。通过绣线菊属叶绿体基因组学研究将有助于厘清该属植物的遗传资源及亲缘关系，有利于进一步发掘其园艺资源。

本研究首次测序组装并分析了绣线菊属中两种重要彩叶绣线菊‘粉霜’和‘黄金喷泉’的叶绿体基因组，发现这两种彩叶绣线菊的叶绿体基因组均为典型的四分体结构，即含有1个LSC、1个SSC和2个IR，并与其他绣线菊属物种的叶绿体基因组结构一致，是植物叶绿体基因组中最常见的结构类型^[34]。‘粉霜’和‘黄金喷泉’叶绿体基因组分别为155 953和155 941 bp，是绣线菊属叶绿体基因组中较小的。此外，本研究通过比较绣线菊属叶绿体基因组，发现除渐尖叶粉花绣线菊缺失psaA基因外，该属的叶绿体基因组在基因组成和GC含量均较稳定。一般而言，光合作用相关的psa基因以往在寄生植物中大量缺失，可能与寄生植物无需进行光合作用有关^[35–37]，而这类基因的丢失现象近年来也在其他植物中发现，可能是由于其他基因(如核基因)发挥了其功能有关^[38]。渐尖叶粉花绣线菊叶绿体基因组丢失psa基因的原因值得后续深入研究。

叶绿体基因组的编码区序列进化速率慢，因而相对保守，而基因内含子及基因间隔区，以及简单重复序列均有较高的变异速率，叶绿体基因组内不同的区域具有不同变异速率的特点非常适用于不同分类阶元的系统发育研究。例如，保守的编码基因适用于科、目及以上分类阶元的系统发育分析^[39]，而其非编码区可用于近缘物种或种下的进化研究^[40]。叶绿体简单重复序列(SSR)是研究植物系统与进化的重要遗传标记，主要为单核苷酸重复序列^[41–42]。在本研究中，‘粉霜’、‘黄金喷泉’及菱叶绣线菊的叶绿体基因组均含有最多的单核苷酸重复序列，与以往的研究结果一致^[41]。此外，本研究还发现4碱基重复序列也较多，为后续开展相关研究奠定了基础。除叶绿体简单重复序列外，本研究使用滑动窗口方法筛选出‘粉霜’和‘黄金喷泉’的7个高变异区域，这些高变位点有4个处于编码区，其余位于基因间隔区。在这些高变位点中，trnH_GUG-psbA和ycf1在其他植物中也具有很高的变异，因此曾被推荐为DNA条形码候选位点^[43–44]。此外，其他片段也具有很高的变异率(特别是rpl32)，说明这些位点在绣线菊属中具有作为DNA条形码位点的潜力，为后续开展相关研究提供了重要研究方向。

本研究对10个有叶绿体基因组数据的绣线菊属物种进行了系统发育分析，结果表明几乎所有的分支节点均有很高的支持率(> 95%)。在以往的研究中，菱叶绣线菊与绣球绣线菊、三裂绣线菊等具有非常近缘的亲缘关系，但彼此的关系并未解决^[45]。本研究支持以上结论，并进一步厘清了这些种的关系，说明绣线菊属叶绿体基因组的蛋白质编码区已积累了较大变异，有充足的系统发育信号，为今后大规模开展相关研究奠定了基础。菱叶绣线菊是麻叶绣球和三裂绣线菊的人工杂交种，本研究由于缺少麻叶绣线菊的样品，并不能使用叶绿体基因组进行杂交方向的判定，后续通过引入麻叶绣线菊的样品将有助于解决上述问题。本研究还发现，菱叶绣线菊、‘粉霜’和‘黄金喷泉’虽有近缘关系，但并未聚成单系，这说明两个菱叶绣线菊园艺品种可能分别通过杂交培育的，也可能是由于菱叶绣线菊(尤其是叶的颜色)经历了强烈的人工选择，从而导致其叶绿体基因突变很快积累，并通过无性繁殖的方式保持了这种变异。叶绿体基因组研究可为我们进一步理解绣线菊属及其彩叶园艺品种的亲缘关系提供大量有用信息，为今后发掘该属更多有价值的园艺资源奠定了基础。