Dengue spatiotemporal patterns in Minas Gerais, Brazil, 2014-2023: regional epidemic forces dominate over the environmental impact of the Brumadinho dam collapse

本研究得出结论，尽管2014年至2023年间巴西米纳斯吉拉斯州的登革热发病率受区域流行病驱动表现出强烈的时空同步性，但2019年布鲁马迪纽大坝溃决对受影响帕拉奥佩巴河流域的登革热传播模式并无显著独立影响。

要点

以下是用简单语言和日常类比对该研究论文的解释。

大局观：洪水与波浪

想象巴西的米纳斯吉拉斯州是一座巨大而繁忙的城市。多年来，这座城市一直面临着一个反复出现的问题：登革热，一种由蚊子传播的疾病。把登革热爆发想象成巨大的海浪，在同一时间席卷整个城市，无论人们住在哪里，都会受到冲击。

2019 年 1 月，城市的一个特定区域发生了一场具体灾难：布鲁马迪纽大坝坍塌。这导致大量泥浆和废物涌入当地河流。科学家们想知道：这次局部的泥石流是否在河流沿岸的城镇中引发了新的、独立的登革热“海啸”，还是说登革热常规的“海浪”像淹没其他人一样淹没了这些城镇？

这项研究旨在通过分析 2014 年至 2023 年的数据来回答这个问题。

他们如何调查（侦探工作）

研究人员像侦探一样比较了两组城镇：

1. “泥浆”组：直接位于大坝坍塌路径上的城镇（帕拉奥佩巴河流域）。
2. “对照组”：位于同一卫生区域但未受泥浆冲击的邻近城镇。

他们逐年查看这些城镇的登革热病例数量。为了进行公平的比较，他们不仅仅统计原始数字（因为有些城镇比其他城镇大）；他们观察的是爆发的模式。这就像比较音乐的节奏，而不仅仅是音量。

他们的发现

结果出奇地清晰：大坝坍塌并没有创造出独特的登革热模式。

以下是他们发现的分解，使用简单的比喻：

• “区域节奏”是主宰：决定登革热何时爆发的最大因素是卫生区域（一个大的行政区域），而不是特定的河流。这就像整个社区可能因为全市范围内的事件而在同一时间听到警报声，而不是因为某户人家发生了火灾。大坝附近的城镇和远处的城镇都随着同样的节拍起舞。
• “巨浪”淹没了“局部涟漪”：研究发现，2016 年、2019 年和 2023 年，整个州都出现了登革热的巨大激增。2019 年的激增发生在大坝坍塌之时，但这是一场巨大的全州性事件，而不是由泥浆引起的局部事件。
  • 类比：想象你正试图听到隔壁传来的低语（大坝的影响），而旁边正在举办一场摇滚音乐会（全州范围的登革热疫情）。由于音乐会太吵，低语根本听不见。研究人员发现，只有这场“音乐会”（全州疫情）在驱动数据。
• 没有特殊的“泥浆区”：当他们根据登革热模式对城镇进行分组时，大坝附近的城镇并没有聚集在一起。它们与其他城镇混在一起。受灾地区没有独特的“登革热特征”。

为什么大坝没有改变情况？

你可能会想，“但是大坝破坏了环境！这难道不会滋生更多蚊子吗？”

研究人员提出了几个原因来解释为什么数据中没有显示出这一点：

1. “州级浪潮”太强大：2019 年的全国性疫情规模如此之大，很可能压倒了大坝可能造成的任何微小的局部变化。
2. 地形：河流附近的许多城镇是农村且分布分散。登革热通常在拥挤的城市中爆发，那里蚊子和人挤在一起。在这些农村地区，“人群”不够密集，大坝造成的环境破坏不足以引发大规模的局部爆发。
3. 数据局限性：该研究依赖于官方报告的患病人数。有时，当灾难发生时，医院会不堪重负，或者人们无法前往医院，因此病例未被报告。这可能掩盖了微小的局部影响，但基于他们确实拥有的数据，没有观察到任何影响。

结论

该研究得出结论，米纳斯吉拉斯州的登革热是由巨大的全州性力量（如天气模式和病毒周期）驱动的，而不是由像布鲁马迪纽大坝坍塌这样的局部环境灾难驱动的。

虽然大坝是一场造成巨大环境破坏的悲剧，但它并没有在受灾城镇登革热爆发的时间或模式上留下可检测的痕迹。冲击这些城镇的登革热“浪潮”与冲击该州其他地区的浪潮相同，其驱动因素远比河流谷地本身要大得多。

技术摘要

技术摘要：2014–2023 年米纳斯吉拉斯州登革热的区域驱动因素

问题陈述
登革热是巴西面临的一项重大公共卫生挑战，其中米纳斯吉拉斯州承受了显著的病例负担。登革热的传播动力学极为复杂，由气候、环境、人口和社会经济因素之间的相互作用所驱动。2019 年 1 月，布鲁马迪纽大坝发生灾难性溃坝，向帕拉奥佩巴河流域释放了约 1200 万立方米的铁矿石尾矿。该事件造成了严重的环境扰动，可能改变了埃及伊蚊的水生繁殖栖息地，并扰乱了供水系统，从而可能改变当地的登革热传播动力学。然而，在同期大规模全国性疫情爆发（特别是 2019 年和 2023 年）以及系统性监测局限性的背景下，将此类灾难引发的局部流行病学信号区分开来，仍是一个实证挑战。本研究旨在调查，与区域对照组相比，布鲁马迪纽大坝溃坝是否导致受影响市镇的登革热时空模式发生了可检测的、独立的变化。

方法学
作者利用巴西法定传染病信息系统（SINAN）的数据，对 2014 年至 2023 年米纳斯吉拉斯州的登革热报告病例进行了生态时空分析。

• 研究人群：分析涵盖了 67 个市镇：26 个被指定为“帕拉奥佩巴流域”（直接/间接受大坝溃坝影响）和 41 个“区域对照”（同一行政卫生区域内未受影响的市镇）。第三组州级对照用于提供描述性背景。
• 数据处理：计算了每千人的年度发病率，并使用市镇特定的 Z 分数进行标准化，以比较独立于绝对量级的时间模式。
• 分析方法：
  • 模式识别：使用基于皮尔逊相关性的层次聚类和非度量多维缩放（NMDS）评估时间相似性。
  • 多变量分析：采用置换多变量方差分析（PERMANOVA）和主成分分析（PCA）来量化变异来源（例如，卫生区域与流域状态），并识别驱动变异的年份。
  • 统计建模：拟合了一个以市镇为随机效应的线性混合效应模型，以检验 2019 年后发病率的变化。该模型将年份、帕拉奥佩巴流域状态、卫生区域及其交互作用作为固定效应。使用边际均值估计了前/后对比。

主要结果

• 时间同步性：登革热发病率在整个州表现出强烈的时间同步性，主要疫情高峰在 2015–2016 年、2019 年和 2023 年同时出现。无论暴露状态如何，这些高峰在所有市镇中均被观察到。
• 变异驱动因素：卫生区域解释了时间发病率曲线变异的 31.5%（$p = 0.001$），而帕拉奥佩巴流域状态未解释显著的变异（$p = 0.998$）。层次聚类识别出五个不同的市镇组，经过多重检验校正后，其中没有任何一组显著富集帕拉奥佩巴流域市镇。
• 多变量驱动因素：PCA 表明，年际变异主要由疫情年份驱动，特别是 2023 年、2019 年和 2016 年。
• 灾难的影响：在线性混合效应模型中，年份是一个显著预测因子（$p < 0.001$），但帕拉奥佩巴流域状态及其与时间的交互作用并不显著。
  • 在非暴露（区域对照）市镇中，2019 年后的发病率显著增加（$p < 0.001$）。
  • 在 2019 年后的帕拉奥佩巴流域市镇中，未检测到统计学显著的变化（$p = 0.088$）。
  • 年份与帕拉奥佩巴流域状态之间的交互作用不显著（$p = 0.197$）。

主要贡献
本研究提供了对重大环境灾难对虫媒传染病传播流行病学影响的稳健、比较评估。通过利用混合效应框架控制卫生区域层面的决定因素，并采用同期对照组，作者将大坝溃坝的特定信号从更广泛的州级疫情力量中分离出来。研究结果挑战了布鲁马迪纽灾难独立改变当地登革热动力学的假设，反而表明在此背景下，区域疫情过程比局部环境扰动占据主导地位。

意义与主张
作者得出结论，研究期间米纳斯吉拉斯州的登革热动力学主要由区域和全州范围的疫情过程驱动，没有统计学上显著的独立效应可归因于布鲁马迪纽大坝溃坝对报告登革热模式的影响。

• 区域驱动因素的主导地位：结果强化了卫生区域作为有意义的流行病学单元的观点，其中共享的决定因素（气候、媒介生态、人口连通性）比局部环境扰动更强烈地调节登革热风险。
• 大规模疫情的混杂：研究强调了当局部灾难效应与大规模、同步的全国性疫情爆发（例如 2019 年和 2023 年疫情）同时发生时，检测局部灾难效应的困难。灾难与重大全国性激增之间的时间重叠起到了关键的混杂作用。
• 公共卫生意义：研究结果强调了使用稳健的比较设计来评估环境灾难后果的必要性。这表明，虽然对媒介栖息地的局部影响可能存在，但在大规模疫情爆发的背景下，这些影响可能过于细微，无法在年度、市镇级别的分辨率中被检测到。作者指出，被动监测漏报（因新冠疫情而加剧）以及缺乏昆虫学数据等局限性可能会掩盖微妙的局部效应，因此需要未来进行具有更高时空分辨率的研究。