- 1Department of Laboratory Medicine, Shenzhen Institute of Translational Medicine, The First Affiliated Hospital of Shenzhen University, Shenzhen Second People’s Hospital, Shenzhen Key Laboratory of Medical Laboratory and Molecular Diagnostics, Shenzhen, China
- 2Shenzhen Yantian District People’s Hospital (Group), Southern University of Science and Technology Yantian Hospital, Shenzhen, China
Purpose: Non-alcoholic fatty liver disease (NAFLD) is the most common liver disease. Phthalates have been suggested to influence the development of NAFLD due to their endocrine-disrupting properties, but studies based on nationally representative populations are insufficient, and existing studies seem to have reached conflicting conclusions. Due to changes in legislation, the use of traditional phthalates has gradually decreased, and the phthalates substitutes is getting more attention. This study aims to delve deeper into how the choice of diagnostic approach influences observed correlations and concern about more alternatives of phthalates, thereby offering more precise references for the prevention and treatment of NAFLD.
Methods: A cohort of 641 participants, sourced from the National Health and Nutrition Examination Survey (NHANES) 2017–2018 database, was evaluated for NAFLD using three diagnostic methods: the Hepatic Steatosis Index (HSI), the US Fatty Liver Indicator (US.FLI), and Vibration Controlled Transient Elastography (VCTE). The urinary metabolite concentrations of Di-2-ethylhexyl phthalate (DEHP), Di-isodecyl phthalate (DIDP), Di-isononyl phthalate (DINP), Di-n-butyl phthalate (DnBP), Di-isobutyl phthalate (DIBP), Di-ethyl phthalate (DEP) and Di-n-octyl phthalate (DnOP) were detected. The association between NAFLD and urinary phthalate metabolites was evaluated through univariate and multivariate logistic regression analyses, considering different concentration gradients of urinary phthalates.
Results: Univariate logistic regression analysis found significant correlations between NAFLD and specific urinary phthalate metabolites, such as Mono-(2-ethyl-5-oxohexyl) phthalate (MEOHP), Mono-(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP), Mono-2-ethyl-5-carboxypentyl phthalate (MECPP), and Mono-(carboxyisoctyl) phthalate (MCiOP), across different diagnostic criteria. In a multivariate logistic regression analysis adjusting only for demographic data, MEOHP (OR = 3.26, 95% CI = 1.19–8.94, p = 0.029), MEHHP (OR = 3.98, 95% CI = 1.43–11.1, p = 0.016), MECPP (OR = 3.52, 95% CI = 1.01–12.2, p = 0.049), and MCiOP (OR = 4.55, 95% CI = 1.93–10.7, p = 0.005) were positively related to NAFLD defined by HSI and VCTE. The correlation strength varied with the concentration of phthalates, indicating a potential dose–response relationship. Adjusting for all covariates in multivariate logistic regression, only MCiOP (OR = 4.22, 95% CI = 1.10–16.2, p = 0.044), as an oxidative metabolite of DINP, remained significantly associated with NAFLD under the VCTE criterion, suggesting its potential role as a risk factor for NAFLD.
Conclusion: This research highlights a significant association between DINP and NAFLD. These findings underscore the need for further investigation into the role of the phthalates substitutes in the pathogenesis of NAFLD and the importance of considering different diagnostic criteria in research.
1 Introduction
Non-alcoholic fatty liver disease (NAFLD) is the most common cause of chronic liver disease and is likely to be the most common indication for liver transplantation by 2030 (1). A meta-analysis of studies published between 1990 and 2015 suggests that the global prevalence of NAFLD is about 25% (2). As of 2019, the global total prevalence was as high as 30.05%, which increased by about 50 percent between 2016–2019, affecting 37% of Americans (3, 4). NAFLD is a liver disease associated with a metabolic disorder characterized by the accumulation of excessive fat in the liver, a condition that is not caused by alcohol consumption (5). The progression of the disease includes steatohepatitis, fibrosis, and eventually, possibly even cirrhosis (6). Multiple factors play an important role in the disease: insulin resistance (7), lipid accumulation, oxidative stress (8), gene mutation (5), sedentary lifestyle and eating patterns (6). In recent years, the relationship between phthalates exposure and diseases has been concerned.
Phthalates, widely used plasticizers in PVC and various plastics, confer flexibility and durability but raise health concerns due to their endocrine-disrupting potential and environmental pervasiveness (9). In daily life, people inevitably encounter phthalates and their metabolites through personal care products, detergents, food packaging, and even some medical devices (10). Ingestion and inhalation are two important routes for phthalates to enter the body, and intake is most often assessed by measuring metabolites in urine (11). Upon entering the human body, phthalates undergo extensive metabolic processing, yielding a multitude of metabolites. We summarize the relationships between the metabolites and their parent materials, as well as their chemical formulas, for the reader’s reference (Supplementary Table S1; Supplementary Figure S1). Phthalates have been linked to a variety of diseases, such as hyperuricemia (12), reproductive health (13), cardiovascular disease risk (14, 15), metabolic syndrome (16) and liver injury (17).
The studies that have been published so far have applied different diagnostic methods to reach controversial conclusions. A study based on the database of a representative population in Korea, The Korean National Environmental Health Survey (KoNEHS), indicated that the higher concentration of MEHHP were positively correlated with the prevalence of NAFLD using HIS diagnostic criteria (18). However, the results of another study based on the US population showed that not only MEHHP, but also MEOHP and MECPP showed a significant positive association with NAFLD employing the HSI diagnostic criteria. Unfortunately, only MEHHP maintained this correlation when they were based on the US.FLI standard (19). Furthermore, a recent study found that under the diagnostic criteria of HSI and US.FLI, there is no correlation between MEHHP, MEOHP, and MECPP with NAFLD. However, based on the VCTE standard, MEHHP is positively correlated with NAFLD (20). Another similar study also employed the VCTE diagnostic criteria, and they found that, after adjusting for all variables, the above three metabolites are not associated with NAFLD (20, 21). Clearly, different researchers applying various diagnostic criteria to assess NAFLD have reached controversial conclusions. Currently, there are no studies that simultaneously apply all three criteria to the same group of subjects for diagnosis and compare how the choice of diagnostic methods affects the evaluation of disease relevance.
Due to changes in legislation, the use of DEHP has gradually decreased, and the phthalates substitutes is getting more attention. For metabolites of DEHP alternatives, MEP can elevate levels of alanine aminotransferase (ALT) and aspartate aminotransferase (AST), enzymes potentially indicative of liver disease (22). Within the European Union (EU), DEHP is recognized as harmful to reproduction and requires clear labeling with a skull-and-crossbones symbol and the word “TOXIC” when in pure form (Regulation of DEHP | Health Care Without Harm1). For drugs and biological products under the CDER’s jurisdiction, the FDA strongly discourages the use of DEHP as an excipient (FDA-2012-D-1135). As a result of these regulations, the use of DEHP is decreasing, but the use of its alternatives is gradually increasing, such as DINP (23). In addition, other DEHP alternatives are also worthy of attention, such as DIDP, DINP, DnBP, DIBP, DEP and DnOP. MnBP, a metabolite of DnBP, is believed to be correlated with NAFLD (24). However, the research on DEHP substitutes and liver diseases is still insufficient, and the research on the correlation of alternative compounds and NAFLD is even less known.
This study utilized nationally representative data from the National Health and Nutrition Examination Surveys (NHANES) database in the United States to explore the association between phthalates and NAFLD under different prediction criteria (including HSI, US.FLI and VCTE). It is recommended that healthcare professionals closely examine the association between different metabolites and disease when selecting predictive models other than the gold standard for evaluating NAFLD, to enhance the accuracy of preventive measures.
2 Materials and methods
2.1 Study design and participants
The National Health and Nutrition Examination Survey (NHANES) is a significant public health database to assess the health and nutritional status of adults and children in the United States. It’s conducted by the National Center for Health Statistics (NCHS), which is part of the Centers for Disease Control and Prevention (CDC). In the process of implementation, researchers conduct laboratory and imaging tests in specially designed and equipped Mobile Examination Centers (MECs). The use of professional collection, processing, and storage protocols ensures that the data obtained from the NHANES database is both accurate and reliable. The analysis of this data not only represents the specific population sampled, but also considering the database employs a complex, multistage probability sampling design to select participants, it is representative of the broader civilian U.S. population.
In this study, we collected data from the NHANES database for the 2017 to 2018 cycles, involving 9,254 participants. Initially, as our focus was on adults, we excluded 3,398 individuals under the age of 18, leaving 5,856 adults. Of these, approximately two-fifths (2,342 participants) underwent testing for indicators used to assess the presence of NAFLD (ALT, AST, BMI, Diabetes, GGT, Waist circumference, Fasting insulin and Fasting glucose). However, only 772 participants were tested for the presence of phthalates and plasticizer metabolites in their urine. Additionally, we excluded patients with liver cancer, autoimmune hepatitis, and viral hepatitis (hepatitis C antibody, hepatitis C RNA positive; hepatitis B surface antigen positive). Furthermore, to eliminate the influence of excessive alcohol consumption on metabolic indicators and NAFLD diagnosis, we excluded 106 participants who were heavy drinkers. Excessive alcohol consumption was characterized by an average consumption exceeding 20 grams per day for men and 10 grams per day for women (25). Ultimately, 641 participants were enrolled (Figure 1).
2.2 Measurement of phthalates in urine
In our study, we evaluated ten major metabolites derived from seven different parent phthalates, including Mono-2-ethyl-5-carboxypentyl phthalate (MECPP), Mono-(2-ethyl-5-hydroxyhexyl) phthalate (MEHHP), Mono-(2-ethyl-5-oxohexyl) phthalate (MEOHP), Mono-(carboxyisononyl) phthalate (MCiNP), Mono-oxoisononyl phthalate (MOiNP), Mono-(carboxyisoctyl) phthalate (MCiOP), Mono-n-butyl phthalate (MnBP), Mono-isobutyl phthalate (MiBP), Mono-ethyl phthalate (MEP), Mono-(3-carboxypropyl) phthalate (MCPP). A thorough summary of each phthalate and its corresponding metabolites is provided in the Supplementary Figure S1. In this study, the urine samples were first processed by enzymolysis of glucosylated analytes, and then quantitatively detected by high performance liquid chromatography-electrospray ionization-tandem mass spectrometry (HPLC-ESI-MS/MS). In addition to this, the collection of these samples is often carried out in multiple locations across the country to ensure that the samples are representative.2
For sample analysis, values not detected due to being below the lower limit of detection (LOD) are populated with the LOD divided by the square root of 2. This approach aligns with the guidelines set by the National Center for Health Statistics (NCHS). A detailed analysis guide can be found at the following website: https://wwwn.cdc.gov/nchs/nhanes/analyticguidelines.aspx#analytic-guidelines.
2.3 Definition of NAFLD
Liver biopsy, traditionally the gold standard for diagnosing NAFLD, faces limitations due to its invasiveness and poor patient acceptance, making it less ideal for regular screening. To address this, researchers have introduced two effective models for disease screening: the HSI criteria and the US.FLI criteria. In the HSI prediction model, the AUC (95%CI) was 0.812 (0.801–0.824), and the cutoff values were 30 and 36. Values below 30, considering factors like ALT, AST, BMI, and gender, imply a minimal risk of NAFLD. In contrast, values above 36 suggest a high likelihood of the condition (26). In the US.FLI prediction model, the AUC (95%CI) was 0.80 (0.77–0.83), and the cutoff values were 30. If US.FLI ≥ 30, it was defined as NAFLD (27). HSI and US.FLI were calculated as follows:
HSI = 8*ALT (IU/L)/AST (IU/L) + body mass index (BMI, kg/m2) + 2 (if female) +2 (if type 2 diabetes).
HSI = 8*ALT (IU/L)/AST (IU/L) + body mass index (BMI, kg/m2) + 2 (if female) +2 (if type 2 diabetes).
VCTE is a specific liver elastography technique that assesses the stiffness of the liver by measuring its response to brief mechanical vibrations, defined by two variables, CAP and LSM. CAP was used to assess hepatic steatosis, and LSM was used to assess hepatic fibrosis. In this study, we referred to the conclusions of Eddowes et al., and defined NAFLD as LSM ≥ 8 kPa and CAP≥274 dB/m. In the VCTE criteria, the AUC (95%CI) was 0.87 (0.82–0.92) (28).
2.4 Definition of covariates
This study selected some of the demographic data, examination data, laboratory data, and variables defined by questionnaire data, including age (>18), gender (male/female), race (Mexican American/ Other Hispanic/ Non-Hispanic White/ Non-Hispanic Black/Non-Hispanic Asian/Other Race), education level (the highest degree obtained is either above or below high school) and other indicators. Smoking status was defined based on the total number of cigarettes smoked in a lifetime and whether they had smoked recently. Recreational physical activity was divided into vigorous activity, moderate activity, sedentary activity. Body mass index (BMI) was categorized into under/normal weight (<25 kg/m2), overweight (25–30 kg/m2), obese (> = 30 kg/m2). In addition, there was another indicator of obesity, waist circumference. Ensuring data accuracy, systolic blood pressure (SBP) measurements were conducted three times, with the average of these values being utilized for analysis. Diabetes was defined as a hemoglobin A1C (HbA1C) concentration > = 6.5% or fasting plasma glucose (FPG) > = 126 mg/dL or a person who had been told by a doctor that he or she had diabetes and was currently taking insulin. Other covariates, such as alanine aminotransferase (ALT), aspartate aminotransferase (AST), gamma-glutamyl transpeptidase (GGT), lactate dehydrogenase (LDH), albumin (ALB), globulin (GLO), were also included in this study. Covariates were selected based on their association with exposure and outcomes, informed by clinical experience and supported by the available literature.
2.5 Statistical analysis
The data in NHANES uses a complex multi-stage probabilistic sampling design. Therefore, this study referred to the methodology of the National Center for Health Statistics (NCHS) to accurately evaluate large population data in the United States, considering the primary sampling units (PSU) and the stratum variable (29). Due to the smallest subpopulation composition of phthalates and their metabolites, this study incorporated the use of subsample weights (wtsb2yr). All statistical analyses were performed using R (version 4.3.2), and Beckman Coulter DxAI platform v2.0.3 Furthermore, to address the skewed distribution, we transformed the urinary phthalate metabolites data using the natural logarithm (ln) method.
Continuous variables were described as median (IQR). Categorical variables were described as n (unweighted) (%). However, the calculation process is weighted analysis. Differences between groups were evaluated using Mann–Whitney U tests for continuous variables, depending on the skewness of the data. Chi-square tests were utilized for categorical variables. Multiple logistic regression analyses were conducted to explore the associations between Phthalates and suspected NAFLD, adjusting for potential confounders. All tests were two-tailed, and a p-value of less than 0.05 was considered statistically significant. Variables included in the logistic regression analysis underwent strict screening procedures. Specifically, due to the high correlation between the two variables (FPG vs. HbA1C; BMI vs. waist circumference, Pearson’s r ≥ 0.6) (Supplementary Figure S2), only one of the indicators was selected for adjustment, and variables with high correlation to the outcome variables were excluded (Supplementary Figure S3). Finally, we identify the variables that need to be adjusted for different models: model 1 was unadjusted; model 2, age and race were adjusted; model 3, age, race, recreational physical activity, SBP and HbA1C were adjusted.
3 Results
3.1 Baseline characteristics
Based on the inclusion and exclusion criteria previously mentioned, a total of 641 participants were enrolled in the study. These participants were evaluated using three different NAFLD prediction criteria: HSI > 36, US.FLI ≥ 30, and VCTE: LSMs >8 kPa and CAP ≥274 dB/m. Unfortunately, only 572 participants underwent VCTE assessment, owing to missing data on median CAP, median stiffness, or elastograph status. Table 1 presents the baseline characteristics. Under the HSI diagnostic criteria, participants with higher levels of obesity, fasting blood glucose, HbA1C, ALT, GGT, lower albumin, and higher globulin were more likely to have NAFLD. An increase in the concentration of phthalates and plasticizers metabolites in urine, specifically MEOHP, MEHHP, MECPP, MCiOP, and MCiNP, May be associated with NAFLD. Under the US.FLI diagnostic criteria, older age was more likely to be associated with the disease (42 vs. 53, p < 0.001). Similar to the HSI group, in addition to significant statistical differences in obesity, fasting blood glucose, HbA1C, ALT, GGT, albumin, and globulin, an increase in SBP and AST (117 vs. 125, p < 0.001; 18 vs. 20, p = 0.01) May also be related to NAFLD. Furthermore, higher concentrations of MnBP, MiBP, MCPP, and MEHHP might be linked to an increased likelihood of the disease. Under the VCTE criteria, trends were similar to those in the US.FLI results, but there were no statistical differences in AST, albumin, and globulin concentrations between the two groups. MEOHP, MEHHP, MCiOP, and MCiNP May be related to NAFLD. Additionally, Table 1 indicates that irrespective of the diagnostic criteria used, factors such as gender, education level, smoking status, and recreational physical activity showed no significant statistical differences, suggesting no association with NAFLD.
3.2 Association between baseline information and NAFLD
Univariate logistic regression analysis was conducted on variables exhibiting statistical differences in demographic, measurement, and laboratory test data as presented in Table 2. The objective was to identify factors potentially associated with outcomes in the same cohort of subjects when subjected to different diagnostic criteria for NAFLD. This analysis aids in the selection of variables for subsequent multivariate logistic regression.
Table 2. Univariate logistic regression analysis for all variables except urinary phthalates levels.
As demonstrated in Table 2, heterogeneous results were obtained when the same cohort underwent diagnosis using varying predictive criteria. Under the HSI criteria, physical measurements and laboratory indicators such as BMI, waist circumference, diabetes status, HbA1C, fasting plasma glucose, ALT, and globulin levels were positively correlated with NAFLD, all showing statistical significance. In contrast, serum albumin levels were inversely related to the disease (OR = 0.28, 95%CI = 0.13–0.60, p = 0.003). Additionally, compared to Mexican Americans, other races showed a negative correlation with the disease. Under the US.FLI criteria, the older the age, the greater the positive association with the disease (OR = 1.02, 95%CI = 1.01–1.03, p < 0.001). Furthermore, GGT was positively correlated with NAFLD (OR = 1.06, 95% CI = 1.01–1.12, p = 0.034), with other results similar to those under the HSI criteria. Under the VCTE criteria, age, race, BMI, diabetes status, waist circumference, systolic blood pressure, fasting plasma glucose, HbA1C, ALT, and GGT were all positively correlated with NAFLD, as detailed in Table 2.
3.3 Association between urinary phthalates and NAFLD
The final metabolites of various phthalates May have different associations with NAFLD under different diagnostic criteria (Table 3). We observed that under the HSI criteria, there were significant positive correlations between NAFLD and the MEOHP, MEHHP, MECPP, and MCiOP (OR = 1.46, 95% CI = 1.11–1.93, p = 0.011; OR = 1.46, 95% CI = 1.10–1.95, p = 0.013; OR = 1.44, 95% CI = 1.04–2.00, p = 0.030; OR = 1.26, 95% CI = 1.01–1.58, p = 0.043). Under the US.FLI criteria, in addition to MEOHP, MEHHP, and MECPP, MnBP and MiBP also showed a positive correlation with NAFLD (OR = 1.29, 95% CI = 1.02–1.61, p = 0.033; OR = 1.26, 95% CI = 1.03–1.54, p = 0.026), while the correlation with MCiOP was not observed. Similarly to the HSI criteria, under the VCTE standards, MEOHP, MEHHP, MECPP, and MCiOP were positively correlated with the disease (OR = 1.47, 95% CI = 1.11–1.94, p = 0.010; OR = 1.45, 95% CI = 1.07–1.96, p = 0.020; OR = 1.53, 95% CI = 1.10–2.13, p = 0.014). Additionally, MCiNP also displayed a statistically significant positive correlation with the disease (OR = 1.36, 95% CI = 1.03–1.79, p = 0.035).
The above results illustrate the independent association of each indicator with NAFLD under different disease prediction criteria. However, it is yet to be ascertained whether these indicators are interrelated, and consequently, if they collectively contribute to the pathogenesis of NAFLD. Furthermore, to ascertain if the effect exhibited concentration-dependence, we classified phthalates into four distinct concentration gradients. Consequently, we conducted a multivariate logistic regression analysis to account for the combined impact of various factors on the outcome.
3.4 Correlation between phthalates and NAFLD based on logistic regression model
Three logistic regression models were constructed: model 1, no covariates were adjusted (Table 4); model 2, age and race were adjusted (Table 5); model 3, age, race, recreational physical activity, systolic blood pressure and HbA1C were adjusted (Table 6). We investigated the relationship between varying concentration gradients of urinary phthalates and NAFLD in a cohort, using different diagnostic criteria. Under the HSI criteria, elevated concentrations of MEP, MCPP, MOiNP, MEOHP, MEHHP, MECPP, MCiOP, and MCiNP exhibited positive correlations with NAFLD, suggesting an increased disease likelihood at higher levels. Specifically, most of this correlation is concentrated in the third quartile (Q3) versus the first (Q1), regardless of model 1 (M1) or model 2 (M2). Odds Ratios for Q3 vs. Q1 were as follows: MEP (M1: 1.93; M2: 1.89), MCPP (M1: 2.27; M2: 2.38), MOiNP (M1: 2.86; M2: 2.93), MEOHP (M1: 2.68; M2: 2.81), MEHHP (M1: 2.75; M2: 2.85), MECPP (M1: 4.19; M2: 4.22), with p-values <0.05. MCiOP and MCiNP also showed significant positive correlation, particularly in Q2 vs. Q1 (Tables 4, 5). Under the US.FLI standard, similar trends were observed. MEP, MnBP, MiBP, MCPP, MEOHP, MEHHP, and MECPP displayed positive correlations with NAFLD in model 1 and model 2. The higher the urinary phthalates, the higher the positive correlation with NALFD. Furthermore, MiBP and MCPP demonstrated significant associations across multiple quartiles, indicating a consistent trend. The VCTE criteria further corroborated these findings. Compounds like MEOHP, MEHHP, MECPP, and MCiOP exhibited positive correlations with NAFLD (Tables 4, 5).
Table 4. Model 1—univariate logistic regression analysis for NAFLD according to urinary phthalates quantile levels.
Table 5. Model 2—multivariate logistic regression analysis for NAFLD according to urinary phthalates levels.
Table 6. Model 3—multivariate logistic regression analysis for NAFLD according to urinary phthalates levels.
However, when all covariates were adjusted for multivariate logistic regression analysis, we found that no matter which prediction method was used for diagnosis (HSI or US.FLI) in our cohort, all the phthalates mentioned above lost their relevance. Only under the criterion of VCTE, a higher concentration of MCiOP in urine was positively correlated with NAFLD (OR = 4.22, 95% CI = 1.10–16.2, p = 0.044), and this correlation was found in models 1,2, and 3 (Table 6). If participants used VCTE to predict disease outcomes, those with more MCiOP in their urine were more likely to develop the disease, suggesting that MCiOP May be a risk factor for NAFLD, although this was cross-sectional study.
4 Discussion
This study, utilizing NHANES database data, reveals that higher concentrations of specific phthalate metabolites correlate positively with NAFLD, with the correlation’s strength varying by diagnostic criteria. Notably, only under VCTE criteria did the Mono-(carboxyisoctyl) phthalate (MCiOP) correlate positively with NAFLD. Based on what we know, this research innovatively compares multiple NAFLD diagnostic methods (HSI>36, US.FLI ≥ 30 and VCTE:LSM ≥ 8 kPa or CAP≥274 dB/m), highlighting the impact of diagnostic method selection on correlations. Furthermore, our analysis extends beyond traditional phthalates and their metabolites, incorporating a broader range of substitutes to provide a more comprehensive understanding of how various metabolites influence disease.
Phthalates May affect NAFLD through multiple pathways. Cell experiments show that phthalates exposure can trigger lipid accumulation and amplify oxidative stress in hepatocytes via the suppression of the JAK2/STAT5 pathway in BRL-3A cells, raising the risk of NAFLD (30). In animal studies, phthalates have been shown to facilitate macrophage polarization by simultaneously activating PPARα and PPARγ. This coordinated activation contributes to the disruption of lipid homeostasis within the liver (31). Moreover, phthalates can also cause liver cell damage by affecting oxidative stress metabolism (32). In population-based studies, researchers utilizing data from The Korean National Environmental Health Survey (KoNEHS) indicated that the third and fourth quartiles of MEHHP were positively correlated with the prevalence of NAFLD among Korean adults (OR = 1.33, 95% CI = 1.00–1.78; OR = 1.39, 95% CI = 1.00–1.92), employing the HSI diagnostic criteria (18). Another analysis, using data from the NHANES 2003–2016 database of a representative US population, based on the HSI diagnostic criteria, showing that MEOHP (OR = 1.56, 95% CI = 1.08–2.24), MEHHP (OR = 1.55, 95% CI = 1.09–2.21), and MECPP (OR = 1.44, 95% CI = 1.06–1.95) are associated with NAFLD. However, only MEHHP (OR = 1.98, 95% CI = 1.32–2.97) was found to be related to NAFLD under the US.FLI diagnostic method (19). Recently, another study included individuals from the NHANES 2017 to 2018 cycles and discovered a significant correlation between MECPP (OR = 2.719, 95% CI = 1.296–5.700) and MEHHP (OR = 2.073, 95% CI = 1.111–3.867) with NAFLD under the VCTE criteria. Nonetheless, this correlation vanished under HSI and US.FLI standards. Similarly, He et al. found that patients with higher urinary DEHP metabolite concentrations had an increased risk of NAFLD in the partially adjusted logistic regression model (OR = 1.22, 95%CI = 1.09–1.36), but lost the association in the fully adjusted model based on the VCTE (20, 21). These contradictory findings indicate that the relationship between phthalate metabolites and NAFLD is influenced by two factors: the selection of the population and the choice of diagnostic methods. In our research, fully adjusted multivariate logistic regression results demonstrated that three metabolites of DEHP (MECPP, MEHHP, and MEOHP) were not associated with NAFLD, similar to the above conclusions. Due to legislative impacts, exposure to DEHP has been decreasing over the years, with its alternatives, such as DINP—a less toxic substitute—increasingly being used (33). Utilizing three predictive methods, HSI, US.FLI, and VCTE, our analysis of various phthalates and their substitutes’ metabolites found that only when the VCTE standard was used for disease diagnosis, one of DINP’s metabolites, MCiOP, could potentially have an adverse effect on the occurrence of NAFLD (OR = 4.22, 95% CI = 1.10–16.2, p = 0.044). The reason for these differences could be partly due to the HSI and US.FLI diagnostic criteria relying on serum laboratory markers, which are influenced by various factors and are not as sensitive as VCTE. VCTE can directly measure liver stiffness, which more accurately reflects the degree of fibrosis associated with NAFLD (34). Another factor could be the differences in inclusion and exclusion criteria between studies, as well as the covariates considered, leading to different logistic regression models being constructed.
DINP emerges as a leading alternative in industrial applications to DEHP (35). Humans are widely exposed to DINP and May affect NAFLD through multiple pathways (36, 37). At present, the understanding of how DINP influences NAFLD is limited. Nonetheless, potential pathways through which DINP May contribute to NAFLD include alterations in lipid metabolism, oxidative stress, inflammation, and cellular communication mechanisms. Firstly, DINP promotes adipocyte differentiation and lipid accumulation in 3 T3-L1 pre-adipocytes through the activation of PPARγ (Peroxisome Proliferator-Activated Receptor Gamma), highlighting its role in adipogenesis and lipid metabolism (38). High doses of DINP can cause lipid metabolism disorders in male mice (39), suggesting that it May be conducive to the formation of fatty liver. Moreover, DINP modulates lipid metabolism by influencing the Endocannabinoid System (ECS), resulting in heightened accumulation of lipids and triglycerides in hepatic tissue (40). Secondly, DINP acts as a toxic agent that causes damage to liver tissues, induces oxidative stress by increasing ROS levels and depleting glutathione (GSH), promotes lipoperoxidation and DNA damage, and elevates pro-inflammatory cytokine levels (41). Moreover, DIDP can also lead to increased serum concentrations of ALT and AST among Balb/c mice, illustrating its potential to induce liver damage (42). Thirdly, the perturbation in ROS homeostasis, induced by DINP, catalyzes the activation of inflammatory responses, precipitating the secretion of inflammatory mediators. Studies have demonstrated that DINP induces an upregulation of IL-1 and TNF-α in Kunming mice (41). This phenomenon is ostensibly linked to the activation of the NF-κB signaling pathway, suggesting a mechanistic pathway through which DINP exerts its pro-inflammatory effects (43). Additionally, dermal exposure DINP enhances pulmonary inflammation via activation of the IL-31/TRPV1 pathway, leading to increased levels of IL-4, IL-5, IL-6, and IL-13, and decreased IFN-gamma in mice (44). DINP can also promote airway inflammation by activating the PI3K/AKT/NF-κB signaling pathway (45). Whether this mechanism exists in the liver remains to be further studied. Finally, DINP activates the MAPK-Erk1/2 pathway and, upon extended exposure, can inhibit gap junctional intercellular communication (GJIC), associated with the progression of fatty liver disease, fibrosis, and cirrhosis (46, 47). Regrettably, there is still a lack of research on human liver cells or tissues, and the mechanisms mentioned May exhibit interspecies variations, requiring further verification. However, it can still be stated that DINP May become a significant underlying factor influencing NAFLD through pathways such as lipid metabolism, oxidative stress, inflammatory responses, and cellular communication.
This study’s limitations are noteworthy. Initially, its cross-sectional approach, based on single-point phthalate level measurements from NHANES data, fails to capture the fluctuations in phthalate exposure. Phthalate metabolites have a short half-life in the body. This means that a single urine sample May only reflect short-term exposure (from a few hours to a few days) and May not accurately represent long-term or cumulative exposure. While a single measurement is still valuable for providing insights into recent exposure and possible associations with health outcomes. Future research should consider employing repeated sampling to obtain a more comprehensive understanding of long-term exposure. Subsequently, in exploring the association between multiple phthalates and NAFLD, this study leverages diverse diagnostic criteria, including HSI, US.FLI, and VCTE, through rigorous participant selection criteria. However, the availability of VCTE data is currently restricted to the 2017–2018 NHANES cycle. The anticipated update of the NHANES database will provide an opportunity for statistical analyses encompassing a broader dataset. Furthermore, although these three diagnostic methods can assess NAFLD, it is important to note that NAFLD progresses through multiple stages, such as Non-Alcoholic Fatty Liver (NAFL), Non-Alcoholic Steatohepatitis (NASH), fibrosis, and cirrhosis. Unfortunately, the current database does not provide sufficient data to allow for a stratified analysis of these stages. We acknowledge this limitation and recognize the need for future research to conduct more detailed studies on the staging of NAFLD to further investigate the association between phthalates and different stages of the disease. In addition, the covariates included in this study, such as smoking status, recreational physical activity, and insulin injection status, were all obtained through questionnaires, which May introduce recall bias. However, it is worth noting that the questionnaire data included in the NHANES database were collected by trained interviewers using standardized questionnaires and interview procedures, ensuring the scientific rigor and reliability of the data. Lastly, HSI, US.FLI, and VCTE are important tools in NAFLD research but are not considered gold standards like liver biopsy, which, despite its accuracy, is invasive and unsuitable for large-scale screenings. These non-invasive methods face limitations in sensitivity and specificity, potentially exaggerating phthalates’ role in NAFLD. However, these methods are acknowledged and have demonstrated utility in clinical diagnosis (26–28). Our study not only confirms a correlation between phthalates and NAFLD but also reveals that the choice of predictive criteria affects this correlation. VCTE is identified as a more accurate method than HSI and US.FLI for assessing liver stiffness (34), and its use has become prevalent for evaluating NAFLD. Our data suggest a notable link between DINP and NAFLD when diagnoses are made using VCTE criteria, which could improve early NAFLD detection and recommend that patients limit their exposure to DINP-containing plastic products to protect their liver health.
Numerous studies have investigated the association between phthalates and NAFLD. In our research, we consolidate existing viewpoints and uncovered the novel influence of different diagnostic methods on study outcomes. Additionally, we reveal that the correlation between DINP and NAFLD is significant, particularly when clinicians use non-invasive methods like the VCTE criteria for disease prediction. This awareness can prompt clinicians to consider environmental and chemical exposures as part of their assessment for NAFLD risk factors. Our study suggests that when utilizing VCTE criteria to assess NAFLD, the content of DINP secondary metabolites in urine May serve as a more potent predictor. It suggests that limiting exposure to DINP-containing products could be a viable preventative measure for at-risk populations. Furthermore, these insights advocate for incorporating routine environmental health screenings into clinical practice, enabling early intervention and personalized patient education. Future research should continue exploring the mechanistic pathways of DINP and validate these findings across diverse populations and longer-term studies.
5 Conclusion
This study establishes a positive link between phthalate exposure, specifically DINP (the parent body of the MCOP), and NAFLD risk in adults, highlighting the influence of diagnostic criteria on this relationship. Given DINP’s potential liver damage risk, it is imperative to conduct future case–control or cohort studies to confirm a causal relationship. Additionally, combining two or more 2-year cycles of the continuous NHANES is beneficial to understand this association better.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.
Author contributions
JZ: Conceptualization, Data curation, Resources, Software, Writing – original draft. QG: Investigation, Methodology, Resources, Visualization, Writing – original draft. DG: Funding acquisition, Supervision, Validation, Writing – review & editing.
Funding
The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was supported by the National Key Research and Development Program of China (2022YFC2302700); National Health Commission Hospital Management Research Institute Project (JCWJ2024ZD01); Clinical Research Special Project of the Medical and Health Technology Development Research Center of the National Health Commission (WKZX2024DN0181); Guangdong Science and Technology Foundation (2023B0101200003, 2022B1111020001, 2021A1515220084, 2020JD076); Shenzhen Science and Technology Foundation (KJZD20230923115359001, JSGG20220301090003004, ZDSYS20210623092001003, JCYJ20200109120205924); Medical-Engineering Interdisciplinary Research Foundation of ShenZhen University (2023YG009); Innovation Project of Guangxi Graduate Education (JGY2023168); and Shenzhen Portion of Shenzhen-Hong Kong Science and Technology Innovation Cooperation Zone (HTHZQSWS- KCCYB-2023060).
Acknowledgments
We extend our deepest gratitude to the National Health and Nutrition Examination Survey (NHANES) for providing access to their invaluable datasets, which made this research possible. Special thanks to our colleagues and the technical staff who offered their expertise and support throughout the study. We also appreciate the constructive feedback from peer reviewers, which significantly enhanced the quality of our manuscript. Lastly, we acknowledge the financial support, which was crucial for the completion of this project.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2024.1407976/full#supplementary-material
Footnotes
References
1. Byrne, CD, and Targher, G. NAFLD: a multisystem disease. J Hepatol. (2015) 62:S47–64. doi: 10.1016/j.jhep.2014.12.012
2. Younossi, ZM, Koenig, AB, Abdelatif, D, Fazel, Y, Henry, L, and Wymer, M. Global epidemiology of nonalcoholic fatty liver disease-Meta-analytic assessment of prevalence, incidence, and outcomes. Hepatology. (2016) 64:73–84. doi: 10.1002/hep.28431
3. Younossi, ZM, Golabi, P, Paik, JM, Henry, A, Van Dongen, C, and Henry, L. The global epidemiology of nonalcoholic fatty liver disease (NAFLD) and nonalcoholic steatohepatitis (NASH): a systematic review. Hepatology. (2023) 77:1335–47. doi: 10.1097/HEP.0000000000000004
4. Persons, PA, Bersoux, S, and Whited, MH. Is your patient at risk for NAFLD? Cleve Clin J Med. (2022) 89:719–27. doi: 10.3949/ccjm.89a.22005
5. Cobbina, E, and Akhlaghi, F. Non-alcoholic fatty liver disease (NAFLD)—pathogenesis, classification, and effect on drug metabolizing enzymes and transporters. Drug Metab Rev. (2017) 49:197–211. doi: 10.1080/03602532.2017.1293683
6. Kanuri, G, and Bergheim, I. In vitro and in vivo models of non-alcoholic fatty liver disease (NAFLD). Int J Mol Sci. (2013) 14:11963–80. doi: 10.3390/ijms140611963
7. Farzanegi, P, Dana, A, Ebrahimpoor, Z, Asadi, M, and Azarbayjani, MA. Mechanisms of beneficial effects of exercise training on non-alcoholic fatty liver disease (NAFLD): roles of oxidative stress and inflammation. Eur J Sport Sci. (2019) 19:994–1003. doi: 10.1080/17461391.2019.1571114
8. Guo, X, Yin, X, Liu, Z, and Wang, J. Non-alcoholic fatty liver disease (NAFLD) pathogenesis and natural products for prevention and treatment. Int J Mol Sci. (2022) 23:15489. doi: 10.3390/ijms232415489
9. Koch, HM, and Calafat, AM. Human body burdens of chemicals used in plastic manufacture. Philos Trans R Soc Lond Ser B Biol Sci. (2009) 364:2063–78. doi: 10.1098/rstb.2008.0208
10. Muscogiuri, G, Barrea, L, Laudisio, D, Savastano, S, and Colao, A. Obesogenic endocrine disruptors and obesity: myths and truths. Arch Toxicol. (2017) 91:3469–75. doi: 10.1007/s00204-017-2071-1
11. Council NR. Phthalates and cumulative risk assessment: The tasks ahead. Washington, DC: The National Academies Press (2008). 208 p isbn:978-0-309-12841-4 Available at: https://nap.nationalacademies.org/catalog/12528/phthalates-and-cumulative-risk-assessment-the-tasks-ahead.
12. Tan, Y, Fu, Y, Yao, H, Wu, X, Yang, Z, Zeng, H, et al. Relationship between phthalates exposures and hyperuricemia in U.S. general population, a multi-cycle study of NHANES 2007-2016. Sci Total Environ. (2023) 859:160208. doi: 10.1016/j.scitotenv.2022.160208
13. Weng, X, Zhu, Q, Liao, C, and Jiang, G. Cumulative exposure to phthalates and their alternatives and associated female reproductive health: body burdens, adverse outcomes, and underlying mechanisms. Environ Sci Technol. (2023) 57:8189–212. doi: 10.1021/acs.est.3c00823
14. Mariana, M, Castelo-Branco, M, Soares, AM, and Cairrao, E. Phthalates’ exposure leads to an increasing concern on cardiovascular health. J Hazard Mater. (2023) 457:131680. doi: 10.1016/j.jhazmat.2023.131680
15. Mariana, M, Feiteiro, J, Verde, I, and Cairrao, E. The effects of phthalates in the cardiovascular and reproductive systems: a review. Environ Int. (2016) 94:758–76. doi: 10.1016/j.envint.2016.07.004
16. Merida, DM, Moreno-Franco, B, Marques, M, Leon-Latre, M, Laclaustra, M, and Guallar-Castillon, P. Phthalate exposure and the metabolic syndrome: a systematic review and meta-analysis. Environ Pollut. (2023) 333:121957. doi: 10.1016/j.envpol.2023.121957
17. Yu, L, Yang, M, Cheng, M, Fan, L, Wang, X, Xu, T, et al. Associations between urinary phthalate metabolite concentrations and markers of liver injury in the US adult population. Environ Int. (2021) 155:106608. doi: 10.1016/j.envint.2021.106608
18. Yang, YJ, Kim, T, and Hong, YP. Urinary phthalate levels associated with the risk of nonalcoholic fatty liver disease in adults: the Korean National Environmental Health Survey (KoNEHS) 2012-2014. Int J Environ Res Public Health. (2021) 18:6035. doi: 10.3390/ijerph18116035
19. Cai, S, Fan, J, Ye, J, Rao, X, and Li, Y. Phthalates exposure is associated with non-alcoholic fatty liver disease among US adults. Ecotoxicol Environ Saf. (2021) 224:112665. doi: 10.1016/j.ecoenv.2021.112665
20. Chen, X, Tian, F, Wu, J, Liu, L, Li, Y, Yu, G, et al. Associations of phthalates with NAFLD and liver fibrosis: a nationally representative cross-sectional study from NHANES 2017 to 2018. Front Nutr. (2022) 9:1059675. doi: 10.3389/fnut.2022.1059675
21. He, Y, Zou, J, Hong, T, and Feng, D. Association between Di-2-ethylhexyl phthalate and nonalcoholic fatty liver disease among US adults: mediation analysis of body mass index and waist circumference in the NHANES. Food Chem Toxicol. (2023) 179:113968. doi: 10.1016/j.fct.2023.113968
22. Milosevic, N, Milanovic, M, Sudji, J, Bosic Zivanovic, D, Stojanoski, S, Vukovic, B, et al. Could phthalates exposure contribute to the development of metabolic syndrome and liver disease in humans? Environ Sci Pollut Res Int. (2020) 27:772–84. doi: 10.1007/s11356-019-06831-2
23. Zota, AR, Calafat, AM, and Woodruff, TJ. Temporal trends in phthalate exposures: findings from the National Health and nutrition examination survey, 2001-2010. Environ Health Perspect. (2014) 122:235–41. doi: 10.1289/ehp.1306681
24. Yang, EJ, Choi, BS, and Yang, YJ. Risk of nonalcoholic fatty liver disease is associated with urinary phthalate metabolites levels in adults with subclinical hypothyroidism: Korean National Environmental Health Survey (KoNEHS) 2012-2014. Int J Environ Res Public Health. (2022) 19:3267. doi: 10.3390/ijerph19063267
25. Chalasani, N, Younossi, Z, Lavine, JE, Charlton, M, Cusi, K, Rinella, M, et al. The diagnosis and management of nonalcoholic fatty liver disease: practice guidance from the American Association for the Study of Liver Diseases. Hepatology. (2018) 67:328–57. doi: 10.1002/hep.29367
26. Lee, JH, Kim, D, Kim, HJ, Lee, CH, Yang, JI, Kim, W, et al. Hepatic steatosis index: a simple screening tool reflecting nonalcoholic fatty liver disease. Dig Liver Dis. (2010) 42:503–8. doi: 10.1016/j.dld.2009.08.002
27. Ruhl, CE, and Everhart, JE. Fatty liver indices in the multiethnic United States National Health and nutrition examination survey. Aliment Pharmacol Ther. (2015) 41:65–76. doi: 10.1111/apt.13012
28. Eddowes, PJ, Sasso, M, Allison, M, Tsochatzis, E, Anstee, QM, Sheridan, D, et al. Accuracy of fibro scan controlled attenuation parameter and liver stiffness measurement in assessing steatosis and fibrosis in patients with nonalcoholic fatty liver disease. Gastroenterology. (2019) 156:1717–30. doi: 10.1053/j.gastro.2019.01.042
29. Chen, TC, Clark, J, Riddles, MK, Mohadjer, LK, and Fakhouri, THI. National Health and nutrition examination survey, 2015-2018: sample design and estimation procedures. Vital Health Stat. (2020) 2:1–35.
30. Zhang, Y, Wang, S, Zhao, T, Yang, L, Guo, S, Shi, Y, et al. Mono-2-ethylhexyl phthalate (MEHP) promoted lipid accumulation via JAK2/STAT5 and aggravated oxidative stress in BRL-3A cells. Ecotoxicol Environ Saf. (2019) 184:109611. doi: 10.1016/j.ecoenv.2019.109611
31. Xu, M, Li, Y, Wang, X, Zhang, Q, Wang, L, Zhang, X, et al. Role of hepatocyte-and macrophage-specific PPARgamma in hepatotoxicity induced by Diethylhexyl phthalate in mice. Environ Health Perspect. (2022) 130:17005. doi: 10.1289/EHP9373
32. Aydemir, D, Aydogan-Ahbab, M, Barlas, N, and Ulusu, NN. Effects of the in-utero dicyclohexyl phthalate and di-n-hexyl phthalate administration on the oxidative stress-induced histopathological changes in the rat liver tissue correlated with serum biochemistry and hematological parameters. Front Endocrinol (Lausanne). (2023) 14:1128202. doi: 10.3389/fendo.2023.1128202
33. Yang, S, Arcanjo, RB, and Nowak, RA. The effects of the phthalate DiNP on reproductiondagger. Biol Reprod. (2021) 104:305–16. doi: 10.1093/biolre/ioaa201
34. Nogami, A, Yoneda, M, Iwaki, M, Kobayashi, T, Kessoku, T, Honda, Y, et al. Diagnostic comparison of vibration-controlled transient elastography and MRI techniques in overweight and obese patients with NAFLD. Sci Rep. (2022) 12:21925. doi: 10.1038/s41598-022-25843-6
35. Wang, Y, Zhu, H, and Kannan, K. A review of biomonitoring of phthalate exposures. Toxics. (2019) 7:21. doi: 10.3390/toxics7020021
36. Silva, MJ, Reidy, JA, Preau, JL Jr, Needham, LL, and Calafat, AM. Oxidative metabolites of diisononyl phthalate as biomarkers for human exposure assessment. Environ Health Perspect. (2006) 114:1158–61. doi: 10.1289/ehp.8865
37. Wittassek, M, and Angerer, J. Phthalates: metabolism and exposure. Int J Androl. (2008) 31:131–8. doi: 10.1111/j.1365-2605.2007.00837.x
38. Zhang, L, Sun, W, Duan, X, Duan, Y, and Sun, H. Promoting differentiation and lipid metabolism are the primary effects for DINP exposure on 3T3-L1 preadipocytes. Environ Pollut. (2019) 255:113154. doi: 10.1016/j.envpol.2019.113154
39. Huang, Y, Sun, F, Tan, H, Deng, Y, Sun, Z, Chen, H, et al. DEHP and DINP induce tissue-and gender-specific disturbances in fatty acid and Lipidomic profiles in neonatal mice: a comparative study. Environ Sci Technol. (2019) 53:12812–22. doi: 10.1021/acs.est.9b04369
40. Forner-Piquer, I, Mylonas, CC, Calduch-Giner, J, Maradonna, F, Gioacchini, G, Allara, M, et al. Endocrine disruptors in the diet of male Sparus aurata: modulation of the endocannabinoid system at the hepatic and central level by Di-isononyl phthalate and Bisphenol a. Environ Int. (2018) 119:54–65. doi: 10.1016/j.envint.2018.06.011
41. Ma, P, Yan, B, Zeng, Q, Liu, X, Wu, Y, Jiao, M, et al. Oral exposure of Kunming mice to diisononyl phthalate induces hepatic and renal tissue injury through the accumulation of ROS. Protective effect of melatonin. Food Chem Toxicol. (2014) 68:247–56. doi: 10.1016/j.fct.2014.03.027
42. Chen, Y, Li, C, Song, P, Yan, B, Yang, X, Wu, Y, et al. Hepatic and renal tissue damage in Balb/c mice exposed to diisodecyl phthalate: the role of oxidative stress pathways. Food Chem Toxicol. (2019) 132:110600. doi: 10.1016/j.fct.2019.110600
43. Liu, J, Han, X, Zhang, T, Tian, K, Li, Z, and Luo, F. Reactive oxygen species (ROS) scavenging biomaterials for anti-inflammatory diseases: from mechanism to therapy. J Hematol Oncol. (2023) 16:116. doi: 10.1186/s13045-023-01512-7
44. Peng, Q, Wu, Y, Li, Y, Lu, C, Yao, R, Hu, S, et al. The IL-31/TRPV1 pathway mediates allergic asthma exacerbated by DINP dermal exposure in OVA-sensitized Balb/c mice. Sci Total Environ. (2024) 912:169613. doi: 10.1016/j.scitotenv.2023.169613
45. Chen, L, Chen, J, Xie, CM, Zhao, Y, Wang, X, and Zhang, YH. Maternal Disononyl phthalate exposure activates allergic airway inflammation via stimulating the phosphoinositide 3-kinase/Akt pathway in rat pups. Biomed Environ Sci. (2015) 28:190–8. doi: 10.3967/bes2015.025
46. Hernandez-Guerra, M, Hadjihambi, A, and Jalan, R. Gap junctions in liver disease: implications for pathogenesis and therapy. J Hepatol. (2019) 70:759–72. doi: 10.1016/j.jhep.2018.12.023
Keywords: NAFLD, phthalates, National Health and Nutrition Examination Survey (NHANES), MCiOP, DINP
Citation: Zou J, Gu Q and Gu D (2024) Association between phthalates exposure and non-alcoholic fatty liver disease under different diagnostic criteria: a cross-sectional study based on NHANES 2017 to 2018. Front. Public Health. 12:1407976. doi: 10.3389/fpubh.2024.1407976
Edited by:
Yanqiang Li, Boston Children’s Hospital and Harvard Medical School, United StatesReviewed by:
Yi Zheng, Merck & Co., Inc., United StatesPingping Xiang, Icahn School of Medicine at Mount Sinai, United States
Copyright © 2024 Zou, Gu and Gu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Dayong Gu, d2FuaG9vZEAxNjMuY29t