The commercial tea from different regions of China was sampled to detect different REEs' concentrations in various teas and analyze their distribution patterns. According to the 2013-2014 China National Food Consumption Survey, tea consumption data were obtained. Combined with the REE content and leaching ratio in tea, the REE exposure risk assessment of tea drinkers was carried out.
2.1 Sampling and determination
From 2012 to 2015, 4326 tea samples were collected from supermarkets and local markets in 24 provinces, autonomous regions and municipalities of China, including green tea, oolong tea, black tea, dark tea, white tea and flowering tea. Since only three samples, the yellow tea was not included in the study.
The concentrations of REEs in tea samples were analyzed in local laboratories in 24 provinces, autonomous regions and municipalities. All laboratories used the same analysis procedure and training was provided before the analysis. The sixteen REEs except for Pm in teas were analyzed following a protocol for inductive coupled plasma-mass spectrometry (ICP-MS) elemental analysis in China National Food Safety Standard GB5009.94-2012. Chemical data were independently verified to ensure data quality and accuracy. The limit of detection (LOD) for the method ranged from 0.003 to 0.6 μg/kg for different tea categories. All the samples with results below the LOD were calculated as 1/2LOD depending on the recommended method by the GEMS/Food (WHO 1995) [20]. The percentage of samples in different tea categories to each REE below LOD was over 60%, so the LOD value (UB, upper bound) was assigned to those non-detectable samples. Otherwise, a half of LOD was assigned to the non-detectable (MB, middle bound).
2.2 Tea consumption
Population intake rates for tea were derived from the China National Food Consumption Survey during 2013-2014. In the survey, 43,386 subjects were selected through stratified multi-stage cluster sampling from 16 provinces, autonomous regions, and municipalities in China. Since there are only 21 individuals in the tea-drinking population under 18, the tea-drinking population in this study only considers adults. Food consumption data were collected with a 24-hour dietary recall method on three discrete days, and body weights were obtained in each corresponding individual. Averages of each individual's food type consumption were calculated on three discrete days to represent chronic intake. Because the concentration of REEs in flavored tea beverages and tea infusion was hard to achieve, the consumption rate of tea consumers was derived from surveyed individuals reporting processed tea consumption. Consumption within the same category of processed tea was summed up for every investigated individual to match each tea category's REE concentration data. Tea high consumer group is defined as the 95th percentile of consumers.
2.3 REE exposure assessment of tea consumers
To estimate the risk of REE exposure attributable to tea consumption in adults, used the following formula to calculate the average daily dose (ADD) value:
ADD =∑([REEs]*IR*LR) / BW/1000
Where ADD= average daily dose (in μg/kg BW/d), [REEs] = mean total REEs concentration (in μg/kg) in the specific category of tea, IR = specific category tea intake rate (g/d), LR=specific category tea leaching rate, and BW = each individual body weight (in kg).
Leaching rate was considered in the respect that people only drink tea infusion. Based on food consumption data, white tea was not consumed by the surveyed individuals. Leaching rates of the five tea categories have to be obtained from the available literature. Because the conversion factor of REEs to REOs (rare earth oxides) is in the range of 1.137 to 1.27(1.228 (Ce-Ce2O3), 1.148 (Dy-Dy2O3), 1.143 (Er-Er2O3), 1.158 (Eu-Eu2O3), 1.153 (Gd-Gd2O3), 1.146 (Ho-H2O3), 1.173 (La-La2O3), 1.137 (Lu-Lu2O3), 1.166 (Nd-Nd2O3), 1.208 (Pr-Pr2O3), 1.534 (Sc-Sc2O3), 1.160 (Sm-Sm2O3), 1.176 (Tb-Tb2O3), 1.142 (Tm-Tm2O3), 1.270 (Y-Y2O3) and 1.139 (Yb-Yb2O3)), we assume that leaching rated designated as REOs can be considered as REEs. Leaching rates from the literature were calculated at the maximum values, not considering tea shape and tea origin. Leaching ratios of dark tea, oolong tea, green tea, flowering tea and black tea were 26.01%, 65.7%, 35%, 33.3% and 7.04%, respectively [19,21-23]. (See Additional file 1)
2.4 Health risk assessment of tea consumers
Hazard quotient (HQ) was applied to judge if the health risk is acceptable. If the value of HQ is less than one, it is believed that there is no significant health risk. If HQ exceeds one, there are possibilities that health risks occur, with a probability increasing as the value of HQ increases [24]. The following equation estimated the HQ:
HQ=ADD /ADI
Where ADI is the acceptable daily dose (μg/kg BW/d).
In 2016, the China Scientific Committee on Food Safety Risk Assessment set the temporary ADI (tADI) of three main REE, lanthanum (51.5μg/kg BW), cerium (645.0μg/kg BW), and yttrium (145.5μg/kg BW) based on the NOAEL of the critical effect in the decreased body weight gain from a 90-day feeding study in rats and application of a safety factor of 200 for inter-species and intra-human variability [25-26]. Besides, the China Scientific Committee selected the lowest tADI of lanthanum as the health-based guidance value of the total REEs following the conservation principle of risk assessment, then derived a group tADI for the total REEs (51.5μg/kg BW). Therefore, we use 51.5μg/kg BW/d as ADI in our study.
2.5 Statistic analysis
REEs concentrations were described by mean, median, interquartile range (IQR), and quantiles. The distributions of total REEs in different categories of tea were positive-skewed, so the median could better reflect the REEs level. Considering the distributions of REEs levels did not follow a normal distribution, Kruskal-Wallis tests were used to assess statistics difference among categories of tea samples and REE exposure of various subgroups, respectively. Moreover, comparisons between multiple groups were carried out by the Wilcoxon test, adjusted by the p.adj function. Statistical significance was two-tailed and set at α=0.05. All statistical analyses and figures were performed using R 3.5.1.