A. R.
Iskandar‡
a,
B.
Titz‡
a,
A.
Sewer‡
a,
P.
Leroy‡
a,
T.
Schneider
a,
F.
Zanetti
a,
C.
Mathis
a,
A.
Elamin
a,
S.
Frentzel
a,
W. K.
Schlage
b,
F.
Martin
a,
N. V.
Ivanov
a,
M. C.
Peitsch
a and
J.
Hoeng
*a
aPMI R&D, Philip Morris Products S.A. (part of the Philip Morris International group of companies), Quai Jeanrenaud 5, CH-2000 Neuchâtel, Switzerland. E-mail: julia.hoeng@pmi.com; Fax: +41 (58)242 2811; Tel: +41 (58)242 2214
bBiology consultant, Max-Baermann-Str. 21, 51429 Bergisch Gladbach, Germany
First published on 29th May 2017
Systems biology combines comprehensive molecular analyses with quantitative modeling to understand the characteristics of a biological system as a whole. Leveraging a similar approach, systems toxicology aims to decipher complex biological responses following exposures. This work reports a systems toxicology meta-analysis in the context of in vitro assessment of a candidate modified-risk tobacco product (MRTP) using three human organotypic cultures of the aerodigestive tract (buccal, bronchial, and nasal epithelia). Complementing a series of functional measures, a causal network enrichment analysis of transcriptomic data was used to compare quantitatively the biological impact of aerosol from the Tobacco Heating System (THS) 2.2, a candidate MRTP, with 3R4F cigarette smoke (CS) at similar nicotine concentrations. Lower toxicity was observed in all cultures following exposure to THS2.2 aerosol compared with 3R4F CS. Because of their morphological differences, a smaller exposure impact was observed in the buccal (stratified epithelium) compared with the bronchial and nasal (pseudostratified epithelium). However, the causal network enrichment approach supported a similar mechanistic impact of CS across the three cultures, including the impact on xenobiotic, oxidative stress, and inflammatory responses. At comparable nicotine concentrations, THS2.2 aerosol elicited reduced and more transient effects on these processes. To demonstrate the benefits of additional data modalities, we employed a newly established targeted mass-spectrometry marker panel to further confirm the reduced cellular stress responses elicited by THS2.2 aerosol compared with 3R4F CS in the nasal culture. Overall, this work demonstrates the applicability and robustness of the systems toxicology approach for in vitro inhalation toxicity assessment.
The pathophysiology of CS-induced effects in the human lung can be studied by obtaining biopsy samples. However, such an approach frequently involves invasive procedures (e.g., bronchoscopy) and is limited to small amounts of samples. Although technologies have improved tremendously, for example, fibered confocal microscopy (“optical biopsy”) can be conducted during bronchoscopy,6,7 the available approaches do not allow to decipher the mechanisms of early toxicity impact on the field of injury.17 Alternatively, toxicological testing frequently relies on the use of animal models. However, strong efforts are being made by the scientific community—and further enforced by various regulatory bodies—to implement the 3Rs principle: to reduce, refine, and replace the use of animals.8 These efforts are equally motivated by the issue of cross-species translatability, in which results obtained from animal studies do not necessarily reflect the potential toxic impact in humans.
Organotypic culture models, which can be derived from human primary cells, are an alternative to investigate inhalation toxicity in human tissues. Unlike submerged monolayer cultures of respiratory epithelial cells, organotypic epithelium cultures are grown at the air–liquid interface and can therefore be directly exposed to CS, aerosols, or nanoparticles on the apical side.9,10 Following exposure to CS, tissue injury occurs not only in the distal lung, but also in the aerodigestive tract. In the oral cavity, CS exposure has been linked to inflammation and cancer.11–14 In response to CS exposure, the methylation and gene expression profiles of oral epithelial cells were similar to those of bronchial epithelial cells.15 In the nasal epithelium, the alteration of genes involved in xenobiotic metabolism—a biosensor of exposure—following smoking was reported to resemble that in the bronchial epithelium.15,16 Increased levels of the cytochrome P450 (CYP) 1A1 and 1B1 genes were also observed in both the nasal and bronchial epithelium of smokers, further supporting the relevance of nasal tissue as a surrogate lung airway epithelium.17 Various human organotypic culture models of the aerodigestive tract are available. Models derived from nasal or bronchial epithelial cells form a fully differentiated epithelium composed of basal, ciliated, and mucus-secreting goblet cells,18,19 while the models derived from buccal epithelial cells develop into a stratified non-keratinizing epithelium.20 Studies have demonstrated that the responses of these human organotypic cultures following CS exposure resemble those observed under in vivo situations.21–24 Recently, we have reported that these human organotypic culture models were useful and relevant to assess the effects of aerosol generated from the Tobacco Heating System (THS)2.2, a candidate modified risk tobacco product (MRTP), as compared with CS.18–20
Here, we report on a systems toxicology meta-analysis of the biological impact of exposure to THS2.2 aerosol and 3R4F CS on three types of organotypic cultures of aerodigestive epithelia (leveraging data from our recent publications: buccal,20 nasal,18 and bronchial19). This work further explores whether the “field of injury” concept15—CS exposure affects similar mechanisms across the entire tissues lining the aerodigestive tract—can be reproduced in these in vitro studies. In the context of such a meta-analysis, this work demonstrates how the systems toxicology framework enables a comprehensive toxicological and mechanistic assessment across tissues and exposure conditions. Overall, the data show that 3R4F CS exposure triggers similar mechanisms in all three tissues, demonstrating an in vitro observation of the “field of injury.” This work also reveals that the reduced levels of harmful and potentially harmful chemicals in the THS2.2 aerosol compared with the 3R4F smoke25 translate into reduced perturbations of critical biological processes including inflammation and cellular stress that were consistently affected by 3R4F CS exposure across the three in vitro human organotypic airway epithelium cultures.
Fig. 2 A series of in vitro studies using human organotypic epithelium cultures. Organotypic culture models recapitulating the human aerodigestive tract lining the “tissue of injury” fields (buccal, bronchial, and nasal) were exposed to 3R4F CS or THS2.2 aerosol at similar nicotine concentrations in an Exposure System (Vitrocell 24/48®). An illustration of the air–liquid interface organotypic culture located in the Base Module of the exposure system is shown. Exposure characterization throughout the study period included measurements of nicotine concentrations (in the CS/aerosol) and deposited carbonyl concentrations (in the phosphate-buffered saline-filled Cultivation Base Module). Various biological endpoints were measured at various post-exposure time points as illustrated. The detailed experimental procedures have been described previously.18–20 ●Three independent exposure-runs were conducted for each item (3R4F and THS2.2); except for those using bronchial cultures in 2014. *Dilution refers to the percentage 3R4F smoke or THS2.2 aerosol diluted with air in the Dilution/Distribution Module of the Exposure System. †Nicotine concentration (mg L−1) refers the corresponding concentration to the specific dilution of smoke/aerosol determined by trapping the diluted smoke/aerosol in the EXtrelut® 3NT column. §The nasal organotypic cultures were reconstituted from the primary nasal epithelial cells of 30 year-old non-smoker male; buccal organotypic cultures were reconstituted from the primary buccal epithelial cells of 46 year-old non-smoker male; and bronchial organotypic cultures were reconstituted from the primary bronchial cells of 28 year-old non-smoker male (except for the first two experimental repetitions in which the bronchial cultures were reconstituted from 23 year-old non-smoker male). NA: not available. |
Soste et al. have demonstrated how targeted proteomic analysis of sentinel proteins, i.e. biological marker proteins, can efficiently capture the activation state of up to 188 biological processes in baker's yeast.40 Guided by this idea, a protein marker panel (sentinel panel) was defined to cover the major effect categories relevant to organotypic exposure studies: xenobiotic metabolism, oxidative stress, metabolic adaptations, unfolded-protein response (UPR), tissue composition changes, barrier function, and senescence (Fig. 9). The list of targeted proteins for these effect categories together with details of the spiked-in peptides is presented in ESI Table 1.†
Total protein was extracted from the cultures using a commercially available sample preparation kit (Biognosys AG, Schlieren, Switzerland). The cells were incubated in denaturing buffer and disrupted by ultrasound treatment. The protein concentration was determined using a Pierce 660 nm protein assay (Thermo Scientific). Fifty microgram of protein was subjected to protein reduction, alkylation, and digestion as described in the manual. Prior to analysis, samples were purified using C18 reversed phase and solid phase extraction plates (Waters). After resuspension in LC-Buffer A (1% acetonitrile, 0.1% formic acid), stable-isotope labeled reference peptides were spiked into the sample for each target of interest. About 0.5 μg of total protein was injected for analysis onto a 15 cm C18 reversed-phase column and analyzed by liquid chromatography coupled with tandem mass spectrometry using an Easy-nano LC 1000 instrument connected online to a Q-Exactive Plus (System 1) or Q-Exactive HF mass-analyzer (System 2, both from Thermo Scientific).
On system 1, peptides were separated using a PepMap RSLC C18 column (50 μm × 15 cm, 2 μm particle size, 100 Å pore-size, Thermo Scientific) with a flow rate of 200 nL min−1 in a gradient starting with 5% LC-Buffer B (95% acetonitrile, 0.1% formic acid) to 28% LC-Buffer B over 50 min followed by a 10 min column wash step at 100% LC-Buffer B. The Q-Exactive Plus mass spectrometer was operated in retention time scheduled PRM mode with a resolution of 17.500, an AGC target value of 1e6, a maximum injection time of 30 ms and a precursor isolation window of 1.2 m/z.
On system 2, peptides were separated using an Acclaim PepMap RSLC C18 column (75 μm × 15 cm, 2 μm particle size and 100 Å pore-size, Thermo Scientific) with a flow rate of 300 nL min−1 in a gradient of 5% LC-Buffer B to 30% LC-Buffer B over 30 min followed by a 10 min column wash step at 100% LC-Buffer B. A Q-Exactive HF mass spectrometer was operated in retention time scheduled PRM mode with a resolution of 30000, an AGC target value of 1e6, a maximum injection time of 60 ms and a precursor isolation window of 1.2 m/z. For retention time scheduling, iRT peptides (Biognosys AG, Schlieren, Switzerland) were used (Escher et al., 2012)95.
Raw files from the PRM acquisition were analyzed with SpectroDive (version 7.5, BiognoSYS AG). Ion chromatograms for the endogenous peptides and the corresponding stable isotope labeled reference peptides were extracted for all the measured transitions (ESI Table 1†) using the software vendor's default settings. For quantification, the area under curve (AUC) intensities of all transitions were summed and the ratios of AUC sums of the endogenous and corresponding reference peptide signals were calculated.
For the statistical analysis, a linear model was fitted for each exposure and the respective sham group, including the experimental repetition as a covariate to account for the pairing between exposure and sham groups.41 The obtained raw p-values (without empirical Bayes moderation, corresponding to a paired t-test) were adjusted across protein markers using the Benjamini–Hochberg false discovery rate (FDR) method. Differentially expressed proteins were defined as those with a FDR-adjusted p-value of <0.05.
In parallel, a gene set analysis (GSA) was performed with pathway maps from the KEGG knowledge base.42 The significance of the gene-set enrichment was assessed using a competitive null hypothesis (Q1) and a self-contained null hypothesis (Q2).43 Whereas Q1 tests for the significance of genes in the set versus those not in the set, Q2 tests for a significant difference between the conditions. With this, we expect Q2 to be more appropriate in the context of comparative toxicity assessment (e.g., to reveal a significant effect on a given gene set compared with air-exposed controls exposure), while Q1 can highlight gene-sets that dominate these responses. Here, the Q1 statistics were calculated with the Camera approach44 and the Q2 statistics with the Roast approach,45 which take the gene correlation structures into account. The resulting p-values were adjusted for multiple hypothesis testing using the Benjamini–Hochberg procedure.
The comparability across the exposure-induced biological impact was examined using clustering with two distinct network-based similarity metrics (amplitude-based and shape-based). The amplitude-based metric considers two exposure conditions with similarly high NPA values to be closely related. To calculate this metric, first all the node-level contributions to the normalized NPA values28 were concatenated across all networks; these values were normalized so that, for all networks, the maximal normalized NPA value across all exposure conditions was 1, and so that all the node contributions from not statistically significantly perturbed networks were 0; finally, the Euclidean distance matrix between the available 24 considered exposure conditions was calculated. The “shape-based” metric specifically compares the relative distributions of the node-level contributions across exposure conditions. Two exposure conditions with similar distributions of their node-level contributions (but not necessarily the same magnitude of non-normalized NPAs) are considered to be closely related. To calculate this metric, it was assumed that the normalized node contributions for any given network and exposure conditions summed to 1; then the corresponding distance matrix was calculated using the Manhattan formula, which is appropriate for comparing concatenated (normalized) distributions. For clustering, an affinity propagation-based approach was applied in two steps to identify clusters and the connections between their “exemplar elements”.46 The defined graph encodes the similarities between exposure conditions within and between the obtained clusters.
Characterization of the 3R4F CS and THS2.2 aerosol in the exposure system (Vitrocell 24/48®) included the measurements of nicotine—previously described47—in the 3R4F CS and THS2.2 aerosol throughout the studies. The measurements were done to assess the reproducibility of the nicotine concentrations for a given dilution of 3R4F CS or THS2.2 aerosol. Fig. 3A shows that reliable nicotine concentrations were achieved for a particular smoke/aerosol dilution, despite some variability (descriptive statistics results are given in ESI Table 3†). Moreover, the concentrations of carbonyls deposited in the Cultivation Base Module were also determined throughout the studies (Fig. 3A) and were found to be comparable across the exposure experiments for the particular doses of 3R4F CS and THS2.2 aerosol (descriptive statistics results are given in ESI Table 3†). The concentrations of the deposited carbonyls remained lower in the PBS samples exposed to THS2.2 aerosol compared with 3R4F CS at all the tested doses although the nicotine concentrations between THS2.2 and 3R4F CS were similar (i.e., the nicotine concentration in the low dose of 3R4F CS was comparable to the low dose of THS2.2; and the high dose of 3R4F CS was comparable to the medium dose of THS2.2). A more extensive characterization of the THS2.2 aerosol has been reported previously.25
Acute 3R4F CS exposure (at 0.13–0.15 mg nicotine per L) elicited a similar cytotoxicity profile (based on the level of adenylate kinase released into the basolateral media upon cell damage) in both the bronchial and nasal cultures: the cytotoxicity increased dose dependently and with the duration of post-exposure (Fig. 3B). Even at higher doses (nicotine concentrations of 0.32–0.51 mg nicotine per L), 3R4F CS elicited a less pronounced cytotoxicity in the buccal cultures compared with the bronchial and nasal cultures following exposure to 3R4F CS at 0.25–0.27 mg nicotine per L. The 3R4F CS-induced cytotoxicity in the buccal cultures increased only slightly (maximum ∼6%) with the post-exposure duration. This cytotoxicity measurement relies on the level of adenylate kinase in the basolateral media; therefore, its level in the apical compartment was not taken into account that could potentially underestimate the cytotoxicity measurement—in particular, for a thick stratified epithelium culture. Nevertheless, a histological assessment of the buccal cultures demonstrated some detachment above the basal cell layer, keratinization/desquamation of the epithelium, and the presence of intracellular granular structures (reported previously20). In contrast, THS2.2 aerosol exposure elicited only minimal cytotoxicity in all three cultures at all the doses and post-exposure time points tested.
An additional functional parameter was inferred from the ciliary beating functionality of the ciliated pseudostratified epithelium airway cultures (bronchial and nasal) following exposure. Mucociliary clearance—a mechanism driven by the coordinated movement (beating) of cilia to transport mucus-containing toxicants in the respiratory tract48—is an initial defense mechanism to clear inhaled toxicants.49Fig. 3C (1st row, Weighted Frequency by the power of the beating signal) shows that 3R4F CS exposure (at 0.25 mg nicotine per L) reduced the frequency of the ciliary beating in both bronchial and nasal cultures immediately after exposure. The weighted frequency remained low at the later post-exposure time points compared with the air control. Fig. 3C (2nd row, Frequency Uniformity) shows that 3R4F CS exposure (at 0.25 mg nicotine per L) interrupted the uniformity of the beating frequency in both bronchial and nasal cultures. Less uniformity of the beating frequency was observed at all post-exposure time points, which suggested that the normal propulsion of the mucus layer was disrupted following 3R4F CS exposure. Both disorganization of the cilia beating frequency and reduction in ciliary beating frequency have been reported in individuals with asthma.48 Moreover, in both bronchial and nasal cultures exposed to 3R4F CS (at 0.25 mg nicotine per L), the areas at which active ciliary beating was detected were reduced (Fig. 3C, 3rd row, Active Area). The smaller active area could be attributed to the loss (shedding) or shortening of the cilia, similarly to what has been observed in nasal biopsies of individuals exposed to high levels of air pollution.50 The power of the detected signal (log FFT power, Fig. 3C, 4th row) further confirmed that 3R4F CS exposure impacted the beating signal immediately after exposure (particularly for the 3R4F CS at 0.25 mg nicotine per L). The overall data further reveal that the cultures have the capacity to recover following exposure: the power of the beating signals eventually returned to levels comparable to the air-exposed samples (at the 24 h and 48 h post-exposure time points, and 72 h for the bronchial cultures). However, the weighted frequency, frequency uniformity, and active area remained low. The observation suggests that despite recovering (regaining power of the beating signals), the beating frequency and the coordinated ciliary movement remained compromised. Moreover, the analysis demonstrated a distinct profile of the ciliary beating function following exposure: a dose-dependent impact of 3R4F CS was observed in the bronchial cultures but not in the nasal cultures. This suggests that the nasal cultures exhibit superior defense capabilities (with regard to mucociliary clearance), which is consistent with the role of the nasal epithelium as the first lines of defense against inhaled pathogens, dusts, and irritants.51 Conversely, exposure to THS2.2 aerosol at all the tested concentrations did not impact the weighted ciliary beating frequency, the uniformity of the frequency, the active area, or the power of the beating signal. The findings were similarly observed in both bronchial and nasal cultures, further supporting the minor cytotoxicity impact of THS2.2 aerosol discussed earlier.
Fig. 4 Mechanistic investigation of the exposure impact based on the transcriptomic data. (A) Barplots showing the number of significantly differentially expressed genes (DEGs) across the exposure conditions (FDR-adjusted p-value <0.05). The heatmaps indicate the expression profiles of the top ten genes (sorted first by the number of significant conditions and then by the mean of the absolute fold-changes). The log2(fold-changes) compared with the respective air control groups are color-coded and the statistical significance level is indicated (FDR-adjusted p-value). (B) Gene set analysis (GSA) was performed with the KEGG gene-set collection using absolute log2(fold-changes) as the gene-level and the mean as the gene-set level statistics. Significance with respect to the treatment effect (Q2, compared with the air control) and dominant effects of individual gene sets (Q1) was assessed with Benjamini-Hochberg based FDR adjustment (FDR adj. p-value <0.05). The numbers of significantly up- and down-regulated gene sets for Q1 and Q2 are shown in the top panel, and the top gene sets, first sorted by the number of significant conditions and then by their average absolute scores, are shown in the bottom panels. (C) The causal network enrichment approach for the analysis of the transcriptomic datasets. For each network category, the relative biological impact factor is shown in radar plots (CFA, Cell Fate; CPR, Cell Proliferation; CST, Cell Stress; IPN, Inflammatory Process Networks). The heatmaps show the network perturbation amplitudes for each network in the collection, across all conditions. Full details of the comparative analyses for all doses of the exposure are given in the ESI Fig. 2.† |
All three organotypic cultures (buccal, bronchial, and nasal) demonstrated a pronounced differential gene expression response upon 3R4F CS exposure at all post-exposure time points (Fig. 4A, top panels). The largest number of differentially expressed genes (DEGs) was observed at 72 h post-exposure to 3R4F CS in the buccal cultures; the largest number of DEGs in the bronchial cultures was observed at 24 h post-exposure following 3R4F CS; whereas the maximum number of DEGs was observed at 4 h post-exposure to 3R4F CS in the nasal cultures. These data suggest a different kinetics of the exposure response across the three cultures. At the similar nicotine concentrations, THS2.2 aerosol elicited a less pronounced change in the DEGs. The maximum number of DEGs was observed 4 h post-exposure to THS2.2 (551 DEGs) in the nasal cultures. In all three cultures, the THS2.2 aerosol-induced DEGs were mostly limited to the earliest post-exposure time point (4 h), suggesting a lower and more transient gene expression response to THS2.2 aerosol exposure than 3R4F CS exposure.
The top ten genes—ranked first by the number of significant conditions and second by their fold-changes—were all involved in either xenobiotic metabolism (CYP1A1, CYP1B1) or the oxidative stress response (TXNRD1, SRXN1, NQO1, HMOX1, SLC7A11, GCLM, LOC344887, MAFF) (Fig. 4A, bottom panels). While 3R4F CS induced persistent changes in these genes at all post-exposure time points, THS2.2 aerosol exposure elicited only transient changes in the expression of these genes, mostly limited to the 4 h post-exposure time point.
Gene-set analysis (GSA) enables aggregation of individual gene responses and links them to specific biological processes. Fig. 4B (upper panels) summarizes the number of significantly affected gene sets of the KEGG gene-set collection.42Fig. 4B (lower panels) shows the profiles for the top 10 affected gene sets. Overall, the GSA-based analysis showed a response pattern similar to that of the DEG profiles: 3R4F CS exposure elicited prominent biological responses across all three cultures at all post-exposure times; in contrast, THS2.2 aerosol resulted in a more transient and weaker response (mostly at 4 h post-exposure time point). The large number of diverse gene sets linked to the 3R4F CS exposure suggests an extensive biological impact of 3R4F CS (the entire GSA results are reported in ESI Fig. 1†). The most consistently and most sturdily affected gene sets across the exposure conditions pointed to the activation of xenobiotic (Metabolism of Xenobiotics by Cytochrome P450), oxidative stress (Glutathione Metabolism), and pro-inflammatory (NF-kappa B Signaling Pathway) responses of the exposed cultures. The top 10 gene sets illustrate the main challenge of GSA, whereby unspecific diverse gene sets are uncovered. Such gene sets often are not directly pertinent to the biological context under investigation. For example, the Pertussis, Legionellosis, and Chemical Carcinogenesis gene sets were also identified in this study. Despite its name, Chemical Carcinogenesis, the genes within this gene set encode xenobiotic metabolism enzymes. In principle, these xenobiotic metabolizing enzymes can also contribute to the transformation of pro-carcinogens into carcinogens,53 but are not directly indicative of an active carcinogenesis process.
As compared with GSA, the causal network enrichment approach offers a more targeted toxicological assessment, in which the relevant biological context is considered [Fig. 1]. In the present work, the exposure-induced perturbation of 29 causal biological network models was analyzed and compared across the buccal, bronchial, and nasal epithelium cultures (with the exception of the Epithelial Mucus Hypersecretion network model, which was not used for the analysis of the buccal cultures as the stratified epithelium model is not comprised of mucus-secreting cells). The analysis demonstrated a similar 3R4F CS-induced perturbation in all four functional network categories for all the three cultures: Cell Stress (CST), Inflammatory Process Network (IPN), Cell Fate (CFT), and Cell Proliferation (CPR) (Fig. 4C, radar plots). The plots demonstrate a pronounced impact of 3R4F CS exposure on the buccal, bronchial, and nasal cultures: the greatest impact was mostly observed at the earliest post-exposure time point and was declining thereafter. However, the buccal cultures’ inflammatory response (modelled in the IPN) following the highest 3R4F CS-induced perturbation was observed at 72 h post-exposure.
Compared with the impact of 3R4F CS, THS2.2 aerosol exposure elicited a lower impact on these network collections. Because the network models are organized hierarchically, further mechanistic investigations can be done at the level of individual networks (Fig. 4C, heatmap panels). For example, in the Cell Stress category, 3R4F CS exposure perturbed the Xenobiotic Metabolism Response, Oxidative Stress, and NFE2L2 Signaling networks in all three cultures. Although THS2.2 aerosol also induced significant perturbations in these networks, the amplitudes of the impact remained well below those observed upon 3R4F CS exposure. The kinetics of the 3R4F CS-induced impact across the three cultures were different. Nasal and bronchial cultures had the highest perturbation at 4 h after 3R4F CS exposure and the perturbations decreased with the duration of post-exposure. Similarly, the buccal cultures were initially perturbed 4 h after 3R4F CS exposure but the impact was lower than that observed in the bronchial and nasal cultures. Thereafter, the perturbations in the buccal cultures generally persisted until the later post-exposure time points—or even exacerbated—unlike those in the bronchial and nasal cultures, e.g., for the Epithelial Innate Immune Activation and Tissue Damage networks.
The analyses of the transcriptomic data based on the gene, gene-set, and causal network level yielded consistent exposure response profiles with increasing statistical and biological confidence. Each of these analyses provided different levels of mechanistic insights. The DEG analysis demonstrated a differential gene expression profile in response to exposure, especially to 3R4F CS, and could identify the most common genes differentially altered following exposure. The GSA approach could link the impact of exposure to various biological processes, including the impact on xenobiotic metabolism and oxidative stress, despite including some annotations that are less relevant in the context of the respiratory system. Finally, the causal network enrichment methodology provided a more concise quantification of the biological and functional processes that are pertinent for respiratory physiology. Together, the results showed that THS2.2 aerosol exposure induced a weaker and more transient biological impact than 3R4F CS at similar nicotine concentrations, as previously observed.18,19,54 THS2.2 aerosol exposure resulted mainly in transient adaptation processes, which is consistent with the observation that less than 10% of the toxicant yields in 3R4F CS are present in THS2.2 aerosol.25
Clustering with the amplitude-based metric clearly reflected the magnitude of the exposure impact for different exposure groups (Fig. 5A). All THS2.2 exposure impact (at all post-exposure time points for all culture types) were clustered together (Fig. 5A, i). However, this cluster also included the later impact of 3R4F CS aerosol exposure on the nasal cultures (48 h and 72 h post-exposure). This observation reflects the overall low and transient network perturbations following THS2.2 aerosol exposure and the rapid recovery of the nasal cultures from the 3R4F CS exposure discussed earlier (i.e., regarding ciliary beating functionality measures (Fig. 3C) and the causal network enrichment analysis (Fig. 4C, radar plots)). The 3R4F CS-induced perturbations at 4 h post-exposure in the bronchial cultures formed an independent cluster (Fig. 5A, ii) and those in the nasal cultures and the 3R4F CS-induced perturbations in the bronchial cultures at 24 h post-exposure formed two interconnected clusters (Fig. 5A, iii). These findings suggest that the degree of impact on the nasal cultures 4 h after 3R4F CS exposure resembles that on the bronchial cultures at 24 h post-exposure, which was consistent with Fig. 4C showing the network perturbation score. Finally, another cluster (Fig. 5A, iv) included the rest of the 3R4F CS-induced impact on the bronchial and nasal cultures at the later post-exposure time points and the 3R4F CS-induced impact on the buccal cultures at all post-exposure time points. This observation resembles the perturbation scores elicited by 3R4F CS exposure shown in Fig. 4C where less impact of exposure in general was observed for the buccal cultures compared with the bronchial and nasal cultures. The overall analysis demonstrated a tissue type specific response following exposure.
To further resolve the similarities of the exposure-induced impact, a cluster analysis was conducted based on the shaped-based metric that relies on the similar pattern of the impacted molecular entities (i.e., network nodes) in the network models, independent of the magnitude (a detailed description is given in the Materials and methods section). If the exposure-induced network perturbations cluster together in this analysis, similar molecular entities were perturbed for the same network model. In contrast, if the exposure impacts do not cluster together, different molecular entities were impacted following exposure, thus the perturbation impacted different network models. Fig. 5B, i shows that the perturbations at 4 h after exposure to 3R4F CS for all the three cultures pointed toward similar perturbed molecular entities and network models. Similarly, the THS2.2 aerosol induced impacts at 4 h post-exposure on all culture types formed one cluster (Fig. 5B, ii). The THS2.2 aerosol induced impacts at the later post-exposure time points clustered together, suggesting that the transient impact of THS2.2 aerosol exposure affecting similar molecular entities and network models. Furthermore, the shaped-based analysis in particular demonstrated tissue-specific clusters for the exposure impact at the later post-exposure time points: those of buccal (Fig. 5B, iv), bronchial (Fig. 5B, v), and nasal (Fig. 5B, vi).
In summary, the clustering analysis based on the exposure-induced network perturbations enabled a high-level (but mechanism focused) meta-analysis of the exposure responses and could distinguish the lower impact observed following THS2.2 aerosol exposure compared with the more pronounced impact of 3R4F CS (using the amplitude-based metric). In addition, the analysis could identify the similarities and differences in the cellular response across different tissues (especially, using the shape-based metric).
Fig. 6 Induction of xenobiotic metabolism and oxidative stress responses across the three organotypic epithelium cultures. (A) Correlation of the perturbation scores of the molecular entities in the Xenobiotic Metabolism Response network across the three organotypic cultures. The color of the heatmap indicates the Spearman correlation coefficients. Scatter plots for each culture pair show the correlation of the activation values of the molecular entities in the Xenobiotic Metabolism Response network. (B) Clustered gene expression matrix for genes of a xenobiotic metabolism gene set (KEGG collection,42 Metabolism of Xenobiotics by Cytochrome P450). The log2 fold-changes compared with the air control are color-coded and FDR-adjusted significance is indicated for the FDR-adjusted p-value <0.05 (x) and <0.01 (*) levels. Gene clustering based on the pair-wise correlation between the fold-changes and the clustering results are shown as a dendrogram (clusters are marked in different colors). (C) As in A, but for the Oxidative Stress network. (D) As in B, but for the oxidative stress response gene set (Reactive Oxygen Species Pathway of the hallmark collection60). |
The correlation plots revealed that the xenobiotic responses to 3R4F CS were well aligned across the three cultures. The perturbation scores of the molecular entities within the Xenobiotic Metabolism Network exhibited high Spearman correlation values across the buccal, bronchial, and nasal cultures (Fig. 6A). The highest correlation values were observed for the 4 h time points following 3R4F CS exposure. Because of the minor (and mainly transient) biological impact of THS2.2 aerosol exposure, the correlation values for the cultures exposed to THS2.2 aerosol were well correlated only at the 4 h post-exposure time point. This is consistent with the observed transient xenobiotic response following THS2.2 aerosol exposure shown in Fig. 4. Furthermore, the expression profiles of genes involved in xenobiotic metabolism (Fig. 6B) showed that genes encoding the cytochrome P450 enzymes CYP1A1 and CYP1B1 were greatly induced in response to 3R4F CS exposure. This observation was consistent in all three organotypic cultures. Moreover, the increased CYP1A1 and CYP1B1 levels only at 4 h post-exposure to THS2.2 aerosol further exemplified the earlier remark regarding the transient response upon exposure to THS2.2 aerosol. Additional xenobiotic metabolism genes—with a similar (transiently altered) profile—were assigned to the orange correlation cluster (i.e., aldehyde dehydrogenases ALDH1A3 and ALDH3A1, aldo-keto reductases such as AKR1C3, and the glucuronosyltransferase UGT1A6). Notably, the expression profiles of some genes in this orange cluster suggested differential xenobiotic responses between the buccal (stratified epithelium) and bronchial/nasal (pseudostratified epithelium) cultures; for example, the epoxide hydrolase EPHX1 gene, which was most prominently up-regulated in the buccal culture following 3R4F CS exposure. The genes within the green cluster were mostly down-regulated in the buccal as compared with the bronchial/nasal cultures; whereas the genes within the purple cluster showed the opposite trend. This analysis further showed that several paralogous genes were assigned to different clusters (e.g. for the green vs. purple cluster: GSTA4 vs. GSTA1/3, ADH5/7 vs. ADH1B/6, ADH3B2 vs. ADH3B1, and CYP2C9 vs. CYP2C8 were assigned). These results suggest tissue-specific differences in the xenobiotic metabolism response, which involve different members of these gene families.
The correlations of the perturbation scores of the molecular entities within the Oxidative Stress network showed a fairly consistent response across the three cultures following 3R4F CS exposure at all post-exposure time points. For THS2.2 aerosol exposure, high correlation values were only observed for the 4 h post-exposure time points (Fig. 6C). A correlation-based clustering of the expression profiles of oxidative stress-responsive genes revealed three main clusters (Fig. 6D): the genes within the main cluster (purple) showed a consistent induction after 3R4F CS exposure for buccal, bronchial, and nasal cultures. This cluster contains genes that represent the main branches of the oxidative stress response, namely: the glutathione system (GCLC, GCLM, MGST1, GSR), the thioredoxin system (thioredoxin reductase 1, TXNRD1), the peroxiredoxin system (PRDX1, SRXN1), quinone NAD(P)H dehydrogenase 1 (NQO1), and metabolic adjustments for NAPDH production (glucose-6-phosphate dehydrogenase (G6PD)). A similar impact on the oxidative stress response pathways was observed previously in rat lung tissue exposed to 3R4F CS56 and in human lung tissue under oxidative stress conditions.57–59 The altered gene expression in this cluster further showed that THS2.2 aerosol induced a transient oxidative stress response that was limited to the 4 h after exposure time point across all three cultures (an exception was observed for NQO1 in the nasal culture). The second cluster (green) contains several genes including the gene encoding catalase (CAT) that were significantly down-regulated only upon 3R4F CS exposure in all three cultures. This is different from a previous study, which found that CS exposure up-regulated the expression of CAT.59 Moreover, following repeated exposure to 3R4F CS, the expression of the CAT gene in human organotypic gingival cultures was found to be up-regulated.54 The third cluster of genes (orange) indicates a distinct oxidative stress response. Genes in this cluster were mainly down-regulated in the buccal (stratified epithelium) but mainly up-regulated in both bronchial and nasal (pseudostratified epithelium) cultures in response to 3R4F CS exposure.
Taken together, both 3R4F CS and THS2.2 aerosol exposure affected the xenobiotic metabolism and oxidative stress response of the three cultures. Both processes were consistently perturbed following 3R4F CS exposure across all three tissue cultures; however THS2.2 aerosol exposure elicited a weak and mostly transient response (occurring only at the early 4 h post-exposure time point). The results also suggest that some specific gene isoforms that are involved in both the xenobiotic and oxidative responses were induced in a tissue-specific manner following the exposures. Such information could identify a distinct response profile across the buccal, bronchial, and nasal epithelium cultures.
Based on the causal network enrichment analysis of the transcriptomic data, Fig. 4 shows that for the three organotypic epithelium cultures, 3R4F CS exposure elicited a common inflammatory response mechanism (as modeled in the Epithelial Innate Immune Activation network). Furthermore, a correlation analysis was conducted, for which the perturbation scores for the molecular entities in this network (see Fig. 1, panel D) were correlated across the three organotypic cultures. Fig. 7A showed that 4 h post-exposure to 3R4F CS, the perturbation scores of the molecular entities of the Epithelial Innate Immune Activation network, across the three tissues, were highly correlated. Higher Spearman correlations were observed for the comparison between the bronchial and nasal cultures (both of which are pseudostratified epithelium) as compared with the correlations between the bronchial/nasal and the stratified buccal epithelium cultures (despite a fairly high correlation at the 4 h post-exposure time point, Fig. 7A). The lower correlation between bronchial/nasal and buccal may reflect the previously reported differences in the immune response between the buccal (stratified epithelium) and bronchial/nasal (pseudostratified epithelium) cultures. For example, differences in the induction of antimicrobial human beta-defensins (hBDs) have been reported: hBDs are induced by TNFα, IL-1, and phorbol 12-myristate 13-acetate, but not by lipopolysaccharides in the oral epithelium; whereas lipopolysaccharides induced the gene expression of hBDs in the tracheal epithelium.62 Another possible reason for this discrepancy could be attributed to the thicker morphology of the buccal cultures, or to the different donors from which the primary cells were isolated for the generation of the culture models. The stratified epithelium may trigger a delayed kinetic response to exposure (discussed earlier). Furthermore, as shown in Fig. 4, THS2.2 aerosol exposure elicited only small alterations in the global gene expression of the buccal, bronchial, and nasal cultures, except at the 4 h post-exposure time point. The perturbation scores of the molecular entities within the Epithelial Innate Immune Activation network model were much lower as compared with that for 3R4F CS exposure (Fig. 4). As a result, the correlation values between the perturbation scores of the molecular entities within the Epithelial Innate Immune Activation were fairly low across the three cultures (and only well correlated at the 4 h post-exposure time point, Fig. 7A). This analysis demonstrates that the causal-network enrichment approach could further identify the exposure-related molecular mechanisms that are common (and different) across the three epithelium cultures.
Fig. 7 Exposure-induced pro-inflammatory responses across the buccal, bronchial, and nasal cultures. (A) Correlations of the perturbation scores of the molecular entities in the Epithelial Innate Immune Activation network across the three organotypic models. (B) Multianalyte profiling data for secreted pro-inflammatory mediators measured at various post-exposure time points (as a cross-sectional sampling) are shown. The fold-changes relative to the respective sham groups are color coded. Grey cells indicate that no measurement was conducted for the specific mediators. ESI Fig. 3† shows the data for all the tested concentrations. (C) Correlations between the secreted mediators and their corresponding mRNAs (based on the transcriptomic data) for CXCL8, MMP1, TGFA, and TIMP-1 are shown. Correlation plots for the complete set of mediators and concentrations are reported in ESI Fig. 4.† The log2 fold-change (relative to the air control) of a secreted mediator at a given post-exposure time point (i.e. 24 h, 48 h, and 72 h) (y-axis) is compared with the average log2 fold-change (relative to the air control) of the respective mRNA for the same and preceding times to account for accumulation of the secreted protein products (x-axis). |
For the analyzed panel of secreted pro-inflammatory mediators, 3R4F CS exposure resulted in shared but also distinct effects across the three culture models – with a higher similarity between bronchial and nasal than buccal cultures (Fig. 7B, ESI Fig. 3†). A generally consistent increase in matrix metalloproteinase 1 (MMP-1) levels was observed at all post-exposure time points following 3R4F CS. The association between CS and the increased expression of MMP-1 has been reported previously.63 In lung epithelial cells, the link between CS and emphysema has been attributed to the alteration in MMP-1 levels.63 Increased levels of the secreted vascular endothelial growth factor alpha (VEGFA) were also observed following 3R4F CS exposure in all three culture models. In contrast, CXCL8 and IL6 levels were decreased following 3R4F CS in the buccal culture, but tended to be increased in the bronchial and nasal cultures. This observation further supports the distinct inflammatory response between the pseudostratified epithelium (bronchial/nasal) and stratified epithelium (buccal) discussed earlier. The effects of THS2.2 aerosol exposure on these pro-inflammatory mediators generally remained lower and less commonly significant than those upon 3R4F CS exposure at these comparable doses. Finally, an integrative analysis of RNA and protein levels revealed good correlations for CXCL8, MMP-1, TGFA, and TIMP-1 levels (Fig. 7C), suggesting that these mediators are regulated predominantly at the transcriptional level.
Across the three cultures, 3R4F CS exposure was linked to a reduced expression of miR-99a, miR-361, miR-149, miR-125b, and miR-100. Various studies have demonstrated the role of miR-99a in growth factor signaling pathways, including the insulin-like growth factor 1 receptor (IGF-1R) and mammalian/mechanistic target of rapamycin (mTOR).65–68 The miR-99a is thought to regulate cell proliferation and differentiation, although many studies were done using cancer cells.66–70 In contrast, the levels of miR-99a in the cultures exposed to THS2.2 aerosol were not significantly different from the air control.
Increased levels of miR-224, miR-192, and miR-132 following 3R4F CS exposure were consistently found in all epithelium cultures (buccal, bronchial, and nasal). A study has shown that miR-132 regulates cell proliferation by suppressing the retinoblastoma protein (RB1).71 Following allergen exposure, the levels of miR-132 in the bronchial brushing samples were increased, which was linked to shedding of bronchial epithelial cells via the repression of cyclin dependent kinase inhibitor 1A (CDKN1A).72 Moreover, following exposure to ozone, increased miR-132 expression in human sputum cells was reported.73 The levels of miR-132 were not impacted by the THS2.2 aerosol exposure at all post-exposure time points. The levels of miR-4521 were regulated differently following exposure to 3R4F CS: decreased expression was observed in the buccal cultures but increased expression was observed in the bronchial and nasal cultures. A study has reported that the single nucleotide polymorphism rs7210250—located near miR-4521 loci—was associated with risk of esophageal adenocarcinoma.74 Tobacco smoking is one of the risk factors of the disease. The observed difference in the directionality of the miR-4521 expression following 3R4F CS across the cultures could be attributed to the different tissue type or the different donor profile from which the organotypic cultures were reconstituted. Nevertheless, the specific role of miR-4521 has not been reported.
Overall, eight out of the ten miRNAs, which were significantly differently expressed compared with the air controls, were similarly altered following 3R4F CS and THS2.2 aerosol exposure in all three cultures (Fig. 8), demonstrating a common regulatory miRNA-profile across the buccal, bronchial, and nasal cultures. Consistent with the other endpoints analyzed in this study, THS2.2 aerosol exposure in general elicited reduced impact on the alterations in miRNAs compared with that following 3R4F CS exposure at comparable nicotine concentrations.
Fig. 9 shows the expression of various proteins that were altered at 48 h after exposure measured using mass spectrometry-based targeted proteomics. Exposure to 3R4F CS was associated with increased levels of proteins involved in the xenobiotic protein response. First, the cytochrome P450 enzymes CYP1A1 and CYP1B1 had the strongest up-regulation following 3R4F CS exposure. Both CYP1A1 and CYP1B1 are known to be regulated upon activation of the arylhydrocarbon receptor following exposure to xenobiotics.75 Second, the levels of the aldehyde dehydrogenase 3A1 (ALDH3A1) protein were increased in the nasal cultures exposed to 3R4F CS compared with the air controls. ALDH3A1 is frequently found to be up-regulated following CS exposure; for example, increased levels of ALDH3A1 were detected in the fluid lining the epithelial cells in smokers.76 Moreover, the abundance of ALDH3A1 in the sputum could distinguish smokers from non-smokers, and thus can serve as a specific marker of the smoking status.77 These observations demonstrate consistently the greater levels of ALDH3A1 after CS exposure both in vivo and in vitro. Finally, the increased expression of AKR1B10, an aldo-keto-reductase, following 3R4F CS exemplifies the culture response against toxic aldehydes and CS exposure, as reported previously.78,79 Overall, the alteration in the protein markers further confirmed that 3R4F CS exposure impacted the xenobiotic metabolism response of the nasal cultures (based on the network enrichment analysis of the transcriptomic data).
Furthermore, the proteomics data also support the 3R4F CS-induced oxidative stress discussed earlier. The protein levels of both the catalytic subunit of the glutamate cysteine ligase (GCLC) and the NAD(P)H-quinone dehydrogenase 1 (NQO1) were significantly increased in the nasal cultures exposed to 3R4F CS compared with the air control (Fig. 7). GCLC is the catalytic subunit of the rate limiting enzyme for glutathione synthesis; thus GCLC is essential for the oxidative stress reaction (which ensues following smoke exposure).58 NQO1 is known as an oxidative stress enzyme, and acts as a quinone reductase. Increased levels of NQO1 following CS exposure in the sputum of smokers has been reported previously.77 Similarly, increased NQO1 levels were observed in the large airway epithelium of smokers.80 In addition, the generation of NADPH is essential for the cellular response against oxidative stress.59,81,82 In this study, the two key enzymes regulating the generation of NADPH—glucose-6-phosphate dehydrogenase (G6PD) and malate enzyme (ME1)—were up-regulated following 3R4F CS exposure. Overall, the proteomics data further confirm that 3R4F CS induced oxidative stress response in the nasal cultures.
These proteomics data could further confirm the findings obtained using the causal network enrichment approach. For example, based on the transcriptomic data, 3R4F CS exposure was linked to a significant perturbation of the Senescence network model (Fig. 4C). Here, the proteomics analysis identified that the expression of p21—a cell cycle regulator known to accumulate in senescent cells83,84—was increased following 3R4F CS (Fig. 9). In addition, the increased expression of stratifin (SFN) and the small proline rich protein 1A (SPRR1A) following 3R4F CS could be linked to the significant perturbation of the Tissue Damage network (Fig. 4C). The expression of SFN is mainly associated with squamous epithelial cells,85 thus suggesting that 3R4F CS exposure might trigger a cellular differentiation in the nasal cultures after exposure. Moreover, SPRR1A plays a role in the barrier function of the epithelium.86 Another isoform of the protein, SPRR3, was also reported to be increased in the sputum and large airway epithelium of smokers.77,80
Consistent with the other data modalities (functional measures, secreted mediator profiles, and transcriptomics), when compared to the impact of 3R4F CS, THS2.2 aerosol exposure resulted in much more limited alterations of protein expression. Only when the nasal cultures were exposed to THS2.2 aerosol at a nicotine concentration three times that of 3R4F CS, increased levels of a subset of the proteins were observed (CYP1B1, AKR1B10, NQO1, ME1, G6PD, and CGN) (Fig. 9). Nevertheless, the degree to which these proteins were altered remained lower than that observed following 3R4F CS exposure. Overall, the results illustrate that an additional data modality—in this present study, the proteomics data—can further support and confirm the observations based on the analysis of transcriptomic data i.e., the reduced impact of THS2.2 aerosol as compared with 3R4F CS was linked to the reduced xenobiotic metabolism and oxidative stress responses.
The present meta-analysis includes functional measurements (cytotoxicity, ciliary beating functionality, and pro-inflammatory mediator profiles) and advanced computational approaches leveraging gene set analyses and causal network enrichment to comprehensively assess the biological impact of 3R4F CS and THS2.2 aerosol exposures on the in vitro human organotypic buccal, bronchial, and nasal cultures. An exposure characterization was included in the assessment to ensure the reproducibility and the consistency of the smoke/aerosol generation throughout the studies. On the mechanistic level, the responses of the three organotypic cultures to 3R4F CS were well aligned, with a prominent engagement of xenobiotic metabolism, oxidative stress, and pro-inflammatory mechanisms. Buccal tissues demonstrated an overall lower sensitivity to the cytotoxic effects of 3R4F CS exposure, which could be attributed to the thicker layer of the buccal culture (a stratified epithelium) and the different donors from which the primary cells were isolated (Fig. 2). The culture responses—based on the transcriptome changes—could further demonstrate that THS2.2 aerosol exposure elicited a more transient response (at the 4 h post-exposure time point) in all the three cultures. The analysis further demonstrated some differences, for example, that specific isoforms of genes or proteins were regulated differently across the three cultures following exposure. The analysis based on the various data modalities further confirmed the overall reduced and more transient biological impact of THS2.2 aerosol exposure, as compared with 3R4F CS. The similarity of study design across the three studies (buccal, bronchial, and nasal) enabled a robust and reproducible analysis, which supported an overall reduced impact of THS2.2 aerosol exposure on the “field of injury” tissues in vitro, compared with that of CS.
In the current studies, these in vitro models were used for the comparative assessment of smoke/aerosol at similar nicotine concentrations, to investigate relative differences in the effects elicited in biologically relevant culture systems. Currently, the in vitro exposure scenario does not directly mimic the exposure situation in the aerodigestive tract – e.g., because of differences in the puffing parameters, the thermodynamic state of the flowing smoke/aerosol, and the geometry of the flow path in vitro and in vivo. Therefore, we have ongoing efforts in our group90,91 that incorporate computational fluid dynamics models to further refine our knowledge and understanding of smoke/aerosol behavior in the human respiratory tract and in vitro systems – to eventually allow for more direct translatability between the culture systems and the in vivo situation.
Overall, the present paper demonstrates that the systems toxicology approach is a robust methodology that can not only detect toxic effects of exposure, but also infer the mechanistic-based toxicity assessment. Moreover, this approach could deduce relevant biological mechanisms, such as cell stress and inflammatory processes—which are relevant for the pathophysiological effects of smoking in humans—thus offering an opportunity for translation to tissue-specific clinical biomarkers.
This work also exemplifies how targeted proteomics can strengthen the findings observed using the causal network enrichment analysis based on transcriptomic data: the xenobiotic metabolism and oxidative stress responses were prominently elicited by 3R4F CS and only slightly impacted by THS2.2 aerosol exposure. In the future, the systems toxicology approach will mature as additional omics data modalities will be incorporated. Such multimodality data will enable investigations into the exposure effects at multiple levels of biological organization. In previous work, our group and others have demonstrated that multi-omics profiling can unravel toxicological mechanisms: an integration of transcriptomics, proteomics, and metabolomics resulted in a comprehensive characterization of the effects of Cyclosporine A and cisplatin on cell stress in renal epithelial cells;92,93 a combination of transcriptomics and metabolomics revealed that reduced oxidative stress was linked to THS2.2 aerosol exposure as compared with exposure to CS in human organotypic gingival cultures;54 and integration of transcriptomics, proteomics, and lipidomics data revealed the impact of candidate MRTP aerosols and CS on lipid metabolism in mouse lungs.82 In the future, it will be pertinent to directly integrate proteomics results with the causal network approach, leading to the generation of multi-scale network models.
The scientific community will play an important role in the success of systems toxicology as an assessment approach. Indeed, making the network models, algorithms, and datasets available to the scientific community is critical to validate the applied methodologies. Motivated by the Tox2194 and the OpenTox projects,95 which incorporate public repositories of large-scale toxicity data, the Causal Biological Network Database (CBN) has been made available to share the Causal Biological Network Models (http://www.causalbionet.com/).29 The Systems Biology Verification project (sbvIMPROVER) (http://sbvimprover.com) was organized to evaluate microarray-based phenotype predictions,96 species translation,97 generation of a comprehensive set of COPD-relevant models,36 and recently blood-based gene expression signatures to classify exposure status between 3R4F CS and MRTP aerosols.98 Furthermore, the INTERVALS platform has been developed to share results from in vivo inhalation studies and in vitro studies in the context of product assessment.99
In the future, the systems toxicology approach—guided by the adverse outcome pathway (AOP) framework100—is envisioned to mature and eventually advance into dynamic AOPs that can be used to quantify and/or predict all the steps of a toxicological response. Such a response is ranging from the initial molecular interactions between the toxicant and the host system (molecular initiating event), to the cellular /organ events, and finally to the organism- and population-level events4 (Fig. 10). An integrative systems toxicology assessment framework can be broadly applied to all relevant areas of toxicology, including food safety, pharmacological drug safety, and occupational safety. Although the present work reports the application of the systems toxicology approach in the context of tobacco product assessment, the same approach has already been shown to be applicable for other toxicity evaluations, such as for nutraceuticals.101 A literature curation platform for causal network models has been established that enables a straight forward expansion to other disease areas, as was demonstrated for atherosclerosis plaque destabilization.102 Indeed, this framework is envisioned to evolve beyond traditional toxicological assessment approaches. For example, incorporating genetics and other personalized risk factors will allow an individual-based safety assessment, similar to personalized medicine.103 Beyond toxicological assessment, this approach also has the potential to generate novel insights into disease mechanisms and support the identification of pharmaceutical interventions, as already proposed for the AOP framework.100 Thus, continuing further on this path, the outlined strategy and approaches are well equipped to make the vision of 21st century toxicology a reality.5,87
Fig. 10 Dimensions of the systems toxicology paradigm. Toxicological effects percolate through biological systems, from the initiating molecular interactions to cellular and tissue responses to population effects (upper panels). Different measurement techniques enable these effects to be quantitatively pursued at all levels. As systems toxicology evolves, additional levels of the exposure effect can be captured in a single mathematical model (lower panel). In this study, the causal network models captured the cellular responses based on transcriptomic measurements, which were complemented by additional functional measurements, e.g. cytotoxicity and cilia beating analysis. The development of dynamic multi-scale models will allow a uniform assessment of all levels of the exposure effects from molecular interactions to tissue/organ responses. Eventually, dynamic adverse outcome pathway (AOP) models will bridge the entire toxicological response, from molecular initiating events to population effects. Adapted from the artwork of Samantha J. Elmhurst (http://www.livingart.org.uk) published previously.4 |
The meta-analysis reported here demonstrates the applicability of the systems toxicology approach to generate comprehensive and consistent data showing the culture response to 3R4F CS and THS2.2 aerosol on several relevant biological mechanisms, including cellular stress and pro-inflammatory responses. The results show consistently across all three in vitro models—buccal, bronchial, and nasal—that THS2.2 aerosol exposure had a considerably reduced and more transient biological impact on these in vitro models compared with equivalent exposures to 3R4F CS.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c7tx00047b |
‡ These authors contributed equally to this work. |
This journal is © The Royal Society of Chemistry 2017 |