生态学报  2017, Vol. 37 Issue (6): 1749-1756

文章信息

韩丽丽 , 于丹婷 , 贺纪正
HAN Lili , YU Danting , HE Jizheng .
土壤病毒生态学研究方法
Research methods for soil viral ecology
生态学报. 2017, 37(6): 1749-1756
Acta Ecologica Sinica. 2017, 37(6): 1749-1756
http://dx.doi.org/10.5846/stxb201511122292

文章历史

收稿日期: 2015-11-12
网络出版日期: 2016-08-02
土壤病毒生态学研究方法
韩丽丽 1, 于丹婷 1,2, 贺纪正 1     
1. 中国科学院生态环境研究中心, 城市与区域生态国家重点实验室, 北京 100085;
2. 中国科学院大学, 北京 100049
摘要: 病毒是地球上最丰富的生物实体,每克土壤中可包含数以亿计的病毒,它不仅影响土壤中其它微生物的群落组成、土壤元素的生物地球化学循环,还会影响土壤微生物的物种进化,甚至影响植物、动物和人体健康。目前人们对土壤中病毒的种类及丰度、分布特征以及功能引起的生态环境效应还知之甚少。在概述病毒生态学研究方法的基础上,对土壤病毒的提取、纯化、定量及分子生态学方法等基本流程进行了比较分析,以期建立一套快速简便、高效稳定的适用于土壤病毒研究的方法,并用于研究土壤病毒的多样性及分布特征,探讨病毒在环境中的生存和传播机制,为土壤病毒的防控及开发利用提供支撑。
关键词: 土壤病毒     病毒生态学     宏病毒组学     分子生物学方法    
Research methods for soil viral ecology
HAN Lili 1, YU Danting 1,2, HE Jizheng 1     
1. State Key Laboratory of Urban and Regional Ecology, Research Center for Eco-Environmental Sciences, Chinese Academy of Sciences, Beijing 100085, China;
2. University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: Viruses are the most abundant biological entities on the planet, and can reach 1010 viral particles per gram of soil. Viruses affect the composition of microbial communities, mediate soil biogeochemical cycles, regulate soil microbial evolution, and threaten the health of plants and animals, including humans. However, there is limited knowledge on the abundance and distribution patterns of viruses and the manner in which they affect ecological processes. This study focuses on the comparative approaches to study soil viruses, including their extraction, purification, quantification, and methods pertaining to molecular ecology. The development of reliable and highly efficient methods to evaluate soil viral particles is desirable, because it aids the study of soil virus diversity and distribution. Therefore, these methods will improve our understanding of viral propagation mechanisms in soils, and will provide information on how to control viruses to ensure adequate public health and safety.
Key words: soil virus     viral ecology     metaviromics     molecular biological methods    

病毒是地球上最丰富的生物实体, 对土壤中营养元素的生物地球化学循环、微生物群落组成以及土壤基因资源的多样性发挥着重要作用。首先, 病毒能够侵染细菌、真菌、蓝细菌等微生物, 影响它们的生理代谢及死亡率, 从而影响生态系统中微生物的群落结构及多样性;其次, 由于病毒侵染导致的细胞裂解释放的有机物影响环境中其它生物的养分利用, 从而影响环境中元素的生物地球化学循环;第三, 噬菌体引起宿主细胞裂解并释放出DNA, 释放到环境中的DNA被其它微生物通过转录方式吸收、整合, 引起微生物的遗传变异, 推动微生物种群的进化。同时, 土壤的高度异质性及其多样的生物环境为病毒提供了丰富而多样的寄生环境, 更利于病毒长期生存和繁殖。

尽管病毒在土壤环境中发挥着重要的作用, 但由于其个体微小、基因组含量低、缺乏通用基因、进化速度快生命周期短等特点, 使得土壤病毒生态学的研究进展缓慢。病毒生态学研究方法经历了表型分析、单基因分析、基因组分析及宏基因组分析等几个重要阶段。如20世纪90年代, 科学家们利用显微技术对病毒和细菌丰度的地域差异进行过大量研究, 发现在撒哈拉及纳米布沙漠这类极端环境中均存在着非常独特的噬菌体类群[1-2]。随着技术的发展, 更多的研究结合了表型和分子生物学方法, 如将透射电子显微镜 (TEM:Transmission Electron Microscope) 和脉冲场凝胶电泳 (PFGE:Pulse Field Gel Electrophoresis)、随机扩增多态性DNA (RAPD:Random Amplified Polymorphic DNA analysis) 方法相结合。这类方法可有效地用于研究单个分离的噬菌体, 也可对噬菌体的群体组成进行研究。这类基于基因组水平的研究方法还可用于评价不同噬菌体在环境中的分布, 如研究发现某些海洋噬菌体基因组只在特定的地点出现, 而其它一些噬菌体分布却很广[3]。有研究利用RAPD技术研究还发现弧菌噬菌体的丰度和地理分布没有直接关系, 而是受生境类型所驱动, 如与海水、沉积物以及无脊椎动物等生境相关[4]

测序技术的出现和日益广泛的应用, 突破了探寻环境微生物黑箱的瓶颈。关于病毒基因组的测序始于1977年, 第1株被测序的病毒是侵染大肠杆菌的ΦΧ174。2002年, Brcitbart等[5]首次通过鸟枪测序法研究海洋病毒的宏基因组。2005年, Edwards等[6]提出了病毒宏基因组的概念。Hurwitz等[7]采用比较基因组学和网络分析方法分析了海洋病毒群落的生态驱动模型。Bolduc等[8]通过网络分析方法揭示了黄石酸性热泉中病毒的群体组成及相互关系。截止至2015年10月NCBI数据库中已收录了4885株病毒的全基因组序列。

与海水及热泉样品相比, 土壤样品具有异质性高、成份复杂等特点, 从土壤中分离提取高纯病毒面临更多的困难, 本文在对病毒生态学研究进展综述的基础上[9], 进一步对病毒生态学研究方法作一简要概述。在对海洋及其它环境病毒的提取、纯化、定量及分子生态学方法进行比较优化的基础上, 重点建立一套快速简便、高效稳定的适用于土壤病毒生态学研究的方法 (图 1), 主要包括病毒样品的富集、定量、核酸提取、指纹图谱分析及宏基因组测序等几个部分。

图 1 土壤病毒生态学研究方法流程 Fig. 1 Research method process of soil viral ecology
1 土壤中病毒的提取及富集 1.1 土壤病毒的提取

病毒对土壤的吸附与病毒种类、土壤pH值、温度、含水量、土壤结构等环境因子密切相关。土壤病毒学研究的第一步是选择合适、高效的提取方法高效地提取出吸附于土壤颗粒表面的病毒粒子。最常用的土壤病毒浸提液有4种[10], 包括10%的牛肉膏、250 mmol/L甘氨酸溶液、10 mmol/L焦磷酸钠和1%的柠檬酸钾溶液。比较几种浸提液对土壤病毒的提取效果, 发现250 mmol/L甘氨酸溶液[11]和7%的牛肉膏提取液[12]最适用于提取污泥中的多种噬菌体, 但是浸提出的病毒溶液难以过滤及进行后续处理, 浸提的病毒溶液粘稠, 无法制片染色, 也不可用荧光显微镜进行计数;1%的柠檬酸钾溶液提取农田土壤病毒效率最高, 并且可通过荧光显微镜直接计数;10 mmol/L焦磷酸钠浸提效率略次于柠檬酸钾浸提液, 缺点是浸提液使用荧光显微镜观察时易形成团块, 无法计数。综上, 对于土壤病毒浸提液以250 mmol/L甘氨酸和1%的柠檬酸钾溶液最为适用, 具体选择还要依据待提取病毒类型、土壤的理化性质而定。

1.2 土壤病毒溶液的分离、纯化及浓缩 1.2.1 切向流过滤法

切向流过滤 (TFF:Tangential-flow filtration) 是指液体流动方向与过滤方向垂直的过滤形式, 这种方法与常规过滤相比, 具有液体流动在过滤介质表面产生剪切力, 减小滤饼层颗粒物的堆积, 不易堵塞, 可保证稳定的过滤速度等优点。TFF主要用于从大量环境样品中分离、浓缩病毒颗粒, 设备及流程见图 2。TFF系统使用中空纤维膜滤柱, 本文以孔径0.2 μm的微滤柱和孔径为100 kDa的超滤柱为例 (图 2), 介绍TFF系统过滤的原理及流程。孔径0.2 μm的微滤柱主要用于过滤去除直径大于0.2 μm的细菌及原生生物等, 而直径小于0.2 μm的病毒则能通过滤柱滤出 (图 2), 通过收集过滤液即可得到较纯的病毒溶液。但此时的病毒溶液浓度很低, 无法用于下一步研究, 因此需要再通过100 kDa的超滤柱浓缩。更换滤柱后, 直径大于100 kDa的病毒被滤柱拦截, 返回到回流端, 而水溶液则能通过滤柱滤出 (图 2), 通过不断回流过滤, 减小病毒溶液体积, 直至残余液体体积小于100 mL。

图 2 TFF过滤系统及流程[13] Fig. 2 TFF filtration system and procedure[13]
1.2.2 聚乙二醇沉淀法

聚乙二醇沉淀法 (PEG:Polyethylene glycol) 主要是利用PEG破坏病毒外壳蛋白质分子表面的水化层而使病毒发生沉淀的方法, 常用于病毒的提纯浓缩, 它可将TFF浓缩的病毒溶液 (50—100 mL) 中的病毒沉淀。首先向浓缩的病毒溶液中加入终浓度为0.5 mol/L的NaCl溶液, 再加入等量的10% PEG6 000, 4℃过夜放置后, 8000 r/min离心30 min, 即可获得病毒沉淀, 再用缓冲液将病毒溶解成1—2 mL, 用于下一步病毒的CsCl密度梯度离心或DNA的提取。

1.2.3 CsCl密度梯度离心法

由于环境中病毒的大小、质量不一, 密度也不一样, 密度梯度离心法既可用于病毒的浓缩, 也可用于病毒颗粒的分离纯化。每种病毒的大小不变, 对应不同介质的密度层, Thurber等[13]在2009年的综述文献中作了较全面的归纳。该方法主要适用于纯病毒的分离, 缺点是仪器设备价格昂贵, 一般实验室难以配备。

2 定量技术 2.1 荧光显微计数

病毒定量目前最常用的方法就是荧光显微 (EFM:Epifluorescence microscope) 计数法, 该方法利用真空抽滤法将稀释的病毒溶液中的病毒颗粒加载在滤膜上, 通过高效的荧光染料染色后, 将荧光信号分子和病毒DNA双链结合, 在荧光显微镜下可观察到发出荧光的病毒DNA, 利用目镜网格尺计数后推算出每克土壤样品中病毒颗粒的数量。常用的荧光标记物有Yo-Pro-1、SYBR Green I和SYBR Gold等嵌入性染料[14], 该方法与透射电镜相比, 价格较低、操作简便且重复性强, 适用于土壤样品中未知病毒颗粒的定量研究。利用EFM直接计数有一点值得注意的是, 由于可能对一些非病毒粒子染色, 会高估病毒的丰度[15]。另外, 与海水病毒样品相比, 土壤中有较高含量的腐殖酸, 易导致较高的荧光背景干扰, 这也是目前土壤病毒EFM计数方法无法避免的缺陷。

2.2 透射电镜分类计数

在分子生物学技术发展以前, 透射电子显微镜 (TEM) 技术是作为病毒生态学研究的最主要手段, 它不仅可应用于病毒定量, 还可对病毒进行分类。具体方法:取50 g土壤样品经浸提、过滤纯化后, 取1 mL的纯化提取物加入9 mL的灭菌去离子水, 在Formvar-coated 650铜丝网网格上旋转, 夹起铜网用滤纸吸取边上多余液体, 立刻将铜网倒扣在1%醋酸双氧铀负染液滴上, 静置1—2 min, 再夹起铜网吸去多余染液, 自然晾干。在透射电镜下先放大4000倍观察全貌, 再逐步放大到40000—85000倍观察病毒形态, 对病毒进行分类计数。这种方法直接准确, 在荧光显微镜下不能观察到的病毒可以通过透射电镜观察到, 且结果精确度高, 但是操作复杂, 对仪器要求较高。

2.3 流式细胞仪计数

流式细胞仪 (Flow Cytometer, FCM) 是以激光为光源、检测生物学颗粒的仪器。流式细胞术是一种能够对单个细胞、分子、生物颗粒进行多参数快速检测的新型分析和分选技术。将待测细胞或微粒进行荧光染色后制成悬液标本, 在一定气体压力下将待测样品压入流动室, 用不含细胞或微粒的缓冲液 (鞘液) 在高压下从鞘液管喷出, 鞘液管入口方向与待测细胞或微粒流呈一定角度, 使鞘液包绕这细胞或微粒高速流动, 形成一个圆形的流速 (鞘流), 待测细胞在鞘液的包裹下单行排列, 一次通过流式细胞仪的检测区域。这种方法测定灵敏度高, 对病毒浓度较高的样品测试结果目前是最准确的, 但无法测试浓度较低的病毒样品[16]。同荧光显微计数方法一样, 流式细胞仪也只能检测被荧光标记的双链DNA病毒, 对单链DNA病毒和RNA病毒仍无法检测。

2.4 病毒计数仪计数

病毒计数仪是为了检测病毒粒子而重新组装的利用特殊液体控制、双重染色标记的流式细胞分析仪。对病毒的定量可以在5 min内完成, 且可以转变成全自动分析平台。主要原理是通过联合染色系统, 利用荧光染料对病毒基因组和包膜蛋白染色, 通过流动室时被激光激发后, 标记的蛋白和核酸分别产生黄色和红色的荧光信号同时通过不同荧光通道被捕捉, 只有同时产生这2种信号才计算为一个完整的病毒颗粒。这种方法相对于噬斑滴度计数法、荧光显微计数法等具有分析速度快、成本低、样品制备简单、灵敏度高等特点, 缺点在于仪器成本高、检测的病毒种类有限, 只能检测人体及动物体内的20多种病毒, 检测结果不能代表环境中所有的病毒。

3 病毒核酸的提取及分子生态学方法 3.1 核酸提取

对富集的病毒浓缩液通过EFM或TEM方法检测病毒纯度, 看是否有其他微生物细胞污染, 对纯化后的病毒提取核酸。裂解病毒前, 使用DNase I降解土壤溶液中游离的DNA, 并使用16S rRNA引物PCR检测无细菌DNA污染后, 进行病毒核酸的提取。病毒核酸提取过程中最重要的是每一步操作都要做到不受病毒及其他微生物的污染。病毒核酸提取有多种试剂盒, 如Mobio公司的PowerViralTM Environmental RNA/DNA Isolation Kit Sample、Qiagen公司的QIAamp Viral RNA Mini Kit、QIAampViral RNA Mini Kit等, 或者使用甲酰胺或CTAB结合酚仿抽提法[13, 17]大量提取, 不同类型土壤适用的方法也不一样, 提取效率也不尽相同。提取的DNA如果浓度较低可通过Phi29DNA聚合酶扩增, RNA样品则需通过反转录酶反转成cDNA后进行后续实验。

3.2 目的基因多样性

病毒生态学研究之所以远远滞后于细菌和真核生物, 最重要的原因之一就是病毒之间缺乏通用的标记基因, 但对于特定种类的病毒, 科学家们还是发现了它们之间存在一些相对保守的标记基因[18], 并通过这些标记基因序列研究特定病毒群体的多样性。常用的标记基因有编码外壳蛋白基因, 壳组装蛋白基因, 尾鞘蛋白基因, 辅助代谢基因以及几种聚合酶基因等。目前研究较多的是T4型噬菌体中编码主要壳体蛋白的g23基因。2005年Filée等[19]比较了16株T4型噬菌体单以g23基因与以基因组中23个基因组合绘制的多基因系统发育树, 发现这2种方法构建的结果非常相似, 证明g23基因可以作为T4型噬菌体类群划分的分子基础, 可以较好地表征复杂环境中T4型噬菌体基因多样性。g20基因编码噬藻体壳组装蛋白, 也已被广泛用于Myoviridae家族病毒多样性的研究, 研究覆盖了稻田土壤, 稻田水体, 海洋及其它多种淡水环境[20-24]。针对不同环境的g20基因多样性研究也发现, 不同环境都有其独特的类群, 说明g20基因的分布与其存在的环境有一定的关系[25]。蓝藻噬菌体中的光合作用相关基因psbApsbD也被广泛应用于蓝藻噬菌体多样性的研究[26-28], psbA基因的系统发育研究表明它能分辨不同环境不同宿主菌, 但却和地理位置没有相关性[29]。其它标记基因如主要的壳粒蛋白基因mcp[30-33]以及磷酸盐饥饿状态相关的phoH基因[34-35]在海洋噬菌体的研究中也有应用, 它们在噬菌体家族中的分布并没有严格的限制。在噬菌体的核质中也有一些和DNA合成相关的基因被作为标记基因, 如Podoviridae家族中T7型噬菌体中的DNA聚合酶基因polA[36-37], Myoviridae家族中的T4型噬菌体的DNA聚合酶基因g43[38], Family B家族中核质大DNA (NCLDVs) 中的DNA聚合酶基因polB[38-39], 以及一些RNA病毒中的RNA聚合酶基因[40-52]

3.3 指纹图谱技术

随着分子生物学技术的快速发展, 更多的病毒学研究引入了指纹图谱技术, 如脉冲场凝胶电泳、随机扩增多态性DNA等, 并将其与TEM方法相结合, 可对噬菌体的群体组成进行研究。这种基于基因组水平的研究方法不仅显示病毒的多样性, 还可结合表型横向比较不同噬菌体在环境中的分布, 如利用TEM和PFGE结合研究海洋中蓝藻噬菌体的基因组, 虽然发现大部分分离的蓝藻噬菌体和长尾噬菌体表型很相似, 但它们的基因组大小则分布在45—125 kb之间[43];利用TEM和RAPD方法研究海洋弧菌噬菌体的时空分布;利用RAPD技术研究还发现弧菌噬菌体的丰度和地理分布没有直接关系, 而是受生境类型所驱动等。

3.4 宏病毒组测序

宏病毒组测序 (Viromics) 是利用宏基因组学研究环境中病毒多样性的高通量测序方法。目前宏基因组测序采用的方法有传统的Sanger/鸟枪法和二代测序技术。应用最为广泛的莫过于第二代测序的三大测试平台, 即罗氏公司的454 GS FLX测序平台、ABI公司的SOLiD测序平台和Illumina公司的Solexa Genome Analyzer测序平台。测序的过程大致如下:提取的病毒DNA (或RNA反转成cDNA), 经质检合格后, 利用超声将DNA打断, 片段化的DNA经过胶回收后, 在两端加接头, 再通过不同的方式将每个片段结合到微珠或芯片上, 形成几百万个同时反应的微型反应池, 每个微型反应池中的DNA片段通过PCR扩增产生的多拷贝作为测序的模板[44];随后通过酶延伸或寡核苷酸连接反应同时对这些模板进行测序。

获得海量的序列数据后, 最困难的问题在于数据分析处理。原始序列经过质检和优化后, 获得的Reads拼接组装后, 利用MetaGene进行基因预测;在此基础上, 可以获得序列组成、物种组成、功能组成和群落特征等重要的宏基因组信息;数据分析的下一个关键步骤是根据相似度或序列组成对序列进行归类, 基于相似度的方法将序列与KEGG、COG、SEED和ACLAME等功能数据库比对, 获得病毒基因的功能注释结果及丰度信息。通过病毒分类或功能基因的丰度差异比较分析、显著性差异特征分析, 为后续研究提供方向。这种方法的局限性在于价格相对昂贵, 较依赖现有数据库, 而病毒的未知序列及开放阅读框数据库比例很高, 大部分的参考序列不能被鉴定。

和传统的病毒研究方法如电子显微镜技术, 基因指纹图谱或靶基因测序等相比, 宏病毒组是唯一一个能反映样点病毒群体多样性和功能的方法[45]。利用宏基因组学技术研究环境中病毒群落的多样性, 避免了病毒中没有通用的标记基因的缺点, 可以用来作为系统发育的指标, 评估环境中病毒的多样性[46]。但环境病毒学的研究还处于起步阶段, 病毒宏基因组学研究目前面临的最大问题是序列注释困难, 大约只有10%的序列能从已知数据库中比对上[47]

本文对已发表的土壤病毒研究文献进行了归纳整理, 表 1列出了不同地区土壤中优势病毒的分布及使用的研究方法。对土壤病毒的研究还是以宏病毒组学方法为主, 比较不同地区的病毒分布, 发现方法对病毒分类的影响尤为重要, 同一实验室研究不同土壤中的病毒分布得出的结果可能相似, 而不同实验室研究同种类型的土壤病毒分布其结果差异则较大, 这可能是由于地理分布的影响, 也可能是方法的选择及操作对提取的病毒存在偏好性。有研究指出Phi29DNA聚合酶对扩增环状DNA具有偏好性, 可能导致具有环状结构的ssDNA被选择性的放大, 这也许是多数土壤中ssDNA病毒占主导优势的一个重要原因[48]。另外, 在已发表的土壤病毒宏基因组研究中, 还未见有RNA病毒的相关报道, 这可能是由于病毒的基因组核酸含量低, 对病毒RNA的提取更为困难, 导致RNA病毒的研究仍存在很多困难, 问题的解决还有待全基因组扩增技术和测序技术的发展。

表 1 不同地区土壤中的优势病毒及研究方法 Table1 Dominant viruses in different areas and research methods
土壤类型
Soil type
地区
Area
代表性病毒
Typical viruses
基因组类型
Genome type
研究方法
Method
参考文献
References
特华拉州土壤
Delaware Soils
美国
USA
有尾噬菌体
Caudovirales
双链DNA
dsDNA
电镜
TEM
[49]
沙漠土壤
Desert
美国
USA
浮游放线菌噬菌体
Actinoplanes ΦAsp2
双链DNA
dsDNA
病毒基因克隆文库
Viral clone libraries
[50]
草原土壤
Prairie
美国
USA
肌尾噬菌体
Myoviridae
热带雨林土壤
Rainforest
秘鲁
Peru
长尾噬菌体
λ-like siphophage
水稻土
Paddy soil
韩国
Korea
微小噬菌体
Microviridae
单链DNA
ssDNA
宏病毒组学
iromics (MDAH)
[48]
环状病毒
Circoviridae
单链DNA
ssDNA
宏病毒组学
iromics (MDAX)
表层土
Surface soil
南极
Antarctica
长尾噬菌体
Siphoviridae
双链DNA
dsDNA
宏病毒组学
iromics
[51]
石下土
Hypolith
长尾噬菌体
Siphoviridae
双链DNA
dsDNA
纳米布沙漠土
Namib desert
南非
Southern Africa
长尾噬菌体
Siphoviridae
双链DNA
dsDNA
宏病毒组学
iromics
[17]
沿海沙地土
Machair
苏格兰
Scotland
微小噬菌体
Microviridae
单链DNA
ssDNA
宏病毒组学
iromics
[52]
棕壤土
Brown earth
微小噬菌体
Microviridae
单链DNA
ssDNA
MDAH:DNA聚合酶加热处理;MDAX:DNA聚合酶不加热处理
4 研究展望

病毒组学方法无疑是病毒生态学研究中最全面、最直接、获得数据量最多的方法。面对如此庞大的信息量, 仅分析病毒基因组信息和病毒基因多样性信息是远远不够的, 在生物信息学水平上, 人们还需要开发更多的分析软件进一步挖掘病毒群体的多样性及其与环境因子之间的相互关系;在进化水平上, 病毒引起基因水平转移及其它转座子的作用机制还有待研究, 微生物基因组和病毒基因组之间的相互作用关系仍存在多种可能性;从技术的发展水平上, 生物信息学和多学科交叉技术的发展, 有助于推动病毒群体动力学和宏病毒组数据数学模型的建立, 这也将是未来病毒生态学研究的重要发展趋势。另外, 目前的研究主要集中在双链DNA病毒, 对于单链DNA病毒及RNA病毒的研究还需依赖方法的发展。最后, 有针对性地开展土壤病毒群体的生态分析, 建立相关技术体系及土壤病毒资源库, 将有助于对有益病毒开发利用及对关键有害病毒进行控制。

参考文献
[1] Prestel E, Salamitou S, DuBow M S. An examination of the bacteriophages and bacteria of the Namib desert. The Journal of Microbiology, 2008, 46(4): 364–372. DOI:10.1007/s12275-008-0007-4
[2] Prigent M, Leroy M, Confalonieri F, Dutertre M, DuBow M S. A diversity of bacteriophage forms and genomes can be isolated from the surface sands of the Sahara Desert. Extremophiles, 2005, 9(4): 289–296. DOI:10.1007/s00792-005-0444-5
[3] Steward G F, Montiel J L, Azam F. Genome size distributions indicate variability and similarities among marine viral assemblages from diverse environments. Limnology and Oceanography, 2000, 45(8): 1697–1706. DOI:10.4319/lo.2000.45.8.1697
[4] Comeau A M, Chan A M, Suttle C A. Genetic richness of vibriophages isolated in a coastal environment. Environmental Microbiology, 2006, 8(7): 1164–1176. DOI:10.1111/emi.2006.8.issue-7
[5] Breitbart M, Salamon P, Andresen B, Mahaffy J M, Segall A M, Mead D, Azam F, Rohwer F. Genomic analysis of uncultured marine viral communities. Proceedings of the National Academy of Sciences of the United States of America, 2002, 99(22): 14250–14255. DOI:10.1073/pnas.202488399
[6] Edwards R A, Rohwer F. Viral metagenomics. Nature Reviews Microbiology, 2005, 3(6): 504–510. DOI:10.1038/nrmicro1163
[7] Hurwitz B L, Westveld A H, Brum J R, Sullivan M B. Modeling ecological drivers in marine viral communities using comparative metagenomics and network analyses. Proceedings of the National Academy of Sciences of the United States of America, 2014, 111(29): 10714–10719. DOI:10.1073/pnas.1319778111
[8] Bolduc B, Wirth J F, Mazurie A, Young M J. Viral assemblage composition in Yellowstone acidic hot springs assessed by network analysis. The ISME Journal, 2015, 9(10): 2162–2177. DOI:10.1038/ismej.2015.28
[9] 韩丽丽, 贺纪正. 病毒生态学研究进展. 生态学报, 2016(16). DOI:10.5846/stxb201501210170
[10] Williamson K E, Wommack K E, Radosevich M. Sampling natural viral communities from soil for culture-independent analyses. Applied and Environmental Microbiology, 2013, 69(11): 6628–6633.
[11] Araujo R M, Lasobras J, Lucena F, Jofre J. Methodological improvements for the recovery of Bacteroides fragilisphages and coliphages from environmental samples. Water Science & Technology, 1993, 27(3/4): 119–122.
[12] Monpoeho S, Maul A, Mignotte-Cadiergues B, Schwartzbrod L, Billaudel S, Ferré V. Best viral elution method available for quantification of enteroviruses in sludge by both cell culture and reverse transcription-PCR. Applied and Environmental Microbiology, 2001, 67(6): 2484–2488. DOI:10.1128/AEM.67.6.2484-2488.2001
[13] Thurber R V, Haynes M, Breitbart M, Wegley L, Rohwer F. Laboratory procedures to generate viral metagenomes. Nature Protocols, 2009, 4(4): 470–583. DOI:10.1038/nprot.2009.10
[14] Suttle C A, Fuhrman J A. Enumeration of virus particles in aquatic or sediment samples by epifluorescence microscopy//Wilhelm S W, Weinbauer M G, ed. Manual of Aquatic Viral Ecology. Suttle CA: ASLO, 2010, 15: 145-153.
[15] Suttle C A. Marine viruses-major players in the global ecosystem. Nature Reviews Microbiology, 2007, 5(10): 801–812. DOI:10.1038/nrmicro1750
[16] Brussaard C P D, Marie D, Bratbak G. Flow cytometric detection of viruses. Journal of Virological Methods, 2000, 85(1/2): 175–182.
[17] Adriaenssens E M, Van Zyl L, De Maayer P, Rubagotti E, Rybicki E, Tuffin M, Cowan D A. Metagenomic analysis of the viral community in Namib Desert hypoliths. Environmental Microbiology, 2015, 17(2): 480–595. DOI:10.1111/1462-2920.12528
[18] Adriaenssens E M, Cowan D A. Using signature genes as tools to assess environmental viral ecology and diversity. Applied and Environmental Microbiology, 2014, 80(15): 4470–5480. DOI:10.1128/AEM.00878-14
[19] Filée J, Tétart F, Suttle C A, Krisch H M. Marine T4-type bacteriophages, a ubiquitous component of the dark matter of the biosphere. Proceedings of the National Academy of Sciences of the United States of America, 2005, 102(35): 12471–12476. DOI:10.1073/pnas.0503404102
[20] Fuller N J, Wilson W H, Joint I R, Mann N H. Occurrence of a sequence in marine cyanophages similar to that of T4 g20 and its application to PCR-based detection and quantification techniques. Applied and Environmental Microbiology, 1998, 64(6): 2051–2060.
[21] Zhong Y, Chen F, Wilhelm S W, Poorvin L, Hodson R E. Phylogenetic diversity of marine cyanophage isolates and natural virus communities as revealed by sequences of viral capsid assembly protein gene g20. Applied and Environmental Microbiology, 2002, 68(4): 1576–1584. DOI:10.1128/AEM.68.4.1576-1584.2002
[22] Dorigo U, Jacquet S, Humbert J F. Cyanophage diversity, inferred from g20 gene analyses, in the largest natural lake in France, Lake Bourget. Applied and Environmental Microbiology, 2004, 70(2): 1017–1022. DOI:10.1128/AEM.70.2.1017-1022.2004
[23] Matteson A R, Loar S N, Bourbonniere R A, Wilhelm S W. Molecular enumeration of an ecologically important cyanophage in a Laurentian Great Lake. Applied and Environmental Microbiology, 2011, 77(19): 6772–6779. DOI:10.1128/AEM.05879-11
[24] McDaniel L D, delaRosa M, Paul J H. Temperate and lytic cyanophages from the Gulf of Mexico. Journal of the Marine Biological Association of the United Kingdom, 2006, 86(3): 517–527. DOI:10.1017/S0025315406013427
[25] 荆瑞勇, KimuraM, 王光华. 不同自然环境下噬藻体g20基因多样性研究进展. 微生物学报, 2013, 53(11): 1149–1157.
[26] Mann N H, Cook A, Millard A, Bailey S, Clokie M. Marine ecosystems: bacterial photosynthesis genes in a virus. Nature, 2003, 424(6950): 741–741.
[27] Millard A, Clokie M R J, Shub D A, Mann N H. Genetic organization of the psbAD region in phages infecting marine Synechococcus strains. Proceedings of the National Academy of Sciences of the United States of America, 2004, 101(30): 11007–11012. DOI:10.1073/pnas.0401478101
[28] Zeidner G, Bielawski J P, Shmoish M, Scanlan D J, Sabehi G, Béjà O. Potential photosynthesis gene recombination between Prochlorococcus and Synechococcus via viral intermediates. Environmental Microbiology, 2005, 7(10): 1505–1513. DOI:10.1111/emi.2005.7.issue-10
[29] Chénard C, Suttle C A. Phylogenetic diversity of sequences of cyanophage photosynthetic gene psbA in marine and freshwaters. Applied and Environmental Microbiology, 2008, 74(17): 5317–5324. DOI:10.1128/AEM.02480-07
[30] Baker A C, Goddard V J, Davy J, Schroeder D C, Adams D G, Wilson W H. Identification of a diagnostic marker to detect freshwater Cyanophages of filamentous Cyanobacteria. Applied and Environmental Microbiology, 2006, 72(9): 5713–5719. DOI:10.1128/AEM.00270-06
[31] Larsen J B, Larsen A, Bratbak G, Sandaa R A. Phylogenetic analysis of members of the Phycodnaviridae virus family, using amplified fragments of the major capsid protein gene. Applied and Environmental Microbiology, 2008, 74(10): 3048–3057. DOI:10.1128/AEM.02548-07
[32] Rowe J M, Fabre M F, Gobena D, Wilson W H, Wilhelm S W. Application of the major capsid protein as a marker of the phylogenetic diversity of Emiliania huxleyi viruses. FEMS Microbiology Ecology, 2011, 76(2): 373–380. DOI:10.1111/fem.2011.76.issue-2
[33] Hopkins M, Kailasan S, Cohen A, Roux S, Tucker K P, Shevenell A, Agbandje-McKenna M, Breitbart M. Diversity of environmental single-stranded DNA phages revealed by PCR amplification of the partial major capsid protein. The ISME Journal, 2014, 8(10): 2093–2103. DOI:10.1038/ismej.2014.43
[34] Goldsmith D B, Crosti G, Dwivedi B, McDaniel L D, Varsani A, Suttle C A, Weinbauer M G, Sandaa R A, Breitbart M. Development of phoH as a novel signature gene for assessing marine phage diversity. Applied and Environmental Microbiology, 2011, 77(21): 7730–7739. DOI:10.1128/AEM.05531-11
[35] Sullivan M B, Huang K H, Ignacio-Espinoza J C, Berlin A M, Kelly L, Weigele P R, DeFrancesco A S, Kern S E, Thompson L R, Young S, Yandava C, Fu R, Krastins B, Chase M, Sarracino D, Osburne M S, Henn M R, Chisholm S W. Genomic analysis of oceanic cyanobacterial myoviruses compared with T4-like myoviruses from diverse hosts and environments. Environmental Microbiology, 2010, 12(11): 3035–3056. DOI:10.1111/emi.2010.12.issue-11
[36] Breitbart M, Miyake J H, Rohwer F. Global distribution of nearly identical phage-encoded DNA sequences. FEMS Microbiology Letters, 2004, 236(2): 249–256. DOI:10.1111/fml.2004.236.issue-2
[37] Lavigne R, Seto D, Mahadevan P, Ackermann H W, Kropinski A M. Unifying classical and molecular taxonomic classification: analysis of the Podoviridae using BLASTP-based tools. Research in Microbiology, 2008, 159(5): 406–514. DOI:10.1016/j.resmic.2008.03.005
[38] Marston M F, Taylor S, Sme N, Parsons R J, Noyes T J E, Martiny J B H. Marine cyanophages exhibit local and regional biogeography. Environmental Microbiology, 2013, 15(5): 1452–1463. DOI:10.1111/emi.2013.15.issue-5
[39] Chen F, Suttle C A. Evolutionary relationships among large double-stranded DNA viruses that infect microalgae and other organisms as inferred from DNA polymerase genes. Virology, 1996, 219(1): 170–178. DOI:10.1006/viro.1996.0234
[40] Culley A I, Lang A S, Suttle C A. High diversity of unknown picorna-like viruses in the sea. Nature, 2003, 424(6952): 1054–1057. DOI:10.1038/nature01886
[41] Culley A I, Lang A S, Suttle C A. Metagenomic analysis of coastal RNA virus communities. Science, 2006, 312(5781): 1795–1798. DOI:10.1126/science.1127404
[42] Culley A I, Steward G F. New genera of RNA viruses in subtropical seawater, inferred from polymerase gene sequences. Applied and Environmental Microbiology, 2007, 73(18): 5937–5944. DOI:10.1128/AEM.01065-07
[43] Hambly E, Tétart F, Desplats C, Wilson W H, Krisch H M, Mann N H. A conserved genetic module that encodes the major virion components in both the coliphage T4 and the marine cyanophage S-PM2. Proceedings of the National Academy of Sciences of the United States of America, 2001, 98(20): 11411–11416. DOI:10.1073/pnas.191174498
[44] 贺纪正, 袁超磊, 沈菊培, 张丽梅. 土壤宏基因组学研究方法与进展. 土壤学报, 2012, 49(1): 155–164.
[45] Thurber R V. Current insights into phage biodiversity and biogeography. Current Opinion in Microbiology, 2009, 12(5): 582–587. DOI:10.1016/j.mib.2009.08.008
[46] Rohwer F, Edwards R. The phage proteomic tree: a genome-based taxonomy for phage. Journal of Bacteriology, 2002, 184(16): 4529–5535. DOI:10.1128/JB.184.16.4529-4535.2002
[47] Mokili J L, Rohwer F, Dutilh B E. Metagenomics and future perspectives in virus discovery. Current Opinion in Virology, 2012, 2(1): 63–77. DOI:10.1016/j.coviro.2011.12.004
[48] Kim K H, Chang H W, Nam Y D, Roh S W, Kim M S, Sung Y, Jeon C O, Oh H M, Bae J W. Amplification of uncultured single-stranded DNA viruses from rice paddy soil. Applied and Environmental Microbiology, 2008, 74(19): 5975–5985. DOI:10.1128/AEM.01275-08
[49] Williamson K E, Radosevich M, Wommack K E. Abundance and diversity of viruses in six Delaware soils. Applied and Environmental Microbiology, 2005, 71(6): 3119–3125. DOI:10.1128/AEM.71.6.3119-3125.2005
[50] Fierer N, Breitbart M, Nulton J, Salamon P, Lozupone C, Jones R, Robeson M, Edward R A, Felts B, Rayhawk S, Knight R, Rohwer F, Jackson R B. Metagenomic and small-subunit rRNA analyses reveal the genetic diversity of Bacteria, Archaea, Fungi, and Viruses in soil. Applied and Environmental Microbiology, 2007, 73(21): 7059–7066. DOI:10.1128/AEM.00358-07
[51] Zablocki O, Van Zyl L, Adriaenssens E M, Rubagotti E, Tuffin M, Cary S C, Cowan D. High-level diversity of tailed phages, eukaryote-associated viruses, and virophage-like elements in the metaviromes of Antarctic soils. Applied and Environmental Microbiology, 2014, 80(22): 6888–6897. DOI:10.1128/AEM.01525-14
[52] Reavy B, Swanson M M, Cock P J, Dawson L, Freitag T E, Singh B K, Torrance L, Mushegian A R, Taliansky M. Distinct circular single-stranded DNA viruses exist in different soil types. Applied and Environmental Microbiology, 2015, 81(12): 3934–3945. DOI:10.1128/AEM.03878-14