Abstract

Inorganic arsenic exposure is ubiquitous, and both exposure and interindividual differences in its metabolism have been associated with cardiometabolic risk. However, the associations of arsenic exposure and arsenic metabolism with the metabolic syndrome (MetS) and its individual components are relatively unknown. We used Poisson regression with robust variance to evaluate the associations of baseline arsenic exposure (urinary arsenic levels) and metabolism (relative percentage of arsenic species over their sum) with incident MetS and its individual components (elevated waist circumference, elevated triglycerides, reduced high-density lipoprotein cholesterol, hypertension, and elevated fasting plasma glucose) in 1,047 participants from the Strong Heart Family Study, a prospective family-based cohort study in American Indian communities (baseline visits were held in 1998–1999 and 2001–2003, follow-up visits in 2001–2003 and 2006–2009). Over the course of follow-up, 32% of participants developed MetS. An interquartile-range increase in arsenic exposure was associated with a 1.19-fold (95% confidence interval: 1.01, 1.41) greater risk of elevated fasting plasma glucose concentration but not with other individual components of the MetS or MetS overall. Arsenic metabolism, specifically lower percentage of monomethylarsonic acid and higher percentage of dimethylarsinic acid, was associated with higher risk of overall MetS and elevated waist circumference but not with any other MetS component. These findings support the hypothesis that there are contrasting and independent associations of arsenic exposure and arsenic metabolism with metabolic outcomes which may contribute to overall diabetes risk.

The metabolic syndrome (MetS) affects one quarter of the global adult population and one third of US adults (1, 2). MetS is characterized by a clustering of abnormalities in waist circumference, triglyceride levels, cholesterol, blood pressure, and glucose levels (3). Persons with MetS have up to a 5-fold greater risk of suffering from coronary heart disease, stroke, or diabetes (2, 4). Efforts are needed to identify and intervene on preventable risk factors for MetS development.

Inorganic arsenic (iAs) is a carcinogen that is associated with increased risk of numerous health effects, including cardiometabolic outcomes (522). After ingestion, iAs is converted to mono- and dimethylated arsenic compounds (monomethylarsonic acid (MMA) and dimethylarsinic acid (DMA)) that are excreted in the urine (2325). Arsenic metabolism is measured by computing relative percentages of urinary iAs, MMA, and DMA over their sum (iAs%, MMA%, and DMA%). Interindividual differences in arsenic metabolism are related to adverse health outcomes after controlling for arsenic exposure. Due to its relatively shorter half-life and rapid excretion through urine as compared with iAs, higher DMA% is considered a more efficient arsenic metabolism profile and protective against arsenic toxicity. However, increasing evidence indicates that the association between arsenic metabolism and disease is more complex. Indeed, while higher urinary MMA% and lower DMA% have been associated with greater risks of arsenic-related cancer (2631) and cardiovascular disease (3234), lower MMA% and higher DMA% have been cross-sectionally associated with higher body mass index (BMI) (35) and prospectively with greater risk of insulin resistance (36), diabetes (16, 3739), and MetS (11). The mechanism behind these contrasting associations is not clear, yet trends have been consistent across exposure levels and ethnicities (40).

MetS is understudied in the context of arsenic and provides a useful opportunity to understand the conflicting associations with metabolic outcomes reported for arsenic metabolism versus the associations reported with cardiovascular disease, because of the clustering of interrelated but distinct components. To our knowledge, no investigators to date have reported on the relationship of arsenic with both MetS and its individual components. We conducted a comprehensive evaluation of the prospective association of arsenic exposure and metabolism with MetS and its components in the Strong Heart Family Study, a family-based cohort study of members of American Indian tribes from Arizona, Oklahoma, and North/South Dakota. The study communities are affected by a high burden of MetS and are exposed to low-to-moderate levels of arsenic from drinking water and food.

METHODS

Study population

The Strong Heart Family Study investigators recruited 2,919 participants in 1998–1999 and 2001–2003. Participants recruited in 1998–1999 (n = 428) had follow-up visits in 2001–2003 and 2006–2009. Participants recruited in 2001–2003 (n = 2,491) had a single follow-up visit in 2006–2009. For this study, we included participants who were free of diabetes at baseline and had donated sufficient urine for arsenic analyses (n = 1,720). We used urinary arsenic level measured from baseline visits and MetS data from 1–2 follow-up visits. Participants missing information on education, smoking, alcohol intake, BMI, estimated glomerular filtration rate, and urinary creatinine concentration (n = 18) were excluded. We further excluded participants missing MetS data (n = 52) and prevalent MetS cases (n = 616), which resulted in 1,047 participants available for incident MetS analyses. For incidence of individual MetS components, we further excluded participants with prevalent cases of each component under analysis (Figure 1).

Selection of study participants for an analysis of the associations between arsenic exposure and the metabolic syndrome (MetS) and its components (elevated waist circumference (WC), elevated triglycerides (TG), reduced high-density lipoprotein cholesterol (HDL-C), hypertension (HT), and elevated fasting plasma glucose (FPG)), Strong Heart Family Study, 1998–2003.
Figure 1.

Selection of study participants for an analysis of the associations between arsenic exposure and the metabolic syndrome (MetS) and its components (elevated waist circumference (WC), elevated triglycerides (TG), reduced high-density lipoprotein cholesterol (HDL-C), hypertension (HT), and elevated fasting plasma glucose (FPG)), Strong Heart Family Study, 1998–2003.

All participants provided informed consent. Study protocols were approved by multiple institutional review boards, participating communities, and the Indian Health Service.

Data collection

Baseline and follow-up visits included biospecimen collection, physical examinations, completion of a food frequency questionnaire, and completion of an interview-administered questionnaire (collecting data on age, sex, education, smoking history, alcohol use, and medical history) (41). Measurement of waist circumference, blood pressure, height, and weight, as well as collection of urine and fasting blood samples, was performed during physical examinations by centrally trained nurses following a standardized protocol (42).

Urinary arsenic

Morning spot urine samples were frozen within 1–2 hours of collection and stored at −70°C. For arsenic analyses, up to 1.0 mL of urine from each sample was transported on dry ice to the Trace Element Laboratory at the University of Graz (Graz, Austria). Methods have been described elsewhere (43). Briefly, iAs, MMA, DMA, and arsenobetaine were measured using high-performance liquid chromatography/inductively coupled plasma mass spectrometry. Limits of detection were 0.1 μg/L for all 4 species. Arsenic species concentrations below the limit of detection (<5% for all species) were imputed as the limit of detection divided by √2. Arsenobetaine levels were low (median, 0.51 μg/L; interquartile range (IQR), 0.34–1.00), confirming infrequent seafood consumption. Urinary creatinine level was analyzed using automated alkaline picrate methodology.

Metabolic syndrome

Fasting triglycerides, cholesterol, and glucose were measured using enzymatic methods. High-density lipoprotein cholesterol was measured by precipitation with heparin and manganese chloride (42). Systolic and diastolic blood pressure measurements were taken on the right arm after 5 minutes’ rest using a Baum mercury sphygmomanometer, and the average of the last 2 of 3 measurements was calculated. Waist circumference was measured at the umbilicus with the participant in a supine position.

MetS was characterized according to the National Cholesterol Education Program Adult Treatment Panel III guidelines (3, 44), as recommended by the Indian Health Service (45), and was defined as 3 or more of the following: elevated waist circumference (≥40 inches (≥102 cm) for men, ≥35 inches (≥89 cm) for women); elevated triglyceride concentration (≥150 mg/dL (or use of medication)); elevated fasting plasma glucose concentration (≥100 mg/dL (or use of medication)); reduced high-density lipoprotein cholesterol concentration (<40 mg/dL for men and <50 mg/dL for women (or use of medication)); and hypertension (systolic blood pressure ≥130 mm Hg or diastolic blood pressure ≥85 mm Hg (or use of medication)).

Other variables

Estimated daily average dietary intake of 1-carbon metabolism micronutrients, including vitamin B2, vitamin B6, and folate, as well as total caloric intake and percentage of kilocalories derived from fat and protein, were measured at baseline through an interviewer-administered Block 119-item food frequency questionnaire (46).

Estimated glomerular filtration rate (mL/minute/1.73 m2; continuous) was calculated from plasma creatinine level, age, and sex using the Chronic Kidney Disease Epidemiology Collaboration formula (47). BMI was calculated as weight in kilometers divided by squared height in meters.

We selected single nucleotide polymorphism rs12768205 on the arsenite methyltransferase gene (AS3MT), identified using the Illumina Cardio-Metabo DNA Analysis BeadChip (MetaboChip; Illumina, Inc., San Diego, California), because it showed the strongest association with arsenic metabolism biomarkers in the Strong Heart Family Study (index single nucleotide polymorphism) (48).

Statistical analysis

Data on arsenic exposure, calculated as the sum of inorganic and methylated arsenic species (∑As) in urine, were right-skewed and log-transformed. Arsenic metabolism was computed by dividing each arsenic metabolite concentration over the sum of those species and multiplying by 100 (iAs%, MMA%, and DMA%). Differences in characteristics between participants with and without MetS were determined using Kruskal-Wallis and χ2 tests for continuous and categorical variables, respectively (Table 1). We assessed the prospective associations of baseline urinary arsenic measures (∑As, iAs%, MMA%, and DMA%) with incident MetS and its individual components. We employed modified Poisson regression with robust variance (49), using generalized estimating equations with an independence working correlation structure to account for family clustering. In the main analyses, we determined the relative risk and 95% confidence interval for overall MetS and each MetS component.

Table 1.

Baseline Characteristics of American Indian Participants According to Incident Metabolic Syndrome Status at Follow-up, Strong Heart Family Study, 1998–2009

Baseline CharacteristicMetS Status at Follow-upP Valuea
Total (n = 1,047)No MetS (n = 709 (67.7%))MetS (n = 338 (32.3%))
Median (IQR)%Median (IQR)%Median (IQR)%
Age, years30.7 (20.6–41.8)28.4 (19.5–40.9)34.0 (24.4–44.4)<0.001
Sex
 Female56.456.755.6
 Male43.643.344.40.80
Education, years
 <1234.736.530.8
 ≥1265.363.569.20.08
Smoking status
 Never smoker43.143.442.3
 Ever smoker17.717.318.3
 Current smoker39.339.239.30.91
Alcohol intake
 Never drinker11.412.49.2
 Ever drinker21.220.223.4
 Current drinker67.467.467.50.20
Body mass indexb27.4 (23.9–31.6)26.1 (22.9–30.1)30.1 (26.5–34.4)<0.001
Urinary creatinine concentration, mg/dL1.5 (1.0–2.1)1.6 (1.0–2.2)1.5 (1.1–2.1)0.42
eGFR, mL/minute/1.73 m2122 (111–134)123 (112–134)120 (108–132)0.02
∑As, μg/L6.5 (4.2–10.8)6.3 (4.2–10.4)7.0 (4.2–11.7)0.21
Biomarkers of arsenic metabolism
 iAs%10.6 (7.5–14.9)10.9 (7.6–15.3)9.9 (7.1–13.5)0.01
 MMA%15.4 (12.2–19.0)16.0 (12.7–19.8)14.2 (11.2–17.6)<0.001
 DMA%73.3 (66.3–78.9)72.1 (65.6–78.2)75.1 (69.1–80.5)<0.001
Components of MetS
 Waist circumference, inchesc36.6 (32.7–40.6)35.4 (31.5–39.0)39.0 (35.4–43.6)<0.001
  Elevated waist circumferenced45.238.160.1<0.001
 Triglycerides, mg/dL104 (80.0–135)96.0 (75.0–125)119.5 (95.2–146.8)<0.001
  Elevated triglyceridese15.912.423.1<0.001
 HDL cholesterol concentration, mg/dLf53.0 (44.0–62.0)54.0 (46.0–63.0)50.0 (42.0–59.8)<0.001
  Reduced HDL concentrationg24.420.931.7<0.001
 Systolic blood pressure, mm Hg115 (108–124)114 (107–122)119 (111–128)<0.001
 Diastolic blood pressure, mm Hg74.0 (67.0–80.0)73.0 (66.0–79.0)76 (69.0–82.0)<0.001
  Hypertensionh18.715.824.90.001
 FPG concentration, mg/dL90.0 (85.0–95.0)90.0 (85.0–95.0)92.0 (86.0–96.0)0.001
  Elevated FPG concentrationi10.510.99.80.67
Baseline CharacteristicMetS Status at Follow-upP Valuea
Total (n = 1,047)No MetS (n = 709 (67.7%))MetS (n = 338 (32.3%))
Median (IQR)%Median (IQR)%Median (IQR)%
Age, years30.7 (20.6–41.8)28.4 (19.5–40.9)34.0 (24.4–44.4)<0.001
Sex
 Female56.456.755.6
 Male43.643.344.40.80
Education, years
 <1234.736.530.8
 ≥1265.363.569.20.08
Smoking status
 Never smoker43.143.442.3
 Ever smoker17.717.318.3
 Current smoker39.339.239.30.91
Alcohol intake
 Never drinker11.412.49.2
 Ever drinker21.220.223.4
 Current drinker67.467.467.50.20
Body mass indexb27.4 (23.9–31.6)26.1 (22.9–30.1)30.1 (26.5–34.4)<0.001
Urinary creatinine concentration, mg/dL1.5 (1.0–2.1)1.6 (1.0–2.2)1.5 (1.1–2.1)0.42
eGFR, mL/minute/1.73 m2122 (111–134)123 (112–134)120 (108–132)0.02
∑As, μg/L6.5 (4.2–10.8)6.3 (4.2–10.4)7.0 (4.2–11.7)0.21
Biomarkers of arsenic metabolism
 iAs%10.6 (7.5–14.9)10.9 (7.6–15.3)9.9 (7.1–13.5)0.01
 MMA%15.4 (12.2–19.0)16.0 (12.7–19.8)14.2 (11.2–17.6)<0.001
 DMA%73.3 (66.3–78.9)72.1 (65.6–78.2)75.1 (69.1–80.5)<0.001
Components of MetS
 Waist circumference, inchesc36.6 (32.7–40.6)35.4 (31.5–39.0)39.0 (35.4–43.6)<0.001
  Elevated waist circumferenced45.238.160.1<0.001
 Triglycerides, mg/dL104 (80.0–135)96.0 (75.0–125)119.5 (95.2–146.8)<0.001
  Elevated triglyceridese15.912.423.1<0.001
 HDL cholesterol concentration, mg/dLf53.0 (44.0–62.0)54.0 (46.0–63.0)50.0 (42.0–59.8)<0.001
  Reduced HDL concentrationg24.420.931.7<0.001
 Systolic blood pressure, mm Hg115 (108–124)114 (107–122)119 (111–128)<0.001
 Diastolic blood pressure, mm Hg74.0 (67.0–80.0)73.0 (66.0–79.0)76 (69.0–82.0)<0.001
  Hypertensionh18.715.824.90.001
 FPG concentration, mg/dL90.0 (85.0–95.0)90.0 (85.0–95.0)92.0 (86.0–96.0)0.001
  Elevated FPG concentrationi10.510.99.80.67

Abbreviations: ∑As, sum of inorganic and methylated arsenic species; DMA, dimethylarsinic acid; eGFR, estimated glomerular filtration rate; FPG, fasting plasma glucose; HDL, high-density lipoprotein; iAs, inorganic arsenic; IQR, interquartile range; MetS, metabolic syndrome; MMA, monomethylarsonic acid.

aP values were determined by Kruskal-Wallis test for continuous variables and χ2 test for categorical variables.

b Weight (kg)/height (m)2.

c 1 inch = 2.54 cm.

d Meets criteria for the elevated waist circumference component of MetS (≥40 inches (≥102 cm) in men and ≥35 inches (≥89 cm) in women).

e Meets criteria for the elevated triglyceride component of MetS (≥150 mg/dL (or use of medication)).

f To convert HDL cholesterol units to mmol/L, multiply values by 0.0259.

g Meets criteria for the reduced HDL cholesterol component of MetS (<40 mg/dL for men and <50 mg/dL for women (or use of medication)).

h Meets criteria for the hypertensive component of MetS (systolic blood pressure ≥130 mm Hg or diastolic blood pressure ≥85 mm Hg (or use of medication)).

i Meets criteria for the elevated FPG component of MetS (≥100 mg/dL (or use of medication)).

Table 1.

Baseline Characteristics of American Indian Participants According to Incident Metabolic Syndrome Status at Follow-up, Strong Heart Family Study, 1998–2009

Baseline CharacteristicMetS Status at Follow-upP Valuea
Total (n = 1,047)No MetS (n = 709 (67.7%))MetS (n = 338 (32.3%))
Median (IQR)%Median (IQR)%Median (IQR)%
Age, years30.7 (20.6–41.8)28.4 (19.5–40.9)34.0 (24.4–44.4)<0.001
Sex
 Female56.456.755.6
 Male43.643.344.40.80
Education, years
 <1234.736.530.8
 ≥1265.363.569.20.08
Smoking status
 Never smoker43.143.442.3
 Ever smoker17.717.318.3
 Current smoker39.339.239.30.91
Alcohol intake
 Never drinker11.412.49.2
 Ever drinker21.220.223.4
 Current drinker67.467.467.50.20
Body mass indexb27.4 (23.9–31.6)26.1 (22.9–30.1)30.1 (26.5–34.4)<0.001
Urinary creatinine concentration, mg/dL1.5 (1.0–2.1)1.6 (1.0–2.2)1.5 (1.1–2.1)0.42
eGFR, mL/minute/1.73 m2122 (111–134)123 (112–134)120 (108–132)0.02
∑As, μg/L6.5 (4.2–10.8)6.3 (4.2–10.4)7.0 (4.2–11.7)0.21
Biomarkers of arsenic metabolism
 iAs%10.6 (7.5–14.9)10.9 (7.6–15.3)9.9 (7.1–13.5)0.01
 MMA%15.4 (12.2–19.0)16.0 (12.7–19.8)14.2 (11.2–17.6)<0.001
 DMA%73.3 (66.3–78.9)72.1 (65.6–78.2)75.1 (69.1–80.5)<0.001
Components of MetS
 Waist circumference, inchesc36.6 (32.7–40.6)35.4 (31.5–39.0)39.0 (35.4–43.6)<0.001
  Elevated waist circumferenced45.238.160.1<0.001
 Triglycerides, mg/dL104 (80.0–135)96.0 (75.0–125)119.5 (95.2–146.8)<0.001
  Elevated triglyceridese15.912.423.1<0.001
 HDL cholesterol concentration, mg/dLf53.0 (44.0–62.0)54.0 (46.0–63.0)50.0 (42.0–59.8)<0.001
  Reduced HDL concentrationg24.420.931.7<0.001
 Systolic blood pressure, mm Hg115 (108–124)114 (107–122)119 (111–128)<0.001
 Diastolic blood pressure, mm Hg74.0 (67.0–80.0)73.0 (66.0–79.0)76 (69.0–82.0)<0.001
  Hypertensionh18.715.824.90.001
 FPG concentration, mg/dL90.0 (85.0–95.0)90.0 (85.0–95.0)92.0 (86.0–96.0)0.001
  Elevated FPG concentrationi10.510.99.80.67
Baseline CharacteristicMetS Status at Follow-upP Valuea
Total (n = 1,047)No MetS (n = 709 (67.7%))MetS (n = 338 (32.3%))
Median (IQR)%Median (IQR)%Median (IQR)%
Age, years30.7 (20.6–41.8)28.4 (19.5–40.9)34.0 (24.4–44.4)<0.001
Sex
 Female56.456.755.6
 Male43.643.344.40.80
Education, years
 <1234.736.530.8
 ≥1265.363.569.20.08
Smoking status
 Never smoker43.143.442.3
 Ever smoker17.717.318.3
 Current smoker39.339.239.30.91
Alcohol intake
 Never drinker11.412.49.2
 Ever drinker21.220.223.4
 Current drinker67.467.467.50.20
Body mass indexb27.4 (23.9–31.6)26.1 (22.9–30.1)30.1 (26.5–34.4)<0.001
Urinary creatinine concentration, mg/dL1.5 (1.0–2.1)1.6 (1.0–2.2)1.5 (1.1–2.1)0.42
eGFR, mL/minute/1.73 m2122 (111–134)123 (112–134)120 (108–132)0.02
∑As, μg/L6.5 (4.2–10.8)6.3 (4.2–10.4)7.0 (4.2–11.7)0.21
Biomarkers of arsenic metabolism
 iAs%10.6 (7.5–14.9)10.9 (7.6–15.3)9.9 (7.1–13.5)0.01
 MMA%15.4 (12.2–19.0)16.0 (12.7–19.8)14.2 (11.2–17.6)<0.001
 DMA%73.3 (66.3–78.9)72.1 (65.6–78.2)75.1 (69.1–80.5)<0.001
Components of MetS
 Waist circumference, inchesc36.6 (32.7–40.6)35.4 (31.5–39.0)39.0 (35.4–43.6)<0.001
  Elevated waist circumferenced45.238.160.1<0.001
 Triglycerides, mg/dL104 (80.0–135)96.0 (75.0–125)119.5 (95.2–146.8)<0.001
  Elevated triglyceridese15.912.423.1<0.001
 HDL cholesterol concentration, mg/dLf53.0 (44.0–62.0)54.0 (46.0–63.0)50.0 (42.0–59.8)<0.001
  Reduced HDL concentrationg24.420.931.7<0.001
 Systolic blood pressure, mm Hg115 (108–124)114 (107–122)119 (111–128)<0.001
 Diastolic blood pressure, mm Hg74.0 (67.0–80.0)73.0 (66.0–79.0)76 (69.0–82.0)<0.001
  Hypertensionh18.715.824.90.001
 FPG concentration, mg/dL90.0 (85.0–95.0)90.0 (85.0–95.0)92.0 (86.0–96.0)0.001
  Elevated FPG concentrationi10.510.99.80.67

Abbreviations: ∑As, sum of inorganic and methylated arsenic species; DMA, dimethylarsinic acid; eGFR, estimated glomerular filtration rate; FPG, fasting plasma glucose; HDL, high-density lipoprotein; iAs, inorganic arsenic; IQR, interquartile range; MetS, metabolic syndrome; MMA, monomethylarsonic acid.

aP values were determined by Kruskal-Wallis test for continuous variables and χ2 test for categorical variables.

b Weight (kg)/height (m)2.

c 1 inch = 2.54 cm.

d Meets criteria for the elevated waist circumference component of MetS (≥40 inches (≥102 cm) in men and ≥35 inches (≥89 cm) in women).

e Meets criteria for the elevated triglyceride component of MetS (≥150 mg/dL (or use of medication)).

f To convert HDL cholesterol units to mmol/L, multiply values by 0.0259.

g Meets criteria for the reduced HDL cholesterol component of MetS (<40 mg/dL for men and <50 mg/dL for women (or use of medication)).

h Meets criteria for the hypertensive component of MetS (systolic blood pressure ≥130 mm Hg or diastolic blood pressure ≥85 mm Hg (or use of medication)).

i Meets criteria for the elevated FPG component of MetS (≥100 mg/dL (or use of medication)).

The associations of arsenic exposure with incident MetS and MetS components are reported per IQR increase and with restricted cubic splines to allow for flexibility in the dose-response. For arsenic metabolism, we conducted 3 types of analyses: 1) conventional model, 2) leave-one-out, and 3) principal component (PC). In conventional models, we evaluated each metabolite separately, reporting the relative risk for each outcome per 5–percentage-point increase in each arsenic species. In leave-one-out models, for each arsenic species, we also included one of the species not being evaluated. For example, when evaluating iAs%, if we include MMA%, the regression coefficient estimates the relative risk associated with substituting DMA% for the equivalent iAs% while holding constant MMA% and ∑As. This leave-one-out approach has been used previously for arsenic metabolism (37). Finally, we evaluated arsenic metabolism using PC analyses, as recently proposed (48, 50). PC analyses were conducted for arsenic species percentages after normalizing them to have a mean value of 0 and a standard deviation of 1. This method is useful for arsenic metabolism because it removes the interdependence between the 3 biomarkers, with the 2 resulting PCs allowing for potentially more biologically meaningful interpretation, by summarizing primary and secondary methylation steps, respectively (48, 50). Although we cannot be certain that the PCs truly reflect these methylation steps, since arsenic metabolism in the body is still being debated, together with leave-one-out models, this modeling strategy contributes to a more comprehensive picture of how the 3 metabolites interact with respect to their association with health outcomes.

We used progressive adjustments for known determinants of arsenic metabolism and MetS. For ∑As models, model 1 adjusted for urinary creatinine concentration. Model 2 further adjusted for age (years; continuous), sex, region (Arizona, Oklahoma, or North/South Dakota), and education (<12 years or ≥12 years (or above/below the appropriate level of schooling if aged <18 years)). Model 3 further adjusted for BMI (<25, 25–29.9, 30–34.9, 35–39.9, or ≥40), smoking status (never, former, or current smoker), alcohol use (never, former, or current drinker), and estimated glomerular filtration rate. For arsenic metabolism models, model 1 additionally adjusted for log ∑As. In models evaluating waist circumference, BMI was not included as an adjustment factor, for consistency with previous literature (51).

We analyzed possible effect modification of the association between arsenic metabolism and MetS by region, education, smoking status, BMI, ∑As, age, B vitamin intake, and AS3MT genotype by adding interaction terms. Interaction P values were obtained using Wald tests for multiple coefficients. Participants who were missing information on AS3MT genotype (n = 12) or vitamin intake (n = 75) were excluded from interaction analyses.

We conducted several sensitivity analyses for arsenic metabolism and MetS models. First, we adjusted for number of MetS events present at baseline (0, 1, or 2), since the greater the number of baseline abnormalities the fewer “new” abnormalities there are to be gained for development of MetS. Second, we adjusted for dietary 1-carbon metabolism vitamin intake and estimated daily caloric intake (both log-transformed). Third, we adjusted for percentage of kilocalories derived from fat and protein. Finally, we adjusted for baseline waist circumference in models evaluating arsenic metabolism and elevated waist circumference.

RESULTS

Participant characteristics

Of 1,047 participants without MetS at baseline, 338 (32%) developed MetS over a median follow-up period of 5.3 (IQR, 4.7–6.5) years. The median age was 30.7 years; 56.4% of participants were female and 43.6% were male. Median urinary ∑As concentration was 6.5 μg/L, and median urinary concentrations of iAs%, MMA%, and DMA% were 10.6%, 15.4%, and 73.3%, respectively (Table 1). Compared with participants who remained free of MetS throughout follow-up, those who developed MetS were older, were more educated, had higher BMI and DMA%, and had lower estimated glomerular filtration rate, iAs%, and MMA% at baseline (P < 0.05). The first arsenic metabolism PC explained 80.1% of the variance in arsenic species and reflected higher DMA% and lower iAs% and MMA%, possibly representing overall methylation to DMA, or the secondary arsenic methylation step (Table 2). The second arsenic metabolism PC explained 19.9% of the variance and reflected higher MMA% and lower iAs%, independent of DMA%, potentially reflecting the primary methylation step.

Table 2.

Principal Components of Urinary Arsenic Species Measured in the Strong Heart Family Study, 1998–2009

VariablePrincipal Component
12
Variance in arsenic species explained, %a80.1 (1.6)19.9 (0.8)
Weight
 iAs%−0.55−0.67
 MMA%−0.530.74
 DMA%0.640.04
VariablePrincipal Component
12
Variance in arsenic species explained, %a80.1 (1.6)19.9 (0.8)
Weight
 iAs%−0.55−0.67
 MMA%−0.530.74
 DMA%0.640.04

Abbreviations: DMA, dimethylarsinic acid; iAs, inorganic arsenic; MMA, monomethylarsonic acid.

a Values are presented as mean (standard deviation).

Table 2.

Principal Components of Urinary Arsenic Species Measured in the Strong Heart Family Study, 1998–2009

VariablePrincipal Component
12
Variance in arsenic species explained, %a80.1 (1.6)19.9 (0.8)
Weight
 iAs%−0.55−0.67
 MMA%−0.530.74
 DMA%0.640.04
VariablePrincipal Component
12
Variance in arsenic species explained, %a80.1 (1.6)19.9 (0.8)
Weight
 iAs%−0.55−0.67
 MMA%−0.530.74
 DMA%0.640.04

Abbreviations: DMA, dimethylarsinic acid; iAs, inorganic arsenic; MMA, monomethylarsonic acid.

a Values are presented as mean (standard deviation).

Associations of arsenic exposure and metabolism with incident MetS

In adjusted models, ∑As was not associated with incident MetS (per IQR increase in ∑As, relative risk (RR) = 1.03, 95% confidence interval (CI): 0.90, 1.18) (Table 3). For arsenic metabolism, in adjusted conventional models, the relative risks of incident MetS per 5% increase in iAs%, MMA%, and DMA% were 0.94 (95% CI: 0.88, 1.01), 0.87 (95% CI: 0.79, 0.95), and 1.07 (95% CI: 1.02, 1.12), respectively. In iAs-fixed leave-one-out models for MMA% and DMA%, the relative risk remained the same for MMA% (RR = 0.87, 95% CI: 0.79, 0.97) but strengthened to 1.14 (95% CI: 1.03, 1.27) for DMA% (Table 3). Consistently, in flexible dose-response analyses, DMA% showed a linear association with MetS when iAs% was included in the model (see Web Figure 1, available at https://academic.oup.com/aje). The association for DMA%, however, was attenuated in the MMA%-fixed leave-one-out model. In both leave-one-out models for iAs% (fixing either MMA% or DMA%), there remained no significant association with MetS. When modeling arsenic metabolism using PC analysis, an IQR increase in the first arsenic metabolism PC (greater DMA%, lower MMA% and iAs%) was associated with a 1.19-fold (95% CI: 1.06, 1.34) increased risk of MetS, while the second arsenic metabolism PC (greater MMA%, lower iAs%, independent of DMA%) was associated with a 0.93-fold (95% CI: 0.84, 1.03) decreased risk (Table 3).

Table 3.

Relative Risk of Incident Metabolic Syndrome According to Increases in the Sum of Inorganic and Methylated Arsenic Species and Biomarkers of Arsenic Metabolism (iAs%, MMA%, DMA%, and Principal Components) (n = 1,047), Strong Heart Family Study, 1998–2009

VariableModel 1aModel 2bModel 3c
RR95% CIRR95% CIRR95% CI
∑As, per IQR increased1.130.99, 1.291.050.92, 1.191.030.90, 1.18
iAs% (per 5% increase)
 Conventional model0.900.84, 0.970.910.85, 0.970.940.88, 1.01
 Leave-one-out models
  MMA% fixed0.970.91, 1.040.970.90, 1.040.980.91, 1.05
  DMA% fixed1.211.04, 1.401.191.02, 1.391.120.96, 1.30
MMA% (per 5% increase)
 Conventional model0.790.72, 0.870.800.72, 0.880.870.79, 0.95
 Leave-one-out models
  iAs% fixed0.810.73, 0.890.810.73, 0.910.870.79, 0.97
  DMA% fixed0.830.71, 0.960.840.72, 0.980.890.77, 1.04
DMA% (per 5% increase)
 Conventional model1.111.06, 1.171.111.06, 1.171.071.02, 1.12
 Leave-one-out-models
  iAs% fixed1.241.12, 1.381.231.10, 1.371.141.03, 1.27
  MMA% fixed1.030.96, 1.101.040.96, 1.111.020.95, 1.10
Arsenic metabolism PCs
 PC 1, per IQR increasee1.321.17, 1.481.321.16, 1.491.191.06, 1.34
 PC 2, per IQR increasef0.900.82, 0.980.910.82, 1.000.930.84, 1.03
VariableModel 1aModel 2bModel 3c
RR95% CIRR95% CIRR95% CI
∑As, per IQR increased1.130.99, 1.291.050.92, 1.191.030.90, 1.18
iAs% (per 5% increase)
 Conventional model0.900.84, 0.970.910.85, 0.970.940.88, 1.01
 Leave-one-out models
  MMA% fixed0.970.91, 1.040.970.90, 1.040.980.91, 1.05
  DMA% fixed1.211.04, 1.401.191.02, 1.391.120.96, 1.30
MMA% (per 5% increase)
 Conventional model0.790.72, 0.870.800.72, 0.880.870.79, 0.95
 Leave-one-out models
  iAs% fixed0.810.73, 0.890.810.73, 0.910.870.79, 0.97
  DMA% fixed0.830.71, 0.960.840.72, 0.980.890.77, 1.04
DMA% (per 5% increase)
 Conventional model1.111.06, 1.171.111.06, 1.171.071.02, 1.12
 Leave-one-out-models
  iAs% fixed1.241.12, 1.381.231.10, 1.371.141.03, 1.27
  MMA% fixed1.030.96, 1.101.040.96, 1.111.020.95, 1.10
Arsenic metabolism PCs
 PC 1, per IQR increasee1.321.17, 1.481.321.16, 1.491.191.06, 1.34
 PC 2, per IQR increasef0.900.82, 0.980.910.82, 1.000.930.84, 1.03

Abbreviations: ∑As, sum of inorganic and methylated arsenic species; CI, confidence interval; DMA, dimethylarsinic acid; iAs, inorganic arsenic; IQR, interquartile range; MMA, monomethylarsonic acid; PC, principal component; RR, relative risk.

a Model 1 adjusted for log total arsenic (except in models where ∑As was the predictor of interest) and urinary creatinine concentration.

b Model 2 further adjusted for sex, region, and education.

c Model 3 further adjusted for alcohol intake, smoking status, kidney function, and body mass index.

d 10.84 μg/L (75th percentile) vs. 4.20 μg/L (25th percentile).

e 1.04 (75th percentile) vs. −0.92 (25th percentile).

f 0.50 (75th percentile) vs. −0.38 (25th percentile).

Table 3.

Relative Risk of Incident Metabolic Syndrome According to Increases in the Sum of Inorganic and Methylated Arsenic Species and Biomarkers of Arsenic Metabolism (iAs%, MMA%, DMA%, and Principal Components) (n = 1,047), Strong Heart Family Study, 1998–2009

VariableModel 1aModel 2bModel 3c
RR95% CIRR95% CIRR95% CI
∑As, per IQR increased1.130.99, 1.291.050.92, 1.191.030.90, 1.18
iAs% (per 5% increase)
 Conventional model0.900.84, 0.970.910.85, 0.970.940.88, 1.01
 Leave-one-out models
  MMA% fixed0.970.91, 1.040.970.90, 1.040.980.91, 1.05
  DMA% fixed1.211.04, 1.401.191.02, 1.391.120.96, 1.30
MMA% (per 5% increase)
 Conventional model0.790.72, 0.870.800.72, 0.880.870.79, 0.95
 Leave-one-out models
  iAs% fixed0.810.73, 0.890.810.73, 0.910.870.79, 0.97
  DMA% fixed0.830.71, 0.960.840.72, 0.980.890.77, 1.04
DMA% (per 5% increase)
 Conventional model1.111.06, 1.171.111.06, 1.171.071.02, 1.12
 Leave-one-out-models
  iAs% fixed1.241.12, 1.381.231.10, 1.371.141.03, 1.27
  MMA% fixed1.030.96, 1.101.040.96, 1.111.020.95, 1.10
Arsenic metabolism PCs
 PC 1, per IQR increasee1.321.17, 1.481.321.16, 1.491.191.06, 1.34
 PC 2, per IQR increasef0.900.82, 0.980.910.82, 1.000.930.84, 1.03
VariableModel 1aModel 2bModel 3c
RR95% CIRR95% CIRR95% CI
∑As, per IQR increased1.130.99, 1.291.050.92, 1.191.030.90, 1.18
iAs% (per 5% increase)
 Conventional model0.900.84, 0.970.910.85, 0.970.940.88, 1.01
 Leave-one-out models
  MMA% fixed0.970.91, 1.040.970.90, 1.040.980.91, 1.05
  DMA% fixed1.211.04, 1.401.191.02, 1.391.120.96, 1.30
MMA% (per 5% increase)
 Conventional model0.790.72, 0.870.800.72, 0.880.870.79, 0.95
 Leave-one-out models
  iAs% fixed0.810.73, 0.890.810.73, 0.910.870.79, 0.97
  DMA% fixed0.830.71, 0.960.840.72, 0.980.890.77, 1.04
DMA% (per 5% increase)
 Conventional model1.111.06, 1.171.111.06, 1.171.071.02, 1.12
 Leave-one-out-models
  iAs% fixed1.241.12, 1.381.231.10, 1.371.141.03, 1.27
  MMA% fixed1.030.96, 1.101.040.96, 1.111.020.95, 1.10
Arsenic metabolism PCs
 PC 1, per IQR increasee1.321.17, 1.481.321.16, 1.491.191.06, 1.34
 PC 2, per IQR increasef0.900.82, 0.980.910.82, 1.000.930.84, 1.03

Abbreviations: ∑As, sum of inorganic and methylated arsenic species; CI, confidence interval; DMA, dimethylarsinic acid; iAs, inorganic arsenic; IQR, interquartile range; MMA, monomethylarsonic acid; PC, principal component; RR, relative risk.

a Model 1 adjusted for log total arsenic (except in models where ∑As was the predictor of interest) and urinary creatinine concentration.

b Model 2 further adjusted for sex, region, and education.

c Model 3 further adjusted for alcohol intake, smoking status, kidney function, and body mass index.

d 10.84 μg/L (75th percentile) vs. 4.20 μg/L (25th percentile).

e 1.04 (75th percentile) vs. −0.92 (25th percentile).

f 0.50 (75th percentile) vs. −0.38 (25th percentile).

In interaction analyses of the association of arsenic metabolism with incident MetS according to subgroup characteristics (Table 4), there were interactions between BMI and DMA% (P for interaction = 0.03) and between the AS3MT rs12768205 polymorphism and MMA% (P for interaction = 0.05). The association between DMA% and MetS was weaker with increasing BMI category, and the inverse association between MMA% and MetS was markedly stronger for participants with the rs12768205 AA genotype.

Table 4.

Relative Risk of Incident Metabolic Syndrome According to Increases in Biomarkers of Arsenic Metabolism (iAs%, MMA%, and DMA%) and Levels of Other Risk Factors (n = 972), Strong Heart Family Study, 1998–2009a

CharacteristicNo. of CasesNo. of ControlsBiomarker of Arsenic Metabolism
iAs%MMA%DMA%
RR95% CIP ValuebRR95% CIP ValueRR95% CIP Value
Age, years0.320.120.18
 ≤30.01193530.990.90, 1.090.930.80, 1.091.020.95, 1.11
 30.1–49.91402400.880.78, 0.990.830.73, 0.931.121.04, 1.20
 ≥50.049710.970.78, 1.210.730.59, 0.891.130.98, 1.30
Sex0.580.310.40
 Male1312870.920.83, 1.030.80.70, 0.921.111.02, 1.19
 Female1773770.960.87, 1.050.880.78, 1.001.060.99, 1.13
Education, years0.440.940.59
 ≤12972410.980.87, 1.090.850.70, 1.021.060.97, 1.15
 >122114230.920.84, 1.010.840.76, 0.941.091.02, 1.16
Smoking status0.970.710.95
 Never smoker1312840.950.84, 1.070.810.70, 0.931.091.00, 1.19
 Ever smoker591190.950.80, 1.130.860.69, 1.061.070.95, 1.20
 Current smoker1182610.930.83, 1.040.880.76, 1.011.071.00, 1.16
Body mass indexc0.110.080.03
 <25.0 (normal)502800.820.66, 1.020.710.57, 0.891.211.06, 1.38
 25.0–29.9 (overweight)1042070.900.79, 1.020.820.70, 0.951.111.02, 1.20
 ≥30.0 (obese)1541771.020.92, 1.120.920.81, 1.051.020.95, 1.08
Folate level, μg0.380.380.28
 <3521433170.980.87, 1.090.880.76, 1.021.050.97, 1.13
 ≥3521653470.910.84, 1.000.810.73, 0.911.111.04, 1.17
Vitamin B6 level, mg0.190.710.30
 <1.71513330.990.89, 1.120.860.76, 0.981.050.98, 1.13
 ≥1.71573310.900.81, 0.990.840.73, 0.941.111.03, 1.17
Vitamin B2 level, mg0.510.410.36
 <1.71543300.960.87, 1.070.890.78, 1.021.050.98, 1.13
 ≥1.71543340.920.84, 1.020.830.73, 0.931.091.03, 1.17
∑As concentration, μg/L0.440.990.79
 ≤4.19771680.940.82, 1.090.840.69, 1.021.080.98, 1.20
 4.20–6.57671871.050.90, 1.250.860.70, 1.061.020.91, 1.14
 6.58–10.83821530.900.77, 1.050.840.70, 1.011.090.99, 1.22
 ≥10.84821560.900.78 1.020.860.74, 1.001.091.00, 1.19
AS3MT genotype0.120.050.16
 GG1543290.980.90, 1.060.900.80, 1.001.040.98, 1.11
 GA1322630.840.73, 0.970.800.69, 0.941.141.05, 1.25
 AA22721.111.75, 0.700.540.35, 0.841.150.89, 1.49
Overall3086640.940.87, 1.010.840.76, 0.931.081.02, 1.14
CharacteristicNo. of CasesNo. of ControlsBiomarker of Arsenic Metabolism
iAs%MMA%DMA%
RR95% CIP ValuebRR95% CIP ValueRR95% CIP Value
Age, years0.320.120.18
 ≤30.01193530.990.90, 1.090.930.80, 1.091.020.95, 1.11
 30.1–49.91402400.880.78, 0.990.830.73, 0.931.121.04, 1.20
 ≥50.049710.970.78, 1.210.730.59, 0.891.130.98, 1.30
Sex0.580.310.40
 Male1312870.920.83, 1.030.80.70, 0.921.111.02, 1.19
 Female1773770.960.87, 1.050.880.78, 1.001.060.99, 1.13
Education, years0.440.940.59
 ≤12972410.980.87, 1.090.850.70, 1.021.060.97, 1.15
 >122114230.920.84, 1.010.840.76, 0.941.091.02, 1.16
Smoking status0.970.710.95
 Never smoker1312840.950.84, 1.070.810.70, 0.931.091.00, 1.19
 Ever smoker591190.950.80, 1.130.860.69, 1.061.070.95, 1.20
 Current smoker1182610.930.83, 1.040.880.76, 1.011.071.00, 1.16
Body mass indexc0.110.080.03
 <25.0 (normal)502800.820.66, 1.020.710.57, 0.891.211.06, 1.38
 25.0–29.9 (overweight)1042070.900.79, 1.020.820.70, 0.951.111.02, 1.20
 ≥30.0 (obese)1541771.020.92, 1.120.920.81, 1.051.020.95, 1.08
Folate level, μg0.380.380.28
 <3521433170.980.87, 1.090.880.76, 1.021.050.97, 1.13
 ≥3521653470.910.84, 1.000.810.73, 0.911.111.04, 1.17
Vitamin B6 level, mg0.190.710.30
 <1.71513330.990.89, 1.120.860.76, 0.981.050.98, 1.13
 ≥1.71573310.900.81, 0.990.840.73, 0.941.111.03, 1.17
Vitamin B2 level, mg0.510.410.36
 <1.71543300.960.87, 1.070.890.78, 1.021.050.98, 1.13
 ≥1.71543340.920.84, 1.020.830.73, 0.931.091.03, 1.17
∑As concentration, μg/L0.440.990.79
 ≤4.19771680.940.82, 1.090.840.69, 1.021.080.98, 1.20
 4.20–6.57671871.050.90, 1.250.860.70, 1.061.020.91, 1.14
 6.58–10.83821530.900.77, 1.050.840.70, 1.011.090.99, 1.22
 ≥10.84821560.900.78 1.020.860.74, 1.001.091.00, 1.19
AS3MT genotype0.120.050.16
 GG1543290.980.90, 1.060.900.80, 1.001.040.98, 1.11
 GA1322630.840.73, 0.970.800.69, 0.941.141.05, 1.25
 AA22721.111.75, 0.700.540.35, 0.841.150.89, 1.49
Overall3086640.940.87, 1.010.840.76, 0.931.081.02, 1.14

Abbreviations: ∑As, sum of inorganic and methylated arsenic species; AS3MT, arsenite methyltransferase gene; CI, confidence interval; DMA, dimethylarsinic acid; iAs, inorganic arsenic; MMA, monomethylarsonic acid; RR, relative risk.

a Models adjusted for log urinary total arsenic, log urinary creatinine, age, sex, region, smoking status, alcohol intake, B vitamin intake, caloric intake, education, body mass index, and kidney function.

bP values for interaction were obtained using Wald tests for multiple coefficients.

c Weight (kg)/height (m)2.

Table 4.

Relative Risk of Incident Metabolic Syndrome According to Increases in Biomarkers of Arsenic Metabolism (iAs%, MMA%, and DMA%) and Levels of Other Risk Factors (n = 972), Strong Heart Family Study, 1998–2009a

CharacteristicNo. of CasesNo. of ControlsBiomarker of Arsenic Metabolism
iAs%MMA%DMA%
RR95% CIP ValuebRR95% CIP ValueRR95% CIP Value
Age, years0.320.120.18
 ≤30.01193530.990.90, 1.090.930.80, 1.091.020.95, 1.11
 30.1–49.91402400.880.78, 0.990.830.73, 0.931.121.04, 1.20
 ≥50.049710.970.78, 1.210.730.59, 0.891.130.98, 1.30
Sex0.580.310.40
 Male1312870.920.83, 1.030.80.70, 0.921.111.02, 1.19
 Female1773770.960.87, 1.050.880.78, 1.001.060.99, 1.13
Education, years0.440.940.59
 ≤12972410.980.87, 1.090.850.70, 1.021.060.97, 1.15
 >122114230.920.84, 1.010.840.76, 0.941.091.02, 1.16
Smoking status0.970.710.95
 Never smoker1312840.950.84, 1.070.810.70, 0.931.091.00, 1.19
 Ever smoker591190.950.80, 1.130.860.69, 1.061.070.95, 1.20
 Current smoker1182610.930.83, 1.040.880.76, 1.011.071.00, 1.16
Body mass indexc0.110.080.03
 <25.0 (normal)502800.820.66, 1.020.710.57, 0.891.211.06, 1.38
 25.0–29.9 (overweight)1042070.900.79, 1.020.820.70, 0.951.111.02, 1.20
 ≥30.0 (obese)1541771.020.92, 1.120.920.81, 1.051.020.95, 1.08
Folate level, μg0.380.380.28
 <3521433170.980.87, 1.090.880.76, 1.021.050.97, 1.13
 ≥3521653470.910.84, 1.000.810.73, 0.911.111.04, 1.17
Vitamin B6 level, mg0.190.710.30
 <1.71513330.990.89, 1.120.860.76, 0.981.050.98, 1.13
 ≥1.71573310.900.81, 0.990.840.73, 0.941.111.03, 1.17
Vitamin B2 level, mg0.510.410.36
 <1.71543300.960.87, 1.070.890.78, 1.021.050.98, 1.13
 ≥1.71543340.920.84, 1.020.830.73, 0.931.091.03, 1.17
∑As concentration, μg/L0.440.990.79
 ≤4.19771680.940.82, 1.090.840.69, 1.021.080.98, 1.20
 4.20–6.57671871.050.90, 1.250.860.70, 1.061.020.91, 1.14
 6.58–10.83821530.900.77, 1.050.840.70, 1.011.090.99, 1.22
 ≥10.84821560.900.78 1.020.860.74, 1.001.091.00, 1.19
AS3MT genotype0.120.050.16
 GG1543290.980.90, 1.060.900.80, 1.001.040.98, 1.11
 GA1322630.840.73, 0.970.800.69, 0.941.141.05, 1.25
 AA22721.111.75, 0.700.540.35, 0.841.150.89, 1.49
Overall3086640.940.87, 1.010.840.76, 0.931.081.02, 1.14
CharacteristicNo. of CasesNo. of ControlsBiomarker of Arsenic Metabolism
iAs%MMA%DMA%
RR95% CIP ValuebRR95% CIP ValueRR95% CIP Value
Age, years0.320.120.18
 ≤30.01193530.990.90, 1.090.930.80, 1.091.020.95, 1.11
 30.1–49.91402400.880.78, 0.990.830.73, 0.931.121.04, 1.20
 ≥50.049710.970.78, 1.210.730.59, 0.891.130.98, 1.30
Sex0.580.310.40
 Male1312870.920.83, 1.030.80.70, 0.921.111.02, 1.19
 Female1773770.960.87, 1.050.880.78, 1.001.060.99, 1.13
Education, years0.440.940.59
 ≤12972410.980.87, 1.090.850.70, 1.021.060.97, 1.15
 >122114230.920.84, 1.010.840.76, 0.941.091.02, 1.16
Smoking status0.970.710.95
 Never smoker1312840.950.84, 1.070.810.70, 0.931.091.00, 1.19
 Ever smoker591190.950.80, 1.130.860.69, 1.061.070.95, 1.20
 Current smoker1182610.930.83, 1.040.880.76, 1.011.071.00, 1.16
Body mass indexc0.110.080.03
 <25.0 (normal)502800.820.66, 1.020.710.57, 0.891.211.06, 1.38
 25.0–29.9 (overweight)1042070.900.79, 1.020.820.70, 0.951.111.02, 1.20
 ≥30.0 (obese)1541771.020.92, 1.120.920.81, 1.051.020.95, 1.08
Folate level, μg0.380.380.28
 <3521433170.980.87, 1.090.880.76, 1.021.050.97, 1.13
 ≥3521653470.910.84, 1.000.810.73, 0.911.111.04, 1.17
Vitamin B6 level, mg0.190.710.30
 <1.71513330.990.89, 1.120.860.76, 0.981.050.98, 1.13
 ≥1.71573310.900.81, 0.990.840.73, 0.941.111.03, 1.17
Vitamin B2 level, mg0.510.410.36
 <1.71543300.960.87, 1.070.890.78, 1.021.050.98, 1.13
 ≥1.71543340.920.84, 1.020.830.73, 0.931.091.03, 1.17
∑As concentration, μg/L0.440.990.79
 ≤4.19771680.940.82, 1.090.840.69, 1.021.080.98, 1.20
 4.20–6.57671871.050.90, 1.250.860.70, 1.061.020.91, 1.14
 6.58–10.83821530.900.77, 1.050.840.70, 1.011.090.99, 1.22
 ≥10.84821560.900.78 1.020.860.74, 1.001.091.00, 1.19
AS3MT genotype0.120.050.16
 GG1543290.980.90, 1.060.900.80, 1.001.040.98, 1.11
 GA1322630.840.73, 0.970.800.69, 0.941.141.05, 1.25
 AA22721.111.75, 0.700.540.35, 0.841.150.89, 1.49
Overall3086640.940.87, 1.010.840.76, 0.931.081.02, 1.14

Abbreviations: ∑As, sum of inorganic and methylated arsenic species; AS3MT, arsenite methyltransferase gene; CI, confidence interval; DMA, dimethylarsinic acid; iAs, inorganic arsenic; MMA, monomethylarsonic acid; RR, relative risk.

a Models adjusted for log urinary total arsenic, log urinary creatinine, age, sex, region, smoking status, alcohol intake, B vitamin intake, caloric intake, education, body mass index, and kidney function.

bP values for interaction were obtained using Wald tests for multiple coefficients.

c Weight (kg)/height (m)2.

In sensitivity analyses including additional adjustments for baseline number of MetS criteria and B vitamins (Web Table 1), as well as percentage of kilocalories derived from fat and protein (data not shown), we observed consistent results.

Associations of arsenic exposure and metabolism with incident MetS components

In adjusted models, ∑As was associated with elevated fasting plasma glucose concentration (RR = 1.19, 95% CI: 1.01, 1.41) (Table 5) but not with any other MetS component. ∑As remained associated with elevated fasting plasma glucose level in flexible dose-response models (Figure 2).

Table 5.

Relative Risk of Incident Components of the Metabolic Syndrome According to Increases in the Sum of Inorganic and Methylated Arsenic Species and Biomarkers of Arsenic Metabolism (iAs%, MMA%, and DMA%), Strong Heart Family Study, 1998–2009

Component, Arsenic Biomarker, and ModelModel 1aModel 2bModel 3c
RR95% CIRR95% CIRR95% CI
Elevated waist circumference (n = 574)d
 ∑As, per IQR increasee0.990.86, 1.140.990.86, 1.140.980.85, 1.13
 iAs%
  Conventional model0.960.88, 1.040.970.89, 1.060.970.89, 1.06
  LOO model (MMA% fixed)1.020.94, 1.121.020.94, 1.121.030.94, 1.12
  LOO model (DMA% fixed)1.291.10, 1.511.251.07, 1.451.251.08, 1.45
 MMA%
  Conventional model0.800.74, 0.880.830.76, 0.910.830.76, 0.91
  LOO model (iAs% fixed)0.800.72, 0.880.820.75, 0.910.820.75, 0.91
  LOO model (DMA% fixed)0.780.66, 0.910.800.69, 0.930.800.69, 0.93
 DMA%
  Conventional model1.081.03, 1.141.071.01, 1.131.071.01, 1.13
  LOO model (iAs% fixed)1.261.14, 1.391.221.10, 1.341.221.10, 1.34
  LOO model (MMA% fixed)0.980.89, 1.060.980.90, 1.060.970.89, 1.06
Elevated triglyceride concentration (n = 881)f
 ∑As, per IQR increasee0.960.81, 1.141.000.84, 1.190.980.82, 1.19
 iAs%
  Conventional model0.950.86, 1.040.930.84, 1.030.930.84, 1.03
  LOO model (MMA% fixed)0.960.87, 1.060.950.86, 1.060.940.84, 1.05
  LOO model (DMA% fixed)0.990.82, 1.181.010.84, 1.220.990.81, 1.19
 MMA%
  Conventional model0.950.85, 1.060.920.82, 1.030.930.83, 1.04
  LOO model (iAs% fixed)0.970.86, 1.090.940.83, 1.060.960.84, 1.08
  LOO model (DMA% fixed)1.010.84, 1.210.990.82, 1.191.010.84, 1.23
 DMA%
  Conventional model1.040.98, 1.101.060.99, 1.131.050.99, 1.13
  LOO model (iAs% fixed)1.030.92, 1.161.060.94, 1.201.050.92, 1.18
  LOO model (MMA% fixed)1.040.94,1.151.050.94, 1.171.060.95, 1.18
Reduced HDL cholesterol concentration (n = 792)g
 ∑As, per IQR increasee1.070.91, 1.250.960.81, 1.140.950.80, 1.14
 iAs%
  Conventional model0.890.81, 0.980.900.82, 0.990.920.84, 1.02
  LOO model (MMA% fixed)0.940.86, 1.040.930.85, 1.030.940.85, 1.04
  LOO model (DMA% fixed)1.110.93, 1.341.040.86, 1.251.000.83, 1.22
 MMA%
  Conventional model0.820.73, 0.930.870.77, 0.980.910.81, 1.03
  LOO model (iAs% fixed)0.850.75, 0.960.900.79, 1.020.940.82, 1.07
  LOO model (DMA% fixed)0.900.75, 1.080.960.80, 1.161.000.82, 1.21
 DMA%
  Conventional model1.111.04, 1.181.091.02, 1.161.061.00, 1.13
  LOO model (iAs% fixed)1.181.04, 1.341.120.98, 1.271.060.93, 1.22
  LOO model (MMA% fixed)1.060.96, 1.161.070.97, 1.181.060.96, 1.17
Hypertension (n = 851)h
 ∑As, per IQR increasee1.110.91, 1.361.050.90, 1.231.030.88, 1.20
 iAs%
  Conventional model1.000.93, 1.090.980.90, 1.061.000.92, 1.08
  LOO model (MMA% fixed)1.000.91, 1.091.000.91, 1.101.010.91, 1.12
  LOO model (DMA% fixed)0.970.80, 1.181.090.89, 1.331.050.85, 1.30
 MMA%
  Conventional model1.020.90, 1.160.920.83, 1.030.960.86, 1.08
  LOO model (iAs% fixed)1.020.89, 1.180.920.81, 1.050.960.84, 1.10
  LOO model (DMA% fixed)1.030.85, 1.250.920.75, 1.120.950.77, 1.17
 DMA%
  Conventional model0.990.93, 1.051.030.98, 1.091.010.96, 1.06
  LOO model (iAs% fixed)0.980.85, 1.121.080.95, 1.231.040.91, 1.19
  LOO model (MMA% fixed)1.000.92, 1.101.000.91, 1.100.990.89, 1.09
Elevated fasting plasma glucose concentration (n = 937)i
 ∑As, per IQR increasee1.481.28, 1.731.251.07, 1.461.191.01, 1.41
 iAs%
  Conventional model1.040.95, 1.131.010.92, 1.101.030.95, 1.12
  LOO model (MMA% fixed)1.070.98, 1.171.040.95, 1.131.030.94, 1.12
  LOO model (DMA% fixed)1.211.01, 1.441.160.97, 1.381.030.86, 1.23
 MMA%
  Conventional model0.920.82, 1.030.910.81, 1.031.010.90, 1.14
  LOO model (iAs% fixed)0.890.78, 1.000.900.79, 1.021.000.88, 1.14
  LOO model (DMA% fixed)0.830.69, 0.990.860.72, 1.030.970.82, 1.16
 DMA%
  Conventional model1.010.95, 1.071.020.96, 1.090.980.93, 1.04
  LOO model (iAs% fixed)1.131.00, 1.271.120.98, 1.271.000.88, 1.13
  LOO model (MMA% fixed)0.930.86, 1.020.970.88, 1.050.970.89, 1.06
Component, Arsenic Biomarker, and ModelModel 1aModel 2bModel 3c
RR95% CIRR95% CIRR95% CI
Elevated waist circumference (n = 574)d
 ∑As, per IQR increasee0.990.86, 1.140.990.86, 1.140.980.85, 1.13
 iAs%
  Conventional model0.960.88, 1.040.970.89, 1.060.970.89, 1.06
  LOO model (MMA% fixed)1.020.94, 1.121.020.94, 1.121.030.94, 1.12
  LOO model (DMA% fixed)1.291.10, 1.511.251.07, 1.451.251.08, 1.45
 MMA%
  Conventional model0.800.74, 0.880.830.76, 0.910.830.76, 0.91
  LOO model (iAs% fixed)0.800.72, 0.880.820.75, 0.910.820.75, 0.91
  LOO model (DMA% fixed)0.780.66, 0.910.800.69, 0.930.800.69, 0.93
 DMA%
  Conventional model1.081.03, 1.141.071.01, 1.131.071.01, 1.13
  LOO model (iAs% fixed)1.261.14, 1.391.221.10, 1.341.221.10, 1.34
  LOO model (MMA% fixed)0.980.89, 1.060.980.90, 1.060.970.89, 1.06
Elevated triglyceride concentration (n = 881)f
 ∑As, per IQR increasee0.960.81, 1.141.000.84, 1.190.980.82, 1.19
 iAs%
  Conventional model0.950.86, 1.040.930.84, 1.030.930.84, 1.03
  LOO model (MMA% fixed)0.960.87, 1.060.950.86, 1.060.940.84, 1.05
  LOO model (DMA% fixed)0.990.82, 1.181.010.84, 1.220.990.81, 1.19
 MMA%
  Conventional model0.950.85, 1.060.920.82, 1.030.930.83, 1.04
  LOO model (iAs% fixed)0.970.86, 1.090.940.83, 1.060.960.84, 1.08
  LOO model (DMA% fixed)1.010.84, 1.210.990.82, 1.191.010.84, 1.23
 DMA%
  Conventional model1.040.98, 1.101.060.99, 1.131.050.99, 1.13
  LOO model (iAs% fixed)1.030.92, 1.161.060.94, 1.201.050.92, 1.18
  LOO model (MMA% fixed)1.040.94,1.151.050.94, 1.171.060.95, 1.18
Reduced HDL cholesterol concentration (n = 792)g
 ∑As, per IQR increasee1.070.91, 1.250.960.81, 1.140.950.80, 1.14
 iAs%
  Conventional model0.890.81, 0.980.900.82, 0.990.920.84, 1.02
  LOO model (MMA% fixed)0.940.86, 1.040.930.85, 1.030.940.85, 1.04
  LOO model (DMA% fixed)1.110.93, 1.341.040.86, 1.251.000.83, 1.22
 MMA%
  Conventional model0.820.73, 0.930.870.77, 0.980.910.81, 1.03
  LOO model (iAs% fixed)0.850.75, 0.960.900.79, 1.020.940.82, 1.07
  LOO model (DMA% fixed)0.900.75, 1.080.960.80, 1.161.000.82, 1.21
 DMA%
  Conventional model1.111.04, 1.181.091.02, 1.161.061.00, 1.13
  LOO model (iAs% fixed)1.181.04, 1.341.120.98, 1.271.060.93, 1.22
  LOO model (MMA% fixed)1.060.96, 1.161.070.97, 1.181.060.96, 1.17
Hypertension (n = 851)h
 ∑As, per IQR increasee1.110.91, 1.361.050.90, 1.231.030.88, 1.20
 iAs%
  Conventional model1.000.93, 1.090.980.90, 1.061.000.92, 1.08
  LOO model (MMA% fixed)1.000.91, 1.091.000.91, 1.101.010.91, 1.12
  LOO model (DMA% fixed)0.970.80, 1.181.090.89, 1.331.050.85, 1.30
 MMA%
  Conventional model1.020.90, 1.160.920.83, 1.030.960.86, 1.08
  LOO model (iAs% fixed)1.020.89, 1.180.920.81, 1.050.960.84, 1.10
  LOO model (DMA% fixed)1.030.85, 1.250.920.75, 1.120.950.77, 1.17
 DMA%
  Conventional model0.990.93, 1.051.030.98, 1.091.010.96, 1.06
  LOO model (iAs% fixed)0.980.85, 1.121.080.95, 1.231.040.91, 1.19
  LOO model (MMA% fixed)1.000.92, 1.101.000.91, 1.100.990.89, 1.09
Elevated fasting plasma glucose concentration (n = 937)i
 ∑As, per IQR increasee1.481.28, 1.731.251.07, 1.461.191.01, 1.41
 iAs%
  Conventional model1.040.95, 1.131.010.92, 1.101.030.95, 1.12
  LOO model (MMA% fixed)1.070.98, 1.171.040.95, 1.131.030.94, 1.12
  LOO model (DMA% fixed)1.211.01, 1.441.160.97, 1.381.030.86, 1.23
 MMA%
  Conventional model0.920.82, 1.030.910.81, 1.031.010.90, 1.14
  LOO model (iAs% fixed)0.890.78, 1.000.900.79, 1.021.000.88, 1.14
  LOO model (DMA% fixed)0.830.69, 0.990.860.72, 1.030.970.82, 1.16
 DMA%
  Conventional model1.010.95, 1.071.020.96, 1.090.980.93, 1.04
  LOO model (iAs% fixed)1.131.00, 1.271.120.98, 1.271.000.88, 1.13
  LOO model (MMA% fixed)0.930.86, 1.020.970.88, 1.050.970.89, 1.06

Abbreviations: ∑As, sum of inorganic and methylated arsenic species; CI, confidence interval; DMA, dimethylarsinic acid; HDL, high-density lipoprotein; iAs, inorganic arsenic; IQR, interquartile range; LOO, leave-one-out; MMA, monomethylarsonic acid; RR, relative risk.

a Model 1 adjusted for log total arsenic (except in models where ∑As was the predictor of interest) and urinary creatinine concentration.

b Model 2 further adjusted for age, sex, region, and education.

c Model 3 further adjusted for alcohol intake, smoking status, kidney function, and body mass index (except in waist circumference models).

d Waist circumference ≥40 inches (≥102 cm) in men and ≥35 inches (≥89 cm) in women.

e 10.84 μg/L (75th percentile) vs. 4.20 μg/L (25th percentile).

f Triglyceride concentration ≥150 mg/dL (or use of medication).

g HDL cholesterol concentration <40 mg/dL for men and <50 mg/dL for women (or use of medication).

h Systolic blood pressure ≥130 mm Hg or diastolic blood pressure ≥85 mm Hg (or use of medication).

i Fasting plasma glucose concentration ≥100 mg/dL (or use of medication).

Table 5.

Relative Risk of Incident Components of the Metabolic Syndrome According to Increases in the Sum of Inorganic and Methylated Arsenic Species and Biomarkers of Arsenic Metabolism (iAs%, MMA%, and DMA%), Strong Heart Family Study, 1998–2009

Component, Arsenic Biomarker, and ModelModel 1aModel 2bModel 3c
RR95% CIRR95% CIRR95% CI
Elevated waist circumference (n = 574)d
 ∑As, per IQR increasee0.990.86, 1.140.990.86, 1.140.980.85, 1.13
 iAs%
  Conventional model0.960.88, 1.040.970.89, 1.060.970.89, 1.06
  LOO model (MMA% fixed)1.020.94, 1.121.020.94, 1.121.030.94, 1.12
  LOO model (DMA% fixed)1.291.10, 1.511.251.07, 1.451.251.08, 1.45
 MMA%
  Conventional model0.800.74, 0.880.830.76, 0.910.830.76, 0.91
  LOO model (iAs% fixed)0.800.72, 0.880.820.75, 0.910.820.75, 0.91
  LOO model (DMA% fixed)0.780.66, 0.910.800.69, 0.930.800.69, 0.93
 DMA%
  Conventional model1.081.03, 1.141.071.01, 1.131.071.01, 1.13
  LOO model (iAs% fixed)1.261.14, 1.391.221.10, 1.341.221.10, 1.34
  LOO model (MMA% fixed)0.980.89, 1.060.980.90, 1.060.970.89, 1.06
Elevated triglyceride concentration (n = 881)f
 ∑As, per IQR increasee0.960.81, 1.141.000.84, 1.190.980.82, 1.19
 iAs%
  Conventional model0.950.86, 1.040.930.84, 1.030.930.84, 1.03
  LOO model (MMA% fixed)0.960.87, 1.060.950.86, 1.060.940.84, 1.05
  LOO model (DMA% fixed)0.990.82, 1.181.010.84, 1.220.990.81, 1.19
 MMA%
  Conventional model0.950.85, 1.060.920.82, 1.030.930.83, 1.04
  LOO model (iAs% fixed)0.970.86, 1.090.940.83, 1.060.960.84, 1.08
  LOO model (DMA% fixed)1.010.84, 1.210.990.82, 1.191.010.84, 1.23
 DMA%
  Conventional model1.040.98, 1.101.060.99, 1.131.050.99, 1.13
  LOO model (iAs% fixed)1.030.92, 1.161.060.94, 1.201.050.92, 1.18
  LOO model (MMA% fixed)1.040.94,1.151.050.94, 1.171.060.95, 1.18
Reduced HDL cholesterol concentration (n = 792)g
 ∑As, per IQR increasee1.070.91, 1.250.960.81, 1.140.950.80, 1.14
 iAs%
  Conventional model0.890.81, 0.980.900.82, 0.990.920.84, 1.02
  LOO model (MMA% fixed)0.940.86, 1.040.930.85, 1.030.940.85, 1.04
  LOO model (DMA% fixed)1.110.93, 1.341.040.86, 1.251.000.83, 1.22
 MMA%
  Conventional model0.820.73, 0.930.870.77, 0.980.910.81, 1.03
  LOO model (iAs% fixed)0.850.75, 0.960.900.79, 1.020.940.82, 1.07
  LOO model (DMA% fixed)0.900.75, 1.080.960.80, 1.161.000.82, 1.21
 DMA%
  Conventional model1.111.04, 1.181.091.02, 1.161.061.00, 1.13
  LOO model (iAs% fixed)1.181.04, 1.341.120.98, 1.271.060.93, 1.22
  LOO model (MMA% fixed)1.060.96, 1.161.070.97, 1.181.060.96, 1.17
Hypertension (n = 851)h
 ∑As, per IQR increasee1.110.91, 1.361.050.90, 1.231.030.88, 1.20
 iAs%
  Conventional model1.000.93, 1.090.980.90, 1.061.000.92, 1.08
  LOO model (MMA% fixed)1.000.91, 1.091.000.91, 1.101.010.91, 1.12
  LOO model (DMA% fixed)0.970.80, 1.181.090.89, 1.331.050.85, 1.30
 MMA%
  Conventional model1.020.90, 1.160.920.83, 1.030.960.86, 1.08
  LOO model (iAs% fixed)1.020.89, 1.180.920.81, 1.050.960.84, 1.10
  LOO model (DMA% fixed)1.030.85, 1.250.920.75, 1.120.950.77, 1.17
 DMA%
  Conventional model0.990.93, 1.051.030.98, 1.091.010.96, 1.06
  LOO model (iAs% fixed)0.980.85, 1.121.080.95, 1.231.040.91, 1.19
  LOO model (MMA% fixed)1.000.92, 1.101.000.91, 1.100.990.89, 1.09
Elevated fasting plasma glucose concentration (n = 937)i
 ∑As, per IQR increasee1.481.28, 1.731.251.07, 1.461.191.01, 1.41
 iAs%
  Conventional model1.040.95, 1.131.010.92, 1.101.030.95, 1.12
  LOO model (MMA% fixed)1.070.98, 1.171.040.95, 1.131.030.94, 1.12
  LOO model (DMA% fixed)1.211.01, 1.441.160.97, 1.381.030.86, 1.23
 MMA%
  Conventional model0.920.82, 1.030.910.81, 1.031.010.90, 1.14
  LOO model (iAs% fixed)0.890.78, 1.000.900.79, 1.021.000.88, 1.14
  LOO model (DMA% fixed)0.830.69, 0.990.860.72, 1.030.970.82, 1.16
 DMA%
  Conventional model1.010.95, 1.071.020.96, 1.090.980.93, 1.04
  LOO model (iAs% fixed)1.131.00, 1.271.120.98, 1.271.000.88, 1.13
  LOO model (MMA% fixed)0.930.86, 1.020.970.88, 1.050.970.89, 1.06
Component, Arsenic Biomarker, and ModelModel 1aModel 2bModel 3c
RR95% CIRR95% CIRR95% CI
Elevated waist circumference (n = 574)d
 ∑As, per IQR increasee0.990.86, 1.140.990.86, 1.140.980.85, 1.13
 iAs%
  Conventional model0.960.88, 1.040.970.89, 1.060.970.89, 1.06
  LOO model (MMA% fixed)1.020.94, 1.121.020.94, 1.121.030.94, 1.12
  LOO model (DMA% fixed)1.291.10, 1.511.251.07, 1.451.251.08, 1.45
 MMA%
  Conventional model0.800.74, 0.880.830.76, 0.910.830.76, 0.91
  LOO model (iAs% fixed)0.800.72, 0.880.820.75, 0.910.820.75, 0.91
  LOO model (DMA% fixed)0.780.66, 0.910.800.69, 0.930.800.69, 0.93
 DMA%
  Conventional model1.081.03, 1.141.071.01, 1.131.071.01, 1.13
  LOO model (iAs% fixed)1.261.14, 1.391.221.10, 1.341.221.10, 1.34
  LOO model (MMA% fixed)0.980.89, 1.060.980.90, 1.060.970.89, 1.06
Elevated triglyceride concentration (n = 881)f
 ∑As, per IQR increasee0.960.81, 1.141.000.84, 1.190.980.82, 1.19
 iAs%
  Conventional model0.950.86, 1.040.930.84, 1.030.930.84, 1.03
  LOO model (MMA% fixed)0.960.87, 1.060.950.86, 1.060.940.84, 1.05
  LOO model (DMA% fixed)0.990.82, 1.181.010.84, 1.220.990.81, 1.19
 MMA%
  Conventional model0.950.85, 1.060.920.82, 1.030.930.83, 1.04
  LOO model (iAs% fixed)0.970.86, 1.090.940.83, 1.060.960.84, 1.08
  LOO model (DMA% fixed)1.010.84, 1.210.990.82, 1.191.010.84, 1.23
 DMA%
  Conventional model1.040.98, 1.101.060.99, 1.131.050.99, 1.13
  LOO model (iAs% fixed)1.030.92, 1.161.060.94, 1.201.050.92, 1.18
  LOO model (MMA% fixed)1.040.94,1.151.050.94, 1.171.060.95, 1.18
Reduced HDL cholesterol concentration (n = 792)g
 ∑As, per IQR increasee1.070.91, 1.250.960.81, 1.140.950.80, 1.14
 iAs%
  Conventional model0.890.81, 0.980.900.82, 0.990.920.84, 1.02
  LOO model (MMA% fixed)0.940.86, 1.040.930.85, 1.030.940.85, 1.04
  LOO model (DMA% fixed)1.110.93, 1.341.040.86, 1.251.000.83, 1.22
 MMA%
  Conventional model0.820.73, 0.930.870.77, 0.980.910.81, 1.03
  LOO model (iAs% fixed)0.850.75, 0.960.900.79, 1.020.940.82, 1.07
  LOO model (DMA% fixed)0.900.75, 1.080.960.80, 1.161.000.82, 1.21
 DMA%
  Conventional model1.111.04, 1.181.091.02, 1.161.061.00, 1.13
  LOO model (iAs% fixed)1.181.04, 1.341.120.98, 1.271.060.93, 1.22
  LOO model (MMA% fixed)1.060.96, 1.161.070.97, 1.181.060.96, 1.17
Hypertension (n = 851)h
 ∑As, per IQR increasee1.110.91, 1.361.050.90, 1.231.030.88, 1.20
 iAs%
  Conventional model1.000.93, 1.090.980.90, 1.061.000.92, 1.08
  LOO model (MMA% fixed)1.000.91, 1.091.000.91, 1.101.010.91, 1.12
  LOO model (DMA% fixed)0.970.80, 1.181.090.89, 1.331.050.85, 1.30
 MMA%
  Conventional model1.020.90, 1.160.920.83, 1.030.960.86, 1.08
  LOO model (iAs% fixed)1.020.89, 1.180.920.81, 1.050.960.84, 1.10
  LOO model (DMA% fixed)1.030.85, 1.250.920.75, 1.120.950.77, 1.17
 DMA%
  Conventional model0.990.93, 1.051.030.98, 1.091.010.96, 1.06
  LOO model (iAs% fixed)0.980.85, 1.121.080.95, 1.231.040.91, 1.19
  LOO model (MMA% fixed)1.000.92, 1.101.000.91, 1.100.990.89, 1.09
Elevated fasting plasma glucose concentration (n = 937)i
 ∑As, per IQR increasee1.481.28, 1.731.251.07, 1.461.191.01, 1.41
 iAs%
  Conventional model1.040.95, 1.131.010.92, 1.101.030.95, 1.12
  LOO model (MMA% fixed)1.070.98, 1.171.040.95, 1.131.030.94, 1.12
  LOO model (DMA% fixed)1.211.01, 1.441.160.97, 1.381.030.86, 1.23
 MMA%
  Conventional model0.920.82, 1.030.910.81, 1.031.010.90, 1.14
  LOO model (iAs% fixed)0.890.78, 1.000.900.79, 1.021.000.88, 1.14
  LOO model (DMA% fixed)0.830.69, 0.990.860.72, 1.030.970.82, 1.16
 DMA%
  Conventional model1.010.95, 1.071.020.96, 1.090.980.93, 1.04
  LOO model (iAs% fixed)1.131.00, 1.271.120.98, 1.271.000.88, 1.13
  LOO model (MMA% fixed)0.930.86, 1.020.970.88, 1.050.970.89, 1.06

Abbreviations: ∑As, sum of inorganic and methylated arsenic species; CI, confidence interval; DMA, dimethylarsinic acid; HDL, high-density lipoprotein; iAs, inorganic arsenic; IQR, interquartile range; LOO, leave-one-out; MMA, monomethylarsonic acid; RR, relative risk.

a Model 1 adjusted for log total arsenic (except in models where ∑As was the predictor of interest) and urinary creatinine concentration.

b Model 2 further adjusted for age, sex, region, and education.

c Model 3 further adjusted for alcohol intake, smoking status, kidney function, and body mass index (except in waist circumference models).

d Waist circumference ≥40 inches (≥102 cm) in men and ≥35 inches (≥89 cm) in women.

e 10.84 μg/L (75th percentile) vs. 4.20 μg/L (25th percentile).

f Triglyceride concentration ≥150 mg/dL (or use of medication).

g HDL cholesterol concentration <40 mg/dL for men and <50 mg/dL for women (or use of medication).

h Systolic blood pressure ≥130 mm Hg or diastolic blood pressure ≥85 mm Hg (or use of medication).

i Fasting plasma glucose concentration ≥100 mg/dL (or use of medication).

Relative risk of the metabolic syndrome and its individual components according to the sum of inorganic and methylated arsenic species (∑As), Strong Heart Family Study, 1998–2009. ∑As was calculated as the sum of urinary inorganic arsenic, monomethylarsonic acid, and dimethylarsinic acid concentrations. The solid black lines represent adjusted relative risks of: incident metabolic syndrome, defined as meeting 3 or more of the individual criteria (A); elevated waist circumference, defined as ≥40 inches (≥102 cm) in men and ≥35 inches (≥89 cm) in women (B); elevated triglyceride level, defined as ≥150 mg/dL (or use of medication) (C); reduced high-density lipoprotein cholesterol level, defined as <40 mg/dL for men and <50 mg/dL for women (or use of medication) (D); hypertension, defined as systolic blood pressure ≥130 mm Hg or diastolic blood pressure ≥85 mm Hg (or use of medication) (E); and elevated fasting plasma glucose level, defined as ≥100 mg/dL (or use of medication) (F). Relative risks were calculated using restricted cubic splines for ∑As with knots at the 10th, 50th, and 90th percentiles of the ∑As distribution and were adjusted for log urinary creatinine, age, sex, body mass index (excluding waist circumference models), educational level, region (Arizona, Oklahoma, or North/South Dakota), smoking status, alcohol intake, and kidney function. The referent was set at the 10th percentile of the ∑As distribution. Dotted lines show 95% confidence intervals. Histograms represent the distribution of log-transformed ∑As values.
Figure 2.

Relative risk of the metabolic syndrome and its individual components according to the sum of inorganic and methylated arsenic species (∑As), Strong Heart Family Study, 1998–2009. ∑As was calculated as the sum of urinary inorganic arsenic, monomethylarsonic acid, and dimethylarsinic acid concentrations. The solid black lines represent adjusted relative risks of: incident metabolic syndrome, defined as meeting 3 or more of the individual criteria (A); elevated waist circumference, defined as ≥40 inches (≥102 cm) in men and ≥35 inches (≥89 cm) in women (B); elevated triglyceride level, defined as ≥150 mg/dL (or use of medication) (C); reduced high-density lipoprotein cholesterol level, defined as <40 mg/dL for men and <50 mg/dL for women (or use of medication) (D); hypertension, defined as systolic blood pressure ≥130 mm Hg or diastolic blood pressure ≥85 mm Hg (or use of medication) (E); and elevated fasting plasma glucose level, defined as ≥100 mg/dL (or use of medication) (F). Relative risks were calculated using restricted cubic splines for ∑As with knots at the 10th, 50th, and 90th percentiles of the ∑As distribution and were adjusted for log urinary creatinine, age, sex, body mass index (excluding waist circumference models), educational level, region (Arizona, Oklahoma, or North/South Dakota), smoking status, alcohol intake, and kidney function. The referent was set at the 10th percentile of the ∑As distribution. Dotted lines show 95% confidence intervals. Histograms represent the distribution of log-transformed ∑As values.

For arsenic metabolism, the strongest association with individual MetS components was for waist circumference (Figure 3, Table 5). In the conventional model for MMA%, the relative risk of elevated waist circumference for a 5% increase in MMA% was 0.83 (95% CI: 0.76, 0.91), with consistent associations in both leave-one-out models. For DMA%, the relative risk of elevated waist circumference increased from 1.07 (95% CI: 1.01, 1.13) for a 5% increase in DMA% in the conventional model to 1.22 (95% CI: 1.10, 1.34) in the iAs%-fixed leave-one-out model. For iAs%, the only association with waist circumference was in the DMA%-fixed leave-one-out model (RR = 1.25, 95% CI: 1.08, 1.45). In the waist circumference sensitivity analysis that additionally adjusted for baseline waist circumference, results were consistent (data not shown). Clear associations were not observed between arsenic metabolism and other MetS components beyond waist circumference, although in conventional models, generally, higher DMA% and lower MMA% were nonsignificantly associated with higher risk of elevated triglyceride levels, reduced high-density lipoprotein cholesterol level, and hypertension (Figure 3, Table 5).

Relative risk of the metabolic syndrome (MetS) and its individual components according to the sum of inorganic and methylated arsenic species (∑As) and per 5% increase in biomarkers of arsenic metabolism, Strong Heart Family Study, 1998–2009. ∑As was calculated as the sum of urinary inorganic arsenic (iAs), monomethylarsonic acid (MMA), and dimethylarsinic acid (DMA) concentrations. Part A shows the relative risk of incident MetS and its individual components (elevated waist circumference (WC; ≥40 inches (≥102 cm) in men and ≥35 inches (≥89 cm) in women), elevated triglyceride (TG) level (≥150 mg/dL (or use of medication)), reduced high-density lipoprotein (HDL) cholesterol level (<40 mg/dL for men and <50 mg/dL for women (or use of medication)), hypertension (HT; systolic blood pressure ≥130 mm Hg or diastolic blood pressure ≥85 mm Hg (or use of medication)), and elevated fasting plasma glucose (FPG) level (≥100 mg/dL (or use of medication)) per interquartile-range increase in ∑As. Points (●) represent relative risk, and vertical lines represent 95% confidence intervals. Parts B–D show the relative risk of MetS and its individual components per 5% increase in percentage of iAs (B), percentage of MMA (C), and percentage of DMA (D). The 3 different types of points in panels B–D reflect the 3 different models used to estimate relative risk: the conventional model (■) and 2 leave-one-out (LOO) models (● and ▲). Relative risks were adjusted for log ∑As (excluding exposure models), log urinary creatinine, age, sex, body mass index (excluding waist circumference models), educational level, region (Arizona, Oklahoma, or North/South Dakota), smoking status, alcohol intake, and kidney function.
Figure 3.

Relative risk of the metabolic syndrome (MetS) and its individual components according to the sum of inorganic and methylated arsenic species (∑As) and per 5% increase in biomarkers of arsenic metabolism, Strong Heart Family Study, 1998–2009. ∑As was calculated as the sum of urinary inorganic arsenic (iAs), monomethylarsonic acid (MMA), and dimethylarsinic acid (DMA) concentrations. Part A shows the relative risk of incident MetS and its individual components (elevated waist circumference (WC; ≥40 inches (≥102 cm) in men and ≥35 inches (≥89 cm) in women), elevated triglyceride (TG) level (≥150 mg/dL (or use of medication)), reduced high-density lipoprotein (HDL) cholesterol level (<40 mg/dL for men and <50 mg/dL for women (or use of medication)), hypertension (HT; systolic blood pressure ≥130 mm Hg or diastolic blood pressure ≥85 mm Hg (or use of medication)), and elevated fasting plasma glucose (FPG) level (≥100 mg/dL (or use of medication)) per interquartile-range increase in ∑As. Points (●) represent relative risk, and vertical lines represent 95% confidence intervals. Parts B–D show the relative risk of MetS and its individual components per 5% increase in percentage of iAs (B), percentage of MMA (C), and percentage of DMA (D). The 3 different types of points in panels B–D reflect the 3 different models used to estimate relative risk: the conventional model (■) and 2 leave-one-out (LOO) models (● and ▲). Relative risks were adjusted for log ∑As (excluding exposure models), log urinary creatinine, age, sex, body mass index (excluding waist circumference models), educational level, region (Arizona, Oklahoma, or North/South Dakota), smoking status, alcohol intake, and kidney function.

DISCUSSION

In American Indian men and women aged ≥14 years from Arizona, Oklahoma, and North/South Dakota, arsenic exposure was associated with increased risk of elevated fasting plasma glucose level but not with MetS or other individual components of MetS. Arsenic metabolism patterns, independent of arsenic exposure, were associated with both incident MetS and elevated waist circumference but not with other components of the syndrome. The relative percentage of MMA appeared to be the main driver behind these associations. The distinct and independent associations of arsenic exposure and arsenic metabolism with MetS and its individual components suggest that these are unrelated phenomena that could be contributing to overall diabetes risk at low levels of arsenic exposure. For arsenic exposure, the association appears to be predominately with hyperglycemia. For arsenic metabolism, the association appears to be mainly with central adiposity. The possible role played by arsenic exposure in hyperglycemia and by arsenic metabolism in central adiposity could be underlying mechanisms for the observed associations between arsenic and diabetes in multiple populations.

MetS has been identified as the driving force behind the global epidemics of diabetes and cardiovascular disease (2). The incidence of metabolic conditions continues to rise, with projections for global diabetes incidence alone to double by 2025 (2). Similar trends are evident in the United States (1, 52). While the National Health and Nutrition Examination Survey does not monitor trends in American Indians, Russell et al. (45) reported a MetS prevalence of up to 63% in an older American Indian cohort (ages 45–74 years). American Indians/Alaska Natives are reported to have the highest age-adjusted prevalence of diabetes in the United States (53). To our knowledge, our study was the first to evaluate the incidence of MetS from adolescence onwards among American Indians. With a median age of 30 years, our baseline MetS prevalence of 36% and incidence of 32% over a median of 5.3 years of follow-up is undoubtedly high. The elevated risks of diabetes and MetS among American Indians as compared with other US populations may be influenced in part by genetic background; however, environmental factors, including disproportionate exposure to arsenic in drinking water, have also been implicated (5458). Our study supports the hypothesis that exposure to environmental contaminants, particularly arsenic, could play a role.

The association between higher ∑As and elevated fasting plasma glucose level supports previous experimental and epidemiologic evidence suggesting that iAs may have a diabetogenic effect. Although the mechanisms are not fully understood, experimental studies have shown that arsenic may induce diabetes by inhibiting insulin signaling and insulin-dependent glucose uptake, with both microRNAs and mitochondria being proposed as potential mechanistic links (21, 22, 59, 60). Epidemiologic studies have found associations between arsenic and diabetes at both high (≥100 μg/L) (22, 6164) and moderate (<100 μg/L) (1618, 38, 6568) exposure levels, although findings for low levels are mixed (37, 6971).

The association between arsenic exposure and MetS, as well as other components of MetS beyond fasting plasma glucose, has been less studied. In contrast to our null findings on the relationship between ∑As and MetS, 2 studies conducted in highly exposed Taiwanese populations found a positive association (11, 12). Similarly, studies conducted in high-exposure regions suggest that arsenic exposure is associated with increased risk of hypertension (14, 20, 32, 72, 73), with mostly null findings in low-exposure populations (16, 19, 74), which is consistent with our results for hypertension. A few studies have suggested an association between arsenic exposure and triglycerides (16, 75) and high-density lipoprotein cholesterol (76, 77), in contrast to our null findings. However, epidemiologic studies evaluating low-to-moderate arsenic and triglyceride levels, high-density lipoprotein cholesterol, or hypertension have been cross-sectional. Finally, the association between arsenic exposure and waist circumference has not been evaluated before. In our study, we found no association.

For arsenic metabolism, our finding that lower MMA% and higher DMA% are prospectively associated with both increased risk of incident MetS and elevated waist circumference is consistent with studies that have evaluated these relationships (11, 35) and related outcomes such as diabetes (16, 37, 39, 68) and BMI (35, 78, 79), although most have been cross-sectional. To our knowledge, our evaluation of arsenic metabolism and MetS in a population exposed to low-to-moderate levels of arsenic is novel, as is our prospective assessment of the association between arsenic metabolism and waist circumference. Further, the results from our leave-one-out models add to prospective evidence showing that greater DMA% due to reductions in MMA% is related to incident diabetes (37).

Mechanistic understanding of the increased risk reported for metabolic outcomes with lower MMA% is limited. MMA and DMA in this study, as in most epidemiologic studies, include pentavalent and trivalent forms. Trivalent forms oxidize rapidly to pentavalent forms in urine, making them indistinguishable (80). Trivalent methylated arsenicals are considered the most toxic metabolites, and MMAIII has been proposed as a mechanism for the associations between higher MMA% and increased risk of cancer, cardiovascular disease, and skin lesions (21, 26, 34). In turn, some studies have suggested that the association between higher DMA% and increased risk of metabolic outcomes may be due to DMAIII (16, 37, 68). Other studies have suggested that confounding by diet may explain these associations. For example, a nutritionally sufficient diet with high fat and protein intake may enhance arsenic metabolism efficiency while also having adverse metabolic effects (39). In our study, however, adjusting for percentage of kilocalories derived from protein and fat did not affect results. Diet has also been suggested to play a role through 1-carbon metabolism (37). One-carbon metabolism, regulated by B vitamins, is responsible for providing methyl groups necessary for methylation reactions in arsenic metabolism and could be driving the association between arsenic metabolism and metabolic outcomes (25). However, we saw no evidence of confounding or effect modification by estimated intake of B vitamins (folate, vitamin B2, and vitamin B6).

In pregnant women, DMA% increases as women progress through pregnancy (81, 82). Because adiposity increases during pregnancy, it has been suggested that adiposity may be driving shifts in arsenic metabolism patterns (35). Our finding that waist circumference had the strongest associations with arsenic metabolism supports the possibility of adiposity’s playing a key role in the pathway between arsenic metabolism and metabolic outcomes. Our association was prospective, and it remained after adjustment for baseline waist circumference, suggesting that adiposity may be a consequence of arsenic metabolism, not a cause; still, we cannot discard the possibility that reverse causality explains the association. Notably, however, we saw significant interactions of DMA% and MetS with BMI, with an increasing strength of association with lower BMI. This could be interpreted as evidence against reverse causality, as we are seeing an association between higher DMA% and MetS among persons with normal BMI at baseline. In a previous study in this cohort, participants with the AS3MT rs12768205 AA genotype had lower MMA% (48). We observed an interaction between MMA% and AS3MT genotype, with risk of MetS being further reduced with increasing MMA% with each additional A allele, suggesting that the inverse association between MMA% and MetS is stronger among persons with genetically lower MMA%.

Our study strengths included a broad age distribution, a well-established cohort, high-quality laboratory methods, and consideration of genetic contributions. As in most epidemiologic studies, we were unable to differentiate between trivalent and pentavalent forms of MMA and DMA. We confirmed that seafood intake in our population was infrequent and adjusted for intake of B vitamins; however, we cannot discount the potential for remaining confounding by food sources. Further, we did not have information on vitamin B12, an important vitamin in the 1-carbon metabolism pathway. Because we did not have repeated measurements of arsenic exposure and metabolism, we could not evaluate how changes in adiposity and fasting glucose were related to changes in arsenic exposure over time; therefore, we cannot totally exclude the possibility of reverse causality, despite the study’s prospective design.

In conclusion, our findings support previous evidence that arsenic has diabetogenic effects even at low-level exposures and may be contributing, in part, to the high burden of diabetes in arsenic-exposed populations. Further, our results suggest that disruption of glucose regulation maybe a key pathway driving the association between low-level arsenic exposure and cardiometabolic outcomes. In addition, we found that arsenic metabolism patterns, specifically lower MMA%, are prospectively associated with increased risk for the development of MetS and elevated waist circumference. These findings support the importance of preventing low-level arsenic exposure and the need to better elucidate mechanisms underlying the contrasting and independent associations of arsenic exposure and metabolism with MetS.

ACKNOWLEDGMENTS

Author affiliations: Department of Environmental Health and Engineering, Bloomberg School of Public Health, Johns Hopkins University, Baltimore, Maryland (Miranda J. Spratlen, Ana Navas-Acien); Department of Environmental Health Sciences, Mailman School of Public Health, Columbia University, New York, New York (Miranda J. Spratlen, Poojitha Balakrishnan, Mary V. Gamble, Ana Navas-Acien); Area of Cardiometabolic and Renal Risk, Clinical Research Foundation of Valencia, Valencia, Spain (Maria Grau-Perez); Department of Statistics and Operational Research, Faculty of Mathematics, University of Valencia, Valencia, Spain (Maria Grau-Perez); Missouri Breaks Industries Research, Inc., Eagle Butte, South Dakota (Lyle G. Best, Joseph Yracheta); Department of Medicine, School of Medicine, Johns Hopkins University, Baltimore, Maryland (Mariana Lazo, Dhananjay Vaidya); Department of Epidemiology, Bloomberg School of Public Health, Johns Hopkins University, Baltimore, Maryland (Mariana Lazo, Dhananjay Vaidya); Department of Analytical Chemistry, Institute of Chemistry, University of Graz, Graz, Austria (Kevin A. Francesconi, Walter Goessler); Texas Biomedical Research Institute, San Antonio, Texas (Shelley A. Cole); MedStar Health Research Institute, Hyattsville, Maryland (Jason G. Umans, Barbara V. Howard); and Department of Medicine, School of Medicine, Georgetown University, Washington, DC (Jason G. Umans, Barbara V. Howard).

This work was funded by National Institute of Environmental Health Sciences grants 1F31ES027796-01, 5T32ES007141-33, R01ES025216, 5P30ES009089, and P42ES010349 and National Heart, Lung, and Blood Institute grants R01-HL090863, R01-HL109315, R01HL109301, R01HL109284, R01HL109282, R01HL109319, U01-HL41642, U01-HL41652, U01-HL41654, U01-HL65520, and U01-HL65521.

Conflict of interest: none declared.

Abbreviations

     
  • ∑As

    sum of inorganic and methylated arsenic species

  •  
  • AS3MT

    arsenite methyltransferase gene

  •  
  • BMI

    body mass index

  •  
  • CI

    confidence interval

  •  
  • DMA

    dimethylarsinic acid

  •  
  • iAs

    inorganic arsenic

  •  
  • IQR

    interquartile range

  •  
  • MetS

    metabolic syndrome

  •  
  • MMA

    monomethylarsonic acid

  •  
  • PC

    principal component

  •  
  • RR

    relative risk

REFERENCES

1

Aguilar
M
,
Bhuket
T
,
Torres
S
, et al. .

Prevalence of the metabolic syndrome in the United States, 2003–2012
.
JAMA
.
2015
;
313
(
19
):
1973
1974
.
2

Alberti
GZ
,
Shaw
J
,
Grundy
SM
. The IDF Consensus Worldwide Definition of the Metabolic Syndrome.

Brussels, Belgium
:
International Diabetes Federation
;
2006
. https://www.idf.org/e-library/consensus-statements/60-idfconsensus-worldwide-definitionof-the-metabolic-syndrome. Accessed May 2, 2018.
3

Grundy
SM
,
Cleeman
JI
,
Daniels
SR
, et al. .

Diagnosis and management of the metabolic syndrome: an American Heart Association/National Heart, Lung, and Blood Institute Scientific Statement
.
Circulation
.
2005
;
112
(
17
):
2735
2752
.
4

Stern
MP
,
Williams
K
,
González-Villalpando
C
, et al. .

Does the metabolic syndrome improve identification of individuals at risk of type 2 diabetes and/or cardiovascular disease?
Diabetes Care
.
2004
;
27
(
11
):
2676
2681
.
5

Paul
DS
,
Hernández-Zavala
A
,
Walton
FS
, et al. .

Examination of the effects of arsenic on glucose homeostasis in cell culture and animal studies: development of a mouse model for arsenic-induced diabetes
.
Toxicol Appl Pharmacol
.
2007
;
222
(
3
):
305
314
.
6

Kuo
CC
,
Moon
K
,
Thayer
KA
, et al. .

Environmental chemicals and type 2 diabetes: an updated systematic review of the epidemiologic evidence
.
Curr Diab Rep
.
2013
;
13
(
6
):
831
849
.
7

Moon
K
,
Guallar
E
,
Navas-Acien
A
.

Arsenic exposure and cardiovascular disease: an updated systematic review
.
Curr Atheroscler Rep
.
2012
;
14
(
6
):
542
555
.
8

Moon
KA
,
Guallar
E
,
Umans
JG
, et al. .

Association between exposure to low to moderate arsenic levels and incident cardiovascular disease: a prospective cohort study
.
Ann Intern Med
.
2013
;
159
(
10
):
649
659
.
9

Tsuji
JS
,
Perez
V
,
Garry
MR
, et al. .

Association of low-level arsenic exposure in drinking water with cardiovascular disease: a systematic review and risk assessment
.
Toxicology
.
2014
;
323
:
78
94
.
10

Sung
TC
,
Huang
JW
,
Guo
HR
.

Association between arsenic exposure and diabetes: a meta-analysis
.
Biomed Res Int
.
2015
;
2015
:
368087
.
11

Chen
JW
,
Wang
SL
,
Wang
YH
, et al. .

Arsenic methylation, GSTO1 polymorphisms, and metabolic syndrome in an arseniasis endemic area of southwestern Taiwan
.
Chemosphere
.
2012
;
88
(
4
):
432
438
.
12

Wang
SL
,
Chang
FH
,
Liou
SH
, et al. .

Inorganic arsenic exposure and its relation to metabolic syndrome in an industrial area of Taiwan
.
Environ Int
.
2007
;
33
(
6
):
805
811
.
13

Wang
W
,
Xie
Z
,
Lin
Y
, et al. .

Association of inorganic arsenic exposure with type 2 diabetes mellitus: a meta-analysis
.
J Epidemiol Community Health
.
2014
;
68
(
2
):
176
184
.
14

Wang
SL
,
Li
WF
,
Chen
CJ
, et al. .

Hypertension incidence after tap-water implementation: a 13-year follow-up study in the arseniasis-endemic area of southwestern Taiwan
.
Sci Total Environ
.
2011
;
409
(
21
):
4528
4535
.
15

Navas-Acien
A
,
Silbergeld
EK
,
Pastor-Barriuso
R
, et al. .

Arsenic exposure and prevalence of type 2 diabetes in US adults
.
JAMA
.
2008
;
300
(
7
):
814
822
.
16

Mendez
MA
,
González-Horta
C
,
Sánchez-Ramírez
B
, et al. .

Chronic exposure to arsenic and markers of cardiometabolic risk: a cross-sectional study in Chihuahua, Mexico
.
Environ Health Perspect
.
2016
;
124
(
1
):
104
111
.
17

Kim
NH
,
Mason
CC
,
Nelson
RG
, et al. .

Arsenic exposure and incidence of type 2 diabetes in southwestern American Indians
.
Am J Epidemiol
.
2013
;
177
(
9
):
962
969
.
18

James
KA
,
Marshall
JA
,
Hokanson
JE
, et al. .

A case-cohort study examining lifetime exposure to inorganic arsenic in drinking water and diabetes mellitus
.
Environ Res
.
2013
;
123
:
33
38
.
19

Li
X
,
Li
B
,
Xi
S
, et al. .

Prolonged environmental exposure of arsenic through drinking water on the risk of hypertension and type 2 diabetes
.
Environ Sci Pollut Res Int
.
2013
;
20
(
11
):
8151
8161
.
20

Jiang
J
,
Liu
M
,
Parvez
F
, et al. .

Association between arsenic exposure from drinking water and longitudinal change in blood pressure among HEALS cohort participants
.
Environ Health Perspect
.
2015
;
123
(
8
):
806
812
.
21

Douillet
C
,
Currier
J
,
Saunders
J
, et al. .

Methylated trivalent arsenicals are potent inhibitors of glucose stimulated insulin secretion by murine pancreatic islets
.
Toxicol Appl Pharmacol
.
2013
;
267
(
1
):
11
15
.
22

Maull
EA
,
Ahsan
H
,
Edwards
J
, et al. .

Evaluation of the association between arsenic and diabetes: a National Toxicology Program workshop review
.
Environ Health Perspect
.
2012
;
120
(
12
):
1658
1670
.
23

Chen
JW
,
Chen
HY
,
Li
WF
, et al. .

The association between total urinary arsenic concentration and renal dysfunction in a community-based population from central Taiwan
.
Chemosphere
.
2011
;
84
(
1
):
17
24
.
24

Hall
MN
,
Liu
X
,
Slavkovich
V
, et al. .

Influence of cobalamin on arsenic metabolism in Bangladesh
.
Environ Health Perspect
.
2009
;
117
(
11
):
1724
1729
.
25

Hall
MN
,
Gamble
MV
.

Nutritional manipulation of one-carbon metabolism: effects on arsenic methylation and toxicity
.
J Toxicol
.
2012
;
2012
:
595307
.
26

Steinmaus
C
,
Bates
MN
,
Yuan
Y
, et al. .

Arsenic methylation and bladder cancer risk in case-control studies in Argentina and the United States
.
J Occup Environ Med
.
2006
;
48
(
5
):
478
488
.
27

Agusa
T
,
Kunito
T
,
Kubota
R
, et al. .

Exposure, metabolism, and health effects of arsenic in residents from arsenic-contaminated groundwater areas of Vietnam and Cambodia: a review
.
Rev Environ Health
.
2010
;
25
(
3
):
193
220
.
28

Chen
YC
,
Guo
YL
,
Su
HJ
, et al. .

Arsenic methylation and skin cancer risk in southwestern Taiwan
.
J Occup Environ Med
.
2003
;
45
(
3
):
241
248
.
29

Yu
RC
,
Hsu
KH
,
Chen
CJ
, et al. .

Arsenic methylation capacity and skin cancer
.
Cancer Epidemiol Biomarkers Prev
.
2000
;
9
(
11
):
1259
1262
.
30

Huang
YK
,
Huang
YL
,
Hsueh
YM
, et al. .

Arsenic exposure, urinary arsenic speciation, and the incidence of urothelial carcinoma: a twelve-year follow-up study
.
Cancer Causes Control
.
2008
;
19
(
8
):
829
839
.
31

Steinmaus
C
,
Yuan
Y
,
Kalman
D
, et al. .

Individual differences in arsenic metabolism and lung cancer in a case-control study in Cordoba, Argentina
.
Toxicol Appl Pharmacol
.
2010
;
247
(
2
):
138
145
.
32

Huang
YK
,
Tseng
CH
,
Huang
YL
, et al. .

Arsenic methylation capability and hypertension risk in subjects living in arseniasis-hyperendemic areas in southwestern Taiwan
.
Toxicol Appl Pharmacol
.
2007
;
218
(
2
):
135
142
.
33

Huang
YL
,
Hsueh
YM
,
Huang
YK
, et al. .

Urinary arsenic methylation capability and carotid atherosclerosis risk in subjects living in arsenicosis-hyperendemic areas in southwestern Taiwan
.
Sci Total Environ
.
2009
;
407
(
8
):
2608
2614
.
34

Chen
Y
,
Wu
F
,
Liu
M
, et al. .

A prospective study of arsenic exposure, arsenic methylation capacity, and risk of cardiovascular disease in Bangladesh
.
Environ Health Perspect
.
2013
;
121
(
7
):
832
838
.
35

Gribble
MO
,
Crainiceanu
CM
,
Howard
BV
, et al. .

Body composition and arsenic metabolism: a cross-sectional analysis in the Strong Heart Study
.
Environ Health
.
2013
;
12
:
107
.
36

Grau-Perez
M
,
Kuo
CC
,
Gribble
MO
, et al. .

Association of low-moderate arsenic exposure and arsenic metabolism with incident diabetes and insulin resistance in the Strong Heart Family Study
.
Environ Health Perspect
.
2017
;
125
(
12
):
127004
.
37

Kuo
CC
,
Howard
BV
,
Umans
JG
, et al. .

Arsenic exposure, arsenic metabolism, and incident diabetes in the Strong Heart Study
.
Diabetes Care
.
2015
;
38
(
4
):
620
627
.
38

Navas-Acien
A
,
Silbergeld
EK
,
Pastor-Barriuso
R
, et al. .

Arsenic exposure and prevalence of type 2 diabetes in US adults
.
JAMA
.
2008
;
300
(
7
):
814
822
.
39

Nizam
S
,
Kato
M
,
Yatsuya
H
, et al. .

Differences in urinary arsenic metabolites between diabetic and non-diabetic subjects in Bangladesh
.
Int J Environ Res Public Health
.
2013
;
10
(
3
):
1006
1019
.
40

Kuo
CC
,
Moon
KA
,
Wang
SL
, et al. .

The association of arsenic metabolism with cancer, cardiovascular disease, and diabetes: a systematic review of the epidemiological evidence
.
Environ Health Perspect
.
2017
;
125
(
8
):
087001
.
41

Yurgalevitch
SM
,
Kriska
AM
,
Welty
TK
, et al. .

Physical activity and lipids and lipoproteins in American Indians ages 45–74
.
Med Sci Sports Exerc
.
1998
;
30
(
4
):
543
549
.
42

North
KE
,
Howard
BV
,
Welty
TK
, et al. .

Genetic and environmental contributions to cardiovascular disease risk in American Indians: the Strong Heart Family Study
.
Am J Epidemiol
.
2003
;
157
(
4
):
303
314
.
43

Scheer
J
,
Findenig
S
,
Goessler
W
, et al. .

Arsenic species and selected metals in human urine: validation of HPLC/ICPMS and ICPMS procedures for a long-term population-based epidemiological study
.
Anal Methods
.
2012
;
4
(
2
):
406
413
.
44

de Simone
G
,
Devereux
RB
,
Chinali
M
, et al. .

Prognostic impact of metabolic syndrome by different definitions in a population with high prevalence of obesity and diabetes: the Strong Heart Study
.
Diabetes Care
.
2007
;
30
(
7
):
1851
1856
.
45

Russell
M
,
de Simone
G
,
Resnick
HE
, et al. .

The metabolic syndrome in American Indians: the Strong Heart Study
.
J Cardiometab Syndr
.
2007
;
2
(
4
):
283
287
.
46

Fretts
AM
,
Howard
BV
,
McKnight
B
, et al. .

Associations of processed meat and unprocessed red meat intake with incident diabetes: the Strong Heart Family Study
.
Am J Clin Nutr
.
2012
;
95
(
3
):
752
758
.
47

Shara
NM
,
Wang
H
,
Mete
M
, et al. .

Estimated GFR and incident cardiovascular disease events in American Indians: the Strong Heart Study
.
Am J Kidney Dis
.
2012
;
60
(
5
):
795
803
.
48

Balakrishnan
P
,
Vaidya
D
,
Franceschini
N
, et al. .

Association of cardiometabolic genes with arsenic metabolism biomarkers in American Indian communities: the Strong Heart Family Study (SHFS)
.
Environ Health Perspect
.
2017
;
125
(
1
):
15
22
.
49

Zou
G
.

A modified Poisson regression approach to prospective studies with binary data
.
Am J Epidemiol
.
2004
;
159
(
7
):
702
706
.
50

Jansen
RJ
,
Argos
M
,
Tong
L
, et al. .

Determinants and consequences of arsenic metabolism efficiency among 4,794 individuals: demographics, lifestyle, genetics, and toxicity
.
Cancer Epidemiol Biomarkers Prev
.
2016
;
25
(
2
):
381
390
.
51

Holvoet
P
,
Lee
DH
,
Steffes
M
, et al. .

Association between circulating oxidized low-density lipoprotein and incidence of the metabolic syndrome
.
JAMA
.
2008
;
299
(
19
):
2287
2293
.
52

Ford
ES
,
Giles
WH
,
Dietz
WH
.

Prevalence of the metabolic syndrome among US adults: findings from the Third National Health and Nutrition Examination Survey
.
JAMA
.
2002
;
287
(
3
):
356
359
.
53

American Diabetes Association
. Statistics about diabetes. Overall numbers, diabetes and prediabetes.

2016
. http://www.diabetes.org/diabetes-basics/statistics/?referrer=https://www.google.com/. Updated December 12, 2016. Accessed February 5, 2018.
54

Aminov
Z
,
Haase
R
,
Carpenter
DO
.

Diabetes in Native Americans: elevated risk as a result of exposure to polychlorinated biphenyls (PCBs)
.
Rev Environ Health
.
2016
;
31
(
1
):
115
119
.
55

Sharma
S
,
Yacavone
M
,
Cao
X
, et al. .

Dietary intake and development of a quantitative FFQ for a nutritional intervention to reduce the risk of chronic disease in the Navajo Nation
.
Public Health Nutr
.
2010
;
13
(
3
):
350
359
.
56

Ballew
C
,
White
LL
,
Strauss
KF
, et al. .

Intake of nutrients and food sources of nutrients among the Navajo: findings from the Navajo Health and Nutrition Survey
.
J Nutr
.
1997
;
127
(
10 suppl
):
2085S
2093S
.
57

Burrows
NR
,
Geiss
LS
,
Engelgau
MM
, et al. .

Prevalence of diabetes among Native Americans and Alaska Natives, 1990–1997: an increasing burden
.
Diabetes Care
.
2000
;
23
(
12
):
1786
1790
.
58

Focazio
MJ
,
Welch
AH
,
Watkins
SA
, et al. . A Retrospective Analysis on the Occurrence of Arsenic in Ground-Water Resources of the United States and Limitations in Drinking-Water-Supply Characterizations. (Water-Resources Investigations Report 99-4279).

Reston, VA
:
US Geological Survey
;
2000
. https://pubs.er.usgs.gov/publication/wri994279. Modified February 2, 2012. Accessed May 2, 2018.
59

Padmaja Divya
S
,
Pratheeshkumar
P
,
Son
YO
, et al. .

Arsenic induces insulin resistance in mouse adipocytes and myotubes via oxidative stress-regulated mitochondrial Sirt3-FOXO3a signaling pathway
.
Toxicol Sci
.
2015
;
146
(
2
):
290
300
.
60

Beck
R
,
Styblo
M
,
Sethupathy
P
.

Arsenic exposure and type 2 diabetes: microRNAs as mechanistic links?
Curr Diab Rep
.
2017
;
17
(
3
):
18
.
61

Tseng
CH
,
Tai
TY
,
Chong
CK
, et al. .

Long-term arsenic exposure and incidence of non-insulin-dependent diabetes mellitus: a cohort study in arseniasis-hyperendemic villages in Taiwan
.
Environ Health Perspect
.
2000
;
108
(
9
):
847
851
.
62

Rahman
M
,
Tondel
M
,
Ahmad
SA
, et al. .

Diabetes mellitus associated with arsenic exposure in Bangladesh
.
Am J Epidemiol
.
1998
;
148
(
2
):
198
203
.
63

Lai
MS
,
Hsueh
YM
,
Chen
CJ
, et al. .

Ingested inorganic arsenic and prevalence of diabetes mellitus
.
Am J Epidemiol
.
1994
;
139
(
5
):
484
492
.
64

Wang
SL
,
Chiou
JM
,
Chen
CJ
, et al. .

Prevalence of non-insulin-dependent diabetes mellitus and related vascular diseases in southwestern arseniasis-endemic and nonendemic areas in Taiwan
.
Environ Health Perspect
.
2003
;
111
(
2
):
155
159
.
65

Gribble
MO
,
Howard
BV
,
Umans
JG
, et al. .

Arsenic exposure, diabetes prevalence, and diabetes control in the Strong Heart Study
.
Am J Epidemiol
.
2012
;
176
(
10
):
865
874
.
66

Feseke
SK
,
St.-Laurent
J
,
Anassour-Sidi
E
, et al. .

Arsenic exposure and type 2 diabetes: results from the 2007–2009 Canadian Health Measures Survey
.
Health Promot Chronic Dis Prev Can
.
2015
;
35
(
4
):
63
72
.
67

Kim
Y
,
Lee
BK
.

Association between urinary arsenic and diabetes mellitus in the Korean general population according to KNHANES 2008
.
Sci Total Environ
.
2011
;
409
(
19
):
4054
4062
.
68

Del Razo
LM
,
García-Vargas
GG
,
Valenzuela
OL
, et al. .

Exposure to arsenic in drinking water is associated with increased prevalence of diabetes: a cross-sectional study in the Zimapán and Lagunera regions in Mexico
.
Environ Health
.
2011
;
10
:
73
.
69

Steinmaus
C
,
Yuan
Y
,
Liaw
J
, et al. .

Low-level population exposure to inorganic arsenic in the United States and diabetes mellitus: a reanalysis
.
Epidemiology
.
2009
;
20
(
6
):
807
815
.
70

Makris
KC
,
Christophi
CA
,
Paisi
M
, et al. .

A preliminary assessment of low level arsenic exposure and diabetes mellitus in Cyprus
.
BMC Public Health
.
2012
;
12
:
334
.
71

Chen
Y
,
Ahsan
H
,
Slavkovich
V
, et al. .

No association between arsenic exposure from drinking water and diabetes mellitus: a cross-sectional study in Bangladesh
.
Environ Health Perspect
.
2010
;
118
(
9
):
1299
1305
.
72

Li
X
,
Li
B
,
Xi
S
, et al. .

Association of urinary monomethylated arsenic concentration and risk of hypertension: a cross-sectional study from arsenic contaminated areas in northwestern China
.
Environ Health
.
2013
;
12
:
37
.
73

Abhyankar
LN
,
Jones
MR
,
Guallar
E
, et al. .

Arsenic exposure and hypertension: a systematic review
.
Environ Health Perspect
.
2012
;
120
(
4
):
494
500
.
74

Jones
MR
,
Tellez-Plaza
M
,
Sharrett
AR
, et al. .

Urine arsenic and hypertension in US adults: the 2003–2008 National Health and Nutrition Examination Survey
.
Epidemiology
.
2011
;
22
(
2
):
153
161
.
75

Waghe
P
,
Sarkar
SN
,
Sarath
TS
, et al. .

Subchronic arsenic exposure through drinking water alters lipid profile and electrolyte status in rats
.
Biol Trace Elem Res
.
2017
;
176
(
2
):
350
354
.
76

Karim
MR
,
Rahman
M
,
Islam
K
, et al. .

Increases in oxidized low-density lipoprotein and other inflammatory and adhesion molecules with a concomitant decrease in high-density lipoprotein in the individuals exposed to arsenic in Bangladesh
.
Toxicol Sci
.
2013
;
135
(
1
):
17
25
.
77

Ettinger
AS
,
Bovet
P
,
Plange-Rhule
J
, et al. .

Distribution of metals exposure and associations with cardiometabolic risk factors in the “Modeling the Epidemiologic Transition Study”
.
Environ Health
.
2014
;
13
:
90
.
78

Su
CT
,
Lin
HC
,
Choy
CS
, et al. .

The relationship between obesity, insulin and arsenic methylation capability in Taiwan adolescents
.
Sci Total Environ
.
2012
;
414
:
152
158
.
79

Gomez-Rubio
P
,
Roberge
J
,
Arendell
L
, et al. .

Association between body mass index and arsenic methylation efficiency in adult women from southwest US and northwest Mexico
.
Toxicol Appl Pharmacol
.
2011
;
252
(
2
):
176
182
.
80

Vahter
M
,
Concha
G
.

Role of metabolism in arsenic toxicity
.
Pharmacol Toxicol
.
2001
;
89
(
1
):
1
5
.
81

Hopenhayn
C
,
Huang
B
,
Christian
J
, et al. .

Profile of urinary arsenic metabolites during pregnancy
.
Environ Health Perspect
.
2003
;
111
(
16
):
1888
1891
.
82

Gardner
RM
,
Nermell
B
,
Kippler
M
, et al. .

Arsenic methylation efficiency increases during the first trimester of pregnancy independent of folate status
.
Reprod Toxicol
.
2011
;
31
(
2
):
210
218
.
This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/about_us/legal/notices)

Supplementary data