Next Article in Journal
Tailoring InSb Nanowires for High Thermoelectric Performance Using AAO Template-Assisted Die Casting Process
Previous Article in Journal
Z-Scheme Heterojunction of SnS2/Bi2WO6 for Photoreduction of CO2 to 100% Alcohol Products by Promoting the Separation of Photogenerated Charges
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Nanomaterials
Volume 12
Issue 12
10.3390/nano12122031
nanomaterials-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Characterization of ENM Dynamic Dose-Dependent MOA in Lung with Respect to Immune Cells Infiltration

by
Angela Serra
1,2,3,†,
Giusy del Giudice
1,2,3,†,
Pia Anneli Sofia Kinaret
4,
Laura Aliisa Saarimäki
1,2,3,
Sarah Søs Poulsen
5,
Vittorio Fortino
6,
Sabina Halappanavar
7,
Ulla Vogel
5 and
Dario Greco
1,2,3,4,*
1
Faculty of Medicine and Health Technology, Tampere University, 33520 Tampere, Finland
2
BioMediTech Institute, Tampere University, 33520 Tampere, Finland
3
Finnish Hub for Development and Validation of Integrated Approaches (FHAIVE), 33520 Tampere, Finland
4
Institute of Biotechnology, University of Helsinki, 00014 Helsinki, Finland
5
National Research Centre for the Working Environment, 2100 Copenhagen, Denmark
6
Institute of Biomedicine, University of Eastern Finland, 70211 Kuopio, Finland
7
Environmental Health Science and Research Bureau, Health Canada, Ottawa, ON K1A 0K9, Canada
*
Author to whom correspondence should be addressed.
†
These authors contribute equally to this work.
Nanomaterials 2022, 12(12), 2031; https://doi.org/10.3390/nano12122031
Submission received: 19 May 2022 / Revised: 8 June 2022 / Accepted: 9 June 2022 / Published: 13 June 2022
Browse Figures
Versions Notes

Abstract

:
The molecular effects of exposures to engineered nanomaterials (ENMs) are still largely unknown. In classical inhalation toxicology, cell composition of bronchoalveolar lavage (BAL) is a toxicity indicator at the lung tissue level that can aid in interpreting pulmonary histological changes. Toxicogenomic approaches help characterize the mechanism of action (MOA) of ENMs by investigating the differentially expressed genes (DEG). However, dissecting which molecular mechanisms and events are directly induced by the exposure is not straightforward. It is now generally accepted that direct effects follow a monotonic dose-dependent pattern. Here, we applied an integrated modeling approach to study the MOA of four ENMs by retrieving the DEGs that also show a dynamic dose-dependent profile (dddtMOA). We further combined the information of the dddtMOA with the dose dependency of four immune cell populations derived from BAL counts. The dddtMOA analysis highlighted the specific adaptation pattern to each ENM. Furthermore, it revealed the distinct effect of the ENM physicochemical properties on the induced immune response. Finally, we report three genes dose-dependent in all the exposures and correlated with immune deregulation in the lung. The characterization of dddtMOA for ENM exposures, both for apical endpoints and molecular responses, can further promote toxicogenomic approaches in a regulatory context.

Graphical Abstract

1. Introduction

When ENMs interact with biological systems, they induce a cascade of molecular events which are mainly dependent on their physicochemical characteristics, such as the size and the surface chemistry [1,2,3].
Cellular responses to ENM exposures are complex, as they are trying to counteract the induced perturbation and achieve system homeostasis. In classical inhalation toxicology, bronchoalveolar lavage (BAL) cell compositions are widely used to clarify the immunopathological phenotypes in the lung [4]. ENMs exert diverse responses both in terms of cell types recruited in the lung and in their numbers [5,6,7,8,9,10]. However, apical endpoints, such as the one derived by BAL counts, do not provide details about the underlying mechanisms and responses to distinct ENMs. Usually, each significant molecular alteration triggered by ENM exposure is considered as the MOA of that compound [11].
Microarray experiments are an effective tool for the characterization of the mechanism of action detected at transcriptome levels (transcriptional mechanism of action, hereafter referred to as tMOA) [12,13,14]. However, the dissection of this response in the molecular events directly triggered by the exposure to the ENM and the ones induced by the biological adaptation process is complicated.
In risk assessment and toxicology, dose-dependent modeling is an established tool to identify points of departure (POD) that can be used as a guidance value denoting the adverse effect levels, such as benchmark dose (BMD), no observed adverse effect level (NOAEL), and lowest observed adverse effect level (LOAEL) [15]. Similar modeling strategies have emerged for the analysis of transcriptomic experiments to identify the genes directly and progressively affected with respect to the reference dose or concentration [16,17,18,19] and to characterize their POD.
The tMOA of an exposure arises from different mechanisms of transcriptional regulation. Therefore, some of these events are observed as dose-dependent, whereas others are consequences of more complex regulatory loops.
The common strategy for the investigation of dose responses is to study their alteration at each time point of exposure separately. However, this strategy might not accurately interpret the kinetics of the molecular adaptation processes [20,21]. The joint dose-time dynamic POD (dPOD) allows us to identify the lowest dose and earliest time points at which a notable change with respect to the controls takes place, and this effect follows a monotonic dose-dependent trend. Moreover, few efforts have been made to link the traditional immune cell quantification and transcriptional responses together [22,23,24,25]. To the best of our knowledge, the correlation between specific portions of dddtMOA and different immune cells has never been systematically investigated.
To overcome these problems, we applied a computational approach [21] to model the dose- and time-dependent alteration patterns of publicly available BAL cell counts and lung transcriptomic responses from multiple ENMs, namely, carbon black (CB), titanium dioxide (TiO2), and two multiwalled carbon nanotubes (MWCNT), at three different time points. Our goal was to characterize the dddtMOA as the portion of the adaptation response directly linked to the ENM, and to further associate its effect with the immune cell infiltrations.
Our approach suggests that the dddtMOA has the potential to be applied in the context of an ENM hazard assessment, further facilitating the implementation of toxicogenomic-based approaches in a regulatory context.

2. Materials and Methods

2.1. Dataset Collection and Preparation

Microarray-based transcriptomics data from lungs of mice exposed, by a single intratracheal instillation, to carbon black (CB), multiwalled carbon nanotubes (MWCNTs) and titanium dioxide nanoparticles were downloaded from the GEO database. Details on the nanoparticles are reported in Table 1.
Data for MWCNTs came from the studies of Poulsen et al. [26] (available on GEO with ID GSE55286), where female C57BL/6 mice were exposed by a single intratracheal instillation to 18, 54, or 162 g/mouse of a short MWCNT (NRCWE-26, 847 ± 102 nm in length) or long MWCNT (NM-401, 4048 ± 336 nm in length). Data for CB nanomaterials came from the studies of Bourdon et al. [27,28] (available on GEO with ID GSE35193), where female C57BL/6 mice were exposed by a single intratracheal instillation to 18, 54, or 162 g/mouse of Printex 90 carbon-black nanoparticles. Data for nano-TiO2 nanomaterials came from the studies of Husain et al. [29] and Saber et al. [30] (available on GEO with ID GSE41041), where female C57BL/6 mice were exposed to rutile nano-TiO2 (primary size of 20.6 nm and surface area of 107.7 m2/g). Experiments related to CB and TiO2 were performed by using the Agilent-014868 Whole Mouse Genome Microarray 4x44K G4122F platform, whereas for the MWCNT data, the Agilent-028005 SurePrint G3 Mouse GE 8x60K Microarray platform was used. Raw mRNA expression values were imported in R v. 3.4.4 using the Limma read.maimages function [31]. Raw data were preprocessed using the eUTOPIA application [32]. First, low-quality probes were removed: only probes with a value higher than 75% quantile of negative control probes in at least 75% of the samples were selected. Afterwards, the expression values were log2 transformed and quantile normalized between arrays. Surrogate technical variables were identified and removed in the CB and TiO2 datasets by using the ComBat method from the R Ensembl gene ID level using annotation files provided by the microarray manufacturer (Agilent Technologies, Santa Clara, CA, USA).
Finally, differential expression between sample groups was evaluated with the Limma package using the corrected batches as covariates for the model. The combination of the three exposure doses (18, 54, and 162 µg) and three time points (1, 3, and 28 days) yields a total of nine comparisons when each dose- and time-point combination is compared to its corresponding control sample. Genes were considered significantly differentially expressed with an absolute log2 fold change > 0.58 and Benjamini–Hochberg-adjusted p-value < 0.05. For the sake of comparison with respect to the dynamic dose-dependent genes, we computed the union of all the differentially expressed genes revealed by the nine comparisons and used that to characterize each ENM.
Moreover, the BAL counts for neutrophils, eosinophils, macrophages, and lymphocytes, present in the lungs after exposure to the same ENMs considered in this study, were retrieved from previous publications [26,28,30] (Table S1).

2.2. Integrated Time and Dose Analysis of MOA

Sample-wise log2 fold changes were calculated between each exposed sample and each of its corresponding control samples. Afterwards, each gene was analyzed utilizing the TinderMIX tool [21]. Briefly, linear- and second-order polynomial models were applied for the fitting of the log2 fold change. For each gene, the optimal model was selected based on the lowest goodness-of-fit p-value. Genes with goodness-of-fit p-values > 0.05 were removed from the analysis. The dose and time ranges were divided into 50 equally distributed bins, and the optimal model of each gene was used to predict their corresponding log2 fold changes. Thus, the gene was represented by an activation map that interpolates the space of time and dose also for the experimental conditions (doses and time) not originally included in the analysis.
From each activation map, a monotonically increasing or decreasing area with respect to the dose, with a fold change greater than the activity threshold (fold change > |1.5|), was identified by analyzing the gradient vector field of the map. If present, this area is marked as responsive, and the gene is deemed altered in a dynamic dose-dependent manner.
Eventually, an activity label was assigned to each gene by dividing both axes of the activation map into three sections and obtaining a grid with 9 cells. The sections of the dose axis were named as “sensitive”, “intermediate”, and “resilient”, whereas for time axis, labels “early”, “middle”, and “late” were assigned. The final dPOD is obtained by identifying the earliest and most sensitive activation and is labeled concatenating the corresponding dose and time labels.

2.3. Integrated Time and Dose Analysis of Cell Counts

For each pair of dose and time points, the cell counts were normalized for the counts in the control sample in the following way:
C C t p x   d s = log 2 c o u n t s t r /   c o u n t s c t r l
where tp is the time point, ds is the dose, tr is the treated samples, and ctrl are the controls. Within the samples and controls of a particular pair of dose and time points, all the possible pairwise comparisons are performed.
Afterwards, the TinderMIX tool was used to identify the dynamic dose dependency of the BAL cell counts. Starting from the normalized cell counts, the linear- and second-order polynomial models were fitted in the joint dose-time space. The optimal model was used to predict an activation map associated with each cell count in the form of contour plots. The optimal model was identified as the one with the lowest goodness-of-fit p-value, among those with p-values < 0.05. If none of the models met the requirement for the p-value, the optimal model was selected as the one with the smallest AIC (Akaike information criterion). Eventually, a dPOD activity label was associated with each cell count in each material.

2.4. Correlation between Gene Expression and Cell Counts

For each material, the correlation between the maps of its dynamic dose-dependent genes and the maps of the cell counts were computed. For each pair of gene and cell types, we divided the map in three different regions based on if both gene and cell were dose-dependent (called region A), only one of the two was dose-dependent (called region B), or none of the two were dose-dependent (called region C). For each portion of the map, we computed a binary vector indicating if the activation was higher (1) or lower (0) than the mean activation value of that portion. The Pearson correlations between the binary vectors were computed. This led to three correlation values for each pair of gene and cell, which were combined by a weighted sum that accounted for 70% of the correlation from region A, 20% of correlation from region B, and 10% of correlation from region C. Since the aim of this study was to identify dynamic-dose dependent alterations induced by the exposure to ENMs, the weighted schema was selected to give more relevance to the area where both the gene and the cell counts were dynamic-dose dependent. However, rather than discarding correlations outside this area, they were considered still relevant to explain interactions between the genes and cell populations. To avoid spurious correlations, gene activity and cell counting maps were considered correlated when the final absolute Pearson correlation value was greater than 0.6.

2.5. Pathway Analyses

The KEGG pathway enrichment analysis was performed with the compareCluster function from the clusterProfiler R package [33]. To compute the combined p-values for the pathways shared across the ENMs, we used the sum of logs method (or Fisher’s method) implemented in the sumlog function of the metap R package [34].

3. Results and Discussion

3.1. Characterization of the dddtMOA

The tMOA of an exposure is commonly defined as the set of significant variations in gene expression observed between the exposed and unexposed samples. Therefore, we first performed a classical differential expression analysis of the dataset including exposure to TiO2, CB, and two MWCNTs of different length with multiple doses (18, 54, and 162 mg) and multiple time points (1, 3, and 28 days).
Particularly, we identified lists of differentially expressed genes for each combination of material, dose, and time. As for the lung inflammation, we retrieved information on BAL counts of macrophages, neutrophils, eosinophils, and lymphocytes in the same experimental conditions [26,28,30] (Table S1).
To characterize the dynamic dose-dependent portion of the tMOA (dddtMOA) of each nanomaterial, we performed an integrated dose and time analysis of the molecular alteration induced on the whole transcriptome and we investigated the distribution of the number of dynamic dose-dependent genes across their dose and time dPOD (Figure 1A, Tables S2–S5).
Both MWCNTs showed the highest number of dose-dependent genes (2484 and 2670 for NM401 and NRCWE26, respectively), whereas the number of genes altered by CB and TiO2 was lower (232 and 298, respectively). Our results suggest that the magnitude of the response could reflect the ENM hazard level. For example, a high aspect ratio and rigid nanomaterials, such as NM401, are usually more persistent in the lung, and their reduced clearance contributes to the toxic profile [35].
In our dose- and time-dependent analysis, the genes were labeled based on where their expression deregulation starts with respect to the time (early, middle, or late) and the dose levels (sensitive, intermediate, or resilient).
As expected, the kinetics of the molecular response in all ENMs immediately start at short time points, having an early dPOD. This might be partially influenced by the single-exposure experimental design, where stronger acute responses can be detected at shorter time intervals, whereas mechanisms associated with recovery can be detected at later time points. However, our results are consistent with the acute adaptive response to ENM exposure [36]. As for the exposure dose, MWCNTs seem to activate most genes at low doses, having a sensitive dPOD (Figure 1A). As discussed before, the dose required to induce the alteration can be linked to the toxicity level of the ENMs. Interestingly, the most represented label in the CB and TiO2 exposures is early resilient. This may indicate that for spherical-shaped nanomaterials, a portion of the response requires higher doses to be evident.
To functionally characterize the dddtMOA, we first compared the differentially expressed genes (DEG) and the genes with a dPOD activity label of the four ENMs, and we found an intersection in all the materials (Figure 1B). Particularly, 185 genes were found for the TiO2, 216 for the CB, and 2347 and 2638 for the NM401 and NRCWE26, respectively. This supports our assumption that the MOA of the different ENMs is only partially dose-dependent, and that monotonic portions of the molecular adaptation response may be directly altered by exposure.
We then investigated the functional profile of the dPOD genes by mapping them onto the KEGG pathways (Tables S6–S9). The functionalities associated with each dddtMOA were ENM specific; however, some functions were enriched by all the studied nanomaterials, as well as by specific subcategories (Figure 2). Among the specific ENM responses, our analysis highlighted pathways associated with cardiac muscle contraction and signaling in cardiomyocytes, suggesting a possible cardiotoxic effect of titanium dioxide (Figure 2D). It was previously shown that exposure to TiO2 nanoparticles both in vitro and in vivo increases cardiac excitability and conduction velocity [37,38]. Similarly, the most significant alterations in NM401 exposure pointed to lysosomal effects (Figure 2G). Straight and rigid carbon nanotubes, such as NM401, are known to induce disruption of cellular structures and vesicles, including the lysosomes [6].
As for shared responses, all spherical-shaped ENMs enriched pathways related to microorganisms and viral infections (Figure 2A). Similarities in the host response between carbon-based ENMs and viruses were already observed [39,40]. Interestingly, pattern-recognition receptors are also commonly present. Toll-like receptors are the sentinels of foreign biomolecules, and they mediate unspecific recognition of damage- and pathogen-associated molecular patterns (DAMP and PAMP, respectively). This may suggest that the exposed cells sense spherical nanoparticles through similar signaling pathways that evolved to recognize pathogens, reinforcing the importance of physicochemical properties of the nanomaterials in shaping the adaptive response. Interestingly, toll-like receptors are known to induce necroptosis, an inflammatory form of programmed cell death, which was enriched in the carbon-black dddtMOA. MWCNTs specifically enrich pathways related to the extracellular matrix and focal adhesion, as well as vascular structures (Figure 2C). MWCNTs are known profibrotic agents, and their hazardous effects on the lung epithelium were already described in detail [40,41,42,43]. Interestingly, carbon-based materials (CNMs) share a direct effect on the complement and coagulation cascade (Figure 2B). CNMs are indeed used in coagulation disorder therapy because of their potential to interact with platelets and induce thrombosis and hemolysis [44].
Finally, all the investigated ENMs shared a core of functions related to inflammation and immunity (Figure 2H). Nanomaterials trigger the immune response and induce activation of inflammatory mediators through various mechanisms. Importantly, NFkB and TNF are very well characterized mediators of nanomaterial adaptation [45,46,47,48]. Similarly, upregulation of a plethora of cytokines was reported for all the studied nanomaterials [45,46].

3.2. Combination of Transcriptional Changes and BAL Cell Count Dose-Time Responses Informs on ENM-Specific Immune Cell Activation

The interaction of ENMs with the immune system is an accepted theory; however, clarifying the cellular mechanism of ENM-mediated toxicity is still considered a major challenge.
We hypothesized that combining the dose-dependent and time-related transcriptional changes with the BAL cell count activation may help to better characterize the molecular basis of the immune cell activation and response.
First, we performed a similar dose- and time-integrated analysis on the BAL cell counts and identified their joint dynamic dose-dependent dPOD (Table S10).
In the case of neutrophils, the dPOD for all the ENMs was determined as sensitive-early, specifying a dynamic dPOD starting at a low dose range and an early time interval and preserved through the postexposure time points. This suggests that the immediate cell infiltration consequent to the insult extends until later time points. Neutrophils have long been known as “the first line of defense”, as well as the initiator of the inflammatory response [49]. In this light, it is expected that neutrophils, as quickly reacting cells, are labeled as early responders.
Macrophages showed different dPOD labels in different ENMs. They were sensitive-early in all the CNMs and intermediate-late in TiO2. Macrophages are equally considered as rapid responders, but unlike neutrophils, which migrate to the site of insult and survive only 1 to 2 days, macrophages are already residing in lung tissue and survive much longer at the site of inflammation. More circulating macrophages might be recruited at later time points or higher doses, thus possibly explaining their varying activation patterns across the different ENMs.
Eosinophils showed sensitive-late dPOD in CB, still suggesting cell infiltration after the exposure but at the longest time points [50].
Lymphocytes also showed a dose-time response only in spherical-shaped ENM. In the case of CB they were labeled as sensitive-early, whereas in TiO2 they were activated later at low doses.
To contextualize the dddtMOA with respect to the activation of immune cells, we identified the sets of genes with a dynamic dose-dependent activation which were also correlated with the immune cell counts (Figure 3).
Macrophages and neutrophils showed the highest number of correlated dPOD genes (Figure 3A). This can be easily explained by their key role in acute responses to exogenous compounds, as well as their extremely plastic activation upon cues received from their immediate microenvironment [51,52].
However, when we analyzed the relative composition, a difference between nanotubes and spherical-shaped nanomaterials arose (Figure 3B). MWCNTs are predominantly correlated with macrophages (about 75% for NM401) and neutrophils (about 60% for NRCWE26). On the contrary, TiO2 mainly induced the response of lymphocytes and eosinophils, and CB of lymphocytes exclusively. Several studies have already shown that ENMs can interact in diverse ways with the immune system [36]. The difference seems to be driven by physicochemical characteristics, such as the shape and the size, or the hydrophobicity level. For example, Liu et al. [53] proved that the diameter in silica–titania nanoparticles directly affects the amplitude of the inflammatory reaction.
In conclusion, correlating the dddtMOA transcriptional profile with dPOD labels derived by cell BAL counts can inform on different cellular mechanisms induced by specific ENMs.

3.2.1. ENM Physicochemical Properties Induce Differences in Immune Activation

To characterize the dynamics of the immune response for each ENM, their individual functional profiles were investigated by mapping the genes onto the KEGG pathways (Figure 4, Tables S11–S14). Two clear patterns emerged from genes correlated with neutrophils and macrophages.
Our set of macrophage-correlated genes (Figure 4A) in TiO2 nanoparticles include several mediators, such as CCL7, CCL8, CCL9, and IL4. IL4 is a powerful inducer of M2 polarization and marker of M2 macrophages. M2 macrophages are known to promote tissue repair and an anti-inflammatory phenotype [54]. Comparable results were previously observed for TiO2 nanotubes, which exert a protective effect on epithelial integrity by inducing M2 macrophages [55,56]. CNMs, instead, even if with different amplitudes, share a common functional core of the macrophage adaptive response, including the expression of IL1, NfkBia, TNF, IL10, and TGF. Both TNF and IL1 belong to the panel of inducers of M1 polarization, whereas IL10 and TGF belong to the M2 panel instead. We previously described that CNMs trigger macrophages towards a hybrid M1/M2 phenotype in a nonacute exposure setup [46]. This evidence points to a different mechanism of polarization of macrophages across ENMs that is possibly dependent on ENM chemistry (Figure 5A).
As for neutrophil-correlated genes (Figure 4B), they were characterized by a different adaptive response in TiO2 and CB as compared to NM401 and NRCWE26.
Genes correlated with neutrophils in TiO2 and CB were monocyte-lineage inducers and chemoattractants related to chemokines and cytokine receptor–ligand activity (Figure 4B). It is well recognized that the successful initiation and resolution of inflammation is mediated by neutrophils. GM-CSF genes, such as Csf2 and IL8, directly promote proliferation and maturation of neutrophils, but also affect macrophages and eosinophils via STAT3 and STAT5, facilitating defense against infections. In fact, most of these genes correlated with neutrophils overlap between the two spherical-shaped ENMs, highlighting that sustainment and homeostasis of neutrophils is a conserved mechanism between them. However, a set of genes was specific for CB, including IL6, Clec4n, CCL12, Ctsk, IL33, Pigr, Vinn1, and Ifit1. IL6 is a crucial mediator since it induces neutrophilia under inflammation. Furthermore, when activated via STAT3, IL6 signaling represents an important event for the termination of the immune response [57,58]. Overall, these results point to the central role of the JAK-STAT pathway in sustaining neutrophil proliferation and production of inflammatory mediators (Figure 5B).
On the contrary, neutrophil-correlated genes in MWCNTs suggest induced neutrophilia through a different mechanism (Figure 5B). NRCWE26 has only a few genes correlated with neutrophils, among which are genes associated with cell proliferation (such as Pol2 and Mcm2/3), and, interestingly, various hemoglobin chains and subunits. Hemoglobin was recently described to have a pro-inflammatory and pro-oxidant effect, and it has become clear that it can trigger neutrophil infiltration by recognition of DAMP or GPCR [59,60]. Multiple studies report that the main effect of neutrophils when triggered through GPCR is sustained proliferation and chemotaxis [61]. NM401 neutrophil-correlated genes enrich the epidermal growth factor (EGF) pathway and multiple pathways related to glutathione transferases. Uddin et al. demonstrated that EGF production by epithelial cells, a well-known factor of epithelial repair, increases neutrophil accumulation, which can promote intervention for acute injury, as well as potentially enhance chronic inflammation in the airways [62].
Overall, these results suggest that the mechanism of neutrophil activation is different across ENMs (Figure 5B). This could be due to the activation of a different panel of membrane receptors, and hence, we hypothesize the role of ENM shape as a critical property in activating a more epithelial or strictly immunological response. Supporting our hypothesis, Danielsen et al. already observed that pulmonary exposure to differently shaped CNMs induced differences in the neutrophil influx in vivo [63]. In addition, CB-induced neutrophil influx and transcription of pro-inflammatory genes in lung tissue were shown to be unaffected by TLR2 and TLR4 status, whereas CNT-induced neutrophil influx was somewhat reduced in TLR4−/− mice as compared to wildtype mice [63].
Finally, genes correlated with lymphocytes and eosinophils (Figure 4C,D) enriched only a few pathways and did not show specific trends associated with ENM properties. A possible explanation could be the small number of genes correlated with eosinophils and lymphocytes.
Of note, even if the role of eosinophils in TiO2 and NRCWE26 has not been investigated in detail in the community, our results highlight genes enriching the IL-17 pathway. Th17-derived cytokines, such as IL-17A and IL-17F, exert distinct effects in airway inflammation, and it was postulated that aberrant IL-17A/F production may drive severe forms of asthma [64].

3.2.2. Dynamic Dose-Dependent Genes Shared by All the ENMs

Increasing evidence suggests that dose-dependent alterations highlight the most biologically significant events [65]. In this light, benchmark dose analysis can explain the MOA of ENM exposures by highlighting relevant responders already at short time points.
The approach applied here revealed 39 genes with dose-dependent alteration in the four ENMs that show, however, a specific dynamic pattern for each ENM (Figure 6).
The genes were clustered into seven groups based on their dynamic dose-dependent dPOD genes. Some clusters show specific labels for individual or grouped ENMs. Interestingly, cluster 5 shows an early alteration in the spherical-shaped nanomaterials, but not in the carbon nanotubes. The Sftpd gene in this group was previously proposed as a biomarker of TiO2 exposure [66]. Cluster 7, instead, shows an early alteration specific to CNMs. The Lrg1 gene in this group is an acute phase protein mainly expressed in neutrophils, involved in wound healing and fibrosis. Notably, it was proposed as a reliable marker of airway inflammation in asthma [67].
The biggest cluster includes genes with early activation across the four nanomaterials. Since the aim of our work was to identify dose-dependent responder genes that also correlate with the phenotype observed in the lungs, we focused our analysis on those that are also correlated to at least one cell population in each ENM. Of these, three genes met our criteria, namely, CCL12, CCL7, and IL1b (Table 2). The three genes mostly show a sensitive-early dPOD, meaning that they are already activated at low doses and early time points. However, CCL7 and CCL12 show a resilient response in TiO2 and CB. Moreover, the genes show a material-specific correlation with the different cell populations (Table 2).
CCL7 is an important chemoattractant of neutrophils, macrophages, and natural killer cells, whose expression allows prompt control and the initiation of inflammation [68]. Importantly, CCL7 is a key player in the resolution of inflammation in lungs, suggesting that it may have a central role in the recovery phase of the exposure [69]. CCL12 binds to the same receptor of CCL7, and these axes seem to play a key role in insult response in the lung. In detail, mice with lung epithelial cell-specific deletion of CCL12 were protected from bleomycin-induced fibrosis and had expression of CCL2 and CCL7 similar to that of control mice treated with bleomycin. Many MWCNTs seem to induce bleomycin-like fibrosis [40]. The protective effect of low levels of CCL12 against fibrosis makes this gene an important factor and possibly a therapeutic target. IL1b is an important marker of inflammation, already established in cancer as a prognostic and therapeutic target [70,71]. Importantly, IL1 stimulates CCL7 expression through the NfkB and MAPK pathways in astrocytes [72]. Even if this relationship has not been proven in the lung, our results suggest an important cytokine–chemokine axis, whose deregulation is potentially predictive of lung toxicity.

4. Conclusions

Even though ENMs are widely spread, clarifying the molecular events and cellular mechanisms of adaptation upon exposure is still considered a challenge [12]. BAL cell counts are commonly used to identify inflammatory events in the lungs. Although they help identify which cell types and in which proportion they contribute to inflammation, they lack the ability to inform on the crosstalk across the different populations, as well as molecular mediators. To overcome this problem, toxicogenomic studies were successfully applied to characterize the MOA of nanomaterials.
Disentangling the direct effects of the exposures from other secondary ones arising from complex regulatory loops is a main challenge in toxicogenomics. Here we assumed that genes whose expression is dose-dependent are directly altered by the exposure at that specific dose. Dose-dependent modeling is indeed an effective method for the identification of transcriptomic alterations that have a monotonic trend across the doses. However, classical dose-dependent analyses are performed at each time point separately, thus limiting the capability of describing the kinetic effects of the exposures.
In this study we performed an integrative analysis to characterize the dynamic dose-dependent portion of the dddtMOA of four ENMs, namely, TiO2, CB, and two MWCNTs. We disentangled the dddtMOA and highlighted specific, functional profiles of each nanomaterial and subcategory. Furthermore, we modeled the commonalities between the dddtMOA and four immune cell populations. We demonstrated that dynamic dose-dependent and differentially expressed genes, that also correlate with the patterns of infiltrations of immune cells in the lungs, are effective in characterizing the direct effects of the exposures. Furthermore, our analyses highlighted specific adaptive responses modulated by different properties of the nanomaterials, such as chemistry and shape. The effect of specific properties on the mediators’ induction has important, potential therapeutic applications.
The integrated dynamic dose-dependent analysis has the potential to highlight key mediators and early responders of exposures. Here we highlighted a panel of 39 genes with dynamic dose-dependent profiles across the nanomaterials. Of these genes, three were correlated to at least one immune cell population, having the potential to inform on lung immunotoxicity.
The present study has several limitations. The recovery after acute toxic response that is highlighted by our analysis might be partially biased by the experimental setup used when performing transcriptomic experiments, which includes only one exposure at the baseline. This could also explain why most of the genes identified are correlated with macrophages and neutrophils, which are known to be early responders against intruders. The availability of a bigger set of ENMs and data from repeated exposures might help in characterizing a chronic adaptive response that would correlate with lymphocytes and eosinophils. Furthermore, the association between physicochemical properties of nanomaterials and the adaptive response may be more robustly detected by including a wider set of nanomaterials.
Nonetheless, the definition of toxicological points of departure for key genes altered by ENM exposures opens new perspectives in the application of transcriptomic data for regulatory purposes.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/nano12122031/s1, Table S1: Table of BAL fluid cell data; Table S2: List of dynamic dose-dependent genes in TiO2; Table S3: List of dynamic dose-dependent genes in CB; Table S4: List of dynamic dose-dependent genes in NM401; Table S5: List of dynamic dose-dependent genes in NRCWE26; Table S6: List of KEGG pathways enriched by the dynamic dose-dependent genes in TiO2; Table S7: List of KEGG pathways enriched by the dynamic dose-dependent genes in CB; Table S8: List of KEGG pathways enriched by the dynamic dose-dependent genes in NM401; Table S9: List of KEGG pathways enriched by the dynamic dose-dependent genes in NRCWE26; Table S10: dynamic dose-dependent labels of the immune cell BAL counts; Table S11: List of KEGG pathways of the dynamic dose-dependent genes correlated with macrophages; Table S12: List of KEGG pathways of the dynamic dose-dependent genes correlated with neutrophils; Table S13: List of KEGG pathways of the dynamic dose-dependent genes correlated with eosinophils; Table S14: List of KEGG pathways of the dynamic dose-dependent genes correlated with lymphocytes.

Author Contributions

Conceptualization, D.G. and U.V.; methodology, A.S.; software, A.S.; validation, A.S., G.d.G. and P.A.S.K.; formal analysis, A.S., G.d.G. and P.A.S.K.; investigation, A.S., G.d.G. and P.A.S.K.; resources, D.G.; data curation, L.A.S., S.S.P., V.F., S.H. and U.V.; writing—original draft preparation, A.S., G.d.G. and P.A.S.K.; writing—review and editing, A.S., G.d.G., L.A.S., S.S.P., V.F., S.H., U.V. and D.G.; visualization, A.S. and G.d.G.; supervision, D.G.; project administration, D.G. and A.S.; funding acquisition, D.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the European Union Horizon 2020 research and innovation programme under grant agreement No 814426 (NanoinformaTIX), and the Academy of Finland (Grant No. 322761).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Publicly available datasets were analyzed in this study. This data can be found here: https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE55286 (accessed on 1 December 2020); https://ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE35193 (accessed on 1 December 2020).

Acknowledgments

The transcriptomic data were generated at Health Canada and were supported by Health Canada’s Genomics Research and Development Initiative.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kinaret, P.; Ilves, M.; Fortino, V.; Rydman, E.; Karisola, P.; Lähde, A.; Koivisto, J.; Jokiniemi, J.; Wolff, H.; Savolainen, K.; et al. Inhalation and Oropharyngeal Aspiration Exposure to Rod-Like Carbon Nanotubes Induce Similar Airway Inflammation and Biological Responses in Mouse Lungs. ACS Nano 2017, 11, 291–303. [Google Scholar] [CrossRef]
  2. Danielsen, P.H.; Knudsen, K.B.; Štrancar, J.; Umek, P.; Koklič, T.; Garvas, M.; Vanhala, E.; Savukoski, S.; Ding, Y.; Madsen, A.M.; et al. Effects of physicochemical properties of TiO2 nanomaterials for pulmonary inflammation, acute phase response and alveolar proteinosis in intratracheally exposed mice. Toxicol. Appl. Pharmacol. 2020, 386, 114830. [Google Scholar] [CrossRef]
  3. Papadiamantis, A.G.; Jänes, J.; Voyiatzis, E.; Sikk, L.; Burk, J.; Burk, P.; Tsoumanis, A.; Ha, M.K.; Yoon, T.H.; Valsami-Jones, E.; et al. Predicting Cytotoxicity of Metal Oxide Nanoparticles using Isalos Analytics Platform. Nanomaterials 2020, 10, 2017. [Google Scholar] [CrossRef]
  4. Welker, L.; Jörres, R.A.; Costabel, U.; Magnussen, H. Predictive value of BAL cell differentials in the diagnosis of interstitial lung diseases. Eur. Respir. J. 2004, 24, 1000–1006. [Google Scholar] [CrossRef] [Green Version]
  5. Hadrup, N.; Rahmani, F.; Jacobsen, N.R.; Saber, A.T.; Jackson, P.; Bengtson, S.; Williams, A.; Wallin, H.; Halappanavar, S.; Vogel, U. Acute phase response and inflammation following pulmonary exposure to low doses of zinc oxide nanoparticles in mice. Nanotoxicology 2019, 13, 1275–1292. [Google Scholar] [CrossRef] [Green Version]
  6. Rahman, L.; Jacobsen, N.R.; Aziz, S.A.; Wu, D.; Williams, A.; Yauk, C.L.; White, P.; Wallin, H.; Vogel, U.; Halappanavar, S. Multi-walled carbon nanotube-induced genotoxic, inflammatory and pro-fibrotic responses in mice: Investigating the mechanisms of pulmonary carcinogenesis. Mutat. Res. Genet. Toxicol. Environ. Mutagen. 2017, 823, 28–44. [Google Scholar] [CrossRef] [Green Version]
  7. Nikota, J.; Banville, A.; Goodwin, L.R.; Wu, D.; Williams, A.; Yauk, C.L.; Wallin, H.; Vogel, U.; Halappanavar, S. Stat-6 signaling pathway and not Interleukin-1 mediates multi-walled carbon nanotube-induced lung fibrosis in mice: Insights from an adverse outcome pathway framework. Part. Fibre Toxicol. 2017, 14, 37. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Rahman, L.; Wu, D.; Johnston, M.; William, A.; Halappanavar, S. Toxicogenomics analysis of mouse lung responses following exposure to titanium dioxide nanomaterials reveal their disease potential at high doses. Mutagenesis 2017, 32, 59–76. [Google Scholar] [CrossRef] [Green Version]
  9. Husain, M.; Kyjovska, Z.O.; Bourdon-Lacombe, J.; Saber, A.T.; Jensen, K.A.; Jacobsen, N.R.; Williams, A.; Wallin, H.; Halappanavar, S.; Vogel, U.; et al. Carbon black nanoparticles induce biphasic gene expression changes associated with inflammatory responses in the lungs of C57BL/6 mice following a single intratracheal instillation. Toxicol. Appl. Pharmacol. 2015, 289, 573–588. [Google Scholar] [CrossRef] [Green Version]
  10. Bengtson, S.; Knudsen, K.B.; Kyjovska, Z.O.; Berthing, T.; Skaug, V.; Levin, M.; Koponen, I.K.; Shivayogimath, A.; Booth, T.J.; Alonso, B.; et al. Differences in inflammation and acute phase response but similar genotoxicity in mice following pulmonary exposure to graphene oxide and reduced graphene oxide. PLoS ONE 2017, 12, e0178355. [Google Scholar] [CrossRef] [Green Version]
  11. Alexander-Dann, B.; Pruteanu, L.L.; Oerton, E.; Sharma, N.; Berindan-Neagoe, I.; Módos, D.; Bender, A. Developments in toxicogenomics: Understanding and predicting compound-induced toxicity from gene expression data. Mol. Omics 2018, 14, 218–236. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Kinaret, P.A.S.; Serra, A.; Federico, A.; Kohonen, P.; Nymark, P.; Liampa, I.; Ha, M.K.; Choi, J.-S.; Jagiello, K.; Sanabria, N.; et al. Transcriptomics in toxicogenomics, part I: Experimental design, technologies, publicly available data, and regulatory aspects. Nanomaterials 2020, 10, 750. [Google Scholar] [CrossRef] [Green Version]
  13. Federico, A.; Serra, A.; Ha, M.K.; Kohonen, P.; Choi, J.-S.; Liampa, I.; Nymark, P.; Sanabria, N.; Cattelani, L.; Fratello, M.; et al. Transcriptomics in toxicogenomics, part II: Preprocessing and differential expression analysis for high quality data. Nanomaterials 2020, 10, 903. [Google Scholar] [CrossRef] [PubMed]
  14. Serra, A.; Fratello, M.; Cattelani, L.; Liampa, I.; Melagraki, G.; Kohonen, P.; Nymark, P.; Federico, A.; Kinaret, P.A.S.; Jagiello, K.; et al. Transcriptomics in toxicogenomics, part III: Data modelling for risk assessment. Nanomaterials 2020, 10, 708. [Google Scholar] [CrossRef] [Green Version]
  15. Sand, S.; Parham, F.; Portier, C.J.; Tice, R.R.; Krewski, D. Comparison of Points of Departure for Health Risk Assessment Based on High-Throughput Screening Data. Environ. Health Perspect. 2017, 125, 623–633. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Phillips, J.R.; Svoboda, D.L.; Tandon, A.; Patel, S.; Sedykh, A.; Mav, D.; Kuo, B.; Yauk, C.L.; Yang, L.; Thomas, R.S.; et al. BMDExpress 2: Enhanced transcriptomic dose-response analysis workflow. Bioinformatics 2019, 35, 1780–1782. [Google Scholar] [CrossRef]
  17. Serra, A.; Saarimäki, L.A.; Fratello, M.; Marwah, V.S.; Greco, D. BMDx: A graphical Shiny application to perform Benchmark Dose analysis for transcriptomics data. Bioinformatics 2020, 36, 2932–2933. [Google Scholar] [CrossRef] [PubMed]
  18. Thomas, R.S.; Allen, B.C.; Nong, A.; Yang, L.; Bermudez, E.; Clewell, H.J.; Andersen, M.E. A method to integrate benchmark dose estimates with genomic data to assess the functional effects of chemical exposure. Toxicol. Sci. 2007, 98, 240–248. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Halappanavar, S.; Rahman, L.; Nikota, J.; Poulsen, S.S.; Ding, Y.; Jackson, P.; Wallin, H.; Schmid, O.; Vogel, U.; Williams, A. Ranking of nanomaterial potency to induce pathway perturbations associated with lung responses. NanoImpact 2019, 14, 100158. [Google Scholar] [CrossRef]
  20. Saarimäki, L.A.; Kinaret, P.A.S.; Scala, G.; del Giudice, G.; Federico, A.; Serra, A.; Greco, D. Toxicogenomics analysis of dynamic dose-response in macrophages highlights molecular alterations relevant for multi-walled carbon nanotube-induced lung fibrosis. NanoImpact 2020, 20, 100274. [Google Scholar] [CrossRef]
  21. Serra, A.; Fratello, M.; Del Giudice, G.; Saarimäki, L.A.; Paci, M.; Federico, A.; Greco, D. TinderMIX: Time-dose integrated modelling of toxicogenomics data. Gigascience 2020, 9, giaa055. [Google Scholar] [CrossRef]
  22. Wither, J.E.; Prokopec, S.D.; Noamani, B.; Chang, N.-H.; Bonilla, D.; Touma, Z.; Avila-Casado, C.; Reich, H.N.; Scholey, J.; Fortin, P.R.; et al. Identification of a neutrophil-related gene expression signature that is enriched in adult systemic lupus erythematosus patients with active nephritis: Clinical/pathologic associations and etiologic mechanisms. PLoS ONE 2018, 13, e0196117. [Google Scholar] [CrossRef]
  23. Poulsen, S.S.; Knudsen, K.B.; Jackson, P.; Weydahl, I.E.K.; Saber, A.T.; Wallin, H.; Vogel, U. Multi-walled carbon nanotube-physicochemical properties predict the systemic acute phase response following pulmonary exposure in mice. PLoS ONE 2017, 12, e0174167. [Google Scholar] [CrossRef] [Green Version]
  24. Billing, A.M.; Knudsen, K.B.; Chetwynd, A.J.; Ellis, L.-J.A.; Tang, S.V.Y.; Berthing, T.; Wallin, H.; Lynch, I.; Vogel, U.; Kjeldsen, F. Fast and robust proteome screening platform identifies neutrophil extracellular trap formation in the lung in response to cobalt ferrite nanoparticles. ACS Nano 2020, 14, 4096–4110. [Google Scholar] [CrossRef]
  25. Saber, A.T.; Mortensen, A.; Szarek, J.; Koponen, I.K.; Levin, M.; Jacobsen, N.R.; Pozzebon, M.E.; Mucelli, S.P.; Rickerby, D.G.; Kling, K.; et al. Epoxy composite dusts with and without carbon nanotubes cause similar pulmonary responses, but differences in liver histology in mice following pulmonary deposition. Part. Fibre Toxicol. 2016, 13, 37. [Google Scholar] [CrossRef] [Green Version]
  26. Poulsen, S.S.; Saber, A.T.; Williams, A.; Andersen, O.; Købler, C.; Atluri, R.; Pozzebon, M.E.; Mucelli, S.P.; Simion, M.; Rickerby, D.; et al. MWCNTs of different physicochemical properties cause similar inflammatory responses, but differences in transcriptional and histological markers of fibrosis in mouse lungs. Toxicol. Appl. Pharmacol. 2015, 284, 16–32. [Google Scholar] [CrossRef]
  27. Bourdon, J.A.; Halappanavar, S.; Saber, A.T.; Jacobsen, N.R.; Williams, A.; Wallin, H.; Vogel, U.; Yauk, C.L. Hepatic and pulmonary toxicogenomic profiles in mice intratracheally instilled with carbon black nanoparticles reveal pulmonary inflammation, acute phase response, and alterations in lipid homeostasis. Toxicol. Sci. 2012, 127, 474–484. [Google Scholar] [CrossRef] [Green Version]
  28. Bourdon, J.A.; Saber, A.T.; Jacobsen, N.R.; Jensen, K.A.; Madsen, A.M.; Lamson, J.S.; Wallin, H.; Møller, P.; Loft, S.; Yauk, C.L.; et al. Carbon black nanoparticle instillation induces sustained inflammation and genotoxicity in mouse lung and liver. Part. Fibre Toxicol. 2012, 9, 5. [Google Scholar] [CrossRef] [Green Version]
  29. Husain, M.; Saber, A.T.; Guo, C.; Jacobsen, N.R.; Jensen, K.A.; Yauk, C.L.; Williams, A.; Vogel, U.; Wallin, H.; Halappanavar, S. Pulmonary instillation of low doses of titanium dioxide nanoparticles in mice leads to particle retention and gene expression changes in the absence of inflammation. Toxicol. Appl. Pharmacol. 2013, 269, 250–262. [Google Scholar] [CrossRef] [Green Version]
  30. Saber, A.T.; Jacobsen, N.R.; Mortensen, A.; Szarek, J.; Jackson, P.; Madsen, A.M.; Jensen, K.A.; Koponen, I.K.; Brunborg, G.; Gützkow, K.B.; et al. Nanotitanium dioxide toxicity in mouse lung is reduced in sanding dust from paint. Part. Fibre Toxicol. 2012, 9, 4. [Google Scholar] [CrossRef] [Green Version]
  31. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef] [PubMed]
  32. Marwah, V.S.; Scala, G.; Kinaret, P.A.S.; Serra, A.; Alenius, H.; Fortino, V.; Greco, D. eUTOPIA: solUTion for Omics data PreprocessIng and Analysis. Source Code Biol. Med. 2019, 14, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Yu, G.; Wang, L.-G.; Han, Y.; He, Q.-Y. clusterProfiler: An R package for comparing biological themes among gene clusters. OMICS 2012, 16, 284–287. [Google Scholar] [CrossRef] [PubMed]
  34. Dewey, M. Metap: Meta-Analysis of Significance Values. 2020. Available online: https://cran.r-project.org/web/packages/metap/metap.pdf (accessed on 1 December 2020).
  35. Luanpitpong, S.; Wang, L.; Rojanasakul, Y. The effects of carbon nanotubes on lung and dermal cellular behaviors. Nanomedicine 2014, 9, 895–912. [Google Scholar] [CrossRef] [Green Version]
  36. Alsaleh, N.B.; Brown, J.M. Immune responses to engineered nanomaterials: Current understanding and challenges. Curr. Opin. Toxicol. 2018, 10, 8–14. [Google Scholar] [CrossRef]
  37. Savi, M.; Rossi, S.; Bocchi, L.; Gennaccaro, L.; Cacciani, F.; Perotti, A.; Amidani, D.; Alinovi, R.; Goldoni, M.; Aliatis, I.; et al. Titanium dioxide nanoparticles promote arrhythmias via a direct interaction with rat cardiac tissue. Part. Fibre Toxicol. 2014, 11, 63. [Google Scholar] [CrossRef] [Green Version]
  38. Husain, M.; Wu, D.; Saber, A.T.; Decan, N.; Jacobsen, N.R.; Williams, A.; Yauk, C.L.; Wallin, H.; Vogel, U.; Halappanavar, S. Intratracheally instilled titanium dioxide nanoparticles translocate to heart and liver and activate complement cascade in the heart of C57BL/6 mice. Nanotoxicology 2015, 9, 1013–1022. [Google Scholar] [CrossRef] [Green Version]
  39. Kinaret, P.A.S.; Ndika, J.; Ilves, M.; Wolff, H.; Vales, G.; Norppa, H.; Savolainen, K.; Skoog, T.; Kere, J.; Moya, S.; et al. Toxicogenomic profiling of 28 nanomaterials in mouse airways. Adv. Sci. 2021, 8, 2004588. [Google Scholar] [CrossRef]
  40. Nikota, J.; Williams, A.; Yauk, C.L.; Wallin, H.; Vogel, U.; Halappanavar, S. Meta-analysis of transcriptomic responses as a means to identify pulmonary disease outcomes for engineered nanomaterials. Part. Fibre Toxicol. 2016, 13, 25. [Google Scholar] [CrossRef] [Green Version]
  41. Poulsen, S.S.; Jackson, P.; Kling, K.; Knudsen, K.B.; Skaug, V.; Kyjovska, Z.O.; Thomsen, B.L.; Clausen, P.A.; Atluri, R.; Berthing, T.; et al. Multi-walled carbon nanotube physicochemical properties predict pulmonary inflammation and genotoxicity. Nanotoxicology 2016, 10, 1263–1275. [Google Scholar] [CrossRef] [Green Version]
  42. Knudsen, K.B.; Berthing, T.; Jackson, P.; Poulsen, S.S.; Mortensen, A.; Jacobsen, N.R.; Skaug, V.; Szarek, J.; Hougaard, K.S.; Wolff, H.; et al. Physicochemical predictors of Multi-Walled Carbon Nanotube-induced pulmonary histopathology and toxicity one year after pulmonary deposition of 11 different Multi-Walled Carbon Nanotubes in mice. Basic Clin. Pharmacol. Toxicol. 2019, 124, 211–227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Gaté, L.; Knudsen, K.B.; Seidel, C.; Berthing, T.; Chézeau, L.; Jacobsen, N.R.; Valentino, S.; Wallin, H.; Bau, S.; Wolff, H.; et al. Pulmonary toxicity of two different multi-walled carbon nanotubes in rat: Comparison between intratracheal instillation and inhalation exposure. Toxicol. Appl. Pharmacol. 2019, 375, 17–31. [Google Scholar] [CrossRef] [PubMed]
  44. Panigrahi, A.R.; Yadav, P.; Beura, S.K.; Singh, S.K. Blood Coagulation System and Carbon-Based Nanoengineering for Biomedical Application. In Advanced Micro- and Nano-Manufacturing Technologies; Springer: Berlin/Heidelberg, Germany, 2022; pp. 279–298. [Google Scholar]
  45. Sharma, M.; Salisbury, R.L.; Maurer, E.I.; Hussain, S.M.; Sulentic, C.E.W. Gold nanoparticles induce transcriptional activity of NF-κB in a B-lymphocyte cell line. Nanoscale 2013, 5, 3747–3756. [Google Scholar] [CrossRef]
  46. Kinaret, P.A.S.; Scala, G.; Federico, A.; Sund, J.; Greco, D. Carbon nanomaterials promote M1/M2 macrophage activation. Small 2020, 16, e1907609. [Google Scholar] [CrossRef] [Green Version]
  47. Rydman, E.M.; Ilves, M.; Koivisto, A.J.; Kinaret, P.A.S.; Fortino, V.; Savinko, T.S.; Lehto, M.T.; Pulkkinen, V.; Vippola, M.; Hämeri, K.J.; et al. Inhalation of rod-like carbon nanotubes causes unconventional allergic airway inflammation. Part. Fibre Toxicol. 2014, 11, 48. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Shvedova, A.A.; Kisin, E.; Murray, A.R.; Johnson, V.J.; Gorelik, O.; Arepalli, S.; Hubbs, A.F.; Mercer, R.R.; Keohavong, P.; Sussman, N.; et al. Inhalation vs. aspiration of single-walled carbon nanotubes in C57BL/6 mice: Inflammation, fibrosis, oxidative stress, and mutagenesis. Am. J. Physiol. Lung Cell. Mol. Physiol. 2008, 295, L552–L565. [Google Scholar] [CrossRef] [Green Version]
  49. Rosales, C.; Lowell, C.A.; Schnoor, M.; Uribe-Querol, E. Neutrophils: Their role in innate and adaptive immunity 2017. J. Immunol. Res. 2017, 2017, 9748345. [Google Scholar] [CrossRef]
  50. Akuthota, P.; Wang, H.B.; Spencer, L.A.; Weller, P.F. Immunoregulatory roles of eosinophils: A new look at a familiar cell. Clin. Exp. Allergy 2008, 38, 1254–1263. [Google Scholar] [CrossRef]
  51. Zhang, X.; Mosser, D.M. Macrophage activation by endogenous danger signals. J. Pathol. 2008, 214, 161–178. [Google Scholar] [CrossRef]
  52. Giese, M.A.; Hind, L.E.; Huttenlocher, A. Neutrophil plasticity in the tumor microenvironment. Blood 2019, 133, 2159–2167. [Google Scholar] [CrossRef]
  53. Liu, Y.; Hardie, J.; Zhang, X.; Rotello, V.M. Effects of engineered nanoparticles on the innate immune system. Semin. Immunol. 2017, 34, 25–32. [Google Scholar] [CrossRef] [PubMed]
  54. Kim, S.Y.; Nair, M.G. Macrophages in wound healing: Activation and plasticity. Immunol. Cell Biol. 2019, 97, 258–267. [Google Scholar] [CrossRef] [PubMed]
  55. Yu, W.P.; Ding, J.L.; Liu, X.L.; Zhu, G.D.; Lin, F.; Xu, J.J.; Wang, Z.; Zhou, J.L. Titanium dioxide nanotubes promote M2 polarization by inhibiting macrophage glycolysis and ultimately accelerate endothelialization. Immun. Inflamm. Dis. 2021, 9, 746–757. [Google Scholar] [CrossRef] [PubMed]
  56. Xu, W.-C.; Dong, X.; Ding, J.-L.; Liu, J.-C.; Xu, J.-J.; Tang, Y.-H.; Yi, Y.-P.; Lu, C.; Yang, W.; Yang, J.-S.; et al. Nanotubular TiO2 regulates macrophage M2 polarization and increases macrophage secretion of VEGF to accelerate endothelialization via the ERK1/2 and PI3K/AKT pathways. Int. J. Nanomed. 2019, 14, 441–455. [Google Scholar] [CrossRef] [Green Version]
  57. Braun, D.A.; Fribourg, M.; Sealfon, S.C. Cytokine response is determined by duration of receptor and signal transducers and activators of transcription 3 (STAT3) activation. J. Biol. Chem. 2013, 288, 2986–2993. [Google Scholar] [CrossRef] [Green Version]
  58. Velazquez-Salinas, L.; Verdugo-Rodriguez, A.; Rodriguez, L.L.; Borca, M.V. The role of interleukin 6 during viral infections. Front. Microbiol. 2019, 10, 1057. [Google Scholar] [CrossRef] [Green Version]
  59. Bozza, M.T.; Jeney, V. Pro-inflammatory Actions of Heme and Other Hemoglobin-Derived DAMPs. Front. Immunol. 2020, 11, 1323. [Google Scholar] [CrossRef]
  60. Rifkind, J.M.; Mohanty, J.G.; Nagababu, E. The pathophysiology of extracellular hemoglobin associated with enhanced oxidative reactions. Front. Physiol. 2014, 5, 500. [Google Scholar] [CrossRef] [Green Version]
  61. Futosi, K.; Fodor, S.; Mócsai, A. Neutrophil cell surface receptors and their intracellular signal transduction pathways. Int. Immunopharmacol. 2013, 17, 638–650. [Google Scholar] [CrossRef] [Green Version]
  62. Uddin, M.; Lau, L.C.; Seumois, G.; Vijayanand, P.; Staples, K.J.; Bagmane, D.; Cornelius, V.; Dorinsky, P.; Davies, D.E.; Djukanović, R. EGF-induced bronchial epithelial cells drive neutrophil chemotactic and anti-apoptotic activity in asthma. PLoS ONE 2013, 8, e72502. [Google Scholar] [CrossRef] [Green Version]
  63. Danielsen, P.H.; Bendtsen, K.M.; Knudsen, K.B.; Poulsen, S.S.; Stoeger, T.; Vogel, U. Nanomaterial- and shape-dependency of TLR2 and TLR4 mediated signaling following pulmonary exposure to carbonaceous nanomaterials in mice. Part. Fibre Toxicol. 2021, 18, 40. [Google Scholar] [CrossRef]
  64. Wang, Y.-H.; Wills-Karp, M. The potential role of interleukin-17 in severe asthma. Curr. Allergy Asthma Rep. 2011, 11, 388–394. [Google Scholar] [CrossRef] [Green Version]
  65. Ji, Z.; LeBaron, M.J.; Schisler, M.R.; Zhang, F.; Bartels, M.J.; Gollapudi, B.B.; Pottenger, L.H. Dose-Response for Multiple Biomarkers of Exposure and Genotoxic Effect Following Repeated Treatment of Rats with the Alkylating Agents, MMS and MNU. Mutagenesis 2016, 31, 297–308. [Google Scholar] [CrossRef] [Green Version]
  66. Okada, T.; Lee, B.W.; Ogami, A.; Oyabu, T.; Myojo, T. Inhalation of titanium dioxide (P25) nanoparticles to rats and changes in surfactant protein (SP-D) levels in bronchoalveolar lavage fluid and serum. Nanotoxicology 2019, 13, 1396–1408. [Google Scholar] [CrossRef]
  67. Honda, H.; Fujimoto, M.; Miyamoto, S.; Ishikawa, N.; Serada, S.; Hattori, N.; Nomura, S.; Kohno, N.; Yokoyama, A.; Naka, T. Sputum Leucine-Rich Alpha-2 Glycoprotein as a Marker of Airway Inflammation in Asthma. PLoS ONE 2016, 11, e0162672. [Google Scholar] [CrossRef]
  68. Ford, J.; Hughson, A.; Lim, K.; Bardina, S.V.; Lu, W.; Charo, I.F.; Lim, J.K.; Fowell, D.J. CCL7 Is a Negative Regulator of Cutaneous Inflammation Following Leishmania major Infection. Front. Immunol. 2018, 9, 3063. [Google Scholar] [CrossRef]
  69. Szymczak, W.A.; Deepe, G.S. The CCL7-CCL2-CCR2 axis regulates IL-4 production in lungs and fungal immunity. J. Immunol. 2009, 183, 1964–1974. [Google Scholar] [CrossRef] [Green Version]
  70. Zhang, J.; Veeramachaneni, N. Targeting interleukin-1β and inflammation in lung cancer. Biomark. Res. 2022, 10, 5. [Google Scholar] [CrossRef]
  71. Cross, L.J.M.; Matthay, M.A. Biomarkers in acute lung injury: Insights into the pathogenesis of acute lung injury. Crit. Care Clin. 2011, 27, 355–377. [Google Scholar] [CrossRef] [Green Version]
  72. Thompson, W.L.; Van Eldik, L.J. Inflammatory cytokines stimulate the chemokines CCL2/MCP-1 and CCL7/MCP-3 through NFkB and MAPK dependent pathways in rat astrocytes [corrected]. Brain Res. 2009, 1287, 47–57. [Google Scholar] [CrossRef] [Green Version]
Nanomaterials 12 02031 g001 550
Figure 1. (A) Number of dose-dependent genes for each ENM for each of the dose-time labels. The labels can be interpreted as follows: sensitive—genes that respond at low doses; intermediate—genes that respond at intermediate doses; resilient—genes that respond at high doses.; early—genes that respond at short time points; middle—genes that respond at intermediate time points; late—genes that respond at late time points. The blue color gradient describes the number of dose-dependent genes: dark blue represents a smaller number of genes identified, whereas light blue represents a higher number of dose-dependent genes. (B) Number of DEG (red) and genes with a dPOD activity label (green) for each ENM.
Figure 1. (A) Number of dose-dependent genes for each ENM for each of the dose-time labels. The labels can be interpreted as follows: sensitive—genes that respond at low doses; intermediate—genes that respond at intermediate doses; resilient—genes that respond at high doses.; early—genes that respond at short time points; middle—genes that respond at intermediate time points; late—genes that respond at late time points. The blue color gradient describes the number of dose-dependent genes: dark blue represents a smaller number of genes identified, whereas light blue represents a higher number of dose-dependent genes. (B) Number of DEG (red) and genes with a dPOD activity label (green) for each ENM.
Nanomaterials 12 02031 g001
Nanomaterials 12 02031 g002 550
Figure 2. Functions of the dddtMOA shared across ENM categories (A,B,C,H) and specific to each ENM (D,E,F,G).
Figure 2. Functions of the dddtMOA shared across ENM categories (A,B,C,H) and specific to each ENM (D,E,F,G).
Nanomaterials 12 02031 g002
Nanomaterials 12 02031 g003 550
Figure 3. (A) Number of dynamic dose-dependent genes correlated with cell counts for each ENM. (B) Proportion of dynamic dose-dependent genes correlated with cell counts for each ENM.
Figure 3. (A) Number of dynamic dose-dependent genes correlated with cell counts for each ENM. (B) Proportion of dynamic dose-dependent genes correlated with cell counts for each ENM.
Nanomaterials 12 02031 g003
Nanomaterials 12 02031 g004 550
Figure 4. KEGG pathways of the dynamic dose-dependent genes for the different ENMs. The different panels represent the pathways enriched by the dose-dependent genes that for each ENM correlate with macrophages (A), neutrophils (B), eosinophils (C), and lymphocytes (D). The size of the dots represent the gene ration whereas the color represent the adjusted p-values of the enrichment analysis. The numbers in brackets below the ENM names represent the number of genes correlated with the cell populations.
Figure 4. KEGG pathways of the dynamic dose-dependent genes for the different ENMs. The different panels represent the pathways enriched by the dose-dependent genes that for each ENM correlate with macrophages (A), neutrophils (B), eosinophils (C), and lymphocytes (D). The size of the dots represent the gene ration whereas the color represent the adjusted p-values of the enrichment analysis. The numbers in brackets below the ENM names represent the number of genes correlated with the cell populations.
Nanomaterials 12 02031 g004
Nanomaterials 12 02031 g005 550
Figure 5. Graphic representation of the effect of ENM psychochemical properties on immune cell populations. In (A), the effect of different-shaped ENMs on neutrophils is reported by showing the different mediators and signaling pathways activated. In (B), ENMs with different chemistry are shown to differently polarize macrophages.
Figure 5. Graphic representation of the effect of ENM psychochemical properties on immune cell populations. In (A), the effect of different-shaped ENMs on neutrophils is reported by showing the different mediators and signaling pathways activated. In (B), ENMs with different chemistry are shown to differently polarize macrophages.
Nanomaterials 12 02031 g005
Nanomaterials 12 02031 g006 550
Figure 6. Genes with a dynamic dose-dependent transcriptomic profile in the four ENMs. Colors indicate the dynamic point of departure of the gene with respect to early, middle, and late time points.
Figure 6. Genes with a dynamic dose-dependent transcriptomic profile in the four ENMs. Colors indicate the dynamic point of departure of the gene with respect to early, middle, and late time points.
Nanomaterials 12 02031 g006
Table
Table 1. Datasets used in this study.
Table 1. Datasets used in this study.
ENM IDTypeLength/Diameter (nm)Surface Area (m2/g)No. SamplesDoses (µg)Time (Day)GSERef.
NM401Multiwalled carbon nanotube4048/671813918/54/1621/3/28GSE55286[26]
NRCWE-026Multiwalled carbon nanotube847/1024513918/54/1621/3/28GSE55286[26]
CB (Printex 90)Carbon black14 (diameter)295–3386718/54/1621/3/28GSE35193[27,28]
TiO2 (L181 UVTitan)Nano-TiO220.6 (diameter)107.76518/54/1621/3/28GSE41041[29,30]
Table
Table 2. Genes that are dose-dependent in all the ENMs and are correlated to at least one cell population. In dPOD column: RE stands for resilient-early, SE stands for sensitive-early. In cell population column: N stands for neutrophils, M stands for macrophages, E stands for eosinophils.
Table 2. Genes that are dose-dependent in all the ENMs and are correlated to at least one cell population. In dPOD column: RE stands for resilient-early, SE stands for sensitive-early. In cell population column: N stands for neutrophils, M stands for macrophages, E stands for eosinophils.
GeneDescriptionENMdPODCell Population
CCL7chemokine (C-C motif) ligand 7TiO2REN, M
CBSEN
NM401SEM
NRCWE26SEN
CCL12chemokine (C-C motif) ligand 12TiO2REE
CBREN
NM401SEN
NRCWE26SEN
IL1binterleukin 1 betaTiO2SEE
CBSEM
NM401SEM
NRCWE26SEM, E
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Serra, A.; del Giudice, G.; Kinaret, P.A.S.; Saarimäki, L.A.; Poulsen, S.S.; Fortino, V.; Halappanavar, S.; Vogel, U.; Greco, D. Characterization of ENM Dynamic Dose-Dependent MOA in Lung with Respect to Immune Cells Infiltration. Nanomaterials 2022, 12, 2031. https://doi.org/10.3390/nano12122031

AMA Style

Serra A, del Giudice G, Kinaret PAS, Saarimäki LA, Poulsen SS, Fortino V, Halappanavar S, Vogel U, Greco D. Characterization of ENM Dynamic Dose-Dependent MOA in Lung with Respect to Immune Cells Infiltration. Nanomaterials. 2022; 12(12):2031. https://doi.org/10.3390/nano12122031

Chicago/Turabian Style

Serra, Angela, Giusy del Giudice, Pia Anneli Sofia Kinaret, Laura Aliisa Saarimäki, Sarah Søs Poulsen, Vittorio Fortino, Sabina Halappanavar, Ulla Vogel, and Dario Greco. 2022. "Characterization of ENM Dynamic Dose-Dependent MOA in Lung with Respect to Immune Cells Infiltration" Nanomaterials 12, no. 12: 2031. https://doi.org/10.3390/nano12122031

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop