Page 7 - 《广西植物》2023年第3期
P. 7

3 期             谭小明等: 广西喀斯特地区毛唇芋兰根内与根际土壤真菌群落组成分析                                            4 0 7

            将样品带回广西中医药大学仙葫校区的实验室ꎬ                              ITS1 区域进行高通量测序ꎮ 测序结果表明ꎬ质控
            随即在流水下冲洗干净根表杂质ꎬ使用滤纸吸干                              后每个样品的测序标签数 均 超 过110 000 条 ( 表
            根表水分后ꎬ迅速采用液氮急冻ꎬ之后置于-80 ℃                           1)ꎻ以测序数据量对应 OTU 数量构建的稀释曲线
            冰箱中保存ꎬ备用ꎮ                                          趋向平坦( 图 1)ꎬ说明测序数据量能够覆盖样品
                 根际土壤采自毛唇芋兰根系表面 0 ~ 1 cm 范                     中绝大多数真菌种类ꎬ测序数据量合理ꎮ 每个样
            围内的土层ꎬ将采集的土壤置于无菌离心管中ꎬ编                             品的覆盖率均达到 0.999( 表 1)ꎬ说明样本中有序

            号分别为 DXTR(大新土壤) 和 LZTR( 龙州土壤)ꎻ                     列没有被测出的概率极低ꎬ测序深度合格ꎮ
            分别取大新毛唇芋兰和龙州毛唇芋兰的主根、走                                  Chao 1 指数和 Ace 指数用于评定样品中群落

            茎(各 3 个重复)置于无菌离心管中ꎬ编号为 DXZG                        的丰富度ꎻ香农指数则可反映出群落分布多样性ꎬ
            ( -1ꎬ-2ꎬ-3)、DXZJ( -1ꎬ-2ꎬ-3)、LZZG( -1ꎬ-2ꎬ           其受到样品群落中物种丰富度以及物种均匀度的
            -3)、LZZJ( -1ꎬ-2ꎬ-3)ꎮ                               影响(庞美霞等ꎬ2021)ꎮ 大新主根( DXZG) Chao 1
            1.2 样品基因组 DNA 的提取和分析                               指数、Ace 指数和香农指数分别为 302. 67 ± 7. 23、
                 样品经冷链送至广州基迪奥生物科技有限公                           316. 04 ± 3. 41、 1. 10 ± 0. 23 ( 图 2 : A)ꎬ 龙 州 主 根
            司ꎬ进行 DNA 的提取和 ITS rDNA 的测序分析ꎻ采                     (LZZG)Chao 1 指数、Ace 指数和香农指数分别为
            用带有 barcode 的特异引物扩增 ITS 的 ITS1 plant               212.72±143.99、215.14±145.47、2.23±0.15ꎮ 龙州主
            区 ( 引 物 序 列 分 别 为 ITS1 F KYO2: 5′ ̄                 根群体中反应群落多样性变化的香农指数( P = 0.
            TAGAGGAAGTAAAAGTCGTAA ̄3′ 和 ITS86R: 5′ ̄             003 5) 显著性上升( 图 2:A)ꎬ而表征物种丰度的
            TTCAAAGATTCGATGATTCAC ̄3′)ꎻ 将 纯 化 后 的               Chao 1 指数(P>0.05)和 Ace 指数(P>0.05)相较于
            PCR 产物(即扩增子) 连接测序接头ꎬ构建测序文                          大新主根群体均有所下降但无明显差异ꎮ 大新走
            库ꎬ上机进行 Illumina MiSeq 高通量测序ꎬ并进行                    茎(DXZJ)Chao 1 指数、Ace 指数和香农指数分别为
            相应的信息分析ꎮ                                           205.41 ± 113.57、210.70 ± 118.47、0.72 ± 0.25( 图 2:
            1.3 数据处理和分析                                        B)ꎬ龙州走茎(LZZJ) Chao 1 指数、Ace 指数和香农
                 过滤低质量 Readsꎬ将双端 Reads 拼接为 Tagꎬ                指数分别为 209.37±47.62、212.68±73.34、0.69±0.11
            对 Tag 进 行 低 质 量 过 滤ꎬ得 到 的 数 据 称 为 Clean            (图 2:B)ꎮ 由此可见ꎬ龙州走茎的真菌物种丰度及
            Tagꎻ使 用 Usearch 软 件 进 行 聚 类 ( Schloss et al.ꎬ      多样性低于大新走茎组ꎬ但差异不显著(P>0.05)ꎮ
            2009)ꎬ去除聚类过程中检测到的嵌合体 Tagꎬ按照                        2.2 根系共生真菌多样性及相关性分析
            97%相 似 性 进 行 可 操 作 分 类 单 元 ( operational               从毛唇芋兰 2 个根际土壤( DXTR、LZTR) 样
            taxonomic unitꎬOTU)聚类分析ꎬ获得 OTU 后ꎬ基于                品、6 个 走 茎 ( DXZJ、 LZZJ) 样 品 和 6 个 主 根
            Effective Tag 进 行 OTU 丰 度 统 计ꎮ 利 用 R 语 言           (DXZG、LZZG)样品内检测到大量 OTUꎮ 首先ꎬ不
            vegan 包进行维恩图分析( Oktalira et al.ꎬ 2019) 和           同土壤样品分析结果表明ꎬDXTR 共获得 411 个
            使用 Qiime 软件( Schloss et al.ꎬ 2009) 进行 alpha        OTUꎬ多于 LZTR 中的 OTU 数量 406ꎻDXTR 的真菌
            多 样 性 指 数ꎬ 包 括 物 种 数 ( Sobs)、 香 农 指 数              香农指数、辛普森指数、Chao 1 指数以及 Pd 指数
            (Shannon 指数)、辛普森指数( Simpson 指数)、丰                  均略高于 LZTR(表 1)ꎬ说明大新毛唇芋兰根际土
            富度指数( Chao 1 指数)、基于丰度的覆盖估计值                        壤真菌群落在丰富度和多样性上均高于龙州毛唇
            (Ace 指数)、谱系多样性指数( Pd 指数) 及覆盖率                      芋兰ꎮ 另外ꎬOTUs ̄Venn 图能够直观反映样本间相
            (Good’s coverage)分析ꎻ同时ꎬ采用统计学方法进                    同及其独特的 OTUꎮ DXTR 与 LZTR 真 菌 OTUs ̄

            行差异特征发现和显著性检验ꎮ                                     Venn 图如图 3 所示ꎬ2 个土壤样本中共获得 705
                                                               个 OTUꎬ它们之间共有的 OTU 为 112 个ꎬDXTR 特
            2  结果与分析                                           有 OTU 为 299 个ꎬLZTR 特有 OTU 为 294 个ꎮ 因

                                                               此ꎬ综合 Alpha 多样性指数及 OTUs ̄Venn 图发现ꎬ
            2.1 测序结果分析                                         大新和龙州的毛唇芋兰根际土壤真菌的组成均较
                 通过对大新县和龙州县的野生毛唇芋兰主根                           为丰富ꎬ但它们在真菌的组成、结构和相对丰度上
            和走茎及根际土壤共 14 个样品的真菌 rDNA 基因                        存在一定的差异ꎮ
   2   3   4   5   6   7   8   9   10   11   12