《Journal of Xenobiotics》：Space-Time Analysis of Burgeoning US Atrial Septal Defect Rates Driven by Cannabis
Albert Stuart Reece and
Gary Kenneth Hulse
编辑推荐：
本研究针对美国日益严峻的房缺（ASD）流行趋势，通过分析全国出生缺陷监测网络（NBDPN）2003-2020年数据，结合药物使用、收入及大麻合法化等多源信息，首次应用时空回归模型揭示了ASD率与大麻使用显著正相关，并证实大麻合法化状态是ASD高发风险因素。结果提示大麻暴露及其合法化政策是推动ASD高发的重要公共卫生问题。
在美国，一种称为房间隔缺损（Atrial Septal Defect, ASD）的先天性心脏病正以前所未有的速度蔓延，部分地区已影响到约1/11.3的新生儿。这种心脏结构的异常不仅威胁儿童健康，其背后驱动因素的未知性更成为公共卫生领域的紧迫课题。近年多国大型研究暗示，日益普遍的大麻使用可能与ASD发病率上升存在关联，然而传统的流行病学方法难以揭示两者在时空维度上的动态关联。为厘清大麻、其合法化政策与ASD流行之间是否存在因果关系，以及是否存在非吸入性暴露途径等问题，研究人员Albert Stuart Reece与Gary Kenneth Hulse在《Journal of Xenobiotics》上发表了此项开创性研究。
研究整合了2003年至2020年间美国27个州、总计17,751,437名活产儿中110,107例ASD病例的数据。通过获取国家出生缺陷预防网络（NBDPN）的ASD率（ASDR）信息，并结合国家药物使用与健康调查（NSDUH）中的香烟、酒精、大麻、镇痛药、可卡因使用数据、美国人口普查收入数据及大麻合法化状态，研究团队首次将时空分析、因果推断工具与多元统计模型结合，系统评估了物质使用、社会经济因素与ASD风险的关联。
关键技术方法
研究核心为空间回归分析与因果推断。利用R软件包splm进行空间回归建模，以处理数据的空间自相关性。通过构建时空序列单变量与双变量地图，可视化ASDR与各物质使用的空间分布及动态变化。为控制混杂偏倚，研究采用了逆概率加权（IPW）的广义回归模型，并使用E值评估未测量混杂因素需达到的强度才能推翻观察到的关联。此外，为验证结果的稳健性，研究进行了非重叠数据集分析、时间滞后分析、剔除高值州（内华达州）的敏感性分析，并利用Belsley-Kuh-Welsch共线性诊断评估模型的多重共线性问题。数据缺失部分（7.6%）通过线性回归插补处理。
研究结果
3.1. 序列地图-图形分析
时序地图显示，ASDR的空间分布模式与大麻使用指标高度一致，而非其他物质（如香烟、酒精、可卡因）。展示2005年至2018年美国（对数）ASDR分布，2018年高发区包括内华达、新墨西哥、肯塔基、田纳西、密歇根、密苏里等州。大麻使用地图（图2）显示，高使用率集中在俄勒冈、内华达、科罗拉多、加利福尼亚等州。双变量地图（图3与图4）更直观揭示，ASDR高发区（紫色/粉红色）与大麻高使用区在空间上高度重叠，尤其是在2014年后的科罗拉多、内华达、俄勒冈、阿拉斯加等州。值得注意的是，中西部肯塔基、田纳西、密西西比、密苏里等州呈现高ASDR但大麻使用率不高的“异常”模式，提示可能存在吸入途径以外的暴露（如环境水体污染）。大麻合法化状态的时空分布（图5）与大麻使用模式高度相似，且在双变量分析中与ASDR显著相关。
3.2. 空间回归
多元空间回归模型经变量筛选后，仅大麻（包括Δ9-THC和大麻二酚CBD）使用与ASDR保持显著正相关（p=0.00002 for cannabidiol），而香烟、酒精、可卡因、镇痛药滥用及收入等协变量均被剔除。交互作用模型中，大麻相关项同样显著。
3.3. 合法化状态
大麻合法化状态与ASDR显著相关。在将合法化状态分为非法、非刑事化、医用、合法四类时，ASDR存在显著组间差异。当二分法（合法 vs. 其余）分析时，合法州的ASDR风险比（RR）为2.73（2.66, 2.80），人群归因分数（AFP）为3.55%，E值点估计为4.90（95% CI下限4.76），远超因果阈值1.25。时空回归中，合法状态作为因子也显示高度显著（p=0.0013）。
3.4. 多重共线性诊断
Belsley条件指数均低于10，表明模型中预测变量间不存在严重多重共线性，回归估计可靠。
3.5. 稳健性分析
剔除ASDR最高的内华达州后，空间回归结果与主分析基本一致。采用末次观察值结转法进行数据插补也证实了主要结论。
3.6. 调查回归
采用逆概率加权（IPW）的调查回归模型，在控制调查设计后，大麻与ASDR的关联依然高度显著（p<0.0001），且E值中位数达13.78，进一步支持关联的稳健性。
3.7. 非重叠数据集
使用2003–2007、2008–2012、2016–2020三个非重叠时间段的数据进行分析，结果仍显示大麻与大麻素与ASDR强烈相关。
3.8. 时间滞后
考虑妊娠生物学，将预测变量滞后一年分析，大麻与大麻素仍与ASDR显著关联。
3.9. 病例确认
ASD病例确认方式（主动、被动、混合）在各州间的差异与ASDR无显著统计学关联，表明报告偏倚并非主要混杂因素。
3.10. 敏感性分析
计算的44个E值中，79.5%的点估计超过9，全部超过因果阈值1.25；27.3%的95%置信区间下限超过9。E值中位数为13.78（IQR: 9.25, 123.62），其95% CI下限中位数为3.31。这些结果表明观察到的关联能抵抗较强的未测量混杂。
结论与讨论
本研究通过多维度分析证实，大麻使用（尤其是大麻二酚CBD、Δ9-THC和大麻萜酚CBG）及其合法化政策是美国ASD流行的重要驱动因素。时空分析显示两者分布高度同步，多元回归中仅大麻相关变量保留显著，且E值普遍较高，强烈提示因果关联。研究还发现中西部农业州（如肯塔基、田纳西）存在“高ASDR-低大麻使用”的异常模式，推测可能与大麻种植导致的环境污染（如水体污染）有关，这为后续研究指明了新方向。
从生物学机制看，心脏发育受多种形态发生素（如Sonic hedgehog、BMP、FGF、Notch、Wnt信号通路）和表观遗传调控，而大麻素已被证明能直接或表观遗传学抑制这些关键通路，从而干扰心房隔等心脏结构的正常形成。这为大麻的致畸性提供了合理的生物学解释。
研究意义深远。首先，ASD作为最常见的先天性心脏病之一，其发病率飙升至近1/11，已成为严重的公共卫生问题。其次，研究首次将时空分析与因果推断结合，系统揭示了大麻暴露与ASD的强关联，并指向大麻合法化政策的潜在公共卫生风险。最后，结果提示大麻的遗传毒性和表观遗传毒性可能具有跨代效应，呼吁社区在制定公共健康政策时，需严格控制大麻素暴露，以保护当代及后代的表观遗传健康。尽管本研究为生态学设计，无法在个体层面断定因果，但其多层次证据的收敛、高E值以及对多种混杂因素的严格控制，使其结论具有重要警示价值，并为后续机制研究与政策干预提供了关键科学依据。
