Hydrogen peroxide treatment induces the transposition of an insertion sequence in Deinococcus radiopugnans DY59

Deinococcus radiopugnans DY59 (formerly Deinococcus swuensis DY59) is a radiation-resistant bacterium isolated from soil. From the 3.5 Mb genomic DNA sequence of strain DY59 (December 2014), 31 insertion sequence (IS) elements of six IS families including IS1, IS4, IS5, IS66, IS630, and IS701 and five unclassified IS elements were detected. Upon induction of oxidative stress with 80 and 100 mM H2O2, the unique ISs of the IS4 family member were actively translocated into a carotenoid biosynthesis gene phytoene desaturase (QR90_10400), resulting in non-pigment phenotypic selection. Therefore, these active transpositions of a specific IS family member were induced by oxidative stress at 80 and 100 mM H2O2. Furthermore, D. radiopugnans DY59 exhibited extremely higher MIC values against H2O2 treatment. To explain this phenomenon, qRT-PCR was conducted to assess the expression levels of catalase and three LysR family regulators. Our findings indicated that the ISDrpg2 and ISDrpg3 elements of the IS4 family were actively transposed into the phytoene desaturase gene by H2O2 treatment via replicative transposition. However, high H2O2 resistance did not originate from H2O2-induced expression of catalase and LysR family regulators.

Prokaryotic genomes exhibit the smallest transposable element insertion sequences (ISs). They have important roles in genomic evolution including the enhancement of gene inactivation and genome plasticity (Touchon and Rocha, 2007;Siguier et al., 2014Siguier et al., , 2015. IS abundance is positively correlated with the frequency of horizontal gene transfer (HGT), genome size, pathogenicity, non-obligatory ecological associations, and human association (Touchon and Rocha, 2007;Vandecraen et al., 2017). Moreover, recent mobilome studies have characterized IS distribution at the genus level (Blesa et al., 2019;Fayad et al., 2019). IS elements are transferred between genomes through all classical mechanisms of HGT (Frost et al., 2005). IS elements generally consist of one or two transposase ORFs and a terminal inverted repeat (TIR) sequence. Transposases are the most prevalent genes in nature. When IS element integrated into genomic DNA, the direct repeat (DR) sequences were produced generally. Unique and random DR sequences in border of IS elements indicated a specific and random sequence recognition and integration. Their mobile nature not only promotes the dissemination of transposable elements within and between genomes but also leads to mutations and rearrangements that can accelerate biological diversification and consequent evolution (Aziz et al., 2010).
In the case of Deinococcus geothermalis (Dgeo), 73 full-length IS elements belonging to nine IS families were distributed across the different molecules of the genome including two mega plasmids. Several IS elements were actively transposed into other genome loci through the copy-and-paste mechanism of ISDge5, ISDge6, and ISDge7 members under hydrogen peroxide (H 2 O 2 )-induced oxidative stress. Furthermore, we reported the occurrence of a colorless phenotype through the loss of carotenoid biosynthesis via IS transposition, which was supported by several lines of evidence (Tian and Hua, 2010;Lee et al., 2019Lee et al., , 2020Ye et al., 2021;Shin et al., 2022). Research on the insertion sequence of bacteria can be conducted from a molecular evolution perspective, including cutting-edge mutagenesis using transposable elements to examine the plasticity of bacterial genomes induced by long-term culture or oxidative stress. It can also serve as a basic study of genome stability for the potential application of bioremediation in the genus Deinococcus species (Gerber et al., 2015;Wright et al., 2017;Consuegra et al., 2021).
In this work, we describe that the genome of D. radiopugnans DY59 has a total of 36 IS elements with six IS families including five unclassified IS elements. Interestingly, strain DY59_IS elements have multiple copies of the IS66, IS630, IS701, IS4, and IS5 families. Nevertheless, when strain DY59 was exposed to 5 kGy of gamma irradiation, non-pigmented mutants were not observed . Here, we evaluated the occurrence of active transposition induced by short-term oxidative stress conditions and its effects on the genomic plasticity in the radiation-resistant bacterium D. radiopugnans DY59. When oxidative stress was induced via hydrogen peroxide treatment with 80 mM, the ISDrpg2 and ISDrpg3 IS elements of the IS4 family were transposed into a carotenoid synthesis gene encoding a phytoene dehydrogenase (QR90_10400), resulting in a non-pigment phenotype. Furthermore, strain DY59 exhibited a high minimum inhibitory concentration (MIC) of 300 mM H 2 O 2 . Thus, the expression levels of catalase and LysR family members (i.e., its putative controlling regulators) were determined by qRT-PCR.

Materials and methods
Bacterial species and culture conditions Deinococcus radiopugnans DY59 was kindly provided by Dr. M. Kim at Seoul Women's University. This strain was deposited into the Korean Collection for Type Cultures (KCTC33033 T ). Deinococcus geothermalis DSM11300 T and D. radiodurans R1 were used as a positive control for the active transposition of ISs. Deinococcus geothermalis and D. radiodurans R1 were cultured on TGY medium containing 1% tryptone, 0.5% yeast extract, and 0.1% glucose with 150 rpm for broth culture at 48 and 30°C, respectively. Strain DY59 is known to grow well in R2A or LB medium. Therefore, the cells were cultured for 3-4 days at 30°C in R2A medium, whereas the cells in the TGY medium required only 2-3 days to grow. All downstream Deinococcus cultures were thus conducted using TGY culture medium.

Selection and determination of active transposition
Deinococcus radiopugnans DY59 cells were typically grown to a maximum optical density (OD 600 ) for 2-3 days. When the culture reached an absorbance of 2.0 or 4.0 at OD 600 as measured by a UV-VIS spectrophotometer, the cells were harvested by centrifugation at 3,000 rpm and resuspended in a 0.9% NaCl solution to an OD 600 of 2.0 for exposure to H 2 O 2 .
Next, the cells were treated with final H 2 O 2 concentration ranges of 80 and 100 mM/200 and 300 or 0 mM as a control and continuously cultured for 1 h at 30°C with 150 rpm shaking. The samples were then directly diluted to 10 −4 -10 −5 with 0.9% NaCl solution and 100 μL of the cell suspensions were spread on TGY agar medium using sterilized glass beads. After discarding the glass beads, the TGY plates were incubated at 30°C for 2-3 days. Non-pigmented colonies were detected and cell culture dilutions were streaked on TGY agar plates for pure cultivation. Non-pigment forming frequency was calculated the ratio of number of non-pigment colony and CFU. We also evaluated whether hypochlorite and gamma irradiation induced the occurrence of IS transposition. Active transposition was detected via PCR amplification of four selected target genes related to carotenoid biosynthesis (Supplementary Table S1).

Detection of IS elements from genome sequence and non-pigment mutants
The D. radiopugnans DY59 genome sequence was obtained from the NCBI database (GenBank accession number NZ_CP010028.1). IS elements were detected from the annotated transposase genes using Frontiers in Microbiology 03 frontiersin.org protein profiling data obtained from the GenBank database. We first extracted sequences from upstream and downstream of a transposase (1 kb region; 3 kb total length). Then, terminal inverted repeat (TIR) sequences and direct repeat (DR) sequences were determined using DNA sequence alignment software such as BLAST, ClustalW, and BioEdit. All detected IS elements were sorted into IS family members. Unfortunately, the IS detection platform ISfinder (https://isfinder. biotoul.fr; Siguier et al., 2006) still does not support the detection and classification of D. radiopugnans DY59 IS elements with calling ID as of September 2022. Therefore, we classified the IS ourselves based on IS structure analysis and nomenclature rules (Supplementary Table S2).

Assessment of various stressors and MIC of antibiotics
Afterward, 10 ml of cultured cells till to an OD 600 of 2.0 and 4.0 were harvested, resuspended with 0.9% NaCl, and separately exposed for 1 h to various H 2 O 2 concentrations or 100 μg/mL sodium hypochlorite treatment or for in total 5 kGy gamma-irradiation (160 Gy/min). Then, to measure the cell viability using colony forming units, samples were diluted to a 10 −5 ratio and spread on TGY agar plates. After 5 kGy gamma-irradiation, DY59 strain was exhibited 0.5% viability comparing to non-radiation. The viability of hypochlorite treatment was not analytical determined. To measure the minimum inhibitory concentration (MIC) of antibiotics, the strain DY59 cells of 100 μL on OD 600 = 4.0 were spread on TGY agar plates, and the disk diffusion assay was conducted with different concentrations of 10 μL antibiotic solutions, such as streptomycin, kanamycin, and ampicillin. The MIC values were determined based on the formation of a clear zone.
qRT-PCR analysis of catalase and three LysR family members Strain DY59 cells were harvested at OD 600 values of 2.0 and 4.0, resuspended to OD 600 of 2.0 in 5 mL 0.9% NaCl, and then treated with 50-150 mM H 2 O 2 for 1 h at 150 rpm. After the H 2 O 2 challenge, the samples were centrifuged at 10,000 rpm for 5 min and washed once more with 0.9% NaCl. The supernatants were discarded and the pellets were stored at −20°C. The cell wall was broken using phenol and DNA digestion was performed using DNase I. RNA was extracted using an RNA prep kit for RNA isolation (RNeasy mini purification kit; Qiagen, Germany). After measuring the extracted RNA concentration, the concentration was normalized to 1,000 ng in an 8 μL volume for all samples. cDNA synthesis was performed with a dNTP mixture and 6-mer random primers using the following protocol: 60°C for 5 min, 4°C for 3 min, 30°C for 10 min, 42°C for 60 min, and finally 95°C for 5 min (PrimeScript™ first strand cDNA Synthesis Kit; TaKaRa, Japan). In the step at 4°C for 3 min, 4 μL 5× buffer, 4.5 μL RNase free water, 1 μL RTase, and 0.5 μL RNase inhibitor were added. qRT-PCR analysis was performed using TB Green ® Premix Ex Taq™ (TaKaRa, Japan) on a Bio-RAD RT-PCR model CFX96™ Optics Module (Bio-RAD, United States). The expression of catalase (QR90_06310) and three LysR family regulators (QR90_13110, QR90_14595, and QR90_15105) was normalized to that of GAPDH, a gene that is constitutively expressed throughout all growth phases. The related expression levels of both catalase and LysR genes were calculated as described in a previous study . Differences between samples were determined via the Student's t-test using the Prism™ software (ver. 8.0). Differences were deemed statistically significant at p < 0.05 (*) and p < 0.01 (**). The primer sequences for qRT-PCR were shown in Supplementary Table S3.

Physiological properties and genetic distribution of carotenoid biosynthesis
Generally, D. radiopugnans DY59 formed reddish-colored colonies after being cultured for 2-3 days on TGY medium at 30°C. The strain also exhibited a stronger viscosity compared to D. geothermalis. Here, we observed different MIC values, which were indicative of different levels of antibiotic resistance. The MIC of streptomycin was less than 25 μg/mL, whereas those of ampicillin, kanamycin, and chloramphenicol were 30, 300, and 60 μg/mlL, respectively ( Figure 1A; Supplementary Figure S1). Therefore, we conclude that strain DY59 displays a resistant phenotype to kanamycin and chloramphenicol, and a tolerant phenotype to streptomycin and ampicillin when compared to the typical working concentrations of antibiotics.
Surprisingly, the results of our H 2 O 2 viability assays indicated that D. radiopugnans DY59 could form colonies after being exposed to H 2 O 2 concentrations of up to 600 mM for 1 h, which is two times higher than the resistance of two well-studied control strains of D. radiodurans and D. geothermalis (300 mM H 2 O 2 ; Figure 1B).
Genes related to carotenoid biosynthesis in genus Deinococcus, which were linked to the phenotypic reddish color of the colonies, were detected via KEGG pathway analysis and genome-wide genomics studies (Tian and Hua, 2010;Lim et al., 2019). We selected four genes that are involved in the carotenoid biosynthesis pathway, which were marked in a simplified pathway schematic ( Figure 2). The selected four carotenoid biosynthesis genes of strain DY59 exhibited amino acid sequence similarity with more than 60% to D. geothermalis genes: QR90_03795, a phytoene synthase, is 69.4% similar to Dgeo_0523; QR90_10400, a phytoene dehydrogenase (desaturase), is 78.9% similar to Dgeo_0524; QR90_14380, a carotenoid hydratase, is 58.6% similar to Dgeo_2309; and QR90_14400, a FAD-dependent oxidoreductase, is 74.5% similar to Dgeo_2306 (Supplementary Table S1). Interestingly, some carotenoid biosynthesis-related genes of D. geothermalis were adjacently clustered in the genome, for examples crtB-crtI and cruC-cruD-cruF-crtO. In contrast, each carotenoid biosynthesis gene in D. radiopugnans was separated.

Distribution of IS elements in the DY59 genome
From the NCBI genome sequence data (submitted at December 2014) and protein information of DY59 (Kim et al., 2015), we first selected transposases and their genomic loci and defined full-length IS elements. The genome of strain DY59 contains six defined IS family members with DDE type transposase (Tpase) including IS1, IS4, IS5, IS66, IS630, and IS701, as well as five unclassified Tpases within a total  Figure S2). In general, Tpase is a major component of IS elements and simplified schematics of the IS structures except for the IS1 family member, which has only a partial fragment, are illustrated in Supplementary Figure S3. The IS1 family member is an IS element (QR90_RS04110) with a 119 aa-length Tpase. In general, the IS1 family IS element is composed of two transposase-related open reading frames (ORFs) with the same transcriptional direction. Therefore, TIR and DR sequences for this IS region of the IS1 family could not be easily determined because it was only a partial fragment.
The IS4 family members had five copies in the genome. The IS elements including transposases were separated into two subtypes. QR90_RS09690, QR90_04660, and QR90_01215 have identical IS elements with a 320 aa-long transposase, a conserved TIR sequence (CTCTGTACCGGACAACT), and DR sequences of variable sequences with 9 nt ("GCCGTGATC, " "GTCACGCCC, " and "GAAAGCAAT, " respectively). QR90_RS10860 and QR90_04150 possess another identical IS element with an identical transposase with 327 aa-length. Both IS elements have a conserved TIR sequence "CTCGGTAGCTGACAACTTCA" and DR sequences "AGATTGAAC" and "CAGGGTCAG. " Thus, these IS4 family members were divided into the two subtypes ISDrpg2 and ISDrpg3, with their Tpase sequences having a 67.19% amino acid identity. The IS5 family members have four copies in the genome with two subtypes elements: QR90_RS06350 with a 277 aa-length Tpase and QR90_RS04355, QR90_RS07350, and QR90_06585 with 265 aa-length Tpase sequences. The three 265 amino acid-long transposase genes are highly similar, with similarity rates of 97.7-100%. However, the Tpase of QR90_RS06350 with a 277 amino acid length exhibited a quite low similarity of 12-12.8% to the 265 aa-length Tpase in the IS5 family. Therefore, the IS5 family members were separated into two subtypes: ISDrpg4 and ISDrpg5. The TIR and DR sequences of three IS elements with a 265 aa length Tpase were "AGGCTG" and "TAG, " respectively. The TIR and DR sequences of QR90_RS06350 were "ACCTCCTGCGAAAGTC" and "TAG, " respectively. The structural schemes for the forward and reverse area of the Tpase are identical, with a 54 and 9 nt distance, respectively. One interesting finding is that the loci of the IS5 family IS elements are positioned close to the IS701 family member, except for QR_ RS06350.
In this study, we identified five unclassified IS members: QR90_08760 with a 270 aa-long Tpase; QR90_RS08625 with a 477 aa-long Tpase; QR90_08735 with a 434 aa-long Tpase; QR90_RS05950 with a 120 aa-long Tpase; and QR90_04880 with an 89 aa-long Tpase.
Next, the active transposition of IS elements was detected on non-pigment phenotypic selection under oxidative stress conditions induced by H 2 O 2 treatment in the present study.

Detection of active transposition on carotenoid biosynthesis by oxidative stress
A total of 24 and one non-pigmented colonies were detected after low concentration (80 and 100 mM) and high concentration (200 and 300 mM) H 2 O 2 treatment for 1 h, respectively. It seems that generation of non-pigmented mutants is less efficient with higher H 2 O 2 concentration.
In the high-concentration treatment, one is a complete non-pigmented strain (w3) and two isolates exhibited a pale reddish color on TGY plates (w1 and w2; Figure 3A). When carotenoid biosynthesis genes were amplified by PCR with the appropriate primers, the complete non-pigmented strain only exhibited an enlarged PCR product from QR90_10400, a phytoene desaturase ( Figure 3B). However, other carotenoid biosynthesis genes were not affected. The ISDrpg3 of the IS4 family was integrated at the 275th nucleotide of QR90_10400 (Figure 4). The TIR sequence of this IS element was "CTCGGTAGCTGACAACTTCA" and the DR sequence was "ACCCGCCCC. " In the low-concentration H 2 O 2 treatment, five out of 24 non-pigmented strains exhibited PCR amplicons with QR90_10400 gene enlargement: two strains with an OD 600 of 2.0 and three strains with an OD 600 of 4.0 at 80 mM H 2 O 2 ( Figure 5A). The frequency of IS transposition in gene QR90_10400 exhibited 5.1 × 10 −4 and 1.8 × 10 −4 from OD 600 of 2.0 and 4.0, respectively (Supplementary Table S4). All five of these IS-integrated mutants exhibited IS4 family transposition ( Figure 5B; Supplementary Figure S4). Other carotenoid-related genes QR90_03795, QR90_14380, and QR90_14400 were not affected ( Figure 5C). The 19 remaining clones did not exhibit any changes in the size of the PCR products, indicating that there was no IS transposition in the selected four carotenoid biosynthesis genes. However, there might be point mutations in carotenoid-related genes or perhaps IS transposition in non-analyzed carotenoid biosynthesisrelated genes.
At an OD 600 of 2.0, one IS element (w3) was integrated in the forward direction at the 401st nt of phytoene desaturase. The DR sequence of this element was "CTTCTTCGA" and the TIR sequence was "CTCGGTAGCTGACAACTTCA. " Another (w7) integrated in the reverse direction at the 617th nt. The DR sequence of this element was "GTAAACGAG" and the TIR sequence was "CTCGGTA GCTGACAACTTCA" (Figure 5B; Supplementary Figure S4A). Both IS elements were belonged to ISDrpg3 of QR90_04150 or QR90_ RS10860 (Supplementary Table S2).
At an OD 600 of 4.0, one IS element (w5) integrated in the reverse direction at the 587th nt of phytoene desaturase. The DR sequence of this element was "AGGCGCTC" and the TIR sequence was "CTCGGTAGCTGACAACTTCA. " Another (w11) integrated in the reverse direction at the 474th nt; its DR sequence was "GCTCGTAGC" and its TIR sequence was identical to that of w5. Both IS elements were identical to ISDrpg3 (Supplementary Figure S4B). The last IS element (w7) integrated in the forward direction at the 22nd nt with Simplified schematic of the carotenoid biosynthesis pathway in DY59 strain from KEGG carotenoid pathway analysis.
DR sequence "CCAGCAGGC" and TIR sequence "CTCT GTACCGGACAACT" (Figure 5B; Supplementary Figure S4C). This IS element was identical to ISDrpg2 of QR90_01215, QR90_04660, or QR90_09690. Therefore, D. radiopugnans DY59 exhibits active transposition of IS4 family members in H 2 O 2 treatment conditions. The PCR detection of five copies of IS4 family members at the location as found in the genome sequence indicates that the active transposition occurred through replicative mode in present (Supplementary Figure S5).

Effects of hypochlorite and gamma-irradiation in IS transposition
The 100 μM sodium hypochlorite does not affected viability of the strain DY59. 5 kGy irradiated strain DY59 exhibits 99.5% reduction of CFU. The selection of non-pigmented colonies was performed via gamma irradiation exposure and sodium hypochlorite treatment. Unlike the lack of pigment production of D. geothermalis, gamma irradiation of total 5 kGy did not induce non-pigment phenotypes in A B

FIGURE 3
Phenotypic diversity in the high concentration of H 2 O 2 treatment (200 mM for w1 and w2 or 300 mM for w3) for 1 h (A) and detection of gene disruption by PCR of genes involved in the carotenoid pathway (B). QR90_03795, Phytoene synthesis; QR90_14400, FAD-dependent oxidoreductase; QR90_14380, carotenoid 1,2-hydratase; QR90_10400, phytoene dehydrogenase. Lanes: M, size marker; 1,5,9: WT;2,6,10: w1;3,7,11: w2;4,8,12: w3. The dotted arrow indicates the wild-type gene PCR products. The solid arrow indicates the IS-integrated PCR product. Detection of IS integration loci and composition of IS element in QR90_10400 phytoene desaturase in w3 mutant obtained after 300 mM H 2 O 2 treatment. The ISDrpg3 of the IS4 family member was integrated at the 275th nt locus of the QR90_10400 gene with reverse transcriptional direction.
Frontiers in Microbiology 07 frontiersin.org the D. radiopugnans DY59 wild-type strain . Moreover, although hypochlorite treatment induced the non-pigment phenotype, the non-pigmented colony did not exhibit IS element integration in the four analyzed carotenoid biosynthesis genes (data not shown).
This may be the same explanation for the absence of transposition of IS in phenotypic changes due to point mutations in four carotenoid biosynthetic genes and the decay of other genes related to pigment formation. Frontiers in Microbiology 08 frontiersin.org

Expression levels of catalase, three LysR family regulators included oxyR and Tpase of the IS4 family by qRT-PCR analysis
To determine the mechanisms underlying the high H 2 O 2 resistance of DY59, the expression levels of catalase and three LysR family regulators including possible oxyR were first measured by qRT-PCR analysis. The chromosome of DY59 strain has a single catalase QR90_06310 with 72.17 and 74.17% amino acid sequence similarity to D. radiodurans catalase KatE1 and D. geothermalis KatE, respectively. After exposure to various concentrations of H 2 O 2 (50, 100, 150 mM, and an unexposed control) for 1 h and growing the cells to OD 600 values of 2.0 and 4.0, the relative expression levels of catalase and LysR family members were measured via qRT-PCR by the basal expression level of OD 600 of 2.0 with unexposed control as 1.0. Unexpectedly, catalase was not dramatically induced at OD 600 2.0 or 4.0 in any of the tested H 2 O 2 concentrations. Nevertheless, there was a 2-fold increase in catalase expression at OD 600 4.0 when the cells were challenged with 50 and 100 mM H 2 O 2 ( Figure 6A). LysR1 QR90_13110 exhibited a more than 6-fold induction at an OD 600 of 4.0 upon exposure to 100 mM H 2 O 2 and LysR3 QR90_15105 exhibited a more than 6-fold induction at an OD 600 of 2.0 in the 100 mM H 2 O 2 condition. Our findings thus indicated that 100 mM H 2 O 2 controlled the expression of LysR1 and LysR3 at two different growth phases ( Figure 6B). However, LysR2 (QR90_14595) was not affected regardless of growth phase and oxidation condition (data not shown). LysR2 was exhibited 70.3 and 75.7% identities of amino acid sequence to the proposed OxyR of D. radiodurans and D. geothermalis, respectively. Therefore, the high H 2 O 2 resistance phenotype of the DY59 strain cannot be attributed to changes in the expression of catalase and its three potential LysR family regulator members, as determined by RNA production levels measured by qRT-PCR analysis.
There are two subtypes of IS4 family which were actively transposed to other loci. Tpase of ISDrpg3 such as QR90_04150 and QR90_10860/ ISDrpg2 such as QR90_01215, QR90_04660, and QR90_09690 exhibited a more than 5-6-fold induction at an OD 600 of 4.0 in the 50 mM H 2 O 2 condition. The related expression levels were reduced at an OD 600 of 4.0 in the 100 and 150 mM H 2 O 2 condition ( Figure 6C). However, both ISDrpg2 and ISDrpg3 Tpases were gradually induced over 2-fold at an OD 600 of 2.0 in the different H 2 O 2 conditions. Thus, IS elements of the IS4 family were induced by the H 2 O 2 and actively transposed into other genomic loci with replicated mode.

Discussion
Bacillus subtilis, a widely known Gram-positive bacterial model organism, lacks IS or any other transposable element excluding several bacteriophages and other remnants of horizontal gene transfer (HGT) events from genome study of wild-type strain. The IS identification platform "ISfinder" revealed four IS families in the B. subtilis genome including IS4Bsu1 (IS4 family), ISBsu1 (IS3 family), ISBsu2 (IS256 family), and ISBsu3 (IS1595 family; Siguier et al., 2006). Interestingly, pathogenic Bacillus species such as B. cereus, B. anthracis, B. thuringiensis, and others have many IS families and Tn3 elements (Fayad et al., 2019). Therefore, the distribution of transposable elements may be associated with pathogenicity. Additionally, the genomes of the pathogenic bacteria may have undergone multiple acquisition of transposable elements through HGT in response to environmental stimuli (Fayad et al., 2019).
Here, we detected active transposition of IS elements under oxidative stress conditions of H 2 O 2 treatment in the radiation-resistant bacterium D. radiopugnans DY59, which was isolated from mountain soil collected in South Korea. The DY59 isolate was phylogenetically clustered near three Deinococcus species from Antarctic marine environments (Hirsch et al., 2004;Kim et al., 2015). The genome of D. radiopugnans DY59 has a total of 36 IS elements encompassing six IS families. Interestingly, these six family members belonged to D. geothermalis IS families that also contained the IS6, IS982, and IS200/605 families . Unfortunately, the IS names cannot be classified based on the IS distribution in the D. radiopugnans genome from "ISfinder" platform (at Sep. 2022;Siguier et al., 2006). The transposase annotation as a key component of IS element was varied among the genome data and between prediction and analysis.
When wild-type DY59 cells were treated with 80 or 300 mM H 2 O 2 , the active transposition of the IS4 family members was only detected on phytoene desaturase (QR90_10400; Figures 4, 5). The active transposition of the IS4 family was also particularly detected in wildtype D. geothermalis upon H 2 O 2 treatment and dielectric bilayer discharge (DBD) plasma radiation Ye et al., 2022).
When the D. radiodurans strains were exposed to gamma irradiation, there were two particular active transposition events on a trimethoprim-resistant selection. One was an IS integrated into a thyA gene; the other was an uvrA gene disruption that resulted in mitomycin-resistant phenotypic selection (Narumi et al., 1997;Mennecier et al., 2006;Pasternak et al., 2010). When the cells were treated with 5 kGy of gamma irradiation, the wild-type D. geothermalis strain exhibited active transposition of IS families, (e.g., IS1 and IS5 family), whereas the wild-type DY59 strain and D. radiodurans strain did not exhibit non-pigment phenotypic mutations .
Insertion sequence transposition has been detected under various environmental stressors, such as nutrient deprivation, temperature changes, metal ion exposure, and oxidative stress caused by UV irradiation, gamma irradiation, and H 2 O 2 treatment (Ohtsubo et al., 2005;Twiss et al., 2005;Kharat et al., 2006;Mijnendonckx et al., 2011;Suzuki et al., 2021). Interestingly, this active transposition of IS elements was found to vary in a species-specific manner. D. geothermalis wildtype and several particular gene-disrupted mutants have been reported to exhibit different types of IS transposition. For example, ISDge3 of the IS1 family was actively transposed under gamma irradiation, whereas ISDge11 of the IS4 family was transposed under H 2 O 2 treatment and DBD plasma radiation. Moreover, ISDge5 and ISDge6 of the IS5 family were found to be transposed in dps-, oxyR-, cystine importer-, and lysRdeficient mutants Ye et al., 2021).
Here, we explored the IS distribution in the genome size of 3.54 Mb of the DY59 strain. Our findings suggested that DY59 exhibits a less complex IS distribution than that of D. geothermalis. In this study, a strict IS element of the IS4 family was actively transposed into a carotenoid biosynthesis gene QR90_10400 phytoene desaturase under H 2 O 2 treatment (Figures 4, 5B; Supplementary Figure S4). Therefore, the wild-type D. radiopugnans DY59 may serve as a suitable model organism for studying the active transposition of unique IS element using a single oxidation inducer such as H 2 O 2 . These experiments provide the opportunity to determine unique IS transposition events in organisms that exhibit genomic plasticity. However, there are still many challenges associated with IS naming Frontiers in Microbiology 09 frontiersin.org and classification from gene annotation, emphasizing the need for more accurate algorithms and additional criteria including machine learning tools for IS identification and assignment, as well as the creation of a global network of research groups working together. The genome sequence information of D. radiopugnans ATCC19172 was updated twice on June 2019 and August 2020, with contig assembly lengths of 4.33 and 4.3 Mb, respectively (NCBI genomes). Deinococcus radiopugnans ATCC19172 has three catalases, FHR04_11220, FHR04_17100, and FHR04_17320, sharing, respectively, 99, 26.6, and 47.6% identity with the single identified catalase (QR90_06310) from strain DY59. In case the published genome sequence of DY59 is incomplete, this strain might Comparison of the gene expression levels of catalase (A), two LysR family regulators (B) and two transposases in ISDrpg2 and ISDrpg3 (C) by qRT-PCR under different final concentrations of H 2 O 2 (0-150 mM). Pair-wise comparisons between experimental groups were conducted via Student's t-test with the Prism™ ver. 8.0 software (*p < 0.05; **p < 0.01; ***p < 0.001).
Frontiers in Microbiology 10 frontiersin.org also possess homologs of FHR04_17100 and FHR04_17320. Thus, we performed qRT-PCR using primer sets corresponding to FHR04_17100 and FHR04_17320 to investigate expression of possible homologs in strain DY59. The possible FHR04_17100 homolog was not induced by the H 2 O 2 treatment and the possible FHR04_17320 homolog exhibited more than 2-fold induction at OD 600 of 2.0 with 50 and 100 mM H 2 O 2 treatment (Supplementary Figure S6). When 150 mM H 2 O 2 treatment at OD 600 of 4.0 was applied, the possible FHR04_17320 homolog exhibited nearby 3-fold induction. In the present additional catalase expression data, the high hydrogen peroxide resistance of DY59 is not explained by induced expression of catalase genes. DY59 strain may have a high level of constitutive catalase enzyme activity, or it may employ other mechanisms for hydrogen peroxide resistance, such as unidentified protective pathways and physiological defense systems using extracellular matrixe components such as extracellular polymeric substances (EPS), proteins, and eDNA, which aid in biofilm formation and enhance tolerance to oxidative stress, and certain transporters (Li et al., 2013;Molina-Santiago et al., 2021). The genome sequence of strain ATCC19172 is 0.8 Mb larger than that of DY59. Although the reported genome sequence of strain DY59 consist of only a chromosome, we consider the possibility that this strain might have one or more additional genome molecules such as plasmids. If the strain DY59 has plasmids, the IS family and number will expand. IS elements are commonly known to exhibit random movement. However, experiments conducted on Deinococcus geothermalis have shown that the transposition of unique IS elements follows a specific pattern. The exact pattern is still unknown, but this experiment provides a starting point for further research. Specifically, it was observed that the IS4 family of IS element only transposed under hydrogen peroxide conditions in D. radiopugnans DY59. Additionally, the ability of Deinococcus species to endure various forms of stress such as radiation, drying, and toxic chemicals is attempted to be explained from a molecular evolutionary perspective through the transposition of IS elements. Further exploration into this research area presents a significant challenge for the future.

Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number (s) can be found in the article/Supplementary material.