en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。
作者简介:

张梦洁(1999-),硕士研究生,研究方向为林木遗传育种,(E-mail)1481436191@qq.com。

通讯作者:

张汉尧,教授,博士生导师,主要从事林木和酵母遗传育种研究,(E-mail)hanyaoz@163.com。

中图分类号:Q943

文献标识码:A

文章编号:1000-3142(2023)10-1849-12

DOI:10.11931/guihaia.gxzw202207030

参考文献
ASAF S, WAQAS M, KHAN AL, et al. , 2017. The complete chloroplast genome of wild rice (Oryza minuta) and its comparison to related species [J]. Front Plant Sci, 8: 304.
参考文献
BEIER S, THIEL T, MUNCH T, et al. , 2017. MISA-web: a web server for microsatellite prediction [J]. Bioinformatics, 33(16): 2583-2585.
参考文献
CAVALIER-SMITH T, 2002. Chloroplast evolution: secondary symbiogenesis and multiple losses [J]. Current Biol, 12(2): R62-R64.
参考文献
CHEN J, 2007. Flora Reipublicae Popularis Sinicae: Vol. 13 [M]. Beijing: Science Press: 351-352. [陈介, 2007. 中国植物志: 第十三卷 [M]. 北京: 科学出版社: 351-352. ]
参考文献
CHEN Y, 2017. Comparative analysis of structural organization and evolution of Chrysanthemum and Ajania chloroplast genomes [D]. Nanjing: Nanjing Agricultural University: 4-16. [陈韵, 2017. 菊属、亚菊属叶绿体DNA的基因组结构与进化分析 [D]. 南京: 南京农业大学: 4-16. ]
参考文献
DENG YW, LIU XZ, ZHUO DL, et al. , 2021. Structural characteristics and phylogenetic relationship of chloroplast genome of Syzygium acuminatissimum [J]. Guangdong Agric Sci, 48(6): 38-46. [邓演文, 刘晓洲, 卓定龙, 等, 2021. 肖蒲桃叶绿体基因组结构特征及系统发育关系 [J]. 广东农业科学, 48(6): 38-46. ]
参考文献
DU L, ZHANG C, LIU Q, et al. , 2018. Krait: an ultrafast tool for genome-wide survey of microsatellites and primer design [J]. Bioinformatics, 34(4): 681-683.
参考文献
FRAZER KA, PACHTER L, POLIAKOV A, et al. , 2004. VISTA: computational tools for comparative genomics [J]. NuclAcids Res, 32(2): W273-W279.
参考文献
HUANG X, CHU GM, ZHENG XK, et al. , 2022. Analysis of codon usage bias and phylogenetic of chloroplast genome in Nymphaea [J]. J Agric Sci Technol, 24(4): 75-84. [黄祥, 楚光明, 郑新开, 等, 2022. 睡莲属叶绿体基因组密码子偏好性及系统发育分析 [J]. 中国农业科技导报, 24(4): 75-84. ]
参考文献
JI KK, SONG XQ, CHEN CG, et al. , 2020. Codon usage profiling of chloroplast genome in Magnoliaceae [J]. J Agric Sci Technol, 22(11): 52-62. [季凯凯, 宋希强, 陈春国, 等, 2020. 木兰科叶绿体基因组的密码子使用特征分析 [J]. 中国农业科技导报, 22(11): 52-62. ]
参考文献
JIN JJ, YW, YANG JB, SONG Y, et al. , 2020. GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes [J]. Genome Biol, 21(1): 241.
参考文献
KIM KJ, LEE HL, 2004. Complete chloroplast genome sequences from Korean ginseng (Panax schinseng Nees) and comparative analysis of sequence evolution among 17 vascular plants [J]. DNA Res, 11(4): 247-261.
参考文献
KUANG DY, WU H, WANG YL, et al. , 2011. Complete chloroplast genome sequence of Magnolia kwangsiensis (Magnoliaceae): implication for DNA barcoding and population genetics [J]. Genome, 54(8): 663-673.
参考文献
KURTZ S, CHOUDHURI JV, OHLEBUSCH E, et al. , 2001. REPuter: the manifold applications of repeat analysis on a genomic scale [J]. Nucl Acids Res, 29(22): 4633-4642.
参考文献
LI XW, GAO HH, WANG YT, et al. , 2013. Complete chloroplast genome sequence of Magnolia grandiflora and comparative analysis with related species [J]. Sci Sin (Vit), 56(2): 189-198.
参考文献
LIANG YN, HUANG YX, WEI SQ, et al. , 2022. Sequence analysis of SSR loci in the transcriptome of purple bud Liupao tea [J]. J Jiangsu Agric Sci, 50(3): 42-48. [梁燕妮, 黄银霞, 魏诗琴, 等, 2022. 紫芽六堡茶转录组SSR位点序列分析 [J]. 江苏农业科学, 50(3): 42-48. ]
参考文献
LIBRADO P, ROZAS J, 2009. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data [J] Bioinformatics, 25(11): 1451-1452.
参考文献
LIEW Y, CHUA KO, YONG H, et al. , 2022 Complete chloroplast genome of Boesenbergia rotunda and a comparative analysis with members of the family Zingiberaceae [J]. Brazil J Bot: 1-14.
参考文献
LIN S, LIU J, HE X, et al. , 2022. Comprehensive comparative analysis and development of molecular markers for Dianthus species based on complete chloroplast genome sequences [J]. Inter J Mol Sci, 23(20): 12567.
参考文献
LIU C, SHI L, ZHU Y, et al. , 2012. CPGAVAS, an integrated web server for the annotation, visualization, analysis, and GenBank submission of completely sequenced chloroplast genome sequences [J]. BMC Genom, 13: 715.
参考文献
LOHSE M, DRECHSEL O, KAHLAU S, et al. , 2013. OrganellarGenomeDRAW—a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets [J]. Nucl Acids Res, 41(W1): W575-W581.
参考文献
LUO DZ, FU RR, GAN YF, et al. , 2020. Chemical constituents of Syzygium grijsii [J]. J Chin Med Mat, 43(5): 1130-1133. [罗德志, 付容容, 甘宇汾, 等, 2020. 轮叶蒲桃的化学成分研究 [J]. 中药材, 43(5): 1130-1133. ]
参考文献
MORTON BR, 1999. Strand asymmetry and codon usage bias in the chloroplast genome of Euglena gracilis [J]. Proc Natl Acad Sci, 96(9): 5123-5128.
参考文献
MU YN, 2020. Phylogenic and taxonomic research of Hypolepis in China based on the complete chloroplast genome [D]. Harbin: Harbin Normal University: 17-36. [穆雨农, 2020. 基于叶绿体全基因组对中国姬蕨属的系统分类学研究 [D]. 哈尔滨: 哈尔滨师范大学: 17-36. ]
参考文献
NIE X, LV S, ZHANG Y, et al. , 2012. Complete chloroplast genome sequence of a major invasive species, crofton weed (Ageratina adenophora) [J]. PLoS ONE, 7(5): e36 869.
参考文献
PARK I, YANG S, KIM WJ, et al. , 2019. Sequencing and comparative analysis of the chloroplast genome of Angelica polymorpha and the development of a novel indel marker for species [J]. Molecules, 24(6): E1038.
参考文献
SHEN XB, ZHU YJ, XU GB, 2021. Distribution characteristics of SSR loci and development of molecular markers in Taxus fuana [J]. J Cent S Univ For Technol, 41(4): 139-147. [申响保, 朱妍洁, 徐刚标, 2021. 密叶红豆杉SSR位点分布特征及分子标记开发 [J]. 中南林业科技大学学报, 41(4): 139-147. ]
参考文献
SUN ZX, AO PX, BI YF, et al. , 2022. Complete chloroplast genome sequence and characteristics analysis of Medica gosativa ‘Deqin’ [J]. Acta Agr Sin, 30(2): 320-328. [孙志轩, 敖平星, 毕玉芬, 2022. ‘德钦’紫花苜蓿叶绿体基因组序列及特征分析 [J]. 草地学报, 30(2): 320-328. ]
参考文献
WANG L, DONG WP, ZHOU SL, 2012. Structural mutations and reorganizations in chloroplast genomes of flowering plants [J]. Acta Bot Boreal-Occident Sin, 32(6): 1282-1288. [王玲, 董文攀, 周世良, 2012. 被子植物叶绿体基因组的结构变异研究进展 [J]. 西北植物学报, 32(6): 1282-1288. ]
参考文献
WANG XQ, SONG WW, XIAO JJ, et al. , 2021. Phylogeny of Myrtales and related groups based on chloroplast genome [J]. Guihaia, 41(1): 68-80. [王雪芹, 宋卫武, 肖建加, 等. 基于叶绿体基因组探讨桃金娘目及其近缘类群的系统发育关系 [J]. 广西植物, 41(1): 68-80. ]
参考文献
WANG YH, XIE YF, ZHANG ZX, et al. The complete chloroplast genome of Sloanea sinensis and the systematic status of Elaeocarpaceae [J]. Guihaia, 42(1): 39-48. [王一麾, 谢宜飞, 张志翔, 等, 2022. 猴欢喜叶绿体全基因组及杜英科系统地位分析 [J]. 广西植物, 42(1): 39-48. ]
参考文献
XIE H, ZHANG L, ZHANG C, et al. , 2022. Comparative analysis of the complete chloroplast genomes of six threatened subgenus Gynopodium (Magnolia) species [J]. BMC genomics, 23(1): 1-14.
参考文献
XU XL, XU LA, HUAN MR, et al. , 2004. Genetic diversity of microsatellites (SSRs) of natural populations of Quercus variabilis [J]. Hereditas (Beijing), 26(5): 683-688.
参考文献
XUE S, SHI T, LUO W, et al. , 2019. Comparative analysis of the complete chloroplast genome among Prunus mume, P. armeniaca, and P. salicina [J]. Hortic Res, 6(1): 626-638.
参考文献
ZHAO KK, LANDREIN S, BARREET R, et al. , 2019. Phylogeographic analysis and genetic structure of an endemic Sino-Japanese disjunctive genus Diabelia (Caprifoliaceae) [J]. Front Plant Sci, 10: 913.
参考文献
ZHAO Q, YU JX, QIN YW, et al. , 2021. Assembly and sequence analysis of the chloroplast genome of Championella sarcorrhiza based on high-throughput sequencing [J]. Chin Tradit Herb Drugs, 52(6): 1744-1750. [赵祺, 余佳兴, 秦宇雯, 等, 2021, 基于高通量测序的菜头肾叶绿体基因组的组装及序列分析 [J]. 中草药, 52(6): 1744-1750. ]
参考文献
ZHAO YF, 2022. Advances in the application of chloroplast genome in botanical research [J]. Biol Teach, 47(3): 83-85. [赵玉芬, 2022. 叶绿体基因组在植物学研究中的应用进展 [J]. 生物学教学, 47(3): 83-85. ]
参考文献
ZUO WM, ZENG Y, YANG CF, et al. , 2019, High-throughput sequencing of the complete chloroplast genome of Rheum tanguticum and its application in species identification [J]. Chin Trad Herb Drugs, 50(22): 5545-5553. [左文明, 曾阳, 杨春芳, 等, 2019. 基于高通量技术的唐古特大黄叶绿体全基因组测序及应用研究 [J]. 中草药, 50(22): 5545-5553. ]
目录contents

    摘要

    轮叶蒲桃(Syzygium grijsii)系桃金娘科(Myrtaceae)蒲桃属(Syzygium)常绿灌木,其开发前景较好,但其叶绿体基因组特征及系统发育关系尚未有相关报道。为弥补轮叶蒲桃基因组学方面的空缺,该文对轮叶蒲桃的叶绿体基因组进行了系统的研究。运用Illumina高通量测序,并在GetOrganelle平台进行完整组装,同时利用组装好的数据分析轮叶蒲桃叶绿体基因组的结构特征和系统发育关系,其中包括轮叶蒲桃叶绿体基因组结构、功能及特征、密码子偏好性分析、叶绿体基因组的比较分析和系统发育的分析。结果表明:(1)轮叶蒲桃叶绿体基因组大小为158591 bp,包含129个基因。其中,rRNA基因8个,tRNA基因37个,蛋白编码基因84个。分析检测到39个重复序列和84个SSR位点。(2)密码子偏好性分析发现轮叶蒲桃叶绿体基因组中末端存在对A/U的偏性,使用最多的是编码亮氨酸的密码子。(3)与近缘种比较,轮叶蒲桃的边界长度保守,边界处的基因种类与多个蒲桃属物种相似;轮叶蒲桃叶绿体基因组在LSC和SSC区变异度较大,有45处0.010<Pi<0.015,核苷酸多样性水平高。(4)系统发育分析表明,轮叶蒲桃与蒲桃(S. jambos)、滇边蒲桃(S. forrestii)和乌墨(S. cumini)亲缘关系最近,其次是马六甲蒲桃(S. malaccense)。因此,该文结果表明轮叶蒲桃的叶绿体基因组结构保守,特征明显,与同属物种相比,其序列具有相似性,核苷酸具有多样性,并且其与多种桃金娘科物种之间存在亲缘关系。该文对轮叶蒲桃叶绿体基因组进行基因组特征和系统发育分析,为轮叶蒲桃的种质资源鉴定和开发利用奠定理论基础。

    Abstract

    Syzygium grijsii, an evergreen shrub of Myrtaceae and Syzygium, has good prospects for exploitation, but its chloroplast genome characteristics and phylogenetic relationships have not been reported. To fill the gap in the genomics of S. grijsii, a systematic study of the chloroplast genome of S. grijsii was conducted in this paper. Illumina high-throughput sequencing was used to perform complete assembly in the GetOrganelle platform, and the assembled data were also used to analyze the structural features and phylogenetic relationships of the S. grijsii chloroplast genome, including the structure, function and characteristics of the S. grijsii chloroplast genome, codon preference analysis, comparative analysis of the chloroplast genome and phylogenetic analyses. The results were as follows: (1) The chloroplast genome of S. grijsii was 158591 bp in size, which contained 129 genes. Among them, 8 were rRNA genes, 37 were tRNA genes and 84 were protein-coding genes. The analysis detected 39 repetitive sequences and 84 SSR loci. (2) Codon preference analysis revealed that there was a bias towards A/U at the ends of the S. grijsii chloroplast genome, and the most used codon was the one encoding leucine. (3) In comparison with closely related species, the boundary length of S. grijsii was conserved, and the gene species at the boundary were similar to several Syzygium species; the S. grijsii chloroplast genome had a high degree of variability in the LSC and SSC regions, with 45 sites of 0.010 < Pi < 0.015 and a high level of nucleotide diversity. (4) Phylogenetic analysis showed that S. grijsii was most closely related to S. jambos, S. forrestii and S. cumini, followed by S. malaccense. The study concluded that S. grijsii had a conserved chloroplast genome structure with distinct features, sequence similarity and nucleotide diversity compared to Syzygium species, and affinity between it and various Myrtaceae species. In this paper, the genomic characterization and phylogenetic analysis of the chloroplast genome of S. grijsii are conducted to lay the theoretical foundation for the identification and exploitation of the germplasm resources of S. grijsii.

  • 叶绿体(chloroplast,cp)是绿色植物进行光合作用的场所,也是绿色植物所特有的细胞器。叶绿体具有相对自主的遗传系统,有自己的遗传物质,能够进行半自主的复制(Xue et al.,2019)。在高等植物中,大多数植物的叶绿体基因组包括一个大单拷贝区(large single copy,LSC)、一个小单拷贝区(small single copy,SSC)和两个反向的重复区(inverted repeats,IRs)(Li et al.,2013)。叶绿体基因组按功能可分为参与叶绿体生物合成的基因、参与叶绿体基因表达的基因、参与光合作用的基因和一些功能未知的开放阅读框(孙志轩等,2022)。叶绿体拥有简单的基因组结构,它包含大量的遗传信息,基因组内碱基交换的速度处于中间水平,分子进化速度在编码区和非编码区的差别很大,所以叶绿体基因组也被广泛利用在系统发育推断和群体遗传学等领域(赵玉芬等,2022)。近年来,随着高通量测序技术越发成熟,近6 000个物种的叶绿体基因组成功被组装和注释。目前,叶绿体基因组分析已经涵盖了大部分被子植物,例如姜科(Zingiberaceae)(Liew et al.,2022)、兰科(Orchidaceae)(Xie et al.,2022)、石竹属(Dianthus)(Lin et al.,2022)等,这些分析为研究物种变异、物种遗传和系统发育关系等方面做出了卓越的贡献。

  • 桃金娘科(Myrtaceae)主要分布于热带和近热带地区,多为乔木或者灌木,约有150属5 600种,我国现有近100种。蒲桃属(Syzygium)是桃金娘科下的一个重要属,全世界约有500种,中国约有72种(陈介,2007)。因为蒲桃属植物外形美观,抗逆性强,所以大多用于城市绿化,但随着近年研究的深入,蒲桃属中的轮叶蒲桃、蒲桃和洋蒲桃等也被发掘出了药用价值。轮叶蒲桃根和叶可用于祛风散寒,活血化瘀,治疗跌打损伤、伤风感冒和风湿性头痛,具有巨大的开发利用价值(罗德志等,2020)。目前,在蒲桃属中,仅有8个物种通过高通量测序获得了完整的叶绿体基因组,学者们对蒲桃属的主要研究重点包括植物化学成分的分离鉴定以及药物研发等方面。一些学者也探讨过桃金娘科的系统发育关系,但只停留于系统发育分析(王雪芹等,2021)。在蒲桃属中,仅有肖蒲桃的叶绿体基因组被详细分析过,但包括轮叶蒲桃在内的其余物种的叶绿体基因组特征和系统发育相关分析等方面的研究成果则尚未见详细报告。

  • 本文首次对轮叶蒲桃进行测序,并对测序结果进行组装分析。分析内容涵盖了基因组结构、功能及特征、重复序列、SSR分析和系统发育分析等多方面,并且首次对轮叶蒲桃及其亲缘物种进行对比分析,对其物种间遗传变异做出了阐述。轮叶蒲桃开发前景较好,叶绿体基因组资源可以为其物种进化提供有效信息,但轮叶蒲桃叶绿体基因组特征和相关系统发育关系尚未有相关报道,影响了其进一步开发利用。本研究通过使用高通量测序技术对轮叶蒲桃的叶绿体基因组进行解码,并以已发表的密切相关物种肖蒲桃的叶绿体基因组序列作为参考进行拼接,得到完整的轮叶蒲桃叶绿体基因组,并在此基础上进行分析,解析轮叶蒲桃与其他桃金娘科植物的系统发育关系,为轮叶蒲桃在个体发育过程中存在的遗传和变异问题及其开发和系统研究提供理论基础。

  • 1 材料与方法

  • 1.1 实验材料、测序和组装

  • 轮叶蒲桃材料采集于福建龙海,用液氮冷冻处理后交与安诺优达生物技术有限公司,在Illumina NovaSeq 6 000平台上测序(中国浙江),测序共获得4.8 G数据量,测序读长为PE 150。通过GetOrganelle(v 1.6.2)进行数据组装(Jin et al.,2020)。组装后用在线软件CPGAVAS 2(http://47.96.249.172:16019/analyzer/view)进行注释(Liu et al.,2012),使用Geneious(v 9.0.2)进行人工校正。测序产生的原始数据(Submission ID: SUB11257174; BioProject ID: PRJNA826323)和组装完成的叶绿体基因组信息(GenBank accession number: ON243764)均上传至NCBI数据库。

  • 1.2 叶绿体基因组结构、功能分析

  • 通过在线软件CPGAVAS2(http://47.96.249.172:16019/analyzer/view)注释轮叶蒲桃叶绿体基因组后,根据注释信息,用Geneious(v 9.0.2)对轮叶蒲桃叶绿体基因组进行基因组成和基因功能分析(穆雨农,2020)。通过在线软件OGDRAW(https://chlorobox.mpimp-golm.mpg.de/OGDraw.html)绘制轮叶蒲桃叶绿体基因组的物理图谱(Lohse et al.,2009)。

  • 1.3 重复序列和SSR位点分析

  • 使用在线软件REPuter(Kurtz et al.,2001)(https://bibiserv.cebitec. uni-bielefeld)对轮叶蒲桃叶绿体基因组重复序列进行分析,其中包括同向、反向、互补和回文重复。软件的参数设置如下:重复序列数目上限为1 000,最小重复长度为30,汉明距离为3(表示一对重复序列的相似度不能小于90%),编辑距离为默认参数。

  • 通过在线软件MISA(v 2.1)检测轮叶蒲桃叶绿体基因组的简单重复序列(Beier et al.,2017)。参数设置如下:单核苷酸、二核苷酸、三核苷酸、四核苷酸、五核苷酸、六核苷酸的重复次数分别设置为10、5、4、3、3、3,其余均为默认。得到的结果文件用Krait(v1.3.3)软件进行手动修正(Du et al.,2018)。

  • 1.4 密码子偏好性分析

  • 通过CodonW(v 1.4.2)软件对轮叶蒲桃叶绿体基因组进行密码子偏好性分析。为了减少误差,本研究选择的蛋白质编码区长度≥300 bp,并且计算出了轮叶蒲桃叶绿体蛋白编码基因的相对同义密码子使用度RSCU值(relative synonymous codon usage),参数均为初始默认值。

  • 1.5 基因组的比较分析

  • 使用在线网站JSHYCloud(http://cloud.genepioneer.com:9929)比较轮叶蒲桃与4个近缘物种,即白果蒲桃(Syzygium album)、马六甲蒲桃(S. malaccense)、香蒲桃(S. odoratum)和肖蒲桃(S. acuminatissimum)的IR边界(LSC/IRb、IRb/SSC、SSC/IRa和IRa/LSC)信息。

  • 通过在线软件mVISTA(https://genome.lbl.gov/vista/index.shtml)中的Shuffle-LAGAN模型比较轮叶蒲桃与以上4个近缘物种基因组序列差异(Frazer et al.,2004)。

  • 通过Geneious(v 9.0.2)对轮叶蒲桃及4个近缘种进行MAFFT比对(陈韵等,2017),再使用DnaSP(DNA Sequence Polymorphism)进行核酸变异度分析(Librado &Rozas,2009),滑动窗口长度为600 sites,步长为200 sites,其余所有参数均为默认,得出Pi值。

  • 1.6 系统发育分析

  • 为阐明轮叶蒲桃与桃金娘科其他物种的系统发育关系,从GenBank数据库(https://www.ncbi.nlm.nih.gov/nuccore)中下载了桃金娘科22个物种的叶绿体基因组序列,采用Geneious(v 9.0.2)软件(MAFFT插件)比对23种叶绿体基因组全序列。使用RAxML(v 8.2.12)构建最大似然进化树,最佳模型为GTR+F+R2,bootstrap值设置为1 000,以推断各节点的支持率。

  • 2 结果与分析

  • 2.1 叶绿体基因组结构、功能及特征

  • 轮叶蒲桃叶绿体基因组全长158 591 bp(图1),共编码129个基因,其中包括8个rRNA基因,37个tRNA基因,84个蛋白编码基因。轮叶蒲桃叶绿体基因组为环状。其中,大单拷贝区的长度为88 046 bp,GC含量为34.8%;小单拷贝区长度为18 409 bp,GC含量为30.9%;单个反向重复区的长度为26 068 bp,GC含量为42.7%,明显高于LSC区和SSC区。在轮叶蒲桃叶绿体基因组中,总GC含量为37.0%,低于AT(63.0%)含量。由图2可知,白色、灰色和黑色区域为外显子,箭头表示基因的方向。轮叶蒲桃包含了3个外显子,其中两个外显子因在反向重复区而分别具有2个拷贝,5′端位于LSC区,3′端位于IR区。此外,轮叶蒲桃叶绿体基因组中有22个含有内含子的基因,其中在19个基因中发现1个内含子,另外3个基因中发现有两个内含子(表1)。其中,rps12作为一种反式剪切基因。

  • 2.2 重复序列和SSR分析

  • 在轮叶蒲桃叶绿体基因组中可以检测到39个大小为30~45 bp的重复序列,包含了1个反向重复,19个回文重复及19个同向重复,未检测到互补重复(图3)。这些重复序列分布广泛,有18个(包含9个同向重复,9个回文重复)位于ycf2、psaBpsaAtrnS-GCUtrnS-UGAycf1基因上,9个跨区域,8个位于基因间区(IGS,intergenic spacer),4个位于内含子区域,4个最长的重复序列(45 bp)全部位于ycf2基因上。

  • 在轮叶蒲桃叶绿体基因组中共检测到84个SSR位点,平均长度(SSR总长度/SSR总计数)为11.34 bp。其中,单核甘酸60个,二核苷酸2个,三核苷酸8个,四核苷酸13个,五核苷酸1个。此外,轮叶蒲桃叶绿体基因组中大部分SSR位点位于IGS区域(64个),少部分位于基因上(20个)。叶绿体中两个二核苷酸(TA、TC),长度均为10 bp,三核苷酸的长度均为12 bp,四核苷酸重复长度仅有一个为16 bp,其余均为12 bp,唯一的五核苷酸长度为15 bp。

  • 2.3 密码子偏好性分析

  • 对轮叶蒲桃叶绿体基因组的41个大于300 bp的蛋白质编码序列进行密码子偏好性分析(表2)。轮叶蒲桃叶绿体基因密码子中,RSCU值>1的密码子有30个,RSCU最高值为1.98,为亮氨酸中的UUA,最低值为0.33,为亮氨酸中的CUG。其中,有29个以A/U结尾,说明轮叶蒲桃叶绿体基因组密码子更偏向以A/U结尾。轮叶蒲桃叶绿体基因的密码子总数为17 506个,其中1 846个(10.54%)为编码亮氨酸的密码子,仅有195个(1.11%)密码子编码半胱氨酸,表明在轮叶蒲桃叶绿体基因组中使用最多的是亮氨酸,使用最少的是半胱氨酸。在所有密码子中,异亮氨酸中的AUU使用次数最多,次数为786,使用次数最少的是半胱氨酸中的UGC,次数为44。

  • 2.4 叶绿体基因组的比较分析

  • 通过在线网站JSHYCloud比较分析了蒲桃属5个物种的LSC、SSC和IR区的边界位置。结果表明这5个物种的IR区长度较保守(26 068~26 472 bp)。结果如图4所示,5个物种LSC/IRb边界在rps19基因内部,并且扩张程度均一致。轮叶蒲桃和马六甲蒲桃的ndhF基因紧挨着IRb/SSC边界,而白果蒲桃、香蒲桃和肖蒲桃的ndhF基因距边界为53~58 bp。5个物种的SSC/IRa边界都位于ycf1基因上,轮叶蒲桃和马六甲蒲桃的ycf1基因朝SSC区分别进入了4 851 bp和4 857 bp,朝IRa区进入了750 bp。白果蒲桃、香蒲桃和肖蒲桃均有4 857 bp位于SSC区、1 188 bp位于IRa区。

  • 图5结果显示(灰色箭头表示基因的方向,纵坐标表示从50%到100%的一致度;UTR表示非翻译区,CNS表示保守非编码序列),以白果蒲桃为参照序列,轮叶蒲桃的基因组与另外4种植物叶绿体基因组之间有所区别。在变异较大的LSC和SSC区,轮叶蒲桃的变异度较另外4种而言较大。5个物种的IR区较为保守,LSC和SSC区中psbM-psbDpsbZ-trnMCAT)、trnMCAT)-atpEycf4-psbEndhF-ccsA基因之间存在较大的序列差异。

  • 使用DnaSP软件进行分析且经过MAFFT比对后的5种桃金娘科植物叶绿体基因组序列结果如图6所示, X轴表示核苷酸的位置, Y轴表示每个窗口的核苷酸多样性,可显著发现IR区的变异度比LSC和SSC区更低。Pi值最高点为0.015,位于LSC区的psbZ-psaB基因间区。图中有2处Pi=0.015的热点,45处0.010<Pi<0.015,可以看出核苷酸多样性水平较高。

  • 图1 轮叶蒲桃叶绿体基因组图谱

  • Fig.1 Chloroplast genome map of Syzygium grijsii

  • 图2 轮叶蒲桃叶绿体基因组的反式剪切基因组结构示意图

  • Fig.2 Schematic diagram of the trans-sheared genome structure in Syzygium grijsii chloroplast genome

  • 表1 轮叶蒲桃叶绿体基因组基因功能总结

  • Table1 Summary of chloroplast genome gene functions in Syzygium grijsii

  • 注:*表示有一个内含子的基因; **表示有两个内含子的基因;括号内的数字表示多拷贝基因的拷贝数量。

  • Note: *indicates the gene with one intron; **indicates the gene with two introns; the number in parentheses indicates number of copies of a multi-copy gene.

  • 2.5 系统发育分析

  • 结果如图7所示,分支上的数字为1 000次自展支持率。系统发育树中具有20个支持率为100%的拓扑结构节点,遗传变异度为0.01,可信度比较高。蒲桃属、桉属(Eucalyptus)、伞房桉属(Corymbia)和杯果木属(Angophora)亲缘关系较近。在蒲桃属中,轮叶蒲桃与滇边蒲桃和乌墨亲缘关系最近,其次是蒲桃。

  • 图3 轮叶蒲桃叶绿体基因组中的SSR类型和分布

  • Fig.3 Types and distribution of SSR in Syzygium grijsii chloroplast genome

  • 3 讨论与结论

  • 轮叶蒲桃叶绿体基因组全长158 591 bp,总GC含量为37.0%,与已发表的桃金娘科物种叶绿体基因组(王雪芹等,2021)相似。轮叶蒲桃完整叶绿体基因组对A/T具有偏好性(推测可能是因为A/T比G/C更容易解链)(左文明等,2019;赵祺等,2021)。另外,轮叶蒲桃完整叶绿体基因组中有129个基因被编码,这些基因中含有蛋白编码基因84个和tRNA基因37个,与邓演文等(2021)的研究结果一致,表明轮叶蒲桃叶绿体基因组相对保守。可以根据这一性状推测出易变性存在于轮叶蒲桃叶绿体基因组序列,并且易发生变异。

  • 有研究者发现,重复序列分析对于分析叶绿体的基因重排、碱基替换、系统发育和基因组进化有重要的作用(Smith et al.,2002; Nie et al.,2012)。轮叶蒲桃叶绿体基因组中包含了3种重复序列类型,大部分位于基因上,其中,最大的四个重复序列均位于ycf2基因上。在轮叶蒲桃的叶绿体基因组中检测到84个SSR位点,随着基序长度的增加,SSR的频率降低,具有多态性。SSR多态性是由重复单元的伸长或缩短引起,广泛用于种群遗传多样性或种群分类分析(Xu et al.,2004; Zhao et al.,2019)。轮叶蒲桃的单核苷酸均为A/T型,与蒲桃属的肖蒲桃一致(邓演文等,2021)。由前人研究可知,从叶绿体基因组中鉴定出的重复序列和SSR可以作为未来木材生态学、进化和追踪研究的宝贵标记(申响保等,2021;梁燕妮等,2022)。

  • 叶绿体基因对于密码子的使用普遍具有偏好性,而通过密码子的偏好性能反映物种对外界环境的选择压力,在叶绿体基因组的研究中发挥着重要作用(Morton et al.,1999; Asaf et al.,2017)。轮叶蒲桃叶绿体基因组中编码亮氨酸的密码子(1 846,10.54%)使用的频率最多,编码半胱氨酸的密码子(195,1.11%)使用的频率最少。分析表明,轮叶蒲桃密码子具有很强的A/U碱基偏性,与睡莲、木兰等(季凯凯等,2020; 黄祥等,2022)被子植物一致。这说明在被子植物间密码子偏好性具有相似性,这种相似性的产生或许与自然选择、遗传变异或者其他方面有关系。

  • 虽然蒲桃属植物叶绿体基因组的部分特征与其他被子植物基本一致(王玲等,2012),但是轮叶蒲桃和蒲桃属其他物种之间亲缘关系的远近尚不明确。蒲桃属的5个物种边界处的基因种类比较保守。由分析可看出,轮叶蒲桃与白果蒲桃的边界区域比较相似,可能是具有较为接近的起源时间(王一麾等,2022)。IR区的收缩和扩张广泛存在于5个蒲桃属植物中,进行进一步的评估可以揭示轮叶蒲桃及相关类群的进化演变过程。

  • 在全序列比对分析中,参试的5个蒲桃属物种间具有较高的相似性,同时也存在几个变异较大的区域,这些区域的差异位点可以应用于种质资源的鉴定、DNA条形码的开发以及后续其他研究(Kim et al.,2004)。从核苷酸多样性分析可以看出,变异位点的Pi值位于0.00 067~0.015之间,有较高水平的核苷酸多样性,LSC区和SSC区是轮叶蒲桃叶绿体基因组中会发生高频变异的区域,这些区间可以作为进化分析的区域(Park et al.,2019),有助于研究轮叶蒲桃的种间关系。

  • 表2 轮叶蒲桃叶绿体基因密码子偏好性分析

  • Table2 Analysis of codon preference in Syzygium grijsii chloroplast genes

  • 注:带下划线所对应的密码子表示被高频使用的密码子(RSCU>1)。

  • Note: The underlined codon indicates a high frequency codon (RSCU>1) .

  • 本研究从NCBI已公开的植物中选取了包括轮叶蒲桃在内的23种植物,通过系统进化分析确定了轮叶蒲桃的进化地位和亲缘关系。轮叶蒲桃、蒲桃、滇边蒲桃和乌墨构成了一个支持率高的单支,与马六甲蒲桃的关系比较密切,亲缘关系较近。在中国植物志分类系统中,轮叶蒲桃、滇边蒲桃和乌墨属小花蒲桃组,这与本研究中构建的系统发育树一致,表明该文聚类分析结果具有可靠性。虽然蒲桃属物种组成丰富,但是目前所有的系统发育数据较少,轮叶蒲桃系统发育关系的研究弥补了部分空缺,但之后依然需要更加深入的分析。

  • 图4 蒲桃属5种植物叶绿体基因组边界分析

  • Fig.4 Analysis of chloroplast genome boundaries in five Syzygium species

  • 通过高通量测序和组装轮叶蒲桃叶绿体样本,获得了轮叶蒲桃完整的叶绿体基因组。该叶绿体基因组的长度为158 591 bp,具有典型的环状结构,含有129个功能基因。叶绿体基因组含有丰富的SSR位点(A/T作为主要成分),并且以A/U结尾的密码子在轮叶蒲桃叶绿体基因组中更容易被使用。通过叶绿体基因组比较分析可以发现,轮叶蒲桃叶绿体基因组的IR边界保守,与蒲桃属其他物种间有较高的相似性,核苷酸水平高。系统发育树支持率比较高,物种间亲缘关系较近。轮叶蒲桃叶绿体基因组特征和系统发育研究对于了解蒲桃属的不同物种非常重要,有助于开发该属的植物资源,并为进一步研究提供理论基础。

  • 图5 蒲桃属5种植物叶绿体基因组序列比对

  • Fig.5 Sequence comparison of five Syzygium species chloroplast genomes

  • 图6 滑动窗口分析5个蒲桃属植物叶绿体基因组核苷酸多样性

  • Fig.6 Sliding window analysis of chloroplast genomic nucleotide diversity in five Syzygium species

  • 图7 基于叶绿体基因组序列用最大似然法构建的23个桃金娘科物种的系统进化树

  • Fig.7 Phylogenetic trees of 23 Myrtaceae species based on chloroplast genome sequences constructed by maximum likelihood (ML) method

  • 参考文献

    • ASAF S, WAQAS M, KHAN AL, et al. , 2017. The complete chloroplast genome of wild rice (Oryza minuta) and its comparison to related species [J]. Front Plant Sci, 8: 304.

    • BEIER S, THIEL T, MUNCH T, et al. , 2017. MISA-web: a web server for microsatellite prediction [J]. Bioinformatics, 33(16): 2583-2585.

    • CAVALIER-SMITH T, 2002. Chloroplast evolution: secondary symbiogenesis and multiple losses [J]. Current Biol, 12(2): R62-R64.

    • CHEN J, 2007. Flora Reipublicae Popularis Sinicae: Vol. 13 [M]. Beijing: Science Press: 351-352. [陈介, 2007. 中国植物志: 第十三卷 [M]. 北京: 科学出版社: 351-352. ]

    • CHEN Y, 2017. Comparative analysis of structural organization and evolution of Chrysanthemum and Ajania chloroplast genomes [D]. Nanjing: Nanjing Agricultural University: 4-16. [陈韵, 2017. 菊属、亚菊属叶绿体DNA的基因组结构与进化分析 [D]. 南京: 南京农业大学: 4-16. ]

    • DENG YW, LIU XZ, ZHUO DL, et al. , 2021. Structural characteristics and phylogenetic relationship of chloroplast genome of Syzygium acuminatissimum [J]. Guangdong Agric Sci, 48(6): 38-46. [邓演文, 刘晓洲, 卓定龙, 等, 2021. 肖蒲桃叶绿体基因组结构特征及系统发育关系 [J]. 广东农业科学, 48(6): 38-46. ]

    • DU L, ZHANG C, LIU Q, et al. , 2018. Krait: an ultrafast tool for genome-wide survey of microsatellites and primer design [J]. Bioinformatics, 34(4): 681-683.

    • FRAZER KA, PACHTER L, POLIAKOV A, et al. , 2004. VISTA: computational tools for comparative genomics [J]. NuclAcids Res, 32(2): W273-W279.

    • HUANG X, CHU GM, ZHENG XK, et al. , 2022. Analysis of codon usage bias and phylogenetic of chloroplast genome in Nymphaea [J]. J Agric Sci Technol, 24(4): 75-84. [黄祥, 楚光明, 郑新开, 等, 2022. 睡莲属叶绿体基因组密码子偏好性及系统发育分析 [J]. 中国农业科技导报, 24(4): 75-84. ]

    • JI KK, SONG XQ, CHEN CG, et al. , 2020. Codon usage profiling of chloroplast genome in Magnoliaceae [J]. J Agric Sci Technol, 22(11): 52-62. [季凯凯, 宋希强, 陈春国, 等, 2020. 木兰科叶绿体基因组的密码子使用特征分析 [J]. 中国农业科技导报, 22(11): 52-62. ]

    • JIN JJ, YW, YANG JB, SONG Y, et al. , 2020. GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes [J]. Genome Biol, 21(1): 241.

    • KIM KJ, LEE HL, 2004. Complete chloroplast genome sequences from Korean ginseng (Panax schinseng Nees) and comparative analysis of sequence evolution among 17 vascular plants [J]. DNA Res, 11(4): 247-261.

    • KUANG DY, WU H, WANG YL, et al. , 2011. Complete chloroplast genome sequence of Magnolia kwangsiensis (Magnoliaceae): implication for DNA barcoding and population genetics [J]. Genome, 54(8): 663-673.

    • KURTZ S, CHOUDHURI JV, OHLEBUSCH E, et al. , 2001. REPuter: the manifold applications of repeat analysis on a genomic scale [J]. Nucl Acids Res, 29(22): 4633-4642.

    • LI XW, GAO HH, WANG YT, et al. , 2013. Complete chloroplast genome sequence of Magnolia grandiflora and comparative analysis with related species [J]. Sci Sin (Vit), 56(2): 189-198.

    • LIANG YN, HUANG YX, WEI SQ, et al. , 2022. Sequence analysis of SSR loci in the transcriptome of purple bud Liupao tea [J]. J Jiangsu Agric Sci, 50(3): 42-48. [梁燕妮, 黄银霞, 魏诗琴, 等, 2022. 紫芽六堡茶转录组SSR位点序列分析 [J]. 江苏农业科学, 50(3): 42-48. ]

    • LIBRADO P, ROZAS J, 2009. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data [J] Bioinformatics, 25(11): 1451-1452.

    • LIEW Y, CHUA KO, YONG H, et al. , 2022 Complete chloroplast genome of Boesenbergia rotunda and a comparative analysis with members of the family Zingiberaceae [J]. Brazil J Bot: 1-14.

    • LIN S, LIU J, HE X, et al. , 2022. Comprehensive comparative analysis and development of molecular markers for Dianthus species based on complete chloroplast genome sequences [J]. Inter J Mol Sci, 23(20): 12567.

    • LIU C, SHI L, ZHU Y, et al. , 2012. CPGAVAS, an integrated web server for the annotation, visualization, analysis, and GenBank submission of completely sequenced chloroplast genome sequences [J]. BMC Genom, 13: 715.

    • LOHSE M, DRECHSEL O, KAHLAU S, et al. , 2013. OrganellarGenomeDRAW—a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets [J]. Nucl Acids Res, 41(W1): W575-W581.

    • LUO DZ, FU RR, GAN YF, et al. , 2020. Chemical constituents of Syzygium grijsii [J]. J Chin Med Mat, 43(5): 1130-1133. [罗德志, 付容容, 甘宇汾, 等, 2020. 轮叶蒲桃的化学成分研究 [J]. 中药材, 43(5): 1130-1133. ]

    • MORTON BR, 1999. Strand asymmetry and codon usage bias in the chloroplast genome of Euglena gracilis [J]. Proc Natl Acad Sci, 96(9): 5123-5128.

    • MU YN, 2020. Phylogenic and taxonomic research of Hypolepis in China based on the complete chloroplast genome [D]. Harbin: Harbin Normal University: 17-36. [穆雨农, 2020. 基于叶绿体全基因组对中国姬蕨属的系统分类学研究 [D]. 哈尔滨: 哈尔滨师范大学: 17-36. ]

    • NIE X, LV S, ZHANG Y, et al. , 2012. Complete chloroplast genome sequence of a major invasive species, crofton weed (Ageratina adenophora) [J]. PLoS ONE, 7(5): e36 869.

    • PARK I, YANG S, KIM WJ, et al. , 2019. Sequencing and comparative analysis of the chloroplast genome of Angelica polymorpha and the development of a novel indel marker for species [J]. Molecules, 24(6): E1038.

    • SHEN XB, ZHU YJ, XU GB, 2021. Distribution characteristics of SSR loci and development of molecular markers in Taxus fuana [J]. J Cent S Univ For Technol, 41(4): 139-147. [申响保, 朱妍洁, 徐刚标, 2021. 密叶红豆杉SSR位点分布特征及分子标记开发 [J]. 中南林业科技大学学报, 41(4): 139-147. ]

    • SUN ZX, AO PX, BI YF, et al. , 2022. Complete chloroplast genome sequence and characteristics analysis of Medica gosativa ‘Deqin’ [J]. Acta Agr Sin, 30(2): 320-328. [孙志轩, 敖平星, 毕玉芬, 2022. ‘德钦’紫花苜蓿叶绿体基因组序列及特征分析 [J]. 草地学报, 30(2): 320-328. ]

    • WANG L, DONG WP, ZHOU SL, 2012. Structural mutations and reorganizations in chloroplast genomes of flowering plants [J]. Acta Bot Boreal-Occident Sin, 32(6): 1282-1288. [王玲, 董文攀, 周世良, 2012. 被子植物叶绿体基因组的结构变异研究进展 [J]. 西北植物学报, 32(6): 1282-1288. ]

    • WANG XQ, SONG WW, XIAO JJ, et al. , 2021. Phylogeny of Myrtales and related groups based on chloroplast genome [J]. Guihaia, 41(1): 68-80. [王雪芹, 宋卫武, 肖建加, 等. 基于叶绿体基因组探讨桃金娘目及其近缘类群的系统发育关系 [J]. 广西植物, 41(1): 68-80. ]

    • WANG YH, XIE YF, ZHANG ZX, et al. The complete chloroplast genome of Sloanea sinensis and the systematic status of Elaeocarpaceae [J]. Guihaia, 42(1): 39-48. [王一麾, 谢宜飞, 张志翔, 等, 2022. 猴欢喜叶绿体全基因组及杜英科系统地位分析 [J]. 广西植物, 42(1): 39-48. ]

    • XIE H, ZHANG L, ZHANG C, et al. , 2022. Comparative analysis of the complete chloroplast genomes of six threatened subgenus Gynopodium (Magnolia) species [J]. BMC genomics, 23(1): 1-14.

    • XU XL, XU LA, HUAN MR, et al. , 2004. Genetic diversity of microsatellites (SSRs) of natural populations of Quercus variabilis [J]. Hereditas (Beijing), 26(5): 683-688.

    • XUE S, SHI T, LUO W, et al. , 2019. Comparative analysis of the complete chloroplast genome among Prunus mume, P. armeniaca, and P. salicina [J]. Hortic Res, 6(1): 626-638.

    • ZHAO KK, LANDREIN S, BARREET R, et al. , 2019. Phylogeographic analysis and genetic structure of an endemic Sino-Japanese disjunctive genus Diabelia (Caprifoliaceae) [J]. Front Plant Sci, 10: 913.

    • ZHAO Q, YU JX, QIN YW, et al. , 2021. Assembly and sequence analysis of the chloroplast genome of Championella sarcorrhiza based on high-throughput sequencing [J]. Chin Tradit Herb Drugs, 52(6): 1744-1750. [赵祺, 余佳兴, 秦宇雯, 等, 2021, 基于高通量测序的菜头肾叶绿体基因组的组装及序列分析 [J]. 中草药, 52(6): 1744-1750. ]

    • ZHAO YF, 2022. Advances in the application of chloroplast genome in botanical research [J]. Biol Teach, 47(3): 83-85. [赵玉芬, 2022. 叶绿体基因组在植物学研究中的应用进展 [J]. 生物学教学, 47(3): 83-85. ]

    • ZUO WM, ZENG Y, YANG CF, et al. , 2019, High-throughput sequencing of the complete chloroplast genome of Rheum tanguticum and its application in species identification [J]. Chin Trad Herb Drugs, 50(22): 5545-5553. [左文明, 曾阳, 杨春芳, 等, 2019. 基于高通量技术的唐古特大黄叶绿体全基因组测序及应用研究 [J]. 中草药, 50(22): 5545-5553. ]

  • 参考文献

    • ASAF S, WAQAS M, KHAN AL, et al. , 2017. The complete chloroplast genome of wild rice (Oryza minuta) and its comparison to related species [J]. Front Plant Sci, 8: 304.

    • BEIER S, THIEL T, MUNCH T, et al. , 2017. MISA-web: a web server for microsatellite prediction [J]. Bioinformatics, 33(16): 2583-2585.

    • CAVALIER-SMITH T, 2002. Chloroplast evolution: secondary symbiogenesis and multiple losses [J]. Current Biol, 12(2): R62-R64.

    • CHEN J, 2007. Flora Reipublicae Popularis Sinicae: Vol. 13 [M]. Beijing: Science Press: 351-352. [陈介, 2007. 中国植物志: 第十三卷 [M]. 北京: 科学出版社: 351-352. ]

    • CHEN Y, 2017. Comparative analysis of structural organization and evolution of Chrysanthemum and Ajania chloroplast genomes [D]. Nanjing: Nanjing Agricultural University: 4-16. [陈韵, 2017. 菊属、亚菊属叶绿体DNA的基因组结构与进化分析 [D]. 南京: 南京农业大学: 4-16. ]

    • DENG YW, LIU XZ, ZHUO DL, et al. , 2021. Structural characteristics and phylogenetic relationship of chloroplast genome of Syzygium acuminatissimum [J]. Guangdong Agric Sci, 48(6): 38-46. [邓演文, 刘晓洲, 卓定龙, 等, 2021. 肖蒲桃叶绿体基因组结构特征及系统发育关系 [J]. 广东农业科学, 48(6): 38-46. ]

    • DU L, ZHANG C, LIU Q, et al. , 2018. Krait: an ultrafast tool for genome-wide survey of microsatellites and primer design [J]. Bioinformatics, 34(4): 681-683.

    • FRAZER KA, PACHTER L, POLIAKOV A, et al. , 2004. VISTA: computational tools for comparative genomics [J]. NuclAcids Res, 32(2): W273-W279.

    • HUANG X, CHU GM, ZHENG XK, et al. , 2022. Analysis of codon usage bias and phylogenetic of chloroplast genome in Nymphaea [J]. J Agric Sci Technol, 24(4): 75-84. [黄祥, 楚光明, 郑新开, 等, 2022. 睡莲属叶绿体基因组密码子偏好性及系统发育分析 [J]. 中国农业科技导报, 24(4): 75-84. ]

    • JI KK, SONG XQ, CHEN CG, et al. , 2020. Codon usage profiling of chloroplast genome in Magnoliaceae [J]. J Agric Sci Technol, 22(11): 52-62. [季凯凯, 宋希强, 陈春国, 等, 2020. 木兰科叶绿体基因组的密码子使用特征分析 [J]. 中国农业科技导报, 22(11): 52-62. ]

    • JIN JJ, YW, YANG JB, SONG Y, et al. , 2020. GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes [J]. Genome Biol, 21(1): 241.

    • KIM KJ, LEE HL, 2004. Complete chloroplast genome sequences from Korean ginseng (Panax schinseng Nees) and comparative analysis of sequence evolution among 17 vascular plants [J]. DNA Res, 11(4): 247-261.

    • KUANG DY, WU H, WANG YL, et al. , 2011. Complete chloroplast genome sequence of Magnolia kwangsiensis (Magnoliaceae): implication for DNA barcoding and population genetics [J]. Genome, 54(8): 663-673.

    • KURTZ S, CHOUDHURI JV, OHLEBUSCH E, et al. , 2001. REPuter: the manifold applications of repeat analysis on a genomic scale [J]. Nucl Acids Res, 29(22): 4633-4642.

    • LI XW, GAO HH, WANG YT, et al. , 2013. Complete chloroplast genome sequence of Magnolia grandiflora and comparative analysis with related species [J]. Sci Sin (Vit), 56(2): 189-198.

    • LIANG YN, HUANG YX, WEI SQ, et al. , 2022. Sequence analysis of SSR loci in the transcriptome of purple bud Liupao tea [J]. J Jiangsu Agric Sci, 50(3): 42-48. [梁燕妮, 黄银霞, 魏诗琴, 等, 2022. 紫芽六堡茶转录组SSR位点序列分析 [J]. 江苏农业科学, 50(3): 42-48. ]

    • LIBRADO P, ROZAS J, 2009. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data [J] Bioinformatics, 25(11): 1451-1452.

    • LIEW Y, CHUA KO, YONG H, et al. , 2022 Complete chloroplast genome of Boesenbergia rotunda and a comparative analysis with members of the family Zingiberaceae [J]. Brazil J Bot: 1-14.

    • LIN S, LIU J, HE X, et al. , 2022. Comprehensive comparative analysis and development of molecular markers for Dianthus species based on complete chloroplast genome sequences [J]. Inter J Mol Sci, 23(20): 12567.

    • LIU C, SHI L, ZHU Y, et al. , 2012. CPGAVAS, an integrated web server for the annotation, visualization, analysis, and GenBank submission of completely sequenced chloroplast genome sequences [J]. BMC Genom, 13: 715.

    • LOHSE M, DRECHSEL O, KAHLAU S, et al. , 2013. OrganellarGenomeDRAW—a suite of tools for generating physical maps of plastid and mitochondrial genomes and visualizing expression data sets [J]. Nucl Acids Res, 41(W1): W575-W581.

    • LUO DZ, FU RR, GAN YF, et al. , 2020. Chemical constituents of Syzygium grijsii [J]. J Chin Med Mat, 43(5): 1130-1133. [罗德志, 付容容, 甘宇汾, 等, 2020. 轮叶蒲桃的化学成分研究 [J]. 中药材, 43(5): 1130-1133. ]

    • MORTON BR, 1999. Strand asymmetry and codon usage bias in the chloroplast genome of Euglena gracilis [J]. Proc Natl Acad Sci, 96(9): 5123-5128.

    • MU YN, 2020. Phylogenic and taxonomic research of Hypolepis in China based on the complete chloroplast genome [D]. Harbin: Harbin Normal University: 17-36. [穆雨农, 2020. 基于叶绿体全基因组对中国姬蕨属的系统分类学研究 [D]. 哈尔滨: 哈尔滨师范大学: 17-36. ]

    • NIE X, LV S, ZHANG Y, et al. , 2012. Complete chloroplast genome sequence of a major invasive species, crofton weed (Ageratina adenophora) [J]. PLoS ONE, 7(5): e36 869.

    • PARK I, YANG S, KIM WJ, et al. , 2019. Sequencing and comparative analysis of the chloroplast genome of Angelica polymorpha and the development of a novel indel marker for species [J]. Molecules, 24(6): E1038.

    • SHEN XB, ZHU YJ, XU GB, 2021. Distribution characteristics of SSR loci and development of molecular markers in Taxus fuana [J]. J Cent S Univ For Technol, 41(4): 139-147. [申响保, 朱妍洁, 徐刚标, 2021. 密叶红豆杉SSR位点分布特征及分子标记开发 [J]. 中南林业科技大学学报, 41(4): 139-147. ]

    • SUN ZX, AO PX, BI YF, et al. , 2022. Complete chloroplast genome sequence and characteristics analysis of Medica gosativa ‘Deqin’ [J]. Acta Agr Sin, 30(2): 320-328. [孙志轩, 敖平星, 毕玉芬, 2022. ‘德钦’紫花苜蓿叶绿体基因组序列及特征分析 [J]. 草地学报, 30(2): 320-328. ]

    • WANG L, DONG WP, ZHOU SL, 2012. Structural mutations and reorganizations in chloroplast genomes of flowering plants [J]. Acta Bot Boreal-Occident Sin, 32(6): 1282-1288. [王玲, 董文攀, 周世良, 2012. 被子植物叶绿体基因组的结构变异研究进展 [J]. 西北植物学报, 32(6): 1282-1288. ]

    • WANG XQ, SONG WW, XIAO JJ, et al. , 2021. Phylogeny of Myrtales and related groups based on chloroplast genome [J]. Guihaia, 41(1): 68-80. [王雪芹, 宋卫武, 肖建加, 等. 基于叶绿体基因组探讨桃金娘目及其近缘类群的系统发育关系 [J]. 广西植物, 41(1): 68-80. ]

    • WANG YH, XIE YF, ZHANG ZX, et al. The complete chloroplast genome of Sloanea sinensis and the systematic status of Elaeocarpaceae [J]. Guihaia, 42(1): 39-48. [王一麾, 谢宜飞, 张志翔, 等, 2022. 猴欢喜叶绿体全基因组及杜英科系统地位分析 [J]. 广西植物, 42(1): 39-48. ]

    • XIE H, ZHANG L, ZHANG C, et al. , 2022. Comparative analysis of the complete chloroplast genomes of six threatened subgenus Gynopodium (Magnolia) species [J]. BMC genomics, 23(1): 1-14.

    • XU XL, XU LA, HUAN MR, et al. , 2004. Genetic diversity of microsatellites (SSRs) of natural populations of Quercus variabilis [J]. Hereditas (Beijing), 26(5): 683-688.

    • XUE S, SHI T, LUO W, et al. , 2019. Comparative analysis of the complete chloroplast genome among Prunus mume, P. armeniaca, and P. salicina [J]. Hortic Res, 6(1): 626-638.

    • ZHAO KK, LANDREIN S, BARREET R, et al. , 2019. Phylogeographic analysis and genetic structure of an endemic Sino-Japanese disjunctive genus Diabelia (Caprifoliaceae) [J]. Front Plant Sci, 10: 913.

    • ZHAO Q, YU JX, QIN YW, et al. , 2021. Assembly and sequence analysis of the chloroplast genome of Championella sarcorrhiza based on high-throughput sequencing [J]. Chin Tradit Herb Drugs, 52(6): 1744-1750. [赵祺, 余佳兴, 秦宇雯, 等, 2021, 基于高通量测序的菜头肾叶绿体基因组的组装及序列分析 [J]. 中草药, 52(6): 1744-1750. ]

    • ZHAO YF, 2022. Advances in the application of chloroplast genome in botanical research [J]. Biol Teach, 47(3): 83-85. [赵玉芬, 2022. 叶绿体基因组在植物学研究中的应用进展 [J]. 生物学教学, 47(3): 83-85. ]

    • ZUO WM, ZENG Y, YANG CF, et al. , 2019, High-throughput sequencing of the complete chloroplast genome of Rheum tanguticum and its application in species identification [J]. Chin Trad Herb Drugs, 50(22): 5545-5553. [左文明, 曾阳, 杨春芳, 等, 2019. 基于高通量技术的唐古特大黄叶绿体全基因组测序及应用研究 [J]. 中草药, 50(22): 5545-5553. ]