多環(huán)芳烴聯(lián)合暴露與肥胖的關(guān)聯(lián)分析

[摘要] 目的 探索多環(huán)芳烴（polycyclic aromatic hydrocarbon，PAH）聯(lián)合暴露與肥胖之間的關(guān)聯(lián)。方法 使用2011—2016年美國國家健康與營養(yǎng)調(diào)查（National Health and Nutrition Examination Survey，NHANES）的數(shù)據(jù)進行調(diào)查分析，以PAH尿液代謝物評估環(huán)境PAH暴露量。使用限制性立方樣條探索PAH暴露與肥胖的非線性劑量–反應關(guān)系。采用加權(quán)分位數(shù)和（weighted quantile sum，WQS）回歸和基于分位數(shù)的g計算（quantile g-computation，Qgcomp）模型評估PAH聯(lián)合暴露對肥胖的影響。結(jié)果 在WQS正向模型中，PAH聯(lián)合暴露與肥胖之間呈正相關(guān)（Plt;0.05）。未在WQS負向模型中發(fā)現(xiàn)顯著相關(guān)性（Pgt;0.05）。進一步使用Qgcomp模型驗證混合物的總體效應，發(fā)現(xiàn)PAH聯(lián)合暴露依舊與肥胖呈正相關(guān)（Plt;0.05）。結(jié)論 PAH聯(lián)合暴露可增加肥胖的風險，PAH聯(lián)合暴露主要歸因于2-羥基萘。

[關(guān)鍵詞] 多環(huán)芳烴；多環(huán)芳烴代謝物；聯(lián)合暴露；肥胖

[中圖分類號] R122.7" """"[文獻標識碼] A """""[DOI] 10.3969/j.issn.1673-9701.2025.19.005

Association analysis of combined exposure to polycyclic aromatic hydrocarbon and obesity

ZHOU Guoxiu， WEI Wanting， YAN Yizhong

Department of Preventive Medicine， School of Medicine， Shihezi University， Shihezi 832003， Xinjiang， China

[Abstract] Objective To explore the association between combined exposure to polycyclic aromatic hydrocarbons （PAH） and obesity. Methods The data from the National Health and Nutrition Examination Survey （NHANES） from 2011 to 2016 were used for investigation and analysis， and the urine metabolites of PAH were used to evaluate the environmental PAH exposure. The nonlinear dose-response relationship between PAH exposure and obesity was explored using restricted cubic splines. The weighted quantile sum （WQS） and quantile g-computation （Qgcomp） models were used to evaluate the impact of combined exposure to PAH on obesity. Results In the WQS positive model， there was a positive correlation between combined exposure to PAH and obesity （Plt;0.05）. No significant correlation was found in the WQS negative model （Pgt;0.05）. The overall effect of the mixture was further verified using the Qgcomp model， and it was found that combined exposure to PAH was still positively correlated with obesity （Plt;0.05）. Conclusion Mixed exposure to PAH increases the risk of obesity， and combined exposure to PAH is mainly attributed to 2-hydroxynaphthalene.

[Key words] Polycyclic aromatic hydrocarbon; Polycyclic aromatic hydrocarbon metabolites; Combined exposure; Obesity

肥胖已成為全球普遍問題，世界衛(wèi)生組織最新調(diào)查顯示，2000—2021年全球肥胖率增長1.5倍^[1]。肥胖可引起高血壓、血脂異常、糖尿病等代謝性疾病，對健康構(gòu)成嚴重威脅^[2-3]。目前，非傳統(tǒng)因素的環(huán)境污染物被越來越多關(guān)注。多環(huán)芳烴（polycyclic aromatic hydrocarbon，PAH）作為長期暴露的空氣污染物在生活中廣泛存在^[4-5]。PAH可通過影響正常內(nèi)分泌功能，刺激脂肪細胞分化，影響體脂分布^[6]。流行病學研究已證明PAH代謝物與肥胖的關(guān)系，但只評估單一PAH代謝物的暴露風險^[7]。污染物聯(lián)合暴露的傳統(tǒng)研究方法是以變量為中心，基于變量之間的相關(guān)性進行分析，但由于共同暴露的化學物存在潛在共線性與復雜的交互作用，污染物與結(jié)局之間可能存在非線性與非加性的關(guān)系，故傳統(tǒng)分析方法面對高階暴露變量數(shù)據(jù)集較為有限，無法判斷變量間相加、協(xié)同或拮抗效應^[8]。采用加權(quán)分位數(shù)和（weighted quantile sum，WQS）回歸和基于分位數(shù)的g計算（quantile g-computation，Qgcomp）模型作為聯(lián)合暴露統(tǒng)計模型，不僅能得到混合效應評估總分，還可識別對總體效應起最大作用的污染物^[9]。本研究使用2011—2016年美國國家健康與營養(yǎng)調(diào)查（National Health and Nutrition Examination Survey，NHANES）的數(shù)據(jù)，旨在探討6～85歲參與者PAH混合暴露后的尿液代謝物與肥胖的關(guān)系，并使用WQS回歸和Qgcomp模型評估聯(lián)合暴露效應，為肥胖的防控提供流行病學證據(jù)。

1" 資料與方法

1.1 "研究人群選擇

本研究使用NHANES 2011—2012年、2013—2014年、2015—2016年的調(diào)查數(shù)據(jù)，包括8033名年齡6～85歲有PAH檢測信息的參與者。排除教育水平（n=476）、家庭貧困與收入比（poverty income ratio，PIR）（n=668）、血清可替寧（n=557）、飲食質(zhì)量（n=1224）數(shù)據(jù)缺失者，最終共有5108名參與者被納入分析。其中有5073名參與者提供體質(zhì)量指數(shù)（body mass index，BMI）用于研究肥胖，有4999名參與者提供腰圍用于研究腹型肥胖，有3249名參與者提供全身脂肪量。

1.2 "肥胖

測量數(shù)據(jù)與問卷數(shù)據(jù)均由專業(yè)人員進行收集。BMI=體質(zhì)量（kg）/身高（m）²；暴露髂嵴用卷尺測量腰圍，精確到0.1cm；采用雙能X線吸收法檢測全身脂肪量（kg）。BMI≥30kg/m²為肥胖^[1]。定義女性腰圍≥88cm、男性腰圍≥102cm為腹型肥胖^[10]。

1.3 "PAH暴露評估

NHANES三分之一的參與者采集PAH尿液代謝物作為PAH的生物標志物，將樣品置于–20℃下儲存運送至國家環(huán)境健康中心進行測試分析。本研究共分析NHANES 3個周期的6種PAH尿液代謝物，分別為1-羥基萘（hydroxynaphthalene，NAP）、2-NAP、2-羥基芴（hydroxyfluorene，F(xiàn)LUO）、3-FLUO、1-羥基菲（1-hydroxyphenanthrene，1-PHEN）和1-羥基芘（1-hydroxypyrene，1-PYR）。6種PAH尿液代謝物的檢出率均gt;90%，對低于檢出限（limit of detection，LOD）的PAH尿液代謝物濃度，可用檢出限除以2的平方根替代。

1.4 "協(xié)變量

分析中涉及的分類變量有性別、年齡、種族、教育程度。血清可替寧可作為主動吸煙和被動吸煙的指征，濃度≥0.015ng/ml表示暴露于煙草煙霧環(huán)境^[11]。PIR可作為家庭收入指標，在本研究中被編碼分為3個等級，分別為lt;1.3、1.3～3.5和gt;3.5?；诿恐苌眢w活動的總小時數(shù)，定義體力為低體力活動（0MET）、中等體力活動（1～48MET）和高強度體力活動（gt;48MET）^[12]。根據(jù)健康飲食指數(shù)–2020得分將其分為低（lt;42分）、中（42～59分）和高（gt;59分）3個級別；得分越高表明飲食質(zhì)量越好。

1.5 "統(tǒng)計學方法

采用R 4.2.3軟件進行統(tǒng)計學分析。在描述性分析中，分類變量使用頻數(shù)和比例描述，比較采用χ²檢驗；正態(tài)分布的連續(xù)變量以均數(shù)±標準差（"）表示，比較采用方差分析。使用限制性立方樣條探索PAH暴露與肥胖的非線性劑量–反應關(guān)系。使用WQS回歸和Qgcomp模型研究PAH混合暴露對肥胖的聯(lián)合影響。Plt;0.05為差異有統(tǒng)計學意義。

本研究中WQS回歸模型將每個PAH代謝物濃度分為四個有序變量（四分位數(shù)），通過計算線性加權(quán)指數(shù)評估6種PAH代謝物濃度對肥胖的混合效應，計算每個PAH代謝物濃度的四分位數(shù)與結(jié)局的關(guān)聯(lián)強度，代表每種PAH代謝物對總體效應指數(shù)的貢獻程度。由于WQS回歸模型可放大未知混雜因素的偏差，在評估污染物與結(jié)局之間的關(guān)系時有方向同質(zhì)性、線性、可加性限制，故本研究引入Qgcomp模型對WQS回歸模型進行修正與補充。該模型可評估不同方向上總體混合效應，同時計算正負權(quán)重。它不僅消除WQS回歸的方向同質(zhì)性，其本身靈活的g計算也增大結(jié)果的穩(wěn)健性^[13]。

2" 結(jié)果

2.1 "參與者的基本特征

本研究納入5108名參與者，平均年齡（37.86±22.09）歲，女性略多于男性，20～59歲年齡段人群約占50%，其基本特征見表1。共納入5073名參與者分析PAH和肥胖的關(guān)系，與非肥胖者相比，肥胖者女性、20～59歲、非西班牙裔白人、高中以上文化水平、暴露于煙草煙霧環(huán)境、中等體力活動、中等飲食質(zhì)量者占比均較高（Plt;0.05），見表2。共納入4999名參與者分析PAH和腹型肥胖的關(guān)系，與非腹型肥胖者相比，腹型肥胖者女性、20～59歲、非西班牙裔白人、高中以上文化水平、中等體力活動的占比均較高（Plt;0.05），見表3。共納入3249名參與者分析PAH和全身脂肪量的關(guān)系，女性、20～59歲、非西班牙裔黑人、高中以上文化水平、暴露于煙草煙霧環(huán)境、中等體力活動者的全身脂肪量更高（Plt;0.05），見表4。

2.2 "非線性關(guān)系分析

2-NAP、2-FLUO和1-PYR與肥胖和腹型肥胖之間均存在顯著的非線性關(guān)系（Plt;0.05）。2-NAP和2-FLUO與全身脂肪量之間存在顯著的非線性關(guān)系（Plt;0.05），見圖1。

2.3 "WQS回歸模型

在WQS回歸正向模型中，PAH聯(lián)合暴露與肥胖之間存在關(guān)聯(lián)（Plt;0.05）。WQS回歸指數(shù)每增加一個四分位數(shù)，肥胖、腹型肥胖、全身脂肪量風險分別增加1.37倍（95%CI：1.24～1.52）、1.31倍（95%CI：1.19～1.45）、5.65倍（95%CI：2.90～11.01），見圖2。WQS回歸可通過權(quán)重代表混合物中各代謝物的重要性，在所有結(jié)局中，PAH的混合暴露效應主要歸因于2-NAP，其在肥胖、腹型肥胖、全身脂肪量的權(quán)重分別為84.60%、82.30%、65.00%，其次是1-PHEN，分別為14.00%、14.50%、34.20%，見表5。

2.4 "Qgcomp模型

進一步使用Qgcomp模型驗證混合物的總體效應。研究發(fā)現(xiàn)PAH聯(lián)合暴露依舊與肥胖呈正相關(guān)，Qgcomp指數(shù)每增加一個四分位數(shù)，肥胖、腹型肥胖、全身脂肪量風險分別增加1.16倍（95%CI：1.05～ 1.27）、1.25倍（95%CI：1.14～1.37）和1.75倍（95%CI：1.00～3.05），見圖3。在所有肥胖指標中，2-NAP依舊是最大的正權(quán)重，3-FLUO為最大負相關(guān)權(quán)重，見表5。

3" 討論

WQS回歸和Qgcomp模型表明，PAH聯(lián)合暴露顯著增加肥胖的風險，全人群中PAH聯(lián)合暴露與肥胖有正向聯(lián)系。研究發(fā)現(xiàn)PAH化合物具有親脂性，易在脂肪組織中積累，血漿和尿液中PAH水平可能只是冰山一角，在脂肪組織樣本中發(fā)現(xiàn)的PAH濃度要高得多^[14]。一項動物實驗發(fā)現(xiàn)，母鼠產(chǎn)前PAH暴露可導致后代體質(zhì)量和脂肪增加，引起環(huán)加氧酶、脂聯(lián)素等決定脂肪細胞分化和功能的基因表達異常^[15]。體外實驗發(fā)現(xiàn)，PAH可發(fā)揮內(nèi)分泌干擾作用，表現(xiàn)為甲狀腺激素受體拮抗性，抑制甲狀腺激素對脂肪的分解代謝作用^[16]。一項關(guān)于PAH與兒童肥胖的NHANES研究顯示，兒童肥胖與PAH（芴、萘、菲和芘）代謝物分子質(zhì)量總和呈正相關(guān)^[17]。韓國一項研究顯示PAH尿液代謝物與韓國成年人群的BMI呈正相關(guān)^[18]。這些與本研究混合暴露結(jié)果一致，但都僅限于將單一的PAH相加；而在本研究中，污染物的非可加性和非線性偏倚得到控制。

構(gòu)建WQS回歸和Qgcomp模型進一步發(fā)現(xiàn)在混合暴露中2-NAP占主要的正權(quán)重。2-NAP作為PAH代謝物具有一定的氧化應激潛力。研究表明2-NAP可通過生成活性氧誘導氧化應激，導致細胞損傷和炎癥反應^[19]。細胞實驗發(fā)現(xiàn)長期暴露于2-NAP可刺激脂肪細胞的脂質(zhì)積累^[20]。Ranjbar等^[21]研究顯示2-NAP可增加成人肥胖風險，與本研究結(jié)果一致。2-NAP是萘最重要的代謝產(chǎn)物，而萘是最簡單的二環(huán)PAH，極易揮發(fā)，常以氣相存在于空氣中，是吸煙和暴露于煙草煙霧環(huán)境最常接觸的PAH^[22-23]。煙草煙霧和PAH對肥胖存在協(xié)同作用，同時暴露可大大增加肥胖風險^[24]。萘是驅(qū)蛾劑和樟腦丸的關(guān)鍵活性物質(zhì)，美國41%的驅(qū)蛾劑產(chǎn)品含萘，棉質(zhì)衣服可吸收樟腦丸中的高濃度萘成為次要暴露來源^[22]。

本研究亦存在不足。首先，橫斷面的研究設計無法推斷PAH和肥胖的因果關(guān)系，需要更大規(guī)模的前瞻設計研究和毒理學實驗驗證。其次，研究結(jié)果的解釋僅限于特定的NHANES周期和有限的PAH代謝物。最后，橫斷面設計無法排除反向因果關(guān)系，建議未來研究采用縱向隊列設計，明確暴露與結(jié)局的時間關(guān)系。

通過WQS回歸和Qgcomp模型發(fā)現(xiàn)PAH共同暴露可對肥胖產(chǎn)生顯著的正向效應，其中2-NAP占主要的正權(quán)重。

利益沖突：所有作者均聲明不存在利益沖突。

[參考文獻]

[1]"" World Health Organization. Obesity and overweight[EB/OL]. （2025-05-07）[2025-06-11]. https：//www.who.int/news-room/ fact-sheets/detail/obesity-and-overweight.

[2]"" 張曉慧， 彭芯， 譚健霞， 等. 重慶市渝北區(qū)老年人TyG指數(shù)及腹型肥胖與初篩高血壓的關(guān)系研究[J]. 中國實驗診斷學， 2024， 28（11）： 1261–1266.

[3]"" 王鳳華， 楊寧， 王雄關(guān)， 等. 老年血脂異?；颊唧w重指數(shù)與高血壓患病風險的關(guān)系[J]. 中華高血壓雜志（中英文）， 2025， 33（4）： 370–377.

[4]"" 程強， 亢潤潤， 張海蕾， 等. 2017—2022年臺州市某城區(qū)空氣PM2.5中重金屬和多環(huán)芳烴分布特征及來源分析[J]. 環(huán)境與職業(yè)醫(yī)學， 2024， 41（10）： 1144–1150.

[5]"" 李文金， 董發(fā)勤， 呂珍珍， 等. 典型工業(yè)城市德陽大氣多環(huán)芳烴的污染特征及致癌風險評估[J]. 環(huán)境與健康雜志， 2024， 41（10）： 847–851.

[6]"" MLYCZY?SKA E， BONGRANI A， RAME C， et al. Concentration of polycyclic aromatic hydrocarbons （PAHs） in human serum and adipose tissues and stimulatory effect of naphthalene in adipogenesis in 3T3-L1 cells[J]. Int J Mol Sci， 2023， 24（2）： 1455.

[7]"" KIM S H， PARK M J， PARK S K. Urinary concentrations of polycyclic aromatic hydrocarbon metabolites and childhood obesity[J]. Heliyon， 2023， 9（9）： e19335.

[8]"" HENDRYX M， LUO J. Latent class analysis of the association between polycyclic aromatic hydrocarbon exposures and body mass index[J]. Environ Int， 2018， 121（Pt 1）： 227–231.

[9]"" FANG L， ZHAO H， CHEN Y， et al. The combined effect of heavy metals and polycyclic aromatic hydrocarbons on arthritis， especially osteoarthritis， in the U.S. adult population[J]. Chemosphere， 2023， 316： 137870.

[10] GRUNDY S M， CLEEMAN J I， DANIELS S R， et al. Diagnosis and management of the metabolic syndrome： An American Heart Association/National Heart， Lung， and Blood Institute Scientific Statement[J]. Circulation， 2005， 112（17）： 2735–2752.

[11] LUO K， LIU J， WANG Y， et al. Associations between organophosphate esters and sex hormones among 6-19-year old children and adolescents in NHANES 2013—2014[J]. Environ Int， 2020， 136： 105461.

[12] AINSWORTH B E， HASKELL W L， HERRMANN S D， et al. 2011 compendium of physical activities： A second update of codes and MET values[J]. Med Sci Sports Exerc， 2011， 43（8）： 1575–1581.

[13] KEIL A P， BUCKLEY J P， O’BRIEN K M， et al. A quantile-based g-computation approach to addressing the effects of exposure mixtures[J]. Environ Health Perspect， 2020， 128（4）： 47004.

[14] SOUSA S， PAíGA P， PESTANA D， et al. Evaluating the impact of polycyclic aromatic hydrocarbon bioaccumulation in adipose tissue of obese women[J]. Chemosphere， 2024， 353： 141673.

[15] YAN Z， ZHANG H， MAHER C， et al. Prenatal polycyclic aromatic hydrocarbon， adiposity， peroxisome proliferator-activated receptor （PPAR） γ methylation in offspring， grand-offspring mice[J]. PLoS One， 2014， 9（10）： e110706.

[16] SUN H， SHEN O X， XU X L， et al. Carbaryl， 1-naphthol and 2-naphthol inhibit the beta-1 thyroid hormone receptor-mediated transcription in vitro[J]. Toxicology， 2008， 249（2-3）： 238–242.

[17] SCINICARIELLO F， BUSER M C. Urinary polycyclic aromatic hydrocarbons and childhood obesity： NHANES （2001-2006）[J]. Environ Health Perspect， 2014， 122（3）： 299–303.

[18] LEE I， PARK H， KIM M J， et al. Exposure to polycyclic aromatic hydrocarbons and volatile organic compounds is associated with a risk of obesity and diabetes mellitus among Korean adults： Korean National Environmental Health Survey （KoNEHS） 2015—2017[J]. Int J Hyg Environ Health， 2022， 240： 113886.

[19] NIE J， LI J， CHENG L， et al. Maternal urinary 2-hydroxynaphthalene and birth outcomes in Taiyuan， China[J]. Environ Health， 2018， 17（1）： 91.

[20] BRIGHT A， LI F， MOVAHED M， et al. Chronic exposure to low-molecular-weight polycyclic aromatic hydrocarbons promotes lipid accumulation and metabolic inflammation[J]. Biomolecules， 2023， 13（2）： 196.

[21] RANJBAR M， ROTONDI M A， ARDERN C I， et al. Urinary biomarkers of polycyclic aromatic hydrocarbons are associated with cardiometabolic health risk[J]. PLoS One， 2015， 10（9）： e0137536.

[22] JIA C， BATTERMAN S. A critical review of naphthalene sources and exposures relevant to indoor and outdoor air[J]. Int J Environ Res Public Health， 2010， 7（7）： 2903–2939.

[23] 李艷寧， 鄧云珺， 楊玲， 等. 被動吸煙對孕婦尿中4種多環(huán)芳烴代謝物濃度的影響[J]. 環(huán)境與職業(yè)醫(yī)學， 2015， 32（6）： 522–527.

[24] KIM H W， KAM S， LEE D H. Synergistic interaction between polycyclic aromatic hydrocarbons and environmental tobacco smoke on the risk of obesity in children and adolescents： The U.S. National Health and Nutrition Examination Survey 2003—2008[J]. Environ Res， 2014， 135： 354–360.

（收稿日期：2025–01–22）

（修回日期：2025–06–13）

通信作者：閆貽忠，電子信箱：Y202403students@163.com