Skip to main content
Download PDF

Comprehensive proteome analysis of nasal lavage samples after controlled exposure to welding nanoparticles shows an induced acute phase and a nuclear receptor, LXR/RXR, activation that influence the status of the extracellular matrix

Clinical Proteomics volume 15, Article number: 20 (2018) Cite this article

Abstract

Background

Epidemiological studies have shown that many welders experience respiratory symptoms. During the welding process a large number of airborne nanosized particles are generated, which might be inhaled and deposited in the respiratory tract. Knowledge of the underlying mechanisms behind observed symptoms is still partly lacking, although inflammation is suggested to play a central role. The aim of this study was to investigate the effects of welding fume particle exposure on the proteome expression level in welders suffering from respiratory symptoms, and changes in protein mediators in nasal lavage samples were analyzed. Such mediators will be helpful to clarify the pathomechanisms behind welding fume particle-induced effects.

Methods

In an exposure chamber, 11 welders with work-related symptoms in the lower airways during the last month were exposed to mild-steel welding fume particles (1 mg/m3) and to filtered air, respectively, in a double-blind manner. Nasal lavage samples were collected before, immediately after, and the day after exposure. The proteins in the nasal lavage were analyzed with two different mass spectrometry approaches, label-free discovery shotgun LC–MS/MS and a targeted selected reaction monitoring LC–MS/MS analyzing 130 proteins and four in vivo peptide degradation products.

Results

The analysis revealed 30 significantly changed proteins that were associated with two main pathways; activation of acute phase response signaling and activation of LXR/RXR, which is a nuclear receptor family involved in lipid signaling. Connective tissue proteins and proteins controlling the degradation of such tissues, including two different matrix metalloprotease proteins, MMP8 and MMP9, were among the significantly changed enzymes and were identified as important key players in the pathways.

Conclusion

Exposure to mild-steel welding fume particles causes measurable changes on the proteome level in nasal lavage matrix in exposed welders, although no clinical symptoms were manifested. The results suggested that the exposure causes an immediate effect on the proteome level involving acute phase proteins and mediators regulating lipid signaling. Proteases involved in maintaining the balance between the formation and degradation of extracellular matrix proteins are important key proteins in the induced effects.

Background

Occupations with work tasks that result in high emissions of airborne particles pose a risk for respiratory health effects according to the national institution of occupational safety and health, occupational safety and health administration, and the world health organization [1,2,3,4]. Welding exposure has been shown to induce a number of negative health effects and symptoms, including bronchitis, siderosis, asthma, metal fume fever, lung function changes, immune suppression, pulmonary infection, and cancer [5, 6]. During welding, base materials and a filler material are fused at high temperature using different methods, e.g. gas metal arc welding and manual metal arc welding and a mixture of gases and agglomerated particles (networks of interconnected particles) are released into the air [7, 8]. The sizes of the agglomerates (typically ~ 100–1000 nm) are in the respirable range, which means that they can be inhaled and exert potentially negative airway effects. The agglomerates are typically built up of primary nanoparticles, which are in the range of 2–70 nm in diameter [9]. The welding fume particles consist of a complex mixture of different metals, e.g. iron, manganese metal oxides [10], and other trace elements depending on the type of welding method and the electrode material, e.g. mild steel or stainless steel [8].

It has been shown that although the exposure levels do not normally exceed current Swedish permissible occupational exposure limits for inorganic respiratory dust (5 mg/m3), there is a high frequency of upper and lower respiratory symptoms among welders [11, 12]. Previous studies have suggested that welding particle exposure induces oxidative stress and an inflammatory response [13, 14]. Elevated levels of tumor necrosis factor-α (TNF-α), interleukin (IL)-1, IL-6, and IL-8 have been observed after welding exposure [15,16,17]. It has also been shown that even at low exposure levels in a controlled human chamber study where no clinical symptoms were recorded, biomarkers of inflammation were still induced [18]. There is a need for a deeper understanding of the mechanisms involved in the responses induced by welding fume particle exposure, and the identification of proteins and pathways associated with welding fume particle exposure is an important step in this regard.

A combination of different mass spectrometry techniques can rapidly and comprehensively identify and quantify protein contents in highly complex samples [19]. During the discovery phase, shotgun proteomics can serve as a valuable tool in rapidly generating a global profile of the protein changes linked to a specific event. Processing the results from the discovery analysis with the bioinformatics-based ingenuity pathway analysis (IPA) can generate hypotheses on the mechanistic pathways that are induced [20]. Measurement of the identified proteins with a targeted LC-coupled selected reaction monitoring (LC–SRM) method offers a better opportunity to validate multiple biomarker candidates simultaneously because it offers high-throughput analysis and is able to provide label-free relative quantification as well as absolute protein quantification.

The aim of this study was to investigate the effects on the proteome level induced by welding fume particles at exposure levels well below the present Swedish permissible occupational exposure limit in welders with work-related lower airway symptoms.

Methods

Materials

Amicon Ultra-0.5 centrifugal filters (cut-off 3 kDa) were purchased from Millipore (Bedford, MA, USA). The micro bicinchoninic acid (BCA) protein assay kit was from Thermo Scientific (Rockford, IL, USA). Trypsin (sequencing grade) was purchased from Promega Corporation (Madison, WI, USA). Calcium chloride, formic acid (FA), hydrochloric acid, and ammonium acetate were purchased from Merck (Darmstadt, Germany). Dithiothreitol (DTT) and iodoacetamide were purchased from Sigma-Aldrich (St. Louis, MO, USA). Acetonitrile (ACN) was purchased from Lab-scan (Dublin, Ireland).

Study subjects

Welders (n = 11) with work-related lower airway symptoms (wheezing, dyspnoea, and/or cough) in the last month were recruited for the study. The inclusion criteria are described in more detail by Dierschke et al. [13]. Briefly, the study subjects were all non-smoking men ranging from 29 to 66 years of age. Baseline medical examinations of the participants were conducted and included standardized questions about past and current respiratory and cardiovascular diseases and allergies. Work-related symptoms (defined as symptoms that recovered during weekends and holidays) from the upper airway at least once a week during the last month were also recorded. Six welders reported work-related nasal symptoms (runny nose, itching/sneezing) during the last month, and three welders showed a positive skin prick test for a standard panel of aeroallergens. Furthermore, lung function (FVC and FEV1 as % of predicted) and methacholine test (PD15) were carried out [18]. For the whole group the lung function was normal: FVC: 99 (72–114) and FEV1: 101 (66–132; % of predicted, median, min–max). Two out of 10 subjects had signs of bronchial hyperreactivity (PD15 < 167 µg methacholine). One welder started coughing during the test, thus it could not be completed. The two hyperreactive welders were among the three subjects with slightly decreased lung function.

Study design

The subjects were exposed in an exposure chamber (22 m3) for 5.5 h with a lunch break for 1 h after 2.5 h of exposure to either filtered air or welding fumes at a concentration of 1 mg/m3 (PM2,5) as described previously [21]. Briefly, particles were generated by gas-metal arc welding in mild steel and collected in a closed chamber from where they were diluted in fresh air via a carbon filter for removal of gases and then fed into the exposure chamber. The composition of the welding fume particles was mainly iron oxides and up to 20% manganese. The primary particle size ranged from 2 to 70 nm, and aggregates with a mean mobility diameter of 160 nm were formed. The exposures were conducted in a double blind fashion and were performed on two Mondays with 1 week in between, making sure that no occupational exposure to welding occurred for at least 48 h prior to the experimental exposure. Nasal lavage fluids were collected on four occasions, twice prior to exposure and twice after, by instilling the nasal cavity with 18 mL of an isotonic saline solution. The first sample was a washout (NL 0, not analyzed), and the second was the baseline sample (NL 1) that was taken 30 min after the wash out. The third sample was collected immediately at the end of the exposure (NL 2), and the fourth was taken at 18–20 h after the end of the exposure (NL 3). An overview of the exposure and sample collection is given in Fig. 1. All samples were stored at − 80 °C until analysis. Two subjects had missing samples and were therefore excluded in the analysis of the pooled samples.

Fig. 1
figure 1

Overview of study design and nasal lavage sampling. The subjects were exposed for 5.5 h to either filtered air or welding fumes at 1 mg/m3 (PM2.5). Nasal lavage fluids were collected on four occasions, twice prior to exposure and twice after. The first sample was a washout (NL 0, not analyzed), the second was the baseline sample (NL 1), which was taken 30 min after the wash out, the third sample was collected immediately at the end of exposure (NL 2), and the fourth sample was taken 18–20 h after the end of the exposure (NL 3). Two subjects had missing samples and were therefore excluded in the analysis of the pooled samples

Full size image

Samples

Mechanistic pathways were first investigated by analyses of pooled samples on a group level to look for indications of global protein changes using a shotgun proteomic hypothesis-generating method. Findings from the shotgun analysis were further investigated by analyses of individual samples using a targeted label-free analysis with relative quantification of all measured proteins and with absolute quantification of five specific proteins that had been indicated by the shotgun analysis to be associated to the exposure (Fig. 2).

Fig. 2
figure 2

Flowchart of the experiment: Nasal lavage sampling before and after exposure to welding fume particles and filtered air. (1) Protein analysis through shotgun proteomics in combination with IPA analysis generates a hypothesis of the mechanisms induced by the welding fume particle exposure. (2) Targeted protein analysis of specific proteins with label-free relative quantification and absolute quantification of MMP9, A1AT, A2MG, APOB, and elastin. (3) Statistical evaluation of the changed protein level with linear mixed model and Wilcoxon signed rank test. (4) Pathway analysis of significantly changed proteins after welding fume particle exposure compared to filtered air exposure

Full size image

The total protein content in each nasal lavage fluid sample was determined using a BCA protein assay kit.

Analysis of pooled samples with discovery shotgun method

Equal amounts of total protein (50 µg) were pooled from all subjects with samples collected at all time points (n = 9), resulting in six pooled samples. The six samples resulted from the samples collected at the three time points (NL 1, NL 2, and NL 3) and for both types of exposure (welding fume exposure or filtered air exposure) (Fig. 1). Samples were reduced with DTT (1.25 mg/mg protein) for 1 h at 60 °C and then alkylated with iodacetamide (2.5 mg/mg protein) for 30 min at room temperature and in darkness. Samples were desalted in centrifugal filters (cut-off 3 kDa) at 4500×g for 2 h and then washed with 300 µl 50 mM ammonium acetate three times followed by centrifugation at 4500×g for 30 min. Proteins were then trypsin digested (1:20 trypsin:protein) overnight at 37 °C with Promega trypsin and 1 mM calcium chloride.

Instrumentation and data analysis

Protein digests (250 ng/injection) were loaded into a nano liquid chromatography system (EASYnLC, Thermo Scientific, Waltham, MA, USA) with a C18 column (100 mm × 0.75 μm, Agilent Technologies, Santa Clara, CA, USA) coupled to a high-resolution mass spectrometer (Orbitrap Velos Pro, Thermo Fisher). Peptides were separated using a 70 min increase from 2 to 40% ACN, supplemented with 0.1% FA, followed by an increase to 90% ACN for 10 min. Peptides were analyzed by a data-dependent acquisition method utilizing collision-induced dissociation for sequencing on an LTQ Velos Orbitrap Pro mass spectrometer (Thermo Scientific). Identification and relative quantification was performed with the MAXQuant software against a human database downloaded from Uniprot (October 31, 2015). Protein levels were determined using label-free quantification based on peptide intensity with at least two peptides (unique and razor, with at least one unique peptide). Relative protein quantification was performed by normalizing each protein level against the protein level found in the baseline sample (the nasal lavage sample collected before exposure) for each exposure, thus the ratios were NL2/NL1 and NL3/NL1, respectively. For a protein to be further evaluated with the targeted method, it had to be detected in at least 50% of the samples (in at least 3 of 6) and have a ratio for being categorized as an increased or decreased protein level of > 1.3 or < 0.8, respectively. When selections according to these criteria were made, a total of 137 welding-associated proteins were identified (Additional file 1). These proteins were further analyzed with IPA (Additional files 2, 3, 4).

Analysis of individual samples with the targeted SRM method

The samples were evaporated and dissolved in 50 mM ammonium acetate to a concentration of 400 µg/mL. Proteins were reduced, alkylated, and trypsin digested. All samples were spiked with isotopically labeled standards from matrix metalloproteinase (MMP)9 (4 fmol/µL), alpha-1-antitrypsin (A1AT) (4 fmol/µL), alpha-2-macroglobulin (A2MG) (4 fmol/µL), apolipoprotein (APO)B (4 fmol/µL), and elastin (120 fmol/µL). The digested samples were desalted on a solid-phase extraction column (Evolute ABN 25 mg packing material, 1 mL column, Biotag, Uppsala, Sweden) by loading the digested sample onto the column according to the instructions from the manufacturer. The peptides were eluted with a gradient of 200 mL 30% ACN, 40% ACN, 70% ACN (and 0.1% FA), and 90% ACN, and the solution was evaporated to dryness (SpeedVac, HT-4X Evaporator, Genevac Inc, USA) and reconstituted in 200 μL 0.1% FA (≈ 1 mg/mL). Samples were analyzed twice (20 μg protein per injection). To avoid systematic bias, all samples originating from the same subject but from different time points were prepared and analyzed at the same time.

Instrumentation

The SRM analyses were conducted using a 5500 QTRAP hybrid triple quadrupole/linear ion trap mass spectrometer equipped with a TurboIonSpray source (Applied Biosystems/MDS Sciex, Framingham, MA, USA) connected to an LC system (UFLCXR; Shimadzu Corporation, Kyoto, Japan; LC–SRM–MS). The peptides were separated on a C18-column (XBridge peptide BEH C18, 130 Å, 3.5 µm particles, 2.1 mm × 100 mm, Waters Corporation, Milford, Massachusetts, USA) in a 32-min run using a gradient flow of 0.4 mL/min with 0.1% FA in water (A) and 0.1% FA in ACN (B) as mobile phases. The LC gradient was increased from 5 B to 10% B for 10 min, from 10 B to 25% B for 15 min, and from 25 B to 95% B for 1 min where it was kept for 2 min. The column was recalibrated to 5% B for 4 min before the next sample injection.

SRM assays

Selection of SRM assays

Addition of nasal lavage protein SRM assays were added to a method already developed in a previous study [22]. In total, SRM assays for 77 proteins were used from the original SRM method, and new SRM assays were developed for an additional 53 proteins (Additional file 5). Unique peptides were chosen for the SRM assays with at least three transitions per peptide and 1–4 assays per protein were used. The tool UniMaP (http://bioserver-1.bioacademy.gr/Bioserver/UniMaP/Tool_index.php) was used to control if all selected peptides were unique for the targeted protein. The added SRM assays were generated by measuring nasal lavage protein digests in full scan and product ion scan with LC–SRM–MS/MS with a dwell time of 10, 4, and 3 ms. Collision energies were generated in Skyline software according to the earlier developed protocol [22]. Robust transitions with high intensity and good reproducibility [less than 20% coefficient of variation (CV)] were further selected.

The final SRM method contained 256 precursors with 811 transitions targeting 130 proteins and 4 in vivo peptide degradation products from collage 4A1 and elastin (the peptide degradation products were previously determined [23, 24]) (Additional file 6).

Selection and optimization of SRM assays for proteins selected for absolute quantification

The proteins MMP9, A1AT, APOB, A2MG, and elastin were selected for absolute quantification. The peptides selected to represent the protein were chosen based on a minimal probability to be missed or partly cleaved or to be biologically or chemically modified [25]. One peptide for each protein was selected, and this peptide (the standard peptide) and the corresponding isotope-labeled peptide (to be used as the internal standard) were purchased from New England Peptide (Gardner, MA, USA).

The collision energy was optimized, and the ion fragments with the highest intensity were selected (Additional file 7). Sensitivity for the selected SRM transitions was assessed by determining the peak area of the quantifier ion of the targeted peptide after injecting 100 fmol of the standard and the internal standards for MMP9, A1AT, A2MG, APOB, and elastin via infusion. The impact of the matrix was assessed by comparing the chromatographic peaks of the transitions in the presence and absence of biological matrix (pooled nasal lavage from 30 volunteers) to assess whether tryptic peptides of nasal lavage proteins interfered with the intended target peptides.

Absolute quantification of five proteins

Calibration curves with and without biological matrix were prepared at concentrations of 30, 23, 15, 7.6, 3.8, 1.5, 0.76, 0.38, 0.076, and 0 fmol/µL. The samples for each calibration point were prepared separately and analyzed twice by LC–MS/MS. All calibration samples were spiked with isotopically labeled standards for MMP9 (4 fmol/µL), A1AT (4 fmol/µL), A2MG (4 fmol/µL), APOB (4 fmol/µL), and elastin (120 fmol/µL) (Additional file 8).

Determination of the limit of detection (LOD)

LOD was determined by calculating the concentration corresponding to three times the standard deviation of the peak area ratios between the peak area at noise level and the peak area of each corresponding internal standard of eight blank samples and thus LOD in this case corresponds to the noise level that prevailed.

Reproducibility of the method

The stability of the sample preparation and the reproducibility of the mass spectrometry analyses were determined by analysis of quality control samples. The samples were composed of pooled nasal lavage samples spiked with the different isotopically labeled peptides MMP9 (4 fmol/µL), A1AT (4 fmol/µL), A2MG (4 fmol/µL), APOB (4 fmol/µL), and elastin (120 fmol/µL) (Additional file 9).

Data processing of individual samples

Data were processed in an automated manner with the Skyline software version 3.5. Peaks were integrated automatically with the algorithm of the Skyline software.

Relative quantification

In total, 256 peptides measuring 130 proteins and 4 in vivo peptide degradation products were analyzed on the individual level with targeted SRM. The isotopically labeled standards were used as global normalizers for systemic error correction. The total peak area of each SRM assay was normalized to the total peak areas of each isotopically labeled reference peptide MMP9, A1AT, APOB, and A2MG, respectively (the isotopically labeled peptide for A2MG was represented by both double (927 m/z) and triple charged (618 m/z) precursors, using each of the precursors as a global normalizer separately). Peptides were thus normalized with all four different isotope-labeled peptides, and data from each normalization step and the non-normalized data were statistically evaluated separately.

Absolute quantification

Absolute quantification of MMP9, A1AT, APOB, A2MG, and elastin was determined by calculating the ratio between the total peak areas of the standard peptide and the corresponding labeled peptide (the internal standard). The concentration of the peptide in the sample was calculated by dividing the ratio by the slope of the calibration curve.

Statistical analysis of individual samples

Protein level was determined in samples obtained from welding fume exposure and from filtered air exposure. To assess the changes at the proteome level, the statistical analyses were conducted on the peptide level and further summarized on the protein level. Data were analyzed with SPSS for Windows 23.0 (SPSS Inc., Chicago, IL, USA). Statistical significance was defined as p < 0.05. The data were statistically analyzed with linear mixed modeling (LMM) and Wilcoxon signed rank test.

Linear mixed modeling

To assess if the changes in protein levels were statistically significant, the data were analyzed with LMM, a method in which several variables are considered simultaneously, e.g. the time course of the study and the involvement of multiple experimental exposure factors [26].

The dependent variables in LMM were the natural logarithmically transformed total peak area normalized to the isotopically labeled reference peptide of MMP9, A1AT, APOB, A2MG, respectively, for the relative quantification. Fixed effects in the model were the two different exposure types (welding fumes and filtered air) and each time point of sampling (immediately after exposure and the day after exposure). Subject was set to a random effect. The sample taken before exposure was set as a covariate. The repeated covariance types chosen were the unstructured covariance structure, the first-order factor analytic covariance structure, and the compound symmetry covariance structure. The assumption of normally distributed residuals was statistically tested with the Shapiro–Wilk test. The normal distribution of the residuals was used to evaluate the validity of the assumptions of the statistical model and was used as a tool for model selection in combination with Akaike information criterion, which is a measure of the relative quality of a statistical model for a given set of data. Data were analyzed according to following model.

$$Log\,\left( {peptide} \right) = Over\,all\,mean + Exposure + Time$$

Peptide changes were shown by the estimated marginal means of the main effect between the exposure term.

The residual was not normally distributed for all peptides for the relative quantification, and the dataset was therefore further analyzed with an additional non-parametric test to additionally check if the significantly changed proteins could be detected with a less sensitive model.

The concentrations obtained by absolute quantification were natural logarithmically transformed and used as the dependent variable. The data were normally distributed, therefore no more statistical analyses were used for the absolute quantification of the data.

Wilcoxon signed rank test

The Wilcoxon signed rank test was conducted for the relative quantification of the different peptides between the different exposure groups at the different time points, calculating the difference between NL1 and NL2 and between NL1 and NL3 for each peptide for each group.

False discovery rate (FDR) analysis

Control for multiple testing (8 tests × 276 peptides = 2208 tests for the LMM main effect test, 6 tests × 2 time points × 276 peptides = 3312 tests for the Wilcoxon signed rank test, and FDR-level: 0.05) was adjusted for with sequential goodness of fit for multiple testing using the SGoF+ algorithm [27,28,29], which resulted in an adjusted significance level 0.030 for the LMM test and 0.036 for the Wilcoxon signed rank test. The multiple test adjustment was conducted with the software SGoF+ (20) (downloaded from http://webs.uvigo.es/acraaj/SGoF.htm) and with the user-defined significance level set at p < 0.05.

Spearman rank correlation

Associations between MMP9 concentration and FEV1% and FVC% were examined using Spearman’s rank correlation. Statistical significance was defined as p < 0.05 (two-sided p values were determined).

IPA analysis

To evaluate common interactions and pathways of the identified proteins with changed levels, the results were further analyzed with IPA software (Ingenuity Systems, Redwood City, CA, USA, www.ingenuity.com). Significantly changed proteins were mapped and compared to known pathways, diseases, functions, and connecting regulators. The data were also evaluated based on upstream regulators with positive z-scores (activating capacity) and negative z-scores (inhibiting capacity). Key regulators are proteins with an earlier established association with the activation or inhibition of many of the identified proteins, but the regulators are not necessarily detected themselves in the analysis. Default settings were used except for species, which was set to human. In addition, only experimentally observed relationships were considered.

Results

Shotgun data analysis of pooled samples

Protein identification and relative quantification

The six pooled samples from the welders were analyzed by shotgun LC–MS/MS. In total, 336 proteins were identified in the pooled samples, whereof 137 proteins remained when the criteria described in the experimental section were applied (Additional file 1). The first observations (NL2/NL1) and second observations (NL3/NL1) were further analyzed with IPA for selection of candidate proteins to be included with targeted SRM (Additional files 2, 3, 4).

Pathway analysis of proteins identified with shotgun analysis and selection of proteins to be further analyzed with SRM

The pathway analysis indicated several pathways associated with the exposure. These were activation of the LXR/RXR pathway in monocytes, the inflammatory defense mechanism in macrophages, the coagulation system, and anti-inflammatory mechanisms, all of which influence the status of the extracellular matrix (ECM) (Additional files 2, 3). From these pathways, five representative proteins were selected for comprehensive quantitative data collection, including MMP9, which was identified as a key connecting mediator (Additional file 3); A1AT and A2MG, two anti-inflammatory proteins associated with MMP9 activity; APOB, an important lipid-loading protein involved in the low density lipoprotein protein complex; and elastin, an important extracellular protein involved in sustaining the elasticity of the ECM.

Proteins with currently unknown biological functions (according to IPA) were excluded from further analysis. In total, 97 proteins remained in the candidate list to be further analyzed with the targeted SRM method on the individual level. To enable a comprehensive monitoring of biological processes, proteins theoretically involved in the above-mentioned pathways were added to the candidate list. Finally, 150 protein candidates and 4 in vivo degradation peptides from collagen 4A1 and elastin remained in the list (Additional file 4). SRM assays were developed for these candidates (Additional files 5, 6).

In order to further investigate the findings from the analysis of pooled samples, individual samples were analyzed with SRM.

SRM analysis of individual samples

Statistical evaluation

Relative quantification of individual samples using LMM

Forty-six significantly changed proteins were found to be associated with welding fume particle exposure compared to the filtered air (Table 1). The beta value reflects the average effect change constant across the individuals between the exposures to welding fumes and filtered air. After the multiple statistical tests were controlled for FDR, 32 proteins were still significantly changed.

Table 1 Significantly changed proteins from the individual samples evaluated by LMM
Full size table

Relative quantification and statistical analysis with the Wilcoxon signed rank test

In total, 56 proteins were significantly changed when using the Wilcoxon signed rank test to compare the changes of peptide levels between the groups at the different time points (Table 2). After adjustment for FDR, 35 proteins still remained significantly changed.

Table 2 Significantly changed proteins from the individual samples evaluated with Wilcoxon signed rank test
Full size table

Absolute quantification of individual samples using LMM

Quantitative data could only be obtained for MMP9, A1AT, and A2MG (Additional file 9). LODs were determined for the proteins in samples without nasal lavage matrix (Table 3).

Table 3 Absolute quantification of individual sample
Full size table

The absolute quantification of the individual samples showed that 90% of all samples had a concentration higher than the LOD for A1AT, 95% of all samples had a concentration higher than the LOD for A2MG, and 88% of all samples had a concentration higher than the LOD for MMP9 (Table 3). The anti-inflammatory proteins A1AT and A2MG were found in higher amounts than MMP9 in all samples (Table 3). The only protein that was significantly changed due to welding fume exposure compared to the filtered air was MMP9. MMP9 had a beta value of 0.77, indicating a decrease in the protein level (Table 3).

MMP9 levels and lung function

The relationships between MMP9 concentration and lung function were evaluated. Combining the data of the welding fume particle and filtered air exposure for MMP9 concentration and the lung function, revealed that the levels of MMP9 were significantly associated with FVC (p = 0.014) and FEV1 (p = 0.014) before exposure and with FVC (p = 0.015) and FEV1 (p = 0.001) after exposure (Fig. 3). MMP9 had no significant association with FVC or FEV1 when separating the welding fume particle and filtered air exposure data in the statistical test for the data assessed before exposure. The levels of MMP9 in nasal lavage samples collected after exposure to welding fume particles were significantly associated with FEV1 (p = 0.005), and a trend toward an association with FVC (p = 0.079) was seen. The welders with the lowest FVC and FEV1 had low levels of MMP9 compared to the rest of the welders in this study (Additional file 10). The welders with the highest FVC and FEV1 during the baseline medical examination had high levels of MMP9 compared to the rest of the welders in this study.

Fig. 3
figure 3

Correlation between the MMP9 concentration (fmol/µL) and lung function (FVC and FEV1%) aCombining the data from both the welding fume exposure day and filtered air exposure day. *Correlation is significant at the 0.05 level (2-tailed). **Correlation is significant at the 0.01 level (2-tailed)

Full size image

Quality of the SRM data

Duplicate analyses of 126 prepared samples were performed by LC-SRM. Based on analyses of all samples, 69% of the SRM assays had a CV less than 20%. When normalizing the total area of each SRM assay to the total area of either of the isotopically labeled peptides belonging to MMP9, A1AT, A2MG, and APOB, 84, 83, 72, 73, and 84%, respectively, had a CV less than 20%.

Pathway analysis of significantly changed proteins from individual samples

Pathway analysis of the proteins with significantly changed levels according to LMM suggested that the proteins had two major pathways in common (Table 4). The first pathway was activation of acute phase response signaling, and this was supported by protein data indicative of induced inflammation, including acute phase proteins such as hemopexin and serum amyloid A1, both of which were upregulated following exposure to welding fumes. The second pathway was activation of LXR/RXR, a nuclear receptor family involved in lipid, inflammation, and cholesterol regulation, and this was supported by decreased levels of MMP9 and by increased levels of paraxonase 1 and clusterin, both of which are proteins found in high-density lipoproteins.

Table 4 Pathways associated with welding fume exposure
Full size table

Based on the significantly changed proteins, a number of key regulators were identified in the IPA analysis (Fig. 4). A key connecting regulator identified in this study was transforming growth factor (TGF)-β, and it was indicated to be inhibited. The following mediators that were detected in our data set supported this finding: MMP9, tissue inhibitor of metalloproteinases 1, fibronectin, and vimentin. Another identified key regulator was IL-6, a pro-inflammatory cytokine, and this was supported by MMP8, serum amyloid A1, interleukin-1 receptor antagonist, and tissue inhibitor of metalloproteinases 1. IL-6 was also associated with TGF-β. TNF and IL-1β were also identified as key regulators in the induced proteome effect. The key regulators and proteins identified in this study have previously been associated with several diseases, including chronic obstructive pulmonary disease, cystic fibrosis, and asthma (Fig. 4).

Fig. 4
figure 4

Network of the significantly changed proteins and the upstream regulators. The significantly changed proteins from the LMM table were associated with the upstream regulators IL6, estrogen receptor, H1FA, IFNG, IL1B, NEUROG1, SPP1, SYVN1, TGFB1, TNF, and STAT3. The canonical pathways connected to these molecular regulators included glucocorticoid receptor signaling. Biomarker enrichment of this network connection showed connections to cystic fibrosis, chronic obstructive pulmonary disease, and asthma. Red = increased measurement, Green = decreased measurement, Orange = indicated activated, Blue = indicated inhibited

Full size image

The significantly changed proteins identified with the Wilcoxon signed rank test also showed activation of the LXR/RXR pathway and activation of acute phase signaling (Additional file 11).

Discussion

Biological findings

In this study, we observed the effects on the proteome level in the nasal region after exposure to welding fumes. Proteins with changed levels were detected following welding fume particle exposure compared to filtered air, although neither clinical symptoms nor specific pro-inflammatory markers in nasal lavage were recorded in the study subjects [18]. This suggests that by applying a more comprehensive proteomic analysis, important biological information might still be gained even if no clinical symptoms can be recorded. Overall, the data analysis suggests that the proteins that differed in expression between welding fume particle and filtered air exposure are primarily involved in acute phase response and lipid regulation.

Changes in the proteome might result from a large number of factors, including metabolism and current health status, and the levels of many proteins with regulatory functions show a diurnal variation [30, 31]. By applying a study design where each study subject is its own control and where protein changes are associated with welding fume particle exposure compared to filtered air exposure, any changed levels of proteins will most likely be linked to the external exposure. Diurnal variation of protein levels was thus automatically taken into account. It is possible that previous exposure might influence the symptoms and sub-clinical responses, but in this project exposures were performed on Mondays making sure that no occupational exposure to welding occurred at least 48 h prior to the experimental exposure, minimizing the short term response interference.

In epidemiological studies of welders, the study subjects are exposed to both the emitted particles and the gases making it impossible to determine the exact cause of any biological response. In this study, the gases and particles were separated, and the study subjects were mostly exposed to the particle fraction. Previous studies have shown that when welding fumes are inhaled by humans, an inflammatory response is induced [15]. The result from the present study shows that the welding fume particles alone might induce an inflammatory response on the proteome level.

The subjects were exposed to particles at levels corresponding to a fifth of the permissible occupational exposure limit level in Sweden. The inhaled dose was, however, probably even less during the experimental exposure than it would have been in a true exposure situation because the welders were sitting still in a chamber without any physically demanding tasks. Even with this limitation, the welding fume particle exposure could still induce a protein response in the nasal region of the symptomatic welders.

The pro-inflammatory cytokines TNF, IL-1β, and IL-6 are suggested in this study to be key regulators connecting the proteins with changed levels, and this supports previous findings regarding inflammation as an important mechanism behind the response induced by welding fume exposure [5, 32]. Our data also indicate that the welding fume particle exposure not only affects the inflammatory response, but also induces imbalances in the regulation of the ECM. There a number of structural proteins that form the basis of connective tissue in the respiratory tract, including collagen I A2, desmoplakin, fibronectin, moesin, tight junction protein ZO-1, and junctional adhesion molecule A. The formation and degradation of connective tissue is highly regulated by the enzymatic activity of proteases. Significantly changed levels of several structural proteins were found following welding fume particle exposure, many with increased levels the day after exposure. Also, significantly changed levels of proteases such as metalloproteinases and a serine protease HTRA1, were detected. MMP8 increased in nasal lavage after exposure, and MMP9 decreased after exposure. MMP8 and MMP9 belong to a family of proteases that degrade ECM proteins. They play a central role in normal tissue structure to maintain the balance between the formation and degradation of ECM proteins and are involved in respiratory tract remodeling [33, 34]. MMPs are in turn inhibited by tissue inhibitor of metalloproteinases [35]. TGF-β was identified in the IPA analysis (Fig. 4), and this cytokine activates fibroblasts to produce collagen [36] and thus plays an important role in tissue repair by regulating the activities of proteins involved in the structural formation of connecting tissue as well as inflammatory responses.

The increased or decreased levels of structural proteins in nasal lavage fluid might be the result of tissue damage induced by the exposure, such that damaged cells release their intracellular content; inflammatory or epithelial cells secreting extracellular components and cleaving the structural proteins; or damage to the tight junctions of the epithelial cells such that the permeability increases.

Previous studies show an association between welding fume exposure and increased risk of developing fibrosis, chronic obstructive pulmonary disease, and asthma [37,38,39], and MMPs have been suggested to play an important role in the pathogenesis of these respiratory diseases (Fig. 4). MMPs can show both pro- and anti-fibrotic activity by degrading ECM proteins and activating various signaling molecules that influence the inflammatory response as well as through inhibition of ECM degradation processes [40,41,42]. Chronic obstructive pulmonary disease and asthma are both associated with airflow restriction and progressive remodeling, which affect the respiratory tract, although the localization and character of the structural changes throughout the respiratory tract wall are different [43]. The remodeling concerns highly composed structural transformation, which affects the airways and leads to significant functional impairment [44].

In normal tissue, the balance between the formation and degradation of ECM proteins is strictly controlled to maintain tissue structure and function [45]. If this balance is disrupted, accumulation of structural connective tissue can lead to fibrotic tissues. Defined by the accumulation of ECM proteins, fibrosis results in scarring and thickening of the affected tissue; it is in essence an exaggerated wound healing response that interferes with normal organ function. In this study, a significant correlation between the concentration of MMP9 in nasal lavage and lung function (FEV1 and FVC%) was detected (Fig. 3). A relationship between MMP9 concentration and pulmonary function, bronchodilator response, and emphysema severity has previously been suggested [46, 47]. Together, these results suggest that the MMP9 level might be associated with the lung function status of each individual welder in this study. In a previous study of welders with and without symptoms, gene expression analysis in induced sputum [48] revealed that genes differentially expressed following welding fume exposure were connected to responses related to tissue damage, immune defense, and inflammation, and one gene that was highly expressed by the asthmatic subjects codes for MMP25, a membrane bound metalloproteinase. Prolonged or frequent exposure to pro-inflammatory agents might promote the development of chronic inflammation with associated damage and transformation of the respiratory epithelium [49, 50].

Also interesting is the presence and increase of the high density lipoprotein-associated proteins serum amyloid A1, paraxonase 1, and clusterin. Knowledge is limited regarding the presence of high density lipoprotein in nasal mucosa, but in the blood stream serum amyloid A1 is an important acute phase protein, and paraxonase 1 is an important antioxidant [51, 52]. Increased paraxonase 1 levels after welding fume exposure might reflect increased oxidative stress in the mucosa and thereby an increased need for antioxidant activities. Regarding clusterin, it has previously been shown that the human clusterin gene is responsive to acute in vivo oxidative stress induced by heavy metals [53]. However, these suggestions have to be further studied because there is a possibility of leakage from the circulation as result of exposure-related tissue damage. A previous study [18] of a symptomatic welding group showed that they had an increasing trend of blood neutrophils after exposure and that they had a tenfold higher level of exhaled breath condensate leukotriene B4 compared to the other groups. Leukotrienes use lipid signaling to convey information order to regulate immune responses, and leukotriene B4 is involved in inflammation [54]. Its primary function is to recruit neutrophils to areas of tissue damage, though it also helps to promote the production of inflammatory cytokines by various immune cells.

Previous studies have shown that most of the upregulated inflammatory response was diminished within 19 h after exposure to welding fumes [55]. This study showed that the inflammatory response that results from welding fume particle exposure is still detectable the day after exposure and that the levels of proteins in the ECM were affected. In parallel, immunosuppressive mediators (A2MG, interleukin-1 receptor antagonist, and uteroglobin) were also detected the day after exposure suggesting that processes attempting to resolve the inflammation had been initiated. If the inflammation is still present the day after exposure and the exposure continues, this might lead to a prolonged inflammatory response with no time to resolve itself.

Methodological considerations

Discovery followed by a targeted mass spectrometry method

More than 300 proteins were identified with the shotgun analysis. Cutoff criteria were used to minimize labor and the time-consuming analysis of redundant proteins. To focus on the proteins that are changed on global level, the median ratio of NL2/NL1 and NL3/NL1 for the increased or decreased proteins was chosen as the cutoff for proteins considered to have a relatively large enough change to detect. The chosen cutoffs were based on the magnitude of change detected for this set of data, which could differ depending on the exposure and the studied group. The combination of the discovery shotgun method and targeted mass spectrometry method (Fig. 2) enabled us to focus the measurements on relevant proteins with potential interest in the induced effects. Another advantage was that crucial proteins that for some reason were not detected with the shotgun method were identified via the pathway analysis in the targeted method and were selected for further analysis.

Absolute quantification in comparison to relative quantification for individual samples

Statistical analysis of protein levels determined with absolute quantification showed that the level of MMP9 was significantly changed after welding exposure, whereas the statistical analysis of protein levels determined with relative quantification showed that MMP9 as well as A2MG were significantly changed. However, absolute quantification was based on only one peptide, which made this method less robust, although the analysis had high reproducibility. These results show that it is important to include several peptides in each assay for quantification. In shotgun MS identification of a protein is usually based on many peptides. For relative quantitative SRM methods as well, many SRM assays with peptide criteria are included in the measurement of a protein giving a high degree of protein coverage. But since mass spectrometry is not inherently quantitative because of differences in the ionization efficiency and/or detectability of the many peptides in a given sample it cannot determine the absolute abundance of proteins in a sample. The intensity of a peak in a mass spectrum is not a good indicator of the amount of the analyte in the sample, although differences in peak intensity of the same analyte between multiple samples accurately reflect relative differences in its abundance. In this study the absolute quantification was based on peptides as opposed to proteins consequently subjected to trypsin digestion and therefore post-digest addition of internal standards was applied. One limitation with such a strategy is that any loss occurring at the digestion step is not compensated for. It has earlier been shown that losses of peptides may occur during trypsin digestion [56]. For absolute quantification, the selected peptide is also chosen based on several stringent criteria, thus leaving very few peptides selectable for absolute quantification but there is no guarantee that this peptide ionizes well. This might leave us with a peptide with low reproducibility and low sensitivity which will have a major effect on the absolute quantification method. Therefore, the label free relative quantification is advantageous in terms of providing reasonable proteome coverage with only little investment and method optimization, but may be limited in terms of accuracy and selectivity.

Verification of protein changes

Different LC–MS techniques were applied in this study to qualitatively and quantitatively determine proteins. A discovery based shotgun LC–MS analysis was used initially to identify proteins on the group level. Then, semi-quantitative and quantitative methods were applied on individual samples to determine and quantify protein levels. The mass spectrometry methodology and especially in combination with LC-separation offers identification and quantification of very high specificity and accuracy. In the SRM assays we included unique peptides, at least three transitions per peptide and at least three assays per protein were used. In addition, an isotope labelled internal standard was included during absolute quantification to further support the specificity. If robust mass spectrometric assays for quantification are available, e.g. as for MMP9 in this study, such methodology may be sufficient to verify proteome changes determined with semi-quantitative mass spectrometry. Should robust quantitative MS assays not be available then intact protein testing using immunoassays may be more suitable to confirm results. However, immunoassay of proteins may suffer from cross-reactivity generating overestimated results and it could be argued that SRM LC–MS is superior to western blot when it comes to accuracy and specificity [57].

Global normalization

The strength of this approach was the use of the isotopically labeled peptides as global normalizers for the peptides that were relatively quantified. This was applied to control for variations in sample preparation, desalting, sample handling, and instrument operation. This type of normalization is rather rough, and there is a risk of losing important findings for peptides that do not vary in the same way as the global isotopically labeled peptide. Using more than one isotopically labeled peptide as a global normalizer increases the chance of finding a good normalization for each peptide.

Statistical evaluation on the peptide level for individual samples

For most of the significantly changed proteins, only one peptide showed a significantly changed level (Tables 1, 3). The reason for this could be that different signal intensities could be generated from different peptides depending on the ionization capacity of that peptide. A weak signal might cause greater variation and therefore a larger degree of uncertainty. Proteins were represented by 1–4 peptides, and when the trends in the changes in peptide levels representing the same protein were in line with each other, which was the case for a majority of the analyzed peptides, the results can be considered reliable.

Thirty proteins could be identified both with LMM and Wilcoxon signed rank test, and 26 proteins were exclusively identified with the Wilcoxon signed rank (Table 3) and 16 proteins were exclusively identified with LMM. To a large extent, the two statistical models identified similar proteins that were significantly changed upon exposure to welding fumes. The LMM method yielded a higher sensitivity (lower p values) for proteins that were significantly changed. LMM can handle several variables in the same statistical analysis, resulting in the need for fewer statistical tests, and thus lower FDR.

Pathway analysis

Pathway analysis identifies common signaling molecules shared between the proteins and if several proteins can be identified in the same pathway, then there is a higher likelihood that this pathway is involved in the biological response. The changed direction of the proteins for the identified canonical pathways to be activated are expected to go in e certain direction when its activated. Our data does not fully go in line with this canonical pathway protein direction. The changed protein levels after an exposure might be driven by a cascade of changed proteins [58, 59], with different protein changes appearing at different time points. In this study, the sampling is taken just a few times after the exposure. This will make the appearance of all the altered proteins impossible to find from the few sampling points.

Sampling

By instilling the nasal cavity with saline solution, cells and secreted mediators as well as intracellular components released because of cell damage in the nose can be collected. If cells were damaged by the sampling procedure and not by the exposure, this would occur at all sampling times and thus no significant change would be observed between the two exposures. However, if an exposure induces some cell damage so that, for example, the permeability of the nasal mucosa were to increase, this could be observed by detecting components from the mucosal tissue.

Conclusion

In conclusion, by applying a comprehensive mass spectrometric work flow a proteome-level response in the nasal upper airway region was demonstrated following mild steel welding fume particle exposure even though no clinical symptoms could be detected. Proteins with significantly changed levels were associated with inflammation and lipid signaling pathways and with changes in the status of the ECM. Several proteins, including MMP8 and MMP9, that were associated with the welding fume particle exposure should be further investigated as biomarkers in future studies.

Abbreviations

ACN:

acetonitrile

MS:

mass spectrometry

FA:

formic acid

LOD:

limit of detection

FDR:

false discovery rate

SRM:

selected reaction monitoring

IL:

interleukin

LC:

liquid chromatography

IPA:

ingenuity pathway analysis

BCA:

bicinchoninic acid

DTT:

dithiothreitol

FVC:

forced vital capacity

FEV:

forced expiratory volume

PD15 :

provocative dose to produce a 15% fall in FEV1

MCH:

methacholine

MMP9:

metalloproteinase 9

A1AT:

alpha-1-antitrypsine

A2MG:

alpha-2-macrogloboline

APOB:

apolipoprotein-B

CV:

coefficient of variation

LMM:

linear mixed model

ECM:

extracellular matrix

LXR/RXR:

liver X receptor/retinoid X receptor

References

  1. WHO. Hazard prevention and control in the work environment: airborne dust. 1999.

  2. OSHA. Occupational safety and health guidance manual for hazardous waste site activities. 1986.

  3. NIOSH. NIOSH pocket guide to chemical hazards. 2007.

  4. WHO. WHO guidelines for indoor air quality: selected pollutants. Copenhagen: WHO Regional office for Europe. 2010.

  5. Antonini JM. Health effects of welding. Crit Rev Toxicol. 2003;33:61–103.

    Article  CAS  PubMed  Google Scholar 

  6. Guha N, Loomis D, Guyton KZ, Grosse Y, El Ghissassi F, Bouvard V, Benbrahim-Tallaa L, Vilahur N, Muller K, Straif K. Carcinogenicity of welding, molybdenum trioxide, and indium tin oxide. Lancet Oncol. 2017;18:581–2.

    Article  CAS  PubMed  Google Scholar 

  7. Berlinger B, Benker N, Weinbruch S, L’Vov B, Ebert M, Koch W, Ellingsen DG, Thomassen Y. Physicochemical characterisation of different welding aerosols. Anal Bioanal Chem. 2011;399:1773–80.

    Article  CAS  PubMed  Google Scholar 

  8. Jenkins NT, Eagar TW. Chemical analysis of welding fume particles. Weld J New York. 2005;84(6):87.

    Google Scholar 

  9. Isaxon C, Pagels J, Gudmundsson A, Asbach C, John A, Kuhlbusch T, Karlsson J, Kammer R, Tinnerberg H, Nielsen J. Characteristics of welding fume aerosol investigated in three Swedish workshops. In: Journal of physics: conference series. IOP Publishing; 2009. p. 012059.

  10. Ali N, Mattsson K, Rissler J, Karlsson HM, Svensson CR, Gudmundsson A, Lindh CH, Jonsson BA, Cedervall T, Karedal M. Analysis of nanoparticle-protein coronas formed in vitro between nanosized welding particles and nasal lavage proteins. Nanotoxicology. 2016;10:226–34.

    Article  CAS  PubMed  Google Scholar 

  11. Hedmer M, Karlsson JE, Andersson U, Jacobsson H, Nielsen J, Tinnerberg H. Exposure to respirable dust and manganese and prevalence of airways symptoms, among Swedish mild steel welders in the manufacturing industry. Int Arch Occup Environ Health. 2014;87:623–34.

    Article  CAS  PubMed  Google Scholar 

  12. Jonsson LS, Tinnerberg H, Jacobsson H, Andersson U, Axmon A, Nielsen J. The ordinary work environment increases symptoms from eyes and airways in mild steel welders. Int Arch Occup Environ Health. 2015;88:1131–40.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Li H, Hedmer M, Wojdacz T, Hossain MB, Lindh CH, Tinnerberg H, Albin M, Broberg K. Oxidative stress, telomere shortening, and DNA methylation in relation to low-to-moderate occupational exposure to welding fumes. Environ Mol Mutagen. 2015;56:684–93.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Antonini JM, Krishna Murthy GG, Brain JD. Responses to welding fumes: lung injury, inflammation, and the release of tumor necrosis factor-alpha and interleukin-1 beta. Exp Lung Res. 1997;23:205–27.

    Article  CAS  PubMed  Google Scholar 

  15. Zeidler-Erdely PC, Erdely A, Antonini JM. Immunotoxicology of arc welding fume: worker and experimental animal studies. J Immunotoxicol. 2012;9:411–25.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Erdely A, Antonini JM, Young S-H, Kashon ML, Gu JK, Hulderman T, Salmen R, Meighan T, Roberts JR, Zeidler-Erdely PC. Oxidative stress and reduced responsiveness of challenged circulating leukocytes following pulmonary instillation of metal-rich particulate matter in rats. Part Fibre Toxicol. 2014;11:34.

    Article  PubMed  PubMed Central  Google Scholar 

  17. Popstojanov R, Antonini JM, Salmen R, Ye M, Zheng W, Castranova V, Fekedulegn DB, Kan H. Alterations in cardiomyocyte function after pulmonary treatment with stainless steel welding fume in rats. J Toxicol Environ Health A. 2014;77:705–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Dierschke K, Isaxon C, Andersson UB, Assarsson E, Axmon A, Stockfelt L, Gudmundsson A, Jönsson BA, Kåredal M, Löndahl J, Pagels J. Acute respiratory effects and biomarkers of inflammation due to welding-derived nanoparticle aggregates. Int Arch Occup Environ Health. 2017;90:451–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Sjöström M, Ossola R, Breslin T, Rinner O, Malmström L, Schmidt A, Aebersold R, Malmström J, Niméus E. A combined shotgun and targeted mass spectrometry strategy for breast cancer biomarker discovery. J Proteome Res. 2015;14:2807–18.

    Article  PubMed  Google Scholar 

  20. Schmidt A, Forne I, Imhof A. Bioinformatic analysis of proteomics data. BMC Syst Biol. 2014;8:S3.

    Article  PubMed  PubMed Central  Google Scholar 

  21. Isaxon C, Dierschke K, Pagels J, Löndahl J, Gudmundsson A, Hagerman I, Berglund M, Wierzbicka A, Assarsson E, Andersson U. A novel system for source characterization and controlled human exposure to nanoparticle aggregates generated during gas–metal arc welding. Aerosol Sci Technol. 2013;47:52–9.

    Article  CAS  Google Scholar 

  22. Mortstedt H, Karedal MH, Jonsson BA, Lindh CH. Screening method using selected reaction monitoring for targeted proteomics studies of nasal lavage fluid. J Proteome Res. 2013;12:234–47.

    Article  PubMed  Google Scholar 

  23. Veidal SS, Karsdal MA, Leeming DJ, Barascuk N, Skjøt-Arkil H, Vassiliadis E. Pathology biomarker assay. In Book Pathology biomarker assay City: Google Patents; 2014.

  24. Veidal SS, Karsdal MA, Nawrocki A, Larsen MR, Dai Y, Zheng Q, Hägglund P, Vainer B, Skjøt-Arkil H, Leeming DJ. Assessment of proteolytic degradation of the basement membrane: a fragment of type IV collagen as a biochemical marker for liver fibrosis. Fibrogenesis Tissue Rep. 2011;4:22.

    Article  CAS  Google Scholar 

  25. Eidhammer I, Barsnes H, Eide GE, Martens L. Computational and statistical methods for protein quantification by mass spectrometry. London: Wiley; 2012. p. 2018–234.

    Google Scholar 

  26. Clough T, Thaminy S, Ragg S, Aebersold R, Vitek O. Statistical protein quantification and significance analysis in label-free LC–MS experiments with complex designs. BMC Bioinform. 2012;13(Suppl 16):S6.

    Article  CAS  Google Scholar 

  27. Carvajal-Rodriguez A, de Una-Alvarez J, Rolan-Alvarez E. A new multitest correction (SGoF) that increases its statistical power when increasing the number of tests. BMC Bioinform. 2009;10:209.

    Article  Google Scholar 

  28. de Una-Alvarez J. The beta-binomial SGoF method for multiple dependent tests. Stat Appl Genet Mol Biol. 2012; 11:Article 14.

  29. Castro-Conde I, Dohler S, de Una-Alvarez J. An extended sequential goodness-of-fit multiple testing method for discrete data. Stat Methods Med Res. 2015;26(5):2356–75.

    Article  PubMed  Google Scholar 

  30. Boukelia B, Fogarty MC, Davison RCR, Florida-James GD. Diurnal physiological and immunological responses to a 10-km run in highly trained athletes in an environmentally controlled condition of 6 °C. Eu J Appl Physiol. 2017;117:1–6.

    Article  Google Scholar 

  31. Gooley JJ, Chua EC-P. Diurnal regulation of lipid metabolism and applications of circadian lipidomics. J Genet Genomics. 2014;41:231–50.

    Article  CAS  PubMed  Google Scholar 

  32. Jarvela M, Kauppi P, Tuomi T, Luukkonen R, Lindholm H, Nieminen R, Moilanen E, Hannu T. Inflammatory response to acute exposure to welding fumes during the working day. Int J Occup Med Environ Health. 2013;26:220–9.

    Article  PubMed  Google Scholar 

  33. Fanjul-Fernández M, Folgueras AR, Cabrera S, López-Otín C. Matrix metalloproteinases: evolution, gene regulation and functional analysis in mouse models. Biochim Biophys Acta. 2010;1803:3–19.

    Article  PubMed  Google Scholar 

  34. Page-McCaw A, Ewald AJ, Werb Z. Matrix metalloproteinases and the regulation of tissue remodelling. Nat Rev Mol Cell Biol. 2007;8:221–33.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Brew K, Nagase H. The tissue inhibitors of metalloproteinases (TIMPs): an ancient family with structural and functional diversity. Biochim Biophys Acta. 2010;1803:55–71.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Meng X-M, Nikolic-Paterson DJ, Lan HY. TGF-[beta]: the master regulator of fibrosis. Nat Rev Nephrol. 2016;12:325–38.

    Article  CAS  PubMed  Google Scholar 

  37. Meldrum M, Rawbone R, Curran AD, Fishwick D. The role of occupation in the development of chronic obstructive pulmonary disease (COPD). Occup Environ Med. 2005;62:212–4.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  38. Koh DH, Kim JI, Kim KH, Yoo SW. Welding fume exposure and chronic obstructive pulmonary disease in welders. Occup Med. 2015;65:72–7.

    Article  Google Scholar 

  39. Hannu T, Piipari R, Tuppurainen M, Nordman H, Tuomi T. Occupational asthma caused by stainless steel welding fumes: a clinical study. Eur Respir J. 2007;29:85–90.

    Article  CAS  PubMed  Google Scholar 

  40. Bonnans C, Chou J, Werb Z. Remodelling the extracellular matrix in development and disease. Nat Rev Mol Cell Biol. 2014;15:786–801.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Van Lint P, Libert C. Chemokine and cytokine processing by matrix metalloproteinases and its effect on leukocyte migration and inflammation. J Leukoc Biol. 2007;82:1375–81.

    Article  PubMed  Google Scholar 

  42. Giannandrea M, Parks WC. Diverse functions of matrix metalloproteinases during fibrosis. Dis Models Mech. 2014;7:193–203.

    Article  CAS  Google Scholar 

  43. Grzela K, Litwiniuk M, Zagorska W, Grzela T. Airway remodeling in chronic obstructive pulmonary disease and asthma: the role of matrix metalloproteinase-9. Arch Immunol Ther Exp (Warsz). 2016;64:47–55.

    Article  CAS  Google Scholar 

  44. James AL, Elliot JG, Abramson MJ, Walters EH. Time to death, airway wall inflammation and remodelling in fatal asthma. Eur Respir J. 2005;26:429–34.

    Article  CAS  PubMed  Google Scholar 

  45. Lu P, Takai K, Weaver VM, Werb Z. Extracellular matrix degradation and remodeling in development and disease. Cold Spring Harb Perspect Biol. 2011;3:a005058.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Koo H-K, Hong Y, Lim MN, Yim J-J, Kim WJ. Relationship between plasma matrix metalloproteinase levels, pulmonary function, bronchodilator response, and emphysema severity. Int J Obstruct Pulm Dis. 2016;11:1129–37.

    Article  CAS  Google Scholar 

  47. Atkinson JJ, Senior RM. Matrix metalloproteinase-9 in lung remodeling. Am J Respir Cell Mol Biol. 2003;28:12–24.

    Article  CAS  PubMed  Google Scholar 

  48. Jönsson LS, Nielsen J, Broberg K. Gene expression analysis in induced sputum from welders with and without airway-related symptoms. Int Arch Occup Environ Health. 2011;84:105–13.

    Article  PubMed  Google Scholar 

  49. Grzela K, Zagórska W, Jankowska-Steifer E, Grzela T. Review paper < br > chronic inflammation in the respiratory tract and ciliary dyskinesia. Cent Eur J Immunol. 2013;38:122–8.

    Article  CAS  Google Scholar 

  50. Haswell LE, Hewitt K, Thorne D, Richter A, Gaca MD. Cigarette smoke total particulate matter increases mucous secreting cell numbers in vitro: a potential model of goblet cell hyperplasia. Toxicol In Vitro. 2010;24:981–7.

    Article  CAS  PubMed  Google Scholar 

  51. Levels JH, Geurts P, Karlsson H, Marée R, Ljunggren S, Fornander L, Wehenkel L, Lindahl M, Stroes ES, Kuivenhoven JA, Meijers JC. High-density lipoprotein proteome dynamics in human endotoxemia. Proteome Sci. 2011;9:34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  52. Karlsson H, Kontush A, James RW. Functionality of HDL: antioxidation and detoxifying effects. In: von Eckardstein A, Kardassis D, editors. High density lipoproteins: from biological understanding to clinical exploitation. Cham: Springer; 2015. p. 207–28.

    Google Scholar 

  53. Alexopoulos EC, Cominos X, Trougakos IP, Lourda M, Gonos ES, Makropoulos V. Biological monitoring of hexavalent chromium and serum levels of the senescence biomarker apolipoprotein J/Clusterin in welders. Bioinorg Chem Appl. 2008. https://doi.org/10.1155/2008/420578.

    PubMed  PubMed Central  Google Scholar 

  54. Fiancette R, Vincent-Fabert C, Guerin E, Trimoreau F, Denizot Y. Lipid mediators and human leukemic blasts. J Oncol. 2011;2011:7.

    Article  Google Scholar 

  55. Wang Z, Neuberg D, Su L, Kim JY, Chen J-C, Christiani DC. Prospective study of metal fume-induced responses of global gene expression profiling in whole blood. Inhal Toxicol. 2008;20:1233–44.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  56. van den Broek I, Smit NP, Romijn FP, van der Laarse A, Deelder AM, van der Burgt YE, Cobbaert CM. Evaluation of interspecimen trypsin digestion efficiency prior to multiple reaction monitoring-based absolute protein quantification with native protein calibrators. J Proteome Res. 2013;12:5760–74.

    Article  PubMed  Google Scholar 

  57. Aebersold R, Burlingame AL, Bradshaw RA. Western blots versus selected reaction monitoring assays: time to turn the tables? Mol Cell Proteomics. 2013;12:2381–2.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. Gabay C, Kushner I. Acute-phase proteins and other systemic responses to inflammation. N Engl J Med. 1999;340:448–54.

    Article  CAS  PubMed  Google Scholar 

  59. Kholodenko BN. Cell signalling dynamics in time and space. Nat Rev Mol Cell Biol. 2006;2006(7):165–76.

    Article  Google Scholar 

Download references

Authors’ contributions

NA: Protein work flow set up. SRM: method development, all proteomic data analysis, result interpretation, manuscript writing. SL: shotgun analysis, commented on the manuscript. HK: shotgun analysis, manuscript writing. AW: chamber exposure study design and logistics, took part in carrying out the exposures (responsible for HRV monitoring and recording). JP: aerosol knowledge, commented on the manuscript. CHI: chamber exposure study design, welding, measurements, and characterization of the welding fume particles. AG: chamber study design, aerosol knowledge, commented on the manuscript. JR: aerosol knowledge, total surface area calculations, commented on the manuscript. JN: Chamber exposure study design, results interpretation, commenting on the manuscript. CL: Quantitative method development, protein work flow set up. MK: Chamber exposure study design, results interpretation, manuscript writing. All authors read and approved the final manuscript.

Acknowledgements

We would like to thank the welders who participated in the study. Furthermore, we thank Eva Assarsson and Ulla BK Andersson for handling the obtained material.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

Not applicable.

Ethics approval and consent to participate

The study was approved by the Ethical Committee at the Medical Faculty at Lund University (H4 160/2007) and written informed consent was obtained from the patints.

Funding

The study was financed by the Swedish Research Council for Health, Working Life and Welfare (FAS 2006-0803), the Swedish Governmental Funding for Clinical Research, and the Skåne County Council’s research and development foundation (ALF) and was performed within the framework of the Competence Centre for Medicine and Technology for Working Life and Society (METALUND) of Lund University that is supported by FORTE.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

  1. Division of Occupational and Environmental Medicine, Lund University, Lund, Sweden

    Neserin Ali, Jörn Nielsen, Christian H. Lindh & Monica Kåredal

  2. Occupational and Environmental Medicine Center, Department of Clinical and Experimental Medicine, Linköping University, Linköping, Sweden

    Stefan Ljunggren & Helen M. Karlsson

  3. Department of Design Sciences, Ergonomics and Aerosol Technology, Lund University, Lund, Sweden

    Aneta Wierzbicka, Joakim Pagels, Christina Isaxon, Anders Gudmundsson & Jenny Rissler

Authors
  1. Neserin Ali
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Stefan Ljunggren
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Helen M. Karlsson
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Aneta Wierzbicka
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Joakim Pagels
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Christina Isaxon
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Anders Gudmundsson
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Jenny Rissler
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Jörn Nielsen
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Christian H. Lindh
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Monica Kåredal
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Neserin Ali.

Additional files

Additional file 1.

Protein identification and mechanism hypothesis generation from pooled samples. The protein list generated from the analysis of pooled samples. Calculated label free quantification ratios from welding fume exposure—after exposure/before exposure (NL2/NL1) and day after exposure/before exposure (NL3/NL1). The proteins in this list were detected in more than half of the samples and had a ratio more than 1.3 or less than 0.8. aNot detected in NL1. bNot detected in NL2 and/or NL3.

Additional file 2.

Canonical pathway results of pooled samples. IPA results of pathways found to be significantly induced when proteins that had a ratio ≥ 1.3 or ≤ 0.8 from additional file 1 were analyzed. After exposure = NL2/NL1. Day after exposure = NL3/NL1.

Additional file 3.

Upstream analysis for observation immediately after exposure and the day after exposure.

Additional file 4.

Flowchart for protein list generation and analyses. A total of 336 proteins were detected with the shotgun analysis. Cutoff criteria were used to minimize labor and the time-consuming analysis of redundant proteins. Mechanism hypotheses were generated with the help of pathway analysis and literature findings.

Additional file 5.

Proteins to be targeted with a SRM method. Proteins with SRM assays from a previously developed method and proteins to be added to the SRM method.

Additional file 6.

Final SRM method for this study. The final SRM method contained 256 precursors with 811 transitions targeted 130 proteins and 4 in vivo peptide degradation products.

Additional file 7.

SRM assay development for absolute quantification. Transitions and CE for the light and heavy peptides.

Additional file 8.

Calibration curves. Calibration curves prepares in nasal lavage matrix: Two sample preparations of the standard curve were made and run during the beginning and the end of the whole runs. They were each injected twice. Calibration curves prepared in MilliQ water matrix: One sample preparation of the standard curve was prepared in MilliQ water as a matrix instead of nasal lavage to confirm that the matrix of the sample was not affecting the quantification of the peptide.

Additional file 9.

Reproducibility and absolute quantification.

Additional file 10.

Medical examination of the welders. Medical examination of the welders with lower airway symptoms during the previous month, including skin prick test positivity for a standard panel of aeroallergens, methacholine test (MCH) positivity, and lung function (FVC and FEV1 as % of predicted) at the baseline examination before the study. MMP9 concentration (fmol/µL) and lung function (FVC and FEV1 as % of predicted) before and after exposure to welding fume particles (A) or filtered air (B). N/A: no data could be obtained due to missing samples.

Additional file 11.

The significantly changed proteins found with the Wilcoxon signed rank test showed that these canonical pathways were activated due to welding fume exposure. Positive z-score: the canonical pathway was predicted to be activated. Negative z-score: the canonical pathway was predicted to be deactivated.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ali, N., Ljunggren, S., Karlsson, H.M. et al. Comprehensive proteome analysis of nasal lavage samples after controlled exposure to welding nanoparticles shows an induced acute phase and a nuclear receptor, LXR/RXR, activation that influence the status of the extracellular matrix. Clin Proteom 15, 20 (2018). https://doi.org/10.1186/s12014-018-9196-y

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12014-018-9196-y

Keywords

Clinical Proteomics

ISSN: 1559-0275