Study population
The 2018 METAL study was performed to investigate the association between exposure to environmental pollutants and diabetic complications in Chinese type 2 diabetic adults. Study participants were recruited from the outpatient clinics of seven communities in Huangpu District and Pudong District, Shanghai, China. Using a simple random sampling method, we invited half of the diabetic patients from the registration platform in each community healthcare clinic. Chinese diabetic citizens ≥18 years old who had lived in their current area for ≥6 months were included. Those with severe communication problems, acute illness, and an unwillingness
to participate were excluded (n=96).
In one of the seven clinics, all 698 registered diabetic individuals who provided urine and urinary phthalates were included. We excluded participants who were missing ultrasound results, urinary creatinine measurements and information on CVD history (n=23). A total of 675 diabetic participants were involved in the final analyses.
The study protocol was approved by the Ethics Committee of Shanghai Ninth People’s Hospital, Shanghai Jiao Tong University School of Medicine. The study protocol conformed to the ethical guidelines of the 1975 Declaration of Helsinki as reflected in a priori approval by the appropriate institutional review committee. Informed consent was obtained from all participants included in the study.
Measurements
A questionnaire about sociodemographic characteristics, medical history, family history, and lifestyle factors was adopted during the interview. The same group of trained and experienced personnel in the SPECT-China study (Wang et al., 2016; Wang et al., 2017) conducted the interviews and clinical examinations, including measurements of weight, height and blood pressure, according to a standard protocol. Body mass index (BMI) was calculated as weight in kilograms divided by height squared in meters squared. Current smoking was defined as having smoked at least 100 cigarettes in one’s lifetime and currently smoking cigarettes (Xu et al., 2013). Insulin resistance was estimated by the homeostasis model assessment index of insulin resistance (HOMA-IR): (fasting insulin [mIU/L])*(FPG [mmol/L])/22.5.
Blood samples were obtained between 6:00 am and 9:00 am after fasting for at least 8 h. Blood was refrigerated immediately after phlebotomy, and after two hours, it was centrifuged, and the serum was aliquoted and frozen in a central laboratory. Serum C-peptide was detected by the chemiluminescence method (Abbott i2000 SR, USA). Glycated hemoglobin (HbA1c) was measured by high-performance liquid chromatography (MQ-2000PT, Medconn, Shanghai, China). Fasting plasma glucose and lipid profiles were performed with a Beckman Coulter AU 680 (Brea, USA).
Carotid atherosclerosis in the common, internal, and bifurcation sites of the bilateral common carotid arteries (CCA) was assessed by the same batch of trained sonographers who were blinded to any clinical conditions of the participants with a Mindray M7 ultrasound system (MINDRAY, Shenzhen, China) with a 10-MHz probe (Wang et al., 2019). They were trained by performing a carotid ultrasound on the same patients before the study began to achieve an interobserver coefficient of variation of less than 10%. As in our previous paper (Wang et al., 2019), the measurement method was based on the consensus statement from the American Society of Echocardiography Carotid Intima-Media Thickness Task Force (Stein et al., 2008). The bilateral mean value of the carotid intima-media thickness (CIMT) was used for analysis. The CCA diameter was measured between the leading edge of the adventitia-media echo of the near wall and the leading edge of the media-adventitia echo of the far wall, based on an average of the end-diastolic diameter measurements 5-10 mm proximal to the carotid bulb.
Measurement of urinary metabolites of phthalates
We used ultra-performance liquid chromatography tandem mass spectrometry (UPLC-MS/MS) to determine 10 urinary phthalate metabolites, including monomethyl phthalate (MMP), monoethyl phthalate (MEP), mono-n-butyl phthalate (MnBP), monoisobutyl phthalate (MiBP), monobenzyl phthalate (MBzP), mono-2-ethylhexyl phthalate (MEHP), mono-2-ethyl-5-oxohexyl phthalate (MEOHP), mono-2-ethyl-5-hydroxyhexyl phthalate (MEHHP), mono-2-ethyl-5-carboxypentyl phthalate (MECPP) and mono-2-carboxymethyl-hexyl phthalate (MCMHP). Among them, MMP, MEP, MnBP, MiBP and MBzP are the metabolites of dimethyl phthalate (DMP), diethyl phthalate (DEP), di-n-butyl phthalate (DnBP), di-iso-butyl phthalate (DiBP) and benzyl butyl phthalate (BBP), respectively, and MEHP, MEOHP, MEHHP, MECPP, and MCMHP are all metabolites of dis (2-ethylhexyl) phthalate (DEHP). We also assessed the coexposure by calculating the micromolar sum of DEHP metabolites (ΣDEHP) and calculated the percentage of total ΣDEHP excreted as MEHP, referred to as %MEHP [%MEHP=MEHP /ΣDEHP] (Joensen et al., 2012).
The method has been described in our previous study with a slight modification (Dong et al., 2018). Briefly, 1 mL of the urine sample was thawed, transferred to a 10 mL glass tube, and incubated with β-glucuronidase (E. coli K 12; Roche, Mannheim, Germany) at 37 °C for 120 min. The sample was subsequently mixed with 1 mL of aqueous 2% (v/v) acetic acid and 100 μL of internal standard (100 μg/L). The mixture was loaded into a PLS column (Dikma, China; 60 mg/3 mL) previously activated with 2 mL of methanol and 2 mL of aqueous 0.5% (v/v) acetic acid. After sample loading, the column was washed with 2 mL of aqueous 0.5% (v/v) acetic acid. Next, 1 mL of methanol was added to elute the metabolites. Finally, the eluate was passed through a 0.2-μm filter and analyzed (2 μL) by an UPLC-MS/MS system integrated by Waters ACQUITY UPLC H-Class (Waters, USA) coupled with an ABSCIEX QTRAP 6500 (AB Sciex Technologies, Framingham, MA, USA). The analytical column was a Waters ACQUITY UPLC BEH C18 Column (1.7 µm, 2.1x50 mm, Waters, USA).
An internal standard method was used to quantify the target metabolite. For every 20 samples, a procedural blank and two matrix-spiked samples at two different spiking concentrations (5 and 15 ng/mL) were processed. The average recoveries and relative standard deviations (RSDs) of the target metabolites respectively ranged from 86.5% to 123.1% and from 0.5% to 12.5% at 5 μg/L and from 75.3% to 98.4% and from 0.7% to 8.8% at 15 μg/L. Sample concentrations of these metabolites were determined after subtraction of blank values. The limits of detection (LODs) were calculated at a signal-to-noise (S/N) ratio of 3 at concentrations of 0.100, 0.080, 0.006, 0.006, 0.060, 0.006, 0.010, 0.020, 0.020 and 0.040 μg/L of MMP, MEP, MnBP, MiBP, MBzP, MEHP, MEOHP, MEHHP, MECPP and MCMHP, respectively (Supplemental Table 1).
Definitions
Hypertension was assessed by systolic blood pressure ≥ 140 mmHg, diastolic blood pressure ≥90 mmHg, or a self-reported prior diagnosis of hypertension by a physician. Dyslipidemia was defined as total cholesterol ≥6.22 mmol/L (240 mg/dL), triglycerides ≥2.26 mmol/L (200 mg/dL), LDL-C≥ 4.14 mmol/L (160 mg/dL), HDL-C< 1.04 mmol/L (40 mg/dL), or a self-reported previous diagnosis of hyperlipidemia by a physician, according to the modified National Cholesterol Education Program-Adult Treatment Panel III.
The outcome CVD was defined as a self-reported diagnosis by a physician and included coronary heart disease, myocardial infarction or stroke. The related question in the questionnaire was “Have you ever been told by a doctor or other healthcare professionals that you have coronary heart disease, myocardial infarction or stroke?” The same question was adopted by another large study in Chinese, where the validation rate reached 91.07% (Lu et al., 2014). Then, the self-reported diagnoses were further verified in the registration platform. Present CCA plaque was identified as focal thickening (≥ 1.5 mm) of the artery wall.
Statistical Analysis
Data analyses were performed using IBM SPSS Statistics, Version 22 (IBM Corporation, Armonk, NY, USA). A P value <0.05 indicated significance (two sided) unless other values were mentioned. Continuous variables were summarized as medians (interquartile range) and categorical variables as percentages (%). Phthalate metabolite concentrations and urine creatinine-adjusted concentrations were presented as geometric means, percentiles and medians.
First, we performed multiple linear (continuous variable outcome) or logistic (categorical variable outcome) regression models to analyze the associations of urine creatinine-adjusted phthalate metabolites with CVD and vascular measurement. The model was adjusted for sex, age, duration of diabetes, BMI, smoking status, hypertension and dyslipidemia. The concentrations of phthalate metabolites were logarithmically transformed to achieve a normal distribution in the analyses. Backward stepwise regression was used, and the criteria for removal were P > 0.1. Bonferroni correction was used to reduce the false discovery rate with multiple comparisons.
Then, we tested the mediation and moderation effects by the SPSS PROCESS macro in an approach with 5000 bootstrap samples (Hayes, 2018). The mediation and moderation models were also adjusted for sex, age, duration of diabetes, BMI, smoking status, hypertension and dyslipidemia if necessary. First, mediation analysis was used to clarify whether exposure X was proposed as influencing outcome Y via an intervening variable M (Figure 1A). In this study, we predicted that “phthalates” impacted “CVD outcome” with “relevant CVD risk factors” as mediator variables. PROCESS was operated using one independent variable (each of the phthalate metabolites), one mediator [HOMA-IR, BMI, systolic blood pressure (SBP), low density lipoprotein (LDL), HbA1c, etc.], and one dependent variable (CVD). Second, in the moderation or interaction analysis, we predicted that exposure to different levels of CVD risk factors and medication usage would moderate the relationship between phthalate metabolites and CVD. PROCESS was operated using one independent variable (each of the phthalate metabolites), one moderator (each CVD risk factor and antidiabetic usage), and one dependent variable (CVD) (Figure 1B).