叶绿体基因组分析技术在中药资源类专业实验教学中的应用.pdf

收稿日期：通信作者：张朝凤，：基金项目：江苏省高等教育教学改革研究课题，编号 ；中国药科大学高层次人才引进经费项目 陈凌云叶绿体基因组分析技术在中药资源类专业实验教学中的应用陈凌云，答李穆，郑文芳，张朝凤中国药科大学中药学院（南京 ）摘要近年来，叶绿体基因组数据被广泛用于解决植物系统发育问题和中药材鉴定，但我国高等院校中药资源类专业现有课程体系却普遍缺少叶绿体基因组学相关理论教学与实践教学环节。中国药科大学中药资源类专业实验教学使用第二代 测序数据，开展叶绿体基因组的组装、注释及系统进化分析的相关教学活动。实践证明，引入叶绿体基因组实验教学，有助于学生深入理解生物信息学在中药资源开发中的应用，培养学生生物统计分析的能力。关键词叶绿体基因组；中药资源类实验；生物信息学 牞 牞 牞 牗 牞 牞 牞 牘 牶 牞 牞 牞 牞 牞 牶 牷 牷 近年来，中药资源相关研究取得了许多重要成果，尤其在药用植物基因组学分析方面。但国内高校中药资源类本专科生教学活动中，植物基因组学实践教学较少或处于空白状态；传统的中药资源相关学习活动以实物辨识和记忆背诵为主，围绕植物传统分类学与本草沿革等展开。为了提升中药资源与开发专业建设水平，促进学生掌握药用植物基因组学相关知识与技能，教学团队在中国药科大学中药资源类专业的实验教学中引入了叶绿体基因组的组装、注释及进化分析等实验内容。教学实践证明，该实验教学的实施有助于学生深入理解生物信息学在中药资源开发中的应用，提升学生生物统计技术的实践能力。?植物基因组序列分析技术的现状 年美国能源部提出了“人类基因组计划”（）草案 ，年 正式启动。年，科学家发表了拟南芥的全基因组序列信息 。截至 年 月，科学家已完成 多个被子植物物种的全基因组测序（：）。植物基因组学已被广泛应用于药用植物鉴定、繁殖育种、分类演化等研究。叶绿体基因组一药学教育 年 第 卷 第 期般是一个双链环状 ，独立于核基因组，在细胞中拷贝数较多，大小在 到 之间，存在于藻类和绝大多数的绿色植物细胞中。叶绿体基因组较为保守，含有大量功能基因，对于物种鉴定、分类、演化研究有着极为重要的作用。目前，科研工作者普遍采用第二代测序平台获得植物基因组 数据，再通过算法提取叶绿体的小片段 序列（），拼接成片段（），再连接成完整的或者间断的长序列（）。随着测序技术的发展，测序准确率提高，测序价格大幅下降，已完成叶绿体全基因组测序的物种数目正在爆发式增长。近年来，科学家相继报道了人参、三七、丹参等重要药用植物的全基因组序列，本草基因组学应运而生 ，对中药学发展产生了重要影响。例如，人参基因组数据的发表促进了人参育种和人参皂苷生物合成研究 ，黄花蒿的全基因组测序不仅发掘了与青蒿素合成相关的基因，而且使得通过异源基因表达增加青蒿素产量成为可能 ，冬虫夏草菌的全基因组测序使得虫草素的生物合成有了可靠的理论基础 。植物基因组技术在中药学与药用植物学研究中发挥着重要作用，?叶绿体基因组分析技术的教学需求叶绿体基因组学在药用植物学与中药鉴定学相关的教学、研究中有着重要作用。目前，国内 药用植物学 教材 ，仍然采用恩格勒分类系统。而基于基因组学的分类系统如 ，相比恩格勒分类系统更为接近植物真实的进化历史，已经被国外教材和科研工作广泛采纳。在药用植物基源鉴定方面，基于基因组学的分子鉴定相比于形态和显微鉴别更加高效、准确，但中药资源类专业本专科教学实践却较少提及相关技术。为了提升中药资源类专业教育教学水平，医药院校需在中药鉴定学、药用植物学等理论课程和实验课程中增加叶绿体基因组分析相关内容。实验的开设一方面有益于培养学生科研思维，提高学生的实验操作技能；另一方面能够充分调动学生的学习积极性，拓宽学生科研视野。当前，部分高校（如北京协和医学院、武汉理工大学、上海中医药大学等）面向本专科生或研究生已开展了基因组学相关课程教学，但大部分中医药类高校尚无植物基因组学相关教学内容。网络技术的普及，大大降低了基因组学实验教学成本，面向本专科学生开展基因组学相关实践教学成为可能。开展基因组学实验教学，需要的计算机条件如下：以上或 以上，以上的可用运行内存空间，以上的可用硬盘空间。处理器的性能和线程数量主要影响运行速度，近十年来市售的服务器几乎都能满足本实验的需求。叶绿体基因组的组装在 系统上完成，注释和系统进化分析可以在 系统，也可以在 系统上完成。笔者采用的服务器品牌和配置如下：、两 颗、内存、操作系统。主要软件为：、。?叶绿体基因组分析技术的教学案例以西南山梗菜的第二代测序数据为实验材料，笔者设计了叶绿体基因组的分析实验。在知识体系和内容上该实验考虑了学生综合运用知识、解决和探索科学问题的能力，让学生了解和掌握植物基因组学知识。如果学生没有接触过 系统，可能较难完成本实验内容。因此，在开展本实验前，笔者建议用 个课时向学生讲授 相关知识，或让学生自学相关知识。对多个叶绿体基因组进行序列比对并构建系统进化树是教学过程中的主要难点，需要教师在讲解这些内容时更为细致。完成本实验全部流程需要大约 小时，约 个课时。西南山梗菜（）属于山梗菜科，半边莲属，又名野烟，小烟草。记载于 滇南本草：“味辛、麻，性温，有大毒（：）。现代药理研究认为，西南山梗菜有消炎止痛、解毒、杀虫止痒的功效。常用于治疗热毒疮痈、跌打损伤、风湿性关节炎等病症。本实验西南山梗菜采集于重庆市梁平区曲水镇聚宝村（年 月，图 ），种植于中国药科大学药用植物园。笔者取西南山梗菜的叶片送天津诺禾致源科技有限公司使用 平台测序，测序数据用于后续叶绿体基因组组装与分析。实验教学基本流程如图 。?实践过程 叶绿体基因组组装实验使用 软件对原始测序数据去除接 头 和 低 质 量 。从 ：？获取并安装。也可以使用 （：药学教育 年 第 卷 第 期 ）安 装 。在 官 网 下 载 的 系统的安装包，上传到服务器后，使用 指令执行即可。安装完成后，使用 安装 软件。图 实验所用的实验材料 西南山梗菜图 实验教学的基本流程 的执行命令如下：指令需要输入经双末端测序得到的 个数据，并指定输出的 个文件的名称。指定名称后，该指令可 以 支 持 多 个 参 数，以 设 定 执 行 的 步 骤。“”为 软件包自带的一个包含接头序列的文件。“”代表切除与接头序列和反向序列匹配的质量值小于 的 ，“”代表切除与接头序列相匹配的 的最小质量值。“”，从 起始开始切除碱基，“”代表切除质量值小于 的碱基。“”，从 末端开始切除碱基，参数含义同上。“”，划窗剪切，会以窗口扫描过整个 ，再按照窗口切除碱基。“”代表窗口的大小为 个碱基，“”代表如果这一条 的窗口的平均质量值低于 ，则删除一整条 。“，设定 的最小长度。“”代表切除总长度低于 的 。此步骤大约需要 分钟，产生的 个 文件将用于下一步分析。采用 软件进行叶绿体基因组的组装，软 件 从 ：下载，并按照说明书安装。也可以在 中使用指令 完成 的安装。采用如下命令执行叶绿体基因组组装：，各项参数，“”和“”表示正向和反向的测序数据，“”表示需要组装的基因组类型，“”代表输出的结果文件夹的名称，“”代表使用的线程数，“”代表组装时的 长度。在多数组装中，“”设定为 和 能取得较好的结果。默认软件参数下，将自动估计组装需要的 数量。可以添加 参数来设置使用 的数据，但需要注意，会将全部的测序数据一次性写入内存，所以此项需要谨慎设置。结果为一个文件夹，其内容如图 。图 输出的结果从上至下，第一个文件夹为 的输出数据，的组装核心采用的是 ，文件药学教育 年 第 卷 第 期夹内含有上述所设定的不同 的组装结果。第二个文件夹为参考基因组库。第三到第四个文件为 产生的完整的、环装的基因组序列。由于同一个植物的叶绿体基因组存在一定差异的拷贝，因此会有含有编号 、的两个序列，挑选其中一个用于后续分析。第五个文件即为最终结果。第六到第九文件为使用 比对到参考基因组后，能成功匹配的 。第十到第十二个文件是没有去除非成环或非主 链序列的组装结果。（如果最终结果未成环，需要用第十到第十二个文件进行进一步分析）。最后一个文件是运行日志，记录了组装的主要过程。将图 后缀为 的文件导入到软件 。（：）是一个可视化的基因组修改软件，能够直观地分析前述组装的结果。使用 导出一段大单拷贝区，一段小单拷贝区和重复拷贝区的两个部分，也可以得到最终组装结果。在测序数据不完整，测序错误、测序污染未去除或去除不完全等情况下，有可能得到类似如图 的结果，是 ，区别于图 的 ，只能得到叶绿体基因组的大片段。图 由 可视化完整的组装结果图 由 可视化不完整的组装结果不完全组装的叶绿体基因组也能用于许多研究，如系统进化分析。笔者尝试了不同的组装参数（如不同的 值，的数量），一般可以得到完整的叶绿体基因组序列。输入同样的 ，执行相同的命令，其结果相同，具有可重复性。叶绿体基因组注释使用 软件进行叶绿体基因组在线注释，网址为 ：。根据提示，选择“”，导入上一步获得的叶绿体基因组序列，并在参考基因组中选择“”作为参考基因组。提交任务，并记录订单号。等待约 分钟，即可查看结果。结果主要包括四部分：文件，详细描述注释得到的基因；文件，描述注释所得到的基因在基因组上位置；一个绘制好的叶绿体基因组圈图；原始数据整理得到的 文件。在叶绿体基因组圈图上可以看到基因的功能和所在位置，也可以将 文件导入其他圈图绘制和美化软件进行自定义操作。除以上四部分内容，还提供了预测基因的蛋白质序列以及较为粗略的 位点分析结果。功能与进化分析对组装结果进行分析，可以得到更多的特征信息，也可以联系生物化学等课程，学习课程中一些概念的实际用途。笔者选择了 重复序列分析、基因组偏好性分析、系统进化分析三个项目。重复序列分析 重复序列分析采用 软件进行。可以鉴定简单重复序列，其提供了在线版和安装版。笔者选择在线版，网址为 ：？。在网页中上传叶绿体基因组全长序列，填入邮箱以及项目名称，一分钟内即可在邮箱中收到分析结果。结果是一个 格式的压缩文件夹，文件夹中包含上传的序列本身，所选择的分析配置，每一个 重复序列的详表文件和一个重复序列的总述文件。在这里笔者设定的项目名称为 ，个文件与上述内容按顺序相互对应，文件可用记事本打开。后缀为 的文件包含了每个 位点的详情，后缀为 的文件包含了 位点的总数和类型等信息。基因组偏好性分析基因组偏好性分析采用 软件进行，软件从 ：获得，将组装得到的叶绿体基因组序列文件和 放置在同一个文件夹。药学教育 年 第 卷 第 期打开 ，可以看到一个命令行窗口（图）。图 的一级菜单输入 ，载入需要分析的序列；输入 ，进入子菜单，用来选择需要计算的内容。设定好后，输入 进行计算，计算完毕后输入 退出，不能直接关闭窗口，否则结果将不会写入文件，无法获得。计算完成后将得到一个 文件和一个 文件，两者均是表格，都能用记事本打开。文件是上文设定的计算内容得到的结果，文件分别列出了每种密码子的占比数量等信息。系统发育分析选取桔梗科中党参、沙参、桔梗等 种药用植物，从 下载这些物种的叶绿体基因组序列，以邻接法构建系统进化树。将下载的序列合并到同一个 文件中。采用 软件进行序列比对，软件可以输入 命令获取。在完成 的安装后，使用指令 即可完成比对。这个过程大约需要 小时。将结果导入 （：）中，选择 ，在弹出的窗口中将其保存为 文件。然后打开 选项中的 ，构建进化树。结果可以在 分钟内得到，可以将结果保存为 文 件，在 网 页 中 进 行 修 饰（：）。这样就得到了西南山梗菜完整的叶绿体基因组，总长度为 。该叶绿体基因组具有典型的环状四分体结构，其中，区长度为 ，区 域 长 度 为 ，两 个 重 复 拷 贝 区 为 。注释基因组得到了 个基因。其中，光合作用相关的基因有 个，自我复制相关的基因有 个，蛋白编码相关的基因 个和未知功能的开放性阅读框 个（图 ）。微卫星序列分析得到了 个 位点。其中，有 个位于 区域，个位于 区域，个位于图 西南山梗菜叶绿体基因组的注释结果 区域。这些 位点中，有 个为单核苷酸 和构成的重复序列、个为 和 组成的重复序列、个由 构成的三核苷酸重复序列、个复合型序列。基因组偏好性分析得到密码子适应指数为 ，表明西南山梗菜叶绿体基因组密码子偏好性较弱，内源基因表达较低。整个基因组的 含量为 ，表明西南山梗菜叶绿体基因组偏好使用碱基 和 。最后，以邻接法构建系统进化树，值越高，表示进化枝的可信度越高。由系统进化树可知，西南山梗菜和 同聚一支，表明 二 者 系 统 关 系 较 近，和 沙 参 、党参 等不在同一支，表明西南山梗菜和后两者系统关系较远（图）。本实验产生的高通量测序数据存储于 （编号：）。涉及的数据处理方法和关键数 据 详 见 ：。?教学效果反馈该实验安排 节理论课（分钟一节课）；节上机课（分钟一节课），进行实际操作练习。理论课结束后，教师向学生发放详细的上机教程。上机课堂上，部分学生在教师讲解后，需要多次协助才能够进行后续的上机练习。造成这种情况的原因，可能是这部分学生没有提前预习。所有学生均上交了教师布置的课后作业。由于大家使用同样的分析数据，大家药学教育 年 第 卷 第 期的作业大体相同，存在部分学生复制他人作业的可能性。为避免抄袭，笔者要求学生将叶绿体序列的名称额外加上自己的姓名和学号。此外，可以向学生发一份全基因组测序的原始数据，学生在课后独立完成全基因组组装，这将帮助有余力的学生巩固知识并提高水平。图 西南山梗菜和同科植物的系统进化树注：系统进化树使用 软件的 方法构建而成，黄花蒿是外类群学生的课后作业是独立完成叶绿体基因组的组装与注释、系统进化分析，提交叶绿体基因组的注释圈图、不同植物的叶绿体基因组的序列比对截图以及最终构建的系统进化树。学生都完成了此作业，大约有三分之一的学生额外上交了叶绿体注释的 文件。此外，有 位接受本课程的学生在做毕业论文设计，内容是药用植物基因组学方面的工作。笔者让接受过该课程的 位学生帮助其他没有接受过此课程的学生使用服务器、做简单的分析。?教学案例的实用性分析首先，为了让学生较好地掌握叶绿体基因组的组装方法，实验需选择合适的物种和高质量的测序序列。笔者曾尝试组装 （杜衡）与 （水葱）的叶绿体基因组，然而两个物种的组装结果都不尽如人意。杜衡的组装结果为相互不能连接的七个序列，使用 将杜衡的序列比对 数据库，发现即使选取最相似的基因组，其结果也相差甚远，多次调整，均未得到完整的叶绿体基因组序列。检查组装过程中，笔者发现许多短小而复杂的序列，推测是测序数据中核基因组过多，即使进行深度筛选也不能完全分离，导致无法组装。水葱的叶绿体基因组组装得到了许多成型的质粒和诸多大小不同的序列，推测可能是样本受到了污染，或者核酸遭到破坏。教师可以在理论教学中阐明有些植物物种的叶绿体基因组较难组装的原因。其次，西南山梗菜的叶绿体基因组具有典型的环状四分体结构，本实验所介绍的方法仅适用于含有标准环状四分体的叶绿体基因组的组装，部分豆科植物和蕨类植物可能并不在此组装流程的适用范围内。由于二代测序的固有问题，组装结果可能出现 和杂乱的分支序列，这种结果可以通过修改 的 参数，使用更多的 或者手动分析结果、手动延长等方式进行修正。也可能存在叶绿体基因组确实没有被测序覆盖，从而无法组装成完整基因组的情况。这种情况可以使用同源植物作为参考基因组进行填补（但可能存在序列错误），或重新测序。为保证教学质量，笔者建议以本文提供的西南山梗菜的二代测序数据开展实践教学（，编号 ）。最后，运行 和 软件会消耗较长时间，且需要至少三倍于测序 的存储空间，因此应该保证计算机能够长时间运行，且有足够的存储空间（）。价格在 万元之间的小型服务器即可满足 名以上学生在学习 操作系统的同时开展叶绿体组装分析等教学活动。也可以借助学校网络数据中心的公共服务器开展本教学活动。学生在学习有关核酸、植物分类学的相关课程内容后，开展本实践教学更为适当。?结语叶绿体基因组分析不只是植物学、中药鉴定学、基因组学研究的有力工具，也是药用植物学、中药鉴药学教育 年 第 卷 第 期定学等课程亟需补充的知识点。借助本专科教育体系，将基因组知识融入教学，既能使教学内容与专业发展技术需求同步，也能使学生更有成就感，为中药资源类相关科技工作夯实基础。参考文献 ，（）：，：陈士林 本草基因组学：中药组学的发展与未来 北京：科学出版社，张景景 基于人参全基因组的多组学数据解析人参皂苷生物合成 武汉：湖北中医药大学，李星，郝鹤，池剑亭，等 利用异源生物生产青蒿素及其前体的研究进展 植物学报，（）：吴娟娟，吴明月，贾惠，等 冬虫夏草中虫草素相关基因的克隆及生物信息学分析 交通医学，（）：，黄宝康 药用植物学（版）北京：人民卫生出版社，刘春生 药用植物学（版）北京：中国中医药出版社，李德铢 中国维管植物科属词典 北京：科学出版社，李德铢 中国维管植物科属志 北京：科学出版社，：，（）：，：，（）：，：，（）：，：，（）：欢 迎 订 阅 年 药学教育 是由教育部主管，中国药科大学主办的高等药学教育理论及实践研究期刊，受到广大医药教育工作者的欢迎。药学教育 是中国学术期刊综合评价数据库统计源期刊，被中国核心期刊（遴选）数据库、国家哲学社会科学期刊数据库（）、万方数据库、中国期刊全文数据库等多家国内大型检索数据库全文收录。本刊旨在探讨各层次药学教育规律，研究药学教育理论，发表药学教育改革成果，结合实际介绍药学教育先进经验。刊物力求突出思想性、学术性、实用性，在药学教育改革中起到宣传、引导、咨询、借鉴、交流作用，是广大药学教育工作者的良朋益友。药学教育 设有教育研究、学科与课程建设、课程思政、教学学术、调研与评估、药学人物等多个栏目，内容丰富，形式新颖，有较大的信息容量。药学教育 为双月刊，逢双月 日出版，大 开本，页，每册定价 元，全年定价 元。欢迎到当地邮局或直接扫描本刊封底二维码订阅，邮发代号：。药学教育 年 第 卷 第 期