BackgroundAnimal studies have shown that exposure to REEs can cause severe liver damage, but evidence from population studies is still lacking. Therefore, we investigated the relationship between REEs concentrations in urine and liver function in the population. MethodsWe conducted a cross-sectional study on 1024 participants in Nanning, China. An inductively coupled plasma mass spectrometer (ICP-MS) was used to detect the concentrations of 12 REEs in urine. The relationship between individual exposure to individual REE and liver function was analyzed by multiple linear regression. Finally, the effects of co-exposure to 5 REEs on liver function were assessed by a weighted sum of quartiles (WQS) regression model and a Bayesian kernel machine regression (BKMR) model. ResultsThe detection rate of 5 REEs, lanthanum (La), cerium (Ce), praseodymium (Pr), neodymium (Nd), and dysprosium (Dy), is greater than 60%. After multiple factor correction, La, Ce, Pr, Nd, and Dy were positively correlated with serum ALP, Ce, Pr, and Nd were positively correlated with serum AST, while Ce was negatively correlated with serum TBIL and DBIL. Both WQS and BKMR results indicate that the co-exposure of the 5 REEs is positively correlated with serum ALP and AST, while negatively correlated with serum DBIL. There were potential interactions between La and Ce, La and Dy in the association of co-exposure of the 5 REEs with serum ALP. ConclusionsThe co-exposure of the 5 REEs was positively correlated with serum ALP and AST, and negatively correlated with serum DBIL.