Processing math: 20%
  • 中文核心期刊
  • CSCD来源期刊
  • 中国科技核心期刊
  • CA、CABI、ZR收录期刊

人参属植物叶绿体基因组特征及其进化的研究

刘潮, 李敏, 任怡园, 钱柏霖, 韩利红

刘潮,李敏,任怡园,等. 人参属植物叶绿体基因组特征及其进化的研究 [J]. 福建农业学报,2022,37(7):886−896. DOI: 10.19303/j.issn.1008-0384.2022.007.009
引用本文: 刘潮,李敏,任怡园,等. 人参属植物叶绿体基因组特征及其进化的研究 [J]. 福建农业学报,2022,37(7):886−896. DOI: 10.19303/j.issn.1008-0384.2022.007.009
LIU C, LI M, REN Y Y, et al. Characteristics and Evolution of Panax Chloroplast Genomes [J]. Fujian Journal of Agricultural Sciences,2022,37(7):886−896. DOI: 10.19303/j.issn.1008-0384.2022.007.009
Citation: LIU C, LI M, REN Y Y, et al. Characteristics and Evolution of Panax Chloroplast Genomes [J]. Fujian Journal of Agricultural Sciences,2022,37(7):886−896. DOI: 10.19303/j.issn.1008-0384.2022.007.009

人参属植物叶绿体基因组特征及其进化的研究

基金项目: 国家自然科学基金项目(32060710、32100010)
详细信息
    作者简介:

    刘 潮(1980-),男,博士,副教授,研究方向:植物系统发育与进化(E-mail:liuchao_80@163.com

    通讯作者:

    韩利红(1981-),女,博士,副教授,研究方向:真菌系统发育与进化(E-mail:hanlihong9527@126.com

  • 中图分类号: R 282

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比率,发现功能未知基因clpPycf1ycf2受正选择作用。系统发育分析显示,屏边三七和三叶参位于基部支系,四倍体人参和西洋参与其他二倍体物种聚在不同支系,三七、竹节参和越南参则亲缘关系较近。
      结论  人参属叶绿体基因组基因数目和顺序一致,基因组结构保守,重复序列数目和类型存在差异,单拷贝区核苷酸多态性高于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.
  • 【研究意义】随着科技的进步,现代农业产业不断往高度集约化与高效种植方向发展,现代化农业设施得到广泛应用。我国设施园艺栽培面积高居世界首位,2016年已达到476.5万hm2,产值超1.46万亿元[1]。农业设施的补光问题日益凸显,人造光源在植物光照中发挥着越来越重要的作用,试验表明LED光源相对传统光源在植物光照方面有明显的优势[1]。以纯人工光为光源的植物工厂代表了设施农业的发展趋势,植物对光吸收的效率决定了设施农业产业的经济效益,因此寻找植物生长所需的最佳复合光谱是关键。植物特征光谱代表了植物在某个生长周期内的最佳光需求[2],若能够依据植物特征光谱精准给光,将有效促进植物产量和品质的提高,从而提高设施农业产业模式的社会经济效益[3]。植物特征光谱是以植物在多种单色光条件下的光合响应曲线为基础进行定义的[2],因此,研究植物在单色光照射条件下的光合响应曲线及其数学模型,对植物光照应用领域具有重要理论价值与实践指导意义。【前人研究进展】当前,在设施农业中普遍采用LED人工光源补光或进行全人工光栽培,对植物在不同光照环境下表现出的相应特性的研究也越来越多[4-14]。植物的生长及其有效成分的积累与光的波长、强度和光周期密切相关[1]。对于以白光为基础的光合响应曲线模型的研究也有许多成果[15-22]。【本研究切入点】目前的研究基本上都采用红绿蓝三色光为光源,此类光的光谱仅为植物光合有效辐射范围的一部分,无法反映植物光合响应的全貌。仅用三色光对植物进行光照培育,可能会导致植物因为部分波长的光的缺失而出现相应的不良生长现象。同时,在对植物光合响应模型的相关研究中,有关植物在单色光条件下——即将光波长作为光合响应重要影响参数的研究有待深入进行。【拟解决的关键问题】以获取植物特征光谱——即植物在一定环境(非光)中进行光合作用最适宜的光谱、光照强度和光照周期为出发点,选择10种具有代表性的植物幼苗作为试验对象,采用自主研制的多种单色光发生器搭建光合作用(速率)测试系统[2],开展植物在多种单色光条件下的光合响应研究,根据单色光条件下光合响应饱和点实测数据,分析、建立具有较高拟合度的光合响应曲线及数学模型,以期发现植物在不同波长光照射下所呈现光合响应特性的一般规律,为最终形成植物特征光谱提供理论依据和数据支撑。

    2021年6月,样本取自南平市森科种苗有限公司苗圃,该苗圃位于延平区西芹镇坑底村,海拔60 m,湿度65%~100%,介于北纬26°15′~26°52′,东经117°50′~118°40′,属中亚热带海洋季风气候,年均气温17.3 ℃(夏季温度30~40 ℃),年降水量1 669 mm。苗圃为露天的普通林地培育树苗,光照较强。在该苗圃取一年生长势良好的10个品种树苗:红椎、马褂木、红花荷、苦槠、乐昌含笑、闽楠、木荷、刨花楠、乳源木莲和油茶。限于篇幅,选择两个具有代表性的树苗红椎和马褂木的测试数据进行分析研究。

    自主研制的多种单色光发生器,该设备由一种双驼峰型光谱的发光装置[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 ∆λ/nm
    420 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
    下载: 导出CSV 
    | 显示表格

    按光合作用测试仪使用的常规要求,选择树苗上第3~4个发育正常的叶片,置于叶室里,采用植物特征光谱试验方法[2]进行测试,即任选多种单色光发生器产生的一种单色光,让其均匀地照射在光合作用测试系统的叶室上,通过光合作用测试系统就可以获得被测试的植物在该单色光此刻对应的光强下的光合作用速率。为了避免杂光的干扰,从出光口到叶室,采取遮光措施。按由小到大的顺序调整该单色光的光合有效辐射强度PAR(以下简称光强I/µmol·m−2·s−1,分别为10,20,30,40,60,80,100,150,200,300,400,……,最大值),获得该单色光在不同光强条件下的光合作用速率,即光合作用速率与光强的关系曲线,通常称之为单色光植物光合作用响应曲线(以下简称光合响应曲线)。

    应用Matlab 2012b软件中的Curve Fitting Tool曲线拟合工具箱(以下简称CFT工具)对所获得的光合作用速率及其对应的光强数据进行曲线拟合,获得每一种单色光的光合响应曲线,该工具箱采用的是非线性最小二乘法的信任区域理论(Trust Region Methods)求解模型参数(以下简称TRM方法)。

    根据有两类不同的单色光光合响应曲线的试验结果,考虑到不同波长的光合响应曲线的差异,分别采用不同的模型来描述。

    对于渐近线型植物单色光光合响应曲线采用以下模型来描述。

    Pn=PnB(1 - e - αλ(I - Ic)) (1)

    式中Pn为净光合速率,I为光强,Ic为补偿点对应的光强,对于特定的树种和特定的波长Ic为常数,PnB为准饱和点光强对应的光合速率,αλ为响应系数。由式(1)可知,当I=Ic时,Pn=0,这与试验测试结果相符。

    对于有拐点的植物单色光光合响应曲线,采用以下模型来描述。

    Pn=ηλ1βλI1+γλI(IIc) (2)

    式中Pn为净光合速率,I为光强,Ic为补偿点对应的光强,ηλ为补偿响应系数,βλ为修正系数,γλ=ηλ/Pmax为补偿比率系数,其中 {P_{\max }} 为植物的最大光合作用速率。由式(2)可知,当 I = {I_c} 时, {P_n} = 0 ,同样与试验测试结果相符。

    对于渐近线型植物光合响应曲线不存在光合作用速率最大值 {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}} 为准饱和点光强。

    对应有拐点光合响应曲线的饱和点光强的确定,就是将式(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}}

    选择22种单色光,采用植物特征光谱的试验方法,进行试验测试,获得的数据应用CFT工具进行模拟,分别获得了10种树苗的22种波长的光合响应曲线。

    试验结果表明,树苗的光合响应曲线有两类:一类是渐近线型光合响应曲线,另一类是有拐点的光合响应曲线。具有代表性的树种分别为红椎和马褂木。红椎的光合响应曲线是渐近线型曲线,这类树种的光合响应曲线没有极值;马褂木的光合响应曲线则是有拐点的曲线,这类曲线有极值。由此可见,同一树苗对不同波长的光合响应曲线是不同的,不同树苗对同一波长的光合响应曲线也明显不同。

    图13还可以看出光强在20~100 µmol·m−2·s−1光合响应曲线接近是线性的,超过100 µmol·m−2·s−1以后光合响应曲线的特异性表现得比较明显,说明植物光合响应的差别在100 µmol·m−2·s−1以后表现得更为显著。

    图  1  红椎的光合响应曲线
    Figure  1.  LRC of C. hystrix
    图  2  马褂木的光合响应曲线
    Figure  2.  LRC of L. chinense
    图  3  红椎和马褂木对3种不同波长的光合响应曲线
    A、B、C分别表示波长为435、600、700 nm。
    Figure  3.  LRCs of C. hystrix and L. chinense exposed to lights of 3 different wavelengths
    A, B and C represent wavelengths of 435 nm, 600 nm and 700 nm respectively.

    表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
    λ/nm
    PnbIc/
    (μmol·m−2·s−1
    Ic′/
    (μmol·m−2·s−1
    αλ/
    (m2·s−1·μmol−1
    R2波长
    Wavelength
    λ/nm
    PnbIc/
    (μ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
    下载: 导出CSV 
    | 显示表格
    表  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 λ/nm
    Isat/
    (μmol·m−2·s−1
    Isat’/
    (μmol·m−2·s−1
    波长
    Wavelength λ/nm
    Isat/
    (μ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
    下载: 导出CSV 
    | 显示表格

    表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
    下载: 导出CSV 
    | 显示表格

    表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 λ/nm
    Isat/
    (μmol·m−2·s−1
    Isat’/
    (μmol·m−2·s−1
    波长
    Wavelength λ/nm
    Isat/
    (μ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
    下载: 导出CSV 
    | 显示表格

    根据表14所确定的饱和点光强 {I_{sat}} 以及试验获取的补偿点光强 {I_c} ,可以分别做出红椎与马褂木的实测饱和点曲线和补偿点曲线,如图4所示。图中饱和点曲线以饱和点光强实测值为主,对于因单色光光强的因素无法测试到饱和点光强的6个点,采用饱和点光强的理论计算值,所绘出的饱和点曲线,不影响获取该植物特征光谱的试验验证。

    图  4  红椎与马褂木的饱和点曲线和补偿点曲线
    Figure  4.  Light saturation and compensation curves of C. hystrix and L. chinense

    每一种植物均可采用植物特征光谱试验方法,获取图4相似的饱和点曲线和补偿点曲线,这是获取植物特征光谱的第一步,以此为依据,在饱和点曲线和补偿点曲线之间选择类似于饱和点曲线的多种光谱(如30种)进行栽培试验,通过比较选择长势最好和产量最高的前5种光谱所栽培的植物实体,进行有效成分的测试与分析,综合评判确定的最佳光谱为植物特征光谱。

    采用植物特征光谱的试验方法,应用植物特征光谱测试系统进行测试,获取了红椎和马褂木两种树苗对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   人参属物种叶绿体基因组长重复序列类型及分布

    A:正向、回文、反向、互补重复数量;B:长重复类型及数量。

    Figure  1.   Types and distribution of long repeats in chloroplast genomes of Panax species

    A: Counts of forward, palindromic, reverse and complementary repeats; B: types and counts of long repeats.

    图  2   人参属物种叶绿体基因组SSR位点类型及分布

    A:单核苷酸、二核苷酸、三核苷酸、四核苷酸、五核苷酸、六核苷酸重复数量;B:SSR类型及数量。

    Figure  2.   Types and distribution of SSRs in chloroplast genomes of Panax species

    A: Counts of mono-, di-, tri-, tetra-, penta- and hexanucleotides; B: types and counts of SSR.

    图  3   人参属物种叶绿体基因组LSC、SSC和IRs边界比较

    Figure  3.   Borders of LSC, SSC, and IRs in chloroplast genomes of Panax species

    图  4   人参属物种叶绿体基因组序列比较

    Figure  4.   Sequences of chloroplast genomes of Panax species

    图  5   人参属物种叶绿体基因组序列多态性分析

    Figure  5.   Nucleotide diversity of chloroplast genomes of Panax species

    图  6   人参属叶绿体蛋白编码基因dN/dS分析

    Figure  6.   dN/dS ratios on protein coding genes of chloroplasts of Panax species

    图  7   基于IQ-TREE软件构建的人参属物种叶绿体基因组系统发育树

    Figure  7.   Phylogenetic tree of chloroplast genes of Panax species based on IQ-TREE

    表  1   人参属物种叶绿体基因组特征

    Table  1   Chloroplast genomes of Panax species

    物种 Species登录号 Accession长度 Length/bpGC含量 GC content/%
    基因组 GenomeLSCSSCIR基因组 GenomeLSCSSCIR
    人参
    P. ginseng
    MK408938 156 241 86 128 18 077 26 018 38.07 36.27 32.18 43.10
    竹节参
    P. japonicus
    KP036469 156 188 86 199 18 013 25 988 38.07 36.29 32.21 43.05
    疙瘩七
    P. japonicus var. bipinnatifidus
    MK408962 156 244 86 186 18 006 26 026 38.06 36.28 32.21 43.05
    珠子参
    P. major
    MN496312 156 402 86 189 18 007 26 103 38.07 36.28 32.22 43.05
    三七
    P. notoginseng
    MK408945 156 319 86 157 18 004 26 079 38.07 36.27 32.23 43.07
    假人参
    P. pseudoginseng
    MW145449 156 074 86 124 18 008 25 971 38.06 36.27 32.22 43.06
    西洋参
    P. quinquefolius
    MK408953 156 070 86 077 17 993 26 000 38.07 36.27 32.22 43.08
    屏边三七
    P. stipuleanatus
    MK408965 156 007 86 083 18 150 25 887 38.04 36.26 32.07 43.08
    三叶参
    P. trifolius
    MF100782 156 157 86 322 18 047 25 894 38.08 36.27 32.26 43.10
    越南参
    P. vietnamensis
    KP036470 155 993 86 178 17 935 25 940 38.05 36.27 32.22 43.02
    野三七
    P. vietnamensis var. fuscidiscus
    MT798585 156 284 86 171 17 971 26 071 38.06 36.28 32.27 42.98
    越南参变种
    P. vietnamensis var. langbianensis
    MT798583 155 984 86 174 17 934 25 938 38.05 36.27 32.21 43.03
    峨眉三七
    P. wangianus
    MK408964 156 189 86 190 17 969 26 015 38.07 36.29 32.21 43.03
    姜状三七
    P. zingiberensis
    MK408969 156 192 86 116 17 966 26 055 38.07 36.31 32.27 43.00
    下载: 导出CSV

    表  2   基于位点模型的人参属叶绿体蛋白编码基因正选择分析

    Table  2   Potential positive selection test on chloroplast genomes of Panax species based on site model

    基因
    Gene
    M1 vs M2M7 vs M8正选择位点
    Positively selected sites
    2ΔlnLP2ΔlnLP
    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.793
    ycf2 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.
    下载: 导出CSV
  • [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)

图(7)  /  表(2)
计量
  • 文章访问数:  587
  • HTML全文浏览量:  170
  • PDF下载量:  25
  • 被引次数: 2
出版历程
  • 收稿日期:  2022-03-30
  • 修回日期:  2022-05-29
  • 录用日期:  2022-03-30
  • 网络出版日期:  2022-08-28
  • 刊出日期:  2022-07-27

目录

/

返回文章
返回