- Research
- Open Access
- Published:
Proteomics analysis of chronic skin injuries caused by mustard gas
BMC Medical Genomics volume 15, Article number: 175 (2022)
Abstract
Sulfur mustard (SM) is an alkylating and forming chemical that was widely used by Iraqi forces during the Iran–Iraq wars. One of the target organs of SM is the skin. Understanding the mechanisms involved in the pathogenesis of SM may help better identify complications and find appropriate treatments. The current study collected ten SM-exposed patients with long-term skin complications and ten healthy individuals. Proteomics experiments were performed using the high-efficiency TMT10X method to evaluate the skin protein profile, and statistical bioinformatics methods were used to identify the differentially expressed proteins. One hundred twenty-nine proteins had different expressions between the two groups. Of these 129 proteins, 94 proteins had increased expression in veterans' skins, while the remaining 35 had decreased expression. The hub genes included RPS15, ACTN1, FLNA, HP, SDHC, and RPL29, and three modules were extracted from the PPI network analysis. Skin SM exposure can lead to oxidative stress, inflammation, apoptosis, and cell proliferation.
Introduction
Sulfur mustard (SM) could be a chemical operations agent that causes blisters on the skin and mucous membranes [1]. Quite 90% of patients exposed to dichloroethyl sulfide showed various skin lesions within the damaged areas. Skin may absorb about 20% of this blistering agent during exposure. About 70% of chemicals can accumulate within the epidermis and the rest within the basement membrane and dermis [2]. SM appears to wreck the skin by disrupting cell proliferation [3]. The Skin symptoms of SM exposure include itching, burning, painless erythema or sunburn, hypopigmentation, hyperpigmentation, and melanoma in exposed and unexposed areas [4, 5]. The mechanism of SM skin symptoms is not fully understood. Acute skin symptoms caused by SM can emerge as itching, erythema, blistering and burning sensation reckoning on the dose and duration of exposure [5, 6]. 2 to 24 h after exposure to SM, only a tiny low number of basal cells show nuclear changes like chromatin loss. Hydropic degeneration of the cytoplasm, and liquefactive necrosis of the epithelium has been reported in acute exposure [7]. The blisters usually appear as small vesicles after about 17 h of exposure to SM, but usually heal within three weeks [8]. In mild cases, the skin lesions could also be limited to erythema that darkens in about 15 days, while the superficial layers of the epidermis become scaly without causing an actual skin defect [9]. Histological studies of skin exposed to SM show vasodilation and neutrophil leakage, indicating the assembly of vasoactive mediators and chemical adsorbents within the damaged areas [10]. Interleukins IL-1α/β, IL-6, IL-8, and TNF-α are released shortly after exposure to SM [11]. The NF-қB pathway and mitogen-activated protein kinases are significantly involved in the the regulation of genes encoding inflammatory cytokines after SM injury [12]. Increased formation of reactive oxygen species and decreased glutathione production occur after exposure to SM [13]. Various metabolites are formed following the reaction of SM with glutathione, which may be detected within the urine [14]. NO signaling pathways also are important in modulating inflammation and necrobiosis in sulfur mustard-exposed skin cells [15]. Although vast amounts of knowledge are obtained in previous research to know the pathophysiology of SM toxicity, many questions remain and more research is needed to answer them [16]. Our main goal during this study was to gauge the expression of different skin proteins in chemical veterans exposed to SM and to search out a possible link between common skin problems including pruritus, eczema, atopic eczema, hyperpigmentation and hypopigmentation and melanoma with changes the expression of skin proteins among these patients.
Materials and methods
Statement
The research ethics committee of Baqiyatallah University of Medical Sciences has approved the project (Approval ID: IR.BMSU.REC.1398.114). Informed consent was obtained from all subjects, and all methods were carried out in accordance with the relevant guidelines and regulations of REC. All experiments and relevant protocols were performed in accordance with the relevant guidelines and regulations of Baqiyatallah University of Medical Sciences.
Patients and controls
After filling in the questionnaire and informed consent, ten patients (SM, moderate to severe) and ten controls who were matched in sex and age were included in the study. Sample inclusion and exclusion criteria are stated in Table 1. The patients had medical records and their mustard injuries were confirmed by expert Dermatologists. Skin biopsy was done the same dermatologist from the inner part of the arm near the axillary fossa. The size of sampling was about 2 * 1 cm2 and the depth of it was hypodermic fat. After cleaning of the sample from blood, it was inserted it in normal saline and then prepared for laboratory examination.
Sample preparation
Skin samples were resuspended in 800 μl of 1% sodium deoxycholate in 100 mM NH4HCO3 was added to the samples and vortexed vigorously; the samples were centrifuged at 5000 RPM. The protein solution was lyophilized and covered with parafilm. According to the manufacturer’s instructions, protein concentrations were determined using the BCA assay (CAT NO:23235 Thermo Scientific™). The disulfide bond of cysteine in the protein was reduced with 10 mM dithiothreitol (DTT) at 37 °C for 1 h, followed by alkylation with 20 mM Ioadoacetamid (IAA) for 45 min in a dark place with room temperature. The remaining IAA in the sample was quenched with 10 mM DTT for 15 min in a dark with room temperature. Proteolysis of samples using Lys-C protease (protein to enzyme ratio 100:1) was performed overnight at 28 °C, followed by trypsin digestion (protein to enzyme ratio 100:1) for 6 h at 37 °C. The pH of the sample was adjusted to approximately 3 using a final concentration of approximately 1% trifluoroacetic acid (TFA) and then desalted using a solid phase extraction disc containing styrene–divinylbenzene, and stepped tips (Empore SDBRPS 47 mm Extraction Disc, SUPLCO). Stage tips were self-packed into pipette tips, peptides were bound to the stage tips, washed with 0.2 percent TFA, and lastly eluted with 80 percent acetonitrile: 5% ammonium hydroxide. The peptides were vacuum centrifuged and then reconstituted in 200 mM HEPES pH 8.8 before being quantified using the Pierce quantitative colorimetric peptide assay (CAT NO:23235 Thermo Scientific™).
Sample labeling
Peptides were labeled with tandem mass tag (TMT) reagent according to the manufacturer's instructions (CAT NO: 90110 Thermo Scientific™). Briefly, anhydrous acetonitrile was added to each vial labeled with TMT, followed by shaking for 5 min and short centrifugation. Aliquots of individual peptide samples were labeled with one of the individual TMT tags (10 tags total). Labeling was performed at room temperature for 1 h with occasional shaking.
To quench excess TMT label in the samples, 5% hydroxylamine was added to each sample, shaken, and incubated at room temperature for 15 min. A label check experiment was performed by mixing 1.5 μl of each individually labeled vacuum-dried TMT sample to ensure that the same amount of total peptide was pooled in all samples prior to sample pooling. Samples were reconstituted in 2 μN, 0.1 μl water and analyzed by LC connected to a mass spectrometer (QExactive, Thermo Fisher, USA). Normalization factors were obtained from label validation experiments, and samples of TMT-labeled peptides were pooled in a 1:1 ratio across all samples and dried under vacuum. Samples were purified by C18 solid phase extraction desalting (SPE, SepPak, Waters) and centrifuged under a vacuum to dryness. The peptide mixture was fractionated into 96 fractions by high pH reverse phase HPLC (HpH), which were then pooled into 17 fractions before LCMS/MS analysis. HpH HPLC fractions of each TMT kit were reconstituted with sample loading buffer (2% acetonitrile, 97.9% water, 0.1% formic acid) and subjected to LC–MS/MS analysis.
Preparation of the samples for MS analysis
The HpH HPLC fractions of each TMT set were reconstituted with sample loading buffer (2% acetonitrile, 97.9% water, 0.1% formic acid) and subjected to LC–MS/MS analysis.
LC–MS/MS analysis
1D data dependent acquisition (DDA) of peptides on QExactive Quadrupole-Orbitrap (QE-classic)
TMT-labeled peptide samples were injected into an in-house packed trap column and desalted with loading buffer. The peptide was eluted from the trap into an in-house packed analytical column with linear gradients of mobile phases A and B (mobile phase B (30%) over 110 min at a flow rate of 300 nL/min across the gradient). The eluent from the trap was separated on the analysis column. The column eluate was fed to the ion source of the mass spectrometer. An electrospray voltage of 2.6 kV was applied through the fluid connection upstream of the column. Peptide precursors from 350 to 1850 m/z were scanned at 70 k resolution using a 1 × 106 target AGC. The 10 strongest ions in the previous survey scan were fragmented by high energy collision dissociation (HCD) using a collision energy of 35 normalized with a separation width of 0.7 m/z. For MS/MS analysis, only precursors with charge levels of + 2 to + 4. In the MS procedure, the minimum signal required for the MS2 trigger was 2.5 × 104, the AGC target value for MS2 was 2 × 105, and the maximum injection time for MS2 was 250 ms. The MS/MS scan resolution was set to 70 k. The dynamic exclusion was set to 90 s.
Protein identification and quantification
The mass spectrometric data files for each sample set were searched using Proteome Discoverer (version 2.1, Thermo Scientific). Uniprot database (181005_UniPr_HUMAN_Revi + Unrevi.fasta, Ref: http://www.uniprot.org) containing 95,106 human proteins including isoforms and unreviewed (Homo sapiens) was used for searching the data. The quantitative ratios were generated using the quantitative values found in channel 126 (control) as the denominator.
Data processing
TMT quantitative proteomic analysis of human samples resulted in the identification of 3586 high confident (Protein, Peptide and PSM FDR < 1%) proteins in skin samples. Density plot and Box Plot analysis of normalized data (Fig. 1) show similar protein ratio distribution (ratio with reference to channel 126, a control sample) across all samples within each respective data set, therefore, a relative abundance comparison was performed. Normalization is accomplished by dividing abundance values by the total intensity of the provided array (i.e., the sum of all abundance values).
Statistical analysis
The differentially abundant proteins were analyzed using the Limma Packages [17], written in R. The Bioconductor software was used to analyze the data, which fits a linear model to each protein and then uses the Bayes technique to the predicted variances to increase power. The results were validated by the ANOVA test separately. Comparison of normalized protein areas between control and patient samples showed a total of 129 differentially abundant proteins (Fold change ≥ 1.2 and Fold change ≤ 0.8, statistically significant differences between patient and control were accepted as p-value < 0.05) in skin samples. The reason for choosing this Fold change was due to the chronic nature of the disease. Selecting a higher threshold makes it difficult to find the altered proteins for patients with chronic disease. Moreover, regarding the limited number of samples, the adaptation of Adjusted p-value for a low number of samples limits the validation of significant proteins.
Gene ontology and pathway enrichment
Gene Ontology (GO) was used to define gene functions in biological process (BP), molecular function (MF) and cellular component (CC) aspects. Enrichr is an online enrichment analysis tools that provides a functional annotation set to understand the biological meaning behind a large list of genes or proteins [18]. The functional enrichment analysis of the statistically significant proteins, including GO and KEGG pathway enrichment analysis [19], was conducted using Enrichr.
PPI network and network analysis
The list of proteins was analyzed using the online website STRING (https://string-db.org/, version 11.5) for protein–protein interaction analysis [20]. Then, the software Cytoscape was used to establish a PPI network. The Network Analyzer in Cytoscape was utilized to calculate node degree [21]. Also, the MCODE plugin was used to perform modular analysis [22], with the parameters set as follows: a Degree Cutoff = 2, Node Score Cutoff = 0.2, K-Core = 2, and Max. Depth = 100.
Results
Identification of significant abundance proteins
This study evaluated the protein profiles of 10 healthy individuals and 10 chemical warfare veterans with a history of SM exposure were evaluated. A total of 129 proteins were identified with significant increasing and decreasing regulatory changes, including 35 proteins with down-regulation and 94 proteins with up-regulation (Additional file 1). Table 2 shows the top 10 down regulated skin proteins and Table 3 show the top 10 up-regulated skin proteins.
Among the 129 proteins in this study that were associated with changes in statistically significant expression, their role and importance in the occurrence of some skin problems such as pruritus, inflammation, eczema, and melanoma observed in chemical warfare exposed to sulfur SM are different.
GO and pathway enrichment analyses
The list of GO analysis and the KEGG pathway are shown in Additional file 2. Figure 2 shows the top 10 enrichment analysis for all significant abundance proteins. The top 10 KEGG pathway of significant abundance proteins are found in Fig. 2a. In terms of BP, the top 10 significant abundance proteins are found in Fig. 2b, As far as MF is found in Fig. 2c. CC is found in Fig. 2d.
PPI network and identification of hub genes
A PPI network of the proteins was constructed using the STRING and Cytoscape tool. The PPI network contained 129 nodes, including 94 upregulated genes and 35 downregulated genes, and 121 edges; (Fig. 3a). The term “degree” in the PPI network means a number of interactions between two genes (nodes). The hub genes of the PPI network were filtered with a cut-off value of degree > 6. As a result, six hub genes were identified. The hub genes included RPS15, ACTN1, FLNA, HP, SDHC, and RPL29. Three modules were extracted from the PPI network, including module 3 with five nodes (i.e., NDUFA6, SDHC, DUFA12, COX6C, and UQCR10) (Fig. 3b), module 2 with five nodes included HBG1, HBB, HBD, HBA1, and HP (Fig. 3c), and module 1 with six nodes of TPM4, MYH11, ACTN1, TPM2, FLNA, and TPM1 (Fig. 3d).
Discussion
Proteomic analysis of skin biopsy samples could aid in the discovery of potential disease biomarkers, as well as provide a better knowledge of disease pathophysiology. In the current study, TMT-based proteomics technologies on the skin of SM-exposed patients were used to show a different pattern compared to healthy controls. This is the first study to evaluate the skin's global proteome profiles in these patients in comparison with the control group.
Our findings showed six hub proteins (RPS15, ACTN1, FLNA, HP, SDHC, and RPL29) that three modules were extracted from the PPI network. Module 1 with five nodes includes NDUFA6, SDHC, DUFA12, COX6C, and UQCR10 (KEGG Pathway Enrichment: Oxidative phosphorylation, p-value: 9.3E-09 and FDR: 1.12E-07), module 2 with five nodes includes HBG1, HBB, HBD, HBA1, and HP, and module 3 with six nodes includes TPM4, MYH11, ACTN1, TPM2, FLNA, and TPM1 (KEGG Pathway Enrichment: Regulation of actin cytoskeleton, p-value: 0.001 and FDR: 0.004) were recovered from the PPI network.
Protein nodes in the first module are associated with mitochondrial electron transport chain (ETC) that includes 4 mitochondrial complexes. Significant alterations in the protein expression of mitochondrial respiratory chain complexes I, II, III, and IV in the SM-exposed patients compared to healthy controls are in line with previous studies [23,24,25,26]. There have been reports of mitochondrial dysfunction as a result of SM exposure [25, 26]. For example, Sourdeval et al. found that SM causes mitochondrial malfunction, which is followed by an increase in reactive oxygen species (ROS) formation and cellular OS [25]. Similar results were reported by Gould et al. for 2-chloroethyl ethylsulfide (CEES), and half mustard [26]. Another study found that mustards increase ROS production by changing cytochrome P450 activity by uncoupling microsomal electron transport at the flavoenzyme reductase [27]. Another mustard derivative, 11β: aniline mustard, alters the ETC at complex I, increasing ROS and OS generation while inhibiting oxygen consumption [28]. Another study reported chemical casualties had a higher pro-oxidant–antioxidant balance (PAB) score in comparison with the healthy control, indicating that they were exposed to OS [29].
ROS, particularly mitochondrial-derived ROS, play a role in a variety of signaling pathways, including eliciting antioxidative signaling, initiating DNA damage responses, affecting iron homeostasis, stimulating apoptosis, and signaling for cell survival and proliferation [30,31,32]. Thus, mitochondria may act as reservoirs of therapeutic targets in SM-exposed patients.
Altered expression and activity of mitochondrial complexes and increased OS in skin diseases similar to those exposed to mustard, including atopic dermatitis and atopic eczema, have also been reported. Previous research has found that nonlesional atopic dermatitis (ADNL) skin had higher expression of genes associated to OS [33]. Leman et al. discovered that increased complexes I and II activity in ADNL keratinocytes (KCs). Enhanced tricarboxylic acid cycle (TCA) turnover was discovered through metabolomics analysis. Mitochondria are important regulators of reactive oxygen species (ROS) production. In ADNL KCs, increased aerobic metabolism resulted in OS. ADNL human epidermal equivalents had a higher mitochondrial function and an increased OS response when compared with the controls. In addition, they discovered that glycolysis in ADNL KCs complements but does not replace mitochondrial metabolism. As a result, aerobic metabolism dominates in ADNL, causing OS [34].
Of note, patients participating in this study, in addition to skin complications caused by SM, also have respiratory problems and hypoxia. Hypoxia affects the mitochondrial shape, such as cristae structure, as well as protein expression, such as HIF proteins [35, 36]. HIF proteins regulate intermediate metabolism by inducing many target genes, such as raising the expression of glycolysis-related proteins and lowering oxygen-dependent pathways via modifying the ETC complex expression and activity [35] that is consistent with enhanced glycolysis and decreased tricarboxylic acid cycle (TCA) turnover was discovered through metabolomics analysis in SM-exposed patients [37, 38]. Subsequently, HIF-1 cooperates with STAT3 on the Haptoglobin (Hp) promoter to increase Hp gene expression (node in the second module). Hp is a significant acute phase protein that binds to extravascular hemoglobin (Hb) and eliminates the harmful Hb quickly via macrophages' CD163 scavenger receptor [39]. As a result, Hp maintain tissues from oxidative damage caused by Hb [40, 41]. Increased level of Hp was shown in plasma SM-exposed patients in comparison with the control group [42]. Patients with skin illnesses, such as psoriasis, have considerably higher Hp mRNA expression in epidermal keratinocytes than controls [43]. Kkeratinocyte-derived Hp may be involved in the skin's downregulation of inflammatory reactions.
In addition to ROS production, mitochondria also are involved within the function of melanocytes and pigments [44]. Hair disorders and pimples are a part of the wide selection of symptoms of mitochondrial disease. Skin problems including scaling, itching and erythema are reported in some mitochondrial diseases like mitochondrial encephalomyopathy and lactic acidosis [45].
Nodes (Tropomyosin 1, Tropomyosin 2, Tropomyosin 4, Myosin Heavy Chain 11, Filamin A and Actinin Alpha 1) in module 3 are important in the protection of the cellular cytoskeleton that all of them increased in SM-exposed patients compared with the control group. Studies showed that in the short term (several hours after SM exposure), changes in these proteins could cause blisters [46, 47], while increased level of these proteins several years after SM exposure is unclear. It can be due to ROS-induced apoptosis or cell proliferation [35].
In addition to the proteins of three modules mentioned above, some of the top 10 proteins were described in the following part. ISG15 proteins are induced by interferon type I and have many functions. It also acts as an extracellular cytokine and modifier of intracellular proteins [48]. Non-conjugated extracellular ISG15 acts as a cytokine to modulate immune responses [49, 50]. Inflammation-induced interferon type 1 can even appear as ulcerative skin lesions within the armpits, groin, and neck [51]. ISG15 appears to act as a carcinogenic protein as well as a tumor suppressor protein [52]. The presence of dermatitis and melanoma of the skin may not be unrelated to elevated levels of this protein. Previous studies have reported that TSR2 overexpression can suppress NF-қB transcription activity in Hep-2 cells [53]. Suppression of NF-κB induces apoptosis [54, 55]. In our study, a rise within the level of Pre-rRNA-processing protein TSR2 homolog was observed within the skin of chemical operations victims. Given the role that these proteins play in carcinogenesis and apoptosis, it will be expected that these proteins also are involved within the development of melanoma. RPS15 proteins, independent of their major role in ribosomal assembly and protein translation, also accomplish several extra-ribosomal functions, including regulating apoptosis, cell cycle arrest, cell proliferation and migration, invasion, and repair of DNA [56]. RPS15 may be associated with expression disorders in some sorts of cancer, including esophageal, lung and skin cancers [57,58,59]. In an exceedingly, Wang et al. indicated that RPS15 is over-expressed in gastric cancer tissues [60]. In our study, a rise in PRS15 ribosomal proteins was observed in chemical warfare victims exposed to SM. COX6C is found within the inner membrane of the mitochondria and is vital within the process of apoptosis [61]. Previous studies have shown that COX6C can cause chronic kidney disease (CKD), diabetes, carcinoma, melanoma and other diseases [62, 63]. Increased expression of COX6C in various invasive organs has been reported to always be observed in damaged tissue cells [64]. In our study, the level of COX6C proteins was significantly increased. Activation of P2RX7 stimulates several pathways, including the discharge of pro-inflammatory cytokines like IL-1 and IL-6, modulation of cell surface receptors, ROS and RNS formation and necrobiosis [65]. Immuno-histochemical studies indicate that P2RX7 within the upper layer of human skin is involved in the death of differentiated keratinocytes [66]. P2RX7 in Allergic dermatitis in humans is related to up-regulation within the basal epidermal layer of inflamed skin of patients with atopic eczema [67]. Increased P2RX7 levels are reported in human melanoma and in several melanoma cell lines [68]. In our study, we also observed a rise in the extent of P2RX7 protein within the skin of chemical operations victims. GIPC1 is involved in the transport of assorted membrane proteins and regulates various cellular processes including cell proliferation, cell surface polarity, cytokines and migration [69].
WWC3 proteins are involved in the regulation of cell proliferation and migration [70]. WWC3 proteins preventing the proliferation of malignant cells are capable of inducing the expression of apoptotic molecules. Within the presence of WWC3, the activity of caspase-3 and caspase-7 increases and WWC3 proteins act as an autophagy regulator in Non-Small Cell carcinoma. [71]. Down-regulation of WWC3, which is related to increased cell proliferation and high metastasis [72]. In our study, WWC3 proteins within the skin of chemical operations victims were associated with a significant down-regulation. MGMT is expressed as a repairing enzyme in any kind of human tissue, the quantity of which is expressed in several tissues. For instance, It has the highest expression within the liver; Excessive expression in epidermis, fibroblasts and melanocytes; It has relatively high within the lungs, kidneys and colon, and lowest within the pancreas, hematopoietic cells, lymph tissues, and brain [73, 74]. In animal models, increasing MGMT levels is to guard against the event of tumors including carcinoma [75, 76]. Therefore, a significant down-regulating MGMT protein in chemical operations victims is expected to be involved in the development of carcinoma. Delta-5-desachurase (FADS1) and delta-6-desachurase (FADS2) are two enzymes, which accountable for the synthesis of highly unsaturated fatty acids Linolic acid (LA) and Omega-6 fatty acid. Among the two essential fatty acids present within the skin, LA plays a crucial role in regulating the function of the skin barrier [77]. Atopic eczema is related to abnormal fatty acid metabolism, Because linolic acid (LA) deficiency in humans and rodents ends up in certain abnormalities within the skin of patients with atopic eczema [78]. In patients with atopic dermatitis, changes within the regulations of some transcription factors like SREBP-1c and PPAR-a have also been observed [79]. These transcription factors play a serious role in the regulation of FADS1 and FADS2 expressions [78]. In our study, FADS2 was related to a significant down-regulation within the skin of chemical operations victims.
Conclusion
According to what has been observed from skin protein changes in chemical veterans exposed to SM, it can be concluded that the increase and decrease in expression of some proteins are significant with skin problems such as itching, scaling, eczema, and atopic dermatitis. Distinct proteins are associated with OS, inflammatory pathways, apoptosis, and cell proliferation.
Availability of data and materials
The datasets generated during and/or analysed during the current study are not publicly available due to confidentiality agreements, supporting data can only be made available to bona fide researchers subject to a non-disclosure agreement but are available from the corresponding author on reasonable request.
References
Evison D, Hinsley D, Rice P. Chemical weapons. BMJ. 2002;324(7333):332–5.
Smith KJ, et al. Sulfur mustard: its continuing threat as a chemical warfare agent, the cutaneous lesions induced, progress in understanding its mechanism of action, its long-term health effects, and new developments for protection and therapy. J Am Acad Dermatol. 1995;32(5):765–76.
Pechura CM, Rall DP. Dermatological effects of mustard agents and lewisite, in veterans at risk: the health effects of mustard gas and lewisite. National Academies Press (US); 1993.
Balali-Mood M, Hefazi M. The pharmacology, toxicology, and medical treatment of sulphur mustard poisoning. Fundam Clin Pharmacol. 2005;19(3):297–315.
Momeni A-Z, et al. Skin manifestations of mustard gas: a clinical study of 535 patients exposed to mustard gas. Arch Dermatol. 1992;128(6):775–80.
Smith W, Dunn M. Medical defense against blistering chemical warfare agents (Reannouncement with new availability information). Army Medical Research Inst. of Chemical Defense, Aberdeen Proving Ground; 1991.
Pirie A. The action of mustard gas on ox cornea collagen. Biochem J. 1947;41(2):185.
Balali-Mood M, Hefazi M. Comparison of early and late toxic effects of sulfur mustard in Iranian veterans. Basic Clin Pharmacol Toxicol. 2006;99(4):273–82.
Emadi SN, et al. External urethral stenosis: a latent effect of sulfur mustard two decades post-exposure. Wiley Online Library; 2009.
Smith HW, Clowes GH, Marshall E. On dichloroethyl sulfide (mustard gas). IV. The mechanism of absorption by the skin. J Pharmacol Exp Ther. 1919;13:1–30.
Sabourin CL, Petrali JP, Casillas RP. Alterations in inflammatory cytokine gene expression in sulfur mustard-exposed mouse skin. J Biochem Mol Toxicol. 2000;14(6):291–302.
Rebholz B, et al. Role of NF-κB/RelA and MAPK pathways in keratinocytes in response to sulfur mustard. J Investig Dermatol. 2008;128(7):1626–32.
Kumar O, Sugendran K, Vijayaraghavan R. Protective effect of various antioxidants on the toxicity of sulphur mustard administered to mice by inhalation or percutaneous routes. Chem Biol Interact. 2001;134(1):1–12.
Black R, Read RW. Biological fate of sulphur mustard, 1,1′-thiobis (2-chloroethane): identification of β-lyase metabolites and hydrolysis products in human urine. Xenobiotica. 1995;25(2):167–73.
Feister AJ. Medical defense against mustard gas: toxic mechanisms and pharmacological implications. 1991.
Kehe K, et al. Molecular toxicology of sulfur mustard-induced cutaneous inflammation and blistering. Toxicology. 2009;263(1):12–9.
van Ooijen MP, et al. Identification of differentially expressed peptides in high-throughput proteomics data. Brief Bioinform. 2018;19(5):971–81.
Xie Z, et al. Gene set knowledge discovery with enrichr. Curr Protoc. 2021;1(3): e90.
Kanehisa M, et al. KEGG: integrating viruses and cellular organisms. Nucleic Acids Res. 2021;49(D1):D545–51.
Szklarczyk D, et al. The STRING database in 2021: customizable protein–protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021;49(D1):D605–12.
Shannon P, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
Bader GD, Hogue CWV. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinform. 2003;4(1):1–27.
Laskin JD, et al. Oxidants and antioxidants in sulfur mustard–induced injury. Ann N Y Acad Sci. 2010;1203(1):92–100.
Sabnam S, et al. CEES-induced ROS accumulation regulates mitochondrial complications and inflammatory response in keratinocytes. Chem Biol Interact. 2020;321: 109031.
Sourdeval M, et al. Inhibition of caspase-dependent mitochondrial permeability transition protects airway epithelial cells against mustard-induced apoptosis. Apoptosis. 2006;11(9):1545–59.
Gould NS, White CW, Day BJ. A role for mitochondrial oxidative stress in sulfur mustard analog 2-chloroethyl ethyl sulfide-induced lung cell injury and antioxidant protection. J Pharmacol Exp Therap. 2009;328(3):732–9.
Brimfield A, et al. Free radical production from the interaction of 2-chloroethyl vesicants (mustard gas) with pyridine nucleotide-driven flavoprotein electron transport systems. Toxicol Appl Pharmacol. 2009;234(1):128–34.
Fedeles BI, et al. Chemical genetics analysis of an aniline mustard anticancer agent reveals complex I of the electron transport chain as a target. J Biol Chem. 2011;286(39):33910–20.
Nobakht M. Gh BF, et al. Pro-oxidant–antioxidant balance in Iranian veterans with sulfur mustard toxicity and different levels of pulmonary disorders. Drug Chem Toxicol. 2016;39(4):362–6.
Li Z, et al. Roles of reactive oxygen species in cell signaling pathways and immune responses to viral infections. Adv Virol. 2017;162(3):603–10.
Redza-Dutordoir M, Averill-Bates DA. Activation of apoptosis signalling pathways by reactive oxygen species. Biochimica et Biophysica Acta -Molecular Cell Research. 2016;1863(12):2977–92.
Hanson ES, Leibold EA. Regulation of the iron regulatory proteins by reactive nitrogen and oxygen species. Gene Expres J Liver Res. 1999;7(4–5):367–76.
Blunder S, et al. Enhanced expression of genes related to xenobiotic metabolism in the skin of patients with atopic dermatitis but not with ichthyosis vulgaris. J Investig Dermatol. 2018;138(1):98–108.
Leman G, et al. Mitochondrial activity is upregulated in nonlesional atopic dermatitis and amenable to therapeutic intervention. J Investig Dermatol. 2022.
Fuhrmann DC, Brüne B. Mitochondrial composition and function under the control of hypoxia. Redox Biol. 2017;12:208–15.
Plecitá-Hlavatá L, Ježek P. Integration of superoxide formation and cristae morphology for mitochondrial redox signaling. Int J Biochem Cell Biol. 2016;80:31–50.
Nobakht BF, et al. NMR spectroscopy-based metabolomic study of serum in sulfur mustard exposed patients with lung disease. Biomarkers. 2017;22(5):413–9.
Ghoochani BFNM, et al. Metabolomics diagnostic approach to mustard airway diseases: a preliminary study. Iran J Basic Med Sci. 2018;21(1):59.
Nielsen MJ, Moestrup SK. Receptor targeting of hemoglobin mediated by the haptoglobins: roles beyond heme scavenging. Blood J Am Soc Hematol. 2009;114(4):764–71.
Lim YK, et al. Haptoglobin reduces renal oxidative DNA and tissue damage during phenylhydrazine-induced hemolysis. Kidney Int. 2000;58(3):1033–44.
Lim S-K, et al. Increased susceptibility in Hp knockout mice during acute hemolysis. Blood. 1998;92(6):1870–7.
Mehrani H, et al. Plasma proteomic profile of sulfur mustard exposed lung diseases patients using 2-dimensional gel electrophoresis. Clin Proteom. 2010;8(1):1–11.
Li P, et al. Localization of haptoglobin in normal human skin and some skin diseases. Int J Dermatol. 2005;44(4):280–4.
Ni-Komatsu L, Orlow SJ. Identification of novel pigmentation modulators by chemical genetic screening. J Investig Dermatol. 2007;127(7):1585–92.
Bodemer C, et al. Hair and skin disorders as signs of mitochondrial disease. Pediatrics. 1999;103(2):428–33.
Mol MA, van den Berg RM, Benschop HP. Proteomic assessment of sulfur mustard-induced protein adducts and other protein modifications in human epidermal keratinocytes. Toxicol Appl Pharmacol. 2008;230(1):97–108.
Mol M, et al. Proteomics as a strategy to study the mechanistic toxicology of sulfur mustard. In: Proceedings of the 2003 Meeting of NA TO TGO04, Medicin Hat, Canada. 2003.
Blomstrom DC, et al. Molecular characterization of the interferon-induced 15-kDa protein. Molecular cloning and nucleotide and amino acid sequence. J Biol Chem. 1986;261(19):8811–6.
Knight E, Cordova B. IFN-induced 15-kDa protein is released from human lymphocytes and monocytes. J Immunol. 1991;146(7):2280–4.
Guilbaud G, et al. Evidence for sequential and increasing activation of replication origins along replication timing gradients in the human genome. PLoS Comput Biol. 2011;7(12): e1002322.
Martin-Fernandez M, et al. Systemic type I IFN inflammation in human ISG15 deficiency leads to necrotizing skin lesions. Cell Rep. 2020;31(6): 107633.
Andersen J, Hassel B. The interferon regulated ubiquitin-like protein, ISG15, in tumorigenesis: friend or foe? Cytokine Growth Factor Rev. 2006;17(6):411–21.
Pasparakis M. Regulation of tissue homeostasis by NF-κB signalling: implications for inflammatory diseases. Nat Rev Immunol. 2009;9(11):778–88.
Morais C, et al. Pyrrolidine dithiocarbamate exerts anti-proliferative and pro-apoptotic effects in renal cell carcinoma cell lines. Nephrol Dial Transplant. 2006;21(12):3377–88.
Morais C, et al. Inhibition of nuclear factor kappa B attenuates tumour progression in an animal model of renal cell carcinoma. Nephrol Dial Transplant. 2010;25(5):1462–74.
Anderson SJ, et al. Ablation of ribosomal protein L22 selectively impairs αβ T cell development by activation of a p53-dependent checkpoint. Immunity. 2007;26(6):759–72.
Wang Q, et al. Cloning and characterization of full-length human ribosomal protein L15 cDNA which was overexpressed in esophageal cancer. Gene. 2001;263(1–2):205–9.
Dang C, et al. Identification of dysregulated genes in cutaneous squamous cell carcinoma. Oncol Rep. 2006;16(3):513–9.
Dmitriev AA, et al. Genetic and epigenetic analysis of non-small cell lung cancer with NotI-microarrays. Epigenetics. 2012;7(5):502–13.
Wang H, et al. Overexpression of ribosomal protein L15 is associated with cell proliferation in gastric cancer. BMC Cancer. 2006;6(1):1–8.
Belevich I, Verkhovsky MI. Molecular mechanism of proton translocation by cytochrome c oxidase. Antioxid Redox Signal. 2008;10(1):1–30.
Zaza G, et al. Downregulation of nuclear-encoded genes of oxidative metabolism in dialyzed chronic kidney disease patients. PLoS ONE. 2013;8(10): e77847.
Masha RT, Houreld NN, Abrahamse H. Low-intensity laser irradiation at 660 nm stimulates transcription of genes involved in the electron transport chain. Photomed Laser Surg. 2013;31(2):47–53.
Tian B-X, et al. Differential expression and clinical significance of COX6C in human diseases. Am J Transl Res. 2021;13(1):1.
Wiley J, et al. The human P2X7 receptor and its role in innate immunity. Tissue Antigens. 2011;78(5):321–32.
Greig AV, et al. Purinergic receptors are part of a signaling system for keratinocyte proliferation, differentiation, and apoptosis in human fetal epidermis. J Investig Dermatol. 2003;121(5):1145–9.
Pastore S, et al. Stimulation of purinergic receptors modulates chemokine expression in human keratinocytes. J Investig Dermatol. 2007;127(3):660–7.
White N, Butler PE, Burnstock G. Human melanomas express functional P2X 7 receptors. Cell Tissue Res. 2005;321(3):411–8.
De Vries L, et al. GIPC, a PDZ domain containing protein, interacts specifically with the C terminus of RGS-GAIP. Proc Natl Acad Sci. 1998;95(21):12340–5.
Wennmann DO, et al. Evolutionary and molecular facts link the WWC protein family to Hippo signaling. Mol Biol Evol. 2014;31(7):1710–23.
Han Q, et al. WW and C2 domain-containing protein-3 promoted EBSS-induced apoptosis through inhibiting autophagy in non-small cell lung cancer cells. J Thorac Dis. 2020;12(8):4205.
Hou J, Zhou J. WWC3 downregulation correlates with poor prognosis and inhibition of Hippo signaling in human gastric cancer. Onco Targets Ther. 2017;10:2931.
Pieper RO. Understanding and manipulating O6-methylguanine-DNA methyltransferase expression. Pharmacol Ther. 1997;74(3):285–97.
Kalal BS, Upadhya D, Pai VR, Chemotherapy resistance mechanisms in advanced skin cancer. Oncol Rev. 2017;11(1).
Zaidi NH, et al. Transgenic expression of human MGMT protects against azoxymethane-induced aberrant crypt foci and G to A mutations in the K-ras oncogene of mouse colon. Carcinogenesis. 1995;16(3):451–6.
Shiraishi A, Sakumi K, Sekiguchi M. Increased susceptibility to chemotherapeutic alkylating agents of mice deficient in DNA repair methyltransferase. Carcinogenesis. 2000;21(10):1879–83.
Hartop P, Prottey C. Changes in transepidermal water loss and the composition of epidermal lecithin after applications of pure fatty acid triglycerides to the skin of essential fatty acid-deficient rats. Br J Dermatol. 1976;95(3):255–64.
Nakamura MT, Nara TY. Structure, function, and dietary regulation of Δ6, Δ5, and Δ9 desaturases. Annu Rev Nutr. 2004;24:345–76.
Schlotter YM, et al. Altered expression of fatty acid desaturases in the skin of dogs with atopic dermatitis. J Dermatol Sci. 2009;54(1):49–52.
Acknowledgements
This research was supported by the Baqiyatallah University of Medical Sciences, Tehran, Iran. Aspects of this research have been facilitated by access to the Macquarie University trading as Australian Proteome Analysis Facility (APAF).
Funding
Not applicable.
Author information
Authors and Affiliations
Contributions
HB and MGH conceived of the presented idea. SHP and VJ carried out the experiment. ASH and SMِِD carried out data gathering. MA performed the computations. BFNMGH helped to revise the manuscript. All authors discussed the results and contributed to the final manuscript. All authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study involving human participants was in accordance with the ethical standards of the institutional and national research committee and with the 1964 Helsinki Declaration and its later amendments or comparable ethical standards. The study was approved by the Institutional Review Board and the Ethical Committee of Baqiyatallah University of Medical Sciences (Approval ID: IR.BMSU.REC.1398.114). Written informed consent was obtained from all study participants.
Consent for publication
Not applicable.
Competing interests
Authors have no conflict of interest to declare.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1
. The list of diffrentially abundant proteins.
Additional file 2
. The list of GO analysis and the KEGG pathway.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Jamshidi, V., Nobakht M. Gh, B.F., Parvin, S. et al. Proteomics analysis of chronic skin injuries caused by mustard gas. BMC Med Genomics 15, 175 (2022). https://doi.org/10.1186/s12920-022-01328-3
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12920-022-01328-3
Keywords
- Proteomics analysis
- Mustard gas
- Chemical injuries
- Systems biology
- Chemical agents