Page 169 - 《广西植物》2023年第10期
P. 169
10 期 赵渊祥等: 猴樟叶绿体基因组特征分析 1 9 2 3
的园林绿化树种和经济树种ꎮ 已有研究主要集中 ncbi.nlm.nih.gov / Blast.cgi) 比对组装的序列ꎬ得到
在其生态特征( 田小琴和韦小丽ꎬ2011)、抗逆性 第二种注释结果ꎮ 然后ꎬ通过对比去除错误注释
(宋芳琳等ꎬ2012)、栽培和育种( 张怡等ꎬ2014)ꎬ 及冗余的注释ꎬ确定多外显子边界ꎬ获得最终的注
组织培养技术(肖祖飞等ꎬ2020) 等方面ꎬ鲜有关于 释ꎮ 最后ꎬ将注释完成的猴樟叶绿体基因组序列
猴樟叶绿体基因组及系统发育方面的研究ꎮ 本研 提交至 NCBI( https: / / www.ncbi.nlm.nih. gov / )ꎬ获
究基于 Illumina 测序平台ꎬ对贵州产地猴樟的叶绿 得登录号(MW381013)ꎮ
体全基因组进行测序、组装和注释ꎬ从叶绿体基因 1.4 叶绿体基因组 SSR 位点分析
组结构、SSR 位点、密码子使用偏性、反向重复区 使用 MISA( v1.0) 软件( Beier et al.ꎬ 2017) 搜
(inverted repeatsꎬIR) 的扩张与收缩以及系统发育 索叶绿体基因组中的 SSR 标记ꎬ设置单核苷酸重
等方面开展研究ꎬ同时结合樟亚科主要属的 23 个 复次数>10ꎬ二核苷酸重复次数>5ꎬ三核苷酸重复
物种构建系统进化树ꎬ拟揭示猴樟叶绿体基因组 次数 > 4ꎬ 四 核 苷 酸、 五 核 苷 酸 和 六 核 苷 酸 重 复
以下信息:( 1) 猴 樟 叶 绿 体 基 因 组 的 基 本 特 征ꎻ 次数>3ꎮ
(2)猴樟叶绿体基因组 SSR 位点及基因密码子偏 1.5 密码子偏好性分析
好性情况ꎻ(3) 猴樟与近缘种的叶绿体基因组 IR 根据 127 个基因的 CDS 序列ꎬ筛选唯一的( 多
区结构差异ꎻ(4)猴樟叶绿体基因组系统发育所属 个拷贝基因选择一个拷贝) 且序列长度大于 300
分支ꎮ bp 的 CDSꎬ并用 CodonW( v.1.4.2) 软件( Wong et
al.ꎬ 2010)估计每个密码子的相对同义密码子使
1 材料与方法 用(relative synonymous codon usageꎬ RSCU)、有效
密码子数(effective number of codonsꎬ ENc) 和密码
1.1 试验材料 子适应指数(codon adaptation indexꎬ CAI)ꎮ
猴樟叶片采自贵州省贵阳市贵州大学南校区 1.6 边界分析
(106°67′ Eꎬ26°43′ N)ꎬ经贵州省森林资源与环境 利 用 近 缘 物 种 新 樟 ( Neocinnamomum
研究中心吴峰教授鉴定为樟科樟属猴樟ꎮ 样本经 delavayi)、 肉 桂 ( Cinnamomum cassia )、 檫 木
液氮速冻处理后保存于-80 ℃ 冰箱ꎮ (Sassafras tzumu)、月桂(Laurus nobilis)和大叶新木
1.2 叶绿体全基因组 DNA 提取与测序 姜子(Neolitsea levinei) 叶绿体基因组ꎬ通过 Irscope
选取叶片组织ꎬ利用植物基因组 DNA 提取试 ( https: / / irscope. shinyapps. io / irapp / ) 在 线 软 件
剂盒( TIANGEN Beijing China) 提取猴樟 DNAꎬ经 (Amiryousefi et al.ꎬ 2018)进行可视化对比ꎬ分析 4
过 1%琼脂糖凝胶电泳和紫外分光光度计检测质 个区域的变化ꎬ特别是 IR 区的扩张和收缩以及
量后利用 Illumina NovaSeq 平台完成高通量测序ꎮ IR / SC 边界基因种类和位置的变化ꎮ
具体实验流程按照 Illumina 公司提供的标准方法 1.7 叶绿体基因组系统发育分析
执行ꎬ将 DNA 片段化后ꎬ对其进行片段纯化、末端 利用猴樟叶绿体基因组和从 NCBI 下载的 23
修复、3′端加 A、连接测序接头ꎬ经 PCR 扩增形成 个樟亚科(Subfam. Lauroideae) 植物叶绿体基因组
测序文库ꎬ文库质检合格后进行测序ꎮ 构建 系 统 发 育 树ꎮ 利 用 MAFFT ( v. 7. 475) 软 件
1.3 叶绿体全基因组的组装与注释 (Nakamura et al.ꎬ 2018)进行多序列比对ꎬ利用 IQ ̄
采用 SPAdes(v3.10.1) 软件( Bankevich et al.ꎬ TREE(v. 2.0.3)软件(Nguyen et al.ꎬ 2015) 构建系
2012)组装叶绿体基因组ꎮ 采用两种方法对叶绿 统发育树ꎬ步长自检 1 000 次ꎮ
体基因组进行注释ꎬ以提高注释的准确性ꎮ 首先ꎬ
2 结果与分析
使用 PRODIGAL(v2.6.3)软件(Hyatt et al.ꎬ 2010)
注释叶绿体的 CDSꎬ使用 HMMER( v3. 1b2) 软件
(Collyda et al.ꎬ 2006)预测 rRNAꎬ使用 ARAGORN 2.1 叶绿体基因组结构
(v.1. 2. 38) 软 件 ( Laslett & Canbackꎬ 2004) 预 测 猴樟叶绿体基因组全长为 152 727 bp( 图 1)ꎬ
tRNAꎮ 其次ꎬ根据 NCBI 上已经公布的近缘物种ꎬ 由图 1 可知ꎬ呈典型的四分体结构ꎬ由 93 605 bp
提取 其 基 因 序 列ꎬ 再 使 用 BLAST ( https: / / blast. 的大单拷贝区(large single copyꎬ LSC)、 18 858 bp