《JOURNAL OF FOOD QUALITY》：Secondary Unsupervised Profiling of Monte Carlo–Derived Health-Risk Outputs for Heavy-Metal Exposure From Edible Mushrooms
编辑推荐：
本研究对Acar等人（2026）先前发表的概率风险评估中得出的情景级健康风险输出进行了二次探索性分析。本手稿没有生成新的蘑菇样本、元素测量、暴露测量或从头蒙特卡洛模拟。源研究量化了五种可食野生蘑菇（Macrolepiota mastoidea、Bovista
本研究对Acar等人（2026）先前发表的概率风险评估中得出的情景级健康风险输出进行了二次探索性分析。本手稿没有生成新的蘑菇样本、元素测量、暴露测量或从头蒙特卡洛模拟。源研究量化了五种可食野生蘑菇（Macrolepiota mastoidea、Bovista aestivalis、Tricholoma fracticum、Helvella sp.和Rhizopogon roseolus）中的镉（Cd）、铅（Pb）、砷（As）和汞（Hg），并生成了成人和儿童暴露情景的概率非致癌风险输出。上游数据集基于五个野外采集的蘑菇标本，每个标本代表一个评估物种，并在源文章中报告了相应的采集地点和采集者编号。在本分析中，可用的蒙特卡洛衍生的Cd、Pb和As记录被重组为一个二次矩阵，包含150,000条Cd、Pb和As特异性情景记录，对应于每个蘑菇物种内每种金属的10,000次蒙特卡洛迭代。Hg未被纳入多变量剖面模型，因为可用的情景文件未提供与聚类工作流所需相同格式的完整Hg特异性风险向量，且Hg在源文章中包含Helvella sp.的一个未检出案例。估计每日摄入量（EDI）、金属特异性目标危害商数（THQ）和累积危害指数（HI/TTHQ）值被标准化，并使用主成分分析（PCA）、均匀流形近似与投影（UMAP）和K-means聚类进行检验。k=3–6的比较聚类表明，k=4提供了一个平衡的探索性解决方案，尽管轮廓系数（0.36–0.40）和戴维斯-布尔丁指数（0.83–0.93）指示中等而非强烈的分离。四个剖面的平均HI值范围为3.27至15.36，第95百分位HI值范围为5.09至41.72；因此，即使是相对最低的剖面也仍高于常规HI=1的非致癌关注阈值。Rhizopogon roseolus与较高累积风险/毒性压力剖面的关联最强（其情景的49.2%），而Macrolepiota mastoidea最常被分配到较低相对累积风险剖面（48.0%）。这些发现应被解释为描述性的、数据驱动的风险剖面，而非因果毒理学类别。
研究人员在论文《JOURNAL OF FOOD QUALITY》中发表了一项基于已有概率风险评估数据的二次探索性分析。研究背景方面，食用菌中重金属暴露的健康风险日益受到关注，传统确定性风险评估方法难以充分反映多因素暴露下的不确定性。蒙特卡洛模拟（Monte Carlo simulation）虽能处理不确定性，但其输出通常汇总为均值或百分位数，可能掩盖不同暴露情景下的风险子结构。现有研究缺乏对多金属、多物种情景级风险输出的系统性无监督分析，因此有必要开展二次剖面分析，以揭示累积风险模式并提高概率风险评估输出的可解释性。研究人员基于Acar等人（2026）已发表的概率风险评估，重新组织其生成的镉（Cd）、铅（Pb）、砷（As）的估计每日摄入量（EDI，Estimated Daily Intake）、目标危害商数（THQ，Target Hazard Quotient）和累积危害指数（HI/TTHQ，Hazard Index/Total Target Hazard Quotient）等情景级输出，共150,000条记录，来自五种土耳其食用野生蘑菇（Macrolepiota mastoidea、Bovista aestivalis、Tricholoma fracticum、Helvella sp.和Rhizopogon roseolus），样本采集于?anakkale和Van省。研究人员使用主成分分析（PCA，Principal Component Analysis）、均匀流形近似与投影（UMAP，Uniform Manifold Approximation and Projection）和K-means聚类进行无监督剖面分析，最终选择k=4的解决方案，并描述了四个探索性风险剖面。主要结论包括：所有四个剖面的平均HI均超过常规非致癌关注阈值（HI=1），其中Rhizopogon roseolus与高累积风险/毒性压力剖面（Profile 2）关联最强（49.2%的情景），Macrolepiota mastoidea与较低相对累积风险剖面（Profile 0）关联最强（48.0%），与源文章中的物种敏感砷（As）和HI/TTHQ模式一致。该研究的意义在于提供了一个可重复的二次分析框架，用于重新利用蒙特卡洛衍生的风险输出，通过数据驱动的方式描述相对风险模式，而非建立因果毒理学分类，从而增强食品安全和公共卫生风险评估中的专家解释能力。
关键技术方法方面，研究人员使用了三种主要方法：首先，PCA用于线性降维以诊断标准化情景矩阵的主方差结构；其次，UMAP用于非线性降维以捕捉局部情景相似性，参数设置为n_neighbors=15、min_dist=0.10、欧氏距离、二维嵌入；最后，K-means聚类在UMAP嵌入空间中进行分组，比较k=3–6后选择k=4。数据来源于Acar等人（2026）的蒙特卡洛模拟输出，包括从土耳其?anakkale和Van省采集的五种食用野生蘑菇的Cd、Pb、As的EDI、THQ和HI/TTHQ值，共150,000条记录。研究人员对数据进行标准化处理后，依次进行PCA、UMAP和K-means聚类，并进行了参数敏感性分析（n_neighbors为10、15、30，min_dist为0.00、0.10、0.30，random_state为0和42），使用轮廓系数（Silhouette score）、戴维斯-布尔丁指数（Davies–Bouldin index）、调整兰德指数（ARI，Adjusted Rand Index）和归一化互信息（NMI，Normalized Mutual Information）评估稳定性。
研究结果部分如下：
**3.1 降维结果（PCA和UMAP）**
通过PCA对标准化情景矩阵进行线性投影，结果显示情景空间存在异质性，部分情景占据更极端区域，但该结果仅作为探索性诊断，而非统计显著分离的证据。UMAP嵌入（n_neighbors=15，min_dist=0.10）提供了局部情景相似性的二维非线性表示，视觉上显示局部分组，但定量指标中等，因此UMAP仅作为描述性组织工具，支持后续剖面构建，而非独立统计证据。
**3.2 基于UMAP的K-means聚类结果**
对k=3–6进行K-means聚类，比较轮廓系数和Davies-Bouldin指数：k=3时轮廓系数最高（0.40），k=5时Davies-Bouldin指数最低（0.83）。保留k=4解决方案，因其避免了k=3的过度合并和k=5或6的过度碎片化，尽管未最大化任一指标。参数敏感性分析显示，轮廓系数范围为0.350–0.399，Davies-Bouldin指数为0.828–0.952，ARI与基线解相比为0.220–0.615，NMI为0.262–0.599，表明聚类结构对参数设置敏感，仅提供中等描述性分离，因此k=4被视为实用探索性总结而非参数不变的潜在毒理学类别。
**3.3 描述性风险剖面总结**
四个剖面（Profile 0–3）的描述性统计显示：平均HI分别为3.27、10.42、15.36和6.95，第95百分位HI分别为5.09、26.96、41.72和16.40，所有剖面均超过HI=1阈值。Profile 2具有最高累积风险，Profile 3具有最高平均金属特异性THQ（4.56）和平均EDI（1.76×10
-3 mg/kg/day），而Profile 0具有最低值。金属特异性THQ列为描述性均值，非直接重构剖面平均HI的算术组分，且未使用推理性检验比较剖面。
**3.4 风险剖面在蘑菇物种中的分布**
通过原始情景计数和百分比分布显示，各蘑菇物种均出现在多个剖面中，表明物种-剖面关系非一对一。Rhizopogon roseolus在Profile 2中情景最多（49.2%），Macrolepiota mastoidea在Profile 0中最多（48.0%），Bovista aestivalis最常分配至Profile 3（31.5%），Helvella sp.在Profile 0、1、3中分布较均匀，Tricholoma fracticum在Profile 0中最多（34.1%）。这些结果与源文章一致，其中Rhizopogon roseolus具有最高砷浓度和HI/TTHQ值，而Macrolepiota mastoidea砷负荷较低。
总结讨论部分，研究人员强调结果应视为二次风险剖面分析而非新健康风险评估。最强的物种-剖面关联（Rhizopogon roseolus与Profile 2）与源文章一致，且与其他食用菌研究（如Chen等、Qin等、Bucurica等、Okumu?等）的多指标解释相符。然而，聚类指标中等，缺乏Hg、铬（Cr）、镍（Ni）、铜（Cu）、锌（Zn）、锰（Mn）、铁（Fe）等元素，以及无外部验证，限制了毒理学广度。参数敏感性分析进一步支持谨慎解释。方法论价值在于探索性组织现有情景输出，而非预测建模或因果推断。
研究结论部分翻译如下：本手稿对继承了Acar等人（2026）的150,000个蒙特卡洛衍生的Cd-Pb-As情景记录进行了二次无监督剖面分析。保留了k=4的UMAP-K-means解决方案作为描述性折衷，尽管聚类指标显示中等分离。平均HI
