时间尺度细菌组成

1、对虾养殖过程中浮游细菌群落在时间尺度上可预测的演替规律裘钱玲琳日，熊金波日，朱建林b,赵群芬a，陈和平我，张德民a*a海洋学院，宁波大学，宁波，315211，中国b建筑工程与环境学院，宁波大学，宁波，315211，中国摘要：微生物群落的空间分布，即群落相似性随空间距离衰减已被广泛接受，但鲜有微生物群落在时间尺度上演替规律的报道。因此，微生物群落在空间和时间尺度上是否具有相似的变化规律仍不明确。本实验利用高通量测序，研究对虾养 殖过程中细菌在时间尺度上的变化规律。 研究发现：细菌群落的相似性随时间显 著（p=0.002）降低，周转率w=0.004，而多样性并没有显著变化。群落结构随时 间不断变化

2、，但连续采样点样品间相似性却较高。 不同池塘细菌群落的时间演替 有着相似的规律，说明细菌群落的时间尺度上动态演替一定程度上是可以预测的。 细菌群落结构的变化与叶绿素 a、总氮、磷酸盐、碳氮比紧密相关。本研究建立 了细菌群落相似性在时间尺度上变化的模型， 表明细菌群落在时间和空间尺度上 具有相似的演替规律。关键词：浮游细菌；时间动态演替；周转率；相似性一时间关系1引言“养殖病害是养殖主体、养殖环境以及周围菌群相互作用的结果”，三者密切相关。对虾养殖过程中微生物群落的动态演替（尤其在健康塘以及病害塘特 定微生物群落）已被证实2-5。这表明浮游细菌在维持虾体健康方面起着重要作 用6, 7，特定的细菌

3、可以区分对虾的健康状态。因此细菌群落时间演替的研究能有力地预测对虾疾病的出现。为了维持高产量，饲料的供应量远远超过对虾的需求量，这使得受纳水体逐 渐富营养化8, 9。然而，资源的分区是微生物的组成、功能结构的关键因素10,11， 可能会导致微生物群落在对虾养殖过程中的演替。近来研究表明，对虾养殖中微生物群落体现出在时间尺度上的演替6, 7, 12。这些研究只停留在利用特定的技术 研究13特定类群（如弧菌）7, 14，并不能全面理解群落在时间尺度上如何加剧富 营养化。尽管已经研究出天然水生生态系统中群落在时间尺度上的动态演替模式 2,5，目前不清楚动态演替在人为操纵的池塘中（一个不断变化的环境）

4、是否具 有可预见性。在空间尺度上微生物群落相似性随距离衰减15, 16（即群落差异性随地理位置 增强），然而相似性与时间的关系（STRs）模式鲜有深入的研究。Jones以及 他的团队17假设水体细菌群落组成在时间尺度以及空间尺度上的变化是同等重 要的。然而目前对于时间尺度上的研究一般针对于系统发育与时间的关系，即在时间尺度上的不同生态系统物种累积丰度19-22。因此，微生物群落在时间尺度上 是否存在类似于空间尺度的时间与相似性之间的关系并不确定。由此，我们用类似空间尺度的数学公式：S = cTw 18，其中标度指数w反映了微生物群落的周转 率。虾塘在放养虾苗之前经过消毒，是一个相对比较清洁的环

5、境，且各虾塘在养 殖全程采用相同的管理，实验在此基础上对微生物群落在时间尺度上的动态演替 进行了详细的研究23。不同的细菌类群在代谢作用上存在本质的不同，由此推 断由于对虾养殖过程中存在资源的分区17, 18，则间接表明细菌群落动态演替较 活跃（存在高周转率）。实验从时间尺度上采集对虾养殖池塘水样，并且应用高 通量测序技术探究浮游细菌群落在时间尺度上的动态演替规律以及找出主要环 境因子，并进一步探究在空间尺度的影响下相似性与时间之间是否存在一定关系。2材料与方法2.1实验设计和水样采集对虾养殖池塘位于宁波瞻岐（北纬 29 32东经121 3）。30个养殖塘大 小一致（约2000 m2），放养密

6、度、饲料类型以及养殖进程等都采用相同的管理。 2012年3月25日，统一放虾苗。从6个不同的养殖塘采集水样，采样距水表面 50cm。水样取自不同时间点（周期42d，从4月29日至6月10日，每次采样时 间点间隔610d），对应于放虾苗后的第 35d、45d、55d、63d、69d、77d，共收 集36个水样（6个养殖塘冷个时间点）。样品利用冰袋降温保存，并在 4h内送 往实验室。- 2 -水温和 pH 利用仪器在距水面 50cm 处现场测量。总有机碳含量( TOC)、 总氮(TN)、总磷(TP)、化学需氧量(COD)、NO3-和PO43-浓度的监测分 析按照标准方法(APHA ,1989) 2

7、4。叶绿素a的测定时，水样经25mm的GF/F 的过滤器过滤，用90%的二甲基甲酰胺在48C恒温下萃取24h,最后使用分光光 度计(UV-1601, Shimadzu)测定水体中的叶绿素 a浓度。2.2 DNA 提取采集大约3L水样提取DNA，水样应预先经尼龙筛(滤径100 ym)过滤，随 后经 0.2 卩山的滤膜(Millipore, Bost on, MA, USA )过滤。使用 Power Soil? DNA Isolation Kit 试剂盒(MO BIO Laboratories, Carlsbad, CA, USA)，按照试剂盒说 明书提取总 DNA。DNA 浓度通过 Na noD

8、rop ND-1000分光光度计(Nan oDrop Technologies, Wilmington, USA)测定，并贮存在-80C冰箱。2.3细菌16S rRNA基因扩增和454测序每个DNA样品各取50ng作为PCR扩增的模板，使用引物对细菌的16SrRNA 基因 19V1V3 高变区进行扩增， 27F：AGAGTTTGATCMTGGCTCAG ，通过罗 氏454高磷酸测序适配器 A 找出11bp大小的barcode标识序列。519R: GWATTACCGCGGCKGCTG，每个样品扩增有特定的 barcode.在50卩体系中 扩增：94C变性30s，55C退火30s，72C延伸30s

9、,进行30个循环，最后一次 72C延伸10min。将每个样品PCR产物合并，采用纯化试剂盒(TaKaRa)纯化， PCR产物采用 Quant-ltTM PicoGreen kit (Invitrogen, Carlsbad, CA)定量分析。等摩尔量 PCR 产物单独放置，用罗氏 FLX454 焦磷酸测序机( Roche Diagnostics Corporation, Branford, CT, USA)，从有 barcode的 27F方向开始读数。2.4 测序数据处理测序结果利用 QIIME 的程序 25进行处理，嵌合体检查和降噪处理。用 PyroNoise筛选出长度为200 450 bp

10、以及小于6的最大均聚物的序列26，并将 带有相同barcode标识序列的序列归为同一样品25。利用uclus*27确定细菌的种 系型，并将比对大于 97%的序列归为同一 OTU (Operational Taxonomic Unit)。 通过Greengenes数据库28分类各种系型。为避免误差，随机选取4500条序列的 子序列进行进一步分析。2.5 统计分析利用单因素方差（SPSS 13.0）分析不同采样点位的a多样性指数。采用非 度量多维尺度（NMDS）的方法评估整体细菌群落结构的差异而产生的B多样性29。非参数多元方差分析（perMANOVA ）分析不同采样时间对于微生物群落变 化的影响

11、。利用PAST计算相似比例，并分析不同采样点主要类群3, 31。为了减 少多重共线性，利用R v.2.11.0中vegan软件包32，基于方差膨胀因子（VIF） 对采用CCA分析主要环境因子与细菌群落之间的关系。2.6时间周转率评估和相似性一时间关系相似性一时间关系模式（STRs）采用幕函数估算，S = cTw,其中w :尺度 指数，表示细菌群落周转率，这从空间尺度指数借鉴而来18。S?rensen指数更着 重于类群分析33，因此利用STRs来评估微生物群落可以发现相似性随着时间推 移累积降低。由于采样时间的设计，每个池塘的分析数据点不是独立，我们分别 将采样点D35 （第一个采样点）与其他采

12、样点（D45，D55，D63，D69，D77） 进行计算细菌群落相似性，再根据采样时间取各自相似性平均值。幕指数 w的 计算式直接利用线性回归的方法，绘制log -log函数：logS = logc + wlogT，用来 评估相似性与时间的关系18, 21。3结果3.1水质化学指标由于是大棚养殖，水体温度以及pH变化相对比较稳定，温度变化范围：27.330.9C，pH变化范围：6.187.05。NO3-和PO43这两个环境因子在养殖塘 变化相对比较明显，TOC在整个对虾养殖周期持续增加（r = 0.730, P .su-l PUS1EJ!UAvtfiM DJ3 EM$ D55D69D77数量以

13、及Shannon指数。此外多样性与采样时间以及环境因子并不存在显著性相 关（r 0.5）。CxictsissmcdDiviMonOOl huKtarnycgieiDtmpmifotaTeri0 Chfant/Jrii Cff(1%）3.3细菌群落结构分析尽管在养殖周期群落丰度以及多样性没有明显的改变，通过非度量多维尺度（NMDS）对样品OTUs排序分析，主要在第一象限分开，这可以清楚得出在我们跟踪的虾塘中浮游细菌类群连续性的变化规律。NMDS轴1与取样时间存在线性显著相关（Pearson相关系数，p = 0.961 P 0.001）（见图2b）。这个模 式可以用perMANOVA方法进一步说明

14、取样时间对于确定细菌群落的组成起着 重要的作用（F = 1.80, P = 0.001）。在连续采样阶段，群落结构并没有发生显著 的变化，例如群落结构在D45和D55以及D63和D69之间并没有显著差别，然而D45和D63以及D63和D77之间存在显著的变化（P 0.05），详见表1。结 果表明，在短期内细菌群落在时间尺度上存在动态演替。a 0.200.150 10rjhxccsu 丄二0.050.00-0.05fl. 10D33V映774-0 .15 -02NMDS axis 140-0.20 -OS50eotoSampling day (days)图2非度量多维尺度（NMDS ）分析注：图

15、a即微生物群落随时间演替；图b即群落结构与采样时间显著正相关 表1基于非参数多元方差（perMANOVA ）分析群落差异性D35-1D45-1D55-1D63-1D69-1D35D450.672D550.0340.069D630.0350.0330.011D690.0330.0340.0040.451D770.0360.0340.0040.0030.073注：加粗值为各采样日期细菌群落存在显著差异性（p 0.05）3.4环境因子与细菌群落的关系利用CCA方差分析筛选出影响微生物群落结构的关键变量，将微生物群落 的分类结构和水体性质以及叶绿素 a与光合作用相联系。最后筛选出 TN（p = 0.0

16、01）， PO43- （ p = 0.030），碳氮比（p = 0.001），叶绿素 a（ p = 0.007） 4 个重 要指标，由CCA图知，群落变化和相关变量显著相关（见图 3）。群落结构在后 期（D63、D77）与前期（D35、D45）是不同的，主要在第一轴分开，与碳氮比 和叶绿素a成正相关，与TN和PO43-成负相关（见图3）。图3各样品与水质参数之间的CCA分析3.5相似性与时间关系为了评估浮游细菌群落的周转率，实验利用相似性一时间的关系（即STRs）， 来评估在人工养殖对虾池塘整个细菌群落的演替。通过细菌群落相似性和采样时 间用log函数转化，并通过线性回归来求出 STR斜率。S

17、TR显著（p = 0.002）， 其中指数w = 0.400 （见图4）。由此可见，细菌群落的变化趋势是可以预测的， 即在监测的各个塘中细菌变化的运动轨迹是相似的（见图5）。此外，累积OTUs和天数存在显著相关（p 1 一 .m_F=s4o.r =0 995./)= 0 0020.00.2040.60.810 1J 1.41.6Days between observations (1og|0 scale)图4细菌群落相似性一时间关系（STRs）a0.20.00 403*02 山 0Q 0 10.2 D30.4PCoA 1 (17J%-0 1-02-03 -Pond 5Pnnd hEW Pond

18、 3-0 24) 10.00.10.20 304PCoA 1图5不同养殖塘水样PCoA分析4 讨论实验最初是为探究细菌群落与对虾发病之间的关系， 因此只收集了对虾发病 高发期（放苗后58 8d）的水样34。为研究微生物组成时间演替，实验的采样 频率相对不严格， 然而，实验已经表明采样频率对于整体的周转率影响很小 （见 图 2b） 20。高密度对虾养殖会逐渐导致富营养化 9，然而浮游细菌群落对于环境微妙的 变化极为敏感 3 。细菌群落随时间存在明显的变化，细菌群落相似性随时间逐渐 降低（见图 4），然而其多样性维持不变。 细菌群落的时间变化与叶绿素 a、TN、 PO43- 和碳氮比（见图 3）是

19、显著相关的，这与特定种群形成的生态龛的原理相 一致35。叶绿素a 一直作为估算浮游植物数量和生物量的代表值36，在很大程度上影响了细菌群落结构的时间演替 37。相近的采样点的微生物群落的系统发育更为接近 2, 20。同样，由细菌群落的 变化反映出菌株水平上也存在一定影响，在连续性采样和其他采样点相比，如 D45和D55比较、D63和D69比较，均比其他采样点的细菌群落更加相近，见 图 2。此外，在采样期间平行演替的模式除外，同一天不同池塘之间发现有更多 相类似的细菌群落（见图 5）。这清楚表明养殖塘中细菌群落结构随时间的变化 不是随机的， 而是高度可预测的。 之后的实验需要验证对虾养殖周期循环

20、的这种 方式。研究表明浮游细菌群落总细胞丰度的变化远小于随时间推移群落组成的变化38, 39。随时间推移群落动态演替以及其功能的研究存在两种机制，即“调整 方案”和“替代方案” 38。本研究为后一种机制，资源环境的连续变化导致不 同类群的浮游细菌组合物形成。 很多研究表明群落结构决定应对环境变化的代谢 响应 40，主要由于细菌群落结构的变化速率和其功能紧密联系，以及微生物组 成影响生态变化的进程 11, 41。因此，了解细菌的动态演替对于更好的预测生态 系统的功能是至关重要的。虾池细菌群落的相似性与时间关系（STRs）指数（w = 0.400）比在多个环境中 细菌连续变化的斜率（ 0.020.

21、03）更大4。这种差异的形成可用以下原因解释， 我们利用焦磷酸测序技术， 检测到的微生物更加丰富， 特别是能检测出一些罕见 的物种。尽管我们的测序深度并没有达到饱和，但Shade等测序深度比我们低- 10 -的多，每个样品测序深度为 4863526ba相关研究表明，样品群落结构的变化 对于测序深度不是特别敏感，例如，每个样品有降噪序列约占90%，进而分析B多样性42。另外，在一定程度上斜率可能由于随机过程而变大43, 44，例如对虾养殖过程中每天进排水的影响。 尽管这样， 类群与时间关系的斜率和在其他环境 中细菌类群累积率相似 20-22。此外，时间尺度上细菌群落动态演替和环境因子显 著相关（

22、图 3），可以预测的演替轨迹表明影响因素是可以确定的而不是随机的， 其在细菌时间演替方面起决定性作用。按照空间尺度模式，实验总结出相似性与时间关系的曲线 （STRs）来评估细 菌群落周转率。 最近研究表明， 随着时间推移环境因子的变化与水生生物空间异 质性类似 17，即推动群落动态变化与生态过程相类似。其实长期以来，都是将 空间代替时间研究微生态演替时序 45, 46。然而在假设稳定生物化学因素等方面， 这种方法也是饱受争议的46。生态研究需要转变连续物种研究的模式（一般 a 多样性不随时间变化），而去研究B多样性，将其作为生态过程研究的基础。总之，这项研究的重要性在于我们提供了说明细菌群落相

23、似性（B多样性）而不是累积群落丰度的时间动态演替直接的证据16, 18。这反映了微生物随距离衰减的相似性的模式。尤其在较短的时间内（42d）,演替轨迹是可以预见的。 浮游细菌群落时间演替是可以高度预测的 45, 46，这与地理化学特点是显著相关 的 46 。这对于人为调节维持虾塘生态功能是一个机会也是一个挑战。综上，我 们研究出细菌群落结构在空间尺度可以预测的模式， 与细菌时空变化的幂次模型 类似。- 14 -参考文献1 Snieszko SF (1974) The effect of environmental stress on outbreaks of infectious diseas

24、es of ?sh. J Fish Biol 6:197 208.-2 Gilbert JA, Field D, Swift P, Newbold L, Oliver A, Smyth T, et al (2009) The seasonalsuccession of microbial communities in the Western English Channel using 16S rDNA-tag pyroseque ncing. Environ Microbiol 11:3132-3139.3 Or A, Shtrasler L, Gophna U (2012) Fine-sca

25、le temporal dynamics of a fragmented lotic microbial ecosystem. Sci Rep 2:207.4 Shade A, Caporaso JG, Handelsman J, Knight R, Fierer N (2013) A meta-analysis of changes in bacterial and archaeal commu ni ties with time. ISME J 7:1493 -1506.5 Wu QL, Hahn MW (2006) High predictability of the seasonal

26、dynamics of a species-likePolynucleobacter population in a freshwater lake. Environ Microbiol 8:1660-1666.6 Lemo nnier H, Courties C, Mug nier C, Torr ton J, Herbla nd A (2010) Nutrie nt and microbial dynamics in eutrophying shrimp ponds affected or unaffected by vibriosis. Mar Pollut Bull 60:402-41

27、1.7 Lucas R, Courties C, Herbland A, Goulletquer P, Marteau AL, Lemonnier H (2010) Etrophication in a tropical pond: Understanding the bacterioplankton and phytoplankton dynamics during a vibriosis outbreak using ?ow cytometric analyses. Aquaculture 310:112-121.8 Ferreira NC, Bonetti C, Seiffert WQ

28、(2011) Hydrological and water quality indices as management tools in marine shrimp culture. Aquaculture 318: 425 -433.9 Ma Z, Song X, Wan R, Gao L (2013) A modi?ed water quality index for intensive shrimp ponds of Litopenaeus vannamei. Ecol Indic 24:287 -293.10 Goldfarb KC, Karaoz U, Hanson CA, Sant

29、ee CA, Bradford MA, Treseder KK, WallensteinMD, Brodie EL (2011) Differential growth responses of soil bacterial taxa to carbon substrates of varying chemical recalcitrance. Front Microbiol 2:1-10.11 Xiong J, He Z, Van Nostrand JD, Luo G, Tu S, Zhou J, Wang G (2012a) Assessing the microbial communit

30、y and functional genes in a vertical soil profile with long-term arsenic contamination. PLoS One 7:e50507.12 Xu H, Zhu M, Jiang Y , Al-Rasheid KS (2010) Temporal species distributions of planktonicprotist communities in semi-enclosed mariculture waters and responses to environmental stress. Acta Oce

31、anol Sin 29:74-83.13 Case M, Leca EE, Leitaoc SN, SantAnna EE, Schwamborn R, De Moraes Junior AT (2008) Plankton community as an indicator of water quality in tropical shrimp culture ponds. Mar Pollut Bull 56:1343 -1352.14 Sung H, Hsu S, Chen C, Ting Y , Chao W (2001) Relationships between disease o

32、utbreak incultured tiger shrimp (Penaeus monodon) and the composition of Vibrio communities in pond water and shrimp hepatopancreas during cultivation. Aquaculture 192:101-110.15 Griffiths RI, Thomson BC, James P, Bell T, Bailey M, Whiteley AS (2011) The bacterialbiogeography of British soils. Envir

33、on Microbiol 13:1642-1654.16 Xiong J, Liu Y , Lin X, Zhang H, Zeng J, Hou J, et al (2012b) Geographic distance and pH drive bacterial distribution in alkaline lake sediments across Tibetan Plateau. Environ Microbiol 14:2457 -2466.17 Jones SE, Cadkin TA, Newton RJ, McMahon KD (2012) Spatial and tempo

34、ral scales ofaquatic bacterial beta diversity. Front Microbiol 3:318.18 Green JL, Holmes AJ, Westoby M, Oliver I, Briscoe D, Danger?eld M, et al (2 004) Spatial scali ng of microbial eukaryote diversity. Nature 432:74750.19 Kim Y, Jeong J, Wells GF, Park H (2013) General and rare bacterial taxa demo

35、nstrating different temporal dynamic patterns in an activated sludge bioreactor. Appl Microbiol Biotech 97:1755 -765.20 Redford A, Fierer N (2009) Bacterial succession on the leaf surface: a novel system for studying successional dynamics. Microb Ecol 58:189 -198.21 Van Der Gast CJ, Ager D, Lilley A

36、K (2008) Temporal scaling of bacterial taxa is in?u encedby both stochastic and deterministic ecological factors. Environ Microbiol 10:1411-1418.22 Wells GF, Park HD, Eggleston B, Francis CA, Criddle CS (2011) Fine-scale bacterial community dynamics and the taxa-time relationship within a full-scale

37、 activated sludge bioreactor. Water Res 45:5476-5488.23 G mez-Consarnau L, Lindh MV , Gasol JM, Pinhassi J (2012) Structuring of bacterioplankton communities by specific dissolved organic carbon compounds. Environ Microbiol 14:2361 -2378.24 APHA, AWW A, WEF. (1990) Standard methods for the examinati

38、on of water and wastewater, l7th ed. Washington. DC.25 Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al (2010) QIIME allows integration and analysis of high-throughput community sequencing data. Nat Meth 7:335-336.26 DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Bro

39、die EL, Keller K, et al (2006a)Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol 72:5069-5072.27 Edgar RC (2010) Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26:2460 -2461.28 DeSantis TZ, Hugenholtz P, Keller

40、 K, Brodie EL, Larsen N, Piceno YM, et al (2006b) NAST:a multiple sequence alignment server for comparative analysis of 16S rRNA genes. Nucleic Acids Res 34:394-399.29 Legendre P, Legendre L (1998) Numerical Ecology. Second English edition. Elsevier Science, Amsterdam, The Netherlands.30 Clarke KR (

41、1993) Non-parametric multivariate analysis of changes in community structure. Aust J Ecol 18:117-143.31 Hammer ?, Harper DA T, Ryan PD (2001) PAST: Paleontological statistics software package for education and data analysis. Palaeontologia Electronica 4:1.32 R Development Core Team (2011) R: A langu

42、age and environment for statistical computing. Vienna, Austria: R 21 Foundation for statistical computing.33 S?rensen TA (1948) A method of establishing groups of equal amplitude in plant sociologybased on similarity of species and its application to analyses of the vegetation of the Danish commons.

43、 Biologiske Skrifter 5:1-34.34 Lemonnier H, Herbland A, Salery L, Soulard B (2006) Summer syndromeLit”opeinaeusstylirostris grow out ponds in New Caledonia: zootechnical and environmental factors. Aquaculture 261:1039 -1047.35 Teeling H, Fuchs BM, Becher D, Klockow C, Gardebrecht A, Bennke CM, et al

44、 (2012) Substrate-controlled succession of marine bacterioplankton populations induced by a phytoplankton bloom. Scienee 336:608 -611.36 Domingues RB, Barbosa A, Galv? o H (2008) Constraints on the use of phytoplankton as a biological quality element within the water framework directive in Portugues

45、e waters. Mar Pollut Bull 56:1389 -395.37 Paver SF, Hayek KR, Gano KA, Fagen JR, Brown CT, et al (2013) Interactions between specific phytoplankton and bacteria affect lake bacterial community succession. Environ Microbiol 15:2489 -2504.38 Comte J, Del Giorgio PA (2011) Composition influences the pathway but not the outcome of the metabolic response of bacterioplankton to resource shifts. PLoS One 6:e25266.39 Cotner JB, Biddanda BA (2002) Small players, large role: Microbial influence on biogeochemical pr