ABSTRACT
Humans are extensively exposed to Quaternary Ammonium Compounds (QACs). QACs are ubiquitously used in medical settings, restaurants, and homes as cleaners and disinfectants. They are also used on food and in personal care products as preservatives. Despite their prevalence, nothing is known about the health effects associated with chronic low-level exposure. Chronic QAC toxicity was recently identified in mice and resulted in developmental and reproductive deficits and altered immune function. Cell based studies show that QACs increase inflammation, disrupt cholesterol synthesis, and decrease mitochondria function. If these studies translate to human toxicity, multiple physiological functions could be affected. QAC concentrations in humans have not been monitored previously. This study tested whether QAC concentrations could be detected in the blood of 43 random volunteers, and whether QAC concentrations were associated with markers of inflammation, mitochondrial function, and cholesterol synthesis in a dose dependent manner. QAC concentrations were detected in 80% of study participants, and were associated with decreased mitochondrial function and an increase in inflammatory cytokines in a dose dependent manner. Cholesterol synthesis pathway intermediaries were generally increased, indicating disruption in cholesterol homeostasis. This is the first study to demonstrate that chronic exposure to QACs results in measurable concentrations in human blood, and to also demonstrate significant correlations between QAC level and meaningful biomarkers related to health.
INTRODUCTION
Humans are extensively exposed to quaternary ammonium compounds (QACs). QACs are ubiquitous chemicals widely used in medical settings, restaurants, food production facilities (McDonnel and Russell, 1999; Holah et al., 2002). They are also used in the home as common cleaners and disinfectants. QACs are frequently added to personal care products and on fresh foods as preservatives. Additionally, QACs, have wide-ranging industrial applications. Despite their prevalence, little is known about the extent of human exposure and the consequences of chronic low-level contact.
Two common QACs are alkyl dimethyl benzyl ammonium chloride (ADBAC, also termed benzalkonium chloride or BAC) and didecyl dimethyl ammonium chloride (DDAC). Over 750,000 tons of both BAC and DDAC are manufactured each year in the US and in Europe (Steichen, 2001). This level of production identifies BAC and DDAC as High Production Volume Chemicals by the US EPA. During production of BAC, a mixture of products is formed. These products contain a benzyl group bound to central nitrogen with alkyl chains composed of differing numbers of carbons units. The average chain length and the distribution of chain lengths vary according to the starting material and the synthetic process used (Danish EPA, 2000). Most processes favor production of chain lengths containing 10 to 18 carbons. This means that BAC solutions contain a mixture of distinct chemicals, each having their own toxicity. Commercial cleaning products often contain mixtures that are predominantly C14 BAC or C12 BAC, followed by C16 BAC. Additionally, many commercial products also contain 2 or more QAC species in combination, such BAC and DDAC. This is important as chemical mixtures can act synergistically or antagonistically to produce an effect that is different from the sum of the individual components. Evaluation of common mixtures is essential in determining chemical exposure risk (Hermens et al., 1985; Faust et al., 2003; Narotsky et al., 2011).
QAC disinfectants are historically viewed as having low toxicity. Accidental overdoses from ingestion of the concentrated mix have resulted in human deaths, but these events are rare (Chataigner et al., 1991; Ellenhorn et al., 1997). Chronic exposure is known to cause asthma and contact dermatitis (Bernstein et al., 1994; Suneja and Belsito, 2008), as well as ocular inflammation and hypersensitivity (Hong and Bielory, 2009). Risk assessments, conducted by the US EPA in 2006, determined that both BAC and DDAC pose little risk to humans (EPA 2006a, 2006b). When the risk assessment is examined in detail, it is clear that the determination was based on studies with decreased body weight as the toxicological endpoint. These studies, using rats and rabbits, found decreased weight gain in adults, neonates, and fetuses at the higher doses (Tyl, 1989; Neeper-Bradley, 1990, 1993). No other specific and subtle toxicological endpoints have been tested. More importantly, no studies have evaluated the extent and magnitude of human exposure, or the possible systemic response to this exposure.
Chronic QAC toxicity has only recently been identified. We have shown using in-vitro studies on both human and mouse cell lines that QACs increase inflammation, disrupt mitochondrial function, alter estrogen signaling, and inhibit cholesterol synthesis (McDonald, 2017; Datta et al., 2017; Herron et al., 2016). In mice, we have shown that a commercial cleaner containing a mixture of BAC+DDAC used to disinfect the vivarium, affected reproduction, embryo development, and immune function (Melin et al., 2014, 2015; Hrubec et al., 2017; McDonald, 2017). We have also shown that BACs can cross the placenta and alter cholesterol and lipid homeostasis in mouse neonatal brains after gestational exposure through the maternal diet (Herron et al., 2019). If these studies translate to human toxicity, many basic physiological functions and disease processes could be altered. The body of evidence from our animal and cell based studies suggests toxicity on both a cellular and organismal level. Extrapolation from animal and in-vitro studies is the foundation for assessing the safety of drugs and chemicals. Thus, there is a high potential for QAC exposure to cause adverse effects in humans. There is, however, a complete data gap on the effects of chronic low-level systemic exposure to QACs. The objective of this study was to determine the extent of human exposure and determine if there is a systemic response associated with exposure.
METHODS
Our cell and mouse-based studies suggested that QACs promote increased proinflammatory cytokines, decreased mitochondrial function, and inhibition of cholesterol synthesis. Therefore, we focused on these potential consequences in this current assessment of human QAC exposure.
Sample Population Characteristics
The blood samples were collected from participants recruited from Blacksburg, VA in the USA. Blacksburg is a small, college town, largely comprised of individuals associated with Virginia Tech, a large university. Blood samples were collected anonymously, and no personal identifying information was collected. Based on visual categorization of age, approximately 2/3rds of the study participants were likely college students, and 1/3rd older adults likely affiliated with the university in some capacity. Inclusion criteria included: non pregnant individuals who were greater than 18 years of age with no history of chronic illness, no acute illness, and previous blood draw in the last two weeks. All research was approved by the IRB (project number 1302366-8) for the protection of human subjects.
Sample Size Calculation
We hypothesize that there is an association between measured toxicological endpoints and QAC concentration in blood. By estimating the correlation coefficient to be 0.5 and setting α=0.05 and power=0.9, at least 31 blood samples were needed to determine associations between blood analytes and QAC concentrations.
Sample Collection
Blood samples were collected from 43 healthy volunteers by a certified phlebotomist following IRB approved protocols to protect human subjects. Blood was collected from each participant into 3 separate heparinized blood tubes. One tube was used for mitochondrial analysis, a second tube was used for sterol and QAC analysis, while the third tube was used for cytokine analysis. Determination of each blood analyte was made blind to the participant’s QAC concentration.
Mitochondrial Function
Mitochondrial physiology was assessed with an XF24 Extracellular Flux Analyzer (Seahorse System, Agilent, Santa Clara, CA) which determines O2 consumption under different test conditions. Briefly, 6M white blood cells were aliquoted onto Seahorse XF24 plates at ∼20,000 cells per well. Oxygen consumption rate (OCR) was determined at basal rate (ATP synthesis inhibited by oligomycin), maximal mitochondrial stimulation from addition of FCCP, and at full mitochondrial inhibition from addition of Antimycin A/rotenone (ORC from Proton Leak).
Sterol Analysis
Analysis of Cholesterol and the cholesterol synthesis precursors: 7-Dehydrocholesterol (7-DHC), 8-Dehydrocholesterol (8-DHC), 7-Dehydrodesmosterol (7-DHD), Lanosterol, Zymosterol, Lathosterol, and Desmosterol, were determined by mass spectrometry. Lipids were extracted from the whole blood and the sterols analyzed by reverse phase LC-MS/MS assays as described previously (Herron et al., 2018).
Cytokine Production
Cytokines were determined by ELISA. C-reactive protein (CRP) was determined in plasma (R&D Systems, Minneapolis, MN) and NF-κB determined in whole blood (ABNOVA, Taipei, Taiwan). The remaining cytokines, IL-6, IL-10, IL-12, and TNFα, were determined in plasma (R&D Systems) under three conditions. The first determination was an unstimulated baseline prior to culture and stimulation. The remaining whole blood was divided into two aliquots. One received no additives while the other was stimulated with addition of 1µg/mL LPS. Both were cultured simultaneously as a whole blood culture for 12 hours at 37°C with 5% CO2. After culture, plasma was collected and assessed for cytokine concentration. There was no difference between the unstimulated baseline and the unstimulated cultured control, thus the unstimulated cultured control was not included in further analysis. IL-10 is generally considered an anti-inflammatory cytokine, but stimulation with LPS alone converts macrophages to the M1 pro-inflammatory phenotype and results in pro-inflammatory IL-10 production.
QAC Determination
QAC levels in the blood were measured to accurately gauge exposure. Blood samples were spiked with known amounts of deuterium-labeled (d7-benzyl) BACs as internal standards (Herron et al., 2016). All samples were extracted by Folch solution (chloroform/methanol = 2/1). After extraction, the samples were re-constituted in the LC solvent and then analyzed by UPLC-MS/MS (Seguin et al., 2019).
QAC Nomenclature
The nomenclature used to describe QACs varies across publications due to the large number of chemicals in the QAC class, as well as the fact that many common QACs have multiple synonyms. For the purposes of this paper, the following nomenclature will be used. QAC is used to refer to both BAC and DDAC together. BAC, with no modifier, is used when speaking of the various BAC alkyl chain lengths together, while the modifiers C10, C12, C14 and C16 BAC are used to designate the specific BAC chemical with the indicated alkyl chain length.
Statistical Analysis
Univariate descriptive statistics were conducted. QAC concentrations were then compared with absolute values of mitochondrial function, sterol synthesis intermediaries, and cytokine production to determine associations between the toxicological endpoints and the blood QAC concentration. As both the toxicological endpoints and QAC concentration measurements are continuous, the nonparametric Spearman’s Rank Correlation (rs) was used to measure associations between all endpoints. Each QAC species was ranked high to low and correlated with the associated analyte value for that sample. Associations were calculated across the ranked samples, increasing by 1 sample (N) in each calculation. For example, 10 individuals contained C12 BAC in their blood. The 10 samples were ranked from high to low. Correlations were determined between C12 BAC and each analyte for the five individuals with the highest C12 BAC. A second correlation was then calculated for individuals with the six highest C12 BAC concentration. The process was repeated until correlations were determined for all 10 individuals containing C12 BAC. Correlation coefficients of 0.5 with a p < 0.05 were used to identify significance. This is an accepted level of association in studies relating concentrations of an environmental contaminant with a pathophysiological outcome (Li et al., 2011; Jenkins et al., 2002; Bárány et al., 2002). While this correlation coefficient may seem low compared to a tightly controlled laboratory study, the normal variability in blood constituents, and variations in exposure level from randomly sampled individuals, necessitates a more relaxed measure of significant association. For means comparison tests, each blood analyte was stratified into two groups for each QAC; samples with no QAC and samples with QAC. Difference in the means of the two groups was determined by a t-test for parametric data, or Wilcoxin Rank Sum Test for non-parametric data. Each blood analyte was also stratified into low, medium and high evenly spaced QAC terciles. Differences between the means of each tercile were compared by ANOVA for parametric data, or Kruskal Wallis test for non-parametric data. Nineteen of the twenty-seven analytes measured (including QACs) exhibited outliers (determined by Statistix Software, Talahasee, FL). Typically only a single high outlier was seen per analyte; four analytes had two outliers. All outliers were greater than two standard deviations from the mean; 73% of identified outliers were greater than three standard deviations from the mean. The identified outliers were not included in the analysis.
RESULTS
Approximately 80% of individuals contained measurable concentrations of QACs in their blood. Blood concentrations of each QAC species are given in supplementary data (Table S1). One-third of individuals had Total QAC concentrations between 10 – 150 nM, which has been shown to have physiological effects in cell culture models (Datta et el., 2017; Herron et al., 2016). Descriptive statistics for the complete dataset are given in Table S1. Data were analyzed in two ways; by means comparison, and by association with each QAC. The means comparison test evaluates differences between concentration groups rather than at the individual level, and includes individuals with no QAC residues. The association test measures the direct linear relationship between two variables from each individual across the sample population. Thus, the association test evaluates data at the individual level in relation to the whole dataset, and does not include individuals without QACs in their blood.
Descriptive statistics for all measured analytes. The varying sample number (N) reflects loss of some samples during the shipping and handling process.
Means Comparison
For analysis, the measured blood concentrations were stratified into the following exposure groups: no QAC, and then low, mid, and high evenly spaced QAC tercile groups. Exceptions were made for C10 and C12 BAC. Only a single individual contained detectable C10 BAC and no further analysis was conducted for this BAC. Furthermore, as only a limited number of individuals contained C12 BAC residues, this BAC was divided into no C12, low C12, and high C12 groups with no mid-range group. Means were then compared between the groups to identify if a dose response was present.
Eleven individuals had C12 BAC residues in their blood; one individual’s value was identified as an outlier (6.5 standard deviations from the mean) and was not included in the analysis. For C12 BAC, significant differences in means were observed for Baseline IL-12 and Stimulated IL-10 (Fig. 1). Both analytes were significantly elevated in individuals with the higher concentrations of C12 BAC in their blood.
Differences in selected analytes by C12 and C14 BAC concentration found in the blood. Individual bar graphs depict the mean ± SE for different analytes. Bars identified by different letters indicate a significant difference (p ≤ 0.05) between groups. Analytes stratified into two groups were compared by either t-test or Wilcoxin Rank Sum Test, while those stratified into four groups were compared by ANOVA or Kruskal Wallis test. Scatterplot correlations were determined by Spearman’s Rank Correlation (rs).
Twenty individuals contained C14 BAC residues in their blood; two individual’s values were identified as outliers (3.8 and 5.5 standard deviations from the mean) and were not included in the analysis. C14 BAC affected IL-6 a common marker of inflammation demonstrating a dose response across the tercile groups (Fig. 1).
Twenty-four individuals contained C16 BAC residues in their blood; one individual’s value was identified as an outlier (6.3 standard deviations from the mean) and was not included in the analysis. C16 BAC affected sterol homeostasis, inflammation, and mitochondrial function (Fig. 2). Individuals with C16 BAC residues in their blood demonstrated an inverted-U nonmonotonic dose response for the cholesterol synthesis intermediary 8-DHC with the highest concentrations found at the low and mid C16 terciles (Fig. 2). The inflammatory markers, Stimulated IL-6 and Stimulated IL-10, demonstrated a J-shaped nonmonotonic dose response with significantly lower blood concentrations in the low and mid C16 BAC terciles. The mitochondrial respiration marker, Maximum OCR from Basal, generally decreased across the exposure groups. Mitochondrial function also demonstrated a significant decrease in individuals with C16 BAC exposure compared to those with C16 BAC in their blood (Fig. 2).
Differences in selected analytes by C16 BAC concentration found in the blood. Individual bar graphs depict the mean ± SE for different analytes. Bars identified by different letters indicate a significant difference (p ≤ 0.05) between groups. Analytes stratified into two groups were compared by either t-test or Wilcoxin Rank Sum Test, while those stratified into four groups were compared by ANOVA or Kruskal Wallis test. Scatterplot correlations were determined by Spearman’s Rank Correlation (rs).
DDAC residues were seen in 32 of the 43 participants; one individual’s value was identified as an outlier (5.9 standard deviations from the mean) and was not included in the analysis. DDAC affected cholesterol biosynthetic pathway intermediates, inflammatory markers, and mitochondrial function (Fig. 3). 8-DHC exhibited a difference between individuals with DDAC in their blood compared to those with no DDAC exposure. Baseline TNFα was more than an order of magnitude higher in individuals with DDAC present in their blood than those without DDAC. Baseline IL-10 and Stimulated IL-6 both exhibited a J-shaped nonmonotonic dose response with the lowest IL-10 concentrations found in the low and mid DDAC groups.
Differences in selected analytes by DDAC concentration found in the blood. Individual bar graphs depict the mean ± SE for different analytes. Bars identified by different letters indicate a significant difference (p ≤ 0.05) between groups. Analytes stratified into two groups were compared by either t-test or Wilcoxin Rank Sum Test, while those stratified into four groups were compared by ANOVA or Kruskal Wallis test. Scatterplot correlations were determined by Spearman’s Rank Correlation (rs).
Among our sample population, 35 of the 43 individuals contained QAC in their blood; two values were identified as outliers (2.9 and 5.9 standard deviations from the mean) and were not included in the analysis. When the analytes were evaluated in relation to the Total QAC concentration, the responses were similar to those seen for the individual QAC that had the greatest influence on blood analyte (Fig. 4). This is not surprising as the Total QAC determination is just the sum of all 4 QAC species of meaningful exposure (C12, C14, C16 BACs, and DDAC). Stimulated IL-6 showed a dose response similar to that seen with C14 BAC, C16 BAC, and DDAC. An additive response was seen across the QAC species for the means comparison test of Maximal ORC from Basal. Mitochondrial respiration was decreased at higher concentrations of both C16 BAC and DDAC, however a greater decline in respiration was seen when the two QACs were combined in Total QAC analysis. A cancellation of effect was seen with 8-DHC. 8-DHC concentrations were increased with both C16 BAC and DDAC, but with Total QAC. This cancellation of effect is likely due to the modest declines in 8-DHC concentrations with C-12 and C14 BAC blunting the overall response.
Differences in selected analytes by Total QAC concentration found in the blood. Individual bar graphs depict the mean ± SE for different analytes. Bars identified by different letters indicate a significant difference (p ≤ 0.05) between groups. Analytes stratified into two groups were compared by either t-test or Wilcoxin Rank Sum Test, while those stratified into four groups were compared by ANOVA or Kruskal Wallis test. Scatterplot correlations were determined by Spearman’s Rank Correlation (rs).
Associations
Significant associations (rs) were observed between a number of analytes and specific QAC species. As expected, the associations varied with the number of samples included in the calculation. All calculated associations are given in Supplemental Tables 2 - 6.
Spearman’s Rank Correlations (rs) for C12 BAC. Blue cells designate r > 0.5 and p ≤ 0.05. Grey cells designate r > 0.5 and p ≤ 0.1. No notable correlations were observed above N=9 and are not included in the table.
Spearman’s Rank Correlations (rs) for C14 BAC. Blue cells designate r > 0.5 and p ≤ 0.05. Grey cells designate r > 0.5 and p ≤ 0.1. No notable correlations were observed above N=15 and are not included in the table.
Spearman’s Rank Correlations (rs) for C16 BAC. Blue cells designate r > 0.5 and p ≤ 0.05. Grey cells designate r > 0.5 and p ≤ 0.1 Blank cells represents removed outliers, and samples lost during shipping. No notable correlations were observed above N=19 and are not included in the table.
Spearman’s Rank Correlations (rs) for DDAC. Blue cells designate r > 0.5 and p ≤ 0.05. Grey cells designate r > 0.5 and p ≤ 0.1. No notable correlations were observed above N=24 and are not included in the table.
Spearman’s Rank Correlations (rs) for Total QAC. Blue cells designate r > 0.5 and p ≤ 0.05. Grey cells designate r > 0.5 and p ≤ 0.1. No notable correlations were observed above N=26 and are not included in the table.
For C12 BAC, the cholesterol pathway constituent Lanosterol demonstrated a strong positive correlation (rs=0.833, p=0.008, N=9) with C12 concentrations between 4.468 and 0.828 nM across 9 of the 10 individuals (Figs. 1 & 5). The increase in the precursor Lanosterol indicates a slowdown in a downstream step of the cholesterol synthesis pathway or a general upregulation of cholesterol production.
Significant associations of C12, C14, and C16 BAC with blood analytes. The number of individuals with C12, C14, and C16 BAC detected in their blood was 10, 18, and 22 respectively. Each QAC species was ranked high to low and correlated with the associated analyte value for that sample. The N is the number of ranked samples included in the calculation. Correlation rs values were determined by the nonparametric Spearman’s Rank Correlation test. Correlation coefficients ≥ 0.5 with a p ≤ 0.05 are shaded blue and considered significant. Correlation coefficients ≥ 0.5 with a p ≤ 0.1 are shaded grey to designate associations with less certainty. Significant correlations were not found beyond an N of 19 for C16 BAC and are not shown.
For C14 BAC, both cholesterol synthesis intermediates and inflammatory markers demonstrated associations (Figs. 1 & 5). 7-DHD demonstrated a strong negative association (rs=-0.825, p=0.002, N=12) which remained significant over a range of N=10 to N=14 (from 14.383 to 3.337 nM C14 BAC) and with less certainty for N=9. Baseline TNFα was significantly correlated (rs=-0.700, p=0.028, N=10) with C14 BAC; however, the association was inconsistently present. On the other hand, Stimulated TNFα was strongly and significantly associated with C14 BAC across much of N=8 to N=14 (from 14.383 to 3.337 nM C14 BAC) with the highest association seen at rs=0.817, p=0.010, N=9 (Fig. 1 & 5).
Associations were observed with C16 BAC only for inflammatory intermediaries (Fig. 2 & 5). Stimulated TNFα, and Stimulated IL-6 were weakly associated with C16 BAC (rs=0.640, p=007, N=17; and rs=0.563, p=0.013, N=19, respectively). The association was significant from C16 BAC concentrations of 6.411 to 0.935 nM. On the other hand IL-10, in both Stimulated and Baseline samples, was positively associated with C16 BAC. Stimulated IL-10 showed strong correlations across the sample population from N=6 to N=19. The strongest correlation was seen with N=8 (rs=0.905, p=0.004) but significant correlations continued through N=19 (from 6.411 to 0.935 nM C16 BAC). Baseline IL-10 was more weakly correlated (rs=0.655, p=0.004). While usually considered an anti-inflammatory cytokine, IL-10 has both pro and anti-inflammatory properties. Stimulation with LPS alone as was done in this study, converts macrophages to the M1 pro-inflammatory phenotype and results in pro-inflammatory IL-10 production. It is not clear whether IL-10 production in the Baseline sample represented a pro-inflammatory or anti-inflammatory stimulus.
Strong positive correlations were found with DDAC and a number of cholesterol synthesis pathway intermediaries including: Cholesterol, Zymosterol, Desmosterol, 7-DHD, and Lanosterol (Fig. 3, Fig. 6). For Cholesterol, the strongest correlation (rs=0.893, p=0.012) was seen with N=7, but significant correlations were seen from N=6 to N=9 (from 16.421 to 7.432 nM). For Zymosterol, the strongest correlation (rs=0.929, p=0.006) was seen with N=7, but significant correlations were seen from N=6 to N=13 (from 16.421 to 5.940 nM). For Desmosterol, the strongest correlation (rs=0.929, p=0.006) was seen with N=7, with significant correlations seen from N=6 to N=9. Associations with less certainty were seen to N=13 (from 16.421 to 7.432 nM). For Lanosterol, the strongest correlation (rs=0.893, p=0.012) was seen with N=7, but significant correlations were seen from N=6 to N=11 (from 16.421 to 6.479 nM). In addition to cholesterol synthesis pathway intermediaries, Maximal Mitochondrial OCR from Basal showed strong negative correlations with DDAC (rs=-0.745, p=0.012) to N=11 (from 16.421 to 6.479 nM) (Figs. 3, Fig. 6, Table S5). Additionally, the inflammatory markers CRP and Stimulated IL-10 were positively associated with DDAC concentration. Thus, while DDAC had the greatest influence on the cholesterol synthesis pathway, mitochondrial and inflammatory markers were affected as well.
Significant associations of DDAC and Total QAC with blood analytes. The number of individuals with DDAC, and with the sum of all QAC species (Total QAC) in their blood was 31 and 33 respectively. Each QAC species was ranked high to low and correlated with the associated analyte value for that sample. The N is the number of ranked samples included in the calculation. Correlation values were determined by the nonparametric Spearman’s Rank Correlation test. Correlation coefficients ≥ 0.5 with a p ≤ 0.05 are shaded blue and are considered significant. Correlation coefficients ≥ 0.5 with a p ≤ 0.1 are shaded grey to designate associations with less certainty. Significant correlations were not found beyond an N of 20 for DDAC and 26 for Total QAC, and are not shown.
Total QAC concentration represents the sum of the 4 individual constituents C12, C14, C16 BACs, and DDAC. Associations of analytes in relation to Total QAC concentration can be enhanced from what is seen with each QAC species if the effects are additive, or potentiated. On the other hand, associations can be diminished if there are divergent responses between individual QAC species, or if a response to a single QAC is diluted by a lack of response to the other QAC components. Increased Total QAC concentrations were strongly and significantly associated (Figs. 4 & 6) with decreased maximal mitochondrial function (rs=-0.706, p=0.006) from N=11 to N=16. Less certain associations were seen from N=8, over a total QAC range from 37.4 to 5.7 nM. While mitochondrial function was decreased with individual QAC species, the greatest inhibition was seen with the Total QAC comparison. This indicates that the effect in the QAC mixture is likely additive. Strong significant associations were also observed for both Baseline and Stimulated IL-10 production (Fig. 4 &6). These significant associations were seen over much of the sample population. While much of this effect came from the response to C16 BAC, the other BAC species were also positively associated with IL-10 production (see supplemental data) resulting in an additive response to Total QAC. The strong associations between C16 BAC and both TNFα, and IL-6 were not as prominent in the Total QAC mixture. This is most likely due to an obscuring of these responses when Total QACs were assessed as a whole.
DISCUSSION
While contradictory evidence is mounting, the widely held belief is that although exposure to BAC and DDAC occurs, little if any accumulates systemically as these chemicals were assumed to pass through the intestine largely unabsorbed. This belief is based on unpublished rodent studies conducted for regulatory purposes (Henderson, 1992; US EPA 2006a, 2006b). For BAC, 90-98% of an orally administered dose was eliminated in the feces, and only 0.03 to 0.5% of the dose was retained in body tissues (US EPA 2006a). For DDAC, 89-99% of an oral dose was eliminated, unabsorbed, in the feces and only 0.003-0.675% of the DDAC dose remained in the tissue (Henderson 1992). However, these studies DDAC and BAC rely upon a single oral dose; tissue accumulation following chronic exposure, has not been assessed.. Furthermore, substantial recovery of QAC in feces does not definitively demonstrate low absorption. For example, approximately 1/3rd of excreted BAC and 40% of excreted DDAC are oxidized metabolites (EPA 2006a; EU, Directive 98/8/EC, 2015); and we have recently shown that BAC is converted to oxidized metabolites by hepatic cytochrome P450 enzymes expressed in the liver (Seguin, et al. 2019). The presence of the oxidized metabolites in the feces suggests orally administered BAC and DDAC were absorbed and reached the liver. Others have also determined that the major excretory route of high molecular weight QACs such C12, C14, and C16 BACs, and DDAC is via bile into feces (Hughes, et al. 1973; Smit, et al. 1998; Song, et al. 2013). Thus, the extent of oxidized QAC metabolites in the feces could in fact reflect a greater than appreciated extent of intestinal absorption, followed by hepatic metabolism and biliary excretion back into feces.
QAC exposure can also occur by other routes. Xue et al. (2004) found much higher blood and tissue concentrations when BAC was aspirated into the lungs, or given parenterally than when dosed orally. This is consistent with the efficient metabolism of BACs by liver cytochrome P450 enzymes following an oral dose, as aspiration or parenterally administration would not undergo hepatic first-pass metabolism. We have observed the same pattern of effect in our reproductive and developmental toxicity studies in terms of routes of exposure. Undosed mice that received ambient exposure (most likely inhalant) from use of QAC disinfectants in the vivarium, showed similar degrees of reproductive and developmental impairment as mice receiving both oral and ambient exposure (Melin et al., 2015; Hrubec et al., 2017). This implies that most of the toxicity results from the ambient/inhalant exposure. This is troubling as ambient exposure more closely mimics human exposure from use of cleaning products. There are virtually no studies investigating the systemic toxicity from inhalant exposure of BAC or DDAC. A single study in rats, evaluating inhalant toxicity of cetylpyridinium chloride, a common but structurally unrelated QAC, found mortality after a single exposure in all dose groups including the lowest dose of 50 mg dust/m3 (Lin et al., 1991). Calculations of the LD50 indicated that systemic toxicity was greater with the inhalant route than other routes. Taken together, these studies show that respiratory tract exposure can result in higher blood concentrations and increased toxicity than oral dosing. Clearly, the human population is exposed, and equally clearly, our results demonstrate that QACs can be absorbed systemically and be detected in human blood, but the main route or routes of exposure are not known.
Approximately 80% of participants in our study had detectable concentrations of BAC and DDAC in their blood; and some individuals contained high concentrations between 10 – 150 nM. The primary source of exposure responsible for the blood residues is not known as QACs are used occupationally as well as in home and public settings. Occupational exposure is an obvious candidate for primary exposure, but daily cleaning of the home in poorly ventilated spaces may result in equivalent uptake of QACs. Additionally, toxicokinetic studies have not been conducted in humans; and thus, the distribution of QACs in the body, rate of metabolism, and eventual clearance is not known. The BAC and DDAC detected in our study population could indicate prior exposure that day, or be the result of low level exposure over months. In 2006, the US EPA risk assessment identified dermal exposure to the concentrated disinfectant, and inhalant exposure to the working concentration of disinfectant as posing the greatest risk to people. This was based on reports of dermal burns for the former and reports of asthma, allergy and hypersensitivity for the latter. The risk assessment did not evaluate internal systemic exposure as it was presumed that BAC and DDAC do not accumulate in the body. Identifying the main source and the route of systemic exposure is critical in order to modify or limit one’s internal dose.
The different QAC species were detected in the following abundance: DDAC > C14 BAC > C16 BAC > C12 BAC > C10 BAC. This roughly follows the proportion of QACs found in most commercial disinfection products; DDAC is usually twice the concentration of the total BAC component, and C12 BAC is frequently double the concentration of C16 BAC. In our blood samples, the average DDAC concentration was only slightly greater than C14 BAC, and C12 BAC was roughly the same concentration as C16 BAC. Single individuals did not necessarily contain all QAC species. Only one individual contained detectable C10 BAC and only 11 individuals contained C12 BAC. DDAC was the most commonly found QAC and was present in 32 individuals. This difference in QAC species distribution indicates either differences in absorption or elimination rates, or that the blood concentrations reflect exposure to multiple formulations and different sources of QAC components in the environment.
While most individuals clustered within a certain range of QAC concentrations, some individuals had a value that was up to an order of magnitude higher than the next closest value. These extreme measurements were not attributable to a single individual, but were found in different individuals. Corresponding measured analytes in these extreme samples were frequently divergent as well. In all cases, neither the QAC outlier nor the associated analyte were included in the statistical analysis. The numbers of sample at these extreme concentrations of QACs were too few to assess a valid response.
Increasing QAC concentrations accompanied significant increases in inflammatory mediators, decreases to mitochondrial function, and decreases in cholesterol synthesis. Analyte response varied with the different alkyl chain lengths of BAC. It is generally held that longer chain QACs have lower toxicity than shorter alkyl chain QACs due to decreased water solubility (Danish EPA, 2000). However, this does not match our previous observations nor the findings of this study. We found that C16 BAC affected more constituents than the other BAC species, despite the fact that C16 BAC was found at lower concentrations in the blood. Overall C16 BAC had the greatest effect on markers of inflammation while DDAC predominantly affected cholesterol synthesis. Mitochondrial function was decreased both with C16 BAC and with DDAC, but demonstrated the strongest response with Total QAC indicating the response was likely additive across the QAC species tested.
In some instances, the means comparison and the association tests gave different results. This is not unreasonable as the two tests measure different aspects of effect. Association tests can tell the concentration range over which there is a potential effect. Ideally associations should be determined on sufficient numbers of individuals in the mid to high exposure range. When only a few individuals at the higher dose range are evaluated, there are too few individuals to establish a true relationship. Conversely, when sufficient individuals at the lowest blood concentrations are included, associations disappear as the QAC concentration falls below the no-observed-effect-level (NOEL). The NOEL indicates the concentration below which the measured QAC has no discernable effect on the particular blood analyte. Chemicals with higher NOELs are less toxic then ones with lower NOELs. From our data, it appears that DDAC has a higher NOEL as across the board the correlations dropped off sooner than for the BAC species.
Means comparison tests can show a dose response across dose groups. Nonmonotonic or hormetic responses are by definition not linear, and thus will not give a linear association across the entire concentration range. This could result in statistically significant associations across only a portion of sample population as was seen for a number of the analytes measured in this study. The idea of a toxicant having a nonmonotonic /hormetic dose response is not new, and is reviewed in detail by Calabrese (2005). First coined in the 1940s, hormetic dose responses demonstrate either an inverted U-shaped curve with decreased responses at low and high doses, or a J-shaped curve with enhanced responses at low and high doses. In an analysis of the published toxicological response literature, Calabrese and Baldwin (2001) estimated that the overall frequency of hormetic dose response across the spectrum of known toxicants is approximately 40%. This includes responses to many well recognized toxicants such as DDT, mercuric chloride, lead, cadmium, chromium, toluene, and others (Calabrese 2005). With regards to QACs, an early study evaluating immunologic effects of BAC, found that histamine release from mast cells followed a nonlinear dose response (Danish EPA, 2000). Additionally, BAC binds to the acetycholine muscarinic receptor, and also the estrogen receptor (Takeuchi et al., 1993; Datta et al., 2017). Many receptor mediated signaling systems, including those involving the muscarinic, and estrogen receptors, have a biphasic nonmomotonic dose response (Calabrese 2005).
Regardless of whether the response is linear or nonmonotonic, our data clearly demonstrated a dose dependent association between blood analytes and the amount of QAC present in the blood. Cytokine production appeared heavily impacted with increases in both baseline and stimulated responses at higher QAC concentrations. Baseline measurements represent the resting/current concentration, while stimulated samples represent the capacity to respond after appropriate signals. The inflammatory effects were particularly striking for TNFα, IL-6 and IL-10. Mitochondrial function was also strongly and significantly inhibited in a dose dependent manner. These changes are similar to those found in our cell and mouse studies. Also corroborating our previous studies, we found that cholesterol synthesis intermediaries were profoundly affected. Increasing concentrations of DDAC were significantly associated with increasing Zymosterol, 7-Dehydrocholesterol, Desmosterol, Lanosterol, and cholesterol while C12 BAC was positively associated with Lanosterol. Previously we demonstrated in-vitro that BAC inhibition of the final step of cholesterol synthesis led to upstream accumulation of the synthesis precursors, particularly 7-DHC and 8-DHC (Herron et al., 2016). However, results in the present study suggest that there is an overall upregulation of cholesterol biosynthesis by DDAC. Interestingly, as C14 BAC increased, 7-DHD decreased indicating that C14 BAC may act directly to inhibit conversion of Zymosterol to 7-DHD. These results suggest that different types of QACs may have different effects on cholesterol biosynthesis.
The correlations found in the present study are strong. The accepted significant correlation coefficient for evaluating environmental effects in human tissues is 0.5 (Li et al., 2011; Jenkins et al., 2002; Bárány et al., 2002). We observed significant correlations from 0.5 to 0.9 for a number of the blood analytes measured. It must be kept in mind however that this is a small pilot study with a limited sample size. No personal information about the study population was collected. As such, confounding influence of external variables such as diet, weight, exercise status, or pharmaceutical use of cholesterol modulating drugs is not known. Even with this caveat, these data are provocative and indicate possible ramifications for health and disease. If the changes observed in this study translate to the population at large, and are associated with pathophysiology, then QAC exposure could be involved in many and varied disease processes. Potentially impacted may be inflammatory diseases such as diabetes, obesity, osteoarthritis, autoimmunity, and cancer; mitochondrial associated diseases such as cardiovascular disease, obesity, and neurodegenerative disease. Disruptions to lipid homeostasis could potentially influence many essential biochemical pathways as cholesterol is required for synthesis of all steroid hormones and is critically involved in cell signaling, cell proliferation, phagocytosis, and cell motility.
In animal studies, we have demonstrated declines in male and female reproduction (Melin et al., 2014, 2015), and severe defects in neural development (Hrubec et al., 2018). We also have found that chronic oral dosing during breeding and gestation resulted in maternal absorption and transfer of C12 BAC and C16 BAC across the placenta to the fetus (Herron et al., 2019). BAC and DDAC also cross the blood-testis barrier to alter spermatogenesis (Melin, 2015). QACs disrupt cell signaling through alterations in sterol and oxysterol formation (Herron et al., 2019) and through receptor mediated G protein Gi/Go signaling (Higashijima et al., 1990). Taken together there is a wealth of evidence suggesting that BAC and DDAC have profound effects on basic biochemical and physiological processes. The data provided by this pilot study indicates that QACs may affect human responses as well.
CONCLUSIONS
This study clearly demonstrated that QAC residues are present in human blood, and that alterations in inflammatory markers, mitochondrial function, and cholesterol synthesis were observed in a dose dependent manner with QAC concentration. Therefore, further studies are critically needed. It is not known if the exposure identified here correlates with disease in the human population. The primary sources of exposure and toxicokinetic properties such as absorption, distribution, metabolism, and elimination in humans is also not known. The limitations of our pilot study, including the small sample size and lack of information on other influential variables prevent us from answering these important questions. As the first study to identify possible systemic effects from chronic exposure in humans, however, the data clearly indicate the potential for health ramifications from QAC exposure.
Data Availability
Data in the manuscript are available
Footnotes
Communicating Author: Terry C. Hrubec, Department of Biomedical Science, E. Via College of Osteopathic Medicine – Virginia, Blacksburg, VA thrubec{at}vt.vcom.edu, 540-231-1702
CITATIONS
