Next Article in Journal
Association between Lipoprotein Levels and Humoral Reactivity to Mycobacterium avium subsp. paratuberculosis in Multiple Sclerosis, Type 1 Diabetes Mellitus and Rheumatoid Arthritis
Previous Article in Journal
Genomic Insights into the Fungal Lignocellulolytic Machinery of Flammulina rossica
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

sRNA OsiA Stabilizes Catalase mRNA during Oxidative Stress Response of Deincoccus radiodurans R1

1
Department of Plant Science, School of Agriculture and Biology, Shanghai Jiao Tong University, Shanghai 200240, China
2
Biotechnology Research Institute, Chinese Academy of Agricultural Sciences, Beijing 100081, China
*
Authors to whom correspondence should be addressed.
Microorganisms 2019, 7(10), 422; https://doi.org/10.3390/microorganisms7100422
Submission received: 3 September 2019 / Revised: 27 September 2019 / Accepted: 30 September 2019 / Published: 8 October 2019
(This article belongs to the Section Environmental Microbiology)

Abstract

:
Deinococcus radiodurans adapts to challenging environments by modulating gene expression in response to oxidative stress. Recently, bacterial small noncoding RNAs (sRNAs) have been presumed to participate in the transcriptional or translational regulation of stress-responsive genes. We found 24 sRNAs that may be involved in the oxidative stress response of D. radiodurans by deep RNA sequencing. Moreover, a typical stress-inducible sRNA, IGR_3053, named OsiA, was predicted to bind to the mRNA of katA, katE, and sodC by the bioinformatics method. An osiA knockout of D. radiodurans displayed increased sensitivity to H2O2 and the decreased catalase activity and total antioxidant activity, suggesting that OsiA probably serves as a regulator in the adaptation to oxidative environments. Further microscale thermophoresis results demonstrated that OsiA can directly bind to the mRNA of katA, sodC, and katE. The stability test result of katA mRNA showed that its half-life was 2 min in the osiA mutant compared with 5 min in the wildtype(wt) strain. Our results indicated that OsiA can enhance the stability of katA mRNA and the activity of KatA and consequently the oxidation resistance of D.radiodurans. We are the first one to explore the super-strong oxidative stress resistance of D.radiodurans at the level of post-transcriptional regulation, and found a new pathway that provides a new explanation for the long-term adaptability of D.radiodurans in extreme environments.

1. Introduction

Deinococcus bacteria are extremophiles known for their superior resistance to ionizing radiation and oxidative stresses [1,2,3,4,5]. Decades of studies have indicated that their extreme resistance is attributed to the protection of oxidized protein degradation and activation of multiple factors and mechanisms for repairing DNA damage [6,7,8,9]. In the DNA repair system of Deinococcus, RecA(Recombinase A), RecFOR system, DNA glucoamylase, and other Deinococcus-specific proteins such as DdrB(DNA damage response protein) and PprA(DNA damage repair protein) show higher repair activity than Escherichia coli in maintaining genomic stability [10]. When facing oxidative stresses, Deinococcus bacteria have developed a more efficient ROS(Reactive Oxygen Species)-scavenging system to defend against oxidative stress in comparison with other stress-sensitive bacteria [5]. This well-equipped system includes catalases (KatA, KatE, SodC, etc.), antioxidant defense proteins (MsrA, MsrB, etc.), specific metabolites such as carotene and carotenoids [9,10,11,12,13], and a significantly higher Mn/Fe concentration ratio. In addition, some stress response transcriptional factors like CarD(RNA polymerase-interacting regulator), DrRRA(DNA-binding response regulator A), IrrE(Ionizing radiation resistance protein), OxyR have also been reported to play significant roles in the stress resistance of D. radiodurans. Aside from these factors, Deinococcus bacteria have developed a conserved regulon, DdrO and PprI (also named as IrrE), that is involved in transcriptional regulation during the oxidative response [14,15]. DdrO and PprI serve as a global regulator pair that regulates a vast array of genes involved in diverse cellular processes under stress conditions, indicating the complexity of the stress tolerance mechanism in Deinococcus [10]. This has prompted us to investigate whether there are also other factors that cooperatively regulate the expression of dozens of genes in Deinococcus radiodurans.
Oxidative stress reaction becomes the secondary of all stress responses on account of excessive ROS [16]. Thus, the key to improving the biological resistance to adverse environments such as high temperature, high salt, drought, and radiation is to remove the excessive ROS. Catalases are major players among the ROS scavenging pathway [16]. The three members, Kat, Sod, and Pod, work in different divisions; regarding Kat, it is capable of reducing peroxide radicals to molecular oxygen and water; superoxide radicals are reduced to peroxide by Sod; and the role of Pod is peroxide detoxification [17]. So far, substantive investigations have been performed with regard to their post-transcriptional regulation and there is still a lot of space left for their regulation among diverse metabolic pathways. A noncoding RNA RyhB has previously been reported to degrade mRNAs encoding iron binding proteins like the mRNA of sodB, which encodes Fe-superoxide dismutase [18]. Nfis, a previously reported sRNA involved in the nitrogen fixation of Pseudomonas stutzeri A1501 were proven to regulate oxidative stress response by binding to the mRNA of the catalase coding gene katB [19,20]. However, there have been no reports on the post-transcriptional regulation of catalases in D.radiodurans. Current studies have focused on DNA and protein levels, few reports have been reported on the regulation of the post-transcriptional level. There is still a lack of adequate explanation for the super extreme environmental adaptability of D.radiodurans.
Small noncoding RNAs (sRNAs) in bacteria are diverse in sequence, and are perceived as crucial elements in gene regulation [21]. sRNAs usually have a length of 50–500 nt, and most regulate mRNA stability or translation efficiency through base-pairing [22]. sRNAs may be classified into three categories according to the location of gene targeting: trans-encoded, cis-internal, and intergenic sRNAs (in a noncoding region between two genes) [23]. The majority of bacterial sRNAs are trans-encoded sRNAs that serve as important regulatory factors [24]. Regulatory sRNAs are frequently induced in stress responses, and they bind to mRNA of the aforementioned genes encoding the transcription factors and proteins to some extent [25], which together contribute to the ability of cells to cope with specific growth conditions such as irradiative or drought conditions. Few studies have investigated the sRNA regulation of D.radiodurans. A total of 144 noncoding RNAs were detected during gamma-irradiated phases of D.radiodurans and 199 noncoding transcripts were detected in the D.radiodurans genome under normal condition; nevertheless, their role in stress resistance has not been studied in depth [26].
A previous study determined the transcriptome of D. radiodurans under 20 mM H2O2 treated for 30 min [13]; interestingly, we found that a relatively obvious phenotype appeared after 80 mM H2O2 treated for 30 min, therefore we determined the transcriptome of D. radiodurans under such conditions. Combined with the transcriptome data with bioinformatics analysis, one sRNA attracted our attention for its binding to the mRNAs of the KatA, KatE, and SodC coding genes. It is speculated that it may be involved in several metabolic pathways of the oxidative stress response of D. radiodurans, so the biological function was studied in detail.

2. Materials and Methods

2.1. Bacteria Strain, Growth Conditions, and Oxidative Stress Treatment

D. radiodurans cells were purchased from China General Microbiological Culture Collection Center(CGMCC, Beijing, China) and were cultured in TGY medium (1% tryptone, 0.5% yeast extract, and 0.1% glucose, pH 5.3, OXIOD, Hampshire, United Kingdom) at 30 °C and shaken at 220 rpm. When the OD (optical density) of the samples reached 0.8 (7 × 106 CFU/mL approximately), the culture was treated with 80 mM H2O2. After 30 min, the bacteria cells were centrifuged at 7000 rpm for 10 min, and the pellets were washed twice with PBS buffer (Phosphate Buffer Saline, 0.02% KH2PO4, 0.29% Na2HPO4·12H2O, 0.8% NaCl, 0.02% KCl, pH 7.5) and finally centrifuged at 12,000× g for 1 min. An untreated cell culture grown under the same condition was used as the control.

2.2. Total RNA Extraction

Total RNA was extracted according to the protocol of the PureLinkTM RNA Mini Kit (Invitrogen, California, USA) and purified decontaminated using a genomic DNA eraser (Tigan, Beijing, China). The RNA concentration was measured using a Qubit RNA Assay Kit in a Qubit 3.0 Fluorometer (Life Technologies, California, USA), and RNA quality and concentration were assessed using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2100 system (Agilent Technologies, California, USA).The RNA sample was diluted to 1000 ng/µL for RNA-Seq.

2.3. Candidate sRNA Screen and the Prediction of sRNA Target Genes

The raw data from the Illumina sequencing platform (https://www.illumina.com/ Illumina/sequencing and array-based solutions for genetic research, David Walt and Larry Bock, April, 1998) [27] was preprocessed to remove the adapter sequences and filter out the low-quality data for QC purposes. The reference D. radiodurans R1 genome and gene annotations were downloaded directly from the NCBI FTP website [28]. The sRNA screen was performed according to the method described by Zhouwei Chen et al. When the ratio between the MEV (mean expression value) of transcripts in the intergenetic region and their adjacent CDS (Coding sequence) was greater than 2, the transcripts were considered to be an sRNA candidate [29]. The sRNA was further validated by PCR. The RNA interaction prediction tool IntaRNA (http://rna.informatik.uni-freiburg.de/IntaRNA/Input.jsp?jobID=4476547&reload=true IntaRNA-RNA-RNA interaction, Busch, 2008) [30] was used to predict sRNA targets online.

2.4. Gene Ontology (GO) and KEGG Enrichment Analysis

Gene ontology (GO) annotations (http://geneontology.org/ The Gene Ontology Resource, Michael Ashburner, May, 2000) [31] were conducted using the AgBaseGOanna package60 (http://agbase.msstate.edu/GOAnna.html Fiona M. McCarthy, 2006) with a BLASTx (https://blast.ncbi.nlm.nih.gov/Blast.cgi?PROGRAM=blastx&PAGE_TYPE=BlastSearch&LINK_LOC=blasthome Basic Local Alignment Search Tool, Jun, 2005) [32] search against the D. radiodurans R1 genome database. GO enrichment of differentially expressed genes was analyzed by using a Perl module (GO::TermFinder) (http://go.princeton.edu/cgi-bin/GOTermFinder Generic Gene Ontology Term Finder, Gavin Sherlock, Nov, 2009) [33]; KOBAS online database (KEGG Orthology Based Annotation System, http://kobas.cbi.pku.edu.cn/ Jianmin, July, 2006) [34]. It was concluded that GO items with corrected p values less than 0.05 were significantly enriched. The KOBAS online database was applied to enrich the DEGs. The KEGG (Kyoto Encyclopedia of Genes and Genomes, https://www.kegg.jp/ KEGG: Kyoto Encyclopedia of Genes and Genomes, Kanehisa Laboratories, 1995) [35] pathways with a modified p value less than 0.05 were considered to be significantly enriched in the tested pathways.

2.5. Microscale Thermophoresis Analysis

The microscale thermophoresis (MST) technique was used to quantify the interactions between sRNAs and the mRNAs transcribed by the target genes [36,37,38]. A T7 transcription kit was first used to synthesize the 5′FAM-labeled mRNA of the target genes (katA, katE, sodC, etc.) in vitro [34]. The labeled mRNAs were then diluted to the fluorescence range of 200–1000 according to the protocol in the user manual of the Monolith NT.115 (Nano temper company, MÜnchen, Germany).
sRNA (50 µg) was used as a template for in vitro transcription, and the resulting synthetic mRNA was dissolved in DPEC(diethylpyrocarbonate)-treated water to an OD = 2. sRNA (10 µL) was first added to a PCR tube, and then serially diluted with DPEC-treated water using a 1:1 ratio for the next 14 tubes. After that, 10 µL of labeled mRNA was mixed with the diluted sRNAs in each tube. After incubation for 30 s, the samples were loaded into MST NT.115 standard glass capillaries. The molecular interactions were tested and the Kd (dissociation equilibrium constant) values were calculated using the NanoTemper software (Nano Temper company, MÜnchen, Germany) package.

3. Results

3.1. Deep Sequencing sRNAs from D. radiodurans R1 under Oxidative Stress

To globally identify small noncoding RNAs differentially expressed in response to oxidative stress in D. radiodurans, we exposed at least three replicate cultures of D. radiodurans to 80 mM H2O2. Total RNA was isolated from untreated cultures (CDR) and H2O2-treated cultures (TDR). sRNAs were sequenced using a strand-specific size-selected (50–500 nt) sRNA library preparation.
After trimming the adapters and discarding the low-quality reads, there were about 6 Gbp reads for each sample to evaluate the gene expression and predict sRNAs. The valid data of CDR and TDR were 79.77% and 77.92%, respectively (Supplementary Table S1). The transcriptional reads were mapped to the D. radiodurans genome to identify sRNAs that are antisense to the structure genes. As shown in Supplementary Figure S1, sRNAs in D. radiodurans genome ranged from 50 to 500 bp, and most of the sRNAs were from 50 to 200 bp. As the length increased, the number of sRNAs was greatly reduced, indicating that sRNAs regulate the expression of target genes mainly via binding with short sequences less than 200 bp.
To confirm the transcriptomic data, we randomly selected 10 sRNAs to validate their expression levels under the oxidative condition by qRT-PCR (quantitative real-time PCR). As shown in Supplementary Figure S1, the expression of 10 sRNAs quickly responded to oxidative treatment and were in accordance with the RNA-Seq analysis.
To make a thorough inquiry of sRNA function in gene regulation, we analyzed the whole genome using computer prediction and experimental methods to identify transcriptional regulatory elements.
Further functional classification of target genes showed that the target genes of 97 sRNAs act as enzymes (42 upregulate and 55 downregulate), the target genes of 32 sRNAs act as kinases (15 upregulate and 17 downregulate), the target genes of 21 sRNAs act as transcription factors (seven upregulate and 14 downregulate), and the target genes of 22 sRNAs are structural proteins (nine upregulate and 13 downregulate) (Table 1). This result indicates that sRNAs can regulate the expression of both structural genes and regulators (kinases and transcription factors) at the transcriptional level. What we found very amusing about this result is that one sRNA may regulate more than one target gene with diverse functions.

3.2. 24 Oxidative Related sRNAs Were Identified by KEGG Analysis

To analyze the function of target genes during oxidative responses, we used the KOBAS database to enrich the KEGG pathway that is probably regulated by sRNAs. The enriched results showed that the target genes regulated by 24 sRNAs were involved in six oxidative stress response relevant metabolism pathways (p < 0.05) including peroxisome function, DNA mismatch repair, glutathione metabolism, carotenoid biosynthesis, the two-component system, and vitamin B6 metabolism (Table 2, Supplementary Tables S2 and S3). It also showed that the enriched pathways were regulated by multiple sRNA regulators. For example, seven sRNAs are involved in the peroxisome pathway, 15 sRNAs are involved in the DNA mismatch repair pathway, and eight sRNAs are involved in glutathione metabolism. Hence, these results indicate that oxidative stress-inducible sRNAs have regulatory functions in the expression of the genes related to robust repair systems and specialized metabolism in D. radiodurans, and the cooperation of different repair systems or stress signaling pathways may be tuned through the same sRNA.
It was noted that one sRNA, OsiA (Oxidative stress-induced sRNA A), which is located in the 358,768 to 358,969 region of chromosome 2, can target five genes related to six pathways including peroxisome, mismatch repair, and two-component system (Table 2). In the peroxisome functional pathway, three genes coding KatA (DR_1998), SodC (DR_A0202), and KatE (DR_A0259) were the target genes regulated by OsiA. Furthermore, the gene coding RplV(DR_0316) in mismatch repair and the gene coding DR_A0354 are probably also regulated by OsiA. These results indicate that OsiA may play a crucial role in transcriptional regulation in D. radiodurans upon oxidative stress.

3.3. Characterization of OsiA, an Oxidative Response sRNA

Northern blotting showed that the expression of OsiA was increased in TDR compared with CDR, which further verified the reliability of the transcriptome data (Figure 1A). qRT-PCR showed that OsiA expression was first increased and later decreased with an increase in H2O2 concentration, and the maximum expression was reached at 40 mM, which suggested that OsiA responded to hydrogen peroxide shock (Figure 1B). To clarify the function of OsiA in the oxidative stress response, we knocked out osiA in D. raidodurans and analyzed the oxidative phenotype of this mutant. Under normal conditions, there was no difference between the WT, ∆osiA, and osiA complementary strain. With the increase of H2O2 concentration, ∆osiA was more sensitive than WT (Figure 1C), and the survival curve demonstrated that ∆osiA could not survive with a high concentration of H2O2 (Figure 1D), while the phenotype and survival fraction of osiA mutant could be recovered after supplementing osiA in the mutant strain, which indicates that OsiA is necessary for the growth of D. radiodurans under oxidative conditions.
To further quantify the oxidative resistance of D. radiodurans in the absence of osiA, catalase activity (Figure 2A) and total antioxidant capacity (Figure 2B) in wild type, osiA mutant, and osiA complementary strains were measured. The results revealed that both catalase activity and total antioxidant capacity were significantly reduced and the compensatory strain can compensate for this result to some extent.

3.4. OsiA Enhances the Stability of Transcripts by Binding to the Target Genes through Multiple Bases

Interestingly, the bioinformatics analysis showed that the binding sequences of OsiA and the transcripts of the three target genes were all located in the same stem-loop (nucleotides 162–175) of OsiA (Supplementary Figure S3).
To investigate whether OsiA binds directly to the target genes, microscale thermophoresis was applied to identify the binding strength between OsiA or its mutants and the corresponding complementary fragment in the mRNA of target genes (katA, katE, sodC). Among the tested genes, katA mRNA had the strongest binding force with OsiA. The interaction of katE mRNA and OsiA was ranked second, while sodC mRNA had the weakest binding with OsiA. According to the predicted binding strength in different combinations, we mutated different bases in OsiA and then assayed the change in binding strength (Figure 3A,C,E). As expected, the binding strength between the mutated OsiA variants and their target genes was completely disrupted (Figure 3B,D,F), suggesting that OsiA strongly binds to the target genes through different bases. Using katA mRNA as an example, when the binding sequence ACACCCGCCCCGAC of OsiA was replaced by the mutated sequence UGUCGGGCCCGCUG, no binding could be detected between OsiA and katA mRNA.
The binding assay showed that OsiA can regulate multiple genes in different cellular antioxidant pathways, indicating that OsiA is required for adaptation to extreme environments and probably has specific regulatory functions in D. radiodurans.
The half-life of katA mRNA in the wildtype strain was 5 min, but this decreased to 2 min with osiA knocked out under the H2O2 treated condition (Figure 4). This result indicates that OsiA regulated katA mRNA at the posttranscriptional level and enhanced the stability of katA mRNA.

4. Discussion

D. radiodurans is a desiccation- and radiation-tolerant bacterium [2,3,4,5]. In recent years, many studies have examined the molecular mechanism of D.radiodurans in adapting to extreme environments [6,7,8,9]. High-resolution RNA-Seq analyses have indicated that the bacterial regulatory factors include transcription factors, protein kinases, and noncoding RNAs. When facing unfavorable environmental conditions, D. radiodurans R1 effectively activates the DNA repair system to protect DNA from breaking, or enhances the expression of antioxidants to protect proteins from degradation [1,10]. In particular, sRNA-mediated transcriptional regulation is likely to finely regulate stress responses to environmental challenges such as oxidative stress.
Peroxisome is a multifunctional organelle involved in the generation and decomposition of ROS, and carotenoid deinoxanthin is an important chemical involved in resistance against extreme oxidative stresses in D. radiodurans [10,11]. Furthermore, a growing number of studies show that two-component systems play an irreplaceable role in the sensing of environmental changes and adapting to adverse environments [39,40]. Vitamin, a kind of small molecular substance, can assist peroxidase to resist UV damage, reduce ROS content, and enhance oxidative stress resistance [41,42]. We identified 24 sRNAs involved in these six metabolic pathways and we will lucubrate them to fully understand how they participate in the oxidative stress and even defense responses of D.radiodurans, and finally construct the post-transcriptional regulatory network of D.radiodurans.
As previously reported, there are many transcriptional regulators involved in the abiotic stress of D.radiodurans such as DrRRA, OxyR, and IrrE. DrRRA was reported to regulate KatA and SodC, and both DrRRA and OxyR have been reported to regulate KatE in oxidative stress response. As shown in our results, OsiA interacts with the mRNA of katA, katE, and sodC, so we can infer that OsiA may have similar functions in global regulation.
The upstream gene (dr_A0243) of osiA encodes a flavohemoprotein that utilizes O2 and NAD(P)H to convert NO into nitrates [43]. The downstream gene (dr_A0244) of osiA encodes a protein that is very similar to the radical S-adenosyl-L-methionine protein, and acts against a variety of toxic agents caused by oxidative stress [44,45]. The remediating functions of the two genes, along with the target gene analysis of OsiA, suggest that OsiA is a pivotal sRNA during the oxidative response in D. radiodurans.
The interaction results showed that all target genes shared complementary sequences with OsiA at 162–175 nt (Supplementary Figure S2B), while it can be seen from the binding analysis that OsiA targets genes by different recognition sites (Figure 3). Future studies are needed to mutate the different sites in OsiA and explore its other targets to elucidate the detailed function of OsiA and how it crosstalks with other regulators during the adaptation of extremophiles to extreme environments.
A series of reports have focused on the extreme resistance mechanism of D.radiodurans. So far, its extreme resistance can be attributed to the following: its own genome structure, antioxidant enzymes, stress protective protein, carotene biosynthesis, high Mn/Fe ratio, and transcriptional regulator. No one has reported on the post-transcriptional regulation of these resistance genes. This study was conducted to explore the post-transcriptional regulation of D.radiodurans and reveal its long term extreme environmental adaptability from a new perspective. In the future, we will further study the post-transcriptional regulation of D.radiodurans, adding a strong note to the blueprint of its complex extreme environmental resistance mechanism.

Supplementary Materials

Supplementary Materials can be found at https://www.mdpi.com/2076-2607/7/10/422/s1. Table S1: Small RNA sequencing and classification in D. radiodurans R1, Table S2: Information of 24 sRNAs involved in oxidative stress response, Table S3: KEGG pathway analysis of oxidative response sRNAs in D. radiodurans R1, Figure S1: The sRNA length distribution in D. radiodurans R1. The sRNA length (nt) is shown on the x-axis and the percentage of sRNA quantity is represented on the y-axis, Figure S2: qRT-PCR validation of 10 sRNAs identified by RNA-seq. Vertical values represent fold change (log2) of small noncoding RNAs (sRNAs) at the transcriptional level after H2O2 treatment, Figure S3: The deduced second structure of OsiA. The black line in the figure represents the common areas on OsiA where each target gene is combined. The color legends represent base pair probabilities.

Author Contributions

Conceptualization, Y.C., K.Z., M.L., J.W.; Methodology, J.H., J.L., R.G., K.Z.; Validation, Y.C., D.X.; Formal analysis, D.X., W.S.; Investigation, Y.C.; Writing—original draft preparation, Y.C.; Writing—review and editing, M.L.; Supervision, Z.Z., W.Z., M.C.

Funding

This work was funded by the Key Research and Development Projects (No. 2018YFA0901003), the Ministry of Agriculture Transgenic Program (No.2019ZX08010-002), and the National Natural Science Foundation of China (Nos.31800061, 31570080, and 31500063). We also appreciate the support of the Agricultural Science and Technology Innovation Program of CAAS.

Acknowledgments

Special thanks go to Min Lin, Kaijing Zuo, and Zhengfu Zhou for their constructive suggestions on the revision of the article and other authors for their excellent technical support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cabiscol, E.; Tamarit, J.; Ros, J. Oxidative stress in bacteria and protein damage by reactive oxygen species. Int. Microbiol. 2000, 3, 3–8. [Google Scholar]
  2. Singh, O.V.; Gabani, P. Extremophiles: Radiation resistance microbial reserves and therapeutic implications. J. Appl. Microbiol. 2011, 110, 851–861. [Google Scholar] [CrossRef] [PubMed]
  3. Li, J.; Li, Q.; Ma, X.; Tian, B.; Li, T.; Yu, J.; Dai, S.; Weng, Y.; Hua, Y. Biosynthesis of gold nanoparticles by the extreme bacterium Deinococcus radiodurans and an evaluation of their antibacterial properties. Int. J. Nanomed. 2016, 11, 5931–5944. [Google Scholar] [CrossRef] [PubMed]
  4. White, O.; Eisen, J.A.; Heidelberg, J.F.; Hickey, E.K.; Peterson, J.D.; Dodson, R.J.; Haft, D.H.; Gwinn, M.L.; Nelson, W.C.; Richardson, D.L.; et al. Genome sequence of the radioresistant bacterium Deinococcus radiodurans R1. Science 1999, 286, 1571–1577. [Google Scholar] [CrossRef] [PubMed]
  5. Cox, M.M.; Battista, J.R. Deinococcus radiodurans - The consummate survivor. Nat. Rev. Microbiol. 2005, 3, 882–892. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, C.; Wang, L.; Li, T.; Lin, L.; Dai, S.; Tian, B.; Hua, Y. A PerR-like protein involved in response to oxidative stress in the extreme bacterium Deinococcus radiodurans. Biochem. Biophys. Res. Commun. 2014, 450, 575–580. [Google Scholar] [CrossRef] [PubMed]
  7. Lin, L.; Dai, S.; Tian, B.; Li, T.; Yu, J.; Liu, C.; Wang, L.; Xu, H.; Zhao, Y.; Hua, Y. DqsIR quorum sensing-mediated gene regulation of the extremophilic bacterium Deinococcus radiodurans in response to oxidative stress. Mol. Microbiol. 2016, 100, 527–541. [Google Scholar] [CrossRef] [PubMed]
  8. Li, M.; Sun, H.; Feng, Q.; Lu, H.; Zhao, Y.; Zhang, H.; Xu, X.; Jiao, J.; Wang, L.; Hua, Y. Extracellular dGMP Enhances Deinococcus radiodurans Tolerance to Oxidative Stress. PLoS ONE 2013, 8, e54420. [Google Scholar] [CrossRef]
  9. Ghosal, D.; Omelchenko, M.V.; Gaidamakova, E.K.; Matrosova, V.Y.; Vasilenko, A.; Venkateswaran, A.; Zhai, M.; Kostandarithes, H.M.; Brim, H.; Makarova, K.S.; et al. How radiation kills cells: Survival of Deinococcus radiodurans and Shewanella oneidensis under oxidative stress. FEMS Microbiol. Rev. 2005, 29, 361–375. [Google Scholar]
  10. Slade, D.; Radman, M. Oxidative Stress Resistance in Deinococcus radiodurans. Microbiol. Mol. Biol. Rev. 2011, 75, 133–191. [Google Scholar] [CrossRef]
  11. Tian, B.; Sun, Z.; Shen, S.; Wang, H.; Jiao, J.; Wang, L.; Hu, Y.; Hua, Y. Effects of carotenoids from Deinococcus radiodurans on protein oxidation. Lett. Appl. Microbiol. 2009, 49, 689–694. [Google Scholar] [CrossRef] [PubMed]
  12. Lu, H.; Chen, H.; Xu, G.; Shah, A.M.U.H.; Hua, Y. DNA binding is essential for PprI function in response to radiation damage in Deinococcus radiodurans. DNA Repair (Amst.) 2012, 11, 139–145. [Google Scholar] [CrossRef] [PubMed]
  13. Chen, H.; Huang, L.; Hua, X.; Yin, L.; Hu, Y.; Wang, C.; Chen, W.; Yu, X.; Xu, Z.; Tian, B.; et al. Pleiotropic effects of RecQ in Deinococcus radiodurans. Genomics 2009, 94, 333–340. [Google Scholar] [CrossRef] [PubMed]
  14. Devigne, A.; Ithurbide, S.; Bouthier de la Tour, C.; Passot, F.; Mathieu, M.; Sommer, S.; Servant, P. DdrO is an essential protein that regulates the radiation desiccation response and the apoptotic-like cell death in the radioresistant Deinococcus radiodurans bacterium. Mol. Microbiol. 2015, 96, 1069–1084. [Google Scholar] [CrossRef] [PubMed]
  15. Blanchard, L.; Guérin, P.; Roche, D.; Cruveiller, S.; Pignol, D.; Vallenet, D.; Armengaud, J.; de Groot, A. Conservation and diversity of the IrrE/DdrO-controlled radiation response in radiation-resistant Deinococcus bacteria. Microbiologyopen 2017, 6, e477. [Google Scholar] [CrossRef] [PubMed]
  16. Das, K.; Roychoudhury, A. Reactive oxygen species (ROS) and response of antioxidants as ROS-scavengers during environmental stress in plants. Front. Environ. Sci. 2014, 2, 53. [Google Scholar] [CrossRef]
  17. Arora, A.; Sairam, R.K.; Srivastava, G.C. Oxidative stress and antioxidative system in plants. Curr. Sci. 2002, 82, 1227–1238. [Google Scholar]
  18. Massé, E.; Vanderpool, C.K.; Gottesman, S. Effect of RyhB small RNA on global iron use in Escherichia coli. J. Bacteriol. 2005, 187, 6962–6971. [Google Scholar] [CrossRef]
  19. Zhan, Y.; Yan, Y.; Deng, Z.; Chen, M.; Lu, W.; Lu, C.; Shang, L.; Yang, Z.; Zhang, W.; Wang, W.; et al. The novel regulatory ncRNA, NfiS, optimizes nitrogen fixation via base pairing with the nitrogenase gene nifK mRNA in Pseudomonas stutzeri A1501. Proc. Natl. Acad. Sci. USA 2016, 113, E4348–E4356. [Google Scholar] [CrossRef]
  20. Zhang, H.; Zhan, Y.; Yan, Y.; Liu, Y.; Hu, G.; Wang, S.; Yang, H.; Qiu, X.; Liu, Y.; Li, J.; et al. The Pseudomonas stutzeri -specific regulatory ncRNA, NfiS, targets the katB mRNA encoding a catalase essential for optimal oxidative resistance and nitrogenase activity. J. Bacteriol. 2019, 201, e00334-19. [Google Scholar] [CrossRef]
  21. Romby, P.; Charpentier, E. An overview of RNAs with regulatory functions in gram-positive bacteria. Cell. Mol. Life Sci. 2010, 67, 217–237. [Google Scholar] [CrossRef] [PubMed]
  22. Citartan, M.; Raabe, C.A.; Hoe, C.H.; Rozhdestvensky, T.S.; Tang, T.H. Bacterial sRNAs: Regulation in Stress. In Stress and Environmental Regulation of Gene Expression and Adaptation in Bacteria; John Wiley & Sons Inc.: Hoboken, NJ, USA, 2016. [Google Scholar]
  23. Pichon, C.; Felden, B. Small RNA gene identification and mRNA target predictions in bacteria. Bioinformatics 2008, 24, 2807–2813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Tang, T.H.; Polacek, N.; Zywicki, M.; Huber, H.; Brugger, K.; Garrett, R.; Bachellerie, J.P.; Hüttenhofer, A. Identification of novel non-coding RNAs as potential antisense regulators in the archaeon Sulfolobus solfataricus. Mol. Microbiol. 2005, 55, 469–481. [Google Scholar] [CrossRef] [PubMed]
  25. Storz, G.; Vogel, J.; Wassarman, K.M. Regulation by Small RNAs in Bacteria: Expanding Frontiers. Mol. Cell 2011, 43, 880–891. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Tsai, C.H.; Liao, R.; Chou, B.; Contreras, L.M. Transcriptional analysis of Deinococcus radiodurans reveals novel small RNAs that are differentially expressed under ionizing radiation. Appl. Environ. Microbiol. 2015, 81, 1754–1764. [Google Scholar] [CrossRef] [PubMed]
  27. Meyer, M.; Kircher, M. Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb. Protoc. 2010, 6. [Google Scholar] [CrossRef]
  28. Hua, X.; Hua, Y. Improved complete genome sequence of the extremely radioresistant bacterium Deinococcus radiodurans R1 obtained using PacBio singlemolecule sequencing. Genome Announc. 2016, 4, e00886-16. [Google Scholar] [CrossRef]
  29. Chen, Z.; Li, L.; Shan, Z.; Huang, H.; Chen, H.; Ding, X.; Guo, J.; Liu, L. Transcriptome sequencing analysis of novel sRNAs of Kineococcus radiotolerans in response to ionizing radiation. Microbiol. Res. 2016, 192, 122–129. [Google Scholar] [CrossRef]
  30. Busch, A.; Richter, A.S.; Backofen, R. IntaRNA: Efficient prediction of bacterial sRNA targets incorporating target site accessibility and seed regions. Bioinformatics 2008, 24, 2849–2856. [Google Scholar] [CrossRef]
  31. Ashburner, M.; Ball, C.A.; Blake, J.A.; Botstein, D.; Butler, H.; Cherry, J.M.; Davis, A.P.; Dolinski, K.; Dwight, S.S.; Eppig, J.T.; et al. Gene ontology: Tool for the unification of biology. Nat. Genet. 2000, 25, 25–29. [Google Scholar] [CrossRef]
  32. Boyle, E.I.; Weng, S.; Gollub, J.; Jin, H.; Botstein, D.; Cherry, J.M.; Sherlock, G. GO::TermFinder—Open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics 2004, 20, 3710–3715. [Google Scholar] [CrossRef] [PubMed]
  33. McCarthy, F.M.; Wang, N.; Magee, G.B.; Nanduri, B.; Lawrence, M.L.; Camon, E.B.; Barrell, D.G.; Hill, D.P.; Dolan, M.E.; Williams, W.P.; et al. AgBase: A functional genomics resource for agriculture. BMC Genom. 2006, 7, 229. [Google Scholar] [CrossRef] [PubMed]
  34. Wu, J.; Mao, X.; Cai, T.; Luo, J.; Wei, L. KOBAS server: A web-based platform for automated annotation and pathway identification. Nucleic Acids Res. 2006, 34, W720–W724. [Google Scholar] [CrossRef] [PubMed]
  35. Ogata, H.; Goto, S.; Sato, K.; Fujibuchi, W.; Bono, H.; Kanehisa, M. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 1999, 27, 29–34. [Google Scholar] [CrossRef] [PubMed]
  36. Beckert, B.; Kedrov, A.; Sohmen, D.; Kempf, G.; Wild, K.; Sinning, I.; Stahlberg, H.; Wilson, D.N.; Beckmann, R. Translational arrest by a prokaryotic signal recognition particle is mediated by RNA interactions. Nat. Struct. Mol. Biol. 2015, 22, 767–773. [Google Scholar] [CrossRef] [PubMed]
  37. Jerabek-Willemsen, M.; André, T.; Wanner, R.; Roth, H.M.; Duhr, S.; Baaske, P.; Breitsprecher, D. MicroScale Thermophoresis: Interaction analysis and beyond. J. Mol. Struct. 2014, 1077, 101–113. [Google Scholar] [CrossRef] [Green Version]
  38. Buddeweg, A.; Sharma, K.; Urlaub, H.; Schmitz, R.A. sRNA41 affects ribosome binding sites within polycistronic mRNAs in Methanosarcina mazei Gö1. Mol. Microbiol. 2018, 107, 595–609. [Google Scholar] [CrossRef] [PubMed]
  39. Kullapanich, C.; Dubbs, J.M.; Mongkolsuk, S. Inactivation of the Agrobacterium tumefaciens ActSR system affects resistance to multiple stresses with increased H2O2 sensitivity due to reduced expression of hemH. Microbiology 2019. [Google Scholar] [CrossRef]
  40. Navasa, N.; Ferrero, M.Á.; Rodríguez-Aparicio, L.B.; Monteagudo-Mera, A.; Gutiérrez, S.; Martínez-Blanco, H. The role of RcsA in the adaptation and survival of Escherichia coli K92. FEMS Microbiol. Lett. 2019, 366, 1–6. [Google Scholar] [CrossRef]
  41. Czégény, G.; Kőrösi, L.; Strid, Å.; Hideg, É. Multiple roles for Vitamin B6 in plant acclimation to UV-B. Sci. Rep. 2019, 9, 1259–1266. [Google Scholar] [CrossRef]
  42. Coutinho, L.G.; de Oliveira, A.H.S.; Witwer, M.; Leib, S.L.; Agnez-Lima, L.F. DNA repair protein APE1 is involved in host response during pneumococcal meningitis and its expression can be modulated by vitamin B6. J. Neuroinflamm. 2017, 14, 243. [Google Scholar] [CrossRef] [PubMed]
  43. Staerck, C.; Gastebois, A.; Vandeputte, P.; Calenda, A.; Larcher, G.; Gillmann, L.; Papon, N.; Bouchara, J.P.; Fleury, M.J.J. Microbial antioxidant defense enzymes. Microb. Pathog. 2017, 110, 56–65. [Google Scholar] [CrossRef] [PubMed]
  44. Caro, A.A.; Cederbaum, A.I. Antioxidant properties of S-adenosyl-L-methionine in Fe2+-initiated oxidations. Free Radic. Biol. Med. 2004, 36, 1303–1316. [Google Scholar] [CrossRef] [PubMed]
  45. He, L.; Ma, X.; Li, Z.; Jiao, Z.; Li, Y.; Ow, D.W. Maize OXIDATIVE STRESS2 homologs enhance cadmium tolerance in Arabidopsis through activation of a putative SAM-dependent methyltransferase gene. Plant Physiol. 2016, 171, 1675–1685. [Google Scholar] [CrossRef]
Figure 1. Stress tolerance analysis of OsiA. (A) Northern blotting identification of OsiA. CDR is the control conditions, while TDR is the H2O2-treated conditions, with the arrow indicating OsiA, while 5S rRNA bands are shown as loading controls. (B) The expression profile of OsiA. ** means that the expression of OsiA is significantly different from that of the wild-type strain, that is, the p value is less than 0.01. (C) The phenotypic analysis of osiA knockout mutant, the osiA complementary strain, and the control under different H2O2-treated conditions. (D) Survival growth analysis of the osiA knockout mutant, osiA complementary strain, and the control under different H2O2-treated conditions.
Figure 1. Stress tolerance analysis of OsiA. (A) Northern blotting identification of OsiA. CDR is the control conditions, while TDR is the H2O2-treated conditions, with the arrow indicating OsiA, while 5S rRNA bands are shown as loading controls. (B) The expression profile of OsiA. ** means that the expression of OsiA is significantly different from that of the wild-type strain, that is, the p value is less than 0.01. (C) The phenotypic analysis of osiA knockout mutant, the osiA complementary strain, and the control under different H2O2-treated conditions. (D) Survival growth analysis of the osiA knockout mutant, osiA complementary strain, and the control under different H2O2-treated conditions.
Microorganisms 07 00422 g001
Figure 2. Catalase activity and total antioxidative activity of the wild type strain and osiA mutant strain. (A) Catalase activity of the wt, ∆osiA, and osiA complementary strain. (B) Total antioxidative activity of the wt, ∆osiA, and osiA complementary strain. ** means the difference is very significant compared with the control group, that is, the p value is less than 0.01.
Figure 2. Catalase activity and total antioxidative activity of the wild type strain and osiA mutant strain. (A) Catalase activity of the wt, ∆osiA, and osiA complementary strain. (B) Total antioxidative activity of the wt, ∆osiA, and osiA complementary strain. ** means the difference is very significant compared with the control group, that is, the p value is less than 0.01.
Microorganisms 07 00422 g002
Figure 3. Microscale thermophoresis of selected sRNAs and their target genes. Red bases mean complementary bases in the binding sequence, the red curve is the fitted combination curve, and the Kd (dissociation equilibrium constant) value is the binding constant of sRNAs and their targets. (A,C,E) are the binding of dr_1998(katA) mRNA, dr_A0202(sodC) mRNA, and dr_A0259(katE) mRNA, respectively, with OsiA-wt. (B,D,F) are the binding of dr_1998 mRNA, dr_A0202 mRNA, dr_A0259 mRNA, respectively, with OsiA-mut.
Figure 3. Microscale thermophoresis of selected sRNAs and their target genes. Red bases mean complementary bases in the binding sequence, the red curve is the fitted combination curve, and the Kd (dissociation equilibrium constant) value is the binding constant of sRNAs and their targets. (A,C,E) are the binding of dr_1998(katA) mRNA, dr_A0202(sodC) mRNA, and dr_A0259(katE) mRNA, respectively, with OsiA-wt. (B,D,F) are the binding of dr_1998 mRNA, dr_A0202 mRNA, dr_A0259 mRNA, respectively, with OsiA-mut.
Microorganisms 07 00422 g003
Figure 4. Halflife of katA mRNA in the wildtype strain and osiA mutant strain.
Figure 4. Halflife of katA mRNA in the wildtype strain and osiA mutant strain.
Microorganisms 07 00422 g004
Table 1. Molecular function of sRNA target genes.
Table 1. Molecular function of sRNA target genes.
FunctionsRNAs (+)sRNAs (−)
Kinases1517
Transcription factor714
Enzymes4255
Structural protein 913
Table 2. Functional annotation of oxidative response sRNAs in the KEGG pathway.
Table 2. Functional annotation of oxidative response sRNAs in the KEGG pathway.
sRNA NamePathway ManeTarget Genes Products
IGR_3053PeroxisomeKatA, SodC
Mismatch repairRplV
Two component systemmethyl-accepting chemotaxis protein
IGR_2205Peroxisomeclass V aminotransferase
Mismatch repairDNA ligase, DNA polymerase III, tau/gamma subunit, DR_2586
Glutathione metabolismglucose-6-phosphate 1-dehydrogenase
Carotenoid metabolismphytoene dehydrogenase
Two component systemNatA, DR_A0009, methyl-accepting chemotaxis-like protein
Vitamin B6 metabolismpyridoxamine 5-phosphate oxidase
IGR_1662PeroxisomeKatE,
Mismatch repairMutT, Ssb, MutS
Carotenoid metabolismphytoene synthase
Two component system5-(carboxyamino)imidazole ribonucleotide mutase
IGR-1449Mismatch repairMutL, DNA helicase II
Glutathione metabolismpotassium-transporting ATPase subunit C
Two component systempotassium-transporting ATPase subunit C, KdpD-related protein
Vitamin B6 metabolismRpsO, PdxT
IGR_1916Peroxisomelong-chain fatty acid--CoA ligase
Mismatch repairDR_1244
Glutathione metabolismValS
IGR_2408Mismatch repair3-oxoacyl-acyl carrier protein reductase
Two component systemmethyl-accepting chemotaxis protein, CheA-related protein
Vitamin B6 metabolismthreonine synthase
IGR_1612Mismatch repairDNA polymerase III subunit epsilon, RuvA
IGR_2590Peroxisomeacyl-CoA synthetase
Mismatch repair3-isopropylmalate dehydratase large subunit
Vitamin B6 metabolismpyridoxamine kinase
IGR_884Mismatch repairXseA
Glutathione metabolismdiaminopimelate decarboxylase
Two component systempotassium-transporting ATPase subunit B
IGR_771Glutathione metabolismnitrogen regulatory protein P-II, leucyl aminopeptidase
Carotenoid biosynthesislycopene cyclase
IGR_2060Two component systempotassium-transporting ATPase subunit A
Vitamin B6 metabolismRuvB, PdxS
IGR_1174Peroxisomeisocitrate dehydrogenase
Glutathione metabolismisocitrate dehydrogenase
IGR_2150Mismatch repairshort chain dehydrogenase/reductase family oxidoreductase
Glutathione metabolismcephalosporin acylase
IGR_2479Mismatch repairDNA polymerase III subunit alpha
Two-component systemsuccinate dehydrogenase, cytochrome subunit
IGR_2389PeroxisomeCu/Zn family superoxide dismutase
IGR_76Mismatch repairsingle-stranded DNA-binding protein
IGR_585Mismatch repairRpsS
IGR_2012Mismatch repairsingle-stranded-DNA-specific exonuclease
IGR_2690Glutathione metabolismArginase, 6-phosphogluconate dehydrogenase-like protein
IGR_373Two-component systemNADH dehydrogenase II
IGR_826Two-component systemSdhB
IGR_1951Two-component systemDNA-binding response regulator
IGR_388Two-component systemPilH

Share and Cite

MDPI and ACS Style

Chen, Y.; Xue, D.; Sun, W.; Han, J.; Li, J.; Gao, R.; Zhou, Z.; Zhang, W.; Chen, M.; Lin, M.; et al. sRNA OsiA Stabilizes Catalase mRNA during Oxidative Stress Response of Deincoccus radiodurans R1. Microorganisms 2019, 7, 422. https://doi.org/10.3390/microorganisms7100422

AMA Style

Chen Y, Xue D, Sun W, Han J, Li J, Gao R, Zhou Z, Zhang W, Chen M, Lin M, et al. sRNA OsiA Stabilizes Catalase mRNA during Oxidative Stress Response of Deincoccus radiodurans R1. Microorganisms. 2019; 7(10):422. https://doi.org/10.3390/microorganisms7100422

Chicago/Turabian Style

Chen, Yun, Dong Xue, Wenjie Sun, Jiahui Han, Jiang Li, Ruyu Gao, Zhengfu Zhou, Wei Zhang, Ming Chen, Min Lin, and et al. 2019. "sRNA OsiA Stabilizes Catalase mRNA during Oxidative Stress Response of Deincoccus radiodurans R1" Microorganisms 7, no. 10: 422. https://doi.org/10.3390/microorganisms7100422

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