麻叶绣线菊叶绿体微卫星特征分析

摘要" 为研究麻叶绣线菊叶绿体基因组简单重复序列（SSR）的特征，本研究基于GenBank数据库中收录的麻叶绣线菊叶绿体基因组序列，利用REPuter软件和MISA软件分析其叶绿体基因组特征、散在重复序列和叶绿体微卫星（cpSSR）特征。结果表明，麻叶绣线菊的叶绿体全基因组序列长度为155 932 bp，是典型的四分体结构；共识别出43个散在重复序列，其中回文重复最多，互补重复最少。其叶绿体全基因组序列上共识别出57个cpSSR位点，主要分布于大单拷贝区。cpSSR重复碱基组成有4种类型，分别为单碱基、二碱基、四碱基和五碱基，其中单碱基重复序列最多，有45个，占重复序列总数的78.95%；以五碱基最少，仅1个，占重复序列总数的1.75%。cpSSR长度区间为10～19 bp，其中10 bp长度最多，占比43.86%；麻叶绣线菊叶绿体cpSSR的拷贝数在3～19，其中以拷贝数为10的序列最多，有25个，占重复序列总数的43.86%。本试验探究了麻叶绣线菊叶绿体微卫星的分布规律和特点，为绣线菊属植物的物种鉴定、群体遗传和系统发育研究提供参考。

关键词" 麻叶绣线菊；叶绿体基因组；散在重复序列；微卫星

中图分类号" S-3"""""" 文献标识码" A"""""" 文章编号" 1007-7731（2025）01-0067-06

DOI号" 10.16377/j.cnki.issn1007-7731.2025.01.014

Characterization analysis of microsatellite in chloroplast of Spiraea cantoniensis Lour.

Abstract" To investigate the characteristics of simple repeat sequences （SSR） in the chloroplast genome of Spiraea cantoniensis Lour.， based on the chloroplast genome sequences of Spiraea cantoniensis Lour. collected from GenBank databases， the chloroplast genome characteristics， scattered repeat sequences， and chloroplast microsatellites（cpSSR） features were analyzed using REPuter and MISA software. The results showed that the complete chloroplast genome sequence was 155 932 bp， and it is a typical tetrad structure; 43 interspersed repeat sequences were identified， with the most palindromic repeats and the least complementary repeats. A total of 57 cpSSR loci were identified on the chloroplast genome sequence， mainly distributed in the LSC region. There were 4 types of cpSSR repeat base compositions， namely mononucleotide， dinucleotide， tetranucleotide， and pentanucleotide. Among them， mononucleotide repeat sequences were the most numerous with 45， which accounted for 78.95% of the total repeat sequences， and pentanucleotide repeat sequences had the least， only 1， which accounted for 1.75% of the total repeat sequences. The length range of cpSSR was 10-19 bp， with the highest length being 10 bp， accounting for 43.86%; the copy number of cpSSR motifs in the chloroplasts of Spiraea cantoniensis Lour. varied from 3 to 19， with the most sequences having a copy number of 10， with 25， accounting for 43.86% of the total repeats. This experiment investigated the distribution patterns and characteristics of chloroplast microsatellites of Spiraea cantoniensis Lour.， providing a reference for species identification， population genetics， and phylogenetic research of Spiraea.

Keywords" Spiraea cantoniensis Lour.; chloroplast genome; interspersed repeat sequences; microsatellite

绣线菊属（Spiraea）是蔷薇科（Rosaceae）绣线菊亚科（Spiraeoideae）落叶类型中的属，在欧洲、亚洲和北美洲均有分布，我国主要分布于云南、甘肃、内蒙古、陕西、四川和山东等地[1]。该属植物具有耐寒、耐旱、适应性强和易于引种栽培等特点，被广泛应用于园林绿化。在城市园林植物造景中，该属植物可以栽植于水岸、石边、山坡和草坪角隅，起到美化景观、净化空气等作用；其花型雅致，叶型多样，且叶色具有丰富的季相变化，可实现夏季赏花、秋季观叶的景观效果[2-3]。麻叶绣线菊（Spiraea cantoniensis Lour.）是绣线菊属灌木，其高可达1.5 m，小枝细瘦，叶片菱状披针形至菱状长圆形；伞形花序，花瓣近圆形或倒卵形，白色，蓇葖果。该植物4—5月开花，7—9月结果[4]。花色鲜艳、叶型多样且株型奇特，具有较高的观赏价值；同时因其萌蘖性强且花期长，在园林绿化中常被用作绿篱，可形成花带并起到阻隔作用。此外，麻叶绣线菊的根、叶和果实具有清热、凉血、祛瘀和消肿止痛等功能，可用于治疗跌打损伤、疥癣等，具有较高的药用价值[5]。

微卫星（Microsatellite）又称简单重复序列（Simple sequence repeats，SSR），是一种常用的以聚合酶链式反应（Polymerase chain reaction，PCR）为基础的DNA分子标记技术，由1～6个重复核苷酸单元组成，其长度在200 bp以内[6]。由于该标记技术具有分辨率高、多态性高以及对DNA模板要求低等优点，被广泛应用于系统分类、遗传多样性分析和种质资源鉴定等研究[7]。叶绿体微卫星又称叶绿体SSR（cpSSR），除具有微卫星标记的优点外，还兼有叶绿体基因组结构简单、单亲遗传等特点[8]。随着数据库中叶绿体基因组数据逐步增加，基于cpSSR分子标记分析种群遗传结构的研究较多，如其在三七[Panax notoginseng （Burkill） F. H. Chen ex C. H. Chow][8]、甘薯[Dioscorea esculenta （Lour.） Burkill][9]、柑橘属（Citrus）[10]、苦参（Sophora flavescens Aiton）[11]和山药（Dioscorea polystachya Turcz.）[12]等植物研究中均有应用。目前，对于麻叶绣线菊的研究大多集中在生理生化、组织培养和植物资源等方面，如杨建伟等[13]研究表明，麻叶绣线菊生长适宜的土壤水分为田间最大持水量的75%，其净光合速率、气孔导度和蒸腾速率在轻度干旱胁迫下较高。刘正兰等[14]阐述了绣线菊植物组织培养技术流程，可在短期内得到大量试管苗。吴小刚[15]对福建省蔷薇科（木本植物）22属148种进行系统调查发现，该地区蔷薇科植物以灌木为主，乔木、抗性树种较少；在园林绿化方面，常以树林、建筑墙角和山石等景物为背景进行成片绿化。该植物分子水平上的研究相对较少。本研究从GenBank数据库中下载了麻叶绣线菊的叶绿体基因组序列数据和注释文件，根据其序列对叶绿体基因组特征、散在重复序列分布，以及微卫星数量、组成和分布特点进行分析，为该属植物物种的品种鉴定、种质遗传多样性分析及序列引物设计提供参考。

1 材料与方法

1.1 试验材料

从GenBank数据库下载麻叶绣线菊的叶绿体基因组完整序列数据和注释文件（GenBank登录号：OR513068），分别以fasta和gb格式保存。

1.2 试验方法

1.2.1 叶绿体基因组的基因构成分析 将获得的叶绿体完整序列数据进行组装、注释后绘制出叶绿体基因组的基本特征。

1.2.2 叶绿体基因组散在重复序列分析 利用在线工具REPuter（https：//bibiserv.cebitec.uni-bielefeld.de/reputer）[16]查找麻叶绣线菊中的散在重复序列，包括正向重复（Forward repeats，F）、反向重复（Reverse repeats，R）、互补重复（Complement repeats，C）和回文重复（Palindromic repeats，P）4种[17]。

1.2.3 叶绿体微卫星分析 使用MISA（https：//webblast.ipk-gatersleben.de/misa/index.php？action=1）对麻叶绣线菊cpSSR进行识别和定位，设置参数：单碱基重复10次以上，二碱基重复5次以上，三碱基重复4次以上，四碱基、五碱基以及六碱基重复3次以上。统计cpSSR分布区域；基于cpSSR碱基类型，分别统计不同重复类型的微卫星长度和个数，并计算其数量、比例和丰度等信息。

2 结果与分析

2.1 麻叶绣线菊叶绿体基因组特征

麻叶绣线菊叶绿体基因组基本特征如表1所示。其序列长度155 932 bp，是叶绿体基因组典型的四分体结构，即1个大单拷贝区（Large single copy，LSC）、1个小单拷贝区（Small single copy，SSC）和2个长度相同、方向相反的反向重复区（Inverted repeats，IR）。其中，LSC区长度84 363 bp，SSC区长度18 883 bp，两个IR区长度均为26 343 bp。查看注释文件可知，该序列共有132个基因，包括84个mRNA基因、37个tRNA 基因、8个rRNA基因和3个假基因。

2.2 麻叶绣线菊叶绿体基因组散在重复序列分布

麻叶绣线菊叶绿体基因组散在重复序列如表2所示。基于REPuter程序，检测到麻叶绣线菊的叶绿体基因组共有43个散在重复序列，其中正向重复（F）14个，反向重复（R）4个，互补重复（C）1个，回文重复（P）24个，占比分别为32.6％、9.3%、2.3%和55.8％。其中30 bp的序列最多，为12个，且4种类型的重复在该长度中均有分布；其次为31、33和34 bp的序列大小，均有5个，32 bp有4个，其余大小的重复序列均有1～2个（图1）。

2.3 麻叶绣线菊cpSSR特征

2.3.1 cpSSR位点分布 麻叶绣线菊叶绿体微卫星分布特征如表3所示。使用MISA在线软件在麻叶绣线菊叶绿体基因组中共识别出57个cpSSR位点（不含复合位点），总长662 bp，占总基因组大小的0.42％。cpSSR分布在基因组292～139 898 bp，平均每2 449 bp出现1个cpSSR位点。结果显示，麻叶绣线菊中，cpSSR在LSC区分布最多，其次分布在SSC区，在IR区分布最少。LSC区中，基因间隔区分布最多，其次分别为外显子区和内含子区。SSC区和IR区中，cpSSR均分布在外显子区。

2.3.2 cpSSR重复碱基组成与丰度 麻叶绣线菊cpSSR重复碱基组成及丰度等信息如表4所示。碱基组成分析结果表明，单碱基重复基序最多，共45个，占cpSSR总量的78.95％，主要为A/T重复，共44个（97.78%），C/G重复1个（2.22％）。其次为四碱基，有6个，占cpSSR总量的10.53％，其中AATA/TATT和TTTA/TAAA在该类型重复中数量最多，分别有2个（33.33％）。再次为二碱基，有5个，占cpSSR总量的8.77％，为AT/TA和TA/TA重复。最后为五碱基，有1个，占cpSSR总量的1.75％，为CAATT/AATTG重复。未发现三碱基、六碱基及以上碱基重复基序。所有类型的cpSSR重复中，单碱基丰度最高，为282.17个/Mb；最低为单碱基中的C/G、四碱基中的AAAT/ATTT、ATAA/TTAT和五碱基中的CAATT/AATTG，丰度均为6.41个/Mb。

2.3.3 cpSSR长度分布 麻叶绣线菊不同类型的cpSSR长度如图2所示，其平均长度12 bp，最长为19 bp，最短为10 bp，其中长度10 bp的cpSSR数量最多，有25个（43.86％），其次为12 bp，有12个（21.05％），最少为14、16和19 bp，各有1个（1.75％）。

2.3.4 cpSSR基序拷贝数分布 不同类型重复序列的基序拷贝数分布如图3所示。麻叶绣线菊cpSSR基序的拷贝数在3～19，不同重复类型的基序拷贝数有较大差异。总体来说，基序越长，则拷贝数越少。单碱基重复序列的拷贝数分布范围较广，在10～19，拷贝数为10的重复序列最多，有25个，占重复序列总数的43.86％。其中A重复7个，T重复12个，G重复1个；拷贝数为14、16和19的序列最少，拷贝数为14的序列仅有1个T重复（1.75％），拷贝数为16的序列仅有1个TTTA重复（1.75％），拷贝数为19的序列仅有1个T重复（1.75％）。

3 结论与讨论

本研究使用REPuter在线程序共检测出麻叶绣线菊叶绿体基因组中有43个散在重复序列，其中回文重复最多，正向重复、反向重复次之，互补重复最少。该分布趋势与同属大部分植物相似，如欧亚绣线菊（Spiraea media Schmidt）、土庄绣线菊（Spiraea ouensanensis H. Lév.）[18]、绣球绣线菊（Spiraea blumei G. Don）、三裂绣线菊（Spiraea trilobata L.）、渐尖叶粉花绣线菊（Spiraea japonica var. acuminata Franch.）、珍珠绣线菊（Spiraea thunbergii Siebold ex Blume）和菱叶绣线菊[Spiraea × vanhouttei （Briot） Carrière][19]等。但在同亚科植物紫叶风箱果（Physocarpus opulifolius ‘Purpurea’）中，其数量分布趋势为正向gt;回文gt;反向gt;互补[20]。

使用MISA软件在总长155 932 bp的麻叶绣线菊叶绿体基因组中检测出57个cpSSR。相对于同属其他植物，麻叶绣线菊cpSSR总数较少。原因除物种差异外，可能与查找SSR数量的参数设置有关，如土庄绣线菊中SSR搜索的数量幅度较宽[18]，会影响搜索数量。石汝杰等[21]研究表明，AT含量与基因组中微卫星的分布有关。

麻叶绣线菊cpSSR分布中，主要以单碱基为主，其中以A/T基序最多，本研究中出现的二、四和五碱基重复类型中也以基序A和T最多，C和G重复较少。Borsch等[22]和赵祺等[23]研究表明，与CG相比，AT只有2个氢键连接，更易断裂，在复制过程中可能出现重组、滑脱等不稳定现象，这些现象一定程度上增加了微卫星形成的概率。同时CG的甲基化现象也使GC重复的概率降低[24]。丰富的AT重复在较多物种中大量存在，如紫叶风箱果[20]、蕨麻[Argentina anserina （L.） Rydb.][25]、蒙古韭（Allium mongolicum Regel）[26]和蝴蝶兰（Phalaenopsis aphrodite Rchb. f.）[27]等。

cpSSR的序列长度对其多态性有一定的影响。麻叶绣线菊cpSSR平均长度12 bp，长度为10 bp的重复序列数最多，其次是12 bp，二者共有37个，占比64.91％。长度大于15 bp的有2个，占比3.51％，说明长度越长，丰度越低。这可能与该物种进化过程中的选择压力有关[24，28]。

叶绿体基因组中的微卫星具有稳定、分子量小的优点，且在种内、种间水平均具有高度多态性。因此，常被广泛应用于群体中的遗传学和进化生物学研究。高通量测序技术的不断发展和普及，极大简化了SSR引物开发工作的流程，减少了工作量。

综上，本研究对麻叶绣线菊叶绿体基因组特征、散在重复序列和叶绿体微卫星（cpSSR）序列特征进行分析，结果表明，麻叶绣线菊的叶绿体全基因组序列长度为155 932 bp，是典型的四分体结构；全基因组序列有43个散在重复序列，其中回文重复最多；在其叶绿体全基因组序列上识别出57个cpSSR位点，主要分布于大单拷贝区。cpSSR重复序列共有4种类型，分别为单碱基、二碱基、四碱基和五碱基，其中单碱基重复序列最多，为45个，占重复序列总数的78.95％。cpSSR长度在10～19 bp，其中10 bp最多，占cpSSR总数的43.86％；麻叶绣线菊cpSSR基序的拷贝数在3～19，其中以拷贝数为10的重复序列最多，有25个，占所有重复序列的43.86％。研究结果为微卫星引物开发、分子育种、系统发育及种质资源利用等研究提供参考。

参考文献

[1] 邢燕红. 绣线菊属植物茶品质成分测定及其资源利用[D]. 呼和浩特：内蒙古农业大学，2010.

[2] 江丹丹. 中国绣线菊属植物的分类学研究[D]. 杨凌：西北农林科技大学，2015.

[3] 王金秋，杨永花，孙朝华，等. 适宜兰州地区栽培的绣线菊特征特性及其园林应用[J]. 寒旱农业科学，2024（2）：193-196.

[4] 刘正兰. 麻叶绣线菊组织培养及植株再生体系的研究[D]. 雅安：四川农业大学，2008.

[5] 姚争争，李斌，杜文鹏，等. 绣线菊属植物化学成分及生物活性研究进展[J]. 中药材，2016，39（4）：934-941.

[6] 田甜，王宏霞. 暴马丁香的叶绿体微卫星特征分析[J]. 分子植物育种，2024，22（8）：2501-2507.

[7] SHAHNAZARI N，NOORMOHAMMADI Z，SHEIDAI M，et al. A new insight on genetic diversity of sweet oranges：caps-SSR and SSR markers[J]. Journal，genetic engineering amp; biotechnology，2022，20（1）：105.

[8] 韩岩. 基于cpDNA分子标记的三七居群遗传多样性研究[D]. 昆明：昆明理工大学，2018.

[9] 王崇，王连军，杨新笋，等. 104个甘薯品种的cpSSR指纹图谱构建及遗传多样性分析[J]. 热带作物学报，2021，42（6）：1549-1556.

[10] 杨程，龚桂芝，彭祝春，等. 基于cpInDel标记和cpSSR标记的柑橘属及近缘属植物亲缘关系[J]. 中国农业科学，2022，55（16）：3210-3223.

[11] 宋芸，张鑫瑞，贺嘉欣，等. 基于叶绿体SSR分子标记的苦参种质资源遗传多样性分析[J]. 作物杂志，2023（1）：30-37.

[12] 张静珍，王连军，雷剑，等. 基于cpSSR标记的山药种质资源DNA指纹图谱构建及遗传多样性分析[J]. 浙江农业学报，2021，33（7）：1222-1233.

[13] 杨建伟，孙桂芳，赵丹，等. 干旱胁迫对麻叶绣线菊光合及荧光参数的影响[J]. 河北农业大学学报，2018，41（3）：49-54.

[14] 刘正兰，刘军，冯秋平. 麻叶绣线菊的组织培养与快速繁殖[J]. 植物生理学通讯，2008，44（3）：527-528.

[15] 吴小刚. 福建省蔷薇科木本植物资源调查分析[J]. 西南林学院学报，2003，23（1）：46-49.

[16] KURTZ S，CHOUDHURI J V，OHLEBUSCH E，et al. REPuter：the manifold applications of repeat analysis on a genomic scale[J]. Nucleic acids research，2001，29（22）：4633-4642.

[17] 赵月梅，杨振艳，赵永平，等. 木犀科植物叶绿体基因组结构特征和系统发育关系[J]. 植物学报，2019，54（4）：441-454.

[18] 宋晴，宿逸然，王铁梅，等. 内蒙古林草过渡区欧亚绣线菊和土庄绣线菊叶绿体基因组比较研究[J]. 北京林业大学学报，2024，46（3）：103-114.

[19] 张巧玲，陈琳，张淑燕，等. 菱叶绣线菊彩叶园艺品种‘粉霜’和‘黄金喷泉’的叶绿体基因组研究[J]. 热带亚热带植物学报，2023，31（4）：510-520.

[20] 黄祥，何梦瑶，王子煊，等. 紫叶风箱果叶绿体基因组特征及绣线菊亚科系统发育分析[J]. 草业学报，2024，33（3）：161-173.

[21] 石汝杰，张元媛，蒋雪梅，等. 铁皮石斛转录组微卫星序列统计分析[J]. 浙江农业科学，2020，61（11）：2241-2244.

[22] BORSCH T，QUANDT D. Mutational dynamics and phylogenetic utility of noncoding chloroplast DNA[J]. Plant systematics and evolution，2009，282（3）：169-199.

[23] 赵祺，余佳兴，秦宇雯，等. 基于高通量测序的菜头肾叶绿体基因组的组装及序列分析[J]. 中草药，2021，52（6）：1744-1750.

[24] 吴东洋，业宁，徐逸卿，等. 狗枣猕猴桃叶绿体基因组微卫星特征分析[J]. 北方园艺，2018（9）：30-35.

[25] 田甜，李军乔，王鑫慈，等. 蕨麻cpSSR序列特征分析[J]. 安徽农业科学，2023，51（12）：94-97.

[26] 王宏霞，李文麒. 蒙古韭叶绿体微卫星特征分析[J]. 北方园艺，2023（3）：1-7.

[27] 张君毅. 蝴蝶兰叶绿体DNA微卫星分析与标记开发[J]. 江西农业学报，2011，23（12）：31-33，37.

[28] 杨芩，付燕，邓洁，等. 蓝莓自花和异花授粉花柱转录组SSR位点分析[J]. 安徽农学通报，2024，30（11）：33-41.