- Open Access
Identification of differentially expressed proteins in the injured lung from zinc chloride smoke inhalation based on proteomics analysis
Respiratory Research volume 20, Article number: 36 (2019)
Lung injury due to zinc chloride smoke inhalation is very common in military personnel and leads to a high incidence of pulmonary complications and mortality. The aim of this study was to uncover the underlying mechanisms of lung injury due to zinc chloride smoke inhalation using a rat model.
Methods: Histopathology analysis of rat lungs after zinc chloride smoke inhalation was performed by using haematoxylin and eosin (H&E) and Mallory staining. A lung injury rat model of zinc chloride smoke inhalation (smoke inhalation for 1, 2, 7 and 14 days) was developed. First, isobaric tags for relative and absolute quantization (iTRAQ) and weighted gene co-expression network analysis (WGCNA) were used to identify important differentially expressed proteins. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were used to study the biological functions of differentially expressed proteins. Then, analysis of lung injury repair-related differentially expressed proteins in the early (day 1 and day 2) and middle-late stages (day 7 and day 14) of lung injury after smoke inhalation was performed, followed by the protein-protein interaction (PPI) analysis of these differentially expressed proteins. Finally, the injury repair-related proteins PARK7 and FABP5 were validated by immunohistochemistry and western blot analysis.
Morphological changes were observed in the lung tissues after zinc chloride smoke inhalation. A total of 27 common differentially expressed proteins were obtained on days 1, 2, 7 and 14 after smoke inhalation. WGCNA showed that the turquoise module (which involved 909 proteins) was most associated with smoke inhalation time. Myl3, Ckm, Adrm1 and Igfbp7 were identified in the early stages of lung injury repair. Gapdh, Acly, Tnni2, Acta1, Actn3, Pygm, Eno3 and Tpi1 (hub proteins in the PPI network) were identified in the middle-late stages of lung injury repair. Eno3 and Tpi1 were both involved in the glycolysis/gluconeogenesis signalling pathway. The expression of PARK7 and FABP5 was validated and was consistent with the proteomics analysis.
The identified hub proteins and their related signalling pathways may play crucial roles in lung injury repair due to zinc chloride smoke inhalation.
Military operations are complex and are always very dangerous because of toxic industrial chemicals and materials such as smoke bombs. Smoke bombs generate large amounts of particles and toxic gases during combustion and explosion. Therefore, lung smoke inhalation injuries are common in military personnel. It has been noted that smoke inhalation is the leading cause of lung injury . In November 2014, 15 soldiers sustained lung injuries during smoke bomb inhalation [2, 3]. In addition, 3 soldiers sustained lung injuries during smoke bomb inhalation, and one of the soldiers died from the severe lung injury in January 2017 .
Generally, lung injury from smoke inhalation is defined as the inhalation of thermal or chemical irritants (such as zinc chloride), with a very high incidence of pulmonary complications and mortality [5,6,7]. A cough, difficulty breathing, nausea, vomiting and abnormal liver function are the main symptoms of lung injuries due to smoke inhalation. Despite many years of intense research, the molecular mechanisms involved in the pathology of lung injury from smoke inhalation are poorly illuminated [1, 2, 5]. Rats are model animals suitable for studying the effects of inhaled smoke on lung injury . To study the potential mechanisms of smoke inhalation on lung injury, we developed a rat model with high doses of zinc chloride smoke inhalation to cause lung injury.
It is noted that isobaric tags for relative and absolute quantization (iTRAQ), mass spectrometry-based quantitative approaches and weighted gene co-expression network analysis (WGCNA) have largely contributed to proteomics investigations . Applying the above techniques, we first analysed the differentially expressed proteins from days 1, 2, 7 and 14 after zinc chloride smoke inhalation in our rat model. Then, we further identified lung injury repair-related, differentially expressed proteins in the early and middle-late stages of lung injury. Our study may be helpful in understanding the potential molecular mechanisms in lung injuries due to smoke inhalation.
Materials and methods
The lung injury rat model due to smoke inhalation
Male Wistar rats with a weight of 200 g were supplied by the Academy of Military Medical Sciences of the PLA and were raised carefully in accordance with the ethical guidelines of the National Institutes of Health on animal care. All experimental procedures were approved by the Academy of Military Medical Sciences of the PLA. The homemade smoke chamber consisted of 5 g of smoke bomb materials (aluminium powder, zinc oxide, hexachloroethane, silicon powder and black powder). The smoke bomb was placed into a special rectangular iron pot with a diameter of 2 cm, a magnesium strip was placed in the pot to ignite, and the smoke bomb was burned for approximately 45 s. The smoke chamber is 1 m long, 0.8 m wide, 0.6 m high and 0.48 m3 in volume. The smoke box allows samples to be taken quickly, and these samples are measured by a gas sampler to determine the concentration of the individual smoke components in the gas. The concentration of the 5-g smoke bombs was 10.625 g.m-3. The gas sampling method was used to determine the concentration of the zinc ions in the smoke generated by the smoke bomb, and a concentration of 4.996 mg/L was determined. The male Wistar rats were randomly divided into the control group (n = 8) or the exposed group (n = 164), and the male Wistar rats in the exposed group were only exposed once on the first day for 5 min, and then the smoke was cleared for 10 min. In this study, the animals were only exposed once to a high dose because all human exposures are high-dose accidents and are not repeated occupational exposures. Lung tissues were collected 1, 2, 7, 14, 21 and 28 days after the initial exposure.
The right, upper lobe of the lung was fixed in 10% formalin, dehydrated, embedded in paraffin and cut into 6 μm sections. After deparaffinization, the tissues were haematoxylin and eosin (H&E) stained or Mallory stained, and the morphological lesions and changes in lung tissues were observed under a light microscope.
Lung tissues were harvested and lysed in RIPA lysis buffer with protease inhibitors for 30 min at 4 °C. Lysates were then centrifuged for 10 min at 12000 rpm at 4 °C, and the insoluble debris was discarded. The protein in the lung tissues was extracted, followed by alkylation treatment. The concentration test of the extracted protein was then performed. Sodium dodecyl sulfate was used for protein electrophoresis [10, 11]. Protein digestion was conducted using the filter-aided sample preparation (FASP) method described in Wisniewsk et al. . After protein enzymolysis, the iTRAQ reagent was applied to mark peptides . iTRAQ labelling was performed by using the iTRAQ reagent-8plex Multiplex Kit (AB Sciex California, USA) according to the manufacturer’s instructions. Strong cation exchange chromatography was used to predissociate a balanced mixture of marked peptides . The pooled samples were analysed by 2D LC-MS/MS. LC–electrospray ionization (ESI)–MS/MS analysis was conducted on a liquid chromatography machine coupled with tandem mass spectrometry (nLC, Thermo Fisher Scientific). A 10-μL aliquot of each fraction was injected for 2D LC-MS/MS analysis. A mass spectrometer (MS), (Orbitrap Fusion, Thermo Scientific) was used to analyse the eluted peptides. The MS data were acquired with the Orbitrap detector by using the fastest mode. For database searching, the differentially expressed proteins were identified with an abundance ratio ≥ 1.5 or an abundance ratio < 0.667 and p-values≤0.05. A heat-map analysis of the differentially expressed protein expression was performed by using pheatmap in R.
Construction of the WGCNA co-expression network
To find significant modules and proteins in the proteome data, the WGCNA package in the R software package (http://www.r-project.org/) was used to construct the co-expression network of the proteins. The rationale behind the WGCNA approach is that highly correlated genes, in terms of their expression values, are very likely to work in the same biological processes and/or pathways so that they can be grouped into a common network module . The type of network is signed, which allows positively or negatively regulated genes to be grouped into different network modules. First, Pearson’s correlation matrices were calculated between pairwise genes . Then, the similarity matrix was transformed into an adjacency matrix (AM) using a power β = 10 based on the scale-free topology criterion described in the WGCNA package documents . Afterwards, the adjacency matrix was transformed into a topological overlap matrix and was, in turn, converted into a dissimilarity topological overlap matrix, from which a dendrogram was mapped via hierarchical clustering. The clusters were obtained from the dendrogram by using the dynamic tree cutting technique . The next step was to further identify the most important modules for analysis. We first calculated the eigengene of each module, which is defined as the first principal component of its representation matrix. Then, we used these eigengenes to compute module-smoke inhalation time associations. The module-smoke inhalation time association for a given module and smoke inhalation time is the correlation between the eigengene of the module and the smoke inhalation time. Lastly, we chose six modules in association with the smoke inhalation times. In addition, to explore the association between the modules and the smoke inhalation times, a correlation analysis was performed. The modules with absolute values of correlation coefficients > 0.8 and p < 0.001 were identified to be smoke inhalation time-associated modules. The Cytoscape software (http://cytoscape.org/) was used for the visualization of the most significant time-associated module with a threshold > 0.5. Cytoscape cores are available for layout, scripting, file formatting and linking the network to databases . Nodes represent proteins, and edges represent interactions between two proteins . The pivotal nodes in the network were identified based on their degrees of connectivity.
Functional analysis of proteins in the turquoise module
To study the biological function of proteins in the turquoise module, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed by using the online software GeneCodis3 (http://genecodis.cnb.csic.es/analysis). The false discovery rate (FDR) value of each enriched pathway was assigned with -log10 transformation. The threshold of the false discovery rate (FDR) < 0.05 was set as the criteria for statistical significance.
Identification of hub proteins in the turquoise module
To further obtain key proteins, we screened hub differentially expressed proteins in the turquoise module. The proteins with Coefficient.Time + Coefficient.turquoise> 0.9, p.Time + p.turquoise< 0.01 were identified to be hub proteins. A heat-map analysis of the expression of hub proteins was performed by using pheatmap in R language.
Lung injury repair-related protein analysis in the early and middle-late stages of zinc chloride smoke inhalation
To identify potential injury repair-related proteins and related signalling pathways in the early (day 1 and day 2) and middle-late (day 7 and day 14) stages of lung injury after smoke inhalation, we analysed the differentially expressed proteins on day 1 and day 2 compared with the control and on day 7 and day 14 compared with the control, respectively. All these key differentially expressed proteins intersected with those in the most significant smoke time-associated module.
Protein-protein interaction analysis
To understand the protein interaction between key differentially expressed proteins in the early (4 differentially expressed proteins) and middle-late (96 differentially expressed proteins) stages of smoke inhalation, the STRING database (http://www.string-db.org/) was utilized to select interacting protein pairs. Through the PPI relationship obtained in the previous step, the protein-protein network (PPI) was visualized by Cytoscape software (http://cytoscape.org/). The degree was analysed by Cytoscape software. In the network, the nodes represent proteins, and the edges represent the interactions between the two proteins.
Immunohistochemistry validation of differentially expressed injury repair-related proteins
In this study, PARK7, a differentially expressed protein on day 15 after smoke inhalation, was selected for immunohistochemistry (IHC) analysis on days 1, 2, 7, 14, 21 and 28. Then, 5 μm thick continuous sections were incubated with an anti-PARK7 antibody produced in rabbit (1:200 dilution; Abcam) followed by goat anti-rabbit immunoglobulin antibody conjugated by horseradish peroxidase (1:200 dilution; Vector). Then, the slides were visualized using diaminobenzidine (DAB) substrate (Vector). The number of positive cells was calculated under the microscope.
Western blot analysis of differentially expressed injury repair-related proteins
In this study, the differentially expression proteins PARK7 and FABP5 after day 14 of smoke inhalation were selected for western blot (WB) analysis on days 1, 2, 7 and 14. Protein extracts were prepared by homogenization with protein lysates. Protein concentrations were quantified with an Enhanced BCA Protein Assay Kit (Thermo Fisher Scientific, Rochester, N.Y., USA). Protein samples were run on a 12% sodium dodecyl sulfate polyacrylamide gel for electrophoresis and were transferred to a polyvinylidene fluoride membrane. The membrane was blocked with non-fat dried milk for 1 h and was incubated with the following primary antibodies: rabbit anti-human FABP5 (1:1000 dilution), rabbit anti-human PARK7 (1:500 dilution) and mouse monoclonal IgG1 to GAPDH (1:3000 dilution). Primary antibody binding was detected using horseradish peroxidase-conjugated secondary antibody (1:10000 dilution). Detection was achieved with the ECL chemiluminescence kit (Pierce). Bands from three separate western blots were analysed by Quantiscan software (Biosoft, Great Shelford, Cambridge, UK).
Immunohistochemistry and western blot statistical analyses were performed using GraphPad Prism (GraphPad Software, La Jolla, CA). Student’s t-test was used to assess statistical significance. Statistical significance was ascribed to a p value < 0.05. Graphs are presented as the mean ± SD. Immunohistochemistry and western blot experiments were repeated independently at least three times.
Histopathology analysis of rat lungs after zinc chloride smoke inhalation
The histopathology analysis results of rat lungs after zinc chloride smoke inhalation are shown in Fig. 1. The control group showed a normal lung structure. On day 1 after smoke inhalation, vascular wall thickening, oedema, a large number of black particles around the lumen, and neutrophilic granulocytes increased. On day 2 after smoke inhalation, exudative pneumonia, a large amount of exudation in the alveolar cavity and infiltration of lymphocytes around the bronchioles were observed. On day 8 after smoke inhalation, bleeding, emphysema, obvious cavities around the blood vessels, fibroblast proliferation, and massive lymphocytes in the alveolar cavity were observed. On day 14 after smoke inhalation, the alveolar cavity was obviously dilated, and scattered inflammatory foci were visible. On day 21 after smoke inhalation, the bronchial wall dilated and developed into an interstitial emphysema. The alveolar sac was extremely dilated, the interstitial sac was obviously widened, and there was a large amount of inflammatory cell infiltration. On day 28 after smoke inhalation, fibrous tissue hyperplasia, pulmonary fibrosis, eosinophilic infiltration of alveolar cavities, and alveolar expansion were observed.
Mallory staining analysis of rat lungs after smoke inhalation
The Mallory method was used for tricolour staining of lung tissues. The Mallory staining results of rat lungs after smoke inhalation are shown in Fig. 2. The dark blue substance was significantly increased in lung tissues on days 1 (B), 2 (C), 7 (D), 14 (E), 21 (F) and 28 (G) after smoke inhalation compared with that in the control group (A), indicating that pulmonary fibrosis occurred after smoke inhalation.
Identification of differentially expressed proteins
In this study, a proteomics analysis was used for differentially expressed protein identification in lung tissues on days 1, 2, 7 and 14 after smoke inhalation. A total of 3377 differentially expressed proteins were identified. Among these, 27 common differentially expressed proteins were obtained on days 1, 2, 7 and 14 after smoke inhalation. A Venn diagram of the differentially expressed proteins on day 1 vs day 0, day 2 vs day 0, day 7 vs day 0 and day 14 vs day 0 is shown in Fig. 3. In addition, 27 common differentially expressed proteins are listed in Table 1. Figure 4 shows the heat map of the 27 common differentially expressed proteins.
WGCNA co-expression network analysis
In the WGCNA co-expression network analysis, a total of 6 modules (green, blue, yellow, brown, turquoise and grey) were identified (Fig. 5). To further explore the correlation between the modules and the smoke inhalation time, a correlation analysis was performed, and the modules with a correlation coefficient > 0.8 and p < 0.001 were considered to be related. The correlation analysis showed that the turquoise module (which had 909 proteins) was associated with smoke inhalation time, and the other modules were not related to smoke inhalation time (p > 0.05), as shown in Table 2. In addition, the visualization of the turquoise module, which consists of 76 nodes and 370 edges, is shown in Fig. 6.
Functional annotation of the turquoise module
According to the GO enrichment analysis (Fig. 7), protein folding, metabolic processes and proteolysis were the most significantly enriched biological processes; the cytoplasm, mitochondria, and cytosol were the most significantly enriched cellular components; nucleotide binding, ATP binding and peptidase activity were the most significantly enriched molecular functions. In addition, based on the KEGG enrichment analysis, proteasome, protein processing in endoplasmic reticulum, lysosome and glycolysis/gluconeogenesis were enriched signaling pathways (Table 3).
Hub protein identification in the turquoise module
To further identify important proteins, we screened the hub proteins in the turquoise module. A total of 55 hub proteins were identified (Table 4). Additionally, the heat map of 55 hub proteins is shown in Fig. 8. In Fig. 8, 55 differentially expressed hub proteins were significantly upregulated on day 7 and day 14 (middle-late stage) after smoke inhalation, which suggests that these proteins may be associated with lung injury repair from smoke inhalation.
Injury repair-related differentially expressed proteins analysis in the early and middle-late stages of zinc chloride smoke inhalation
The histopathology results showed obvious lung injury on day 1 and day 2 after zinc chloride smoke inhalation. Therefore, it is necessary to find smoke inhalation lung injury repair-related differentially expressed proteins and related signalling pathways in the early and middle-late stages of smoke inhalation. A total of 4 differentially expressed proteins, including Myl3, Ckm, Adrm1 and Igfbp7, were identified after overlapping the early stage of injury with the turquoise module after smoke inhalation. In addition, a total of 96 proteins were identified after overlapping the middle-late stage of injury from smoke inhalation with the turquoise module. A Venn diagram of related differentially expressed proteins in the early and middle-late stages of smoke inhalation is shown in Fig. 9. In addition, after KEGG enrichment analysis of these 96 differentially expressed proteins, we found that the lysosome, glycolysis/gluconeogenesis pathway (involving Eno3 and Tpi1), pentose phosphate pathway, amino sugar and nucleotide sugar metabolism, dust collection duct acid secretion, fructose and mannose metabolism, starch and sucrose metabolism and citrate cycle (TCA cycle) were the enriched signalling pathways (Table 5).
The PPI network of identified key differentially expressed proteins in the early stage (4 differentially expressed proteins were identified after overlapping the early stage of smoke inhalation with the turquoise module) and middle-late stage (96 differentially expressed proteins were identified after overlapping the middle-late of smoke inhalation with the turquoise module) of smoke inhalation was determined by Cytoscape software. There were 61 nodes and 205 edges in the network (Fig. 10). Interestingly, Eno3 (degree = 20), Gapdh (degree = 19), Acly (degree = 17), Umps (degree = 17), Tpi1 (degree = 15), Pygm (degree = 15), Tnni2 (degree = 15), Actn3 (degree = 15), and Acta1 (degree = 15) were proteins with a high degree of differential expression.
Immunohistochemistry validation of PARK7
In this study, PARK7, a differentially expressed protein after day 14 of smoke inhalation, was selected for IHC analysis on days 1, 2, 8, 14, 21 and 28 (Fig. 11). The brown colour represents the positive staining of PARK7. The expression of PARK7 was mainly located in the cytoplasm and the nucleus. Moreover, the density of PARK7 staining in the nucleus increased over time after smoke inhalation. The strong positive staining signal (especially in the control group) was located in the bronchial epithelial cells and the alveolar epithelial cells. The validation result of PARK7 was consistent with the results of the proteomics analysis.
Western blot analysis of PARK7 and FABP5
In this study, the differentially expressed proteins PARK7 and FABP5 on day 14 after smoke inhalation were selected for WB analysis on days 1, 2, 7 and 14 (Fig. 12). The expression of PARK7 was highest in the control group and decreased with the smoke inhalation time. The expression of FABP5 was highest on day 14. The validation result is consistent with the proteomics analysis.
Since the first case of smoke bomb inhalation injury was reported in 1945, cases of smoke-induced inhalation lung injuries have been continuously reported . Zinc chloride is often used for battle, military exercises, and fire-fighting training . Zinc chloride aerosols are highly hygroscopic in the respiratory tract. Both the inhalation of the smoke and zinc chloride lead to pulmonary oedema, alveolitis in the early phase, late-phase pulmonary fibrosis, and often fatal acute respiratory distress syndrome in confined spaces [20, 22]. Smoke inhalation injury has a complex pathogenesis, and there is no specific treatment. Current animal models for studying inhaled lung injuries mainly include rats, ovines and rabbits [23,24,25]. Rats are easy to breed and are easy to develop into model animals for studying various diseases. Zhu et al. reported that the rat model of inhaled smoke is particularly suitable for studying the effects of long-term inhalation of smoke on lung injury . With high-throughput analysis, the emergence of profiles has become an effective method to identify differentially expressed genes and proteins in a variety of diseases. These profiles help to explore pathogenesis and to develop biomarkers. A high-throughput sequencing analysis has suggested that circRNAs were differentially expressed during acute lung injury induced by smoke inhalation in rat lung tissue . To our knowledge, this is the first time that differential proteins have been identified after an inhaled lung injury in rats using proteomic analysis. Here, we established a model of rat zinc chloride smoke inhalation to define proteins associated with inhaled lung injury. In this study, we explored the histopathology of rat lungs after zinc chloride smoke inhalation. On day 1 of smoke inhalation, lung injury occurred. On day 7 after smoke inhalation, injury repair began. On day 14 after smoke inhalation, repair to the injury was obvious. This suggested that smoke inhalation had an effect on the lungs. In addition, after the WGCNA co-expression network analysis, we finally obtained a total of 55 differentially expressed hub proteins. Interestingly, 55 hub proteins were upregulated on days 7 and 14 after smoke inhalation, which suggests that these proteins may be associated with injury repair. Based on this, we further analysed the differentially expressed proteins related to lung injury repair in the early (day 1 and day 2) and middle-late (day 7 and day 14) stages after smoke inhalation. Finally, a total of 4 proteins, including Ckm, Adrm1, Igfbp7 and Myl3, were identified in the early stages of lung injury after smoke inhalation.
Creatine kinase muscle (Ckm) is an early, sensitive tissue leakage biomarker of skeletal muscle injury . Adhesion regulating molecule 1 (Adrm1) encodes the integral plasma membrane protein involved in proteolysis and deubiquitination [27,28,29,30]. It has been reported that Adrm1 is downregulated in ageing-associated idiopathic pulmonary fibrosis . Furthermore, downregulation of Adrm1 at the mRNA and protein levels leads to slower wound-healing . Insulin-like growth factor binding protein 7 (Igfbp7), a cell cycle arrest biomarker, is released in the earliest stages of acute kidney injury, which has been suggested to be an important predictor of acute kidney injury following cardiac surgery [33,34,35]. We hypothesize that Ckm, Adrm1 and Igfbp7 may also be involved in the early repair of lung injury by zinc chloride smoke inhalation.
Myosin light polypeptide 3 (Myl3) is a biomarker for monitoring muscle injury and recovery . It has been found that Myl3 is significantly downregulated in spinal cord injuries . In addition, it is a causal gene for hypertrophic cardiomyopathy . It has been demonstrated that Myl3 is a plasma marker of cardiac injury . In this study, we found that Myl3 was differentially expressed in the early stages of lung injury after smoke inhalation. It has been found that cardiac dysfunction results from the inhalational of Cl2; this exposure may result in pulmonary hypertension due to severe lung injury . In addition, cardiac output varies in patients with acute lung injury with or without right ventricular dilatation . We speculate that Myl3 may be involved in lung injury repair in the early stage after smoke inhalation.
Based on the KEGG enrichment analysis, the proteasome, protein processing in endoplasmic reticulum, the lysosome and glycolysis/gluconeogenesis were enriched signalling pathways and organelles. Aldehyde dehydrogenase 7A1 (ALDH7A1), Enolase 3 (Eno3) and triosephosphate isomerase 1 (Tpi1) were enriched in the glycolysis/gluconeogenesis signalling pathways. ALDH7A1 is a member of the ALDH superfamily. The primary function of these enzymes is the detoxification of endobiotic and xenobiotic aldehyde compounds into their corresponding weak carboxylic acids. Giacalone et al. reported that ALDH7A1 is the potential target of ALDH enzymes in the treatment of lung cancer . Eno3 is a highly expressed protein in the lungs . In guinea pig lungs, Eno3 is involved in glycolysis/gluconeogenesis . It has been reported that Eno3 is upregulated in traumatic brain injury . The expression of triosephosphate isomerase 1 (Tpi1) is low in smoking-induced emphysema lung models . It is worth mentioning that Tpi1 is involved in the glycolysis or gluconeogenesis pathway in acute lung injury . The lungs depend on circulation for glucose acquisition. Therefore, the glycolysis/gluconeogenesis signalling pathway is important in the physiological activities of the lungs. Our results suggest that ALDH7A1, Eno3 and Tpi1 may be involved in glucose metabolism in the middle-late stage of lung injury repair.
According to the PPI, we found highly expressed proteins, including Gapdh, Acly, Tnni2, Acta1, Actn3, Pygm, Eno3 and Tpi1. It is noted that all these proteins were differentially expressed in the middle-late stage of injury after smoke inhalation. It has been shown that glyceraldehyde-3-phosphate dehydrogenase (Gapdh) has an anti-inflammatory function in the lipopolysaccharide-induced severe acute lung injury mouse model . ATP citrate lyase (Acly), a lipid biosynthesis transcript, plays an important role in chronic obstructive pulmonary disease . It has been pointed out that Acly is downregulated in acrolein-induced acute lung injury . Troponin I, skeletal fast 2 (Tnni2) plays a crucial role in lung veins . In rats, Tnni2 is involved in angiogenesis in the early stage of skin incision wounds . Actin, alpha 1 skeletal muscle (Acta1) is a representative marker of myocardial adaptive response to damage . Significantly, Acta1 is a differentially regulated gene in nickel-induced acute lung injury . Muscle glycogen phosphorylase (Pygm) is involved in glycogen metabolism . The homozygous mutation in the Pygm gene leads to glycogen metabolism disorder in patients with acute kidney injury .
Parkinson’s disease 7 (Park7) is autosomal recessive with an early onset and is a metabolic enzyme. The expression of Park7 is downregulated in lung development . It has been found that Park7 is downregulated in asbestos-exposed epithelial and mesothelial lung cell lines . In this study, we found that Park7 was upregulated on day 14 after zinc chloride smoke inhalation. In addition, the immunohistochemistry results showed that the expression of PARK7 was highest in the control group and decreased with time after zinc chloride smoke inhalation.
In addition to the proteins listed above, we also found some differentially expressed proteins (such as Lta4h and Fabp5) in the middle-late stages of injury after zinc chloride smoke inhalation. It is noted that Lta4h is a hub protein in the turquoise module of WGCNA analysis (Table 4). Leukotriene A4 hydrolase (Lta4h) polymorphisms are associated with baseline lung function . It has been found that Lta4h is related to chronic obstructive pulmonary disease . The epidermal fatty acid binding protein 5 (Fabp5) is expressed in lung tissues, especially in alveolar macrophages . A previous study indicated that Fabp5 was downregulated in the airway epithelial cells of smokers with chronic obstructive pulmonary disease . In addition, a microarray experiment showed the expression of Fabp5 in a mouse lung during nickel-induced lung injury . Our WB results showed that the expression of FABP5 was the highest on day 14, which was consistent with the proteomics analysis. In addition, 10 differentially expressed proteins, including Gapdh, Acly, Tnni2, Acta1, Actn3, Pygm, Eno3, Tpi1, Lta4h, Fabp5 and the glycolysis/gluconeogenesis signalling pathway, may play roles in the middle-late stage of zinc chloride smoke inhalation lung injury repair. However, there are limitations to our study. The biological functions of the identified differentially expressed proteins were not thoroughly investigated. Some in vivo or in vitro experiments, such as cell experiments, are needed.
We first analysed differentially expressed proteins on days 1, 2, 7 and 14 after zinc chloride smoke inhalation in rats. Then, we further identified lung injury repair-related differentially expressed proteins in the early and middle-late stages of lung injury. The identified key proteins and related signalling pathways may play crucial roles in smoke inhalation lung injury repair.
False discovery rate
Ematoxylin and eosin
Isobaric tags for relative and absolute quantization
Kyoto Encyclopedia of Genes and Genomes
Weighted gene co-expression network analysis
Zhang F, Li MY, Lan YT, Wang CB. Imbalance of Th17/Tregs in rats with smoke inhalation-induced acute lung injury. Sci Rep. 2016;6:21348.
Xie F, Zhang X, Xie L. Prognostic value of serum zinc levels in patients with acute HC/zinc chloride smoke inhalation. Medicine (Baltimore). 2017;96:e8156.
Cao L, Zhang XG, Wang JG, Wang HB, Chen YB, Zhao DH, Shi WF, Xie LX. Pulmonary function test findings in patients with acute inhalation injury caused by smoke bombs. J Thorac Dis. 2016;8:3160–7.
El Idrissi A, van Berkel L, Bonekamp NE, Dalemans DJ, van der Heyden MA. The toxicology of zinc chloride smoke producing bombs and screens. Clin Toxicol (Phila). 2017;55:167–74.
Rehberg S, Maybauer MO, Enkhbaatar P, Maybauer DM, Yamamoto Y, Traber DL. Pathophysiology, management and treatment of smoke inhalation injury. Expert Rev Respir Med. 2009;3:283–97.
Shirani KZ, Pruitt BA Jr, Mason AD Jr. The influence of inhalation injury and pneumonia on burn mortality. Ann Surg. 1987;205:82–7.
Suzuki M, Aikawa N, Kobayashi K, Higuchi R. Prognostic implications of inhalation injury in burn patients in Tokyo. Burns. 2005;31:331–6.
Zhu F, Qiu X, Wang J, Jin Y, Sun Y, Lv T, Xia Z. A rat model of smoke inhalation injury. Inhal Toxicol. 2012;24:356–64.
Ong SE, Mann M. Mass spectrometry-based proteomics turns quantitative. Nat Chem Biol. 2005;1:252–62.
Liu X, Wang J, Gao L, Liu H, Liu C. iTRAQ-based proteomic analysis of neonatal kidney from offspring of protein restricted rats reveals abnormalities in Intraflagellar transport proteins. Cell Physiol Biochem. 2017;44:185–99.
Li Y, Wang X, Ao M, Gabrielson E, Askin F, Zhang H, Li QK. Aberrant Mucin5B expression in lung adenocarcinomas detected by iTRAQ labeling quantitative proteomics and immunohistochemistry. Clin Proteomics. 2013;10:15.
Wisniewski JR, Zougman A, Nagaraj N, Mann M. Universal sample preparation method for proteome analysis. Nat Methods. 2009;6:359–62.
Guo Z, Liu X, Li M, Shao C, Tao J, Sun W. Differential urinary glycoproteome analysis of type 2 diabetic nephropathy using 2D-LC-MS/MS and iTRAQ quantification. J Transl Med. 2015;13:371.
Alaei S, Sadeghi B, Najafi A, Masoudi-Nejad A. LncRNA and mRNA integration network reconstruction reveals novel key regulators in esophageal squamous-cell carcinoma. Genomics. 2018;111:76-89.
Chen PF, Wang F, Nie JY, Feng JR, Liu J, Zhou R, Wang HL, Zhao Q. Co-expression network analysis identified CDH11 in association with progression and prognosis in gastric cancer. Onco Targets Ther. 2018;11:6425–36.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Liu S, Xie F, Xiang X, Dong S, Qu K, Lin T. Identification of differentially expressed genes, lncRNAs and miRNAs which are associated with tumor malignant phenotypes in hepatoblastoma patients. Oncotarget. 2017;8:97554–64.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Zhai K, Kong X, Liu B, Lou J. Bioinformatics analysis of gene expression profiling for identification of potential key genes among ischemic stroke. Medicine (Baltimore). 2017;96:e7564.
Hsu HH, Tzao C, Chang WC, Wu CP, Tung HJ, Chen CY, Perng WC. Zinc chloride (smoke bomb) inhalation lung injury: clinical presentations, high-resolution CT findings, and pulmonary function test results. Chest. 2005;127:2064–71.
Greenfield RA, Brown BR, Hutchins JB, Iandolo JJ, Jackson R, Slater LN, Bronze MS. Microbiological, biological, and chemical weapons of warfare and terrorism. Am J Med Sci. 2002;323:326–40.
Homma S, Jones R, Qvist J, Zapol WM, Reid L. Pulmonary vascular lesions in the adult respiratory distress syndrome caused by inhalation of zinc chloride smoke: a morphometric study. Hum Pathol. 1992;23:45–50.
Ye Z, Liu X, Yang Y, Zhang X, Yu T, Li S, Feng Y, Luo G. The differential expression of novel circular RNAs in an acute lung injury rat model caused by smoke inhalation. J Physiol Biochem. 2018;74:25–33.
Lange M, Hamahata A, Enkhbaatar P, Esechie A, Connelly R, Nakano Y, Jonkam C, Cox RA, Traber LD, Herndon DN, Traber DL. Assessment of vascular permeability in an ovine model of acute lung injury and pneumonia-induced Pseudomonas aeruginosa sepsis. Crit Care Med. 2008;36:1284–9.
Fu ZH, Guo GH, Xiong ZF, Liao X, Liu MZ, Luo J. Early anticoagulation therapy for severe burns complicated by inhalation injury in a rabbit model. Mol Med Rep. 2017;16:7375–81.
Vlasakova K, Lane P, Michna L, Muniappa N, Sistare FD, Glaab WE. Response of novel skeletal muscle biomarkers in dogs to drug-induced skeletal muscle injury or sustained endurance exercise. Toxicol Sci. 2017;156:422–7.
Lamerant N, Kieda C. Adhesion properties of adhesion-regulating molecule 1 protein on endothelial cells. FEBS J. 2005;272:1833–44.
Cherix N, Froquet R, Charette SJ, Blanc C, Letourneur F, Cosson P. A Phg2-Adrm1 pathway participates in the nutrient-controlled developmental response in Dictyostelium. Mol Biol Cell. 2006;17:4982–7.
Hamazaki J, Iemura S, Natsume T, Yashiroda H, Tanaka K, Murata S. A novel proteasome interacting protein recruits the deubiquitinating enzyme UCH37 to 26S proteasomes. EMBO J. 2006;25:4524–36.
Qiu XB, Ouyang SY, Li CJ, Miao S, Wang L, Goldberg AL. hRpn13/ADRM1/GP110 is a novel proteasome subunit that binds the deubiquitinating enzyme, UCH37. EMBO J. 2006;25:5742–53.
Hao X, Du Y, Qian L, Li D, Liu X. Upregulation of long noncoding RNA AP003419.16 predicts high risk of agingassociated idiopathic pulmonary fibrosis. Mol Med Rep. 2017;16:8085–91.
Fejzo MS, Anderson L, von Euw EM, Kalous O, Avliyakulov NK, Haykinson MJ, Konecny GE, Finn RS, Slamon DJ. Amplification target ADRM1: role as an oncogene and therapeutic target for ovarian Cancer. Int J Mol Sci. 2013;14:3094–109.
Kashani K, Al-Khafaji A, Ardiles T, Artigas A, Bagshaw SM, Bell M, Bihorac A, Birkhahn R, Cely CM, Chawla LS, et al. Discovery and validation of cell cycle arrest biomarkers in human acute kidney injury. Crit Care. 2013;17:R25.
Kellum JA, Chawla LS. Cell-cycle arrest and acute kidney injury: the light and the dark sides. Nephrol Dial Transplant. 2016;31:16–22.
Wang Y, Zou Z, Jin J, Teng J, Xu J, Shen B, Jiang W, Zhuang Y, Liu L, Luo Z, et al. Urinary TIMP-2 and IGFBP7 for the prediction of acute kidney injury following cardiac surgery. BMC Nephrol. 2017;18:177.
Liu Y, Wang Y, Teng Z, Zhang X, Ding M, Zhang Z, Chen J, Xu Y. DNA microarray analysis in screening features of genes involved in spinal cord injury. Med Sci Monit. 2016;22:1571–81.
Marian AJ, Salek L, Lutucuta S. Molecular genetics and pathogenesis of hypertrophic cardiomyopathy. Minerva Med. 2001;92:435–51.
Thoonen R, Ernande L, Cheng J, Nagasaka Y, Yao V, Miranda-Bezerra A, Chen C, Chao W, Panagia M, Sosnovik DE, et al. Functional brown adipose tissue limits cardiomyocyte injury and adverse remodeling in catecholamine-induced cardiomyopathy. J Mol Cell Cardiol. 2015;84:202–11.
Zaky A, Bradley WE, Lazrak A, Zafar I, Doran S, Ahmad A, White CW, Dell'Italia LJ, Matalon S, Ahmad S. Chlorine inhalation-induced myocardial depression and failure. Physiol Rep. 2015;3:e12439.
Cepkova M, Kapur V, Ren X, Quinn T, Zhuo H, Foster E, Matthay MA, Liu KD. Clinical significance of elevated B-type natriuretic peptide in patients with acute lung injury with or without right ventricular dilatation: an observational cohort study. Ann Intensive Care. 2011;1:18.
Giacalone NJ, Den RB, Eisenberg R, Chen H, Olson SJ, Massion PP, Carbone DP, Lu B. ALDH7A1 expression is associated with recurrence in patients with surgically resected non-small-cell lung carcinoma. Future Oncol. 2013;9:737–45.
Wu J, Zhou D, Deng C, Wu X, Long L, Xiong Y. Characterization of porcine ENO3: genomic and cDNA structure, polymorphism and expression. Genet Sel Evol. 2008;40:563–79.
Jain R, Dey B, Tyagi AK. Development of the first oligonucleotide microarray for global gene expression profiling in Guinea pigs: defining the transcription signature of infectious diseases. BMC Genomics. 2012;13:520.
Wu P, Zhao Y, Haidacher SJ, Wang E, Parsley MO, Gao J, Sadygov RG, Starkey JM, Luxon BA, Spratt H, et al. Detection of structural and metabolic changes in traumatically injured hippocampus by quantitative differential proteomics. J Neurotrauma. 2013;30:775–88.
Uh ST, Koo SM, Jang AS, Park SW, Choi JS, Kim YH, Park CS. Proteomic differences with and without ozone-exposure in a smoking-induced emphysema lung model. Korean J Intern Med. 2015;30:62–72.
Li Y, Jin S, Lei L, Pan Z, Zou X. Deciphering deterioration mechanisms of complex diseases based on the construction of dynamic networks and systems analysis. Sci Rep. 2015;5:9283.
Takaoka Y, Goto S, Nakano T, Tseng HP, Yang SM, Kawamoto S, Ono K, Chen CL. Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) prevents lipopolysaccharide (LPS)-induced, sepsis-related severe acute lung injury in mice. Sci Rep. 2014;4:5204.
Li J, Zhao P, Yang L, Li Y, Tian Y, Li S, Bai Y. Integrating 3-omics data analyze rat lung tissue of COPD states and medical intervention by delineation of molecular and pathway alterations. Biosci Rep. 2017;37:BSR20170042.
Fabisiak JP, Medvedovic M, Alexander DC, McDunn JE, Concel VJ, Bein K, Jang AS, Berndt A, Vuga LJ, Brant KA, et al. Integrative metabolome and transcriptome profiling reveals discordant energetic stress between mouse strains with differential sensitivity to acrolein-induced acute lung injury. Mol Nutr Food Res. 2011;55:1423–34.
Millino C, Sarinella F, Tiveron C, Villa A, Sartore S, Ausoni S. Cardiac and smooth muscle cell contribution to the formation of the murine pulmonary veins. Dev Dyn. 2000;218:414–25.
Kameyama H, Udagawa O, Hoshi T, Toukairin Y, Arai T, Nogami M. The mRNA expressions and immunohistochemistry of factors involved in angiogenesis and lymphangiogenesis in the early stage of rat skin incision wounds. Leg Med (Tokyo). 2015;17:255–60.
Aztatzi-Aguilar OG, Uribe-Ramirez M, Arias-Montano JA, Barbier O, De Vizcaya-Ruiz A. Acute and subchronic exposure to air particulate matter induces expression of angiotensin and bradykinin-related genes in the lungs and heart: angiotensin-II type-I receptor as a molecular target of particulate matter exposure. Part Fibre Toxicol. 2015;12:17.
Mallakin A, Kutcher LW, McDowell SA, Kong S, Schuster R, Lentsch AB, Aronow BJ, Leikauf GD, Waltz SE. Gene expression profiles of Mst1r-deficient mice during nickel-induced acute lung injury. Am J Respir Cell Mol Biol. 2006;34:15–27.
Yavari A, Bellahcene M, Bucchi A, Sirenko S, Pinter K, Herring N, Jung JJ, Tarasov KV, Sharpe EJ, Wolfien M, et al. Mammalian gamma2 AMPK regulates intrinsic heart rate. Nat Commun. 2017;8:1258.
Zhao X, Li A, Soni M, Muriello MJ, Jones CH, Whittier WL. McArdle disease: a "pediatric" disorder presenting in an adult with acute kidney injury. CEN Case Rep. 2017;6:156–60.
Chen L, Wilson R, Bennett E, Zosky GR. Identification of vitamin D sensitive pathways during lung development. Respir Res. 2016;17:47.
Nymark P, Lindholm PM, Korpela MV, Lahti L, Ruosaari S, Kaski S, Hollmen J, Anttila S, Kinnula VL, Knuutila S. Gene expression profiles in asbestos-exposed epithelial and mesothelial lung cell lines. BMC Genomics. 2007;8:62.
Via M, De Giacomo A, Corvol H, Eng C, Seibold MA, Gillett C, Galanter J, Sen S, Tcheurekdjian H, Chapela R, et al. The role of LTA4H and ALOX5AP genes in the risk for asthma in Latinos. Clin Exp Allergy. 2010;40:582–9.
Tulah AS, Parker SG, Moffatt MF, Wardlaw AJ, Connolly MJ, Sayers I. The role of ALOX5AP, LTA4H and LTB4R polymorphisms in determining baseline lung function and COPD susceptibility in UK smokers. BMC Med Genet. 2011;12:173.
Ghelfi E, Karaaslan C, Berkelhamer S, Akar S, Kozakewich H, Cataltepe S. Fatty acid-binding proteins and peribronchial angiogenesis in bronchopulmonary dysplasia. Am J Respir Cell Mol Biol. 2011;45:550–6.
Gally F, Chu HW, Bowler RP. Cigarette smoke decreases airway epithelial FABP5 expression and promotes Pseudomonas aeruginosa infection. PLoS One. 2013;8:e51784.
Bein K, Wesselkamper SC, Liu X, Dietsch M, Majumder N, Concel VJ, Medvedovic M, Sartor MA, Henning LN, Venditto C, et al. Surfactant-associated protein B is critical to survival in nickel-induced injury in mice. Am J Respir Cell Mol Biol. 2009;41:226–36.
Availability of data and materials
The dataset supporting the conclusions of this article is included within the article.
Ethics approval and consent to participate
This study was approved by the ethical committee of Chinese PLA General Hospital.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Xie, X., Zhao, J., Xie, L. et al. Identification of differentially expressed proteins in the injured lung from zinc chloride smoke inhalation based on proteomics analysis. Respir Res 20, 36 (2019). https://doi.org/10.1186/s12931-019-0995-0
- Smoke inhalation
- Lung injury
- Differentially expressed proteins
- Signalling pathways