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
   164   165   166   167   168   169   170   171   172   173   174