EUR
en
Electronic nicotine delivery systems (ENDS) are devices that use an electric current instead of tobacco combustion to vaporize nicotine-containing e-liquids for a user to inhale. Since their commercial release in 2003, ENDS have been marketed as smoking cessation tools and safer alternatives to conventional cigarettes because they emit lower, yet detectible levels of carcinogens, ultrafine particles, and other lung irritants. With the advent of concealable USB-like designs and flavors such as fruit, menthol, and candy, ENDS rapidly became the most popular nicotine delivery product used by adolescents in the U.S. From 2011 to 2015, ENDS use increased by 900% in U.S. middle schoolers and high schoolers, prompting the U.S. Food and Drug Administration (FDA) and the U.S. Surgeon general to declare a youth vaping epidemic.
In response to the vaping epidemic, the National Academies of Sciences, Engineering, and Medicine (NASEM) and the U.S. Surgeon General published reports calling for research to characterize ENDS aerosol constituents, which are largely unknown due to the lack of U.S. FDA regulation before 2016 and no current testing standards for acceptable emissions. Since then, volatile organic compounds (VOCs), heavy metals, and ultrafine particles have been detected in ENDS emissions across multiple devices and e-liquids. Importantly, ENDS have been shown to release higher concentrations of ultrafine particles relative to conventional cigarettes, and at levels exceeding health-based limits. This is concerning given ultrafine particles can penetrate deep into the respiratory tract and act as carriers for toxic chemicals to enter the blood stream. Additionally, the Flavor and Extract Manufacturers Association (FEMA) has identified over 1000 flavorings commonly used in e-liquids that may pose a respiratory hazard due to possible volatility and irritant properties. Not surprisingly, International Agency for Research on Cancer (IARC) group 1 and 2 carcinogens, including acrolein, acetaldehyde, formaldehyde, and benzene, have been detected in ENDS emissions. Several of these VOCs are also associated with respiratory infections and lung diseases such as chronic obstructive pulmonary disease (COPD) and asthma. Therefore, regardless of how ENDS compare to conventional cigarettes, ENDS emit potentially hazardous toxicants, making their impact on respiratory health a critical question of utmost importance to public safety.
In contrast to primary exposure, little is known about the health effects of inhaling secondhand ENDS aerosols, which are created when an ENDS user exhales near bystanders or non-users. This is particularly concerning because bystanders are exposed to many of the same potentially toxic ultrafine particles, VOCs, and metals that are present in firsthand emissions according to the World Health Organization (WHO). Importantly, ENDS usage in indoor workplaces and public enclosed places is unregulated in several countries, as well as the majority of U.S. states, enhancing the risk of unintentional exposure. Using the U.S. EPA unacceptable risk measurements for long-term exposure, the WHO recently concluded formaldehyde and acetaldehyde concentrations presented excess lifetime cancer risk, and concentrations of particulate matter (PM) increased mortality risk in exposed populations. Moreover, indoor particulate matter emission rates from secondhand ENDS emissions exceeded those from secondhand cigarette smoke in some studies, and secondhand particulate matter at vaping conventions reached a level of 12 times higher than the U.S. EPA 24-h standard. Still, several gaps remain in the literature on ENDS secondhand emissions. For example, most studies to-date use smoking machines and fail to distinguish the impact of different devices and vaping topographies on indoor air quality. Because the physicochemical profile of firsthand aerosols varies across ENDS devices, e-liquids, and vaping topographies, future studies are needed that make these same comparisons for secondhand emissions in diverse populations, especially as ENDS devices continue to evolve. These gaps are a prevailing barrier to our understanding of the short- and long-term human health impact of ENDS on bystanders.
To address these knowledge gaps, the present study characterizes the levels of particles and VOCs present in firsthand and secondhand emissions from ten established ENDS devices and users, respectively. We characterized the impact of puffing topography on emissions and found the ENDS users puff volume and flow rate are powerful mediators of firsthand and secondhand particle and chemical exposures. Using computational modeling, we showed substantial particle deposition and VOC inhalation risks for ENDS users, youth, and adult bystanders. Furthermore, untargeted high-resolution metabolomic analysis of saliva from ENDS users revealed firsthand emission exposures elicited unique profiles of lung injury biomarkers associated with asthma and lung cancer pathogenesis. Taken together, these results demonstrate the importance of vaping behavior such as puff volume in mediating particle and chemical exposures, as well as metabolic and respiratory health outcomes in individuals exposed to firsthand and secondhand ENDS emissions.
This human subject study was approved by Advarra Institutional Review Board (IRB Protocol No. H22081). All methods performed in this study, including recruitment, informed consent, and human data collection, was carried out in accordance with the approved IRB protocol (No. H22081). We utilized targeted social media strategies in conjunction with a recruitment agency to recruit adult participants aged 18–35 who could speak, read, and write in English who have the presence of at least 16 natural teeth with at least 2 abutting teeth in each sextant. This age group was selected as individuals older than age 35 may be more likely to have used other tobacco products throughout the course of their lifetime, increasing the difficulty of elucidating the specific contribution of ENDS usage to disease endpoints. In addition, participants were excluded if were currently using a nicotine patch, had active substance abuse, are currently using tobacco products, were pregnant or lactating, had a history of diabetes, HIV infection, or were currently prescribed immunosuppressant medications or antibiotics or if another household member was enrolled in the study. Eligible individuals provided informed consent and were scheduled for in-person visits with the study team.
A total of 19 participants (nine never established ENDS users and ten established ENDS users) were consented and enrolled in the study. Never established ENDS users were defined as individuals who have used less than 100 tobacco products in their lifetime. Survey results collected in our study, show the most recent use of ENDS products was one year prior to the study for two participants with the remaining never established ENDS users never having used ENDS products in their lifetime. Established ENDS users were defined as individuals who have not used combustible or any other form of tobacco in the past 90 days and exclusively used ENDS daily. These groups are defined based on current CDC definitions and were used in previous investigations. Recent evidence indicates cessation of conventional tobacco use, results in reduced inflammation within 30 days. Groups were composed of equal numbers of men and women and participants were not included/excluded based on socioeconomic status, race, or poverty status. Upon their visit, saliva samples were collected, puffing topography assessments were completed, and exhaled emission characterization were assessed as described further below.
The Clinical Research Support System (CReSS; Pocket; Borgwaldt-KC, Richmond, VA) device was used to collect real-time puffing topography metrics, including puff volume (mL/puff), flow rate (mL/s), peak flow (mL/puff), inter-puff interval (s), puff duration (s), and time of peak (s) for each ENDS user. For completion of the puffing topography assessment, each participant’s own ENDS device was attached to the CReSS Pocket via a universal device adaptor and participants were instructed to take 10 puffs from their own device as they normally would. Collected real-time puffing topography information was utilized to establish puffing topography profiles for each participant and these parameters were further applied to produce firsthand exposure scenarios from their purchased device via the EAGS platform as described below (IES Technologies, Morgantown, WV).
To investigate differential levels of ENDS aerosol exposure in real-life scenarios, three distinct exposure scenarios (baseline, firsthand, and secondhand) were employed. Baseline, firsthand, and secondhand exposure scenarios were evaluated by characterizing the emissions from never established ENDS users (n = 9), device emissions (n = 10), and established ENDS users (n = 10), respectively. In the baseline scenario, never established ENDS users generated 20 breaths into a 1 m 3 stainless steel environmental exposure chamber. In the secondhand exposure scenario, ENDS users inhaled the aerosol generated by their personal ENDS devices attached to a CReSS Pocket device, as described above, and exhaled emissions via a sterile mouthpiece into a 6m 3 chamber, providing a total of 10 puffs for analysis. To determine device emissions, personal ENDS devices from each user were purchased and utilized to generate firsthand emission. ENDS user puffing topography profiles, were applied to the EAGS platform, controlled by the IES software (IES Technologies, Morgantown, WV) allowing for the production of firsthand emissions from the users device which match their individual distinct puffing profile. The emissions from each scenario were characterized for particulates and VOC, with normalization based on the number of puffs for further analysis. All participant experiments were conducted with the environmental exposure chamber at static state with a minimal air exchange rate and supply air flow was applied only to compensate sampling air flow.
Particle sizes and concentrations present in generated aerosols were monitored for each scenario (e.g., baseline, firsthand, secondhand) by using the Scanning Mobility Particle Sizer (SMPS; models TSI 3080 and 3786, TSI Inc., Shoreview, MN) for 7–300 nm sized particles and Optical Particle Sizer (OPS; TSI 3330, TSI Inc., Shoreview, MN) for 0.3–10 µm sized particles. Considering the short time duration of the aerosol generation process, time resolution of the measurements for both pieces of equipment was set at 1 min. Total particle count, total mass, and geometric mean diameter (GMD) of the emitted aerosol were monitored and calculated for each scenario. For the particle size distribution analysis, data from SMPS and OPS has been merged to obtain entire particle count distribution (dN/d log Dp).
The multiple-path particle dosimetry (MPPD; version 3.04, Applied Research Associates Inc., Albuquerque, NM) was used to determine potential particle inhalation risks and particle deposition estimates. For the analysis, the Stochastic model was selected to mimic real-life ENDS usage and exposure with input parameters as follows: constant exposure with oral breathing scenario; breathing frequency 12 per minute; tidal volume 625 mL; exposure time duration set at 6 h per day, 7 days per week for 52 weeks. Particle emission data collected from each scenario and participant, including count mean diameter, geometric standard deviation, and mass concentration (mg/m 3) were also utilized as inputs to determine the average deposited dose within the airways. Based on particle emission information inputted in the system, the average deposited dose was calculated and compared between each exposure scenario (firsthand and secondhand) and individual puff volume-based groups (e.g., low, medium, high) as further detailed in the results.
For each of the exposure scenarios, generated emissions within the exposure chamber were collected on Tenax ® TA (300 mg) 60/80 mesh sorbent tube and DNPH cartridge for VOC and aldehyde analysis, respectively. Both VOC and aldehyde samples were collected for 60 min at controlled flow rates of 0.2 LPM and 0.5 LPM, respectively, by using mass flow controllers (MFC; Alicat Scientific, Tucson, Arizona). Collected VOC and aldehyde samples were analyzed with thermal desorption gas chromatography/mass spectrometry (TD-GC/MS) and high-performance liquid chromatography (HPLC) at an ISO/IEC 17025 accredited indoor air quality laboratory. The data collected from the TD-GC/MS were evaluated to ensure meeting the guidelines set for identification for the limit of quantification (LOQ).
Participants were instructed to refrain from consumption of food and liquids (aside from water) for 3 h prior to their visit. Upon arrival to the facility, participants were asked to rinse their mouth with water and after 10 min and to provide a saliva sample via passive drool technique. A total of 2 mL of saliva was collected from each participant and stored in labeled sterile cryovials with the addition of Halt™ Protease and Phosphatase Inhibitor cocktail (1% v/v, ThermoScientific, Waltham, MA) and stored at − 80°C for downstream untargeted high-resolution HPLC–MS metabolic analysis as previously described. To process for metabolic profiling, 750 μL of cold methanol was added to 300 μL of saliva. Samples were and centrifuged at 16,000×_g_ for 8 min followed by collection of supernatants. Vacuum concentrated and dried fractions of saliva samples were reconstituted in 60 μL of HPLC diluent (95% water, 5% acetonitrile, 0.1% formic acid) followed by sonication for 5 min and centrifugation at 16,000×_g_ for 8 min. The supernatants were then randomly transferred to a HPLC autosampler vial and run at one time. Separations were performed on an Agilent 1290 system (Palo Alto, CA), with a mobile phase flow rate of 0.45 mL/min. The metabolites were assayed using a Waters HSS T3 column (1.8 µm, 2.1 × 100 mm), where the mobile phases were A (0.1% formic acid in ddiH 2 O) and B (0.1% formic acid in acetonitrile). The calculated relative standard deviation for retention time was 0.27%; for area 8.4%, and for height 10.9%, all within the accepted maximum tolerance of relative standard deviation of 30% for metabolomic studies using mass spectrometry. These metabolite alterations were previously assessed in a separate companion study to evaluate ENDS-induced oral health effects. In our current study these metabolite alterations have been reanalyzed via a distinct bioinformatic approach (described below) to examine exposure and inhalation effects associated with ENDS usage.
Peak deconvolution, integration, and statistical analysis were performed within MS-DIAL (version 4.7). Peak annotations were performed either using the MSP metabolomics MS/MS library, based on authentic standards or made against Kyoto Encyclopedia of Genes and Genomes and the Human Metabolome Database, with the presence of primary ion enforced. Mass tolerances were 0.005 Da for MS1 and 0.01 Da for MS2 or 10 ppm maximum tolerance in _mz_. During alignment, a peak had to be present in at least 60% of at least one sample group and with a height greater than 5× of the corresponding peak in the blanks to be retained. xmsPANDA was used for feature selection based on _p_ values calculated using one-way ANOVA (_limma_) and linear regression model (_lmreg_) based on puff volume. Significant features with _p_< 0.05 from one-way ANOVA and linear regression model were then selected for one-way hierarchical cluster analysis (HCA) based on Pearson’s correlation coefficient as distance metric and Ward´s linkage was used as the method to create HCA. Partial least squares-discriminant analysis (PLS-DA) was also performed using MetaboAnalyst. Significant metabolic features (_limma p_< 0.05 and _lmreg p_< 0.05) were selected for metabolic pathway enrichment analysis using mummichog 2.6.0. This statistical approach has been utilized in previous evaluations employing untargeted metabolomics.
Annotated significant metabolic features (_limma p_< 0.05 and _lmreg p_< 0.05) were subject to univariate and multivariate receiver operating characteristic (ROC) analysis using MetaboAnalyst. Metabolic features were first normalized by median, log 10-transformed, and auto scaled. Monte-Carlo cross validation (MCCV) with balanced sub-sampling was employed to generate ROC curves. Area under ROC curve (AUROC), T-statistics with false discovery rate (FDR), and Log 2 fold change (FC) were calculated. Linear support vector machine (SVM) was used for model classification and feature ranking. By default, two-thirds of the participants were used for feature selection and model evaluation, and the remaining were used for model validation and feature selection. Subjects were assigned to either never established ENDS user or ENDS user group for model testing and validation. Selected features were also used to reconstruct a genome-scale metabolic network using OmicsNet 2.0. Recon3D database was used for metabolic reconstruction.
Top ten VOCs commonly shared by firsthand and secondhand emissions with at least 70% non-missing values were selected for further analyses. To further investigate the metabolic association across different data domains, network analyses were performed for saliva metabolome, selected VOCs, particulate matter, and puff topography (puff volume) from the same exhaled or device samples (n = 10) using xMWAS based on the partial least-squares (PLS) regression method. Correlation threshold was set at |r| > 0.64 and _p_< 0.05 for community detection. Network was visualized using Cytoscape (version 3.9.1) with Perfuse Force OpenCL layout. Metabolic features identified from network analysis were further used for pathway enrichment analysis using mummichog 2.6.0.
Consumer Exposure Model (CEM) developed by the U.S. EPA was employed to calculate the risk assessment from the Acute Dose Rate (ADR), Chronic Dose Rate (CDR), and Lifetime Average Daily Dose (LADD) for both firsthand and secondhand exposure scenarios. The ADR is observed to be exposure where the inhalation rate for the individual occurs during short-term exposure, 24 h/day, scenarios while the CDR uses an inhalation rate based on continued long-term exposure, 24 h/day/year, scenarios. Based on the finding from a recent literature, an average daily puff number used by ENDS user was inputted as 132 puffs in this analysis. Furthermore, the risk calculation for firsthand utilized near-field exposure zone information, to obtain an accurate estimation of the dosage that an active ENDS user is experiencing when using the device. Top ten VOCs with at least 70% non-missing values were selected to calculate the ADR and the CDR for both firsthand and secondhand exposure scenarios. The risk assessment was performed for several age groups (adult, youth, child and infant) represented by the average body weight and inhalation rates for each group. The room size and airflow exchange rate were set to 50 m 3 and 0.45 (per hour), respectively, based on environmental parameters from the CEM. Cancer risk (CR) was calculated for formaldehyde and acetaldehyde using CR = LADD × SF (Eq.1), where SF is the slope factor obtained from U.S. EPA’s Integrated Risk Information System (IRIS). SF for formaldehyde and acetaldehyde is 0.0455 (mg/kg bw/d)-1 and 0.0077 (mg/kg bw/d)-1, respectively.
All the exposure scenarios include at least three replicates throughout the study. A principal component analysis (PCA) was carried out for all puffing topography profiles of ENDS users to determine the variable most representative of exposure. The results of which suggest puff volume as the mediating factor, as described in further detail below. A paired t-test was performed to compare participle emissions and deposition results between firsthand and secondhand exposures. Further, one-way ANOVA with post-hoc testing using Fisher’s LSD was performed to determine differences across puff volume groups. Sample size and replicates are noted in each figure and table, results are presented as means ± standard error with _p_< 0.05 considered statistically significant. All statistical analyses were carried out using GraphPad Prism (version 9.5.0).
This human subject study was approved by Advarra Institutional Review Board (IRB Protocol No. H22081).
The general study participant characteristics (9 never established ENDS users and 10 established ENDS users) are summarized. An approximately equal proportion of male and female participants were recruited. No significant difference between the ages of never established ENDS users and ENDS users was detected.
Puffing topography assessments were carried out in 10 recruited established ENDS users. To determine a potential mediator of exposure, PCA analysis was performed. Puff volume represented the highest variance and exhibited a strong correlation compared to other variables. Based on these findings, the 10 established ENDS users were stratified into three groups based on their puff volume; low (n = 3; < 90 mL/puff), medium (n = 4; 90–150 mL/puff), and high (n = 3; > 150 mL/puff). Details of average combined puffing topography profiles for low, medium, and high puff volume groups are described along with comprehensive details of their devices used for firsthand and secondhand exposure scenarios. It is important to note that device format or brand, participant sex, or age were not taken into consideration in the determination of puff volume-based stratification. This stratification strategy was employed in a companion study investigating oral health effects of ENDS.
Particle concentrations and sizes were monitored for all three exposure scenarios (baseline, firsthand, and secondhand). The variations in number of particle emission, mass concentration, GMD, and particle size distribution between low (n = 3), medium (n = 4), and high (n = 3) groups from each scenario are summarized. The level of particle number emission of firsthand emissions was significantly greater compared to secondhand (_p_< 0.01) and minimal levels were present in baseline emissions (_p_< 0.0001). Moreover, an increasing trend in number emission was noted across low, medium, and high groups in both firsthand and secondhand scenarios, although significance was only noted within firsthand groups. Likewise, increased level of mass concentration was also noted in the medium and high puff volume groups compared to the low puff volume group. Furthermore, a significant particle size reduction between firsthand and secondhand emissions was noted, especially in low (_p_< 0.05) and medium (_p_< 0.01) groups, although no significance was noted between firsthand and secondhand particle sizes in the high group.
The particle distribution of firsthand and secondhand scenarios was compared based on their puff volume. In the low group comparison, firsthand showed its highest peak at 79.86 nm, while secondhand showed its two peaks at 20.35 and 36.19 nm. Both firsthand and secondhand in the low group were positively skewed towards more fine particles according to SMPS data, while OPS data showed a peak in firsthand at 375 nm. Similarly, in the medium group comparison, firsthand showed its major peak at 65.52 nm, second peak at 186.01 nm with the highest peak at 550 nm. Unlike the secondhand low group, secondhand medium group showed its highest peak at 42.55 nm with the second highest peak at 116.52 nm. Within the high puff volume group comparison, firsthand also showed a bimodal distribution with one peak at 71.69 nm and another at 375 nm. Unlike other puff volume groups, in the high group comparison, particles emitted from secondhand only showed a single highest peak at 38.89 nm.
Overall, an increased level of number and mass emission was noted as the puff volume increased. Across different puff volume-based groups in both firsthand and secondhand, the bimodal size distribution was noted with a positively skewed trend toward more ultrafine particles. In particular, the particles ranging from 10 to 100 nm showed an increased level of particle number concentrations in higher puff volume groups (medium and high) compared to the low group, suggesting an increased level of PM exposure with greater puff volume. Additionally, OPS data consistently showed the highest peak in firsthand vape in the range from 300 to 600 nm, which could be an artifact due to the water content of the aerosols increasing light scattering. Furthermore, firsthand showed a slight shift towards bigger-sized particles with greater puff volume, while secondhand remained within the range between 10 and 100 nm size.
The MPPD-calculated inhalation dosimetry results (total deposited mass, total deposited surface area, total deposited dose, deposition fraction) of firsthand and secondhand are summarized. Overall, increased levels of total deposited mass and surface area were noted in firsthand compared to secondhand, while both groups showed an increase with greater puff volume. However, the difference between firsthand and secondhand was not significant in any of puff volume-based groups. Likewise, the MPPD-calculated deposited dose was also higher in firsthand compared to secondhand, especially for the low (_p_< 0.05) and medium (_p_< 0.01) groups, showing an increase with greater puff volume. However, a greater total deposition fraction was noted in secondhand than firsthand across all puff volume-based groups with the significance noted for medium (_p_< 0.0001) and high groups. Specifically, tracheobronchial (TB) and pulmonary alveoli (P) regions showed a greater deposition fraction compared to head region in all puff volume-based groups. In particular, secondhand medium group showed significantly higher regional deposition fraction in TB (_p_< 0.0001) and P (_p_< 0.0001) compared to firsthand medium group.
The firsthand and secondhand emissions were cross-referenced to identify the common VOCs between both cases. The top 20 VOC emissions were collected by identifying the VOCs that are commonly shared by user exhalation and device emission and have at least 30% non-missing values. The top 20 VOCs have been observed as chemical components of the e-liquids or by-products from their decomposition and consumption from the vaporization process of the devices. In general, the concentration of VOCs from firsthand emission increases as the puff volume increases, especially in 1,2-propanediol, hexanal, nonanal, and toluene, while secondhand showed an increase in decanal, formaldehyde, 3-piperidinol, 1-ethyl-, 1,2-propanediol, and acetaldehyde. The secondhand emissions are much less than the firsthand emissions for almost all the VOCs, except for a few bio-active VOCs.
Over 50,000 metabolic features were yielded from Agilent 6545 Q-TOF mass spectrometer with positive ionization mode after data filtering using xmsPANDA. One-way ANOVA using _limma_ showed 10,095 features were significantly altered (_p_< 0.05). Additionally, linear regression analysis using _lmreg_ revealed 2096 significant features (_p_< 0.05), with 1160 positively associated and 916 negatively associated. A total of 263 metabolic features were finally selected by identifying the overlapping features from _limma_ and _lmreg_ analyses. One-way HCA of 263 metabolic features selected from one-way ANOVA and linear regression analyses clearly showed two clusters, with 53 metabolites in cluster I (negative association) and 210 metabolites in cluster II (positive association). PLS-DA of selected metabolic features showed a clear separation between ENDS users and never established ENDS users, as well as a noted trajectory from low to high puff volumes among ENDS users. Together with the one-way HCA, PLS-DA demonstrated the vaping behavioral response was largely driven by ENDS participants with a high puff volume. Pathway enrichment analysis was carried out using _mummichog_ for the 263 selected metabolic features. 263 significant features were putatively mapped to 15 unique empirical compounds, resulting perturbation in major metabolic pathways include energy metabolisms [pyruvate (_p_ = 2 × 10–5), glycolysis and gluconeogenesis (_p_ = 0.01)], amino acids [glycine, serine, alanine and threonine metabolism (_p_ = 0.0075), β-alanine metabolism (_p_ = 0.024)], and lipid regulation [carnitine shuttle (_p_ = 0.014), propanoate metabolism (_p_ = 0.019), prostaglandin formation (_p_ = 0.041), glycerophospholipid metabolism (_p_ = 0.047)]. Noticeably, mean saliva levels of L-palmitoylcarnitine, lactate, methylglyoxal, and 6-hydroxyhexnaoate, were significantly lower than never established ENDS users, while leukotriene D4 and 3-hydroxyquinine had higher mean levels in ENDS users with an increasing trajectory from low to high puff volumes.
To evaluate the potential application of those 15 metabolic features with identified empirical compound ID as biomarker for ENDS users, we performed univariate and multivariate ROC analyses using MCCV. Methylglyoxal, butyrate, 4-hydroxyphenylacetaldehyde, and 6-hydroxyhexanoate were the top 4 metabolites with AUROC > 0.9 and mean importance measure > 2.9. MCCV generated a model with two features with an AUC of 0.95 (95% CI 0.79–1) and a prediction accuracy of 84.4%. As the number of features increase, the better performance of MCCV was observed, with highest AUC of 1 (95% CI 1–1) and prediction accuracy of 91.3% when top 10 metabolic features were used for classification.
We next utilized the 15 metabolic features to reconstruct genome-scale metabolic network. One out of two subnetworks that had disease correlations were identified. A total of 10 metabolites were associated in the identified subnetwork, including pregnenolone sulfate, 4-hydroxyphenylacetaldehyde, methylglyoxal, butyrate, l-palmitoylcarnitine, cholesterol sulfate, leukotriene D4, taurolithocholate, retinoate, and lactate. Metabolite-disease network analysis showed several disease-related pathways following the genome-scale reconstruction. Gastric cancer and lung diseases including asthma (_p_ = 0.000187), non-small cell lung cancer (_p_ = 0.0178), and small cell lung cancer (_p_ = 0.00781) were identified to be directly associated with the input queries of metabolites. Other disease types including sensory dysfunction, neurological diseases, endocrine and metabolic diseases, and circulatory dysfunctions were associated with the genes from the genome-scale metabolic network reconstruction.
Integrative network analysis was conducted to further investigate the potential correlation among four different datasets, including puff topography (puff volumes), particle characterization, VOCs emission, and saliva metabolome. A total of 5822 metabolic features were identified and clustered into seven clusters. Puff volume was clustered into cluster 7, with 1,2-propanediol, hexanal, nonanal, toluene from firsthand emission, and 1238 metabolic features. Pathway enrichment analysis using the 1238 metabolic features from cluster 7 generated eight pathways including _N_-glycan biosynthesis (_p_ = 0.012), porphyrin metabolism (_p_ = 0.015), glycosphingolipid biosynthesis (_p_ = 0.02), bile acid biosynthesis (_p_ = 0.022), leukotriene metabolism (_p_ = 0.031), prostaglandin formation (_p_ = 0.031), blood group biosynthesis (_p_ = 0.031), and arachidonic acid metabolism (_p_ = 0.034). The pathways associated with particle characterization are shown. The mass concentration and size of the particle emitted from firsthand exposure were clustered into cluster 5, with the acetaldehyde from the secondhand exposure and 564 metabolic features. Pathway enrichment analysis revealed two pathways, including alkaloid biosynthesis (_p_ = 0.007) and squalene and cholesterol biosynthesis (_p_ = 0.012) were significantly associated. The total particle count from firsthand emission was clustered into cluster 2, with decanal and heptanal also from firsthand emission and 2082 metabolic features. A total of nine pathways, mainly focused on lipid-regulated pathways, were identified. The mass concentration, particle size, and total particle counts from secondhand emission all grouped into cluster 6, with 1-ethyl-3-piperidinol from firsthand and secondhand emission, as well as 637 metabolic features. Seven pathways, including major vitamin and cofactor metabolisms and two lipid-regulated pathways, were significantly associated with cluster 6.
The top ten VOCs shared between firsthand and secondhand emissions were distributed into all seven clusters. Toluene, decanal, nonanal, hexanal, xylenes, and heptanal had similar total abundance across firsthand and secondhand emissions. 1-Ethyl-3-piperidinol, 1,2-propanediol, acetaldehyde, decanal, formaldehyde had higher abundance in firsthand emission. Heptanal from firsthand emission had highest centrality score, followed by decanal from secondhand and firsthand emission. Pathway enrichment analysis using the metabolic features from each cluster indicate xylenes and toluene had the most differential association between firsthand and secondhand exposure. Xylenes and toluene from firsthand exposure were mainly associated with vitamins and cofactors, lipid regulation, and glycans, whereas in secondhand exposure they were only associated with sialic acid metabolism. 1,2-Propanediol, nonanal, hexanal, and heptanal had more pathway associated from secondhand emission and were in the same cluster. 1-Ethyl-3-piperidinol and decanal from firsthand and secondhand emission were in the same cluster and had the exact pathways associated. Overall, majority of the selected VOCs were associated with inflammatory response related pathways, for example, arachidonic acid, leukotriene, prostaglandins, and bile acids, metabolisms.
Risk assessment using the CEM model revealed almost all ten selected VOCs from secondhand vaping emissions could pose significantly higher acute exposure dose and chronic exposure dose, and therefore higher health risks to younger age groups. Importantly, a puff behavior-dependent exposure risk was identified, showing increasing ADR and CDR as the puff volumes increase. Decanal, formaldehyde, 1-ethyl-3-piperidinol, and 1,2-propanediol from secondhand emission were the major VOCs driving the puff volume-dependent pattern. Higher exposure risk was found associated with high puff volume group for acetaldehyde. Notably, we also identified a clear puff volume dependent exposure risk for several VOCs from firsthand emission, including 1,2-propanediol, hexanal, nonanal, and toluene, which is consistent with the puff volume dependent increase in each corresponding VOC. While the majority of VOC emission and exposure risk were decreased in secondhand exposure scenario in comparison to firsthand emission, several VOCs had higher exposure and risk level in secondhand exposure. Formaldehyde exposure and estimated health risks were substantially increased in secondhand exposure scenario for medium and high puff group. Acetaldehyde exposure and estimated health risks were increased in secondhand exposure scenario for high puff group, whereas xylenes exposure and associated health risks were increased in secondhand exposure scenario for medium puff group. Toluene exposure and estimated health risks increased in secondhand exposure scenario for low and medium puff group. Based on the modeled lifetime average daily exposure levels, the cancer risk of formaldehyde from secondhand emissions showed a puff volume-dependent increase with 1.13 × 10–5 for medium puff group and 1.03 × 10–5 for high puff group, both were significantly higher than 8.51 × 10–6 from low puff volume group. The trend was opposite for formaldehyde from firsthand emission with a cancer risk of 8.37 × 10–5 for low puff group, 4.5 × 10–6 for medium puff group, and 6.83 × 10–6 for high puff volume group. Lower cancer risks were found with acetaldehyde, with highest of 4.97 × 10–6 for low puff volume group from firsthand emission and 1.35 × 10–6 for high puff group from secondhand emission. Leukotriene D4 level in saliva increased with higher formaldehyde level and higher acetaldehyde level, though not statistically significant. In contrast, retinoate level in saliva decreased with higher formaldehyde level and higher acetaldehyde level.
In this study, a comprehensive investigation of potential secondhand vaping exposure risks in an indoor environment was performed based on the physical and chemical characterization of the emitted aerosols along with its impact on various metabolic processes. Current studies of potential impact of secondhand vaping on indoor air quality are primarily focusing on compositional analysis of the emission, while this study is one of the few that demonstrates real-life personal exposure scenarios by conducting a comparative analysis between firsthand and secondhand with puffing topography profile-matched data and metabolomic outcomes associated with each exposure scenario. Moreover, the calculated vital capacity (VC) and forced vital capacity (FVC) did not show significant differences among the low, medium, and high puff volume groups. The comparable lung functions across different puff
Bookmark
Daniel Féau processes personal data in order to optimise communication with our sales leads, our future clients and our established clients.
This site is protected by reCAPTCHA and the Google Privacy Policy and Terms of Service apply.