- Open Access
Smoking induces sex-specific changes in the small airway proteome
Respiratory Research volume 22, Article number: 234 (2021)
Cigarette smoke triggers many cellular and signaling responses in the lung and the resulting inflammation plays a central role in smoke-related lung diseases, such as COPD. We explored the effects of smoking on the small airway proteome in samples obtained by collection of exhaled particles with the aim to identify specific proteins dysregulated by smoking.
Exhaled particles were obtained from 38 current smokers, 47 former smokers and 22 healthy controls with the PExA method. 120 ng of sample was collected from individual subjects and analyzed with the SOMAscan proteomics platform. General linear model-based statistics were performed.
Two hundred and three proteins were detected in at least half of 107 total samples. Active smoking exerted a significant impact on the protein composition of respiratory tract lining fluid (RTLF), with 81 proteins altered in current smokers compared to never smokers (p < 0.05, q < 0.124). Among the proteins most clearly discriminating between current and never smokers were sRAGE, FSTL3, SPOCK2 and protein S, all of them being less abundant in current smokers. Analysis stratified for sex unveiled sex differences with more pronounced proteomic alterations due to active smoking in females than males. Proteins whose abundance was altered by active smoking in women were to a larger extent related to the complement system. The small airway protein profile of former smokers appeared to be more similar to that observed in never smokers.
The study shows that smoking has a strong impact on protein expression in the small airways, and that smoking affects men and women differently, suggesting PExA sampling combined with high sensitivity protein analysis offers a promising platform for early detection of COPD and identification of novel COPD drug targets.
More than 8 million deaths a year around the world are attributable to tobacco use . Smoking is the most important risk factor for chronic obstructive pulmonary disease (COPD), with 15–36% of smokers developing COPD [2, 3]. The gaseous and particulate matter (tar) phase of cigarette smoke contain approximately 4700 different toxic substances, many of which are highly reactive oxygen/nitrogen species (RONS) . Furthermore, a single puff of cigarette smoke is estimated to contain more than 1014 oxidizing molecules  and smoking one cigarette exposes the human respiratory tract to 10,000–40,000 μg of particulate matter (PM). These particles have a mean diameter of < 1 μm which allows a high rate of deposition in the distal parts of the human lung , especially in the transition zone of small airways. Free radicals, aldehydes, ketones and other species found in cigarette smoke may induce oxidative damage in many types of biological macromolecules, compromising their structure and/or function . The inhalation of cigarette smoke triggers many cellular and signalling responses in the lung and the resulting inflammation plays a central role in smoke-related lung diseases . The mechanisms by which smoking contributes to inflammation and tissue damage are, however, not completely understood . Several proteomic studies applied to sputum, lung tissue and bronchoalveolar lavage (BAL) fluid from never, former and current smokers [7, 9] have identified smoke-induced alterations in proteins involved in the response to oxidative stress and inflammation as well as in extracellular matrix organization and wound healing [9, 10]. The complement system has been proposed an important driver of this inflammation  and its role in the pathogenesis of COPD is emerging [12, 13]. Several complement factors are synthesized locally in the lungs by the alveolar type II cells , providing a basis for local complement system activation and initiation of inflammation. The activation of the complement system leads to opsonisation, phagocytosis and recruitment of inflammatory cells to the infected or injured area . Further, exposure to cigarette smoke has been reported to activate the alternative complement pathway , whereas inhibition of the alternative pathway is protective in preclinical models of lung injury .
In addition, analysis of lung epithelium gene expression has shown that whereas most smoking-related changes reverse following smoking cessation, some appear to be permanent, even in long-term former smokers [17,18,19].
It has also been reported that smoke-related pathogenic mechanisms in COPD may differ between women and men [20,21,22,23], with women appearing more susceptible to the adverse effects of smoking as compared to men . Studies have shown that female smokers have a faster annual lung function decline  and a greater loss of FEV1 per pack-year of smoking , report a higher level of dyspnoea and are more likely to have bronchial hyperresponsiveness than men . There is also a higher prevalence of females among patients with severe, early-onset COPD . The increased susceptibility to cigarette smoke in women may have several aetiologies. First, there may be a genetic predisposition for smoking-related lung damage that is sex-specific . Second possible explanation may be a dose-dependent effect. The airways of women are smaller and thus each cigarette represents a proportionately greater exposure . Third, there may be hormonally mediated differences in tobacco-smoke metabolism, resulting in increased oxidative stress in female lungs .
The small airways may contribute substantially to the pathogenesis of COPD as they represent a “silent zone” where disease can progress over many years without causing symptoms and being detected . The respiratory tract lining fluid (RTLF) is a protective interface between the respiratory epithelium and the external environment, and in small airways it consists mainly of lung surfactant, but also contains other biomolecules released from small airway epithelial cells [32, 33]. Changes in the small airway RTLF may reflect biological processes in this important lung compartment, which can be studied by collecting particles in exhaled air, a novel sampling method allowing non-invasive retrieval of biological material from the small airways [33,34,35,36]. The molecular composition of exhaled particles (PEx) has been explored previously and 124 different proteins were identified by LC/MS in pooled samples. A comparison of the identified PEx proteins with published bronchoalveolar lavage (BAL) proteomic data showed a high degree of overlap, with 103 (83%) of the PEx proteins having previously been detected in BAL .
In the present study we explored the effects of smoking on the small airway proteome in samples obtained by the collection of particles in exhaled air with the aim to identify dysregulation of specific proteins associated with smoking and to determine whether such changes are reversible following smoking cessation.
Study design and participants
One hundred and seven subjects with a median age of 61 years, (51% were women), were examined in 2016 and 2017. All subjects were recruited from our previous studies or by an advertisement in a daily paper. The cohort consisted of 38 current smokers (CS), 47 former smokers (FS) and 22 healthy never smokers (NS). The inclusion criterion for never smoking controls was post-bronchodilation FEV1/FVC > 0.70.
Subjects were defined as current smokers if they had smoked on a regular daily basis for at least one year at the time of examination. Former smokers were defined as those that had not smoked in the last 12 months but had smoked on a regular daily basis before that. Those who had never smoked on a regular basis were classified as never smokers.
Participants provided written informed consent prior to the measurements and the Regional Ethics Committee at the University of Gothenburg approved the study (442-17 and 390-06).
All subjects filled out a questionnaire on medical history, symptoms and use of medication and a detailed smoking history was obtained, as well as information about second-hand exposure to cigarette smoke. Pack-years were calculated by multiplying the average number of packs smoked per day with the duration of smoking in years. Blood samples were obtained and analysed for hsCRP and white blood cell differential.
All subjects were instructed to withdraw from short-acting bronchodilators and long-acting bronchodilators at least 6 h and 24 h prior to the examination, respectively.
Spirometry was performed using a Spirare spirometer (Spirare, Stockholm, Sweden) before and after bronchodilation with 400 µg of salbutamol in accordance with the ATS/ERS criteria . Forced vital capacity (FVC), forced expired volume in one second (FEV1) an FEV1/FVC ratio were expressed as a percentage of the reference value (% pred) according to Brisman et al. (note the corrigendum) .
Exhaled particles (PEx) were collected using the PExA instrument (PExA AB, Gothenburg, Sweden), as previously described [32, 39]. Study subjects breathed via a mouthpiece and a two-way, non-re-breathing valve into the instrument which consists of a thermostated box (36 °C) containing an optical particle counter (Grimm Aerosol Technik GmbH & Co, Ainring, Germany) and an impactor (Dekati Ltd, Tampere, Finland). The measured particle sizes cover diameters between 0.41 and 2.98 µm. Subjects inhaled HEPA-filtered air for a minimum of three breaths before the sampling in order to remove particles from ambient air. All subjects wore a nose clip throughout the procedure. A standardized breathing maneuver was used [35, 40], starting with an exhalation at normal flow rate to residual volume, breath holding for 5 s, followed by a maximal inhalation to total lung capacity, immediately followed by a normal exhalation to functional residual capacity. Exhalation flow was measured by an ultrasonic flow meter (OEM flow sensor; Spiroson-AS, Medical Technologies, Zürich, Switzerland), enabling visualization of the expiratory flow and volume . Between breathing maneuvers, the subject breathed particle-free air tidally for 30 to 60 s. Each sampling session continued until 120 ng of exhaled particles were collected. After collection the sample holder was transferred to a clean air room and the substrate was cut out with a scalpel from the sample holder area and placed in Millipore Ultrafree-MC LH Centrifugal Filter insert (FC30LH25) and stored at − 80 °C for subsequent extraction and SOMAscan analysis. True blank samples were generated by applying the same sample handling procedure as for real samples but without collecting PExA sample from the study subjects.
SOMAscan analysis and processing of data
PEx sample preparation is described in detail in Additional file 1. Prior to SOMAscan analysis, the volume of sample buffer was adjusted to reach the same concentration of PEx in all samples in order to normalize the samples for the differences in the collected amount of PEx.
SOMAscan (SomaLogic Inc, Boulder, USA) is an aptamer-based proteomics platform that uses slow off-rate modified DNA aptamers (SOMAmers) as high affinity protein capture reagents to simultaneously quantify more than 1300 human proteins in all types of protein extracts. The SOMAscan assay quantitatively transforms the proteins present in a biological sample into a specific SOMAmer-based DNA signal, expressed as relative fluorescent units (RFUs), which is directly proportional to the amount of target protein in the initial sample .
Intra-plate and inter-plate normalization were performed by SomaLogic according to their SOMAscan assay good laboratory practice (GLP) data quality-control procedures. Sample data was first normalized to remove hybridization variation within a run (Hybridization Control Normalization), followed by Calibration Normalization to remove assay differences between runs. SomaLogic’s data normalization procedures have been described previously in more details .
Limit of detection (LOD) was calculated as the mean plus 3 standard deviations based on two blank samples. Proteins with RFU values > LOD in more than 50% of the samples were considered for further analyses.
Statistical analyses of the protein data were performed using general linear model-based statistics (Qlucore Omics Explorer 3.6 software, Qlucore AB, Lund, Sweden). SOMAscan data was log2 transformed before the analysis to achieve normal distribution. General linear model (GLM), with each variable normalized to mean 0 and variance of 1, was used to determine differences in protein abundance between NS, FS and CS and all the analyses were adjusted for the investigator performing the PExA measurements, as well as the age of the subjects. Statistical analyses are described in more details in Additional file 1. Group comparisons of SOMAscan data were considered hypothesis free and proteins with p-value < 0.05 were considered to be of interest in this explorative study. Linear regression was used to assess correlations between pack-years, time since smoking cessation and protein levels. Statistical analysis of clinical and demographic data was performed using IBM SPSS Statistics for Windows, version 26 (IBM Corp., Armonk, N.Y., USA) with the significance level set to p < 0.05.
Demographic and clinical characteristics
Female former smokers had significantly lower post-FEV1/FVC and post-MMEF/FVC ratio compared to never smokers, while no significant difference was observed between female current and never smokers. The results also indicated that post-FEV1/FVC and post-MMEF/FVC ratio was lower in male current and former smokers compared to never smokers. The number of pack-years in former smokers was significantly lower compared to current smokers both in males and females. A higher number of blood neutrophils was observed in CS compared to NS and FS in males. Similar pattern was observed in females, although the difference did not reach statistical significance (Table 1).
Effects of smoking on small airway proteome
All in all, 203 proteins were detected in at least 50% of the 107 samples (Table S1, Additional file 2). A comparison of proteins detected in our samples with previously published SOMAscan BAL proteomic data showed a high degree of overlap, with 180 (89%) proteins having previously been detected in BAL with the SOMAscan analysis . Using GLM-based statistics and adjusting for age and variance introduced by the investigator revealed statistically significant differences in 120 proteins between at least two of the groups (q < 0.05). As demonstrated by clustering analysis based on smoking status, these proteins mainly distinguished CS from FS and NS (Fig. 1).
Linear regression (adjusted for age) in FS and CS revealed a significant correlation between the relative abundance of 70 proteins and number of pack-years (p < 0.05, ∣R∣ ≥ 0.22).
Current smokers versus never smokers and sex differences
Active smoking exerted a significant impact on the protein composition of respiratory tract lining fluid (RTLF), with relative abundance of 81 proteins altered in CS compared to NS (p < 0.05, q < 0.124) (Table S2, Additional file 2). The majority of the differentially abundant proteins (62 proteins) were less abundant in CS compared to NS as exemplified by soluble receptor for advanced glycation end products (sRAGE), follistatin-related protein 3 (FSTL3), testican-2 (SPOCK2) and protein S (PROS1), which were among the proteins most clearly distinguishing CS from NS (Fig. 2). The proteins are described in more details in Table 2. This observation was true for both females and males (Fig. 3 and Table 3). We also observed significant differences in the effect of smoking on the RTLF protein profiles between females and males. When stratifying for sex, the relative abundance of 58 proteins was significantly altered in female CS as compared to NS (p < 0.05, q ≤ 0.167), while 27 proteins were significantly altered in male CS (p < 0.05, q ≤ 0.370; 8 with q < 0.2).
Former smokers versus never smokers
When comparing FS to NS, the levels of 31 proteins were altered, with all proteins being less abundant in FS (p < 0.05, q ≤ 0.31; 14 with q < 0.2). Stratifying for sex, 7 proteins were altered in female FS (p < 0.05, q ≤ 0.96) and 9 proteins were altered in male FS (p < 0.05, q ≤ 0.77) as compared to NS. Although p-values were significant in the analysis stratified by sex, none of the proteins in females and only one protein in males were found to be significantly altered after correction for multiple testing (q-value < 0.2).
Effects of smoking cessation on the protein profile of small airways
To determine whether the changes in protein profile caused by active smoking persist after smoking cessation, proteins altered in CS were compared to those remaining altered in FS compared to NS (Table 4). Among the 62 proteins less abundant in CS as compared to NS, four (MRC1, CD55, ST2, AK1) remained significantly less abundant in FS, while the levels of the rest of these proteins were more similar to the levels observed in NS. Of the 19 proteins more abundant in CS as compared to NS, none remained more abundant in FS as compared to NS, twelve of these proteins were found to be significantly less abundant in FS as compared to NS (Table 4), and seven proteins were present at the same levels in FS as in NS.
To further study the effect of smoking cessation on the small airway protein profile, the relationship between relative protein abundance and time since smoking cessation was assessed using linear regression adjusted for age. The relative abundance of 20 proteins in males and 12 proteins in females correlated with time since smoking cessation (p < 0.05). When only looking at proteins that were decreased in CS compared to NS, the abundance of 9 proteins increased with time after smoking cessation in male FS as compared to 6 proteins in female FS (Table 5).
Several complement factors were less abundant in female CS compared to NS (not observed in male). Of these, C3 and C5 were significantly less abundant in female CS (Fig. 4). No significant differences were observed in anaphylatoxin C5a between CS and NS in females (p = 0.891) or males (p = 0.672), but a higher C5a/C5 ratio was observed in female CS (log2(fold change) = 0.63, p = 0.006). Complement factors C2, C6, C7 and C9 were less abundant in female CS compared to NS, as well as factor D which was significantly less abundant in both female and male CS (Table 3). Factor H (Fig. 4) and factor I (log2(fold change) = − 0.43, p = 0.035) were also less abundant in female CS, and complement decay-accelerating factor (CD55) was less abundant in female and male CS compared to NS (Fig. 4 and Table 3).
In the present study we used a novel, non-invasive method to obtain biological samples from the small airways enabling assessment of the effects of smoking on the small airway proteome. Clear differences in the small airways protein profiles between never, former and current smokers were observed. These alterations were mainly detected in current smokers, while the protein profile of former smokers appeared to have returned to that observed in never smokers, with the exception of four proteins (MRC1, CD55, ST2, AK1) that remained downregulated. Further stratification revealed sex-associated differences with more pronounced alterations in female current smokers than in males. Based on observed changes in protein levels, the complement system was identified as being a biological pathway strongly affected by current smoking, particularly in female smokers, and a model of the small airway inflammation in smokers was extrapolated (Fig. 5).
The levels of 81 proteins were significantly altered in current smokers as compared to never smokers, and the proteins most clearly distinguishing current smokers from never smokers were soluble receptor for advanced glycation end products (sRAGE), testican-2 (SPOCK2), follistatin-related protein 3 (FSTL3) and protein S (PROS1). sRAGE has anti-inflammatory properties and deficiency of circulating sRAGE is associated with various human diseases. Decreased levels of plasma sRAGE have been previously observed in smokers with or without COPD, as compared to never smokers, and the levels were significantly correlated to lung function [44, 45]. In line with these observations, we observed significantly decreased levels of sRAGE in the lining fluid from small airways in current smokers, whereas, the levels of sRAGE in former smokers appeared to be similar to those in never smokers. These findings are consistent with decreased sRAGE levels being a marker of deficient inflammatory control, as previously suggested .
Levels of FSTL3, a protein structurally and functionally related to follistatin, were decreased in current smokers. Follistatin has anti-inflammatory properties and has been seen to be decreased in cigarette smoke-exposed human bronchial epithelial cells and administration of follistatin was found to attenuate cigarette smoke-induced airway inflammation in mice .
To our knowledge, differential abundance of testican-2 (SPOCK2) has not been previously reported to be associated with smoking or COPD. SPOCK2 is known to modulate matrix metalloproteinases expression and activation  with potential implications for small airway remodelling in context of COPD . We did observe a slightly lower abundance of MMP-2 (also known as gelatinase A) in current smokers as compared to never smokers and a significant positive correlation between the levels of SPOCK2 and MMP-2 (not shown). Previous studies of the association between MMP-2 and smoking or COPD show conflicting results, some reporting increased MMP-2 expression in COPD patients [50,51,52], whereas others have observed decreased small airway levels of MMP-2 in COPD patients . However, other metalloproteinases more often associated with COPD, such as MMP-9 and MMP-12, were not detected in this study.
Protein S is involved in the inhibition of coagulation, clearance of apoptotic cells and inhibition of inflammation. It binds and localizes C4BP, an important regulator of complement activation, to apoptotic cells where C4BP can down-regulate complement activation and therefore inhibit inflammation at the surface of apoptotic cells. We observed that protein S was significantly less abundant both in female and male current smokers as compared to never smokers. These findings suggest that the removal of apoptotic cells is compromised in the small airways of active smokers and are consistent with previous studies reporting the role of decreased clearance of apoptotic cells in the pathogenesis of COPD . Furthermore, due to a decreased protein S in active smokers, the localization of C4BP to apoptotic cells may be diminished, which may result in excessive complement activation and inflammation in the small airways.
Our results show that active smoking in women causes the most pronounced alterations in the proteins related to the complement system. Female current smokers had decreased levels of the majority of the measured complement factors as well as decreased levels of complement factors H and I, key regulators of the alternative pathway (AP) preventing its excessive activation. Factor H also binds to apoptotic cells to limit the inflammatory potential of complement . Furthermore, decreased levels of another important inhibitor of complement activation, CD55, were seen in the small airways of both female and male current smokers. Unlike the classical and the lectin pathways that are generally activated by recognition of exogenous material, the AP is constitutively active at a low level under normal conditions , and full activation is triggered by the presence of factor C3b and factor B which form a complex that eventually leads to the formation of AP C3 convertase that cleaves C3 into C3a and C3b, starting a positive feedback loop needed to commence AP activation. During inflammation, when cells are damaged, C3b is formed by the classical/lectin pathway, making the AP a powerful amplification mechanism of complement activation. Therefore, AP requires continuously active control mechanisms to maintain homeostasis . Dysfunction of these regulatory proteins has been identified as a cause of several diseases , and decreased abundance of factors H and I in current smokers could potentially lead to chronic over-activation of the AP, which could be an important initiating step in the pathogenesis of small airway inflammation. As observed in current smokers, depletion of complement factors C3 and C5 as well as the concomitant increase in the pro-inflammatory mediator C5a (as C5a/C5 ratio), supports this hypothesis. The observed changes in the C5a/C5 ratio, but not in C5a levels, could be due to the nonspecificity of the SOMAscan aptamer for C5a, as it also binds to C5 albeit with approximately a three-fold lower affinity. A similar pattern was seen for C3a (as C3a/C3 ratio), however, SOMAscan aptamers for C3a are known to also bind to C3 with similar affinity and therefore the results could be influenced by the changes in C3. It has also been shown that exposure of C3 to cigarette smoke extract produces a functionally modified form of the molecule capable of activating the AP .
Heat shock protein 90 (HSP90) and translationally-controlled tumor protein (also known as fortilin or TCTP) were among the few proteins to be increased in current smokers as compared to never smokers. HSP90 plays an important role in the “unfolded protein response” (UPR) , a compensatory response elicited by impaired protein folding due to oxidative stress. The primary function of the UPR is to reduce the accumulation of aberrantly folded proteins and promote cell survival by increasing the expression of genes involved in protein chaperoning and folding, translation and degradation [59, 60]. However, when cells are damaged beyond repair, the UPR promotes elimination of afflicted cells through apoptosis . Previous studies have also reported an increase in the expression of heat shock proteins in the alveolar epithelial cells exposed to cigarette smoke extract . TCTP, a pro-survival molecule, has also been shown to play an important role in the UPR where it protects cells from the apoptotic cell death . Our finding of increased HSP90 and TCTP levels in current smokers as compared to never smokers are consistent with previous studies suggesting cigarette smoke causes protein structural changes and impairs their folding which in turn results in the increased expression of the proteins involved in the UPR in order to promote cell survival.
The majority of alterations due to active smoking appeared to be reversible after smoking cessation. The former smokers had quit smoking at least 12 months prior to the start of the study, and time from smoking cessation ranged from one to 45 years, with a median of 11 years in females and 19 years in males. Only four proteins (MRC1, CD55, ST2 and AK1) remained less abundant in former smokers as compared to never smokers, indicative of permanent dysregulation. Reduced expression of MRC1 has been previously associated with reduced removal of apoptotic cells in patients with COPD, which may lead to lung tissue damage and COPD progression due to release of cytotoxic products from apoptotic cells . ST2 is a receptor for IL-33. The membrane-bound ST2 binds IL-33, inducing pro-inflammatory immune responses, and the soluble form functions as a decoy receptor neutralizing IL-33 activity and is therefore considered to have an anti-inflammatory function [64, 65]. Persistent alterations of these proteins in smokers after smoking cessation could help explain why former smokers remain at an elevated risk for COPD despite quitting smoking. The proteins induced in current smokers were all less abundant in former smokers, some to levels even lower than observed in never smokers. These findings are consistent with an earlier observation that most cigarette smoke-related changes in the lung proteome are reversible .
Finally, we observed a surprisingly small overlap of proteins altered due to smoking in females and males, suggesting smoking affects women and men differently and by distinct mechanisms. It has previously been proposed that smoking may have more deleterious effects on lung health in women, and that female smokers may be more susceptible to developing COPD [24, 27]. Sex-specific differences in the BAL cell proteome between healthy smokers and smokers with COPD with more pronounced alterations in females have been reported [20, 23]. When comparing never smokers and smokers with normal lung function, one study reported a significant impact of smoking on the BAL cell proteome with alterations of more than 500 proteins (representing 15 molecular pathways) due to smoking. However, the majority of these alterations were sex-independent .
This study highlights alterations in the composition of the respiratory tract lining fluid caused by smoking, however, there are several limitations that need to be considered when interpreting the results. First of all, the sample size is small and an independent validation cohort is needed to confirm our findings. We also observed variability depending on the investigator performing the PExA sampling, and as a result the statistical analysis had to be adjusted accordingly. Never, former and current smokers were, however, examined in a random order and a systematic error therefore seems unlikely. Another potential cofounder was that never smokers were on average younger than the former smokers and current smokers, and therefore we adjusted the analysis for age. Choosing to include proteins in the analysis with RFU values > LOD in more than 50% of samples could be considered another possible limitation. However, due to the exploratory nature of the study, we chose to select 50% as the limit instead of the more commonly used 70–80%, to open the analytical window and explore proteins that could potentially be missed otherwise. For that same reason and also due to the small sample size, a more inclusive approach was taken and the changes were identified as significant at p-value < 0.05.
The observation that the majority of proteins with altered levels were less abundant in current smokers may be due to decreased protein expression in response to smoking. Another possibility that cannot be excluded is that smoke-induced changes in protein structure may compromise the affinity for the SOMAscan aptamers, leading to weaker signals. However, when comparing proteins detected in current smokers and never smokers, the overlap was significant, indicating that smoking status does not influence the SOMAscan analysis.
Another possible limitation is that the SOMAscan data is normalized based on the total amount of PEx and reported in relative values, which means that reported fold change values may not reflect the true magnitude of differences in the concentration of the differentially abundant proteins. Others ways to normalize the data were considered, such as normalization based on the total protein mass or to a housekeeping protein, but as PEx samples contain a small amount of biological material, it is difficult to measure the total protein concentration in order to use that for normalization. Additionally, the exhaled particle is a novel matrix and knowledge about potential housekeeping proteins that could be used for normalization is limited. Therefore, the data was instead normalized by diluting each sample prior to the SOMAscan analysis to equal concentrations of biological material as reflected by the total mass of the exhaled particles in each individual sample.
It will also be interesting to replicate and extend in a more severe COPD cohort to address whether the small airway proteome changes we observed in current smokers are present also in COPD patients and correlate with disease severity.
In conclusion, the study shows that smoking has a strong impact on protein expression in the small airways, and that smoking affects men and women differently. The observed protein alterations consistent with complement pathway activation in female smokers indicate that PExA proteomics can identify novel COPD disease processes as well as novel biomarkers. Thus, PExA sampling combined with high sensitivity protein analysis offers a promising platform for early detection of COPD and identification of novel COPD drug targets.
Availability of data and materials
The dataset used and analysed during the current study are available from the corresponding author on reasonable request.
Chronic obstructive pulmonary disease
Reactive oxygen/nitrogen species
Respiratory tract lining fluid
Forced expiratory volume in 1 s
Forced vital capacity
Maximal mid-expiratory flow
High-sensitivity C-reactive protein
Relative fluorescent units
Limit of detection
Unfolded protein response
World Health Organization: WHO report on the global tobacco epidemic 2019: offer help to quit tobacco use. https://www.who.int/teams/health-promotion/tobacco-control/who-report-on-the-global-tobacco-epidemic-2019&publication=9789241516204; 2019.
Rennard SI, Vestbo J. COPD: the dangerous underestimate of 15%. Lancet. 2006;367:1216–9.
Løkke A, Lange P, Scharling H, Fabricius P, Vestbo J. Developing COPD: a 25 year follow up study of the general population. Thorax. 2006;61:935.
Borgerding M, Klus H. Analysis of complex mixtures–cigarette smoke. Exp Toxicol Pathol. 2005;57(Suppl 1):43–73.
Pryor WA, Stone K. Oxidants in cigarette smoke. Radicals, hydrogen peroxide, peroxynitrate, and peroxynitrite. Ann N Y Acad Sci. 1993;686:12–27.
Ghio AJ, Hilborn ED, Stonehuerner JG, Dailey LA, Carter JD, Richards JH, Crissman KM, Foronjy RF, Uyeminami DL, Pinkerton KE. Particulate matter in cigarette smoke alters iron homeostasis to produce a biological effect. Am J Respir Crit Care Med. 2008;178:1130–8.
Colombo G, Clerici M, Giustarini D, Portinaro NM, Aldini G, Rossi R, Milzani A, Dalle-Donne I. Pathophysiology of tobacco smoke exposure: recent insights from comparative and redox proteomics. Mass Spectrom Rev. 2014;33:183–218.
Steiling K, Kadar AY, Bergerat A, Flanigon J, Sridhar S, Shah V, Ahmad QR, Brody JS, Lenburg ME, Steffen M, Spira A. Comparison of proteomic and transcriptomic profiles in the bronchial airway epithelium of current and never smokers. PLoS ONE. 2009;4:12.
Titz B, Sewer A, Schneider T, Elamin A, Martin F, Dijon S, Luettich K, Guedj E, Vuillaume G, Ivanov NV, et al. Alterations in the sputum proteome and transcriptome in smokers and early-stage COPD subjects. J Proteomics. 2015;128:306–20.
Mossina A, Lukas C, Merl-Pham J, Uhl FE, Mutze K, Schamberger A, Staab-Weijnitz C, Jia J, Yildirim AÖ, Königshoff M, et al. Cigarette smoke alters the secretome of lung epithelial cells. Proteomics. 2017;17:1600243.
Sun S, Wang H, Zhao G, An Y, Guo Y, Du L, Song H, Qiao F, Yu H, Wu X, et al. Complement inhibition alleviates paraquat-induced acute lung injury. Am J Respir Cell Mol Biol. 2011;45:834–42.
Pandya PH, Wilkes DS. Complement system in lung disease. Am J Respir Cell Mol Biol. 2014;51:467–73.
Marc MM, Korosec P, Kosnik M, Kern I, Flezar M, Suskovic S, Sorli J. Complement factors C3a, C4a, and C5a in chronic obstructive pulmonary disease and asthma. Am J Respir Cell Mol Biol. 2004;31:216–9.
Strunk RC, Eidlen DM, Mason RJ. Pulmonary alveolar type II epithelial cells synthesize and secrete proteins of the classical and alternative complement pathways. J Clin Invest. 1988;81:1419–26.
Ekdahl KN, Mohlin C, Adler A, Aman A, Manivel VA, Sandholm K, Huber-Lang M, Fromell K, Nilsson B. Is generation of C3(H2O) necessary for activation of the alternative pathway in real life? Mol Immunol. 2019;114:353–61.
Kew RR, Ghebrehiwet B, Janoff A. Cigarette smoke can activate the alternative pathway of complement in vitro by modifying the third component of complement. J Clin Investig. 1985;75:1000–7.
Zhang L, Lee JJ, Tang H, Fan Y-H, Xiao L, Ren H, Kurie J, Morice RC, Hong WK, Mao L. Impact of smoking cessation on global gene expression in the bronchial epithelium of chronic smokers. Cancer Prev Res (Philadelphia, PA). 2008;1:112–8.
Beane J, Sebastiani P, Liu G, Brody JS, Lenburg ME, Spira A. Reversible and permanent effects of tobacco smoke exposure on airway epithelial gene expression. Genome Biol. 2007;8:R201.
Hijazi K, Malyszko B, Steiling K, Xiao X, Liu G, Alekseyev YO, Dumas Y-M, Hertsgaard L, Jensen J, Hatsukami D, et al. Tobacco-related alterations in airway gene expression are rapidly reversed within weeks following smoking-cessation. Sci Rep. 2019;9:6978.
Yang M, Kohler M, Heyder T, Forsslund H, Garberg HK, Karimi R, Grunewald J, Berven FS, Nyrén S, Magnus Sköld C, Wheelock ÅM. Proteomic profiling of lung immune cells reveals dysregulation of phagocytotic pathways in female-dominated molecular COPD phenotype. Respir Res. 2018;19:39.
Naz S, Bhat M, Ståhl S, Forsslund H, Sköld CM, Wheelock ÅM, Wheelock CE. Dysregulation of the tryptophan pathway evidences gender differences in COPD. Metabolites. 2019;9:212.
Naz S, Kolmert J, Yang M, Reinke SN, Kamleh MA, Snowden S, Heyder T, Levänen B, Erle DJ, Sköld CM, et al. Metabolomics analysis identifies sex-associated metabotypes of oxidative stress and the autotaxin-lysoPA axis in COPD. Eur Respir J. 2017;49:6.
Kohler M, Sandberg A, Kjellqvist S, Thomas A, Karimi R, Nyrén S, Eklund A, Thevis M, Sköld CM, Wheelock ÅM. Gender differences in the bronchoalveolar lavage cell proteome of patients with chronic obstructive pulmonary disease. J Allergy Clin Immunol. 2013;131:743–51.
Sansores RH, Ramírez-Venegas A. COPD in women: susceptibility or vulnerability? Eur Respir J. 2016;47:19–22.
Gan WQ, Man SFP, Postma DS, Camp P, Sin DD. Female smokers beyond the perimenopausal period are at increased risk of chronic obstructive pulmonary disease: a systematic review and meta-analysis. Respir Res. 2006;7:52.
Prescott E, Bjerg AM, Andersen PK, Lange P, Vestbo J. Gender difference in smoking effects on lung function and risk of hospitalization for COPD: results from a Danish longitudinal population study. Eur Respir J. 1997;10:822–7.
Ohar J, Fromer L, Donohue JF. Reconsidering sex-based stereotypes of COPD. Prim Care Respir J. 2011;20:370–8.
Silverman EK, Weiss ST, Drazen JM, Chapman HA, Carey V, Campbell EJ, Denish P, Silverman RA, Celedon JC, Reilly JJ, et al. Gender-related differences in severe, early-onset chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2000;162:2152–8.
Han MK, Postma D, Mannino DM, Giardino ND, Buist S, Curtis JL, Martinez FJ. Gender and chronic obstructive pulmonary disease: why it matters. Am J Respir Crit Care Med. 2007;176:1179–84.
Sin DD, Cohen SB, Day A, Coxson H, Paré PD. Understanding the biological differences in susceptibility to chronic obstructive pulmonary disease between men and women. Proc Am Thorac Soc. 2007;4:671–4.
Hogg JC, Pare PD, Hackett TL. The contribution of small airway obstruction to the pathogenesis of chronic obstructive pulmonary disease. Physiol Rev. 2017;97:529–52.
Larstad M, Almstrand AC, Larsson P, Bake B, Larsson S, Ljungstrom E, Mirgorodskaya E, Olin AC. Surfactant protein A in exhaled endogenous particles is decreased in chronic obstructive pulmonary disease (COPD) patients: a pilot study. PLoS ONE. 2015;10:e0144463.
Bredberg A, Gobom J, Almstrand AC, Larsson P, Blennow K, Olin AC, Mirgorodskaya E. Exhaled endogenous particles contain lung proteins. Clin Chem. 2012;58:431–40.
Behndig AF, Mirgorodskaya E, Blomberg A, Olin A-C. Surfactant Protein A in particles in exhaled air (PExA), bronchial lavage and bronchial wash - a methodological comparison. Respir Res. 2019;20:214.
Almstrand AC, Bake B, Ljungstrom E, Larsson P, Bredberg A, Mirgorodskaya E, Olin AC. Effect of airway opening on production of exhaled particles. J Appl Physiol. 1985;2010(108):584–8.
Larsson P, Bake B, Wallin A, Hammar O, Almstrand A-C, Lärstad M, Ljungström E, Mirgorodskaya E, Olin A-C. The effect of exhalation flow on endogenous particle emission and phospholipid composition. Respir Physiol Neurobiol. 2017;243:39–46.
Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, Crapo R, Enright P, van der Grinten CP, Gustafsson P, et al. Standardisation of spirometry. Eur Respir J. 2005;26:319–38.
Brisman J, Kim JL, Olin AC, Toren K, Bake B. Spirometric reference equations for Swedish adults. Clin Physiol Funct Imaging. 2017;37:640–5.
Almstrand AC, Ljungstrom E, Lausmaa J, Bake B, Sjovall P, Olin AC. Airway monitoring by collection and mass spectrometric analysis of exhaled particles. Anal Chem. 2009;81:662–8.
Holmgren H, Gerth E, Ljungstrom E, Larsson P, Almstrand AC, Bake B, Olin AC. Effects of breath holding at low and high lung volumes on amount of exhaled particles. Respir Physiol Neurobiol. 2013;185:228–34.
Gold L, Ayers D, Bertino J, Bock C, Bock A, Brody EN, Carter J, Dalby AB, Eaton BE, Fitzwater T, et al. Aptamer-based multiplexed proteomic technology for biomarker discovery. PLoS ONE. 2010;5:e15004.
Candia J, Cheung F, Kotliarov Y, Fantoni G, Sellers B, Griesman T, Huang J, Stuccio S, Zingone A, Ryan BM, et al. Assessment of variability in the SOMAscan assay. Sci Rep. 2017;7:14248.
DeBoer EM, Wagner BD, Popler J, Harris JK, Zemanick ET, Accurso FJ, Sagel SD, Deterding RR. Novel application of aptamer proteomic analysis in cystic fibrosis bronchoalveolar lavage fluid. Proteomics. 2019;13:1800085.
Iwamoto H, Gao J, Pulkkinen V, Toljamo T, Nieminen P, Mazur W. Soluble receptor for advanced glycation end-products and progression of airway disease. BMC Pulm Med. 2014;14:68.
Pouwels SD, Klont F, Kwiatkowski M, Wiersma VR, Faiz A, van den Berge M, Horvatovich P, Bischoff R, ten Hacken NHT. Cigarette smoking acutely decreases serum levels of the chronic obstructive pulmonary disease biomarker sRAGE. Am J Respir Crit Care Med. 2018;198:1456–8.
Smith DJ, Yerkovich ST, Towers MA, Carroll ML, Thomas R, Upham JW. Reduced soluble receptor for advanced glycation end-products in COPD. Eur Respir J. 2011;37:516.
Verhamme FM, Bracke KR, Amatngalim GD, Verleden GM, Van Pottelberge GR, Hiemstra PS, Joos GF, Brusselle GG. Role of activin-A in cigarette smoke-induced inflammation and COPD. Eur Respir J. 2014;43:1028.
Liu G, Ren F, Song Y. Upregulation of SPOCK2 inhibits the invasion and migration of prostate cancer cells by regulating the MT1-MMP/MMP2 pathway. PeerJ. 2019;7:e7163.
Churg A, Zhou S, Wright JL. Matrix metalloproteinases in COPD. Eur Respir J. 2012;39:197.
Mahor D, Kumari V, Vashisht K, Galgalekar R, Samarth RM, Mishra PK, Banerjee N, Dixit R, Saluja R, De S, Pandey KC. Elevated serum matrix metalloprotease (MMP-2) as a candidate biomarker for stable COPD. BMC Pulm Med. 2020;20:302.
Baraldo S, Bazzan E, Zanin ME, Turato G, Garbisa S, Maestrelli P, Papi A, Miniati M, Fabbri LM, Zuin R, Saetta M. Matrix metalloproteinase-2 protein in lung periphery is related to COPD progression. Chest. 2007;132:1733–40.
Chen Y, Chen P, Hanaoka M, Droma Y, Kubo K. Enhanced levels of prostaglandin E2 and matrix metalloproteinase-2 correlate with the severity of airflow limitation in stable COPD. Respirology. 2008;13:1014–21.
Gosselink JV, Hayashi S, Elliott WM, Xing L, Chan B, Yang L, Wright C, Sin D, Paré PD, Pierce JA, et al. Differential expression of tissue repair genes in the pathogenesis of chronic obstructive pulmonary disease. Am J Respir Crit Care Med. 2010;181:1329–35.
Yun JH, Henson PM, Tuder RM. Phagocytic clearance of apoptotic cells: role in lung disease. Expert Rev Respir Med. 2008;2:753–65.
Martin M, Leffler J, Smolag KI, Mytych J, Bjork A, Chaves LD, Alexander JJ, Quigg RJ, Blom AM. Factor H uptake regulates intracellular C3 activation during apoptosis and decreases the inflammatory potential of nucleosomes. Cell Death Differ. 2016;23:903–11.
Noris M, Remuzzi G. Overview of complement activation and regulation. Semin Nephrol. 2013;33:479–92.
Thurman JM, Holers VM. The central role of the alternative complement pathway in human disease. J Immunol. 2006;176:1305.
Marcu MG, Doyle M, Bertolotti A, Ron D, Hendershot L, Neckers L. Heat shock protein 90 modulates the unfolded protein response by stabilizing IRE1alpha. Mol Cell Biol. 2002;22:8506–13.
Kelsen SG, Duan X, Ji R, Perez O, Liu C, Merali S. Cigarette smoke induces an unfolded protein response in the human lung: a proteomic approach. Am J Respir Cell Mol Biol. 2008;38:541–50.
Jorgensen E, Stinson A, Shan L, Yang J, Gietl D, Albino AP. Cigarette smoke induces endoplasmic reticulum stress and the unfolded protein response in normal and malignant human lung cells. BMC Cancer. 2008;8:229.
Pinkaew D, Chattopadhyay A, King MD, Chunhacha P, Liu Z, Stevenson HL, Chen Y, Sinthujaroen P, McDougal OM, Fujise K. Fortilin binds IRE1α and prevents ER stress from signaling apoptotic cell death. Nat Commun. 2017;8:18.
Somborac-Bačura A, Rumora L, Novak R, Rašić D, Dumić J, Čepelak I, Žanić-Grubišić T. Differential expression of heat shock proteins and activation of mitogen-activated protein kinases in A549 alveolar epithelial cells exposed to cigarette smoke extract. Exp Physiol. 2018;103:1666–78.
Hodge S, Hodge G, Ahern J, Jersmann H, Holmes M, Reynolds PN. Smoking alters alveolar macrophage recognition and phagocytic ability: implications in chronic obstructive pulmonary disease. Am J Respir Cell Mol Biol. 2007;37:748–55.
Liew FY, Pitman NI, McInnes IB. Disease-associated functions of IL-33: the new kid in the IL-1 family. Nat Rev Immunol. 2010;10:103–10.
Liew FY. Cigarette smoke resets the Alarmin IL-33 in COPD. Immunity. 2015;42:401–3.
Yang M, Kohler M, Heyder T, Forsslund H, Garberg HK, Karimi R, Grunewald J, Berven FS, Magnus Sköld C, Wheelock ÅM. Long-term smoking alters abundance of over half of the proteome in bronchoalveolar lavage cell in smokers with normal spirometry, with effects on molecular pathways associated with COPD. Respir Res. 2018;19:40.
Robinson AB, Stogsdill JA, Lewis JB, Wood TT, Reynolds PR. RAGE and tobacco smoke: insights into modeling chronic obstructive pulmonary disease. Front Physiol. 2012;3:89.
Ahn N, Kim WJ, Kim N, Park HW, Lee SW, Yoo JY. The Interferon-inducible proteoglycan testican-2/SPOCK2 functions as a protective barrier against virus infection of lung epithelial cells. J Virol. 2019;93:1.
Bartholin L, Guindon S, Martel S, Corbo L, Rimokh R. Identification of NF-kappaB responsive elements in follistatin related gene (FLRG) promoter. Gene. 2007;393:153–62.
Webb JH, Blom AM, Dahlbäck B. Vitamin K-dependent protein S localizing complement regulator C4b-binding protein to the surface of apoptotic cells. J Immunol. 2002;169:2580–6.
Kask L, Trouw LA, Dahlbäck B, Blom AM. The C4b-binding protein-protein S complex inhibits the phagocytosis of apoptotic cells. J Biol Chem. 2004;279:23869–73.
Ali Assad N, Sood A. Leptin, adiponectin and pulmonary diseases. Biochimie. 2012;94:2180–9.
Gao N, Wang Y, Zheng C-M, Gao Y-L, Li H, Li Y, Fu T-T, Xu L-L, Wang W, Ying S, Huang K. β(2)-Microglobulin participates in development of lung emphysema by inducing lung epithelial cell senescence. Am J Physiol Lung Cell Mol Physiol. 2017;312:L669–77.
Picard D. Heat-shock protein 90, a chaperone for folding and regulation. Cell Mol Life Sci. 2002;59:1640–8.
Lee CH, Goag EK, Lee SH, Chung KS, Jung JY, Park MS, Kim YS, Kim SK, Chang J, Song JH. Association of serum ferritin levels with smoking and lung function in the Korean adult population: analysis of the fourth and fifth Korean National Health and Nutrition Examination Survey. Int J Chron Obstruct Pulmon Dis. 2016;11:3001–6.
Yamasaki K. Eeden SFv: lung macrophage phenotypes and functional responses: role in the pathogenesis of COPD. Int J Mol Sci. 2018;19:582.
Kambara K, Ohashi W, Tomita K, Takashina M, Fujisaka S, Hayashi R, Mori H, Tobe K, Hattori Y. In vivo depletion of CD206+ M2 macrophages exaggerates lung injury in endotoxemic mice. Am J Pathol. 2015;185:162–71.
We would like to thank the research nurses and the biomedical technicians at Occupational and environmental medicine unit at University of Gothenburg for clinical work and technical support. We would also like to extend our appreciation to all subjects for their participation in the present study.
Open access funding provided by University of Gothenburg. This study was funded by Forte (Grant No. 2016-00639), Swedish Heart and Lung Foundation (Grant No. 21080209) and supported by the Centre for Allergy Research Highlights Asthma Markers of Phenotype (ChAMP) consortium, funded by the Swedish Foundation for Strategic Research, the Karolinska Institute, AstraZeneca & Science for Life Laboratory Joint Research Collaboration, and the Vårdal Foundation as well as the Swedish Research Council and the Linnaeus university (Grant Nos. 2016-2075-5.1 and 2016-04519).
Ethics approval and consent to participate
The Regional Ethics Committee at the University of Gothenburg approved the study (442-17 and 390-06) and participants provided written informed consent prior to the measurements.
Consent for publication
ACO is a chair-holder and a board member of PExA AB (www.PEXA.se). JÖ reports personal fees from PExA AB during the conduct of the study and was employed by PExA AB while writing the manuscript, but not during the planning and completion of the study. HKO and BG are full-time employees of AstraZeneca.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
A detailed description of PEx sample preparation and analysis.
List of 203 proteins with RFU values > LOD in more than 50% of the samples that was used for further analysis. Table S2. Proteins significantly altered in current smokers as compared to never smokers in the joint model.
About this article
Cite this article
Kokelj, S., Östling, J., Georgi, B. et al. Smoking induces sex-specific changes in the small airway proteome. Respir Res 22, 234 (2021). https://doi.org/10.1186/s12931-021-01825-6
- Small airways
- Respiratory tract lining fluid
- Exhaled particles