多学科风险评估方法能够识别主要水道内的退化段落：以瓜达尔基维尔河为例

《Journal of Contaminant Hydrology》：A multidisciplinary approach to risk assessment enables the identification of deteriorated sections within a major waterway: The Guadalquivir River as a case study

　　

迄今为止，很少有研究评估过伊比利亚半岛南部主要水文流域——瓜达尔基维尔河的水质。在此，研究人员首次提出了一项综合性研究，整合了物理化学表征、细菌群落评估以及沿其整个水道的质量指数。水质指数（Water Quality Index, WQI）显示，下游水质逐渐且持续恶化，其中10%的站点显示为一般水质，出现了富营养化迹象（约400?mg/m³叶绿素）和微生物污染（大肠杆菌Escherichia coli，高达3000?CFU/100?mL）。中下游段落显示出有机碳、金属、悬浮固体和浊度的显著变化，以及受农业活动影响的高硝酸盐（高达19?mg?L^?1）和氨（高达1.4?mg?L^?1）水平。水质还可能受到支流影响的负面影响。由于微生物无处不在并发挥着不可或缺的环境作用，研究人员评估了水微生物组组成、结构和功能性的变化。宏条形码分析揭示了显著的变化，尤其是在上游河段，尽管该区域普遍存在与良好水质相关的参数。在该段落，α-多样性（α-diversity）显著下降，微gAMBI生态指数（microgAMBI ecological index）值非常高，以及潜在的功能变化揭示了环境退化的迹象。此外，预测到的“ABC转运蛋白”（ABC transporters）和“双组分系统”（two-component system）的高普遍性表明，由少数耐受性机会微生物（尤其是丛毛单胞菌科Comamonadaceae）驱动的对不利条件的响应。该区域以广泛的橄榄单一栽培为特征，传统上一直暴露于杀生物剂的污染。在实际意义上，该研究证实，尽管广泛认可的指数可能报告较高的整体水质，但补充性的微生物组分析可能有助于揭示影响整个生态系统的未检测到的危害。这项多学科研究可作为设计风险图以识别需采取行动的“热点”（hotspots）的模型，从而改善重要水资源的管理。

本研究以西班牙南部伊比利亚半岛主要水道——瓜达尔基维尔河为案例，探讨了多学科风险评估方法在识别水资源退化段落的潜力。研究背景指出，作为安达卢西亚地区经济命脉的瓜达尔基维尔河，长期受农业（尤其是橄榄单一栽培）、城市及工业活动的多重影响，尽管欧水框架指令要求实现良好生态状态，现有水质评估多依赖单一理化指标，缺乏对微生物群落结构及功能的整合分析，难以全面揭示隐蔽性生态风险。研究人员旨在通过物理化学表征、水质指数计算及细菌宏条形码技术的综合应用，沿河全程系统评估水质退化程度，并识别关键影响段落。

研究人员于2021年11月至12月在瓜达尔基维尔河沿岸20个站点（G1-G20）采集水样，开展了以下研究：物理化学变量分析（包括pH、溶解氧、总有机碳、硝酸盐、氨、叶绿素、金属等19项指标），计算基于美国国家卫生基金会的水质指数（WQI），并通过16S rRNA基因扩增子测序（V8高变区）对细菌群落进行组成、多样性及生态指数评估；同时利用Tax4Fun软件预测微生物功能通路，结合冗余分析（RDA）和偏最小二乘判别分析（PLS-DA）探究理化变量与微生物组变化的关联。主要结论包括：WQI显示下游水质逐渐恶化，中下游段落（G16-G18）为“一般”等级，富营养化（叶绿素高达~400?mg/m³）和微生物污染（大肠杆菌Escherichia coli高达3000?CFU/100?mL）显著；上游河段（G2-G7）虽理化指标良好，但α-多样性显著下降、微gAMBI生态指数达“差”至“坏”等级，功能预测显示ABC转运蛋白和双组分系统高度表达，指示耐受性机会菌（如丛毛单胞菌科Comamonadaceae）主导的应激响应，与橄榄单作长期暴露于杀生物剂相关。该研究的重要意义在于，证实了广泛使用的综合指数（如WQI）可能掩盖早期生态退化，而补充微生物组分析能为绘制风险地图、优先治理“热点”区域提供多学科支撑。论文发表在《Journal of Contaminant Hydrology》。

关键的技术方法包括：（1）沿瓜达尔基维尔河20个站点（G1-G20）的现场水样采集与分析，样本来自安达卢西亚海洋科学研究所（CSIC）和科尔多瓦大学；（2）物理化学表征测定pH、溶解氧、总有机碳、硝酸盐、氨、叶绿素、金属等19项变量；（3）基于美国国家卫生基金会的水质指数（WQI）计算，整合溶解氧、总大肠菌群、pH、生化需氧量等9个加权参数；（4）16S rDNA扩增子测序（V8区域）分析细菌群落组成、α-多样性（Shannon-Wiener指数、Pielou均匀度）及微gAMBI生态指数；（5）Tax4Fun功能预测代谢通路，结合K-means聚类和Pearson相关分析；（6）冗余分析（RDA）和偏最小二乘判别分析（PLS-DA）评估理化变量与微生物组的关联。

研究结果：
3.1 物理化学表征与水质评价：通过K-means聚类将变量分为5个簇，上游（G1-G9）以高溶解氧、碱度和氧化还原电位为主；中游（G11-G16）硝酸盐、叶绿素a/b及大肠杆菌水平显著升高，富营养化趋势明显，与农业和城市输入相关；G16-G18段落悬浮固体、浊度、铁和锰浓度激增，与侵蚀和污水排放有关；下游河口（G19-G20）电导率、硬度和主要离子骤升，指示海水入侵和除草剂（如敌草隆、特丁津）输入。WQI沿河逐渐下降，G16-G17达“一般”下限，但整体未超标，表明仅靠WQI可能低估生态风险。

3.2 细菌群落分类学分析：16S rRNA测序鉴定出140个细菌科，优势门为变形菌门（Proteobacteria，>70%）和放线菌门（Actinobacteria，>30%）。α-多样性（Shannon-Wiener指数）在上游G2-G7和河口G19-G20显著低于源头G1（p<0.001）。微gAMBI指数显示，上游G2-G7段生态状态从“良好”恶化至“差”或“坏”，中游G10-G15略有恢复，但河口再次恶化。相关性分析表明α-多样性与溶解氧和总有机碳呈负相关，微gAMBI正相关，提示物理化学变量驱动微生物组变化。

3.3 细菌组功能分析：Tax4Fun预测的50条最丰富代谢通路分成3个簇。簇A（如核糖体、氧化磷酸化、肽聚糖生物合成）在中下游丰度增加，反映代谢活性增强；簇B（如ABC转运蛋白、双组分系统、苯甲酸降解）在上游G2-G7高度富集，指示对胁迫的响应和芳香族化合物降解潜力；簇C（如鞭毛组装、细菌趋化性、谷胱甘肽代谢）在河口丰度上升，关联渗透胁迫和低氧环境。Pearson相关分析显示，微gAMBI与谷胱甘肽代谢正相关，与氧化磷酸化和核苷酸切除修复负相关，说明高环境胁迫优先驱动解毒和生存功能。

3.4 全局分析：冗余分析（RDA）以细菌科组成为响应变量，理化变量解释总变异的59.1%（RDA1 40.9%，RDA2 18.2%）。电导率、温度、悬浮固体和特丁津区分了河口段（G16-G20），假单胞菌科（Pseudomonadaceae）在该段占优；高溶解氧和总有机碳解释上游G2-G7段微gAMBI上升和α-多样性下降，丛毛单胞菌科（Comamonadaceae）丰度突出；总大肠菌群和硝酸盐区分中下游（G9-G15），纤维单胞菌科（Cellulomonadaceae）增加。偏最小二乘判别分析（PLS-DA）显示样本沿河道呈渐进式曲线演化，源头G1明显分离，中上游（G2-G15）与下游（G18-G20）形成两个主要簇，G16-G17为过渡段。

讨论总结：研究证实，尽管广泛认可的指数（如WQI）可能报告较高整体水质，但补充微生物组分析能揭示未检测到的危害。上游G2-G7段理化指标良好但生物指标差，可能源于长期橄榄单作累积的农药影响；中下游受农业和城市输入及支流影响明显；功能预测表明ABC转运蛋白和双组分系统富集是耐受菌应对胁迫的关键机制。该研究为多学科风险评估模型提供了范式，但需扩展至支流、水库及不同气象条件的采样以增强代表性。

结论翻译：这项对瓜达尔基维尔河的创新性研究证明了物理化学、微生物学和宏条形码系统表征在水生生态系统中的适用性。尽管需要覆盖更广泛气象条件的进一步研究，但水质指数（WQI）分析显示沿河水质逐渐恶化，可能源于农业和城市活动（铁、碳、氮、大肠杆菌）释放的污染物，或汇入主流的小支流。物理化学分析强调了管理流域扩散和点源污染的重要性，特别是在中下游段落。同时，不同分类水平细菌微生物的鉴定显示，沿河细菌α-多样性显著下降（通常与环境胁迫相关），并且无论是在上游河段（G2-G9）还是最后的内河段（G19-G20），生态状态（microgAMBI指数）均出现恶化。在河口，代谢途径上调（涉及氨基酸代谢）被解读为对渗透胁迫的适应。微生物组潜在功能性的变化揭示了环境退化迹象，特别是在G2–G7段落。因此，推测与有机芳香族化合物代谢相关途径及其他反映对不利环境条件响应的途径被诱导。在这一背景下，选择了一个非常有限的耐受性和机会性微生物群落。丛毛单胞菌科（Comamonadaceae）成为最丰富的科，包含具有高度降解和解毒环境污染物的成员。橄榄园是迄今为止最广泛的作物。高密度的橄榄种植（在大多数地点几乎是单一栽培），结合地形的陡峭度，使得该区域的瓜达尔基维尔河流域极易受到杀虫剂和除草剂的污染。该段落生物指数（microgAMBI、α-多样性）指示的更差生态状态可能反映了长期累积影响，而水质指数（WQI）未能完全检测到该影响。这一事实强调了结合物理化学和生物整合分析（如本研究所述）以全面评估水生生态系统生态状态的必要性。该信息对于主管部门就该水资源做出知情决策至关重要。然而，沿河水质评估应扩展到包括支流、水库和变化的气候条件，以及实施创新的基于传感器的系统以早期检测污染事件。