Assessment of Grouped Weighted Quantile Sum Regression for Modeling Chemical Mixtures and Cancer Risk

The Association Between a Neighborhood Adverse Childhood Experiences Index and Body Mass Index Among New York City Youth

Background: The role of neighborhood factors in the association between adverse childhood experiences (ACEs) and body mass index (BMI) has not been widely studied. A neighborhood ACEs index (NAI) captures neighborhood environment factors associated with ACE exposure. This study examined associations between BMI and an NAI among New York City (NYC) youth. An exploratory objective examined the NAI geographic distribution across NYC neighborhoods. Methods: Data for students attending NYC public general education schools in kindergarten-12th grade from 2006-2017 (n = 1,753,867) were linked to 25 geospatial datasets capturing neighborhood characteristics for every census tract in NYC. Multivariable hierarchical linear regression tested associations between BMI and the NAI; analyses also were conducted by young (<8 years), school age (8-12 years), and adolescent (>12 years) subgroups. In addition, NAI was mapped by census tract, and local Moran's I identified clusters of high and low NAI neighborhoods. Results: Higher BMI was associated with higher NAI across all sex and age groups, with largest magnitude of associations for girls (medium NAI vs. low NAI: unstandardized β = 0.112 (SE 0.008), standardized β [effect size]=0.097, p < 0.001; high NAI vs. low NAI: unstandardized β = 0.195 (SE 0.008), standardized β = 0.178, p < 0.001) and adolescents (medium NAI vs. low NAI: unstandardized β = 0.189 (SE 0.014), standardized β = 0.161, p < 0.001, high NAI vs. low NAI: unstandardized β = 0.364 (SE 0.015), standardized β = 0.334, p < 0.001 for adolescent girls; medium NAI vs. low NAI: unstandardized β = 0.122 (SE 0.014), standardized β = 0.095, p < 0.001, high NAI vs. low NAI: unstandardized β = 0.217 (SE 0.015), standardized β = 0.187, p < 0.001 for adolescent boys). Each borough of NYC included clusters of neighborhoods with higher and lower NAI exposure, although clusters varied in size and patterns of geographic dispersion across boroughs. Conclusions: A spatial index capturing neighborhood environment factors associated with ACE exposure is associated with higher BMI among NYC youth. Findings complement prior literature about relationships between neighborhood environment and obesity risk, existing research documenting ACE-obesity associations, and the potential for neighborhood factors to be a source of adversity. Collectively, evidence suggests that trauma-informed place-based obesity reduction efforts merit further exploration as potential means to interrupt ACE-obesity associations.

Prenatal per- and polyfluoroalkyl substances (PFAS) and maternal oxidative stress: Evidence from the LIFECODES study

BACKGROUND

Per- and polyfluoroalkyl substances (PFAS) are synthetic chemicals linked to adverse pregnancy outcomes. Although their underlying biological mechanisms are not fully understood, evidence suggests PFAS may disrupt endocrine functions and contribute to oxidative stress (OS) and inflammation.

OBJECTIVE

We examined associations between early pregnancy PFAS exposure and OS biomarkers, exploring potential effect modifications by fetal sex and maternal race.

METHODS

We used data from 469 LIFECODES participants with measured plasma PFAS (median 10 weeks gestation) and repeated measures (median 10, 18, 26, and 35 weeks gestation) of urinary OS biomarkers [8-iso-prostaglandin-F2α (8-isoprostane) and 8-hydroxydeoxyguanosine (8-OHdG)]. Protein damage biomarkers (chlorotyrosine, dityrosine, and nitrotyrosine) were additionally measured in plasma from a subset (N = 167) during the third visit. Associations between each PFAS and OS biomarkers were examined using linear mixed-effects models and multivariable linear regressions, adjusting for potential confounders, including maternal age, race, education level, pre-pregnancy BMI, insurance status, and parity. Effect modifications were evaluated by including an interaction term between each PFAS and fetal sex or maternal race in the models.

RESULTS

We observed significant positive associations between PFOS and 8-isoprostane, with a 9.68% increase in 8-isoprostane levels (95% CI: 0.10%, 20.18%) per interquartile range increase in PFOS. In contrast, PFUA was negatively associated [9.32% (95% CI: -17.68%, -0.11%)], while there were suggestive positive associations for MPAH and PFOA with 8-isoprostane. The associations of several PFAS with 8-OHdG varied by fetal sex, showing generally positive trends in women who delivered females, but negative or null in those who delivered males. No significant effect modification by maternal race was observed.

CONCLUSIONS

This study provides evidence linking PFAS exposure to OS during pregnancy, with potential sex-specific effects of certain PFAS on 8-OHdG. Further research should explore additional OS/inflammatory biomarkers and assess the modifying effects of dietary and behavioral patterns across diverse populations.

A review of common statistical methods for dealing with multiple pollutant mixtures and multiple exposures

Traditional environmental epidemiology has consistently focused on studying the impact of single exposures on specific health outcomes, considering concurrent exposures as variables to be controlled. However, with the continuous changes in environment, humans are increasingly facing more complex exposures to multi-pollutant mixtures. In this context, accurately assessing the impact of multi-pollutant mixtures on health has become a central concern in current environmental research. Simultaneously, the continuous development and optimization of statistical methods offer robust support for handling large datasets, strengthening the capability to conduct in-depth research on the effects of multiple exposures on health. In order to examine complicated exposure mixtures, we introduce commonly used statistical methods and their developments, such as weighted quantile sum, bayesian kernel machine regression, toxic equivalency analysis, and others. Delineating their applications, advantages, weaknesses, and interpretability of results. It also provides guidance for researchers involved in studying multi-pollutant mixtures, aiding them in selecting appropriate statistical methods and utilizing R software for more accurate and comprehensive assessments of the impact of multi-pollutant mixtures on human health.

Association between plasma metal exposure and health span in very elderly adults: a prospective cohort study with mixture statistical approach

BACKGROUND

Metals have been linked to a diverse spectrum of age-related diseases; however, the effects of metal exposure on health span remains largely unknown. This cohort study aims to determine the association between plasma metal and health span in elder adults aged ≥ 90 years.

METHODS

The plasma concentrations of seven metals were measured at baseline in 300 elder adults. The end of the health span (EHS) was identified as the occurrence of one of eight major morbidities or mortality events. We used Cox regression to assess hazard ratios (HR). The combined effects of multiple metal mixtures were estimated using grouped-weighted quantile sum (GWQS), quantile g-computation (Q-gcomp), and Bayesian kernel machine regression (BKMR) methods.

RESULTS

The estimated HR for EHS with an inter-quartile range (IQR) increment for selenium (Se) was 0.826 (95% confidence interval [CI]: 0.737-0.926); magnesium (Mg), 0.806 (95% CI: 0.691-0.941); iron (Fe), 0.756 (95% CI: 0.623-0.917), and copper (Cu), 0.856 (95% CI: 0.750-0.976). The P for trend of Se, Mg, and Fe were all < 0.05. In the mixture analyses, Q-gcomp showed a negative correlation with EHS (P = 0.904), with the sum of the negative coefficients being -0.211.

CONCLUSION

Higher plasma Se, Mg, and Fe reduced the risk of premature end of health span, suggesting that essential metal elements played a role in health maintenance in elder adults.

Prenatal and early life exposure to air pollution and the risk of severe lower respiratory tract infections during early childhood: the Espoo Cohort Study

BACKGROUND

There is inconsistent evidence of the effects of exposure to ambient air pollution on the occurrence of lower respiratory tract infections (LRTIs) in early childhood. We assessed the effects of individual-level prenatal and early life exposure to air pollutants on the risk of LRTIs in early life.

METHODS

We studied 2568 members of the population-based Espoo Cohort Study born between 1984 and 1990 and living in 1991 in the City of Espoo, Finland. Exposure assessment was based on dispersion modelling and land-use regression for lifetime residential addresses. The outcome was a LRTI based on data from hospital registers. We applied Poisson regression to estimate the incidence rate ratio (IRR) of LTRIs, contrasting incidence rates in the exposure quartiles to the incidence rates in the first quartile. We used weighted quantile sum (WQS) regression to estimate the joint effect of the studied air pollutants.

RESULTS

The risk of LRTIs during the first 2 years of life was significantly related to exposure to individual and multiple air pollutants, measured with the Multipollutant Index (MPI), including primarily sulphur dioxide (SO2), particulate matter with a dry diameter of up to 2.5 µm (PM2.5) and nitrogen dioxide (NO2) exposures in the first year of life, with an adjusted IRR of 1.72 per unit increase in MPI (95% CI 1.20 to 2.47). LRTIs were not related to prenatal exposure.

CONCLUSIONS

We provide evidence that ambient air pollution exposure during the first year of life increases the risk of LRTIs during the first 2 years of life. SO2, PM2.5 and NO2 were found to contribute the highest weights on health effects.

Urinary neonicotinoid concentrations and obesity: A cross-sectional study among Chinese adolescents

Epidemiological and toxicological studies on neonicotinoids and obesity have been relevant to adults and young children, but data are limited in adolescents. This study aimed to examine the association between urinary neonicotinoid concentrations and obesity measures among Chinese adolescent. A total of 524 urine samples from 300 boys (11.3-16.1 years) and 224 girls (12.1-15.8 years) were collected to detect the concentrations of eleven neonicotinoids. Generalized linear regression, weighted quantile sum regression (WQS) and Bayesian kernel machine regression (BKMR) were used to estimate covariate-adjusted associations between detectable neonicotinoids and ten indicators of obesity. Nitenpyram concentration was associated with increased body mass index z-score (β = 0.170, 95% CI: 0.041, 0.299) and greater odds of being general obesity (OR = 2.46, 95% CI: 1.11, 5.46). N-desmethyl- acetamiprid concentration was associated with an increase in waist-to-height ratio (β = 0.102, 95% CI: 0.029, 0.176) and waist-to-hip ratio (β = 0.083, 95% CI: 0.011, 0.155). The concentrations of clothianidin (OR = 2.06, 95% CI: 1.10, 3.88) and flonicamid (OR = 2.39, 95% CI: 1.07, 5.32) were associated with greater odds of being abdominal obesity. In contrast, the concentrations of imidacloprid (OR = 0.35, 95% CI: 0.14, 0.88) and thiacloprid (OR = 0.28, 95% CI: 0.08, 0.99) were associated with lower odds of being general obesity. The estimates of general obesity and abdominal obesity increased significantly when concentrations of neonicotinoids mixture were at or above the 55th and 65th percentiles, respectively, compared to the 50th percentile concentration. Sex modified the association between nitenpyram and clothianidin and the risk of obesity with a positive association among boys, and a nonsignificant inverse association among girls. The findings suggest that these associations may be mixed and sex-specific.

Ambient air pollutants and breast cancer stage in Tehran, Iran

This study aimed to examine the impacts of single and multiple air pollutants (AP) on the severity of breast cancer (BC). Data of 1148 diagnosed BC cases (2008-2016) were obtained from the Cancer Research Center and private oncologist offices in Tehran, Iran. Ambient PM10, SO2, NO, NO2, NOX, benzene, toluene, ethylbenzene, m-xylene, p-xylene, o-xylene, and BTEX data were obtained from previously developed land use regression models. Associations between pollutants and stage of BC were assessed by multinomial logistic regression models. An increase of 10 μg/m3 in ethylbenzene, o-xylene, m-xylene, and 10 ppb of NO corresponded to 10.41 (95% CI 1.32-82.41), 4.07 (1.46-11.33), 2.89 (1.08-7.73) and 1.08 (1.00-1.15) increase in the odds of stage I versus non-invasive BC, respectively. Benzene (OR, odds ratio = 1.16, 95% CI 1.01-1.33) and o-xylene (OR = 1.18, 1.02-1.38) were associated with increased odds of incidence of BC stages III & IV versus non-invasive stages. BC stage I and stage III&IV in women living in low SES areas was associated with significantly higher levels of benzene, ethylbenzene, o-xylene, and m-xylene. The highest multiple-air-pollutants quartile was associated with a higher odds of stage I BC (OR = 3.16) in patients under 50 years old. This study provides evidence that exposure to AP is associated with increased BC stage at diagnosis, especially under premenopause age.

Association between plasma and dietary trace elements and obesity in a rural Chinese population

Trace elements may play an important role in obesity. This study aimed to assess the plasma and dietary intake levels of four trace elements, Mn, Cu, Zn and Se in a rural Chinese population, and analyse the relationship between trace elements and obesity. A cross-sectional study involving 2587 participants was conducted. Logistic regression models were used to analyse the association between trace elements and obesity; restricted cubic spline (RCS) models were used to assess the dose-response relationship between trace elements and obesity; the weighted quantile sum (WQS) model was used to examine the potential interaction of four plasma trace elements on obesity. Logistic regression analysis showed that plasma Se concentrations in the fourth quartile (Q4) exhibited a lower risk of developing obesity than the first quartile (Q1) (central obesity: OR = 0·634, P = 0·002; general obesity: OR = 0·525, P = 0·005). Plasma Zn concentration in the third quartile (Q3) showed a lower risk of developing obesity in general obesity compared with the first quartile (Q1) (OR = 0·625, P = 0·036). In general obesity, the risk of morbidity was 1·727 and 1·923 times higher for the second and third (Q2, Q3) quartiles of dietary Mn intake than for Q1, respectively. RCS indicated an inverse U-shaped correlation between plasma Se and obesity. WQS revealed the combined effects of four trace elements were negatively associated with central obesity. Plasma Zn and Se were negatively associated with obesity, and dietary Mn was positively associated with obesity. The combined action of the four plasma trace elements had a negative effect on obesity.

Identification of dietary components in association with abdominal aortic calcification

The precise impact of dietary components on vascular health remains incompletely understood. To identify the dietary components and their associations with abdominal aortic calcification (AAC), the data from NHANES was employed in this cross-sectional study. The LASSO method and logistic regression were utilized to identify dietary components that exhibited the strongest association with AAC. Grouped WQS regression analysis was employed to evaluate the combined effects of dietary components on AAC. Furthermore, principal component analysis was employed to identify the primary dietary patterns in the study population. The present analysis included 1862 participants, from whom information on 35 dietary macro- and micronutrient components was obtained through 24-hour dietary recall interviews. The assessment of AAC was performed utilizing dual-energy X-ray absorptiometry. The LASSO method identified 10 dietary components that were associated with AAC. Total protein, total fiber, vitamin A, and β-cryptoxanthin exhibited a negative association with AAC. Compared to the first quartile, the adjusted odds ratios (95% CIs) for the highest quartile were 0.59 (0.38, 0.93), 0.63 (0.42, 0.93), 0.59 (0.41, 0.84), and 0.68 (0.48, 0.94), respectively. Grouped WQS regression demonstrated a positive association between the lipid group and AAC (aOR: 1.29; 95% CI: 1.12, 1.50), while the proteins and phytochemical group exhibited a negative association with AAC (aOR: 0.69; 95% CI: 0.58, 0.82). For the dietary pattern analysis, high adherence to the plant-based pattern (aOR: 0.62; 95% CI: 0.44, 0.88) was associated with a lower risk of AAC, whereas the caffeine and theobromine pattern (aOR: 1.73; 95% CI: 1.25, 2.41) was associated with a higher risk of AAC. The findings of this study indicate that adopting a dietary pattern characterized by high levels of protein and plant-based foods, as well as reduced levels of fat, may offers potential advantages for the prevention of AAC.

Association between Vitamin A and E Forms and Prostate Cancer Risk in the Singapore Prostate Cancer Study

PURPOSE

This study aimed to assess associations between forms of vitamin A and E (both individually and collectively) and the risk of prostate cancer, as well as identify potential effect modifiers.

METHODS

Utilizing data from the Singapore Prostate Cancer Study, a hospital-based case-control study, we measured the serum concentrations of 15 different forms of vitamins A and E in 156 prostate cancer patients and 118 control subjects, using a high-performance liquid chromatography technique. These forms included retinol, lutein, zeaxanthin, α-cryptoxanthin, β-cryptoxanthin, α-carotene, β-carotene, lycopene, ubiquinone, δ-tocopherol, γ-tocopherol, α-tocopherol, δ-tocotrienol, γ-tocotrienol, and α-tocotrienol. The odds ratio and 95% confidence interval for associations between vitamin A and E and prostate cancer risk were estimated using logistic regression models after adjustment for potential confounders. The analyses were further stratified by smoking and alcohol consumption status. The mixture effect of micronutrient groups was evaluated using weighted quantile sum regression.

RESULTS

Higher concentrations of retinol, lutein, α-carotene, β-carotene, ubiquinone, α-tocopherol, δ-tocotrienol, γ-tocotrienol, and α-tocotrienol were significantly and positively associated with overall prostate cancer risk. Among ever-smokers, associations were stronger for lutein, β-cryptoxanthin and β-carotene compared with never-smokers. Among regular alcohol drinkers, associations were stronger for lutein, β-cryptoxanthin, ubiquinone, γ-tocotrienol and α-tocotrienol compared with non-regular alcohol drinkers. Retinol and α-tocotrienol contributed most to the group indices 'vitamin A and provitamin A carotenoids' and 'vitamin E', respectively.

CONCLUSIONS

Several serum vitamin A and E forms were associated with prostate cancer risk, with significant effect modification by smoking and alcohol consumption status. Our findings shed light on prostate cancer etiology.

Modeling historic environmental pollutant exposures and non-Hodgkin lymphoma risk

Environmental exposures to chemicals are suspected risk factors for non-Hodgkin lymphoma (NHL), but few studies have assessed historic environmental risk factors. In this study, we estimated the associations between NHL and 1) historic environmental pollutant emissions from the Risk Screening Environmental Indicators (RSEI) model, which uses a database from the Environmental Protection Agency of toxic release emissions to air, water, and land, and 2) chemical mixtures measured in house dust (groups of PCBs, PAHs, and two mixtures of pesticides) for study participants enrolled in the NCI-SEER population-based case-control study (1998-2000) at four SEER centers - Detroit, Iowa, Los Angeles County, and Seattle. We assigned 11 years of annual temporally-varying historic environmental exposure scores by intersecting residential locations from participants' residential histories with a fine grid from the RSEI model and by performing inverse distance weighting between facilities releasing specific carcinogenic chemicals and residential locations for spatially-precise exposure assignments. We used Bayesian index low-rank kriging multiple membership models to identify important lag times for RSEI scores, cumulative effects of RSEI scores, and specific carcinogenic chemical releases into the environment. We found a significant positive association between RSEI scores and NHL at the maximum time lag of 11 years (OR = 1.17, 95% CI (1.06, 1.32)) and a significant cumulative RSEI score effect (OR = 1.30, 95% CI (1.02, 1.84)) for long-term residents in Detroit, where benzene and trichloroethylene were the most important chemicals driving this association. Additionally, we identified significant inverse associations for two study centers and time lags that did not persist in cumulative exposure models. Large weights for dichloromethane and pentachlorophenol in models of cumulative exposure also support evidence for their association with NHL risk. These results underscore the importance of considering historic and cumulative environmental exposures and using residential histories for diseases with long latency periods such as NHL.

Prenatal Exposure to Legacy and Alternative Per- and Polyfluoroalkyl Substances and Neuropsychological Development Trajectories over the First 3 Years of Life

The neurotoxic effects of prenatal exposure to per- and polyfluoroalkyl substances (PFAS) on offspring animals are well-documented. However, epidemiological evidence for legacy PFAS is inconclusive, and for alternative PFAS, it is little known. In this investigation, we selected 718 mother-child pairs from the Chinese Maoming Birth Cohort Study and measured 17 legacy and alternative PFAS in the third-trimester serum. Neuropsychological developments (communication, gross motor function, fine motor function, problem solving ability, and personal-social skills) were assessed at 3, 6, 12, 18, 24, and 36 months using the Ages and Stages Questionnaires 3rd edition. Trajectories of each subscale were classified into persistently low and persistently high groups via group-based trajectory modeling. Logistic regression and grouped weighted quantile sum were fitted to assess the potential effects of individual PFAS and their mixtures, respectively. Higher linear PFHxS levels were associated with elevated odds for the persistently low trajectories of communication (OR = 1.73; 95% CI: 1.12, 2.66) and problem solving ability (OR = 2.11; 95% CI: 1.14, 3.90). Similar findings were observed for linear PFOS, 1m-PFOS, PFDA, PFDoDA, PFUnDA, and legacy PFAS mixture. However, no association was observed for alternative PFAS and their mixture. We provided insights into the longitudinal links between prenatal legacy/alternative PFAS exposure and neuropsychological development trajectories over the first 3 years of life.

Concurrent dietary intake to nitrate, thiocyanate, and perchlorate is negatively associated with hypertension in adults in the USA

We aimed to comprehensively evaluate the association of urinary nitrate, thiocyanate, and perchlorate metabolites with hypertension among a nationally representative sample of the US adult population. This cross-sectional study investigated data from 15,717 adults aged more than 20 years obtained from the National Health and Nutritional Examination Survey (NHANES) for the years 2005-2016. In the survey, urinary levels of nitrate, thiocyanate, and perchlorate were measured using ion chromatography combined with electrospray tandem mass spectrometry. Blood pressure was calculated as the mean of three measurements. Hypertension was defined as (a) systolic BP ≥130 and/or diastolic BP ≥80 mmHg and/or (b) self-report. Multivariate logistic regression and weighted quantile sum (WQS) regression models were applied to estimate the association between exposure to multiple inorganic anions and hypertension. Restricted cubic spline (RCS) regressions were fitted to discern the potential relationship between the anion exposure and hypertension. These innovation methods used to support our results. Overall, 7533 (49.95%) people with and 7638 (50.35%) without hypertension were included in this study. In the multivariable-adjusted logistic regression models, urinary nitrate (P < 0.001) and perchlorate (P < 0.001) were independently negatively associated with increased occurrence of hypertension, while urinary thiocyanate was insignificantly associated with hypertension (P = 0.664). The WQS regression index showed that, in combination, the three inorganic anions mixture were negatively correlated with hypertension (adjusted OR 0.89; 95% CI 0.83-0.95, P < 0.001). Urinary nitrate was the most heavily weighted component in the hypertension model (weight = 0.784). RCS regression demonstrated that nitrate (nonlinearity P = 0.205) and perchlorate (nonlinearity P = 0.701) were linearly associated with decreased occurrence of hypertension. Concurrent exposure to nitrate, thiocyanate, and perchlorate is associated with a decreased risk of hypertension, with the greatest influence coming from nitrate probably; urinary specific thiocyanate alone had an insignificant association with hypertension.

Preconception exposure to perfluoroalkyl and polyfluoroalkyl substances and couple fecundity: A couple-based exploration

Numerous studies have examined the adverse health effects of perfluoroalkyl and polyfluoroalkyl substances (PFAS) but it remains unclear whether preconception exposure to PFAS affects couple fecundity. This prospective preconception cohort study with 936 Chinese couples aimed to comprehensively assess the effects of PFAS on couple fecundity [measured by the time to pregnancy (TTP)] and infertility (i.e., TTP > 12 menstrual cycles) with a focus on the effects of partner-specific exposure and joint-effects of couple-based exposure. Twenty-five PFAS were quantified in plasma from each partner, including seven branched isomers, two chlorinated polyfluoroalkyl ether sulfonic acids, four emerging PFAS replacements [i.e., 6:2 fluorotelomer phosphate diester (6:2 diPAP) and three short-chain alternatives: perfluoro-n-butanoicacid, perfluorobutane sulfonate and perfluoroheptanoic acid (PFHpA)]. Using a two-phase regression approach composed of elastic net regression and principal component analysis, we found that exposure to 6:2 diPAP and PFHpA rather than legacy PFAS in women and the couple-based exposure patterns characterized by high level of female 6:2 diPAP were significantly associated with reduced couple fecundity, which was independent of the adjustment of co-exposed PFAS homogenous from both partners. For example, a ln unit increase in female 6:2 diPAP was associated with 15 % [fecundity odds ratio (FOR) = 0.85, 95 %CI: 0.76, 0.96)] lower odds of couple fecundability (i.e., longer TTP) and 45 % increased risk of infertility [OR = 1.45 (95 %CI: 1.16, 1.81)], respectively. While most PFAS in men were not associated with couple fecundity, certain PFAS (e.g., perfluorohexane sulfonic acid) in men were negatively associated with infertility risk. However, the combined effects of PFAS mixture in couples were nonsignificant. Our findings suggest that PFAS in men and women may exert different impacts on couple fecundity. Preconception exposure to 6:2 diPAP and PFHpA in women may have the potential to impair couple fecundity.

Public healthcare costs associated with long-term exposure to mixtures of persistent organic pollutants in two areas of Southern Spain: A longitudinal analysis

BACKGROUND

Polychlorinated biphenyls and organochlorine pesticides are persistent organic pollutants (POPs) that had been banned or restricted in many countries, including Spain. However, their ubiquity still poses environmental and human health threats.

OBJECTIVE

To longitudinally explore public healthcare costs associated with long-term exposure to a mixture of 8 POPs in a cohort of residents of two areas of Granada Province, Southern Spain.

METHODS

Longitudinal study in a subsample (n = 385) of GraMo adult cohort. Exposure assessment was performed by analyzing adipose tissue POP concentrations at recruitment. Average primary care (APC) and average hospital care (AHC) expenditures of each participant over 14 years were estimated using the data from their medical records. Data analyses were performed by robust MM regression, weighted quantile sum regression (WQS) and G-computation analysis.

RESULTS

In the adjusted robust MM models for APC, most POPs showed positive beta coefficients, being Hexachlorobenzene (HCB) significantly associated (β: 1.87; 95% Confidence interval (95%CI): 0.17, 3.57). The magnitude of this association increased (β: 3.72; 95%CI: 0.80, 6.64) when the analyses were restricted to semi-rural residents, where β-HCH was also marginally-significantly associated to APC (β: 3.40; 95%CI: -0.10, 6.90). WQS revealed a positive but non-significant mixture association with APC (β: 0.14; 95%CI: -0.06, 0.34), mainly accounted for by β-HCH (54%) and HCB (43%), that was borderline-significant in the semi-rural residents (β: 0.23; 95%CI: -0.01, 0.48). No significant results were observed in G-Computation analyses.

CONCLUSION

Long-term exposure to POP mixtures might represent a modifiable factor increasing healthcare costs, thus affecting the efficiency of the healthcare systems. However, and owing the complexity of the potential causal pathways and the limitations of the present study, further research is warranted to fully elucidate ascertain whether interventions to reduce human exposure should be considered in healthcare policies.

Association of urinary or blood heavy metals and mortality from all causes, cardiovascular disease, and cancer in the general population: a systematic review and meta-analysis of cohort studies

Amounting epidemiological evidence has shown detrimental effects of heavy metals on a wide range of diseases. However, the effect of heavy metal exposure on mortality in the general population remains unclear. The primary objective of this study was to clarify the associations between heavy metals and mortality from all causes, cardiovascular disease (CVD), and cancer based on prospective studies. We comprehensively searched Pubmed, Embase, and Web of Science electronic databases to identify studies published from their inception until 1 March 2022. Investigators identified inclusion criteria, extracted study characteristics, and assessed the methodological quality of included studies according to standardized guidelines. Meta-analysis was conducted if the effect estimates of the same outcome were reported in at least three studies. Finally, 42 original studies were identified. The results of meta-analysis showed that cadmium and lead exposure was significantly associated with mortality from all causes, CVD, and cancer in the general population. Moderate evidence suggested there was a link between arsenic exposure and mortality. The adverse effects of mercury and other heavy metals on mortality were inconclusive. Epidemiological evidence for the joint effect of heavy metal exposure on mortality was still indeterminate. In summary, our study provided compelling evidence that exposure to cadmium, lead, and arsenic were associated with mortality from all causes, CVD, and cancer, while the evidence on other heavy metals, for example mercury, was insignificant or indeterminate. Nevertheless, further prospective studies are warranted to explore the joint effects of multiple metal exposure on mortality.

The Intersection of Neighborhood Environment and Adverse Childhood Experiences: Methods for Creation of a Neighborhood ACEs Index

This study evaluated methods for creating a neighborhood adverse childhood experiences (ACEs) index, a composite measure that captures the association between neighborhood environment characteristics (e.g., crime, healthcare access) and individual-level ACEs exposure, for a particular population. A neighborhood ACEs index can help understand and address neighborhood-level influences on health among individuals affected by ACEs. Methods entailed cross-sectional secondary analysis connecting individual-level ACEs data from the Philadelphia ACE Survey (n = 1677) with 25 spatial datasets capturing neighborhood characteristics. Four methods were tested for index creation (three methods of principal components analysis, Bayesian index regression). Resulting indexes were compared using Akaike Information Criteria for accuracy in explaining ACEs exposure. Exploratory linear regression analyses were conducted to examine associations between ACEs, the neighborhood ACEs index, and a health outcome—in this case body mass index (BMI). Results demonstrated that Bayesian index regression was the best method for index creation. The neighborhood ACEs index was associated with higher BMI, both independently and after controlling for ACEs exposure. The neighborhood ACEs index attenuated the association between BMI and ACEs. Future research can employ a neighborhood ACEs index to inform upstream, place-based interventions and policies to promote health among individuals affected by ACEs.

Association and mediation analyses among multiple metals exposure, plasma folate, and community-based impaired estimated glomerular filtration rate in central Taiwan

BACKGROUND

Chronic kidney disease (CKD) is increasing, with heavy metal exposure an important risk factor. Additionally, the antioxidant folic acid has been studied for reducing blood arsenic levels and related tissue damage. Therefore, we explored the association and mediation effects among various heavy metal levels in blood, plasma folate, other CKD risk factors, and impaired estimated glomerular filtration rate (eGFR).

METHODS

We constructed a community-based cross-sectional study from the Human Biomonitoring and Environmental Health Program in central Taiwan. A total of 1643 participants had lived locally for > 5 years, > 40 years old, and completely received health examinations and biospecimen collections. Impaired eGFR was defined as one single eGFR < 60 mL/min/1.73 m². Plasma folate and metal levels in blood were determined, as well as urinary 8-hydroxy-2'-deoxyguanosine as an oxidative stress marker. Generalized weighted quantile sum (WQS) regression analysis was used to calculate a WQS score, reflecting overall body-burden of multiple metals (arsenic, cadmium, chromium, nickel, and lead) in blood.

RESULTS

Impaired eGFR was identified in 225 participants. Participants with high WQS scores had increased risk of impaired eGFR (odds ratio = 1.67; 95% confidence interval [CI]: 1.34, 2.07). Of five metals, arsenic, lead, and cadmium were weighted highly in impaired eGFR. Participants with high WQS and folate insufficiency (< 6 ng/mL) had 2.38-fold risk of impaired eGFR compared to those with low WQS and high folate (≥6 ng/mL) (95% CI: 1.55, 5.17). Similar increased 4.16-fold risk of impaired eGFR was shown in participants with high WQS and uric acid levels (95% CI: 2.63, 6.58). However, there were no significant WQS-folate (p = 0.87) or WQS-uric acid (p = 0.38) interactions on impaired eGFR risk. As a mediator, uric acid contributed 24% of the association between WQS score and impaired eGFR risk (p < 0.0001). However, no mediation effect of plasma folate was observed.

CONCLUSION

WQS analysis could be applied to evaluate the joint effects of multiple metals exposure. High WQS scores may influence impaired eGFR risk through increased uric acid levels. A large-scale and prospective cohort study is necessary to validate these results and demonstrate any causal relationship.

Imputation of Below Detection Limit Missing Data in Chemical Mixture Analysis with Bayesian Group Index Regression

There is growing scientific interest in identifying the multitude of chemical exposures related to human diseases through mixture analysis. In this paper, we address the issue of below detection limit (BDL) missing data in mixture analysis using Bayesian group index regression by treating both regression effects and missing BDL observations as parameters in a model estimated through a Markov chain Monte Carlo algorithm that we refer to as pseudo-Gibbs imputation. We compare this with other Bayesian imputation methods found in the literature (Multiple Imputation by Chained Equations and Sequential Full Bayes imputation) as well as with a non-Bayesian single-imputation method. To evaluate our proposed method, we conduct simulation studies with varying percentages of BDL missingness and strengths of association. We apply our method to the California Childhood Leukemia Study (CCLS) to estimate concentrations of chemicals in house dust in a mixture analysis of potential environmental risk factors for childhood leukemia. Our results indicate that pseudo-Gibbs imputation has superior power for exposure effects and sensitivity for identifying individual chemicals at high percentages of BDL missing data. In the CCLS, we found a significant positive association between concentrations of polycyclic aromatic hydrocarbons (PAHs) in homes and childhood leukemia as well as significant positive associations for polychlorinated biphenyls (PCBs) and herbicides among children from the highest quartile of household income. In conclusion, pseudo-Gibbs imputation addresses a commonly encountered problem in environmental epidemiology, providing practitioners the ability to jointly estimate the effects of multiple chemical exposures with high levels of BDL missingness.

Residential exposure to carbamate, organophosphate, and pyrethroid insecticides in house dust and risk of childhood acute lymphoblastic leukemia

BACKGROUND

Self-reported residential use of pesticides has consistently been associated with increased risk of childhood leukemia. However, these studies were limited in their ability to identify specific insecticide active ingredients that were associated with risk.

OBJECTIVE

We used household carpet dust measurements of 20 insecticides (two carbamate, 10 organophosphate, two organochlorine, and six pyrethroid) as indicators of exposure and evaluated associations with the risk of childhood acute lymphoblastic leukemia (ALL).

METHODS

We conducted a population-based case-control study of 252 ALL cases diagnosed from 1999 to 2007 and 306 birth certificate controls from 35 counties in Central and Northern California. Carpet dust was collected at a second interview (2001-2007) for cases who had not moved since diagnosis (comparable reference date for controls) using a specialized vacuum cleaner in the room where the child spent most of their time or from the household vacuum. Insecticides were categorized as detected (yes/no), or as tertiles or quartiles of their distributions among controls. We calculated odds ratios (OR) and 95% confidence intervals (CI) using unconditional logistic regression adjusting for demographic characteristics, interview year, and season of dust collection.

RESULTS

Permethrin, chlorpyrifos, diazinon, and carbaryl were the most frequently detected insecticide active ingredients. When we compared the highest quartile to the lowest or to non-detections, there was no association with ALL for permethrin (OR Q4 vs. Q1 = 0.81; 95% CI 0.50-1.31), carbaryl (OR Q4 vs. non-detects = 0.61, 95% CI 0.34-1.08) or chlorpyrifos (OR Q4 vs. Q1 = 0.60; 95% CI 0.36-1.00). The highest quartile of diazinon concentration was inversely associated with risk in the single pesticide model but without a monotonic exposure-response (p-trend = 0.14). After adjusting for other common insecticides, the OR was not significant (OR Q4 vs. Q1 = 0.58; 95% CI 0.33-1.05). None of the other insecticides were associated with risk.

CONCLUSION

Our results should be interpreted within the limitations of the case-control study design including the use of a single post-diagnosis dust sample and restriction to residentially stable participants, which may have resulted in selection bias. Although difficult to implement, additional studies with assessment of exposure to insecticide active and non-active ingredients are necessary to elucidate the role of these common exposures in childhood leukemia risk.

Bayesian Group Index Regression for Modeling Chemical Mixtures and Cancer Risk

There has been a growing interest in the literature on multiple environmental risk factors for diseases and an increasing emphasis on assessing multiple environmental exposures simultaneously in epidemiologic studies of cancer. One method used to analyze exposure to multiple chemical exposures is weighted quantile sum (WQS) regression. While WQS regression has been demonstrated to have good sensitivity and specificity when identifying important exposures, it has limitations including a two-step model fitting process that decreases power and model stability and a requirement that all exposures in the weighted index have associations in the same direction with the outcome, which is not realistic when chemicals in different classes have different directions and magnitude of association with a health outcome. Grouped WQS (GWQS) was proposed to allow for multiple groups of chemicals in the model where different magnitude and direction of associations are possible for each group. However, GWQS shares the limitation of WQS of a two-step estimation process and splitting of data into training and validation sets. In this paper, we propose a Bayesian group index model to avoid the estimation limitation of GWQS while having multiple exposure indices in the model. To evaluate the performance of the Bayesian group index model, we conducted a simulation study with several different exposure scenarios. We also applied the Bayesian group index method to analyze childhood leukemia risk in the California Childhood Leukemia Study (CCLS). The results showed that the Bayesian group index model had slightly better power for exposure effects and specificity and sensitivity in identifying important chemical exposure components compared with the existing frequentist method, particularly for small sample sizes. In the application to the CCLS, we found a significant negative association for insecticides, with the most important chemical being carbaryl. In addition, for children who were born and raised in the home where dust samples were taken, there was a significant positive association for herbicides with dacthal being the most important exposure. In conclusion, our approach of the Bayesian group index model appears able to make a substantial contribution to the field of environmental epidemiology.