Next Article in Journal
Synthesis, Physicochemical Characterization using a Facile Validated HPLC Quantitation Analysis Method of 4-Chloro-phenylcarbamoyl-methyl Ciprofloxacin and Its Biological Investigations
Next Article in Special Issue
Computational Exploration of Minimum Energy Reaction Pathway of N2O Formation from Intermediate I of P450nor Using an Active Center Model
Previous Article in Journal
Prioritizing Endangered Species in Genome Sequencing: Conservation Genomics in Action with the First Platinum-Standard Reference-Quality Genome of the Critically Endangered European Mink Mustela lutreola L., 1761
Previous Article in Special Issue
Complement System Inhibitory Drugs in a Zebrafish (Danio rerio) Model: Computational Modeling
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 24
Issue 19
10.3390/ijms241914817
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

New Insights into Radio-Resistance Mechanism Revealed by (Phospho)Proteome Analysis of Deinococcus Radiodurans after Heavy Ion Irradiation

by
Shihao Liu
1,
Fei Wang
1,
Heye Chen
1,
Zhixiang Yang
2,
Yifan Ning
1,2,
Cheng Chang
1 and
Dong Yang
1,*
1
State Key Laboratory of Proteomics, Beijing Proteome Research Center, National Center for Protein Sciences (Beijing), Beijing Institute of Lifeomics, Beijing 102206, China
2
College of Life Sciences, Hebei University, Baoding 071002, China
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(19), 14817; https://doi.org/10.3390/ijms241914817
Submission received: 29 August 2023 / Revised: 23 September 2023 / Accepted: 27 September 2023 / Published: 1 October 2023
(This article belongs to the Special Issue Advanced Research in Prediction of Protein Structure and Function)
Browse Figures
Review Reports Versions Notes

Abstract

:
Deinococcus radiodurans (D. radiodurans) can tolerate various extreme environments including radiation. Protein phosphorylation plays an important role in radiation resistance mechanisms; however, there is currently a lack of systematic research on this topic in D. radiodurans. Based on label-free (phospho)proteomics, we explored the dynamic changes of D. radiodurans under various doses of heavy ion irradiation and at different time points. In total, 2359 proteins and 1110 high-confidence phosphosites were identified, of which 66% and 23% showed significant changes, respectively, with the majority being upregulated. The upregulated proteins at different states (different doses or time points) were distinct, indicating that the radio-resistance mechanism is dose- and stage-dependent. The protein phosphorylation level has a much higher upregulation than protein abundance, suggesting phosphorylation is more sensitive to irradiation. There were four distinct dynamic changing patterns of phosphorylation, most of which were inconsistent with protein levels. Further analysis revealed that pathways related to RNA metabolism and antioxidation were activated after irradiation, indicating their importance in radiation response. We also screened some key hub phosphoproteins and radiation-responsive kinases for further study. Overall, this study provides a landscape of the radiation-induced dynamic change of protein expression and phosphorylation, which provides a basis for subsequent functional and applied studies.

1. Introduction

D. radiodurans, a radiation-resistant bacterium, was discovered in 1956 in canned meat after exposure to ionizing radiation. Subsequent studies have shown that the bacteria are highly resistant to a variety of extreme environments including ionizing radiation, such as extreme desiccation, low/high temperatures, oxidative stress, and antibiotics [1,2,3]. When exposed to ionizing radiation, the cell wall structure of bacteria tends to be destroyed [4]. However, the SlpA protein and its unique Hpi protein in D. radiodurans form a thick S-layer, enhancing the rigidity of the cell wall and further improving its resistance ability [5,6]. The nucleoid region of D. radiodurans exhibits a ring-like structure and can still restore its normal morphology within hours after exposure to 15,000 Gy of ionizing radiation [7]. Different radiation qualities also show different biological effects. In human A549 cells, 2Gy of carbon ions caused the same survival ratio as 6Gy of X-rays [8]. This phenomenon seems to be more pronounced in D. radiodurans. In the face of heavy ion radiation, D. radiodurans demonstrates the phenomenon of “overkill” [9]. Under 80 Gy of 12C6+ heavy ion irradiation, the fatality rate of the bacterium reaches nearly half, then decreases with increasing doses, reaching a minimum at 160 Gy, and subsequently recovers. For comparison, the fatality rate of E. coli was up to 95% at 5 Gy and continued to be at a high level when the was dose increased [10].
Previous studies have discovered that D. radiodurans possesses DNA repair proteins which are similar to those found in most other prokaryotic organisms, including proteins in base excision repair (Ung, Mug, AlkA2, MutM, NTH, etc.), double strand break repair (UvrA, RecA, RecD2, RecO, RecR, RecF, RecQ, etc.), mismatch repair (MutS, MutL, UvrD, XseA, RecJ, SSB, DpoⅢ, Lig, etc.), nucleotide excision repair (UvrA, UvrB, UvrC, UvrD, Mfd, DpoⅠ, etc.) [11,12,13,14,15]. But, compared to radiation-sensitive species, these proteins in D. radiodurans may play some more effective roles in DNA repair. Deinococci-specific proteins, such as PprA, DdrA, DdrB, DdrC, DdrD, have also been confirmed to be highly expressed when exposed to irradiation, helping to maintain genome integrity [16,17,18]. Besides that, D. radiodurans also encodes two proteins named ‘DNA protection during starvation protein’ (Dps-1, Dps-2); both can store iron and manganese in ionic states and protect DNA from oxidative damage [19,20,21]. For the oxygen radicals that have already formed, the catalases in D. radiodurans and the abundant H-Mn2+ complex play a major role in their removal [1,22].
Post-translational modifications cannot only affect the structure, activity, stability, and cellular localization of a single protein, but also regulate the function and activity of a series of proteins. Succinylation proteomics studies revealed the regulatory role of lysine succinylation in D. radiodurans, indicating the important role of PTMs in its extreme resistance [23]. Protein phosphorylation (mainly on serine/threonine/tyrosine), as one of the most common and essential PTMs in eukaryotes, can rapidly regulate protein activity and achieve rapid adaptation to the environment [24,25]. Phosphorylation may also be widespread in bacteria and affect a variety of life activities, including radiation response. RqkA, a DNA damage responsive serine/threonine protein kinase, plays a pivotal role in DNA repair and cell cycle in D. radiodurans [26,27]. The phosphorylation of Tyr-77 and Thr-318 had been shown to regulate the activity of DrRecA, thereby increasing its affinity for dsDNA [28]. FtsZ protein is involved in bacterial cell division, and the phosphorylation of DrFtsZ at S235 and S335 causes the loss of its function and reduces the cell tolerance to the gamma ray [29]. However, a smaller number of phosphosites have been identified in bacteria compared to eukaryotes. Even in Escherichia coli, a common model organism, fewer than 1600 high-confidence class I phosphosites (localization probability ≥ 0.75) have been identified [30]. Only a few of them have functional annotations [31]. The profiling and in-depth functional analysis of the dynamic phosphoproteome of D. radiodurans will help us to understand the response mechanisms to high doses of radiation.
In this study, we used label-free proteomic technology to detect the dynamic changes of the (phospho)proteome of D. radiodurans under different heavy ion doses and different time points after irradiation. The dynamic changing patterns and functional interpretations were explored. The pathways and the key hub phosphoproteins, as well as kinases related to irradiation resistance were screened for further studies.

2. Results

2.1. Global Analysis for (Phospho)Proteome of D. radiodurans after Heavy Ion Irradiation

In order to detect the dynamic changes of protein levels and phosphorylation modification levels after D. radiodurans irradiation, we selected three radiation doses (20 Gy, 80 Gy, 160 Gy) [10] and two time points after irradiation (0h, 2h), with a control group at each time point. High-resolution LC-MS/MS was used for peptide identification and relatively accurate protein quantification was obtained (Figure 1A). We quantified 2359 proteins, more than all published data to date [10,32,33] (Figure 1B), with more than 2000 proteins identified per sample on average (Figure 1C, Table S1A). In addition, we identified 1110 class I phosphosites among 568 phosphoproteins (Table S1B), including 587 (53%) serine sites and 407 (37%) threonine sites (Figure 1C). Almost all phosphoproteins have been observed and quantified in the proteome (Figure S1). The violin plot shows the quantification, median, quartile, and expression distribution of each experiment, with no significant difference in each of the 16 measurements of the proteome and phosphoproteome, demonstrating the stability of our dataset (Figure S2A,B). Using PCA (principal component analysis) to assess the variability between different samples of the proteome, we found that samples of the same type cluster together, while samples at different time points or irradiated with different doses could be separated clearly, indicating the repeatable variation between different sample groups (Figure 1D). Spearman correlation analysis found that the correlations between the two repeats were 0.89–0.92, further indicating that the experiment had good repeatability (Figure S2C).
Of all the proteins identified, 1553 (~66%) were significantly changed after irradiation (two-way ANOVA, FDR < 0.05, Figure 1E, Table S1). The counts of upregulated proteins (FDR < 0.05, log2FC > 1) were generally higher than those of downregulated ones (FDR < 0.05, log2FC < 1), and this trend was more obvious with the increase in radiation doses and time points (Figure 2A, Figures S2 and S3). A total of 940 proteins were upregulated in at least one radiation state. Compared with three previous irradiation datasets of D. radiodurans, 33 upregulated proteins were identified in at least three datasets, 32 of which were identified in our data (Figure S5). In contrast, the phosphorylation level of 255 (~23%) sites significantly changed after irradiation, less than the proportion of differential proteins (Figure 1E, Table S1A,B). However, we found an interesting phenomenon that the upregulation of the phosphorylation level (median log2FC = 5.22) was much higher than the level of variation in protein abundance (median log2FC = 1.23) (Figure 1F). This means that a smaller number of differential phosphosites are sensitive response sites for irradiation in D. radiodurans. Their functional role in radiation resistance needs further in-depth analysis.

2.2. Dose and Stage Dependence of D. radiodurans Irradiation Response

We conducted a hierarchical clustering analysis of the 1553 differently expressed proteins (DEPs) after being exposed to irradiation (Figure 2B). The analysis revealed that for samples collected immediately after radiation (0h), there was a clear distinction between DEPs at different doses. However, in samples collected 2h after irradiation, this distinction was not apparent (Figure 2B). The number of upregulated proteins was generally higher than that of the downregulated ones, and this trend was more obvious with an increase in doses (Figure 2B, Figures S3 and S5). Statistics on the number of upregulated proteins at different doses also showed that the samples at 0h after irradiation were significantly differential at each dose, and 178 (23.5%) of the three doses were upregulated together (Figure 2C). Meanwhile, the number of co-upregulated proteins reached 240 (42%) in the samples collected 2h after irradiation (Figure 2D). These results reflect the characteristics of the early and late radiation responses of D. radiodurans.
Since the 0h samples have their own characteristics at each dose, function enrichment analysis was performed for the proteins upregulated at each state alone and the proteins upregulated in all doses. Gene ontology (GO) analysis showed that 20 Gy-specific upregulated proteins tend to participate in the cytokinin biosynthetic process (e.g., DR_0471) and trehalose biosynthetic process (e.g., DR_0464, TreZ). TreZ is a key protein in trehalose synthesis in bacteria [34], and trehalose acts as a protector against various stress conditions [35]. In addition, heat shock protein binding entries (e.g., DnaJ) were also enriched. With the increase in irradiation dose, the expressions of proteins related to transmembrane transport and tetrapyrrole metabolism were upregulated. Proteins related to ligase activity were upregulated at various radiation doses, including the proteins related to the ligation between tRNA and amino acids (e.g., GlyQS, PurD, GatC, ProS), and tRNA methylation (e.g., RtcB) [36,37,38,39]. These processes are related to accurate translation, which is important for cell repair after irradiation.
From the viewpoint of different stages after irradiation, the upregulated proteins at different time points showed obvious different functional characteristics. Proteins related to ATP-dependent DNA metabolism and phosphorylation were immediately upregulated after irradiation, which suggests the significance of phosphorylation for the D. radiodurans radiation response. It is noteworthy that membrane proteins and transmembrane transport proteins were upregulated during irradiation recovery (2h). Receptors on the cell membrane can sense external radiation stimuli and initiate a series of intracellular reactions through signal transduction pathways [40]. Transporters in the membrane help to maintain cell homeostasis and are involved in maintaining membrane integrity [41]. For example, LrgA and LrgB can form protein complexes that exert anti-apoptotic functions by maintaining the integrity of cell membranes [42]. Proteins related to signal transduction and nucleic metabolism were upregulated immediately after irradiation and remained in an upregulated state after 2h. Specifically, sustained upregulated proteins were enriched in the positive regulation of signal transduction, such as TOR signaling and nucleotide biosynthetic process. Protein domain enrichment analysis showed similar results to GO functional analysis (Figure S6).

2.3. KEGG Pathway Analysis of Differentially Expressed (Phospho)Proteins

The pathway data of D. radiodurans were downloaded from the KEGG website (https://www.kegg.jp) [43] and matched with the (phospho)proteome (Table S3). Most of the proteins (91.7%) in the pathways were quantitatively identified. Several proteins in the RNA polymerase (e.g., RpoC, RpoZ) and RNA degradation (e.g., Eno, Rny, GroEL) were upregulated, as well as the level of phosphorylation modification (Table 1). In the antioxidant pathways, such as alpha-linolenic acid metabolism, phenylalanine metabolism, and riboflavin metabolism, the levels of protein expression and phosphorylation were all upregulated to different extents (Figure 3A). These pathways may help to clear radiation-generated ROS and maintain cellular homeostasis [44,45,46]. Notably, 20 of the 29 identified proteins in the oxidative phosphorylation pathway were upregulated, including multiple subunits of NADH dehydrogenase and V-type ATP synthase, and cytochrome c oxidase (Table S3B). ATP produced by the mitochondrial respiratory chain is necessary for maintaining the energy supply of cells and for normal physiological functions and metabolic activities. The expressions of many proteins in the pathways of homologous recombination (RecFOR pathway) and mismatch repair were upregulated in at least one state after irradiation (Figure 3B). Due to the lack of RecB/C/D protein, the RecFOR pathway plays a key role in the recombination repair of D. radiodurans [47]. The key proteins RecF/R in this pathway were upregulated in all radiation states, with the highest fold change up to more than 35. The fold change of RecO in ‘160Gy, 0h’ also reached about 5. The activation of the core proteins in the mismatch repair pathway, MutS and MutL, helps chromosomes to reshape quickly to maintain genome stability.

2.4. Dynamic Patterns of Phosphorylation and Functional Interpretation

Preceding analyses revealed that the upregulation of phosphorylation is significantly higher than the upregulation of protein level, suggesting that phosphorylation may play an important role in the irradiation response. In order to better understand the distribution of 255 differential phosphosites, Fuzzy c-means clustering was used to divide them into four clusters, which showed distinct differences (Figure 4A). Cluster 1 (n = 56) phosphosites were upregulated in the early stage and recovered in the late stage, cluster 2 (n = 66) phosphosites were not significantly changed in the early stage and upregulated in the late stage, and cluster 3 (n = 68) phosphosites showed continuous upregulation, while cluster 4 (n = 65) phosphosites showed downregulation. Motif analysis showed that RS/RT motifs were enriched in different expression phosphosites (Figure S7). For these significantly changed phosphoproteins, we downloaded the protein–protein interaction (PPI) data from the STRING database and constructed a PPI network [48]. A total of 25 core regulatory proteins were screened (Figure 4B). DR_0183, with the highest score, has the glutamine amido transferases class II domain, participating in GABA biosynthesis, which has been reported to develop a significant response to radiation stress in plants [49,50,51]. In addition, a series of proteins, including RpoB, a core enzyme for RNA synthesis, have phosphosites with different dynamic patterns, which may be regulated by different kinases and play different parts in their functional activities.
According to their functions, differential phosphoproteins are mainly divided into nine classes, involved in different vital activities (Figure 4C). There are multiple metabolic proteins, which participate in glycometabolism, lipid metabolism, and the metabolism of other nutrients, and these account for a considerable proportion. In the proteins located in the cell wall, S460 in Deinococcus-specific HPI protein showed late upregulation. Another protein SlpA was upregulated late at T243, and the other three sites were downregulated (Y170, S241, T661). The phosphorylation of ribosomal proteins (RpsN, RpsL, RplB, RplY) and heat shock proteins (GroEL, GroES, DnaK) also showed upregulation. Previous studies have suggested that there is also a relationship between protein disorder and stress response [52]. Thus, we calculated the disorder degree for each amino acid residue and classified the proteins into different disorder categories (Table S2). FtsKL and DR-specific protein DdrD were identified as disordered proteins, and their upregulation sites (S905 in FtsKL, cluster 2; S170 in DdrD, cluster 3) are located in disordered regions. A further analysis revealed that the phosphosites identified by D. radiodurans tended to be located in disordered regions (Figure S8), similar to the situation in eukaryotes [53].

2.5. Comparative Analysis of Phosphosites and Protein Expression

To determine whether changes in phosphorylation are due to the changes in protein levels, we analyzed the consistency of the two omics using log2FC standard (|log2FC| ≤ 1 means no change). The proteins were classified using “both up”, “both down”, “both no change”, “pro > pho”, and “pho > pro” (Figure 5A). Each radiation state was calculated separately (Figure 5B). If one protein has multiple phosphosites, and the dynamic pattern of these sites and the protein level meet different classification criteria, the protein will be counted multiple times. We found that the two omics were inconsistent in the majority of cases, i.e., more than 50% in each irradiation state. The class of “both up” increased after 2 h irradiation, but still maintained a lower level. This suggests that radiation stress specifically regulates proteins primarily at a single omics level.
Functional analysis was performed for proteins inconsistent in the two omics levels and for proteins that were both upregulated in the two omics (Figure 5C). The “pro > pho” class was enriched in lipid metabolism and protein folding. The “pho > pro” class proteins were related to nucleic acid metabolism, translation, and the cellular response to an abiotic stimulus. These terms are more closely related to stress response and repair. For example, we were surprised to find that the protein level of RpoB did not change significantly in any state, but there were significant changes in its six phosphosites (Figure 5D). This suggests that phosphorylation is essential for its function after irradiation.

2.6. Irradiation-Sensitive Proteins Containing Kinase Domains

Studies on the kinome of D. radiodurans are lacking at present. In this study, we identified 79 proteins containing kinase domains, 52 of which were sensitive to irradiation, including 35 proteins that were upregulated in at least one state (Figure 6, Table S4). In the annotation of its domains and functions, we found 20 protein kinases, 11 nucleotide kinases, and four lipid kinases. An analysis of the homologous proteins of these kinases revealed that DR_A0332 (a non-specific serine/threonine protein kinase) is a protein kinase only present in bacteria. We also identified six site-specific kinases; the phosphorylation modification levels of these kinases were upregulated, while protein levels did not change significantly. This also included three glucokinases (DR_2635, PckA, GlgC), two lipid kinases (GGlpK, DR_1363), and a nucleotide kinase (DR_A0020).

3. Discussion

In this study, high-precision mass spectrometry and label-free quantification techniques were used to study the dynamic changes in the protein and phosphorylation level in D. radiodurans after heavy ion irradiation. In total, we quantified more than 2300 proteins and 1100 high-confidence phosphosites, which is the first phosphorylation dataset and the largest proteome dataset of D. radiodurans to date. Compared with previous studies, our criteria were more stringent and the results were better (Table S1D). In total, 74.9% of the upregulated proteins were unique in our dataset (Figure S4) [10,32,33]. Based on the integrated analysis of this high-quality dynamic (phospho)proteome dataset, several new insights into radio-resistance mechanisms were observed.
First of all, compared with protein abundance, the phosphorylation level has a much more significant upregulation after irradiation (Figure 1F). This phenomenon suggests that protein phosphorylation is more sensitive to irradiation than the expression level of proteins. A previous study also found this trend in the CVB3-virus-infected HeLa cells [54], but this is the first time that irradiated samples are being analyzed. Notably, the difference in the upregulation degree between the protein expression level and phosphorylation is much higher in D. radiodurans of this study. Thus, deciphering the dynamic changes in phosphorylation after irradiation is of great significance to fully understand the mechanism of extreme ability of radio tolerance in D. radiodurans.
Second, the dynamic change in the (phospho)proteome at different time points after irradiation could provide valuable information to understand the mechanisms of the immediate response and further repair process. However, prior to this study, there were no studies focusing on this question. Through the identification and in-depth analysis of samples at different time points after irradiation, we found that the radiation response of D. radiodurans is time-dependent. The functions of upregulated proteins at different time points were significantly different (Figure 2F). Early upregulated proteins were involved in stress-related functions, while late upregulated proteins were enriched in repair-related functions. For the significantly changed phosphoproteins, we distinguished four distinct expression patterns according to the dynamic changes at different time points (Figure 4A). Further detailed functional studies will provide more specific clues for the mechanism of irradiation response and repair.
Third, the significantly upregulated (phospho)proteins identified in this study could serve as potential radiation-responsive key molecules for further functional research. In the list of these proteins, a variety of known radiation-responsive proteins, including IrrE, PprA, DdrA, DdrD, PolA, PolB, UvrB, UvrC, RecF, RecO, RecR, are involved, proving the reliability of our research strategy. For the new potential key proteins identified, we further analyzed the possible interaction networks they may form and their potential functions. For example, a PPI network was constructed to screen core phosphoproteins. Three phosphoproteins, DR_0183 (large subunit of glutamate synthase), DR_0269 (a DR-specific highly disordered protein), and Eno (Enolase), were the top ones with the highest connectivity.
Fourth, we found that more than half of the radiation-responsive proteins showed inconsistent changes in protein expression levels and phosphorylation levels (Figure 5A,B). Among them, RpoB, a key protein in the RNA synthesis pathway, has no significant change in its protein expression level, but its expression is upregulated after irradiation at multiple phosphosites. RpoB protein can bind to rifampicin to inhibit transcriptional initiation. Previous studies evaluated the mutation of this gene by analyzing rifampicin-resistant clones to estimate the extent of DNA damage caused by radiation from the space station [55], but the link between its phosphorylation and radiation resistance was lacking in previous studies.
Finally, the results of this study may have potential application value. Some radiation-activated pathways, such as the mismatch repair pathway, RNA metabolic pathway, and antioxidant pathways (alpha-linolenic acid metabolism, phenylalanine metabolism, riboflavin metabolism), also exist in humans. As an essential fatty acid, alpha-linolenic acid has anti-inflammatory, anticancer and anti-metabolic syndrome properties [56]. Riboflavin, a type of vitamin B2, can play a protective role in nervous system diseases by regulating the mitochondria energy production–oxidative stress pathway [57]. Our study provides support for the potential anti-radiation applications of these edible organic compounds.
In summary, this study comprehensively analyzed the (phospho)proteome dynamic changes induced by heavy ion radiation from the aspects of differential expression, function, pathway, protein–protein interaction, subcellular localization, etc., which provide a fundamental framework for a system understanding of this radiation-tolerant organism. In particular, our discovery of core phosphoproteins, radiation-responsive kinases, and multiple radiation-responsive pathways will help us to better understand the underlying mechanisms and to build a foundation for further applied research.

4. Materials and Methods

4.1. Strain and Culture Conditions

D. radiodurans strains were purchased from the China General Microbiological Culture Collection Center (No. 1.3828, CGMCC, Beijing, China). The strains were cultured in TGY liquid medium (1% tryptone, 0.5% glucose, and 0.1% yeast extract) at 30 °C with shaking at 200 r/min.

4.2. Pre-Irradiation Preparation, Irradiation, and Post-Irradiation Treatment

D. radiodurans strains were cultivated to the early stationary phase (OD600 = 1.5). Samples were collected at 4000 r/min, at 4 ℃ for 3 min, and were washed three times with 1× PBS, and then finally re-suspended with 1× PBS (400 μL). The samples were subjected to heavy ion irradiation doses of 20 Gy, 80 Gy, and 160 Gy, respectively, at room temperature with a dose rate of 80 Gy/min (HeavyIon Research Facility in Lanzhou (HIRFL), Institute of Modern Physics, Chinese Academy of Sciences). In addition, a control group without exposure was set up. Immediately after irradiation, 1/2 of the samples were collected at 13,300 r/min for 5 min, and stored at −80 ℃. The remaining 1/2 samples were left at room temperature, and were collected after irradiation for 2 h.

4.3. Protein Lysis and Digestion for (Phospho)Proteomics Analysis

A total of 8M UA solution containing 1× protease and phosphatase inhibitor (Thermo, USA) was added to the bacterial precipitation for cell lysis and placed on ice for 30 min. Then, it was crushed on the ice with ultrasonic (SCIENTZ, China), and centrifuged at 4 ℃ at 14,000 r/min (Thermo, USA) for 5 min. Then, the precipitation and the surface oil were discarded, the supernatant was taken, and the protein concentration was determined using the BCA kit (CWBIO, China). In total, 1.5 mg protein was added to DTT (final concentration was 10 mM), bathed in water at 56 ℃ for 30 min, lowered to room temperature, added into 1 mol/L IAA solution, shielded from light for 30 min, and centrifuged at 14,000 r/min for 10 min. The protein supernatant was transferred to two 10 kD ultrafiltration tubes (Millipore, USA) for loading, and centrifuged for 14,000 r/min for 15 min. Then, 8M urea solution was added, followed by vortex oscillation, centrifuging at 14,000 r/min for 15 min 3 times, and the waste liquid was discarded. Next, 50 mM ammonium bicarbonate was added, vortexed for 5~10 s, centrifuged for 15 min 3 times at 14,000 r/min, and then the waste liquid was discarded. The new casing was replaced, 50 mM ammonium bicarbonate and 1 μg/μL trypsin (mass ratio 1/50, Promega, China) were added, shaken, and mixed, and finally digested at 37 ℃ for 12 h. This was centrifuged at 14,000 r/min for 10 min and the peptide solution was collected. Then, 1/10 of the peptide solution was taken in a new centrifuge tube for whole-protein sRP-HPLC identification, the remaining liquid was acidified with 0.1% formic acid, and then all peptide segments were frozen dry (Eppendorf, German) and stored at −80 ℃.
The lyophilized peptide was dissolved in a solution containing 1M lactic acid, 70% acetonitrile, and 5% TFA (Sigma, USA). The peptide was adsorbed on TiO2 (GL science, Japan) on a suspension apparatus. The peptide was washed with 30% acetonitrile, 0.5% TFA solution, 80% acetonitrile, and 0.5% TFA solution successively 3 times. The 2%, 5%, 8%, 10%, and 40% acetonitrile solutions were gradient-eluted once each, combined into 3 fractions, acidified, desalted, and identified via high-precision mass spectrometry (Thermo Scientific ™ Orbitrap Exploris 480) [58].

4.4. Chromatography and Mass Spectrometry

The 500 ng peptides were loaded onto 20 cm columns packed in-house with C18 particles (1.9 μm). The poly-peptide samples were dissolved in 10 μL of solvent A (0.1% FA in water), injected with 5μL, loaded onto the pre-column at a flow rate of 3 μL/min on the EASY-nano-LC chromatography system, and then separated on the column at a flow rate of 300 nL/min. The gradient was as follows: solvent B (0.1% FA in 80% ACN) increased linearly from 7% to 12% (0 min–10 min); solvent B rose linearly from 12% to 45% (10 min–70 min); and then solvent B was at 95% (55 min–90 min). The mass spectrum data were collected using FAIMS Pro, equipped with specific parameters that were set as follows: The spray voltage of ion source was set to 2.1 kV. FAIMS Pro adopts two voltages of −45 V and −65 V, and each cycle time was set to 1.3 s. Data-dependent acquisition with the MS method was used in one full scan (350–1500 m/z, R = 60,000 at 120 m/z, maximum injection time 50 ms), first with a target of 3 × 106 ions, followed by data-dependent MS/MS scans with higher-energy collision dissociation (AGC target 7.5e4 ions, maximum injection time at 22 ms, isolation window 1.6 m/z, normalized collision energy 27%, and R = 15,000 at 120 m/z). A dynamic exclusion of 20 s was enabled.

4.5. Database Searching, Quality Control and Data Normalization

The (phospho)proteome data were searched using Maxquant (version 2.0.3.0). The “match between runs” parameter was selected. The iBAQ (proteome) and intensity (phosphoproteome) data were normalized via median normalization. For the calculation of FDR values, we used the R function “anova” for processing normalized data. The Prcomp function was used for principal component analysis. Missing values were interpolated with global minimum values if necessary.

4.6. GO and KEGG Annotation

The GO annotation of D. radiodurans was downloaded from the UniProt database, in which there are 2103 annotated genes, and we traced the results, which finally corresponded to a total of 3186 GO terms. Pathway information was downloaded from the KEGG website and matched with the UniProt accession. Annotation information for all genes is listed in Table S2.

4.7. PPI Network Analysis

The STRING database was used to download the interaction network of all differential phosphoproteins. We selected the default parameters in the basic setting. Then, we drew the network diagram using Cytoscape [59]. Since some proteins have multiple differential phosphosites, we used different colors to indicate the cluster that the sites belonged to.

4.8. Protein Structure and Evolutionary Age Annotation

Protein domains were identified using the HMMER-scan program to search the Pfam databases. All Pfam-A domain entries were used to identify domains in each protein [60]. Finally, we obtained 2490 proteins with domain annotations. The domain information for all proteins is shown in Table S2.
Intrinsically disordered regions were identified using SPOTD (default parameters). SPOTD is used to calculate the structural disorder value of each amino acid residue [61,62]. An amino acid residue with a high structural disorder value, which exceeds the default threshold (0.426), was considered to be disordered. After we obtained the results, the structural disorder ratio (SDR, i.e., the percentage of disordered residues) was calculated and classified into 3 grades: 0–10%, highly structured; 10–30%, moderately disordered; 30–100%, highly disordered [63]. Then, we matched the result with phosphorylation site to determine whether the site was disordered or not.
To determine the evolutionary age of each gene in D. radiodurans, a phylostratigraphy approach [64,65,66] was used. We selected representative species with complete genome sequences in the UniProt database. In each genus, the species with the highest number of protein sequences was selected. In addition, all species of the Deinococcus genus were selected to determine whether the gene was genus-specific or species-specific. In total, 3394 species remained. We performed reciprocal alignment of all D. radiodurans protein sequences with the protein sequence of the other species using Diamond BLASTP [67]. The matches with both e-values less than 1E-5 were regarded as homologous genes. At last, we determined the gene age of each gene, and classified them into 4 grades: Cellular organism, Bacteria, Deinococcota, and Deinococcus radiodurans. Finally, we identified 498 Deinococcota-specific proteins and 164 Deinococcus radiodurans-specific proteins.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms241914817/s1.

Author Contributions

Conceptualization, D.Y. and S.L.; methodology, S.L.; software, F.W., Z.Y., H.C., Y.N. and C.C.; validation, D.Y.; formal analysis, S.L., F.W., Z.Y., H.C., Y.N. and C.C.; investigation, S.L.; resources, D.Y.; data curation, S.L.; writing—original draft preparation, S.L.; writing—review and editing, D.Y.; visualization, S.L.; supervision, D.Y.; project administration, D.Y.; funding acquisition, D.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (32270711), the Beijing Nova Program (20220484230), the fund project in the technology field of the basic strengthening plan (2019-JCJQ-JJ-165), and the State Key Laboratory of Proteomics (SKLP-K202004).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed in this study are included in this published article and the Supplementary Information Files. Supplementary Figures and Tables are available at IJMS online.

Acknowledgments

Thanks to the proteomics and bioinformatics platforms at the National Center for Protein Sciences (Beijing) for support in proteomics data generation and large-scale data processing.

Conflicts of Interest

The authors declare that they have no competing interest.

References

  1. Lim, S.; Jung, J.-H.; Blanchard, L.; de Groot, A. Conservation and diversity of radiation and oxidative stress resistance mechanisms in Deinococcus species. FEMS Microbiol. Rev. 2018, 43, 19–52. [Google Scholar] [CrossRef] [PubMed]
  2. Bauermeister, A.; Hahn, C.; Rettberg, P.; Reitz, G.; Moeller, R. Roles of DNA repair and membrane integrity in heat resistance of Deinococcus radiodurans. Arch. Microbiol. 2012, 194, 959–966. [Google Scholar] [CrossRef] [PubMed]
  3. Mogul, R.; Bol’Shakov, A.; Chan, S.L.; Stevens, R.M.; Khare, B.N.; Meyyappan, M.; Trent, J.D. Impact of Low-Temperature Plasmas on Deinococcus radiodurans and Biomolecules. Biotechnol. Prog. 2003, 19, 776–783. [Google Scholar] [CrossRef] [PubMed]
  4. Peltek, S.; Meshcheryakova, I.; Kiseleva, E.; Oshchepkov, D.; Rozanov, A.; Serdyukov, D.; Demidov, E.; Vasiliev, G.; Vinokurov, N.; Bryanskaya, A.; et al. E. coli aggregation and impaired cell division after terahertz irradiation. Sci. Rep. 2021, 11, 20464. [Google Scholar] [CrossRef] [PubMed]
  5. Shukla, S.K.; Manobala, T.; Rao, T.S. The role of S-layer protein (SlpA) in biofilm-formation of Deinococcus radiodurans. J. Appl. Microbiol. 2022, 133, 796–807. [Google Scholar] [CrossRef] [PubMed]
  6. Rothfuss, H.; Lara, J.C.; Schmid, A.K.; Lidstrom, M.E. Involvement of the S-layer proteins Hpi and SlpA in the maintenance of cell envelope integrity in Deinococcus radiodurans R1. Microbiology 2006, 152, 2779–2787. [Google Scholar] [CrossRef]
  7. Levin-Zaidman, S.; Englander, J.; Shimoni, E.; Sharma, A.K.; Minton, K.W.; Minsky, A. Ringlike Structure of the Deinococcus radiodurans Genome: A Key to Radioresistance? Science 2003, 299, 254–256. [Google Scholar] [CrossRef]
  8. Winter, M.; Dokic, I.; Schlegel, J.; Warnken, U.; Debus, J.; Abdollahi, A.; Schnölzer, M. Deciphering the Acute Cellular Phosphoproteome Response to Irradiation with X-rays, Protons and Carbon Ions. Mol. Cell. Proteom. 2017, 16, 855–872. [Google Scholar] [CrossRef]
  9. Mehnati, P.; Morimoto, S.; Yatagai, F.; Furusawa, Y.; Kobayashi, Y.; Wada, S.; Kanai, T.; Hanaoka, F.; Sasaki, H. Exploration of ‘Over Kill Effect’ of High-LET Ar- and Fe-ions by Evaluating the Fraction of Non-hit Cell and Interphase Death. J. Radiat. Res. 2005, 46, 343–350. [Google Scholar] [CrossRef]
  10. Gao, Y.; Li, N.; Zhou, Y.; Zhang, Z.; Zhang, Y.; Fan, P.; Zhou, H.; Zhang, T.; Chang, L.; Gao, H.; et al. iTRAQ-based proteomic analysis of Deinococcus radiodurans in response to 12C6+ heavy ion irradiation. BMC Microbiol. 2022, 22, 264. [Google Scholar] [CrossRef]
  11. Slade, D.; Lindner, A.B.; Paul, G.; Radman, M. Recombination and Replication in DNA Repair of Heavily Irradiated Deinococcus radiodurans. Cell 2009, 136, 1044–1055. [Google Scholar] [CrossRef] [PubMed]
  12. Pagès, V. Single-strand gap repair involves both RecF and RecBCD pathways. Curr. Genet. 2016, 62, 519–521. [Google Scholar] [CrossRef] [PubMed]
  13. Morimatsu, K.; Kowalczykowski, S.C. RecQ helicase and RecJ nuclease provide complementary functions to resect DNA for homologous recombination. Proc. Natl. Acad. Sci. USA 2014, 111, E5133–E5142. [Google Scholar] [CrossRef]
  14. Zahradka, K.; Slade, D.; Bailone, A.; Sommer, S.; Averbeck, D.; Petranovic, M.; Lindner, A.B.; Radman, M. Reassembly of shattered chromosomes in Deinococcus radiodurans. Nature 2006, 443, 569–573. [Google Scholar] [CrossRef] [PubMed]
  15. Seck, A.; De Bonis, S.; Stelter, M.; Ökvist, M.; Senarisoy, M.; Hayek, M.R.; Le Roy, A.; Martin, L.; Saint-Pierre, C.; Silveira, C.M.; et al. Structural and functional insights into the activation of the dual incision activity of UvrC, a key player in bacterial NER. Nucleic Acids Res. 2023, 51, 2931–2949. [Google Scholar] [CrossRef] [PubMed]
  16. Banneville, A.-S.; de la Tour, C.B.; De Bonis, S.; Hognon, C.; Colletier, J.-P.; Teulon, J.-M.; Le Roy, A.; Pellequer, J.-L.; Monari, A.; Dehez, F.; et al. Structural and functional characterization of DdrC, a novel DNA damage-induced nucleoid associated protein involved in DNA compaction. Nucleic Acids Res. 2022, 50, 7680–7696. [Google Scholar] [CrossRef] [PubMed]
  17. Agapov, A.A.; Kulbachinskiy, A.V. Mechanisms of stress resistance and gene regulation in the radioresistant bacterium Deinococcus radiodurans. Biochemistry 2015, 80, 1201–1216. [Google Scholar] [CrossRef]
  18. Selvam, K.; Duncan, J.R.; Tanaka, M.; Battista, J.R. DdrA, DdrD, and PprA: Components of UV and Mitomycin C Resistance in Deinococcus radiodurans R1. PLoS ONE 2013, 8, e69007. [Google Scholar] [CrossRef]
  19. Santos, S.P.; Mitchell, E.P.; Franquelim, H.G.; Castanho, M.A.R.B.; Abreu, I.A.; Romão, C.V. Dps from Deinococcus radiodurans: Oligomeric forms of Dps1 with distinct cellular functions and Dps2 involved in metal storage. FEBS J. 2015, 282, 4307–4327. [Google Scholar] [CrossRef]
  20. Reon, B.J.; Nguyen, K.H.; Bhattacharyya, G.; Grove, A. Functional comparison of Deinococcus radiodurans Dps proteins suggests distinct in vivo roles. Biochem. J. 2012, 447, 381–391. [Google Scholar] [CrossRef]
  21. Grove, A.; Wilkinson, S.P. Differential DNA Binding and Protection by Dimeric and Dodecameric forms of the Ferritin Homolog Dps from Deinococcus radiodurans. J. Mol. Biol. 2005, 347, 495–508. [Google Scholar] [CrossRef] [PubMed]
  22. Gaidamakova, E.K.; Sharma, A.; Matrosova, V.Y.; Grichenko, O.; Volpe, R.P.; Tkavc, R.; Conze, I.H.; Klimenkova, P.; Balygina, I.; Horne, W.H.; et al. Small-Molecule Mn Antioxidants in Caenorhabditis elegans and Deinococcus radiodurans Supplant MnSOD Enzymes during Aging and Irradiation. mBio 2022, 13, e0339421. [Google Scholar] [CrossRef]
  23. Zhou, C.; Dai, J.; Lu, H.; Chen, Z.; Guo, M.; He, Y.; Gao, K.; Ge, T.; Jin, J.; Wang, L.; et al. Succinylome Analysis Reveals the Involvement of Lysine Succinylation in the Extreme Resistance of Deinococcus radiodurans. Proteomics 2019, 19, e1900158. [Google Scholar] [CrossRef]
  24. Li, W.; Han, X.; Lan, P. Emerging roles of protein phosphorylation in plant iron homeostasis. Trends Plant Sci. 2022, 27, 908–921. [Google Scholar] [CrossRef] [PubMed]
  25. Sharma, K.; D’Souza, R.C.J.; Tyanova, S.; Schaab, C.; Wiśniewski, J.R.; Cox, J.; Mann, M. Ultradeep Human Phosphoproteome Reveals a Distinct Regulatory Nature of Tyr and Ser/Thr-Based Signaling. Cell Rep. 2014, 8, 1583–1594. [Google Scholar] [CrossRef] [PubMed]
  26. Sharma, D.K.; Siddiqui, M.Q.; Gadewal, N.; Choudhary, R.K.; Varma, A.K.; Misra, H.S.; Rajpurohit, Y.S. Phosphorylation of deinococcal RecA affects its structural and functional dynamics implicated for its roles in radioresistance of Deinococcus radiodurans. J. Biomol. Struct. Dyn. 2019, 38, 114–123. [Google Scholar] [CrossRef]
  27. Maurya, G.K.; Modi, K.; Banerjee, M.; Chaudhary, R.; Rajpurohit, Y.S.; Misra, H.S. Phosphorylation of FtsZ and FtsA by a DNA Damage-Responsive Ser/Thr Protein Kinase Affects Their Functional Interactions in Deinococcus radiodurans. mSphere 2018, 3, e00325–e18. [Google Scholar] [CrossRef]
  28. Rajpurohit, Y.S.; Bihani, S.C.; Waldor, M.K.; Misra, H.S. Phosphorylation of Deinococcus radiodurans RecA Regulates Its Activity and May Contribute to Radioresistance. J. Biol. Chem. 2016, 291, 16672–16685. [Google Scholar] [CrossRef]
  29. Chaudhary, R.; Mishra, S.; Maurya, G.K.; Rajpurohit, Y.S.; Misra, H.S. FtsZ phosphorylation brings about growth arrest upon DNA damage in Deinococcus radiodurans. FASEB BioAdvances 2022, 5, 27–42. [Google Scholar] [CrossRef]
  30. Potel, C.M.; Lin, M.-H.; Heck, A.J.R.; Lemeer, S. Widespread bacterial protein histidine phosphorylation revealed by mass spectrometry-based proteomics. Nat. Methods 2018, 15, 187–190. [Google Scholar] [CrossRef]
  31. Garcia-Garcia, T.; Poncet, S.; Derouiche, A.; Shi, L.; Mijakovic, I.; Noirot-Gros, M.-F. Role of Protein Phosphorylation in the Regulation of Cell Cycle and DNA-Related Processes in Bacteria. Front. Microbiol. 2016, 7, 184. [Google Scholar] [CrossRef]
  32. Xiong, Y.; Wei, L.; Xin, S.; Min, R.; Liu, F.; Li, N.; Zhang, Y. Comprehensive Temporal Protein Dynamics during Postirradiation Recovery in Deinococcus radiodurans. Oxidative Med. Cell. Longev. 2022, 2022, 1622829. [Google Scholar] [CrossRef] [PubMed]
  33. Ott, E.; Kawaguchi, Y.; Özgen, N.; Yamagishi, A.; Rabbow, E.; Rettberg, P.; Weckwerth, W.; Milojevic, T. Proteomic and Metabolomic Profiling of Deinococcus radiodurans Recovering After Exposure to Simulated Low Earth Orbit Vacuum Conditions. Front. Microbiol. 2019, 10, 909. [Google Scholar] [CrossRef] [PubMed]
  34. Stracke, C.; Meyer, B.H.; Hagemann, A.; Jo, E.; Lee, A.; Albers, S.-V.; Cha, J.; Bräsen, C.; Siebers, B. Salt Stress Response of Sulfolobus acidocaldarius Involves Complex Trehalose Metabolism Utilizing a Novel Trehalose-6-Phosphate Synthase (TPS)/Trehalose-6-Phosphate Phosphatase (TPP) Pathway. Appl. Environ. Microbiol. 2020, 86, 24. [Google Scholar] [CrossRef] [PubMed]
  35. Sakaguchi, M. Diverse and common features of trehalases and their contributions to microbial trehalose metabolism. Appl. Microbiol. Biotechnol. 2020, 104, 1837–1847. [Google Scholar] [CrossRef]
  36. Lu, R.-S.; Li, P.; Qiu, Y.-X. The Complete Chloroplast Genomes of Three Cardiocrinum (Liliaceae) Species: Comparative Genomic and Phylogenetic Analyses. Front. Plant Sci. 2017, 7, 2054. [Google Scholar] [CrossRef]
  37. Caserta, E.; Liu, L.-C.; Grundy, F.J.; Henkin, T.M. Codon-Anticodon Recognition in the Bacillus subtilis glyQS T Box Riboswitch. J. Biol. Chem. 2015, 290, 23336–23347. [Google Scholar] [CrossRef]
  38. Chakravarty, A.K.; Subbotin, R.; Chait, B.T.; Shuman, S. RNA ligase RtcB splices 3′-phosphate and 5′-OH ends via covalent RtcB-(histidinyl)-GMP and polynucleotide-(3′)pp(5′)G intermediates. Proc. Natl. Acad. Sci. USA 2012, 109, 6072–6077. [Google Scholar] [CrossRef]
  39. Zhang, Y.; Morar, M.; Ealick, S.E. Structural biology of the purine biosynthetic pathway. Cell. Mol. Life Sci. 2008, 65, 3699–3724. [Google Scholar] [CrossRef]
  40. Fujimura, S.; Mio, K.; Ohkubo, T.; Arai, T.; Kuramochi, M.; Sekiguchi, H.; Sasaki, Y.C. Diffracted X-ray Tracking Method for Measuring Intramolecular Dynamics of Membrane Proteins. Int. J. Mol. Sci. 2022, 23, 2343. [Google Scholar] [CrossRef]
  41. Cai, L.; Zhao, D.; Hou, J.; Wu, J.; Cai, S.; Dassarma, P.; Xiang, H. Cellular and organellar membrane-associated proteins in haloarchaea: Perspectives on the physiological significance and biotechnological applications. Sci. China Life Sci. 2012, 55, 404–414. [Google Scholar] [CrossRef] [PubMed]
  42. Yang, Y.; Jin, H.; Chen, Y.; Lin, W.; Wang, C.; Chen, Z.; Han, N.; Bian, H.; Zhu, M.; Wang, J. A chloroplast envelope membrane protein containing a putative LrgB domain related to the control of bacterial death and lysis is required for chloroplast development in Arabidopsis thaliana. New Phytol. 2011, 193, 81–95. [Google Scholar] [CrossRef] [PubMed]
  43. Kanehisa, M.; Furumichi, M.; Sato, Y.; Ishiguro-Watanabe, M.; Tanabe, M. KEGG: Integrating viruses and cellular organisms. Nucleic Acids Res. 2020, 49, D545–D551. [Google Scholar] [CrossRef] [PubMed]
  44. Xu, Q.; Cao, Y.; Zhong, X.; Qin, X.; Feng, J.; Peng, H.; Su, Y.; Ma, Z.; Zhou, S. Riboflavin protects against heart failure via SCAD-dependent DJ-1–Keap1–Nrf2 signalling pathway. Br. J. Pharmacol. 2023, 1–23. [Google Scholar] [CrossRef]
  45. Alam, S.-I.; Kim, M.-W.; Shah, F.A.; Saeed, K.; Ullah, R.; Kim, M.-O. Alpha-Linolenic Acid Impedes Cadmium-Induced Oxidative Stress, Neuroinflammation, and Neurodegeneration in Mouse Brain. Cells 2021, 10, 2274. [Google Scholar] [CrossRef]
  46. Feng, L.; Li, W.; Liu, Y.; Jiang, W.-D.; Kuang, S.-Y.; Wu, P.; Jiang, J.; Tang, L.; Tang, W.-N.; Zhang, Y.-A.; et al. Protective role of phenylalanine on the ROS-induced oxidative damage, apoptosis and tight junction damage via Nrf2, TOR and NF-κB signalling molecules in the gill of fish. Fish Shellfish. Immunol. 2016, 60, 185–196. [Google Scholar] [CrossRef]
  47. Bentchikou, E.; Servant, P.; Coste, G.; Sommer, S. A Major Role of the RecFOR Pathway in DNA Double-Strand-Break Repair through ESDSA in Deinococcus radiodurans. PLOS Genet. 2010, 6, e1000774. [Google Scholar] [CrossRef]
  48. Szklarczyk, D.; Kirsch, R.; Koutrouli, M.; Nastou, K.; Mehryary, F.; Hachilif, R.; Gable, A.L.; Fang, T.; Doncheva, N.T.; Pyysalo, S.; et al. The STRING database in 2023: Protein–protein association networks and functional enrichment analyses for any sequenced genome of interest. Nucleic Acids Res. 2022, 51, D638–D646. [Google Scholar] [CrossRef]
  49. Jiao, C.; Gu, Z. iTRAQ-based proteomic analysis reveals changes in response to UV-B treatment in soybean sprouts. Food Chem. 2018, 275, 467–473. [Google Scholar] [CrossRef]
  50. Prasad, S.M.; Kumar, S.; Parihar, P.; Singh, A.; Singh, R. Evaluating the combined effects of pretilachlor and UV-B on two Azolla species. Pestic. Biochem. Physiol. 2016, 128, 45–56. [Google Scholar] [CrossRef]
  51. Huang, G.; Wang, L.; Zhou, Q. Lanthanum (III) Regulates the Nitrogen Assimilation in Soybean Seedlings under Ultraviolet-B Radiation. Biol. Trace Element Res. 2012, 151, 105–112. [Google Scholar] [CrossRef] [PubMed]
  52. Boothby, T.C.; Tapia, H.; Brozena, A.H.; Piszkiewicz, S.; Smith, A.E.; Giovannini, I.; Rebecchi, L.; Pielak, G.J.; Koshland, D.; Goldstein, B. Tardigrades Use Intrinsically Disordered Proteins to Survive Desiccation. Mol. Cell 2017, 65, 975–984. [Google Scholar] [CrossRef] [PubMed]
  53. Ramasamy, P.; Vandermarliere, E.; Vranken, W.F.; Martens, L. Panoramic Perspective on Human Phosphosites. J. Proteome Res. 2022, 21, 1894–1915. [Google Scholar] [CrossRef]
  54. Giansanti, P.; Strating, J.R.P.M.; Defourny, K.A.Y.; Cesonyte, I.; Bottino, A.M.S.; Post, H.; Viktorova, E.G.; Ho, V.Q.T.; Langereis, M.A.; Belov, G.A.; et al. Dynamic remodelling of the human host cell proteome and phosphoproteome upon enterovirus infection. Nat. Commun. 2020, 11, 4332. [Google Scholar] [CrossRef]
  55. Fujiwara, D.; Kawaguchi, Y.; Kinoshita, I.; Yatabe, J.; Narumi, I.; Hashimoto, H.; Yokobori, S.-I.; Yamagishi, A. Mutation Analysis of the rpoB Gene in the Radiation-Resistant Bacterium Deinococcus radiodurans R1 Exposed to Space during the Tanpopo Experiment at the International Space Station. Astrobiology 2021, 21, 1494–1504. [Google Scholar] [CrossRef] [PubMed]
  56. Yuan, Q.; Xie, F.; Huang, W.; Hu, M.; Yan, Q.; Chen, Z.; Zheng, Y.; Liu, L. The review of alpha-linolenic acid: Sources, metabolism, and pharmacology. Phytotherapy Res. 2021, 36, 164–188. [Google Scholar] [CrossRef] [PubMed]
  57. Fila, M.; Chojnacki, C.; Chojnacki, J.; Blasiak, J. Nutrients to Improve Mitochondrial Function to Reduce Brain Energy Deficit and Oxidative Stress in Migraine. Nutrients 2021, 13, 4433. [Google Scholar] [CrossRef]
  58. Stejskal, K.; de Beeck, J.O.; Dürnberger, G.; Jacobs, P.; Mechtler, K. Ultrasensitive NanoLC-MS of Subnanogram Protein Samples Using Second Generation Micropillar Array LC Technology with Orbitrap Exploris 480 and FAIMS PRO. Anal. Chem. 2021, 93, 8704–8710. [Google Scholar] [CrossRef]
  59. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of Biomolecular Interaction Networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  60. Wang, H.; Zhong, H.; Gao, C.; Zang, J.; Yang, D. The Distinct Properties of the Consecutive Disordered Regions Inside or Outside Protein Domains and Their Functional Significance. Int. J. Mol. Sci. 2021, 22, 10677. [Google Scholar] [CrossRef]
  61. Hanson, J.; Paliwal, K.K.; Zhou, Y. Accurate Single-Sequence Prediction of Protein Intrinsic Disorder by an Ensemble of Deep Recurrent and Convolutional Architectures. J. Chem. Inf. Model. 2018, 58, 2369–2376. [Google Scholar] [CrossRef] [PubMed]
  62. Gao, C.; Ma, C.; Wang, H.; Zhong, H.; Zang, J.; Zhong, R.; He, F.; Yang, D. Intrinsic disorder in protein domains contributes to both organism complexity and clade-specific functions. Sci. Rep. 2021, 11, 2985. [Google Scholar] [CrossRef] [PubMed]
  63. Gsponer, J.; Futschik, M.E.; Teichmann, S.A.; Babu, M.M. Tight Regulation of Unstructured Proteins: From Transcript Synthesis to Protein Degradation. Science 2008, 322, 1365–1368. [Google Scholar] [CrossRef] [PubMed]
  64. Trigos, A.S.; Pearson, R.B.; Papenfuss, A.T.; Goode, D.L. Altered interactions between unicellular and multicellular genes drive hallmarks of transformation in a diverse range of solid tumors. Proc. Natl. Acad. Sci. USA 2017, 114, 6406–6411. [Google Scholar] [CrossRef] [PubMed]
  65. Domazet-Lošo, T.; Tautz, D. A phylogenetically based transcriptome age index mirrors ontogenetic divergence patterns. Nature 2010, 468, 815–818. [Google Scholar] [CrossRef]
  66. Domazet-Lošo, T.; Brajković, J.; Tautz, D. A phylostratigraphy approach to uncover the genomic history of major adaptations in metazoan lineages. Trends Genet. 2007, 23, 533–539. [Google Scholar] [CrossRef]
  67. Buchfink, B.; Reuter, K.; Drost, H.-G. Sensitive protein alignments at tree-of-life scale using DIAMOND. Nat. Methods 2021, 18, 366–368. [Google Scholar] [CrossRef]
Ijms 24 14817 g001
Figure 1. Experimental design, quality control and overview of (phospho)proteome data. (A) Experimental design for (phospho)proteome analysis of heavy ion irradiated D. radiodurans. Strains were treated with 3 doses of radiation from 20 to 160 Gy. Samples were collected at two time points (0 h, 2 h) after irradiation. The color scheme for each dose remained essentially the same throughout the paper. (B) A comparison with the count of proteins identified by D. radiodurans proteomics published in the last 3 years [10,32,33]. (C) The number of proteins and phosphosites identified as a whole and the average number identified per sample. (D) PCA analysis of the global proteome samples. Different colors/shapes represent different doses/time points. (E) Fraction of proteins and phosphosites dynamically regulated after multiple dose irradiation, determined by using a two-way ANOVA test (FDR < 0.05). (F) Distribution of the scale of changes of significantly regulated proteins and phosphosites, showing that the differences in phosphorylation are generally more extensive than those at protein levels.
Figure 1. Experimental design, quality control and overview of (phospho)proteome data. (A) Experimental design for (phospho)proteome analysis of heavy ion irradiated D. radiodurans. Strains were treated with 3 doses of radiation from 20 to 160 Gy. Samples were collected at two time points (0 h, 2 h) after irradiation. The color scheme for each dose remained essentially the same throughout the paper. (B) A comparison with the count of proteins identified by D. radiodurans proteomics published in the last 3 years [10,32,33]. (C) The number of proteins and phosphosites identified as a whole and the average number identified per sample. (D) PCA analysis of the global proteome samples. Different colors/shapes represent different doses/time points. (E) Fraction of proteins and phosphosites dynamically regulated after multiple dose irradiation, determined by using a two-way ANOVA test (FDR < 0.05). (F) Distribution of the scale of changes of significantly regulated proteins and phosphosites, showing that the differences in phosphorylation are generally more extensive than those at protein levels.
Ijms 24 14817 g001
Ijms 24 14817 g002
Figure 2. Integration analysis of irradiation sensitive proteins. (A) The number of up- or downregulated proteins at different irradiation doses and time points. (B) Heatmap of z-scored protein abundance (iBAQ value) of the DEPs. (C,D) Venn diagram of the upregulated proteins at 3 radiation doses, and 0h (C) or 2h (D) after irradiation. (E,F) Function enrichment analysis of 0h upregulated proteins with different doses (E), and of upregulated proteins at different time points (F). The red scale represents over-representation, while the blue scale represents under-representation.
Figure 2. Integration analysis of irradiation sensitive proteins. (A) The number of up- or downregulated proteins at different irradiation doses and time points. (B) Heatmap of z-scored protein abundance (iBAQ value) of the DEPs. (C,D) Venn diagram of the upregulated proteins at 3 radiation doses, and 0h (C) or 2h (D) after irradiation. (E,F) Function enrichment analysis of 0h upregulated proteins with different doses (E), and of upregulated proteins at different time points (F). The red scale represents over-representation, while the blue scale represents under-representation.
Ijms 24 14817 g002
Ijms 24 14817 g003
Figure 3. Irradiation response-related pathways. (A) The changes in all pathways of D. radiodurans after irradiation. The ring diagram from outside to inside is as follows: pathway class (a), pathway subclass (b), proportion of proteins identified (c), proportion of upregulated proteins (d), proportion of phosphoproteins (e), and proportion of proteins with upregulated phosphorylation level (f) (Table S3A). (B) Homologous recombination (RecFOR) and mismatch repair pathways. Red color indicates upregulation; (p) means phosphorylation in the protein. The differential expressions of proteins in the pathways are labeled in the table below.
Figure 3. Irradiation response-related pathways. (A) The changes in all pathways of D. radiodurans after irradiation. The ring diagram from outside to inside is as follows: pathway class (a), pathway subclass (b), proportion of proteins identified (c), proportion of upregulated proteins (d), proportion of phosphoproteins (e), and proportion of proteins with upregulated phosphorylation level (f) (Table S3A). (B) Homologous recombination (RecFOR) and mismatch repair pathways. Red color indicates upregulation; (p) means phosphorylation in the protein. The differential expressions of proteins in the pathways are labeled in the table below.
Ijms 24 14817 g003
Ijms 24 14817 g004
Figure 4. Integration analysis of phosphoproteins and phosphosites. (A) Fuzzy c-means clustering of 4 dynamic patterns of differential phosphosites. (B) PPI network of interactions between proteins with differential phosphosites. Displayed network includes 25 core proteins and 85 suggested proteins that have at least one interaction with other protein(s). (C) Overview of major functional classification of proteins with differential phosphosites.
Figure 4. Integration analysis of phosphoproteins and phosphosites. (A) Fuzzy c-means clustering of 4 dynamic patterns of differential phosphosites. (B) PPI network of interactions between proteins with differential phosphosites. Displayed network includes 25 core proteins and 85 suggested proteins that have at least one interaction with other protein(s). (C) Overview of major functional classification of proteins with differential phosphosites.
Ijms 24 14817 g004
Ijms 24 14817 g005
Figure 5. Comparative analysis of protein expression and phosphorylation modification. (A) Classification of protein expression and phosphorylation levels after irradiation. (B) The proportion of different classes in different irradiation states. (C) Functional enrichment analysis of various classified proteins. (D) Dynamic changes in modification levels of site-specific protein RpoB.
Figure 5. Comparative analysis of protein expression and phosphorylation modification. (A) Classification of protein expression and phosphorylation levels after irradiation. (B) The proportion of different classes in different irradiation states. (C) Functional enrichment analysis of various classified proteins. (D) Dynamic changes in modification levels of site-specific protein RpoB.
Ijms 24 14817 g005
Ijms 24 14817 g006
Figure 6. The classification of proteins upregulated in at least one state after irradiation which have a kinase domain. Blue to red color gradient denotes downregulation to upregulation compared with unirradiated samples at this time point. Fold changes were log2 transformed. The text and color on the right, respectively, indicate the name and category of the kinase domain corresponding to the protein. The protein labeled with ※ is a bacteria-specific protein.
Figure 6. The classification of proteins upregulated in at least one state after irradiation which have a kinase domain. Blue to red color gradient denotes downregulation to upregulation compared with unirradiated samples at this time point. Fold changes were log2 transformed. The text and color on the right, respectively, indicate the name and category of the kinase domain corresponding to the protein. The protein labeled with ※ is a bacteria-specific protein.
Ijms 24 14817 g006
Table 1. RNA metabolism-related upregulated proteins or differential phosphosites.
Table 1. RNA metabolism-related upregulated proteins or differential phosphosites.
Uniprot IDProtein NameThe Upregulated Experiment Groups Dose (Gy)-Time (h)Phosphosites
Amino Acid-Cluster
Q9RV39A-adding tRNA nucleotidyltransferase160-0NULL
Q9RR60Enolase20-0, 80-0, 160-0,20-2, 80-2, 160-2T252-4, Y256-1
Q9RSR1Polyribonucleotide nucleotidyltransferaseNULLT750-2
Q9RY23Chaperone protein DnaK20-0, 80-0, 160-0,20-2, 80-2, 160-2S113-2, T133-3
Q9RWQ9Chaperonin GroEL20-0, 80-0, 160-0,20-2, 80-2, 160-2T80-2, T268-2
Q9RSJ6RNA polymerase subunit alphaNULLY28-2
Q9RVV9RNA polymerase subunit betaNULLS53-4, Y98-3, T364-2, S606-1, S670-3, Y1107-4
Q9RVW0RNA polymerase subunit beta’80-0, 160-0, 20-2, 160-2S339-4
Q9RRJ6RNA polymerase subunit omega80-0, 160-0, 160-2NULL
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liu, S.; Wang, F.; Chen, H.; Yang, Z.; Ning, Y.; Chang, C.; Yang, D. New Insights into Radio-Resistance Mechanism Revealed by (Phospho)Proteome Analysis of Deinococcus Radiodurans after Heavy Ion Irradiation. Int. J. Mol. Sci. 2023, 24, 14817. https://doi.org/10.3390/ijms241914817

AMA Style

Liu S, Wang F, Chen H, Yang Z, Ning Y, Chang C, Yang D. New Insights into Radio-Resistance Mechanism Revealed by (Phospho)Proteome Analysis of Deinococcus Radiodurans after Heavy Ion Irradiation. International Journal of Molecular Sciences. 2023; 24(19):14817. https://doi.org/10.3390/ijms241914817

Chicago/Turabian Style

Liu, Shihao, Fei Wang, Heye Chen, Zhixiang Yang, Yifan Ning, Cheng Chang, and Dong Yang. 2023. "New Insights into Radio-Resistance Mechanism Revealed by (Phospho)Proteome Analysis of Deinococcus Radiodurans after Heavy Ion Irradiation" International Journal of Molecular Sciences 24, no. 19: 14817. https://doi.org/10.3390/ijms241914817

APA Style

Liu, S., Wang, F., Chen, H., Yang, Z., Ning, Y., Chang, C., & Yang, D. (2023). New Insights into Radio-Resistance Mechanism Revealed by (Phospho)Proteome Analysis of Deinococcus Radiodurans after Heavy Ion Irradiation. International Journal of Molecular Sciences, 24(19), 14817. https://doi.org/10.3390/ijms241914817

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop