Skip to content

Advertisement

  • Research
  • Open Access

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

  • 1Email authorView ORCID ID profile,
  • 2,
  • 2,
  • 3,
  • 3,
  • 3,
  • 3,
  • 3,
  • 1,
  • 1 and
  • 1
Clinical Proteomics201815:20

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

  • Received: 28 July 2017
  • Accepted: 2 May 2018
  • Published:

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.

Keywords

  • Chamber study
  • Welding fume particles
  • Nasal lavage
  • Effects
  • Proteomics
  • Mass spectrometry
  • Pathways

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 [14]. 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 [1517]. 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
Fig. 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

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
Fig. 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

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 [2729], 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

AC

Gene name

Protein names

Beta value

CI 95% (min–max)

FDR

Sig. peptides

Nr. peptides

Increased

P04040

CATA

Catalase

2.42

(1.46–4.02)

 

1

2

P02766

TTHY

Transthyretin

1.38

(1.02–1.86)

f

1

2

P15924

DESP

Desmoplakin

1.37

(1.05–1.77)

 

1

2

A6NC86

PINLY

Phospholipase A2 inhibitor

1.35

(1.09–1.66)

 

1

1

P80188

NGAL

Neutrophil gelatinase-associated lipocalin

1.35

(1.02–1.78)

f

2

3

P18510

IL1RN

Interleukin-1 receptor antagonist protein

1.30

(1.02–1.67)

f

1

1

Q8TDL5

LPLC1

BPI fold-containing family B member 1

1.29

(1.03–1.61)

f,f

2

2

P0DJI8

SAA1

Serum amyloid A-1 protein

1.28

(1.00–1.62)

f

1

2

P01023

A2MG

Alpha-2-macroglobulin

1.26

(1.10–1.45)

 

1

4

P02790

HEMO

Hemopexin

1.26

(1.05–1.51)

 

1

3

Q08380

LG3BP

Galectin-3-binding protein

1.26

(1.10–1.44)

 

1

2

P01834

IGKC

Immunoglobulin kappa constant

1.24

(1.04–1.48)

 

2

2

P05155

IC1

Plasma protease C1 inhibitor

1.23

(1.02–1.50)

f

1

3

P08603

CFAH

Complement factor H

1.23

(1.03–1.48)

 

1

1

P02760

AMBP

Protein AMBP

1.23

(1.04–1.46)

 

1

2

Q92743

HTRA1

Serine protease HTRA1

1.23

(1.10–1.37)

 

1

2

P22894

MMP8

Neutrophil collagenase

1.22

(1.08–1.38)

 

1

1

P07737

PROF1

Profilin-1

1.20

(1.00–1.43)

f

1

1

P27169

PON1

Serum paraoxonase

1.20

(1.11–1.29)

 

1

3

P10909

CLUS

Clusterin

1.19

(1.11–1.28)

 

2

3

Q96DA0

ZG16B

Zymogen granule protein 16 homolog B

1.18

(1.03–1.37)

 

1

1

P04792

HSPB1

Heat shock protein beta-1

1.16

(1.01–1.34)

f

1

2

Decreased

P07602

SAP

Prosaposin

0.92

(0.88–0.97)

 

1

3

Q9GZX6

IL22

Interleukin-22

0.88

(0.79–0.99)

f

1

3

P01034

CYTC

Cystatin-C

0.88

(0.79–0.97)

 

1

1

Q9Y624

JAM1

Junctional adhesion molecule A

0.87

(0.77–0.97)

 

1

2

Q13822

ENPP2

Ectonucleotide pyrophosphatase

0.85

(0.74–0.97)

 

1

1

O75223

GGCT

Gamma-glutamylcyclotransferase

0.82

(0.73–0.93)

 

1

3

P37802

TAGL2

Transgelin-2

0.82

(0.70–0.95)

 

1

1

Q14118

DAG1

Dystroglycan

0.81

(0.66–1.00)

f

1

2

P04080

CYTB

Cystatin-B

0.81

(0.71–0.92)

 

1

1

P02751

FINC

Fibronectin

0.80

(0.66–0.97)

 

1

3

P14780

MMP9

Matrix metalloproteinase-9

0.78

(0.62–0.97)

f

2

3

P31944

CASPE

Caspase-14

0.77

(0.62–0.95)

 

2

2

P05164

PERM

Myeloperoxidase

0.75

(0.62–0.92)

f

2

2

Q14019

COTL1

Coactosin-like protein

0.72

(0.58–0.89)

 

1

2

P06727

APOA4

Apolipoprotein A-IV

0.71

(0.54–0.93)

 

1

3

P01033

TIMP1

Metalloproteinase inhibitor 1

0.71

(0.57–0.88)

 

1

2

O75594

PGRP1

Peptidoglycan recognition protein 1

0.70

(0.54–0.90)

 

3

3

P29401

TKT

Transketolase

0.69

(0.50–0.94)

 

1

2

P04406

G3P

Glyceraldehyde-3-phosphate dehydrogenase

0.68

(0.57–0.81)

 

1

1

O75556

SG2A1

Mammaglobin-B

0.67

(0.47–0.96)

f

1

3

P02749

APOH

Beta-2-glycoprotein 1

0.65

(0.49–0.86)

f

2

2

P01008

ANT3

Antithrombin-III

0.53

(0.31–0.89)

 

1

1

Q6UWW0

LCN15

Lipocalin-15

0.43

(0.25–0.72)

 

1

2

P08670

VIME

Vimentin

0.72

(0.54–0.95)

 

1

2

Relative quantification of individual samples generated this protein list of the significantly changed proteins after welding fume exposure compared to filtered air, and these were evaluated by LMM, p < 0.05. The beta value represents the average change of the protein level after welding fume exposure

f Peptides not significantly changed after multiple test correction (FDR), p > 0.03. Sig. peptides number of significantly changed peptides. Nr. Peptides number of peptides used in the SRM method representing the specific protein

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

AC

Gene name

Protein name

After exposure

Day after exposure

Significant peptides

FDR

Proteins increased immediately after exposure

P02765

FETUA

Alpha-2-HS-glycoprotein

Increase

 

1

f

Proteins increased immediately after exposure and the day after exposure

A6NC86

PINLY

Phospholipase A2 inhibitor

Increase

Increase

1

 

P01019

ANGT

Angiotensinogen

Increase

Increase

1

–f

P01834

IGKC

Immunoglobulin kappa constant

Increase

Increase

2

 

P01857

IGHG1

Immunoglobulin heavy constant gamma 1

Increase

Increase

1

 

P02768

ALBU

Serum albumin

Increase

Increase

2

 

P26038

MOES

Moesin

Increase

Increase

1

 

Q92743

HTRA1

Serine protease HTRA1

Increase

Increase

1

 

Proteins increased the day after exposure

P01023

A2MG

Alpha-2-macroglobulin

 

Increase

2

 

P01859

IGHG2

Immunoglobulin heavy constant gamma 2

 

Increase

1

f

P02760

AMBP

Protein AMBP

 

Increase

1

f

P02787

TRFE

Serotransferrin

 

Increase

1

f

P02790

HEMO

Hemopexin

 

Increase

1

 

P04040

CATA

Catalase

 

Increase

1

 

P04792

HSPB1

Heat shock protein beta-1

 

Increase

1

 

P07737

PROF1

Profilin-1

 

Increase

1

f

P08123

CO1A2

Collagen alpha-2

 

Increase

1

f

P08603

CFAH

Complement factor H

 

Increase

1

 

P15924

DESP

Desmoplakin

 

Increase

1

f

P18510

IL1RA

Interleukin-1 receptor antagonist protein

 

Increase

1

f

P27169

PON1

Serum paraoxonase

 

Increase

1

 

P31947

1433S

14-3-3 protein sigma

 

Increase

1

 

P36955

PEDF

Pigment epithelium-derived factor

 

Increase

1

 

Q04118

PRB3

Basic salivary proline-rich protein 3

 

Increase

1

f

Q07157

ZO1

Tight junction protein ZO-1

 

Increase

1

f

P11684

UTER

Uteroglobin

 

Increase

2

 

Proteins decreased immediately after exposure but increased the day after exposure

P07602

SAP

Prosaposin

Decrease

Increase

2

f,f

Q96IY4

CBPB2

Carboxypeptidase B2

Decrease

Increase

2

 

P01877

IGHA2

Immunoglobulin heavy constant alpha 2

Decrease

Increase

1

f

Proteins decreased immediately after exposure 

O75556

SG2A1

Mammaglobin-B

Decrease

 

1

 

P01008

ANT3

Antithrombin-III

Decrease

 

1

 

P01034

CYTC

Cystatin-C

Decrease

 

2

 

P01036

CYTS

Cystatin-S

Decrease

 

1

 

P04406

G3P

Glyceraldehyde-3-phosphate dehydrogenase

Decrease

 

1

 

P0DMV9

HS71B

Heat shock 70 kDa protein 1B

Decrease

 

1

 

P13796

PLSL

Plastin-2

Decrease

 

1

 

P14780

MMP9

Matrix metalloproteinase-9

Decrease

 

1

 

P15311

EZRI

Ezrin

Decrease

 

1

 

P15502

ELN

Elastin

Decrease

 

1

 

P37802

TAGL2

Transgelin-2

Decrease

 

1

 

Q13822

ENPP2

Ectonucleotide pyrophosphatase

Decrease

 

1

 

Q14019

COTL1

Coactosin-like protein

Decrease

 

1

 

Q14508

WFDC2

WAP four-disulfide core domain protein 2

Decrease

 

1

 

Q9GZZ8

LACRT

Extracellular glycoprotein lacritin

Decrease

 

1

 

Q9Y624

JAM1

Junctional adhesion molecule A

Decrease

 

1

 

P04080

CYTB

Cystatin-B

Decrease

 

1

 

Q6UXB2

VCC1

C-X-C motif chemokine 17

Decrease

 

1

 

Q6P5S2

CF058

Protein LEG1 homolog

Decrease

 

2

f,f

Proteins decreased immediately after exposure and the day after exposure

O75594

PGRP1

Peptidoglycan recognition protein 1

Decrease

Decrease

2

 

P05164

PERM

Myeloperoxidase

Decrease

Decrease

2

 

P06727

APOA4

Apolipoprotein A-IV

Decrease

Decrease

1

f

Proteins decreased the day after exposure

P02749

APOH

Beta-2-glycoprotein 1

 

Decrease

1

f

P02751

FINC

Fibronectin

 

Decrease

1

 

P08311

CATG

Cathepsin G

 

Decrease

1

 

P20061

TCO1

Transcobalamin-1

 

Decrease

1

f

Q6UWW0

LCN15

Lipocalin-15

 

Decrease

1

f

Relative quantification of individual samples generated this list of the significantly changed proteins after welding fume exposure compared to filtered air, and these were statistically evaluated with the Wilcoxon signed rank test both immediately after exposure and the day after exposure, p < 0.05

f Peptides not significantly changed after multiple test correction (FDR), p > 0.036

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

Protein

Mean slope in NL

Mean slope in MilliQ water

LOD (fmol/µL)

Exposure

Mean; median; (max–min) concentration (fmol/µL)

Statistical evaluation

Before exposure

After exposure

Day after exposure

Beta-value

95% CI

p value

MMP9

0.291

0.275

0.043

A

0.26; 0.14; (0.02–0.61)

0.14; 0.13; (0.02–0.37)

0.19; 0.13; (0.02–0.36)

0.77

0.63–0.93

0.015

B

0.16; 0.15; (0.02–0.42)

0.16; 0.11 (0.03–0.52)

0.18; 0.16; (0.08–0.4)

   

A1AT

0.153

0.158

0.205

A

0.51; 0.37; (0.23–1.28)

0.50; 0.37; (0.23–1.6)

0.44; 0.52; (0.14–0.94)

1.20

0.77–1.87

0.369

B

0.67; 0.51; (0.3–2.12)

0.52; 0.41; (0.21–1.84)

0.56; 0.45; (0.2–1.48)

   

A2 MG

0.184

0.153

0.075

A

0.52; 0.43; (0.09–0.13)

0.42; 0.20; (0.08–1.32)

0.35; 0.29; (0.05–0.95)

0.95

0.73–1.23

0.669

    

B

0.49; 0.32; (0.16–1.13)

0.47; 0.24; (0.09–1.55)

0.48; 0.51; (0.11–1.44)

   

APOB

0.164

0.177

0.356

A/B

NA

NA

NA

Elastin

0.432*

0.339

0.339

A/B

NA

NA

NA

Mean slope of the calibration curve for each quantitative protein and the calculated LOD (determined from blank samples)

NL nasal lavage, A welding fume particle exposure, B filtered air exposure, NA no quantitative data could be obtained

* Big variation between the slopes

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
Fig. 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)

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

Symbol

Protein

AC

Expr. fold change

Expected

LXR/RXR activation (positive z-score 1.26)

AMBP

Alpha-1-microglobulin/bikunin precursor

P02760

1.23

Up

APOA4

Apolipoprotein A4

P06727

− 1.41

Up

APOH

Apolipoprotein H

P02749

− 1.54

Up

CLU

Clusterin

P10909

1.19

Up

HPX

Hemopexin

P02790

1.26

Up

IL1RN

Interleukin 1 receptor antagonist

P18510

1.30

Down

MMP9

Matrix metallopeptidase 9

P14780

− 1.29

Down

PON1

Paraoxonase 1

P27169

1.20

Up

SAA1

Serum amyloid A1

P0DJI8

1.28

Up

TTR

Transthyretin

P02766

1.38

Up

Acute phase response signaling (positive z-score 0.378)

A2M

Alpha-2-macroglobulin

P01023

1.27

Up

AMBP

Alpha-1-microglobulin/bikunin precursor

P02760

1.23

Down

APOH

Apolipoprotein H

P02749

− 1.54

Down

FN1

Fibronectin 1

P02751

− 1.25

Up

HPX

Hemopexin

P02790

1.26

Up

IL1RN

Interleukin 1 receptor antagonist

P18510

1.30

Up

SAA1

Serum amyloid A1

P0DJI8

1.28

 

SERPING1

Serpin family G member 1

P05155

1.24

 

TTR

Transthyretin

P02766

1.38

Down

The significantly changed proteins found with the LMM statistical method showed that these canonical pathways were activated due to welding fume exposure. AC: accession number. Expected: The changed direction of the protein for this canonical pathway to be activated are expected to go in e certain direction when its activated

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
Fig. 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

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 [3739], 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 [4042]. 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

Declarations

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.

Open AccessThis 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.

Authors’ Affiliations

(1)
Division of Occupational and Environmental Medicine, Lund University, Lund, Sweden
(2)
Occupational and Environmental Medicine Center, Department of Clinical and Experimental Medicine, Linköping University, Linköping, Sweden
(3)
Department of Design Sciences, Ergonomics and Aerosol Technology, Lund University, Lund, Sweden

References

  1. WHO. Hazard prevention and control in the work environment: airborne dust. 1999.Google Scholar
  2. OSHA. Occupational safety and health guidance manual for hazardous waste site activities. 1986.Google Scholar
  3. NIOSH. NIOSH pocket guide to chemical hazards. 2007.Google Scholar
  4. WHO. WHO guidelines for indoor air quality: selected pollutants. Copenhagen: WHO Regional office for Europe. 2010.Google Scholar
  5. Antonini JM. Health effects of welding. Crit Rev Toxicol. 2003;33:61–103.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.Google Scholar
  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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedGoogle Scholar
  20. Schmidt A, Forne I, Imhof A. Bioinformatic analysis of proteomics data. BMC Syst Biol. 2014;8:S3.View ArticlePubMedPubMed CentralGoogle 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.View ArticleGoogle 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.View ArticlePubMedGoogle 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.Google Scholar
  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.View ArticleGoogle 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.View ArticleGoogle 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.View ArticleGoogle Scholar
  28. de Una-Alvarez J. The beta-binomial SGoF method for multiple dependent tests. Stat Appl Genet Mol Biol. 2012; 11:Article 14.Google Scholar
  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.View ArticlePubMedGoogle 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.View ArticleGoogle Scholar
  31. Gooley JJ, Chua EC-P. Diurnal regulation of lipid metabolism and applications of circadian lipidomics. J Genet Genomics. 2014;41:231–50.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedPubMed CentralGoogle Scholar
  36. Meng X-M, Nikolic-Paterson DJ, Lan HY. TGF-[beta]: the master regulator of fibrosis. Nat Rev Nephrol. 2016;12:325–38.View ArticlePubMedGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.View ArticleGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedGoogle Scholar
  42. Giannandrea M, Parks WC. Diverse functions of matrix metalloproteinases during fibrosis. Dis Models Mech. 2014;7:193–203.View ArticleGoogle 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.View ArticleGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.View ArticleGoogle Scholar
  47. Atkinson JJ, Senior RM. Matrix metalloproteinase-9 in lung remodeling. Am J Respir Cell Mol Biol. 2003;28:12–24.View ArticlePubMedGoogle 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.View ArticlePubMedGoogle 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.View ArticleGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.PubMedPubMed CentralGoogle Scholar
  54. Fiancette R, Vincent-Fabert C, Guerin E, Trimoreau F, Denizot Y. Lipid mediators and human leukemic blasts. J Oncol. 2011;2011:7.View ArticleGoogle 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.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedGoogle 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.View ArticlePubMedPubMed CentralGoogle Scholar
  58. Gabay C, Kushner I. Acute-phase proteins and other systemic responses to inflammation. N Engl J Med. 1999;340:448–54.View ArticlePubMedGoogle Scholar
  59. Kholodenko BN. Cell signalling dynamics in time and space. Nat Rev Mol Cell Biol. 2006;2006(7):165–76.View ArticleGoogle Scholar

Copyright

© The Author(s) 2018

Advertisement