Heart rate variation and human body burdens of environmental mixtures in the Cree First Nation communities of Eeyou Istchee, Canada

•This article examines contaminant mixtures and HRV in an Indigenous population.•Findings advance the understanding of mixture effects associations with HRV.•The findings helped identify novel components associated with HRV. Heart rate variability (HRV) is a measure of cardiac autonomic regulation that examines the variation in beat-to-beat fluctuations in heart rate. While many exposure-based studies have examined the effects of single or similar groups of contaminants on HRV parameters, none have examined the association between complex environmental mixtures, including organic and elemental contaminants, and HRV. Using data collected from the Multi-Community Environment-and-Health Study in Eeyou Istchee (Quebec, Canada), we assessed HRV in two time domain measures: root mean square of successive differences (RMSSD) and standard deviation of the N-N (RR) intervals (SDNN); and in three frequency domains: high frequency (HF), low-frequency (LF), and very-low frequency (VLF) in 443 participants. We first examined mixture effects of nineteen organic and metal contaminants in blood using principal component analysis (PCA) and a multivariable general linear regression on HRV responses, adjusting for age, sex, body mass index, smoking status, and kidney disease covariates. We subsequently assessed HRV outcome response variables using Bayesian kernel machine regression (BKMR) to further examine individual contaminant contribution and overall mixture effects. In the PCA, a significant positive association was observed between RMSSD and principal component (PC) axis 2, which was highly positively-loaded for nickel and moderately negatively-loaded for mercury. A negative association between SDNN and PC-1, which was highly positively-loaded for all PCBs (polychlorinated biphenyls) and organochlorines and moderately positively-loaded for mercury, was observed. Additionally, a significant and positive association was observed between PC-2 and SDNN and a significant and negative association between PC-3 (negatively loaded for cadmium) and LF. Associations with contaminants were not observed for HF or VLF. BKMR results suggest that trans-nonachlor and cis-nonachlor are primarily responsible for reductions in HRV; however, their contributions to HRV deficits was non-significant when examined as a complete mixture. While PCA results suggested that organochlorines, mercury, nickel, and cadmium may all play a role in altering various HRV, subsequent BKMR analysis illustrate