Page 26 - 《广西植物》2026年第9期
P. 26
1 4 6 2 广 西 植 物 46 卷
RSCU)分析ꎬ所有参数均设置为默认值ꎬ最后得出 trnS ̄GCU 5 个基因均有 2 个拷贝ꎮ 丽江云杉叶绿
丽江云杉叶绿体基因组密码子使用情况ꎬRSCU> 体基因组中共鉴定出 14 个含有内含子的基因ꎬ分
1、 = 1、<1 分别代表该密码子使用频率较高、无偏 别为 atpF、petB、petD、rpoC1、rpl2、rpl16、ycf3、rps12
好性、较低ꎮ 8 个 蛋 白 编 码 基 因 和 trnA ̄UGC、 trnG ̄GCC、 trnI ̄
1.6 系统发育分析 GAU、trnK ̄UUU、 trnL ̄UAA、 trnV ̄UAC 6 个 tRNA 基
从 NCBI 下载了已公布的所有 16 种云杉属植 因ꎮ 除 ycf3 基因序列含有 2 个内含子以外ꎬ其余
物 材 料ꎬ 以 及 两 个 外 类 群 银 杉 ( Cathaya 基因均只有 1 个内含子ꎮ
argyrophyllaꎬ NC _ 014589 ) 和 白 皮 松 ( Pinus 2.2 重复序列分析
bungeanaꎬNC_028421) 进行系统发育重建ꎮ 采用 SSR 一般由 1 ~ 6 个核苷酸重复多次构成ꎬ真
PhyloSuite v1.2.2 软件(Zhang et al.ꎬ 2020) 进行分 核生物基因组中普遍存在ꎬ被广泛应用于系统发
析ꎬ提取到 72 个共有的蛋白编码基因ꎬ选择该软 育、谱系地理、遗传图谱及遗传多样性等方面的研
件 中 内 置 的 MAFFT ̄v7. 505 ( Katoh & Standleyꎬ 究ꎮ 在丽江云杉叶绿体基因组中ꎬ共检测到 29 个
2013)分别进行序列比对ꎬ比对策略选用 FFT ̄NS ̄ SSRꎬ单核苷酸重复 22 个ꎬ其中 20 个均是 A / T 类
1ꎬ之 后 串 联 为 一 个 序 列ꎬ 利 用 PartitionFinder2 型的重复ꎬ仅有 2 个是 G / C 重复ꎬ双核苷酸重复 7
(Lanfear et al.ꎬ 2017)分别选择适用于最大似然法 个ꎬ全部为 AT / TA 重复ꎬ没有检测到三、四、五等其
和贝叶斯法的最优替换模型ꎬ最大似然法采用 IQ ̄ 他类型的核苷酸重复( 图 2)ꎮ 从所在区域来看ꎬ
TREE( Nguyen et al.ꎬ 2015) 模块ꎬ自展值设置为 23 个 位 于 基 因 间 隔 区 ( intergenic spacer regionꎬ
1 000进行系统发育重建ꎬ贝叶斯法采用 MrBayes IGS)ꎬ2 个位于 rpoC1 内含子ꎬ3 个分别位于 trnG ̄
v3.2( Ronquist et al.ꎬ 2012) 进 行 计 算ꎬ马 尔 科 夫 GCC、rpl16 和 ycf3 内含子区域ꎬ1 个位于 psbZ 基因
链-蒙特卡罗( Markov Chain Monte CarloꎬMCMC) 的蛋白编码序列区( coding sequence regionꎬCDS)ꎮ
运行 1 000 000 代ꎬ每 1 000 代取样 1 次ꎬ2 线程运 这说明丽江云杉叶绿体基因组 SSR 序列更偏好
行ꎬ链设置为 4 条ꎬBurnin 比率值设置为 0.25ꎬ得 A / T 碱基ꎬ并且大多数位于 IGSꎮ 这些 SSR 标记将
出的两个树文件利用 iTOL( https: / / itol.embl.de / ) 为未来研究云杉属物种遗传标记和遗传多样性提
(Letunic & Borkꎬ 2007)在线网址进行美化ꎮ 供可靠的信息ꎮ
丽江云杉叶绿体基因组中包含 29 个串联重复
2 结果与分析 序列(TRS)(图 2)ꎮ 串联重复单元长度为 12 ~ 99
bpꎬ重复次数为 2 ~ 9 次ꎬ总长度为 25 ~ 270 bpꎮ 最
2.1 丽江云杉叶绿体基因组基本特征 短的一段位于 trnV ̄UAC 内含子区域ꎬ是一段长度
GetOrganelle 软件拼接得到的 15 条丽江云杉 为 12 bp 的序列且重复了 2.1 次ꎻ最长的一段为
叶绿体全基因组序列长度一致ꎬ均为 124 207 bpꎬ 270 bpꎬ位于 ycf1 基因内部ꎬ是一段长度为 99 bp
根据 reads mapping 结果ꎬ显示 15 条序列均有较多 的序列且重复了 2.7 次ꎮ 就分布位置而言ꎬ11 个
的 reads 覆盖ꎬ说明丽江云杉叶绿体基因组结构存 位于 IGSꎬtrnR ̄CCG-rbcL 间区包含 4 段 TRSꎬ其余
在同分异构体形式ꎬ结构比较多样ꎮ 选择与参考 位置均只有 1 个ꎻ2 个位于内含子区( intron regionꎬ
序列比对覆盖度较高的一条提交至 GenBank(登录 Intron)ꎬ分别是 petB 和 trnV ̄UAC 基因的内含子区
号为 PV543637)并进行后续分析ꎮ 丽江云杉叶绿 域ꎻ16 个位于 CDSꎬ其中 ycf1 和 ycf3 基因内部分别
体基因组 GC 含量为 38.74%ꎬ呈圆环状ꎬ缺失了一 有 10 个和 3 个ꎬrpoC1 和 rps8 分别包含 2 个和 1 个
段 IR 区序列ꎬ不是典型的四分体结构(图 1)ꎬ但是 TRSꎮ 总体而言ꎬycf1 基因是富含串联重复序列的
丽江云杉叶绿体基因组中存在 4 段短的反向重复 区域ꎮ
序列ꎬ其长度分别为 651、540、353、280 bpꎬ分别包 2.3 密码子偏好性分析
含了 psbI、ycf12、trnI ̄CAU 及 trnH ̄GUG 基因ꎮ 丽江 密码子使用偏好性结果(图 3)显示ꎬ在丽江云
云杉叶绿体基因组共有 113 个基因ꎬ其中包括 74 杉中ꎬ存在 30 个高频密码子( RSCU > 1)ꎬ其中以
个蛋白质编码基因、35 个 tRNA 基因及 4 个 rRNA U、A、G 结尾的分别有 16 个、13 个、1 个ꎬ说明丽江
基因(表 1)ꎬ其中 psbI、ycf12、trnH ̄GUG、trnI ̄CAU、 云杉叶绿体基因组的密码子偏好以 U 或 A 碱基结

