Characteristics and Evolution of Panax Chloroplast Genomes
-
摘要:目的 对人参属(Panax L.)物种叶绿体基因组特征及其系统发育进行研究,为我国人参属资源的遗传学研究和开发利用提供理论依据。方法 基于14种人参属植物叶绿体基因组序列,利用生物信息学软件,对叶绿体基因组特征、序列重复、结构变异、基因进化和系统发育进行分析。结果 人参属物种叶绿体基因组均为典型的四分体结构,包含114个unique基因。长重复序列主要为回文重复和正向重复,30~39 bp的重复序列最多,SSR大多为A/T重复,单核苷酸重复是最丰富的类型。人参属叶绿体基因组未发生基因重排,反向重复区(IR)与单拷贝区边界高度保守,鉴定的12个核苷酸高度可变热点中7个位于大单拷贝(LSC)区,5个位于小单拷贝(SSC)区。根据dN/dS比率,发现功能未知基因clpP、ycf1和ycf2受正选择作用。系统发育分析显示,屏边三七和三叶参位于基部支系,四倍体人参和西洋参与其他二倍体物种聚在不同支系,三七、竹节参和越南参则亲缘关系较近。结论 人参属叶绿体基因组基因数目和顺序一致,基因组结构保守,重复序列数目和类型存在差异,单拷贝区核苷酸多态性高于IR区,正选择基因可能与物种的生态适应性有关。Abstract:Objective Characteristics and phylogeny of chloroplast genomes of the medicinally and economically valuable species in Panax genus were studied.Methods Using bioinformatics software, the properties, repeats, structural variation, evolution, and phylogeny of the genomes of chloroplasts from 14 ginseng species were analyzed.Results The genomes consisted of typical quadripartite structure with 114 unique genes. The long repeats in them were mainly of palindromic and forward types with a length between 30 bp and 39 bp. The simple sequence repeats were largely A/T type and most abundantly mononucleotides. No gene rearrangement occurred in the genomes was observed. The boundary between the inverted repeat region and the single copy region was highly conserved. Of the 12 regions with highly variable nucleotides, 7 were in the large and 5 in the small single copy region. Indicated by the dN/dS ratios, the positive selection could occur on clpP, ycf1, and ycf2 with unknown functions. The phylogenetic analysis showed that P. stipulenatus and P. trifolius were in the basal lineage, the tetraploid P. ginseng and P. quinquefolius separated from other diploid species, while P. notoginseng, P. japonicus, and P. vietnamensis closely related.Conclusion The chloroplasts of the ginseng species examined were basically same in number and order of genomes, conservative in structure, but divert in number and type of repeats. The nucleotide polymorphism of the chloroplasts was higher in single copy region than inverted repeat regions. The positive selection genes identified in the study might result from the ecological adaptation of these Panax species.
-
Keywords:
- Panax L. /
- chloroplast genome /
- repeat /
- variation /
- phylogeny
-
0. 引言
【研究意义】随着科技的进步,现代农业产业不断往高度集约化与高效种植方向发展,现代化农业设施得到广泛应用。我国设施园艺栽培面积高居世界首位,2016年已达到476.5万hm2,产值超1.46万亿元[1]。农业设施的补光问题日益凸显,人造光源在植物光照中发挥着越来越重要的作用,试验表明LED光源相对传统光源在植物光照方面有明显的优势[1]。以纯人工光为光源的植物工厂代表了设施农业的发展趋势,植物对光吸收的效率决定了设施农业产业的经济效益,因此寻找植物生长所需的最佳复合光谱是关键。植物特征光谱代表了植物在某个生长周期内的最佳光需求[2],若能够依据植物特征光谱精准给光,将有效促进植物产量和品质的提高,从而提高设施农业产业模式的社会经济效益[3]。植物特征光谱是以植物在多种单色光条件下的光合响应曲线为基础进行定义的[2],因此,研究植物在单色光照射条件下的光合响应曲线及其数学模型,对植物光照应用领域具有重要理论价值与实践指导意义。【前人研究进展】当前,在设施农业中普遍采用LED人工光源补光或进行全人工光栽培,对植物在不同光照环境下表现出的相应特性的研究也越来越多[4-14]。植物的生长及其有效成分的积累与光的波长、强度和光周期密切相关[1]。对于以白光为基础的光合响应曲线模型的研究也有许多成果[15-22]。【本研究切入点】目前的研究基本上都采用红绿蓝三色光为光源,此类光的光谱仅为植物光合有效辐射范围的一部分,无法反映植物光合响应的全貌。仅用三色光对植物进行光照培育,可能会导致植物因为部分波长的光的缺失而出现相应的不良生长现象。同时,在对植物光合响应模型的相关研究中,有关植物在单色光条件下——即将光波长作为光合响应重要影响参数的研究有待深入进行。【拟解决的关键问题】以获取植物特征光谱——即植物在一定环境(非光)中进行光合作用最适宜的光谱、光照强度和光照周期为出发点,选择10种具有代表性的植物幼苗作为试验对象,采用自主研制的多种单色光发生器搭建光合作用(速率)测试系统[2],开展植物在多种单色光条件下的光合响应研究,根据单色光条件下光合响应饱和点实测数据,分析、建立具有较高拟合度的光合响应曲线及数学模型,以期发现植物在不同波长光照射下所呈现光合响应特性的一般规律,为最终形成植物特征光谱提供理论依据和数据支撑。
1. 材料与方法
1.1 试验材料
2021年6月,样本取自南平市森科种苗有限公司苗圃,该苗圃位于延平区西芹镇坑底村,海拔60 m,湿度65%~100%,介于北纬26°15′~26°52′,东经117°50′~118°40′,属中亚热带海洋季风气候,年均气温17.3 ℃(夏季温度30~40 ℃),年降水量1 669 mm。苗圃为露天的普通林地培育树苗,光照较强。在该苗圃取一年生长势良好的10个品种树苗:红椎、马褂木、红花荷、苦槠、乐昌含笑、闽楠、木荷、刨花楠、乳源木莲和油茶。限于篇幅,选择两个具有代表性的树苗红椎和马褂木的测试数据进行分析研究。
1.2 试验仪器
自主研制的多种单色光发生器,该设备由一种双驼峰型光谱的发光装置[23],它产生24~30种单色光,测试试验选择22种单色光,其峰值波长及其对应的半峰全宽如表1所示,该装置具有任意选择单色光的波长并调节其光强的功能,可以根据需要选择单色光和相应的光强。CIRAS-3光合作用测试仪(美国,PP SYSTEMS公司),能够测试在一定光强条件下的植物光合作用速率、气孔导度和蒸腾速率等重要的物理参数。由这两台设备组成普通光植物特征光谱测试系统[2]。
表 1 多种单色光发生器产生的22种单色光Table 1. Twenty-two monochromatic lights produced by light generating device波长
Wavelength λ/nm半峰全宽
Full width of
half peak ∆λ/nm波长
Wavelength λ/nm半峰全宽
Full width of
half peak ∆λ/nm420 13.9 575 22.6 435 17.9 600 20.3 445 12.9 620 17.9 450 19.4 630 17 465 15.6 650 20.5 475 25 660 11.8 500 24.2 675 26.3 520 20.5 690 26.6 535 18.7 700 33.7 550 22.4 730 28.2 565 15.7 750 28.8 1.3 试验方法
按光合作用测试仪使用的常规要求,选择树苗上第3~4个发育正常的叶片,置于叶室里,采用植物特征光谱试验方法[2]进行测试,即任选多种单色光发生器产生的一种单色光,让其均匀地照射在光合作用测试系统的叶室上,通过光合作用测试系统就可以获得被测试的植物在该单色光此刻对应的光强下的光合作用速率。为了避免杂光的干扰,从出光口到叶室,采取遮光措施。按由小到大的顺序调整该单色光的光合有效辐射强度PAR(以下简称光强I/µmol·m−2·s−1,分别为10,20,30,40,60,80,100,150,200,300,400,……,最大值),获得该单色光在不同光强条件下的光合作用速率,即光合作用速率与光强的关系曲线,通常称之为单色光植物光合作用响应曲线(以下简称光合响应曲线)。
1.4 数据处理的工具与方法
应用Matlab 2012b软件中的Curve Fitting Tool曲线拟合工具箱(以下简称CFT工具)对所获得的光合作用速率及其对应的光强数据进行曲线拟合,获得每一种单色光的光合响应曲线,该工具箱采用的是非线性最小二乘法的信任区域理论(Trust Region Methods)求解模型参数(以下简称TRM方法)。
1.5 光合响应曲线的模型
根据有两类不同的单色光光合响应曲线的试验结果,考虑到不同波长的光合响应曲线的差异,分别采用不同的模型来描述。
对于渐近线型植物单色光光合响应曲线采用以下模型来描述。
Pn=PnB(1 - e - αλ(I - Ic)) (1) 式中
Pn 为净光合速率,I 为光强,Ic 为补偿点对应的光强,对于特定的树种和特定的波长Ic 为常数,PnB 为准饱和点光强对应的光合速率,αλ 为响应系数。由式(1)可知,当I=Ic 时,Pn=0 ,这与试验测试结果相符。对于有拐点的植物单色光光合响应曲线,采用以下模型来描述。
Pn=ηλ1−βλI1+γλI(I−Ic) (2) 式中
Pn 为净光合速率,I 为光强,Ic 为补偿点对应的光强,ηλ 为补偿响应系数,βλ 为修正系数,γλ=ηλ/Pmax 为补偿比率系数,其中{P_{\max }} 为植物的最大光合作用速率。由式(2)可知,当I = {I_c} 时,{P_n} = 0 ,同样与试验测试结果相符。1.6 渐近线型光合响应曲线的准饱和点光强计算方法
对于渐近线型植物光合响应曲线不存在光合作用速率最大值
{P_{\max }} ,在植物光照中,当光强达到一定值时,植物的光合作用速率随着光强增加其增量开始减少,当光强达到某一上限值时,再增加光强,对植物的生长没有实际意义,在纯人工光植物工厂中,植物光照的光强通常在200~400 µmol·m−2·s−1[3]。因此,有必要根据植物栽培的实际来确定合理的光强上限值,该值就称为准饱和点光强。可将式(1)求一阶导数得:{P'_n} = {\alpha _\lambda }{P_{nB}}{e^{ - {\alpha _\lambda }\left( {I - {I_c}} \right)}} (3) 采用TRM方法进行求解,可获得相关的参数,即a=
{P_{nB}} ,b={\alpha _\lambda } ,c={I_c} 。如果确定{p'_n} 值,即可由式(3)求出{I'_{sat}} ,即{p'_n} 值确定后,由式(3)计算出的理论计算值{I'_{sat}} 为准饱和点光强。1.7 有拐点植物光合响应曲线的饱和点光强计算方法
对应有拐点光合响应曲线的饱和点光强的确定,就是将式(2)求一阶导数,可得:
P_n^{\prime}=\frac{\eta_\lambda\left(1-\beta_\lambda I\right)-\eta_\lambda \beta_\lambda\left(I-I_c\right)}{1+\gamma_\lambda I}-\frac{\eta_\lambda \beta_\lambda\left(1-\beta_\lambda I\right)\left(I-I_c\right)}{\left(1+\gamma_\lambda I\right)^2} (4) 当式(4)等于0时,就可获得饱和点光强
{I'_{sat}} 的表达式,即{I'_{sat}} = \frac{{\sqrt {{{1 + {{{\gamma _\lambda }} \mathord{\left/ {\vphantom {{{\gamma _\lambda }} {{\beta _\lambda }}}} \right. } {{\beta _\lambda }}} + \left( {{\beta _\lambda } + {\gamma _\lambda }} \right){I_c}{\gamma _\lambda }} \mathord{\left/ {\vphantom {{1 + {{{\gamma _\lambda }} \mathord{\left/ {\vphantom {{{\gamma _\lambda }} {{\beta _\lambda }}}} \right. } {{\beta _\lambda }}} + \left( {{\beta _\lambda } + {\gamma _\lambda }} \right){I_c}{\gamma _\lambda }} {{\beta _\lambda }}}} \right. } {{\beta _\lambda }}}} - 1}}{{{\gamma _\lambda }}} (5) 根据试验测试获取马褂木数据,采用TRM方法进行求解,得到马褂木对各波长对应的
{\eta _\lambda } 补偿响应系数、{\beta _\lambda } 修正系数、{\gamma _\lambda } 补偿比率系数和补偿点光强{I_c} 。再由式(5)就可以获取马褂木的饱和点光强理论计算值
{I'_{sat}} 。2. 结果与分析
2.1 光合响应曲线及其特点
选择22种单色光,采用植物特征光谱的试验方法,进行试验测试,获得的数据应用CFT工具进行模拟,分别获得了10种树苗的22种波长的光合响应曲线。
试验结果表明,树苗的光合响应曲线有两类:一类是渐近线型光合响应曲线,另一类是有拐点的光合响应曲线。具有代表性的树种分别为红椎和马褂木。红椎的光合响应曲线是渐近线型曲线,这类树种的光合响应曲线没有极值;马褂木的光合响应曲线则是有拐点的曲线,这类曲线有极值。由此可见,同一树苗对不同波长的光合响应曲线是不同的,不同树苗对同一波长的光合响应曲线也明显不同。
从图1~3还可以看出光强在20~100 µmol·m−2·s−1光合响应曲线接近是线性的,超过100 µmol·m−2·s−1以后光合响应曲线的特异性表现得比较明显,说明植物光合响应的差别在100 µmol·m−2·s−1以后表现得更为显著。
2.2 渐近线型光合响应曲线的准饱和点光强
表2为红椎的光合响应曲线对应22个单色光的
{P_{nB}} ,{\alpha _\lambda } ,{I_c} 参数值。从表2可以看出,各种单色光的准饱和点光强对应的最大光合作用速率Pnb是不同的,该点构成的谱线有6个峰值,分别是:3.50、3.95、4.47、3.51、3.51、5.64,对应的波长分别为445、535、565、630、565、750 nm;补偿点的实测值与模型的理论值的相关系数的平方值R2大部分在0.99以上或只有少数接近0.99,说明该模型应用于补偿点是成功的。表3为红椎的准饱和点光强的实测值{I_{sat}} 与{p'_n} = 0.001 时的理论计算值{I'_{sat}} 的比较。从表3可以看出有8个波长没有对应的实测值,无法测出准饱和点的实测值,即按照设定的连续3档光强值,测定的光合作用速率的读数误差小于1%,则判定中间挡光强值为准饱和点的光强值,不满足该条件的测试结果则不予采用,即没有对应的实测值。其他确认的实测值与理论值均在设定值档位差值(±100)的范围之内,因此,可认为实测值与理论值基本上是吻合的,即式(1)模型是正确的。表 2 红椎的光合响应曲线对应22个单色光的参数值Table 2. Measurements related to 22 monochromatic lights on LRC of C. hystrix波长
Wavelength
λ/nmPnb Ic/
(μmol·m−2·s−1)Ic′/
(μmol·m−2·s−1)αλ/
(m2·s−1·μmol−1)R2 波长
Wavelength
λ/nmPnb Ic/
(μmol·m−2·s−1)Ic′/
(μmol·m−2·s−1)αλ/
(m2·s−1·μmol−1)R2 420 3.45 5 12 0.015 0.995 2 575 3.14 15 14 0.01 0.995 8 435 3.47 10 13 0.013 0.996 4 600 3.45 10 15 0.013 0.990 9 445 3.5 10 14 0.012 0.993 8 620 3.45 10 14 0.015 0.990 1 450 3.09 10 14 0.013 0.993 6 630 3.51 15 14 0.013 0.993 465 3.12 15 17 0.013 0.996 8 650 3.44 15 15 0.016 0.991 7 475 2.83 10 12 0.013 0.997 2 660 3.37 10 12 0.015 0.998 500 3.13 15 12 0.011 0.995 8 675 3.28 10 11 0.014 0.996 3 520 3.51 15 13 0.013 0.985 4 690 2.81 10 6 0.007 4 0.988 7 535 3.95 15 16 0.012 0.993 1 700 3.33 25 25 0.007 8 0.991 5 550 3.31 15 17 0.011 0.994 7 730 3.59 30 24 0.007 8 0.987 4 565 4.47 15 12 0.006 9 0.998 750 5.64 25 20 0.002 2 0.987 5 表 3 红椎的准饱和点光强实测值{I_{sat}} 与理论计算值{I'_{sat}} 的比较({p'_n} = 0.001 )Table 3. Measured{I_{sat}} and theoretical{I'_{sat}} at quasi-saturation point of C. hystrix({p'_n} = 0.001 )波长
Wavelength λ/nmIsat/
(μmol·m−2·s−1)Isat’/
(μmol·m−2·s−1)波长
Wavelength λ/nmIsat/
(μmol·m−2·s−1)Isat’/
(μmol·m−2·s−1)420 —— 273 575 —— 347 435 300 311 600 300 313 445 300 316 620 200 284 450 300 293 630 200 301 465 300 299 650 200 264 475 200 287 660 400 276 500 —— 328 675 400 287 520 —— 301 690 500 416 535 —— 338 700 500 441 550 300 342 730 —— 449 565 —— 506 750 —— 1154 2.3 有拐点植物光合响应曲线的饱和点光强
表4显示:相邻波长对应的补偿响应系数比较接近,整个有效辐射波段最大的差值为32×10−3,修正系数也是如此,整个有效辐射波段的最大差值为2.09×10−3,补偿比率系数的差别比较大,这是因为它是补偿响应系数与最大的光合作用速率的比值,由于各波长对应的最大光合作用速率的值差别很大所造成的。从表4还可看出各波长对应的补偿点光强的实测值与模型理论计算值吻合的比较好,相关系数的平方值接近0.99,说明该模型能够准确描述补偿点随波长变化的规律。
表 4 马褂木的补偿响应系数、修正系数、补偿比率系数和补偿点光强Table 4. Compensation response coefficient, correction coefficient, compensation ratio coefficient, and light intensity at compensation point of L. chinense波长
Wavelength
λ/nmηλ
(10−3)βλ
(10−3)γλ
(10−3)Ic/
(μmol·m−2·s−1)Ic’/
(μmol·m−2·s−1)R2 波长
Wavelength
λ/nmηλ
(10−3)βλ
(10−3)γλ
(10−3)Ic/
(μmol·m−2·s−1)Ic’/
(μmol·m−2·s−1)R2 420 37 1.3 1.8 20 21 0.982 7 575 42 0.86 2 10 17 0.991 1 435 44 0.52 3.9 20 20 0.991 6 600 52 0.62 3.3 10 16 0.990 4 445 61 0.48 6.2 10 15 0.992 7 620 57 0.2 6.2 20 18 0.996 2 450 56 0.43 7.5 15 17 0.988 630 66 0.38 5.3 10 16 0.984 5 465 39 1.1 1.1 20 18 0.992 9 650 55 0.61 3.1 15 16 0.994 7 475 62 0.47 5.6 10 18 0.991 6 660 57 0.9 2.6 10 15 0.995 3 500 46 1.8 9.1×10−8 10 18 0.994 9 675 67 0.31 6.5 10 14 0.996 5 520 36 1.8 8.4×10−7 20 21 0.988 8 690 36 0.35 2.5 10 11 0.994 4 535 41 2 0.78 15 14 0.995 700 36 0.33 2.9 25 33 0.993 550 34 1.4 0.28 20 19 0.998 3 730 25 0.92 1.7×10−7 15 25 0.996 8 565 38 2.4 6.3×10−7 20 16 0.991 8 750 35 1.8 0.99 20 18 0.999 1 表5为马褂木饱和点光强的实测值
{I_{sat}} 与理论计算值{I'_{sat}} 的比较。由表5可知:饱和点光强的实测值与模型理论计算值的一致性程度比较高,说明式(5)能够准确地描述饱和点光强随波长变化的规律,即式(2)理论模型是成功的。表 5 马褂木饱和点光强的实测值与理论计算值比较({p'_n} = 0 )Table 5. Measured and theoretical values of light intensity at saturation point of L. chinense ({p'_n} = 0 )波长
Wavelength λ/nmIsat/
(μmol·m−2·s−1)Isat’/
(μmol·m−2·s−1)波长
Wavelength λ/nmIsat/
(μmol·m−2·s−1)Isat’/
(μmol·m−2·s−1)420 400 324 575 —— 429 435 600 521 600 500 478 445 600 465 620 800 808 450 500 472 630 700 573 465 300 385 650 500 494 475 500 491 660 400 390 500 300 281 675 600 594 520 —— 292 690 700 755 535 —— 249 700 900 780 550 —— 352 730 —— 553 565 —— 218 750 —— 252 2.4 饱和点曲线和补偿点曲线
根据表1~4所确定的饱和点光强
{I_{sat}} 以及试验获取的补偿点光强{I_c} ,可以分别做出红椎与马褂木的实测饱和点曲线和补偿点曲线,如图4所示。图中饱和点曲线以饱和点光强实测值为主,对于因单色光光强的因素无法测试到饱和点光强的6个点,采用饱和点光强的理论计算值,所绘出的饱和点曲线,不影响获取该植物特征光谱的试验验证。2.5 获取植物特征光谱的方法
每一种植物均可采用植物特征光谱试验方法,获取图4相似的饱和点曲线和补偿点曲线,这是获取植物特征光谱的第一步,以此为依据,在饱和点曲线和补偿点曲线之间选择类似于饱和点曲线的多种光谱(如30种)进行栽培试验,通过比较选择长势最好和产量最高的前5种光谱所栽培的植物实体,进行有效成分的测试与分析,综合评判确定的最佳光谱为植物特征光谱。
3. 讨论与结论
采用植物特征光谱的试验方法,应用植物特征光谱测试系统进行测试,获取了红椎和马褂木两种树苗对22种波长的一系列光强对应的光合作用速率数据,采用CFT工具箱对这些数据进行曲线拟合,分别获得两个树苗的22条单色光的光合响应曲线,相关系数的平方值R2都在0.99以上,说明拟合度很好。同时采用TRM方法求解式(1)和式(2)两个模型对应的参数,进而获得渐近线型光合响应曲线树苗(红椎)的准饱和点光强和有拐点的光合响应曲线(马褂木)饱和点光强的理论计算值,绝大部分的数据比较接近,说明公式(1)可以作为描述渐近线型光合响应曲线类植物光合作用速率的一个重要数学模型。而从表4的数据看,饱和点光强的实测值与理论计算值的吻合程度没有表2的数据那么好,但是,大部分的数据比较接近,说明式(2)是描述有拐点的光合响应曲线类植物光合作用速率的一个比较好的数学模型。补偿点光强的实测值与模拟理论值很接近,而且随波长的变化也不大,由此可以得出结论:植物的补偿点光强与波长的相关性很小,说明在光合有效辐射的波长范围内,任何波长的光只要光强达到该植物的补偿点光强,则该植物的光合作用就能够启动。准饱和点光强对应的光合速率
{P_{nB}} 随波长的变化也很小,响应系数{\alpha _\lambda } 随波长的不同而变,是波长的函数。补偿响应系数{\eta _\lambda } 、修正系数{\beta _\lambda } 、补偿比率系数{\gamma _\lambda } 均随波长的变化而不同,体现出不同植物光合作用规律的特点。综上所述,通过植物特征光谱试验方法及其普通光植物特征光谱测试系统可以获得{P_{nB}} 、{\alpha _\lambda } 、{\eta _\lambda } 、{\beta _\lambda } 、{\gamma _\lambda } 等与光合特性有关又能体现植物自身特点的参数,说明该方法及其系统是研究植物光合特性的普遍方法及通用装备。从试验数据拟合出的单色光植物光合响应曲线中,发现了树苗的光合响应曲线可以分成两类的现象,可能对树种的分类有帮助,其机理有待进一步研究。致谢:感谢南平市森科种苗有限公司薛华同志的支持和帮助。
-
表 1 人参属物种叶绿体基因组特征
Table 1 Chloroplast genomes of Panax species
物种 Species 登录号 Accession 长度 Length/bp GC含量 GC content/% 基因组 Genome LSC SSC IR 基因组 Genome LSC SSC IR 人参
P. ginsengMK408938 156 241 86 128 18 077 26 018 38.07 36.27 32.18 43.10 竹节参
P. japonicusKP036469 156 188 86 199 18 013 25 988 38.07 36.29 32.21 43.05 疙瘩七
P. japonicus var. bipinnatifidusMK408962 156 244 86 186 18 006 26 026 38.06 36.28 32.21 43.05 珠子参
P. majorMN496312 156 402 86 189 18 007 26 103 38.07 36.28 32.22 43.05 三七
P. notoginsengMK408945 156 319 86 157 18 004 26 079 38.07 36.27 32.23 43.07 假人参
P. pseudoginsengMW145449 156 074 86 124 18 008 25 971 38.06 36.27 32.22 43.06 西洋参
P. quinquefoliusMK408953 156 070 86 077 17 993 26 000 38.07 36.27 32.22 43.08 屏边三七
P. stipuleanatusMK408965 156 007 86 083 18 150 25 887 38.04 36.26 32.07 43.08 三叶参
P. trifoliusMF100782 156 157 86 322 18 047 25 894 38.08 36.27 32.26 43.10 越南参
P. vietnamensisKP036470 155 993 86 178 17 935 25 940 38.05 36.27 32.22 43.02 野三七
P. vietnamensis var. fuscidiscusMT798585 156 284 86 171 17 971 26 071 38.06 36.28 32.27 42.98 越南参变种
P. vietnamensis var. langbianensisMT798583 155 984 86 174 17 934 25 938 38.05 36.27 32.21 43.03 峨眉三七
P. wangianusMK408964 156 189 86 190 17 969 26 015 38.07 36.29 32.21 43.03 姜状三七
P. zingiberensisMK408969 156 192 86 116 17 966 26 055 38.07 36.31 32.27 43.00 表 2 基于位点模型的人参属叶绿体蛋白编码基因正选择分析
Table 2 Potential positive selection test on chloroplast genomes of Panax species based on site model
基因
GeneM1 vs M2 M7 vs M8 正选择位点
Positively selected sites2ΔlnL P 值 2ΔlnL P 值 accD 2.39 3.0×10-1 2.47 2.9×10-1 139 Q, 0.505; 141 S, 0.517; 144 R, 0.528; 209 A, 0.856; 264 E, 0.504; 364 I, 0.507; 372 L, 0.526; 494 Q, 0.855 atpF 3.37 1.9×10-1 3.42 1.8×10-1 14 W, 0.541; 77 Q, 0.537; 158 G, 0.854 cemA 7.48 2.4×10-2 7.48 2.4×10-2 14 A, 0.547; 22 W, 0.538; 54 K, 0.536; 102 D, 0.533; 103 R, 0.544; 118 V, 0.973*; 205 F, 0.531; 208 W, 0.540 clpP 13.40 1.2×10-3 13.40 1.2×10-3 38 P, 0.990*; 39 V, 0.622; 40 A, 0.935; 41 S, 0.612; 138 G, 0.629; 160 Q, 0.626 matK 19.22 6.7×10-5 19.25 6.6×10-5 83 R, 0.994** ndhD 0.77 6.8×10-1 0.99 6.1×10-1 22 F, 0.830; 41 I, 0.829 ndhF 19.81 5.0×10-5 20.29 3.9×10-5 17 P, 0.597; 123 R, 0.592; 133 T, 0.583; 223 F, 0.959*; 421 I, 0.964*; 463 Q, 0.591; 483 R, 0.569; 499 H, 0.552; 502 A, 0.965*; 511 Q, 0.591; 512 M, 0.526; 532 K, 0.559; 534 I, 0.962*; 572 N, 0.543; 598 D, 0.573; 606 F, 0.548; 621 L, 0.630; 641 I, 0.549; 642 G, 0.597; 657 G, 0.594; 670 A, 0.963*; 677 I, 0.532; 703 S, 0.548; 737 Y, 0.964*; 743 F, 0.537; 745 D, 0.958*; 746 L, 1.000** rbcL 15.94 3.5×10-4 16.10 3.2×10-4 225 L, 0.632; 226 Y, 0.969*; 230 A, 0.635; 270 I, 0.866; 281 A, 0.644; 328 S, 0.999**; 353 F, 0.584; 439 R, 0.968*; 472 V, 0.609; 474 I, 0.621; 478 V, 0.714 rpl2 0 1 0 1 99 K, 0.651; 265 D, 0.650; 267 L, 0.652 rpoA 14.53 7.0×10-4 14.56 6.9×10-4 34 L, 0.506; 256 E, 0.991**; 264 N, 0.999**; 272 S, 0.533; 320 A, 0.849; 337 L, 0.537 rpoB 0.90 6.4×10-1 0.90 6.4×10-1 37 Y, 0.505; 165 R, 0.502; 193 R, 0.502; 295 I, 0.502; 365 T, 0.504; 459 S, 0.502; 695 K, 0.506; 735 T, 0.504; 935 D, 0.832 rpoC2 0 1 0 1 85 Q, 0.736 rps12 0 1 0.02 9.9×10-1 25 R, 0.502 ycf1 75.81 0 75.78 0 205 K, 1.000**; 206 Y, 0.999**; 288 I, 0.876; 532 T, 0.803; 533 K, 0.956*; 604 R, 0.801; 624 R, 0.825; 668 I, 0.779; 685 W, 0.837; 699 V, 0.824; 701 Q, 0.802; 706 V, 0.814; 720 V, 0.856; 722 T, 0.814; 723 D, 0.808; 731 R, 0.982*; 733 K, 0.809; 734 I, 0.996**; 736 L, 0.855; 737 I, 0.803; 751 N, 0.760; 809 A, 0.828; 906 N, 0.768; 914 E, 0.795; 915 L, 0.788; 929 L, 0.835; 957 K, 0.979*; 1030 E, 0.979*; 1116 K, 0.804; 1172 M, 0.731; 1186 Q, 0.798; 1206 Q, 0.798; 1502 Q, 0.754; 1601 I, 1.000**; 1616 D, 0.773;
1631 R, 0.814; 1652 R, 0.795; 1665 V, 0.738; 1685 E, 0.978*; 1710 I, 0.979*; 1739 R, 0.793ycf2 14.95 5.7×10-4 14.95 5.7×10-4 30 R, 0.767; 56 R, 0.770; 103 L, 0.765; 126 S, 0.765; 248 S, 0.769; 285 L, 0.765; 320 L, 0.766; 332 R, 0.764; 432 V, 0.765; 466 Y, 0.765; 573 I, 0.764; 598 L, 0.766; 606 M, 0.759; 759 F, 0.763; 850 S, 0.965*; 943 E, 0.764; 1037 P, 0.766; 1050 F, 0.765; 1064 S, 0.763; 1074 I, 0.965*; 1093 L, 0.767; 1101 N, 0.964*; 1171 W, 0.767; 1173 S, 0.765; 1245 P, 0.766; 1265 D, 0.762; 1309 D, 0.762; 1344 R, 0.764; 1422 L, 0.965*; 1431 F, 0.763; 1522 D, 0.767; 1550 Y, 0.766; 1656 N, 0.763; 1717 S, 0.768; 1799 R, 0.769; 2045 E, 0.764 *和**分别表示正选择位点概率> 95%和> 99%。
* and ** indicate that the probabilities of positively selected sites are > 95% and > 99%, respectively. -
[1] WEN J, ZIMMER E A. Phylogeny and biogeography of Panax L. (The Ginseng Genus, Araliaceae): Inferences from ITS sequences of nuclear ribosomal DNA [J]. Molecular Phylogenetics and Evolution, 1996, 6(2): 167−177. DOI: 10.1006/mpev.1996.0069
[2] NGUYEN V B, LINH GIANG V N, WAMINAL N E, et al. Comprehensive comparative analysis of chloroplast genomes from seven Panax species and development of an authentication system based on species-unique single nucleotide polymorphism markers [J]. Journal of Ginseng Research, 2020, 41(1): 135−144. DOI: 10.1016/j.jgr.2018.06.003
[3] NGUYEN V B, PARK H S, LEE S C, et al. Authentication markers for five major Panax species developed via comparative analysis of complete chloroplast genome sequences [J]. Journal of Agricultural and Food Chemistry, 2017, 65(30): 6298−6306. DOI: 10.1021/acs.jafc.7b00925
[4] LI S, WANG P, YANG W Z, et al. Characterization of the components and pharmacological effects of mountain-cultivated ginseng and garden ginseng based on the integrative pharmacology strategy [J]. Frontiers in Pharmacology, 2021, 12: 659954. DOI: 10.3389/fphar.2021.659954
[5] SHIN B K, KWON S W, PARK J H. Chemical diversity of ginseng saponins from Panax ginseng [J]. Journal of Ginseng Research, 2015, 39(4): 287−298. DOI: 10.1016/j.jgr.2014.12.005
[6] CHOI H I, KIM N H, LEE J, et al. Evolutionary relationship of Panax ginseng and P. quinquefolius inferred from sequencing and comparative analysis of expressed sequence tags [J]. Genetic Resources and Crop Evolution, 2013, 60(4): 1377−1387. DOI: 10.1007/s10722-012-9926-3
[7] YANG Z J, LIU G Z, ZHANG G H, et al. The chromosome-scale high-quality genome assembly of Panax notoginseng provides insight into dencichine biosynthesis [J]. Plant Biotechnology Journal, 2021, 19(5): 869−871. DOI: 10.1111/pbi.13558
[8] LIU C, CHEN H H, TANG L Z, et al. Plastid genome evolution of a monophyletic group in the subtribe Lauriineae (Laureae, Lauraceae)[J]. Plant Diversity, 2022 , 44 (4): 377 − 388. DOI: 10.1016/j.pld.2021.11.009
[9] 刘潮, 韩利红, 彭悦, 等. 黄丹木姜子叶绿体基因组特征分析 [J]. 南方农业学报, 2022, 53(1):12−20. DOI: 10.3969/j.issn.2095-1191.2022.01.002 LIU C, HAN L H, PENG Y, et al. Characteristics of chloroplast genome of Litsea elongata(wall. ex nees)Benth. et hook. F [J]. Journal of Southern Agriculture, 2022, 53(1): 12−20.(in Chinese) DOI: 10.3969/j.issn.2095-1191.2022.01.002
[10] 刘潮, 韩利红, 代小波, 等. 辣椒属叶绿体基因组特征及进化 [J]. 热带作物学报, 2022, 43(3):447−454. DOI: 10.3969/j.issn.1000-2561.2022.03.002 LIU C, HAN L H, DAI X B, et al. Characteristics and phylogenetics of the complete chloroplast genomes of Capsicum species [J]. Chinese Journal of Tropical Crops, 2022, 43(3): 447−454.(in Chinese) DOI: 10.3969/j.issn.1000-2561.2022.03.002
[11] YU J J, FU J, FANG Y P, et al. Complete chloroplast genomes of Rubus species (Rosaceae) and comparative analysis within the genus [J]. BMC Genomics, 2022, 23(1): 32. DOI: 10.1186/s12864-021-08225-6
[12] 刘潮, 唐利洲, 韩利红. 四川山胡椒叶绿体基因组特征及山胡椒属系统发育 [J]. 林业科学, 2021, 57(12):167−174. LIU C, TANG L Z, HAN L H. Characterization of the chloroplast genome of Lindera setchuenensis and phylogenetics of the genus Lindera [J]. Scientia Silvae Sinicae, 2021, 57(12): 167−174.(in Chinese)
[13] LIU C K, YANG Z Y, YANG L F, et al. The complete plastome of Panax stipuleanatus: Comparative and phylogenetic analyses of the genus Panax (Araliaceae) [J]. Plant Diversity, 2018, 40(6): 265−276. DOI: 10.1016/j.pld.2018.11.001
[14] 岳杰. 四种人参属植物叶绿体全基因组结构及其系统发育分析[D]. 汉中: 陕西理工大学, 2021 YUE J. Structural and phylogenetic analyses of the complete chloroplast genomes of four species within the genus Panax linn. [D]. Hanzhong: Shaanxi University of Technology, 2021. (in Chinese)
[15] 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. DOI: 10.1093/nar/29.22.4633
[16] BEIER S, THIEL T, MÜNCH T, et al. MISA-web: A web server for microsatellite prediction [J]. Bioinformatics, 2017, 33(16): 2583−2585. DOI: 10.1093/bioinformatics/btx198
[17] AMIRYOUSEFI A, HYVÖNEN J, POCZAI P. IRscope: An online program to visualize the junction sites of chloroplast genomes [J]. Bioinformatics (Oxford, England), 2018, 34(17): 3030−3031. DOI: 10.1093/bioinformatics/bty220
[18] FRAZER K A, PACHTER L, POLIAKOV A, et al. VISTA: Computational tools for comparative genomics[J]. Nucleic Acids Research, 2004, 32(Web Server issue): W273-W279.
[19] KATOH K, ROZEWICKI J, YAMADA K D. MAFFT online service: Multiple sequence alignment, interactive sequence choice and visualization [J]. Briefings in Bioinformatics, 2019, 20(4): 1160−1166. DOI: 10.1093/bib/bbx108
[20] ROZAS J, FERRER-MATA A, SÁNCHEZ-DELBARRIO J C, et al. DnaSP 6: DNA sequence polymorphism analysis of large data sets [J]. Molecular Biology and Evolution, 2017, 34(12): 3299−3302. DOI: 10.1093/molbev/msx248
[21] XU B, YANG Z H. PAMLX: A graphical user interface for PAML [J]. Molecular Biology and Evolution, 2013, 30(12): 2723−2724. DOI: 10.1093/molbev/mst179
[22] MINH B Q, SCHMIDT H A, CHERNOMOR O, et al. IQ-TREE 2: New models and efficient methods for phylogenetic inference in the genomic era [J]. Molecular Biology and Evolution, 2020, 37(5): 1530−1534. DOI: 10.1093/molbev/msaa015
[23] GANLEY A R D, KOBAYASHI T. Monitoring the rate and dynamics of concerted evolution in the ribosomal DNA repeats of Saccharomyces cerevisiae using experimental evolution [J]. Molecular Biology and Evolution, 2011, 28(10): 2883−2891. DOI: 10.1093/molbev/msr117
[24] ZEB U, DONG W L, ZHANG T T, et al. Comparative plastid genomics of Pinus species: Insights into sequence variations and phylogenetic relationships [J]. Journal of Systematics and Evolution, 2020, 58(2): 118−132. DOI: 10.1111/jse.12492
[25] ZHU B, QIAN F, HOU Y F, et al. Complete chloroplast genome features and phylogenetic analysis of Eruca sativa (Brassicaceae) [J]. PLoS One, 2021, 16(3): e0248556. DOI: 10.1371/journal.pone.0248556
[26] LI D M, YE Y J, XU Y C, et al. Complete chloroplast genomes of Zingiber montanum and Zingiber zerumbet: Genome structure, comparative and phylogenetic analyses [J]. PLoS One, 2020, 15(7): e0236590. DOI: 10.1371/journal.pone.0236590
[27] XU K W, LIN C X, LEE S Y, et al. Comparative analysis of complete Ilex (Aquifoliaceae) chloroplast genomes: Insights into evolutionary dynamics and phylogenetic relationships [J]. BMC Genomics, 2022, 23(1): 203. DOI: 10.1186/s12864-022-08397-9
[28] KASHI Y, KING D G. Simple sequence repeats as advantageous mutators in evolution [J]. Trends in Genetics, 2006, 22(5): 253−259. DOI: 10.1016/j.tig.2006.03.005
[29] HAN C Y, DING R, ZONG X Y, et al. Structural characterization of Platanthera ussuriensis chloroplast genome and comparative analyses with other species of Orchidaceae [J]. BMC Genomics, 2022, 23(1): 84. DOI: 10.1186/s12864-022-08319-9
[30] YAN J W, LI J H, YU L, et al. Comparative chloroplast genomes of Prunus subgenus Cerasus (Rosaceae): Insights into sequence variations and phylogenetic relationships [J]. Tree Genetics & Genomes, 2021, 17(6): 50.
[31] 马孟莉, 张薇, 孟衡玲, 等. 草果叶绿体基因组特征及系统发育分析 [J]. 中草药, 2021, 52(19):6023−6031. DOI: 10.7501/j.issn.0253-2670.2021.19.025 MA M L, ZHANG W, MENG H L, et al. Characterization and phylogenetic analysis of the complete chloroplast genome of Amomum tsao-ko [J]. Chinese Traditional and Herbal Drugs, 2021, 52(19): 6023−6031.(in Chinese) DOI: 10.7501/j.issn.0253-2670.2021.19.025
[32] KIM K J. Complete chloroplast genome sequences from Korean ginseng (Panax schinseng nees) and comparative analysis of sequence evolution among 17 vascular plants [J]. DNA Research, 2004, 11(4): 247−261. DOI: 10.1093/dnares/11.4.247
[33] LI R, MA P F, WEN J, et al. Complete sequencing of five Araliaceae chloroplast genomes and the phylogenetic implications [J]. PLoS One, 2013, 8(10): e78568. DOI: 10.1371/journal.pone.0078568
[34] TROFIMOV D, CADAR D, SCHMIDT-CHANASIT J, et al. A comparative analysis of complete chloroplast genomes of seven Ocotea species (Lauraceae) confirms low sequence divergence within the Ocotea complex [J]. Scientific Reports, 2022, 12: 1120. DOI: 10.1038/s41598-021-04635-4
[35] YANG J, CHIANG Y C, HSU T W, et al. Characterization and comparative analysis among plastome sequences of eight endemic Rubus (Rosaceae) species in Taiwan [J]. Scientific Reports, 2021, 11: 1152. DOI: 10.1038/s41598-020-80143-1
[36] NIU Y F, GAO C W, LIU J. Comparative analysis of the complete plastid genomes of Mangifera species and gene transfer between plastid and mitochondrial genomes [J]. PeerJ, 2021, 9: e10774. DOI: 10.7717/peerj.10774
[37] GARUD N R, POLLARD K S. Population genetics in the human microbiome [J]. Trends in Genetics:TIG, 2020, 36(1): 53−67. DOI: 10.1016/j.tig.2019.10.010
[38] JI Y H, LIU C K, YANG Z Y, et al. Testing and using complete plastomes and ribosomal DNA sequences as the next generation DNA barcodes in Panax (Araliaceae) [J]. Molecular Ecology Resources, 2019, 19(5): 1333−1345. DOI: 10.1111/1755-0998.13050
[39] KIM N H, CHOI H I, KIM K H, et al. Evidence of genome duplication revealed by sequence analysis of multi-loci expressed sequence tag-simple sequence repeat bands in Panax ginseng Meyer [J]. Journal of Ginseng Research, 2014, 38(2): 130−135. DOI: 10.1016/j.jgr.2013.12.005
-
期刊类型引用(0)
其他类型引用(2)