Microsoft Word xb 刘开辉.doc

Similar documents
Microsoft Word tb 赵宏宇s-高校教改纵横.doc

<D2BDC1C6BDA1BFB5CDB6C8DAD7CAB8DFB7E5C2DBCCB3B2CEBBE1C3FBB5A52E786C7378>

PCA+LDA 14 1 PEN mL mL mL 16 DJX-AB DJ X AB DJ2 -YS % PEN

Microsoft Word 聂雪梅.doc

[1] Nielsen [2]. Richardson [3] Baldock [4] 0.22 mm 0.32 mm Richardson Zaki. [5-6] mm [7] 1 mm. [8] [9] 5 mm 50 mm [10] [11] [12] -- 40% 50%

~ 10 2 P Y i t = my i t W Y i t 1000 PY i t Y t i W Y i t t i m Y i t t i 15 ~ 49 1 Y Y Y 15 ~ j j t j t = j P i t i = 15 P n i t n Y

填 表 说 明 1. 本 表 所 填 数 据 截 至 时 间 为 2011 年 9 月 底 2. 本 表 请 用 A4 纸 双 面 打 印, 加 盖 学 校 公 章 后 上 报 3. 表 内 所 填 数 据 请 学 校 认 真 核 实, 确 保 准 确 无 误 4. 建 设 完 成 情 况 对 照

2 : 111 ( Pinctada maxima), 13, 11,,,,,, 15, 25 cm, cm, 4 5 kg 1-2],,, 3-7] 24] 25-26],,, ] min, 1, 5 m 4 m, 20 m , 9-22],,, (

Dan Buettner / /

中国轮胎商业网宣传运作收费标准

* CUSUM EWMA PCA TS79 A DOI /j. issn X Incipient Fault Detection in Papermaking Wa

TGF-β AngⅡ B SD ~ 220g SPF. SCXK No SYXK ~ 25 40% ~ 70% OR37G-C

~ a 3 h NCEP ~ 24 3 ~ ~ 8 9 ~ km m ~ 500 m 500 ~ 800 m 800 ~ m a 200

現代學術之建立 陳平 美學十五講 淩繼堯 美學 論集 徐複觀 書店出版社 的方位 陳寶生 宣傳 敦煌文藝出版社 論集續篇 徐複觀 書店出版社 莊子哲學 王博 道家 的天方學 沙宗平 伊斯蘭教 周易 經傳十

Microsoft Word - 1--齐继峰_new_.doc

物理学报 Acta Phys. Sin. Vol. 62, No. 14 (2013) 叠 [4]. PET 设备最重要的部件就是探测器环, 探测 备重建图像具有减少数据插值的优势. 器环的性能直接影响 PET 的成像能力. 探头与探头 之间得到的符合直线叫做投影线. 所有的投影线在

林学 园艺

untitled

untitled

畜牧 动物医学 蚕 蜂

MHz 10 MHz Mbps 1 C 2(a) 4 GHz MHz 56 Msps 70 MHz 70 MHz 23 MHz 14 MHz 23 MHz 2(b)

240 生 异 性 相 吸 的 异 性 效 应 [6] 虽 然, 心 理 学 基 础 研 [7-8] 究 已 经 证 实 存 在 异 性 相 吸 异 性 相 吸 是 否 存 在 于 名 字 认 知 识 别 尚 无 报 道 本 实 验 选 取 不 同 性 别 的 名 字 作 为 刺 激 材 料, 通

在 培 养 职 前 卓 越 化 学 教 师 的 院 校, 会 编 一 本 过 去 称 作 化 学 教 学 论 实 验, 现 在 拟 为 卓 越 化 学 教 师 教 育 实 验 教 学 研 究 的 教 材 各 院 校 对 这 门 课 程 所 给 的 学 时 不 太 一 样, 但 都 是 围 绕 实 验

Vol. 22 No. 4 JOURNAL OF HARBIN UNIVERSITY OF SCIENCE AND TECHNOLOGY Aug GPS,,, : km, 2. 51, , ; ; ; ; DOI: 10.

1 GIS 95 Y = F y + (1 F) (1) 0 0 Y0 kg/hm 2 /day F y 0 y c kg/hm 2 /day [12] y m 20 kg/hm 2 /hour Y = cl cn ch G [ F( y ) T m yo + (2) (1 F)(

CONTENTS 目 录 大 事 要 闻 2 两 学 一 做 5 科 研 动 态 8 交 流 合 作 14 队 伍 建 设 17 人 才 培 养 19 运 行 管 理 21

United Nations ~ ~ % 2010

a b

相 关 报 道 近 年 来 临 床 上 运 用 多 种 针 灸 疗 法 治 疗 肥 胖 及 其 引 起 的 并 发 症, 如 高 脂 血 症 糖 尿 病 高 血 压 及 痛 经 等, 均 获 得 满 意 疗 效. caused by sedentary lifestyle and genetic f

Scoones World Bank DFID Sussex IDS UNDP CARE DFID DFID DFID 1997 IDS

室内设计2015年第5期.indd

H 2 SO ml ml 1. 0 ml C 4. 0 ml - 30 min 490 nm 0 ~ 100 μg /ml Zhao = VρN 100% 1 m V ml ρ g

Mixtions Pin Yin Homepage

Microsoft Word 定版

Microsoft Word doc

km km mm km m /s hpa 500 hpa E N 41 N 37 N 121

Microsoft Word 李强.doc

% % % 1. 39% % 2000 ~ % 1

/MPa / kg m - 3 /MPa /MPa 2. 1E ~ 56 ANSYS 6 Hz (a) 一阶垂向弯曲 (b) 一阶侧向弯曲 (c) 一阶扭转 (d) 二阶侧向弯曲 (e) 二阶垂向弯曲 (f) 弯扭组合 2 6 Hz

29期xx(copy)

s 15 contact s W Si 3 N μm N / m 4 AFM 95% AFM WPOJ UPOJ WPOJ UPOJ 260 ~ 280 nm WPOJ-DS UPOJ-DS 1 cm 1 cm DEAE-ce

助 剂 改 善 其 止 血 效 果 1 实 验 1.1 原 料 和 试 剂 家 蚕 蛹 经 过 提 取 蛹 油 蛋 白 质 后 剩 余 的 残 渣 ( 主 要 成 分 为 蛹 皮 ), 烘 干 除 杂 粉 碎 后 待 用 ; 壳 聚 糖 ( 成 都 市 科 龙 化 工 试 剂 厂 ), 脱 乙 酰

2 ( 自 然 科 学 版 ) 第 20 卷 波 ). 这 种 压 缩 波 空 气 必 然 有 一 部 分 要 绕 流 到 车 身 两 端 的 环 状 空 间 中, 形 成 与 列 车 运 行 方 向 相 反 的 空 气 流 动. 在 列 车 尾 部, 会 产 生 低 于 大 气 压 的 空 气 流

1556 地 理 科 学 进 展 30 卷 他 关 于 农 村 住 房 结 构 与 抗 震 性 能 的 研 究, 则 多 是 从 工 程 抗 灾 的 角 度, 研 究 某 种 构 造 类 型 的 房 屋, 力 图 找 到 传 统 房 屋 的 结 构 失 误 和 新 建 房 屋 中 存 在 的 问 [

资源 环境 生态 土壤 气象

Microsoft Word - 8 期中文目次.doc

CYTG866.PDF

第二部分

AFLP % 1. 1 DNA 5 RAPD DNA 0. 1g % ~ ph 6 10 DNA 9 40μl TE SSR STS 20ng /L

<4D F736F F D20A6DBB54DACECBEC7B7A7BDD7A655B667A46ABAF42E646F63>

720 () 2009,,, , ( ) (6 8 ), ( 5) ( 6).,,119 E, 150 km., 25. 5, 25,., 5b. 3 ROMS 3. 1 ROMS, S,.,,,, [ 9 ]. 1/ 32 ( ),25, 18 s, 180 s

Technical Acoustics Vol.27, No.4 Aug., 2008,,, (, ) :,,,,,, : ; ; : TB535;U : A : (2008) Noise and vibr

1

(1) (2) (IVI) (2001) (IVI) 50% ~8% 1~30cm (IVI) Study on the Plant Succession of Slopeland Landslide Areas Following H

j.sd

1 SARS SARS SARS SARS SARS SARS S SARS SARS SARS 2 3 DNA SARS 8 NatureScienceEmerging Infectious DiseasesJournal of Virology NatureScienceLance

%

(CIP) : /. :, (/ ) ISBN T S H CI P (2006) CH IJIASH EN GXIAN G YINSHI WEN H U A Y U CHENGY U 1

5 : 7 [11-12] ( ) [13-14] Diessel [12] TPI/GI [15-16] (K3) 1 K1 50% K2 K3 K7 K2 K3 (20%~30%) K2 K10 1 Fig.1 Lithotypes of major coal seams 2.2 1

标题

1 SARS SARS SARS SARS SARS SARS S SARS SARS SARS 2 3 DNA SARS NatureScienceEmerging Infectious DiseasesJournal of Virology NatureScienceLa

Revit Revit Revit BIM BIM 7-9 3D 1 BIM BIM 6 Revit 0 4D 1 2 Revit Revit 2. 1 Revit Revit Revit Revit 2 2 Autodesk Revit Aut

戊 酸 雌 二 醇 片 联 合 宫 颈 注 射 利 多 卡 因 用 于 绝 经 后 妇 女 取 环 的 临 床 效 果 评 价 陆 琴 芬, 等 371 Keywords groups, no removal difficulties and failure, was statistically s

[1-3] (Smile) [4] 808 nm (CW) W 1 50% 1 W 1 W Fig.1 Thermal design of semiconductor laser vertical stack ; Ansys 20 bar ; bar 2 25 Fig

FROM THE EDITOR

<4D F736F F D D DBACEC0F25FD0A3B6D4B8E55F2DB6FED0A32D2D2DC8A5B5F4CDBCD6D0B5C4BBD8B3B5B7FBBAC52E646F63>

,, , 1 (,2006) %, 1. 47,, %, 4. 5, 84 %(,2008a,2008b,2009),,,,, : %, %,?,,,,,,,,, (20

16 31, %, 15 % 1949 ( 1), : 4 8, ;, , , 5, , 5. 05, 400kg,

KJWJ01 Fig. 1 Geological map of the Kajiwa landslide PDH01 60 ~ 70m 80 ~ 90m 80m 2 930m 556m 641m 598m m m 3 2 Ⅰ m 2

, GC/MS ph GC/MS I

(CIP) : /. :, (/ ) ISBN T S H CI P (2006) XIANGPIAOWANLI JIUW ENH UA YU CH ENGYU

标题

30 學 術 論 文 二 復 旦 內 部 圍 繞 鬥 爭 目 標 的 紛 爭 bk bl bm bn

220 20, Vol. 32, No. 6 食 品 科 学 分 析 检 测 荡 体 系 被 用 来 检 测 和 鉴 别 组 成 极 其 复 杂 的 中 药 2007 年 [7-20] 张 泰 铭 等 提 出 了 中 药 非 线 性 化 学 指 纹 图 谱 的 概 念, 并 对 其 形 成 的 原

Microsoft Word 马 姜-出3355

f 2 f 2 f q 1 q 1 q 1 q 2 q 1 q n 2 f 2 f 2 f H = q 2 q 1 q 2 q 2 q 2 q n f 2 f 2 f q n q 1 q n q 2 q n q n H R n n n Hessian

3期


诺贝尔生理学医学奖获奖者

Syngene GBOXHR Retsch MM400 KQ DE Eppendorf Table 1 The sources of Astragali Radix and Hedysarum Radix No. 1 AMM2011_NMGTY Astr

1 引言

<4D F736F F D F B0E6B8DFB1BBD2FDD6B8CAFDC7B0D1D42E646F63>

doc

Microsoft Word 叶然_new_.doc

: 3 : , , ,, , :, ( Ravallion Chen, 2004 ; Khan Riskin, 2005,Yue et al, 2006),, (,2003 ;,200

ti2 guan4 bo1 bo5 huai4 zheng4 hong1 xi2 luo2 ren4

ZUBAN.dvi

Microsoft Word - 刘 慧 板.doc

标题

Microsoft Word 任 辉_new_.doc

LaDefense Arch Petronas Towers 2009 CCTV MOMA Newmark Hahn Liu 8 Heredia - Zavoni Barranco 9 Heredia - Zavoni Leyva

36(4) (2004) Journal of Soil and Water Conservation, 36(4) (2004) earthworms dig soil on surface and prevent plants to grow. But until D

#4 ~ #5 12 m m m 1. 5 m # m mm m Z4 Z5

一次辽宁暴雨过程的诊断及风场反演分析

(Microsoft Word - 11-\261i\256m\253i.doc)

untitled

标题

I (02) (0

1 119 Clark 1951 Martin Harvey a 2003b km 2

<4D F736F F D20B2CCA3BA4542B2A1B6BE BFB9CCE5D3EBB1C7D1CAB0A9B7D6C6DAB5C4B9D8CFB52E646F63>

Transcription:

微生物学报 Acta Microbiologica Sinica 2017, 57(9): 1314-1322 DOI: 10.13343/j.cnki.wsxb.20170026 Extremophile Diversity 极端微生物多样性 高通量测序分析云南腾冲热海热泉真菌多样性 刘开辉 1, 丁小维 1*, 张波 1, 唐小飞 1, 肖敏 2, 鲜文东 2, 李文均 1 陕西理工大学生物科学与工程学院, 陕西汉中 723001 2 中山大学生命科学大学院, 有害生物控制与资源利用国家重点实验室, 广东广州 510275 2* 摘要 : 目的 本研究揭示云南腾冲热海热泉真菌多样性及群落分布格局, 探讨其理化因子对真菌群落结构的影响 方法 利用 Illumina HiSeq2500 高通量测序平台对腾冲热海热泉沉积物宏基因组 ITS 基因进行测序, 并进行生物信息分析 结果 从 5 个热泉样品中共检测到 343484 条有效序列, 包括 5 个真菌门,20 个纲,66 个目 右姐妹泉 (JMQR) 左姐妹泉(JMQL) 蛤蟆嘴泉(M) 桥泉(QQ) 及鼓明泉 (GMQP) 分别以 Agaricales Eurotiales Capnodiales 和 Hypocreales 等为优势目 在属水平上, 共获得 365 个属, 从 JMQR 中检测到 212 个属, 以裂褶菌属 (Schizophyllum) 为最优势 ; 从 JMQL 中挖掘到 197 个属, 以青霉属 (Penicillium) 为最优势类群 ; 从 M 和 QQ 中分别获得 222 个和 270 个属, 均以枝孢属 (Cladosporium) 为最优势 ; 从 GMQP 中发现 179 个属, 以侧齿霉属 (Engyodontium) 丰度最高 NH 4 含量 温度及 ph 影响不同优势真菌的分布, 其中以 ph 与优势类群 (OTU>1%) 结构变化显著性最高 (P=0.05) 结论 云南腾冲热海高温热泉蕴藏着极其丰富的真菌物种, 其不同样品真菌分布具有差异性,pH 可能是影响热泉真菌群落分布的重要因素之一 关键词 : 腾冲热海热泉, 真菌多样性, 高通量测序 真菌是全球生物多样性的重要组成部分, 并在物质循环 能量流动等生态过程中发挥重要作用 据估计, 地球上可能至少存在 510 多万种真菌 [1], 这些资源无疑是人类最宝贵的财富之一 真菌几乎存在于地球的各类生态环境中, 如土壤 水体 粪肥 动植物体内及大气中 由于纯培养技术的不足, 人类所认识的真菌多样性非常有限, 环境中大量的真菌物种仍处于 未知 的状态 由 于全球环境的变化, 真菌物种多样性可能在锐减 为了发现 新的 真菌资源, 一方面利用免培养技 术如 DGGE [2] 高通量测序等探测环境微生物, 另一方面可挖掘极端环境如盐碱湖 [3] 深海 [4] [5] 热泉等环境中蕴藏的真菌资源 在地球上广泛分布着热泉 火山口及堆肥等 基金项目 : 国家自然科学基金 (31100017); 陕西理工大学科研项目 (SLGQD16-06) * 通信作者 Tel/Fax:86-20-84111727;E-mail: 丁小维,dxw518@163.com; 李文均,liwenjun3@mail.sysu.edu.cn 收稿日期 :2017-01-13; 修回日期 :2017-02-19; 网络出版日期 :2017-03-17

刘开辉等 微生物学报, 2017, 57(9) 1315 极端高温环境 这些环境中的嗜热微生物在生命 进化历程中, 由于高温环境的选择作用, 形成独 特的群落结构 遗传特性及生理代谢机制, 因具 有重要的科学和应用价值而成为微生物学研究的 热点之一 [6] 近年来, 研究人员利用高通量测序对 热泉中的细菌 古菌群落结构进行分析, 并对其 微生物群落结构的影响因素进行分析 [7 9] Colman [9] 等报道了美国黄石公园的浮游生物群落结构, 发现 ph 和温度影响上述微生物的群落组成 Hou [10] 等系统研究云南腾冲热海热泉沉积物细菌 古 菌多样性, 认为其优势群体可能受 ph 和温度影 [11] 响 Jiang 等发现除温度外, 腾冲热海姊妹泉中 的硫 砷含量及可溶性氧气均影响细菌群落结构 目前, 利用高通量测序研究热泉微生物主要集中在 细菌 古菌方面, 而有关嗜热真菌仍未见报道 云南腾冲县蕴藏着丰富的地热资源, 境内热 泉遍布各地, 其中以热海区热泉最具代表性 本 研究通过高通量测序分析该地区热泉沉积物真菌 群落格局及物种多样性, 探测到极其丰富的真菌 物种, 并分析影响真菌群落结构的因素, 为再认 识真菌适应高温的 边界 及其资源的发掘提供理 论依据 1 材料和方法 1.1 样品采集 热泉样品于 2014 年 6 月分别采自云南腾冲热 海区桥泉 (QQ) 左姐妹泉 (JMQL) 右姐妹泉 (JMQR) 鼓明泉(GMQP) 及蛤蟆嘴泉 (M) 底沙土沉积物, 并置放于 50 ml 无菌离心管中 采集的样 品迅速放入液氮中冷冻, 用干冰盒带回实验室, 置于 80 C 冰箱中保存, 用于宏基因组 DNA 的提 取 为了探究热泉一些地化因素对真菌分布的影 响, 本研究选择的热泉区其 ph 温度 S 2 及 NH 4 参数均呈现出明显的梯度特征, 并对上述因子进 行原位测定, 数据如表 1 所示 1.2 宏基因组 DNA 提取 约 10.0 g 热泉沙土样品经充分混匀, 并于液 氮中反复研磨 各称取 0.5 g 经研磨的样品, 利用 Fast DNA Spin Kit for Soil (MP Biomedicals) 土壤 DNA 提取试剂盒, 依说明书抽提并纯化热泉沉积 物微生物基因组总 DNA 每个样品的 DNA 完成 3 次重复抽提, 合并后用 1% 的琼脂糖凝胶电泳检测 所提取的基因组 DNA 1.3 ITS 基因的扩增及 Illumina HiSeq 高通量测序 由于真菌 ITS2 区在揭示属间 属内及近似种 的关系中具有十分突出的种类覆盖度和良好的多 样性深度, 能较准确探测出环境样品真菌的群落 结构, 因此本研究选用特异性引物 ITS2-F:GCATC GATGAAGAACGCAGC;ITS2-R:ATATGTAGGA TGAAGAACGYAGYRAA 用于扩增真菌 ITS2 区, 用于扩增不同样品的 ITS2-F 引物融合了 7 个碱基 的 barcode 序列, 以区别不同来源的 ITS 序列 PCR 反应条件 :94 C 5 min;94 C 1 min,53 C 45 s, 72 C 1 min,35 个循环 ;72 C 5 min PCR 反应 体系 (25 μl) : 2 PCR buffer 12.5 μl, dntps (10 mmol/l) 0.5 μl, 基因组 DNA 100 ng,its2-f/r 表 1. 热泉样品地化信息 Table 1. Geochemical information of hot spring samples Sample ph T/ C c(s 2 )/(μg/l) c(nh 4 )/(mg/l) QQ 6.22 75.6 4 0.97 M 7.88 75.1 202 0.13 JMQL 9.00 90.4 4990 1.26 GMQP 9.05 84.5 3480 1.19 JMQR 9.09 83.4 3010 1.66

1316 Kaihui Liu et al. Acta Microbiologica Sinica, 2017, 57(9) (50 μmol/l) 各 1 μl,taq (5 U/μL) 0.5 μl 及 8.5 μl 的 ddh 2 O 同一样品的 PCR 产物合并后经电泳检测, 通过胶回收试剂盒进行纯化 对 PCR 产物进行定量并均一化, 构建 ITS 测序文库, 利用 Illumina HiSeq2500 平台进行高通量测序 1.4 数据分析利用 Mothur 软件分别去除长度小于 200 bp 序列 平均质量分值低于 25 的序列 含有模糊碱基 化学嵌合体及 barcode 序列中含有错配碱基的序列以获得有效序列 (tags), 并进行去冗余处理, 从中挑选出 unique tag 序列 [12] 基于 Naïve Bayesian 方法的分类器 RDP classifier 对序列进行物种注解, 选择 Confidence Threshold (bootstrapping) 参数为 0.5 进行物种分类 根据 97% 的序列相似度对有效序列进行同源比对, 并对操作分类单元 (Operational Taxonomic Units,OTUs) 进行聚类, 计算香农指数 (Shannon Index) ACE 指数 Chao1 指数及文库覆盖率 (Coverage) 根据 OTU 的表达谱数据, 利用 R 语言绘制热图 利用 Canoco 软件分析热泉环境因子对真菌群落的影响, 利用 Bray-Curtis 分析环境因子对真菌群落影响的显著性 2 结果和分析 2.1 测序结果去除不合格序列后,5 个热泉样品中共获得 343484 条有效序列 ( 表 2), 其中从 GMQP JMQL JMQR M 及 QQ 中分别获得 80026 60813 82345 86815 和 33485 条有效序列, 基于 97% 相似度其 OTUs 数量分别为 3003 3669 2689 4099 和 4198, 组装出来的片段中位的片段长度 N50 为 344 360 bp,ace 指数在 10200.5 18936.6,npShannon 指数在 表 2. 热泉样品测序数据统计分析 Table 2. Sequencing data of hot spring samples Sample Tags OTUs N50 ACE npshannon Coverage/% GMQP 80026 3003 356 11268.8 4.1 97.9 JMQL 60813 3669 355 13757.0 4.1 96.5 JMQR 82345 2689 352 10200.5 3.6 98.0 M 86815 4099 360 14869.2 5.6 97.2 QQ 33485 4198 344 18936.6 5.6 92.2 3.6 5.6 之间, 文库覆盖率在 92.2% 98.0%, 说明测得的 OTU 能基本反映所研究的腾冲热海热泉样品真菌的群落组成 2.2 腾冲热泉真菌多样性及群落格局分析利用 RDP classifier 分类器对 OTU 依次进行门 (Phylum) 纲(Class) 目(Order) 科(Family) 及属 (Genus) 水平上的物种分类注释 从云南腾冲热海 5 个热泉样品中共检测到 5 个真菌门, 包括子囊菌门 (Ascomycota) 担子菌门(Basidiomycota) 接合菌类 (Zygomycota) 壶菌门(Chytridiomycota) 及球囊菌门 (Glomeromycota), 其中以子囊菌为优势菌群, 约占 OTUs 总数的 83.6%, 担子菌次之约为 16.1%, 其余类群丰度合计约为 0.1% 从上述样品中检测到 20 个纲, 其中座囊菌纲 (Dothideomycetes) 散囊菌纲 (Eurotiomycetes) 粪壳菌纲(Sordariomycetes) 伞菌纲 (Agaricomycetes) 及酵母纲 (Saccharomycetes) 等为优势纲, 分别占 OTUs 的 30% 26% 22% 13% 和 4% 除了未知真菌群体外, 从 5 个热泉中共检测到 66 个目, 在图 1-A 中, 从 JMQR 中获得 47 个目,Agaricales Eurotiales Pleosporales Capnodiales 及 Trichosphaeriales 为优势类群, 分别占该热泉真菌 OTU 的 27% 25% 21% 6% 及 5%; 在 JMQL 泉中获得 46 个目, 优势类群如 Eurotiales (60%) Hypocreales (8%) Pleosporales (7%) Capnodiales (6%) Agaricales (6%); 热泉 M 中共发现 45 个目, actamicro@im.ac.cn

刘开辉等 微生物学报, 2017, 57(9) 1317 优势类群 Capnodiales 占 28%, 次之为 Eurotiales (16%) Hypocreales (15%) Pleosporales (13%) 及 Saccharomycetales (8%); 在热泉 QQ 中共检测到 60 个目, 其中超过 5% 的优势类群如 Capnodiales (28%) Pleosporales (19%) Eurotiales (13%) Hypocreales (5%) 及 Agaricales (5%); 热泉 GMQP 中发现 41 个目, 以 Hypocreales 为优势群体 (45%), 次之为 Pleosporales (15%) Eurotiales (14%) Capnodiales (9%) 从 5 热泉样品中共检测 365 个真菌属, 不同样品在真菌属水平上分部格局如在图 1-B 所示, 热泉 JMQR 中获得 212 个属, 其中优势类群包括裂褶菌属 (Schizophyllum) 占 27% 青霉属 (Penicillium) 22% 茎点霉属(Phoma) 10% 链格孢属 (Alternaria) 9% 黑孢属(Nigrospora) 5% 及枝 孢属 (Cladosporium) 5%;JMQL 热泉中共检测到 197 个属, 优势属如 Penicillium 占 61% 侧齿霉属 (Engyodontium) 6% Schizophyllum 为 6% 及 Cladosporium 为 5%; 在 M 中获得 222 个属, 以 Cladosporium (21%) Engyodontium (13%) Aspergillus (11%) Saccharomyces (7%) Penicillium (5%) 及 Devriesia (5%) 为优势属 ; 在 QQ 中获得 270 个属, 优势属如 Cladosporium (21%) Engyodontium (13%) 曲霉(Aspergillus) 占 11% 酵母属 (Saccharomyces) 占 7% Penicillium (5%) 及 Devriesia 为 5%; 从 GMQP 中发现 179 个属, 优势属包括 Engyodontium Penicillium 拟盾壳霉属 (Paraconiothyrium) 及 Cladosporium, 分别占微生物 OTU 的 46% 13% 6% 和 5% 上述研究结果表明在高温热泉环境中蕴藏着极其丰富多样的真菌类群 图 1. 热泉样品真菌类群目 (A) 和属 (B) 水平分布图 Figure 1. Order-(A), genus-level (B) distribution of fungi in hot samples from Tengchong, Yunnan. At each rank, fungal sequences that were not identified into anyone in databases were Unclassified ; rare groups shown as others.

1318 Kaihui Liu et al. Acta Microbiologica Sinica, 2017, 57(9) 2.3 腾冲热泉真菌群落结构相似度分析根据 OTU 的丰度数据, 热图可在属水平上将不同丰度的 OTU 分块聚类, 并根据颜色梯度反映不同样品的物种群落相似性 差异性及物种聚类关系展示样本关系 本研究选择不同热泉样的丰度高的前 30 个真菌属, 利用 R 语言绘制热图 ( 图 2),JMQL 和 QQ 类聚为一个分支, 其他 3 个热泉样 GMQP JMQR 及 M 分别同该分支类聚 JMQL 和 QQ 分别以 Penicillium 和 Cladosporium 丰度最高,GMQP 以 Engyodontium 丰度最高, JMQR 和 M 分别以 Schizophyllum 和 Cladosporium 最丰富 除了上述最丰富类群在不同热泉样品中 显著富集外, 其余多数真菌属在不同样品中的丰度具有相似性 2.4 环境因子对热泉真菌群落结构的影响为了分析腾冲热海热泉土壤沉积物理化因子对真菌群落格局的影响, 本研究进行典型对应分析, 结果如图 3 所示 热泉沉积物 S 2 NH 4 温度及 ph 矢量箭头之间夹角为锐角, 表明这 4 种理化因子之间可能具有协同效应 ; 样品 JMQR JMQL 和 M 同 NH 4 相关性最高,QQ 和 GMQP 同 ph 相关性最高 ;Alternaria Nigrospora Penicillium Phoma 和 Schizophyllum 同 NH 4 含量变化相关性最高,Cladosporium 同温度变化相关性最高, Figure 2. 图 2. 不同热泉样品真菌的热图分析 Heatmap of different samples from hot springs in Tengchong, Yunnan. actamicro@im.ac.cn

刘开辉等 微生物学报, 2017, 57(9) 1319 Figure 3. 图 3. 真菌群落与环境因子的典型对应分析 RDA showing associations between environmental variables and fungal genera. Aspergillus Saccharomyces 和 Engyodontium 同样 品 ph 值相关性最高 为进一步分析热泉环境因子对真菌群落结构 的影响, 本研究选择热泉 2 个关键的环境变量 ph 和温度分别同热泉样品中 OTU 丰度大于 1% 的所有 [13] 优势属进行 Bray-Curtis 分析, 结果如图 4 所示, 样品 ph 对优势属群落结构的影响显著性最高 (P=0.05), 其相关性 r=0.491, 然而, 温度对 OTU 丰度大于 1% 的群体结构影响不显著 (P=0.751), 表明 ph 可能是影响腾冲热海热泉真菌群落分布的重要因素之一 3 讨论 微生物在自然界中常以群落的形式存在, 承担着地球物质生物化学循环的重要角色, 因此微生物群落结构及功能特征的探究已成为微生物生态学研究的热点 长期以来, 由于纯培养技术的局限性使我们对微生物种类 丰度及其随环境变化的驱动规律的认识很不全面, 对群落生态功能及其代谢机制方面了解更少 [14] 近年来, 随着高通量测序技术的不断更新, 微生物分子生态学的 图 4. ph 和真菌优势群体的 Bray-Curtis 分析 Figure 4. Bray-Curtis dissimilarity of ph and fungal dominant genera. 研究方法和研究途径也在不断变化 尤其是 Illumina Hiseq 2500 高通量测序平台的成熟和普 及, 使我们能够对环境微生物进行深度测序分析,

1320 Kaihui Liu et al. Acta Microbiologica Sinica, 2017, 57(9) 能灵敏地探测出环境真菌群落组成随环境因子时 空变化而发生的响应 因此, 该技术在微生物群 落结构及多样性研究方面深受学者认可 [15 17] 本研究通过高通量测序分析腾冲热海 5 个热泉真菌群落结构及物种多样性, 共获得 343484 条有效序列, 经分析获得 5 个真菌门 20 个纲, 以子囊菌门为最优势类群 Mouchacca [18] 认为嗜热真菌主要包括 Mucorales Eurotiales Sphaeriales Hyphomycetes 和 Agonomycete 这 5 个目, 然而我们从腾冲热海样品中共检测到 66 个真菌目 365 个属, 以 Schizophyllum Penicillium 及 Phoma 属为 JMQR 热泉的优势类群 ;Penicillium 则主导 JMQL 热泉的真菌群落结构 ;Cladosporium 在 M 热泉中丰度最高,Engyodontium 和 Aspergillus 次之 ;GMQP 泉以 Engyodontium 为最优势真菌 除了丰度最高的类群外,OTU>1% 的所有优势属在不同热泉样品中的组成及其丰度也存在差异, 暗示环境变量驱动热泉真菌的群落格局, 这与环境因子影响热泉原核微生物群落结构的结论一致 [8 11] Redman 和 Chen 等研究发现 Penicillium 和 Aspergillus 是地热生态环境中可培养数目最多的微生物 [5,19] [20] Pan 等采用纯培养技术从腾冲热海 3 个碱性热泉中分离到 8 个真菌属, 如 Chaetomium Coprinopsis Malbranchea Rhizomucor 和 Talaromyces 等, 以 Talaromyces 为最优势群体 其中 Talaromyces 和 Chaetomium 在该研究的 5 个高温热泉中均被检测到, 而 Coprinopsis 仅在 JMQL JMQR 和 QQ 热泉中存在, 但这些真菌的丰度均较低 本研究表明高温热泉中蕴藏着极其丰富 多样的真菌物种资源 大量研究表明, 环境中某些理化因子例如 ph 温度和无机盐等影响微生物群落结构 [21 23] [9] Colman 等发现 ph 和温度是影响黄石公园的浮 游生物群落结构的主要因素 而 ph 温度 硫含量及可溶性氧则影响云南腾冲热泉不同细菌优势群体的分布 [10 11] 除了影响细菌分布外,pH 也驱动真菌的群落结构 [3] 在本研究中, 热泉样品的优势属 Cladosporium 分布受温度变化影响,Alternaria Nigrospora Penicillium Phoma 和 Schizophyllum 同 NH 4 含量有关,Aspergillus Saccharomyces 和 Engyodontium 受 ph 驱动 此外,Bray-Curtis 分析显示 OTU 丰度大于 1% 所有优势属同样品 ph 显著性最高 (P=0.05), 表明 ph 是驱动腾冲热泉真菌分布的一个重要因素 事实上,pH 影响环境有机质的可溶性和酶活性, 影响微生物对有机质分解 转运及其他养分的利用能力 [24], 进而影响微生物生长 繁殖, 导致不同种属的群落组成及丰度不同 [25] 总之, 利用 Illumina Hiseq 2500 高通量测序技术研究云南腾冲热海热泉沉积物中真菌群落结构, 全面揭示出在高温热泉中存在极其丰富 多样的真菌类群, 并探究影响真菌分布的因素, 发 现 ph 温度及 NH 4 含量可能是影响热泉不同真菌群落结构 物种丰度的重要环境因子, 该研究工作将为热泉嗜热真菌资源的开发利用提供理论依据 参考文献 [1] Blackwell M. The fungi: 1, 2, 3... 5.1 million species? American Journal of Botany, 2011, 98(3): 426 438. [2] van Hannen EJ, Mooij W, van Agterveld MP, Gons HJ, Laanbroek HJ. Detritus-dependent development of the microbial community in an experimental system: qualitative analysis by denaturing gradient gel electrophoresis. Applied and Environmental Microbiology, 1999, 65(6): 2478 2484. [3] Liu KH, Ding XW, Wang HF, Zhang XM, Hozzein WN, Wadaan MAM, Lan A, Zhang B, Li WJ. Eukaryotic microbial communities in hypersaline soils and sediments from the alkaline hypersaline Huama Lake as revealed by 454 pyrosequencing. Antonie van Leeuwenhoek, 2014, 105(5): 871 880. [4] Qu J, Liu KH, Ding XW, Deng BW, Chen WQ, Guo QL, Tian actamicro@im.ac.cn

刘开辉等 微生物学报, 2017, 57(9) 1321 XP, Zhang S, Li WJ. Fungal diversity and enzyme activities in marine sediments in the South China Sea. Acta Microbiologica Sinica, 2014, 54(5): 552 562. (in Chinese) 曲佳, 刘开辉, 丁小维, 邓百万, 陈文强, 郭庆兰, 田新朋, 张偲, 李文均. 南海局部海洋沉积物中真菌多样性及产酶活性. 微生物学报, 2014, 54(5): 552 562. [5] Chen KY, Huang D, Liu CC. The mycoflora of hot spring soil in northern Taiwan. Taiwania, 2003, 48(3): 203 211. [6] Xie SC, Yin HF. Progress and perspective on frontiers of geobiology. Science China Earth Sciences, 2014, 44(6): 1072 1086. 谢树成, 殷鸿福. 地球生物学前沿 : 进展与问题. 中国科学 : 地球科学, 2014, 44(6): 1072 1086. [7] Li HZ, Yang QH, Li J, Gao H, Li P, Zhou HY. The impact of temperature on microbial diversity and AOA activity in the Tengchong geothermal field, China. Scientific Reports, 2015, 5: 17056. [8] Xie W, Zhang CL, Wang JX, Chen YF, Zhu YQ, de la Torre JR, Dong HL, Hartnett HE, Hedlund BP, Klotz MG. Distribution of ether lipids and composition of the archaeal community in terrestrial geothermal springs: impact of environmental variables. Environmental Microbiology, 2015, 17(5): 1600 1614. [9] Colman DR, Feyhl-Buska J, Robinson KJ, Fecteau KM, Xu H, Shock EL, Boyd ES. Ecological differentiation in planktonic and sediment-associated chemotrophic microbial populations in Yellowstone hot springs. FEMS Microbiology Ecology, 2016, 92(9): fiw137. [10] Hou WG, Wang S, Dong HL, Jiang HC, Briggs BR, Peacock JP, Huang QY, Huang LQ, Wu G, Zhi XY, Li WJ, Dodsworth JA, Hedlund BP, Zhang CL, Hartnett HE, Dijkstra P, Hungate BA. A comprehensive census of microbial diversity in hot springs of Tengchong, Yunnan Province China using 16S rrna gene pyrosequencing. PLoS One, 2013, 8(1): e53350. [11] Jiang Z, Li P, Jiang DW, Dai XY, Zhang R, Wang YH, Wang YX. Microbial community structure and arsenic biogeochemistry in an acid vapor-formed spring in Tengchong geothermal area, China. PLoS One, 2016, 11(1): e0146331. [12] Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and Environmental Microbiology, 2009, 75(23): 7537 7541. [13] Galand PE, Casamayor EO, Kirchman DL, Lovejoy C. Ecology of the rare microbial biosphere of the Arctic Ocean. Proceedings of the National Academy of Sciences of the United States of America, 2009, 106(52): 22427 22432. [14] Liu KL, Wang JQ, Bu DP, Li D, Yu P, Zhao SG. Current progress in approaches to the study of structure and function diversities of environmental microbial communities. Acta Ecologica Sinica, 2010, 30(4): 1074 1080. (in Chinese) 刘开朗, 王加启, 卜登攀, 李旦, 于萍, 赵圣国. 环境微生物群落结构与功能多样性研究方法. 生态学报, 2010, 30(4): 1074 1080. [15] Zhao S, Zhou N, Zhao ZY, Zhang K, Tian CY. Endophytic bacterial diversity and dynamics in root of Salicornia europaea estimated via high throughput sequencing. Acta Microbiologica Sinica, 2016, 56(6): 1000 1008. (in Chinese) 赵帅, 周娜, 赵振勇, 张科, 田长彦. 基于高通量测序分析盐角草根部内生细菌多样性及动态规律. 微生物学报, 2016, 56(6): 1000 1008. [16] Jin L, Gao XM, Du JH, Wang HX, Guan LM, Wang JC, Wei GW, Qiu XZ. Peat bacterial diversity and community structure in Gahai Lake wetlandin Gan nan. Microbiology China, 2016, 43(11): 2396 2404. (in Chinese) 靳亮, 高学梅, 杜建华, 王洪秀, 关丽梅, 王金昌, 魏国汶, 邱小忠. 尕海湖湿地泥炭细菌多样性分析. 微生物学通报, 2016, 43(11): 2396 2404. [17] Chen L, Xu JM, Feng YZ, Wang JT, Yu YJ, Brookes PC. Responses of soil microeukaryotic communities to short-term fumigation-incubation revealed by MiSeq amplicon sequencing. Frontiers in Microbiology, 2015, 6: 1149. [18] Mouchacca J. Thermophilic fungi: biodiversity and taxonomic status. Cryptogamie Mycologie, 1997, 18(1): 19 69. [19] Redman RS, Litvintseva A, Sheehan KB, Henson JM, Rodriguez R. Fungi from geothermal soils in Yellowstone National Park. Applied and Environmental Microbiology, 1999, 65(12): 5193 5197. [20] Pan WZ, Huang XW, Wei KB, Zhang CM, Yang DM, Ding JM, Zhang KQ. Diversity of thermophilic fungi in Tengchong Rehai National Park revealed by ITS nucleotide sequence analyses. The Journal of Microbiology, 2010, 48(2): 146 152. [21] Kan JJ, Clingenpeel S, Dow CL, McDermott TR, Macur RE, Inskeep WP, Nealson KH. Geochemistry and mixing drive the spatial distribution of free-living Archaea and Bacteria in Yellowstone Lake. Frontiers in Microbiology, 2016, 7: 210. [22] Gao PK, Tian HM, Wang YS, Li YS, Li Y, Xie JX, Zeng B, Zhou JF, Li GQ, Ma T. Spatial isolation and environmental factors drive distinct bacterial and archaeal communities in different types of petroleum reservoirs in China. Scientific Reports, 2016, 6: 20174. [23] Garcia-Pichel F, Loza V, Marusenko Y, Mateo P, Potrafka RM. Temperature drives the continental-scale distribution of key microbes in topsoil communities. Science, 2013, 340(6140): 1574 1577. [24] Persson T, Lundkvist H, Wirén A, Hyvönen R, Wessén B. Effects of acidification and liming on carbon and nitrogen mineralization and soil organisms in mor humus. Water, Air, and Soil Pollution, 1989, 45(1): 77 96. [25] Bååth E, Berg B, Lohm U, Lundgren B, Lundkvist H, Rosswall T, Söderström B, Wirén A. Effects of experimental acidification and liming on soil organisms and decomposition in a Scots pine forest. Pedobiologia, 1980, 20(2): 85 100.

1322 Kaihui Liu et al. Acta Microbiologica Sinica, 2017, 57(9) High-throughput sequencing to reveal fungal diversity in hot springs of Rehai at Tengchong in Yunnan Kaihui Liu 1, Xiaowei Ding 1*, Bo Zhang 1, Xiaofei Tang 1, Min Xiao 2, Wendong Xian 2, Wenjun Li 2* 1 School of Biological Science and Engineering, Shaanxi University of Technology, Hanzhong 723001, Shaanxi Province, China 2 State Key Laboratory of Biocontrol, School of Life Sciences, Sun Yat-Sen University, Guangzhou 510275, Guangdong Province, China Abstract: [Objective] Fungal diversity and community structures were revealed in Rehai hot springs in Tengchong, Yunnan, and effects of sediment geochemical properties on fungal distribution were analyzed. [Methods] Fungal ITS genes were detected in the sediment metagenomes by high-throughput pyrosequencing and bioinformatic analysis. [Results] Total 343484 valid reads from 5 hot springs were obtained, and these reads included 5 phyla, 20 classes, and 66 orders. The hot springs Jiemeiquan right (JMQR), Jiemeiquan left (JMQL), Hamazui (M), Qiaoquan (QQ) and Gumingquan (GMQP) were dominated by communities such as Agaricales, Eurotiales, Capnodiales, and Hypocreales. A total of 365 genera were discovered in the hot springs studied. There were total 212 genera detected in JMQR, 197 in JMQL, 222 in M, 270 in QQ, and 179 in GMQP. Schizophyllum was a dominant genus of fungal community in JMQR, and Penicillium was highly frequent in JMQL. Cladosporium dominanted fungal community in M and QQ, and Engyodontium was highly rich in GMQP. The geochemical variables such as NH 4, temperature, and ph, exerted effects on distribution of different dominant fungi in the hot springs. The ph in the hot springs had significant effect on fungal community composition. [Conclusion] Rehai hot springs in Tengchong, Yunnan, harbor surprisingly fungal diversity. Fungi in different samples had distinct community structures. The ph in the hot springs could be one of the factors determining the fungal community. Keywords: Tengchong hot spring, fungal diversity, pyrosequencing ( 本文责编 : 张晓丽 ) Supported by the National Natural Science Foundation of China (31100017) and by the Shaanxi University of Technology Project (SLGQD16-06) * Corresponding author. Tel/Fax: 86-20-84111727; E-mail: Xiaowei Ding, dxw518@163.com; Wenjun Li, liwenjun3@mail.sysu.edu.cn Received: 13 January 2017; Revised: 19 February 2017; Published online: 17 March 2017 李文均, 中山大学生命科学学院百人计划 / 珠江学者特聘教授, 博士生导师 1996 年于河南师范大学生物系获学士学位 ;1999 年于云南大学生命科学学院获硕士学位 ;2002 年毕业于中国科学院沈阳应用生态研究所获博士学位 ;1999 2014 年曾在云南大学云南省微生物研究所工作 现任国际原核微生物系统学委员会 (nternational Committee on Systematics of Prokaryotes,ICSP) 国际执委 ; 伯杰氏国际系统微生物学会 (BISMiS) 创始会员 兼任中国青年科技工作者协会常务理事 中国微生物学会微生物教学工作委员会副主任委员 中国微生物学会普通微生物专业委员会委员 中国微生物学会分子微生物学与生物工程专业委员会委员 微生物资源专业委员会委员 Antonie van Leeuwenhoek Frontiers in Microbiology BioMed Research International The Scientific World Journal Journal of Arid Land BISMiS Bulletin 微生物学报 微生物学杂志 等杂志编委及多个国内外期刊的通讯评审专家 主要从事高温 高盐等极端环境微生物分类及系统学, 微生物生态学, 极端酶学研究, 实验室主页 :http://www.liactlab.com/index.html actamicro@im.ac.cn