长柄扁桃Amygdalus pedunculata，中文名野樱桃、毛樱桃，多年生灌木，耐寒抗旱，为蔷薇科李亚科桃属扁桃亚属植物，主要分布于中国西北部与毛乌素沙地，是生产生物柴油、活性炭、高品质食用油、蛋白粉、甘油、苦杏仁苷等产品的重要原料和我国重点发展的12种木本油料作物之一；长柄扁桃油是我国新资源食品原料^[1-4]。长柄扁桃属于重点保护植物、濒危植物^[5-8]和药食同源植物^[9-10]。研究发现长柄扁桃种仁含有苦杏仁苷^[11]和多种不饱和脂肪酸^[12]。苦杏仁苷具有抑制肿瘤和炎症、组织器官纤维化等多种药理功能^[13]。不饱和脂肪酸可缓解动脉粥样硬化并有降低血脂作用^[14-15]。长柄扁桃生态效用显著，具有防风固沙稳定荒漠地区脆弱生态系统的作用^{[2-3, 6, 16]}。研究长柄扁桃叶绿体基因组成和系统地位对保护和开发应用濒危植物有重要的学术价值。植物是生产者，叶绿体为光合作用场所，主要存在于植物绿色部位^[17]。叶绿体基因组中有光合作用、能量代谢、同化或未知等多种类型基因^[18-20]。叶绿体基因组是研究高等植物物种形成、种质资源遗传多样性和基因工程技术等方面理论基础和实验材料。长柄扁桃属于蔷薇科桃属Amygdalus，近缘植物的叶绿体基因组多有被报道^[21-28]。人参Panax ginseng^[29]、丹参Salvia miltiorrhiza^[30]叶绿体基因组已被成功组装。叶绿体基因组序列在揭示物种形成起源、物种之间亲缘关系方面始终具有重要价值^[31-33]。长柄扁桃是濒危保护植物、优良木本油料植物、药食同源植物和园林观赏植物，具有十分重要的经济价值和应用价值。本研究以原料植物长柄扁桃为材料，通过二代测序，取得其叶绿体基因NCBI登录号MG602257^[28]。从叶绿体基因序列结构特征、密码子组成分析和IR/SC边界等方面对长柄扁桃的叶绿体基因组进行对比分析，基于全叶绿体基因组数据构建进化树，以期为濒危物种长柄扁桃叶绿体基因工程、分子遗传变异、种质资源开发提供一定理论基础和实验依据。

试验材料为毛乌素沙地自然分布区的长柄扁桃野生生长植株。采集长柄扁桃新鲜幼嫩叶片，于Illumina平台建库高通量测序，在NCBI网站(http://www.ncbi.nlm.nih.gov/)注册得到A. pedunculata叶绿体基因序列登录号。

运用MEGA5.5^[34]、BioEdit^[35]、Geneious^[36]、Mafft7.037^[37]、IRscope (https://irscope.shinyapps.io/irapp/)和CodonW1.4.2软件^[38-40]对NCBI中长柄扁桃叶绿体基因组序列特征进行分析。

采用MEGA5.5与Geneious软件，以Kimura-2为参数基于1 000次Bootstrap value构建邻接法(Neighbor-joining，NJ)系统进化树。

本次高通量测序获得注册号MG602257，与NCBI中登录号NC_037850为同一序列。

将MG602257、KY101153、MG869261^[27]长柄扁桃叶绿体基因组结构与序列情况进行分析比对(图 1–2与表 1)。从图 1中可观察到，长柄扁桃的叶绿体基因组呈现典型的4分区域结构，即大单拷LSC、SSC以及两个大小相同反向重复区IR。3条长柄扁桃叶绿体基因基本特征信息见表 1。长柄扁桃叶绿体基因组总长度MG602257 < MG869261 < KY101153，分别为157 851 bp、157 873 bp和157 948 bp。叶绿体基因4分区域结构中SSC与IR长度较LSC长度小，LSC长度往往最大。长柄扁桃叶绿体基因组LSC长度MG602257 < MG869261 < KY101153，分别为86 052 bp、86 074 bp和86 144 bp。SSC长度变化差异不大，MG602257=MG869261 < KY101153，分别为19 029 bp、19 029 bp和19 056 bp。IRa与IRb长度相同，长柄扁桃叶绿体基因组IRa长度MG602257与MG869261相同，但均大于总长度第一的KY101153，分别为26 385 bp、26 385 bp和26 374 bp。长柄扁桃叶绿体基因总长度序列差22–97 bp，在LSC区长度变化最大，IR区长度变化最小，3个叶绿体基因组中SSC区与IR区长度相对稳定；长柄扁桃的序列A/T/C/G含量以及AT含量数值相同；编码的基因数目129–131个，外显子数目和tRNA、rRNA基因数目在MG602257、MG869261和KY101153无差异。从长柄扁桃现有3条叶绿体基因组长度变化分析来看，其总长度变化规律与分区变化规律无相关性。LSC区变化相对更多，但SSC区与IR区长度保守，在MG602257与MG869261两个叶绿体基因组中两个区域长度均相同。

图 1 长柄扁桃的叶绿体基因组4分区域结构 Fig. 1 Four segments in the chloroplast genomes of A. pedunculata.

图选项

图 2 长柄扁桃的叶绿体基因组LSC和SSC、IR边界比对 Fig. 2 Comparison of the junction between IR, LSC, SSC of the chloroplast genomes from A. pedunculata.

图选项

长柄扁桃3个叶绿体基因组的IR/SSC区的边界如图 2所示。长柄扁桃在IR区边界处情况具有相同点，即MG602257、MG869261和KY101153在JLB/JSB/JSA/JLA四个边界有相同的基因类型和排列顺序，分别为LSC/IRb的rpl22、rps19和rpl2，LSC/IRa的rpl2、rps19和trnH，IRb/SSC的ycf1，SSC/IRa的ndhF。长柄扁桃rps19基因横跨IRb/LSC边界区，rps19基因在MG602257、MG869261和KY101153总长279 bp，长度相同；MG602257和MG869261的rps19基因均在JLB区进入IRb区的长度为187 bp；而在KY101153中，rps19基因进入IRb区的长度为182 bp；即rps19基因在MG602257与KY101153中，表现不同；rps19基因在IR区占据比例较LSC区多，rps19基因长度总差值在3个叶绿体基因组为0，但在IRa区中扩张差值有5 bp。长柄扁桃ycf1基因在MG602257、MG869261和KY101153横跨IRb/SSC边界区，ycf1基因长度不同，总长分别为5 657、5 588、5 684 bp；ycf1基因长度大小与总叶绿体基因组总长变化趋向相同，可能具有正相关性。在LSC/IRa边界JSA中，ndhF基因在MG602257、MG869261中总长度相同均为2 237 bp，表现相同，均扩张到IRa区18 bp；ndhF基因在KY101153总长度小一些，为2 228 bp，扩张到IRa区2 bp。总体上，MG602257与MG869261在边界基因相似性大于KY101153。

核苷酸替换情况和插入缺失在物种进化历程中有十分重要的影响，是分子变化和演化的基础。对长柄扁桃叶绿体基因组的3个序列突变情况进行分析(图 3–5)，共检测到149个突变位点，其中93次为替换事件，56次为插入缺失事件，未发现序列中有小倒位事件。如图 3所示，替换事件发生的主要区域在LSC区的基因间隔区，占总体数目42%，为39次；其次是SSC区的蛋白编码区，为13次。在SC区中总共发生了82次，占88.2%，IR区中发生替换11次，占11.8%。说明SC区较IR区易发生序列变化，IR区较SC区保守。KY101153序列中有89次替换事件，发生次数远大于MG869261，因此，在长柄扁桃3个序列中，MG869261可能与MG602257较KY101153有更高同源性。图中显示转换较颠换发生次数多。KY101153序列中转换发生53次，比颠换多了17次。在SSC区中CDS发现13次替换事件，其中转换发生9次，为颠换发现次数的2.25倍，说明长柄扁桃叶绿体序列中，转换较颠换发生机率较大。

图 3 长柄扁桃的叶绿体基因组替换情况 Fig. 3 The total numbers of substitutions in the chloroplast genomes of A. pedunculata. IGS: intergenic spacer; Intron: intron region; CDS: protein coding region.

图选项

图 4 长柄扁桃的叶绿体基因组不同序列插入缺失事件分布情况 Fig. 4 The distribution of indels in the chloroplast genomes of A. pedunculata.

图选项

图 5 长柄扁桃的叶绿体基因组插入缺失类型与数目 Fig. 5 Types and numbers of indels in the chloroplast genomes of A. pedunculata.

图选项

长柄扁桃3个序列中插入缺失情况不同。如图 4所示，KY101153、MG869261与MG602257发生次数分别为44、5、7，对应百分比分别为79%、9%、12%。KY101153中Indels发生在序列的LSC区、SSC区和IR区；MG869261与MG602257全部分布在LSC区。Indels主要发生在序列的LSC区，发生次数为47，达到84%。图 5显示了长柄扁桃的叶绿体基因组插入缺失类型与数目。共有11种类型，即1、2、3、4、5、6、7、14、15、20、23、27 bp。由图 4可以看到，1 bp为主要类型，其中插入发生数目为12个，缺失发生数目为19个。14 bp至23 bp类型则发生数目均只有1次。

通过Mega5.5、CodonW软件对长柄扁桃叶绿体基因组的密码子组成进行分析。RSCU > 1，说明该密码子的使用频率较高；RSCU=1，该密码子无使用偏好性^[40]。结果表明(图 6–7)，MG602257序列中有52 617个密码子，UUU编码苯丙氨酸，有2 344个，是数目最多的密码子；CGC编码精氨酸，有262个，是数目最少的密码子；编码亮氨酸的密码子数目最多，达到5 135个，占全部密码子总量的9.76%；编码色氨酸的密码子有691个，数目最少，仅占总量1.31%。MG602257密码子RSCU > 1的类型有31个，其中以A/T碱基结尾的有26个，占总数目的83.87%，表明编码长柄扁桃的密码子偏向以A/T结尾。KY101153共有52 649个密码子，在所有密码子中，UUU最多，有2 293个；而CGC最少，只有255个；UAA是最常用的终止密码子，其数量为1 381个，高于UAG (789个)和UGA (1 005个)。MG869261共有52 624个密码子，在所有密码子中，同样是UUU数目最多，有2 319个；而CGC类型的密码子最少，只有255个；UAA数目为1 415，高于UGA和UAG的数量。长柄扁桃3个叶绿体序列的密码子使用情况在图 6与图 7折线图中有相似的变化弯曲，显示出高度保守性。

图 6 长柄扁桃叶绿体基因组密码子的数目对比 Fig. 6 The codon usage numbers of the chloroplast genomes of A. pedunculata.

图选项

图 7 长柄扁桃叶绿体基因组密码子的RSCU值对比 Fig. 7 The relative synonymous codon usage of the chloroplast genomes of A. pedunculata. RSCU: relative synonymous codon usage.

图选项

长柄扁桃和同属物种榆叶梅是多倍体植物，分类地位一直有争议^[41-44]。在分子系统学中，高等植物叶绿体基因是研究植物分类中可依赖的材料。为阐明长柄扁桃在扁桃亚属上的系统进化位置和与近缘种的亲缘关系，本研究选取NCBI中全部已有的长柄扁桃叶绿体基因组，加上已公布的共9个物种的13条叶绿体基因组序列为研究对象构建进化树。选取的植物物种叶绿体基因组具体情况如表 2所示。以MEGA5.5和Geneious软件构建NJ进化树如图 8所示。图 8中结果显示长柄扁桃MG602257、MG869261和KY101153三条序列聚为一支，KY101153与榆叶梅聚为一小支，显示长柄扁桃和榆叶梅亲缘关系相近，研究结果与董山平等学者结论一致^{[4, 41-44]}。

图 8 基于叶绿体全基因组的长柄扁桃NJ系统发育树 Fig. 8 The NJ phylogenetic tree was constructed based on the chloroplast genomes in NCBI. Numbers above branches are bootstrap support and only values larger than 80% were shown.

图选项

长柄扁桃是濒危的重点规划发展的木本油料植物，是生产生物柴油和高档食用油及医药原料来源物种，具有十分重要的开发价值。本研究通过二代高通量测序技术得到了NCBI登录号为MG602257的长柄扁桃叶绿体全基因序列。运用生物信息学方法，对NCBI中长柄扁桃所有的叶绿体基因组进行结构、序列特征和IR/SC区的边界情况对比分析。基于长柄扁桃全部叶绿体基因序列，构建NJ系统进化树。研究发现长柄扁桃MG602257叶绿体基因组全长157 851 bp，分为经典4个区域(1个LSC、2个IR和1个SSC)，3个长柄扁桃叶绿体序列长度相差最少在22 bp，最多达97 bp。IR与SSC序列长度变化较少。长柄扁桃所有叶绿体基因组注释基因类型中，外显子数目相同，均为32个。3个长柄扁桃叶绿体序列中SC区较IR区易发生序列变化，IR区较SC区保守。替换事件发生为93次，主要区域在LSC区的基因间隔区，占总体数目42%；其次是SSC区的蛋白编码区，为13次。在SC区中总共发生了82次，占88.2%。高等植物叶绿体基因组可独立构建可信的植物系统发育进化树。本研究构建的系统发育树显示长柄扁桃与榆叶梅物种亲缘关系较近。由于当前关于长柄扁桃以及其近缘种的叶绿体全基因组数据十分有限，因此，需要获得更多相关叶绿体基因组序列才能后续更细致深入研究。长柄扁桃是濒危物种，同时是优良的生产原料物种，应加强对其生物多样性和遗传多样性的保护。本次研究为长柄扁桃物种保护、系统演化、叶绿体基因工程和种质资源开发应用方面提供相关实验数据和理论依据。