S1 File -

Dietary uptake is the main pathway of exposure to polycyclic aromatic hydrocarbons (PAHs). However, there is no data regarding the pollution and health risks posed by PAHs in Lilium davidii var. unicolor. We measured the concentrations of 16 PAHs in lily bulbs from Lanzhou; analyzed the bioaccumulation, sources, and pollution pathways of PAHs; assessed the influence of baking on PAH pollution in the bulb; and assessed the cancer risks associated with PAH exposure via lily consumption. The total PAH concentrations in raw bulbs were 30.39–206.55 μg kg-1. The bioconcentration factors of total PAHs ranged widely from 0.92 to 5.71, with a median value of 2.25. Pearson correlation analysis revealed that the octanol-water partition coefficients and water solubility values played important roles in the bioaccumulation of naphthalene, fluorene, phenanthrene, pyrene, and fluoranthene in the raw bulb by influencing PAH availability in soil. Correlation analysis and principal component analysis with multivariate linear regression indicated that biomass and wood burning, coal combustion, diesel combustion, and petroleum leakage were the major sources of PAHs in the raw bulbs. The paired t-test showed that the PAH concentrations in the baked bulbs were higher than those in the raw bulbs. PAH compositions in lily bulb changed during the baking process. Baked bulbs exhibited a higher cancer risk than raw bulbs. Local adults had low carcinogenic risks from consuming lily bulbs. This study fills the knowledge gap about PAH pollution and the related health risks of PAHs in the Lanzhou lily.