《Journal of Hazardous Materials Advances》:Association between exposure to brominated flame retardants and thyroid function in NHANES 2007–2012
编辑推荐:
本研究首次采用广义线性模型(GLM)、主成分分析(PCA)、加权分位数和(WQS)回归及分位数G计算(QGC),在NHANES中评估溴化阻燃剂(BFRs)对甲状腺功能的单独及联合影响。由PBB-153驱动的模式(PC2)与游离三碘甲状腺原氨酸(FT3
本研究首次采用广义线性模型(GLM)、主成分分析(PCA)、加权分位数和(WQS)回归及分位数G计算(QGC),在NHANES中评估溴化阻燃剂(BFRs)对甲状腺功能的单独及联合影响。由PBB-153驱动的模式(PC2)与游离三碘甲状腺原氨酸(FT3)、总三碘甲状腺原氨酸(TT3)和总甲状腺素(TT4)呈稳健负相关。WQS回归揭示BFR混合物与FT3、游离甲状腺素(FT4)和TT3显著负相关,QGC进一步证实与FT4的显著负相关。
溴化阻燃剂(BFRs)是一类用于满足建筑材料、家具及电子产品防火安全标准的卤代有机化合物,人类主要通过摄入脂质丰富食物、粉尘摄入或吸入而暴露。多溴联苯(PBBs)和多溴联苯醚(PBDEs)是常用的添加型阻燃剂,具有极高环境持久性和长人体半衰期(如PBB-153约30年)。尽管多种BFRs已被《斯德哥尔摩公约》列为持久性有机污染物(POPs),其健康风险仍受关注。甲状腺激素,包括游离或总甲状腺素(FT
4/TT
4)和三碘甲状腺原氨酸(FT
3/TT
3),对发育、代谢和分化至关重要,其微小变化即可敏感反映化学干扰。动物研究提示BFRs可干扰甲状腺功能,但人群证据不一致且多基于单一暴露模型,难以反映复合暴露现实。为此,研究人员整合广义线性模型(GLMs)、主成分分析(PCA)、加权分位数和(WQS)回归和分位数G计算(QGC)四种方法,于美国国家健康与营养调查(NHANES)2007–2012数据中同时评估单独与联合的BFR暴露对甲状腺功能的影响。
该研究为横断面分析,最终纳入4,064名具有完整BFRs、甲状腺功能指标及关键协变量的成年人。血清中6种BFRs(PBB-153和PBDE-28、-47、-99、-100、-153)采用自动化液液萃取结合同位素稀释气相色谱-高分辨质谱法测定;甲状腺功能指标包括促甲状腺激素(TSH)、FT
3、FT
4、TT
3和TT
4。统计分析中,GLMs评估单个BFR与各甲状腺激素的关联,并经Benjamini-Hochberg法控制错误发现率(FDR);PCA将6种BFRs降维为共暴露模式;WQS和QGC用于估计BFR混合物的联合效应并识别主要贡献化学物。此外还进行了性别和年龄分层敏感性分析。
**基线特征(Baseline Features)**:研究初始纳入30,442名受试者,排除26,333名数据缺失者和45名孕妇后,4,064人进入最终分析。年龄、性别、种族、BMI、教育程度、饮酒、吸烟、甲状腺疾病及用药、心血管疾病、糖尿病和高血压均与一项或多项甲状腺参数显著相关。
**BFRs分布与相关(Distribution and Correlation of BFRs)**:分析显示PBDE-47与PBDE-99(r=0.95)和PBDE-100(r=0.92)高度正相关,PBDE-28与PBDE-47(r=0.89)相关;PBB-153与各PBDEs相关性较弱(r=0.15–0.23),提示其作为独立暴露模式。
**单个BFRs与甲状腺功能关联(Associations between Individual BFRs and Thyroid Function)**:经调查加权多元线性回归和FDR校正后,PBB-153与FT
3(β=-0.048;95%CI -0.070至-0.025;q=0.0095)及TT
3(β=-1.940;95%CI -3.056至-0.823;q=0.0349)显著负相关。在名义显著性水平下,PBDE-28与FT
3、FT
4、TT
3负相关,PBDE-99与TT
3负相关,PBB-153与FT
4负相关,而PBDE-47与FT
4正相关。
**暴露模式与甲状腺功能(PCA, WQS, QGC分析)**:PCA提取三个主成分,共解释94.4%变异。PC1(66.2%)代表PBDE-28、-47、-99、-100共暴露;PC2(16.2%)由PBB-153驱动;PC3(12.0%)由PBDE-153驱动。FDR校正后,PC2与FT
3(β=-0.062)、TT
3(β=-3.185)和TT
4(β=-0.098)显著负相关;PC3与FT
3(β=0.066)和TT
3(β=2.180)正相关,但与TT
4(β=-0.164)负相关;PC1无显著关联。WQS回归显示BFR混合物与FT
3(β=-0.042)、FT
4