×

污水处理厂微生物动力学的稀疏贝叶斯层次向量自回归模型。 (英语) Zbl 07708466号

摘要:污水处理厂(WWTP)的正常运行依赖于在众多相互竞争的微生物之间保持微妙的平衡。对其中复杂的相互作用网络有一个详细的了解,这不仅对最大限度地提高当前的操作效率至关重要,而且对新处理技术的有效设计也至关重要。宏基因组学通过对现有微生物DNA序列的分析,提供了对这些动态系统的深入了解。通过序列聚类推断出独特的分类群,形成操作分类单元(OTU),并从相应的序列计数中获得每个分类群的丰度估计。本研究中的数据包括来自污水处理厂活性污泥(AS)池的每周OTU计数以及相应的化学和环境(CE)协变量测量值。由于观测数据的高维性和稀疏性,直接将模型拟合到OTU数据是非常具有挑战性的。因此,第一步是使用基于季节性阶段的聚类方法将OTU聚合成12个微生物群落或“箱”。根据CE协变量的多元线性回归,假设12个箱子中的平均丰度随时间变化。然后使用一阶向量自回归(VAR)模型对平均值的偏差进行建模,该模型是对常用的广义Lotka-Volterra(gLV)模型的线性近似。通过在层次贝叶斯框架中进行推理,假设微生物群落之间的相互作用具有稀疏性,该框架使用VAR模型的自回归系数矩阵的收缩先验。在选择常规马蹄铁用于生物应用之前,通过分析模拟数据集来探索不同的收缩先验。发现氨和化学需氧量与几个箱子呈正相关,pH值与一个箱子呈正关系。这些结果得到了生物学文献中的发现的支持。还发现了一些负面的相互作用。这些新的生物学发现表明,不同箱子中的OTU可能在争夺资源,并且这些关系很复杂。尽管VAR模型比gLV模型简单,但它仍然能够为污水处理厂的微生物动力学提供有价值的见解。

MSC公司:

62-08 统计问题的计算方法
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] 阿赫勒格比,D.F。;比利奥,M。;Casarin,R.,结构向量自回归过程的贝叶斯图形模型,J.Appl。经济。,31, 357-386 (2016)
[2] 阿拉维,M。;关闭,S。;卡亚,M。;Spieck,E.,《温度影响活性污泥中亚硝酸盐氧化细菌的种群结构》,《环境》。微生物。代表,118-190(2009)
[3] 贝当古,M。;Girolama,M.,层次模型的哈密尔顿蒙特卡罗,(Upadhyay,S.K.;Singh,U.;Dey,D.K.;Loganathan,A.,《贝叶斯方法及其应用的当前趋势》(2015),CRC出版社:佛罗里达州博卡拉顿CRC出版社),79-101,第4章
[4] 布奇,V。;Tzen,B。;李,N。;西蒙斯,M。;塔努埃,T。;Bogart,E.,MDSINE:微生物群时间序列分析的微生物动力系统推理引擎,《基因组生物学》。,17, 121 (2016)
[5] 邦吉,J。;A.威利斯。;Walsh,F.,《微生物多样性研究中物种数量的估算》,年。修订状态申请。,1, 427-445 (2014)
[6] 卡彭特,B。;Gelman,A。;医学博士霍夫曼。;Lee,D。;古德里奇,B。;贝当古,M。;布鲁贝克,医学硕士。;Guo,J。;李,P。;Riddell,A.,Stan:概率编程语言,J.Stat.Softw。,76, 1-32 (2017)
[7] 卡瓦略,C.M。;波尔森,N.G。;Scott,J.G.,《通过马蹄铁处理稀疏性》,J.Mach。学习。第573-80号决议(2009年)
[8] 卡瓦略,C.M。;新几内亚波尔森。;Scott,J.G.,《稀疏信号的马蹄形估计器》,Biometrika,97,465-480(2010)·Zbl 1406.62021号
[9] 柯蒂斯,T.P。;头颅,I.M。;Graham,D.W.,《同行评论:工程生物学的理论生态学》,环境。科学。技术。,37, 64-70 (2003)
[10] Cydzik-Kwiatkowska,A。;杰琳斯卡,M.,《全面废水处理系统中的细菌群落》,《世界微生物杂志》。生物技术。,32, 66 (2016)
[11] Daims,H。;Lebedeva,E.V。;Pjevac,P。;Han,P。;Herbold,C。;Albertsen,M.,《硝化细菌的完全硝化作用》,《自然》,528,504-509(2015)
[12] 大坝,P。;丰塞卡。;Konstantinidis,K.T。;Voit,E.O.,门多塔湖复杂微生物集合种群的动态模型,npj系统。生物应用。,第2条,第16007页(2016年)
[13] 大卫,洛杉矶。;马特拉纳,A.C。;弗里德曼,J。;麻省理工学院Campos-Baptita。;M.C.布莱克本。;佩罗塔,A。;埃尔德曼,S.E。;Alm,E.J.,《宿主生活方式在日常时间尺度上影响人类微生物群》,《基因组生物学》。,15,R89(2014)
[14] Dugan,P.R.,Zoogloea属,(Starr,M.P.;Stolp,H.;Trüper,H.G.;Balows,A.;Schlegel,H.G,《原核生物:细菌的栖息地、分离和鉴定手册》(1981),斯普林格·柏林-海德堡:斯普林格尔·柏林-海德堡-柏林,海德堡),764-770
[15] 艾勒,A。;Heinrich,F。;Bertilsson,S.,湖泊浮游细菌分类群之间的相干动力学和关联网络,ISME J.,6330-342(2012)
[16] 浮士德,K。;拉提,L。;Gonze,D。;德沃斯,W.M。;Raes,J.,《宏基因组学与时间序列分析:揭示微生物群落动力学》,Curr。操作。微生物。,25, 56-66 (2015)
[17] Fisher,C.K。;Mehta,P.,使用稀疏线性回归从宏基因组时间序列中识别人类肠道微生物组中的关键物种,PLoS ONE,9,1-10(2014)
[18] Gefang,D.,用于VAR收缩的贝叶斯双自适应弹性网拉索,国际预测杂志。,30, 1-11 (2016)
[19] Gelfand,A.E。;Fuentes,M。;Guttorp,P。;Diggle,P.,《空间统计手册》,查普曼和霍尔/CRC现代统计方法手册(2010年),泰勒和弗朗西斯·Zbl 1188.62284号
[20] E.I.乔治。;McCulloch,R.E.,《通过吉布斯抽样选择变量》,美国统计协会,88,881-889(1993)
[21] 吉本斯,S.M。;科尼,S.M。;Smillie,C.S。;Alm,E.J.,《人类肠道微生物组的两种动态机制》,《公共科学图书馆·计算》。生物学,13,1-20(2017)
[22] Girolma,M。;Calderhead,B.,Riemann流形Langevin和Hamilton Monte Carlo方法,J.R.Stat.Soc.B,73,123-214(2011)·Zbl 1411.62071号
[23] 戈亚尔,A。;杜宾基纳,V。;Maslov,S.,稳定婚姻问题解释的微生物群落中的多个稳定状态,ISME J.,12,2823-2334(2018)
[24] Granger,C.W.J.,《用计量经济学模型和交叉谱方法研究因果关系》,《计量经济学》,第37期,第424-438页(1969年)·Zbl 1366.91115号
[25] Imhoff,J.F.,《红杆菌》(Bergey’s Manual of Systematics of Archaea and细菌系统学手册(2015),美国癌症学会),1-4
[26] 考尔,A。;曼达尔,S。;Davidov,O.,《存在多余零的微生物组数据分析》,Front。微生物。,8, 2114 (2017)
[27] van Kessel,M.A.H.J。;斯佩斯,D.R。;阿尔伯森,M。;尼尔森,P.H。;Op den Camp,H.J.M。;Kartal,B。;Jetten,M.S.M。;Lücker,S.,《单个微生物的完全硝化》,《自然》,528555-559(2015)
[28] 科赫,H。;吕克尔,S。;阿尔伯森,M。;基特辛格,K。;Herbold,C。;Spieck,E.,《来自Nitrospira属普遍存在的亚硝酸盐氧化细菌的扩展代谢多功能性》,Proc。国家。阿卡德。科学。,112, 11371-11376 (2015)
[29] 科诺普卡,A。;林德曼,S。;Fredrickson,J.,《微生物群落的动力学:揭示识别原理的机制》,ISME J.,9,1488-1495(2015)
[30] Lee,K.H。;库尔,B。;Moscicki,A.B。;巴斯特,B。;Starr,J.,《多元零膨胀模型的贝叶斯变量选择:微生物计数数据的应用》,Biostat。(2018)
[31] Lotka,A.J.,《物理生物学要素》,科学。掠夺。特温特。世纪(1919-1933)、21、341-343(1926)
[32] McIlroy,S.J。;桑德斯,A.M。;阿尔伯森,M。;Nierychlo,M。;麦克罗伊,B。;Hansen,A.A.,MiDAS:活性污泥微生物现场指南,数据库,2015(2015)
[33] 梅赫拉尼,M。;索博特卡,D。;Kowal,P。;西塞尔斯基,S。;Makinia,J.,氮螺旋菌在脱氮系统中的发生和作用,生物资源。技术。,303,第122936条pp.(2020)
[34] 米切尔·T·J。;Beauchamp,J.J.,线性回归中的贝叶斯变量选择,美国统计协会,832,1023-1032(1988)·Zbl 0673.62051号
[35] 穆尼尔,J。;莫内,C。;Vallaeys,T。;Arditi,R。;萨古,A。;Hélias,A。;Irlinger,F.,奶酪微生物群落中的微生物相互作用,应用。环境。微生物。,74, 172-181 (2008)
[36] Neal,R.M.,使用哈密顿动力学的MCMC,(Brooks,S.;Gelman,A.;Jones,G.;Meng,X.L.,马尔可夫链蒙特卡罗手册。马尔可夫链蒙特卡罗手册,Chapman&Hall/CRC现代统计方法手册(2011)),113-162·兹比尔1229.65018
[37] O.Ovaskainen。;蒂霍诺夫,G。;邓森,D。;格里坦,V。;Engen,S。;Sther,B。;Abrego,N.,物种丰富的社区中物种之间的相互作用是如何构建的?一种分析时间序列数据的新方法皇家社会学学士。科学。,284,第20170768条pp.(2017)
[38] Piironen,J。;Vehtari,A.,马蹄铁和其他收缩前的稀疏信息和规则化,Electron。J.Stat.,11,5018-5051(2017)·Zbl 1459.62141号
[39] Pujalte,M.J。;Lucena,T。;Ruvira,医学硕士。;Arahal,D.R。;Macián,M.C.,《红杆菌科》(The family rhodobacteraceae)(原核生物:α蛋白杆菌和β蛋白杆菌(The Pronuclees:Alphaproteobacteria and Betaproteobaceria)(2014),斯普林格-柏林-海德堡出版社),439-512
[40] 新墨西哥州Shchegolkova。;克拉斯诺夫,G.S。;Belova,A.A。;德米特里夫,A.A。;Kharitonov,S.L。;Klimina,K.M。;梅尔尼科娃,N.V。;Kudryavtseva,A.V.,不同废水成分处理厂活性污泥的微生物群落结构,Front。微生物。,7, 90 (2016)
[41] RStan:R接口到Stan(2020),R包版本2.19.31
[42] R.R.斯坦因。;布奇,V。;北卡罗来纳州图桑。;巴菲,C.G。;Rätsch,G。;Pamer,E.G。;桑德,C。;Xavier,J.B.,《时间序列推断的生态建模:肠道微生物群的动力学和稳定性洞察》,《公共科学图书馆·计算》。生物学,9,1-11(2013)
[43] Volterra,V.,从数学角度考虑物种丰度的波动,《自然》,118558-560(1926)
[44] 瓦格纳,M。;Rath,G。;Koops,H.P。;Flood,J。;Amann,R.,《污水处理厂硝化细菌的原位分析》,《水科学》。技术。,34, 237-244 (1996)
[45] 王,Q。;Garrity,G.M。;Tiedje,J.M。;Cole,J.R.,《用于将rRNA序列快速分配到新细菌分类法中的朴素贝叶斯分类器》,Appl。环境。微生物。,73, 5261-5267 (2007)
[46] Wetzel,R.G.,《氮循环》,(Wetzell,R.G,《湖沼学》(2001),学术出版社:圣地亚哥学术出版社),205-237
[47] 威廉姆斯,K.P。;Sobral,B.W。;Dickerman,A.W.,Alphaproteobacteria的健壮种树,J.细菌学。,4578-4586 (2007)
[48] 夏,Y。;Sun,J。;Chen,D.G.,《微生物组数据与R的统计分析》(2018),Springer Singapore:Springer新加坡·Zbl 1407.62024号
[49] 徐,L。;Xu,X。;D.孔。;顾,H。;Kenney,T.,随机广义Lotka-Volterra模型及其在学习微生物群落结构中的应用(2020年)
[50] 张杰。;Lin,W.,logistic正态多项式模型的可缩放估计和正则化,生物计量学,751098-1108(2019)·Zbl 1448.62197号
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。