Page 94 - 《广西植物》2024年第10期
P. 94

1 8 9 6                                广  西  植  物                                         44 卷
            叶片与健叶之间的微生物群和化学因子差异ꎬ但                              1.3 数据质控及群落结构与多样性分析
            未对不同病害程度的野火病烟叶做出进一步分                                   下 机 数 据 经 Cutadapt ( https:/ / cutadapt.
            析ꎮ 烟草野火病不但在田间为害烟叶ꎬ而且在烟                             readthedocs.io / en / v1.9.1 / )软件剔除、过滤低质量序
            叶采收后至烘干前仍可继续为害ꎬ导致叶片病斑                              列ꎬ以及 UCHIME 算法去除嵌合体序列ꎬ最终得到
            继续扩大变成枯焦状而失去烘烤价值ꎮ 因此ꎬ亟                             有效 序 列 ( clean reads)ꎮ 使 用 UPARSE ( http:/ /
            须阐明其病害流行程度与叶际微生物群落的动态                              www.drive5.com/ usearch / )软件在相似度≥97%的条
            变化ꎮ 鉴于此ꎬ本研究以不同发病程度的烟草野                             件下对 clean reads 进行 OTU( operational taxonomic
            火病烟叶为研究对象ꎬ采用高通量测序技术ꎬ拟探                             units) 聚类ꎬ筛选每个 OTU 的代表性序列ꎬ并基于
            讨以下问题:(1)不同程度发病烟叶病斑组织与健                            Mothur 算法将 OTU 代表性序列与 Unite(7.2) 数据
            康组织中微生物群落结构与多样性差异ꎻ(2) 不同                           库和 Silva132 数据库比对、注释物种分类(Edgar &
            发病程度病斑组织微生物群落的变化ꎻ(3) 不同发                           Robertꎬ 2013ꎻ Urmas et al.ꎬ 2013)ꎮ 利用 R 语言工
            病程度烟叶组织微生物的功能变化ꎮ 旨在为认识                             具绘 制 物 种 群 落 柱 状 图、 主 成 分 分 析 ( principal
            烟草野火病的演变规律提供科学依据ꎮ                                  component analysisꎬ PCA) 图和 Venn 分析图等ꎮ 利
                                                               用 QIIME(https:/ / docs.qiime2.org)按样本最小序列
            1  材料与方法                                           数进行样本序列抽平后ꎬ计算样品微生物群落 α 多
                                                               样性指数( Shannon、Simpson、Chao1 等)ꎻ使用 DPS
            1.1 样品采集                                           软件(http:/ / www.dpsw.cn / download.html)进行不同
                 2019 年 6 月ꎬ在种植品种为云烟 87 的贵州省                   组别物种相对丰度的差异性分析ꎮ 采用 Tax4Fun 程
            施秉 县 烟 叶 产 区 ( 107° 51′ 36″—108° 28′ 51″ E、        序语言及 FUNGuild 和 KEGG 数据库分析真菌、细菌
            26°49′47″—27°21′15″ N)采集移栽 60 d 后发生烟               群落的功能类群(孙美丽等ꎬ2023)ꎮ
            草野火病的烟叶样品ꎬ根据病斑面积占比分为轻
                                                               2  结果与分析
            度(5% ~ 10%)、中度(11% ~ 20%) 和重度(21% ~
            50%)叶片(图 1)ꎬ依次记为 A、B、Cꎮ 用无菌剪刀
            剪取样品叶片的病斑组织( b) 和健康组织( j)ꎬ分                        2.1 测序深度评估及数据质控
            别混合均匀置于 50 mL 无菌离心管中ꎬ每个样品 3                            稀释性曲线结果显示ꎬ叶际细菌的 OTU 数在
            管ꎬ每管 5 gꎮ 样品采集后放入低温保存箱ꎬ并迅                          测序数据量为 2 300 时曲线趋于平坦(图 2:A)ꎬ叶
            速带回实验室ꎬ置于-80 ℃ 冰箱保存、备用ꎮ                            际真菌的 OTU 数在测序数据量为 40 000 时曲线
            1.2 基因组 DNA 提取、PCR 扩增及高通量测序                        趋于平坦( 图 2:B)ꎬ表明测序数据已包含所测样
                                  ®
                 参 照 Fast DNA        SPIN 试 剂 盒 ( MP           品中绝大多数的物种信息ꎬ测序结果合理ꎮ
            Biomedicals)使用说明ꎬ提取叶际微生物基因组总                           原始测序数据质控后ꎬ平均每组样品的细菌和
            DNAꎬ 经 超 微 量 分 光 光 度 计 NanoDrop 2000               真菌分别获得 80 092 条和 70 147 条高质量序列(表
            (Thermo Fisher Scientific)检测质量与浓度后ꎬ用无              1)ꎮ 细菌与真菌原始测序数据的 GenBank Accession

                                         ̄1
            菌水将 DNA 稀释至 1 ngμL 备用ꎮ                           Number 分别为 PRJNA1025832、PRJNA1025370ꎮ
                 细菌 16S rRNA 基因和真菌转录间隔区 ITS1                   2.2 叶际细菌与真菌群落结构
            区分别采用引物 515F / 806R(Yuan et al.ꎬ 2020)和            2.2.1 门水平叶际群落结构  健康与病斑组织的
            ITS5 ̄1737F / ITS2 ̄2043R(Wu et al.ꎬ 2019) 进行扩       细菌群落包括变形菌门( Proteobacteria)、厚壁菌门

            增ꎮ PCR 扩增体系及程序参照 Apprill 等(2005)ꎮ                  (Firmicutes)和拟杆菌门( Bacteroidetes)ꎬ其中病斑
            扩增 产 物 经 2% 浓 度 琼 脂 糖 凝 胶 电 泳 检 测、                 组织变形 菌 门 的 相 对 丰 度 高 于 健 康 组 织 ( 图 3:
            GeneJET 胶 回 收 产 物 纯 化 试 剂 盒 ( Thermo Fisher        A)ꎮ 轻度、中度和重度病叶病斑组织的优势细菌
            Scientific)纯化后送至北京诺禾致源科技股份有限                       门为变形菌门( 相对丰度分别为 49.44%、72.39%
            公 司ꎬ 利 用 Ion Plus Fragment Library 试 剂 盒           和 61.92%)ꎬ其次为厚壁菌门( 相对丰度分别为
            (Thermo Fisher Scientific ) 构 建 文 库 并 在 Ion        0.01%、1.54%和 0.73%)和拟杆菌门( 相对丰度分
            S5TMXL 平台上进行测序ꎮ                                    别为 0.06%、0.19%和 0.06%) ꎮ 健康组织优势细
   89   90   91   92   93   94   95   96   97   98   99