Next Article in Journal
Immunodominant B-Cell Linear Epitope on the VP1 P Domain of a Feline Norovirus Cat Model
Previous Article in Journal
It Is Time to Strengthen the Malaria Control Policy of the Democratic Republic of Congo and Include Schools and School-Age Children in Malaria Control Measures
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interrelation between Stress Management and Secretion Systems of Ralstonia solanacearum: An In Silico Assessment

1
Department of Food Science and Human Nutrition, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
2
Department of Medical Zoology, Kyung Hee University School of Medicine, Seoul 02447, Korea
3
Medical Research Center for Bioreaction to Reactive Oxygen Species and Biomedical Science Institute, School of Medicine, Graduate School, Kyung Hee University, Seoul 02447, Korea
4
Department of Biological Sciences, Texas Tech University, Lubbock, TX 79409, USA
5
Department of Biology, Clarkson University, Potsdam, NY 13699, USA
6
Department of Botany, M.U.C. Women’s College, University of Burdwan, Bardhaman 713104, West Bengal, India
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Pathogens 2022, 11(7), 730; https://doi.org/10.3390/pathogens11070730
Submission received: 19 April 2022 / Revised: 18 June 2022 / Accepted: 24 June 2022 / Published: 27 June 2022
(This article belongs to the Section Bacterial Pathogens)

Abstract

:
Ralstonia solanacearum (Rs), the causative agent of devastating wilt disease in several major and minor economic crops, is considered one of the most destructive bacterial plant pathogens. However, the mechanism(s) by which Rs counteracts host-associated environmental stress is still not clearly elucidated. To investigate possible stress management mechanisms, orthologs of stress-responsive genes in the Rs genome were searched using a reference set of known genes. The genome BLAST approach was used to find the distributions of these orthologs within different Rs strains. BLAST results were first confirmed from the KEGG Genome database and then reconfirmed at the protein level from the UniProt database. The distribution pattern of these stress-responsive factors was explored through multivariate analysis and STRING analysis. STRING analysis of stress-responsive genes in connection with different secretion systems of Rs was also performed. Initially, a total of 28 stress-responsive genes of Rs were confirmed in this study. STRING analysis revealed an additional 7 stress-responsive factors of Rs, leading to the discovery of a total of 35 stress-responsive genes. The segregation pattern of these 35 genes across 110 Rs genomes was found to be almost homogeneous. Increasing interactions of Rs stress factors were observed in six distinct clusters, suggesting six different types of stress responses: membrane stress response (MSR), osmotic stress response (OSR), oxidative stress response (OxSR), nitrosative stress response (NxSR), and DNA damage stress response (DdSR). Moreover, a strong network of these stress responses was observed with type 3 secretion system (T3SS), general secretory proteins (GSPs), and different types of pili (T4P, Tad, and Tat). To the best of our knowledge, this is the first report on overall stress response management by Rs and the potential connection with secretion systems.

1. Introduction

Ralstonia solanacearum (Rs) is a group of aerobic, non-spore-forming, Gram-negative plant pathogens that belong to the Burkholderiaceae family of the β-proteo-bacteria group. Rs is a “bacterial species complex” consisting of four monophyletic phylotypes [1]. This pathogen is designated as the world’s deadliest bacterial plant pathogen because of its wide-range host specificity (>200 plant species), broad geographical distribution, abundance and persistence in soil, root invasion, tissue-specific tropism, multiplication, and colonization (>109 C.F.U. per g fresh weight) in xylem vessels [2,3]. Rs is also considered “quarantine organisms”, “bioterrorism”, and “double usage agents” by different regulatory authorities in the USA and Europe [4]. On the other hand, due to wide host diversity, Rs may be an attractive model to investigate the question of adaptation to the host environment.
The stress response in bacteria involves processes that enable bacteria to survive in adverse environmental conditions. The genetic makeup helps bacteria to sense the changing environment and to react accordingly by modulating the expression of certain gene cascades and protein activities. In particular, stress response management is critical for pathogens to survive inside the host and to defeat the host’s primary immunity. The stress response in bacteria is a multi-network phenomenon, and different sets of genes responsible for different types of stresses (such as acidic, alkaline, desiccation, osmotic, nutrient deficiency, heat, cold, etc.) are expressed accordingly with the help of effector molecules to maintain cellular integrity [5]. According to Hengge-Aronis (2002), the bacterial stress response can be defined as the sequential alteration of the expression of several genes and their protein levels to cope with extreme conditions [6]. Like all other pathogenic bacteria, Rs also must face stresses during the first physical contact with the host. It is already well known that plant roots secrete reactive oxygen species (ROS) for their own protection in the rhizosphere environment, which may impart host-associated stress to Rs during the first interaction with the plant [7,8]. Inside the host, different types of signals (e.g., pH and temperature) and nutrient deficiency accelerate the adaptive stress responses in Rs, along with the expression of virulence genes. For example, the Dps gene was reported to contribute to oxidative stress tolerance during tomato plant infection by Rs [7]. Flores-Cruz and Allen (2011) reported several ROS-scavenging enzymes and the OxyR gene in the genome of Rs, which play a critical role during stress conditions or exposure to extreme environments [9]. Fang et al. (2016) classified bacterial stress responses during their infection of hosts; however, the detailed mechanism of stress management in Rs has not been fully explored to date. Notably, the correlation between stress responses and the secretion systems of Rs needs a deeper understanding.
The identification and distribution of Rs proteins responsible for regulating resistance to different types of host-environment stress conditions could be important factors in pathogenicity. In the present investigation, an in silico approach was taken to identify and monitor the distribution of stress-responsive genes and their relationship with bacterial secretion systems of Rs.

2. Results

2.1. Mining and Categorization of Rs Stress-Responsive Genes

Based on experimental records, the initial number of query candidates for bacterial stress-responsive genes (SRGs) was 47 (Table 1).
However, based on orthologs present in Rs, the numbers were reduced to 28 (Table S1). The protein sequences of these 28 SRGs were then used as a secondary query for pBLAST across available annotated Rs genomes (Table S2). A protein–protein interaction network (PIN) was studied within Rs genomes using STRING V11, taking these 28 SRGs as an input file. In string analysis, increasing interactions (no. of edges expected: 39; found: 88) resulted in a total no. of 35 nodes (representing genes, average node degree: 5.03), which indicated 35 SRGs within Rs (Table 2). Thus, in addition to the initially identified 28 SRGs, 7 new SRGs (total of 35 SRGs) were discovered from STRING analysis (Table S3).
The cluster analysis (k-means clustering) of these 35 SRGs revealed five distinct clusters indicating five different types of stress responses: membrane stress response (MSR), osmotic stress response (OSR), oxidative stress response (OxSR), nitrosative stress response (NxSR), and DNA damage stress response (DdSR) (Figure 1).
However, a PIN was not observed in the case of Histone-like protein H-NS (ortholog of K03746, BBJ49629.1), associated with temperature-dependent stress regulation. Interestingly, DdSR was found in the center among these five clusters, showing independent PIN connections with the other four clusters.

2.1.1. MSR

The present analysis revealed that eight SRGs (including nlpD ortholog CBJ51536.2) contributed to the MSR cluster (Figure 1). The cluster is composed of outer membrane metallopeptidase (NlpD), periplasmic serine endoprotease (DegQ), RNA polymerase σ (RpoE), transcription factor (LexA) repressor, and binary systems (RseA and RseB; RecA and RecX). In Figure 1, the PIN between alternative σE (RpoE and RseA), degQ, and two-component systems Rcs (RcsB and RcsC) is distinct. The MSR cluster is linked to the other four clusters via the DdSR cluster. In particular, a PIN was found among rpoE, degQ, nlpD, recA, and lexA of MSR with rpoA, rpoB, rpoC, and rpoS of DdSR.

2.1.2. OSR

Six SRG candidates (including ompR orthologs CBJ35556.1 and CBJ51104.1) were found to construct the cluster OSR (Figure 1). In the present investigation, a strong PIN was recorded among the osmotic stress regulator ompR and envZ of Rs. Along with ompR, a strong PIN was also connected to the two-component system Bae (BaeR and BaeS). The OSR cluster relates to the DdSR cluster via ompR and envZ.

2.1.3. OxSR

The present analysis revealed that seven SRGs (including Dps ortholog CBJ50185.1) contributed to the OxSR cluster (Figure 1). Three orthologs of the enzyme catalase (katA, katB, and katG) and two orthologs of the enzyme superoxide dismutase (sodB and sodC) were identified within Rs genomes (Table 2). Besides multiple and redundant ROS-scavenging enzymes, PINs within oxidative stress regulators such as Dps and oxyR were also identified (Figure 1). In Figure 1, it appears that KatA, KatB, KatG, SodB, and SodC work together in an interactive fashion to provide protection against oxidative stress. The carbon and nitrogen starvation regulator Dps (Figure 1) interacts with stress-responsive enzymes (katE and sodB) and another stress regulator (oxyR). String analysis also suggests that oxyR also interacts with other stress-responsive enzymes (katE and sodB) and stress regulators (Dps) (Figure 1). The PIN between the OxSR cluster and the DdSR cluster is mediated via sodB, kata, katB, and oxyR.

2.1.4. NxSR

Three SRG candidates (nsrR, norB, and aniA) were found to construct the cluster NxSR (Figure 1). Two probable candidates (polA and fur ortholog CBJ50125.1) were not included in this cluster due to the lack of evidence (known interactions) of the involvement of these two SRGs in the nitrosative stress response. A single ortholog (norB) was identified for the enzyme nitric-oxide reductase (Table 2). The PIN obtained from the STRING server indicates interactions among these three SRG candidates (Figure 1). However, a single PIN was observed between the clusters NxSR and DdSR, mediated through norB.

2.1.5. DdSR

The present investigation discovered eight to ten (including polA and fur ortholog CBJ50125.1) SRG candidates that constitute the DdSR cluster (Figure 1). This cluster is composed of DNA polymerase I (polA), different RNA polymerase subunits (rpoA for α, rpoC for β, rpoS for σ, and rpoZ for ω), transcription factors (dskA and fnr), and kinase (gmk). RpoC is a known amino acid starvation regulator, and Fur is a known metal stress regulator, whereas FNR is considered an oxygen sensor. The DdSR cluster is the central cluster of all SRG candidates, and all other clusters are independently linked to it via PIN.

2.2. Homogeneous Distribution of SRGs across Rs Genomes

Among 35 identified SRGs (Table 2), 27 were found to be present in 110 fully annotated Rs genomes taken from the NCBI database. The heatmap analysis with the BLAST results of these 35 SRGs also indicates the same (Figure 2).
We also constructed a cluster dendrogram of the 110 experimental Rs genomes based on the presence and absence of SRGs identified in the present study. We observed most of the strains (92) to be on a straight line (Figure S1), indicating no significant differences in the pattern of stress responses among the different strains of Rs.

2.3. SRGs and Secretion Systems of Rs

For secretion system components, we used the KO numbers published in our previous study with Bradyrhizobium spp. [48] as reference data, and the results are provided in Table S4. String analysis of stress-responsive genes in relation to different secretion systems (type secretion systems including T1SS, T2SS, T3SS, T4SS, T5SS, and T6SS; general secretory protein system or GSPS; and pili systems including T4P, Tad, and Tat) of Rs was performed. A strong PIN with SRGs was found for three different secretion systems, viz. T3SS, GSPS, and pili systems (T4P, Tad, and Tat).

2.4. SRGs and T3SS

Initially, all 35 SRGs (found in this study, Table 2) and 14 T3SS candidate genes (from our previous study) were used as the input file for the String analysis. However, for a better understanding of the PIN between SRGs and T3SS, the input file was manually curated. The final output file (no. of edges expected: 39; found: 99) resulted in a total no. of 39 nodes (representing genes, average node degree: 5.08). Among these 39 nodes, 8 represent T3SS candidate genes (hrcC, hrcJ, hrcN, fliF, fliG, fliH, flhA, and flgJ) and are present as a single cluster (Figure 3).
The candidates for the T3SS cluster are provided in Table S4. In Figure 4, out of 35 SRGs (except for h-ns, fur, ompR, and dnr), 31 were found to be distributed into the five predetermined clusters (Figure 3). However, a PIN was found to be present between T3SS (via flgJ and fliG) and MSR (via rpoE and nlpD) and DdSR (via polA) clusters of SRGs (Figure 3).

2.5. SRGs and GSPS

At first, all 35 SRGs (found in this study, Table 2) and 9 GSPS candidate genes (from our previous study) were used as an input file for the String analysis. The final output file (no. of edges expected: 54; found: 146) resulted in a total no. of 40 nodes (representing genes, average node degree: 7.3). Among these 40 nodes, 9 represent GSPS candidate genes (secA, secB, secD, secE, sceF, secG, secY, ftsY, and ffh) and are present as a single cluster (Figure 4). Information on the candidate genes in the GSPS cluster is provided in Table S4.
Among 35 SRGs (except for h-ns, fur, ompR, and dnr), 31 were distributed into the five predetermined clusters (Figure 4). There was no direct PIN between GSPS and DdSR, GSPS and NxSR, and GSPS and OSR. In contrast, a PIN was recorded between the GSPS cluster (most of the candidate genes) and the MSR cluster (via nlpD and recA). A PIN was also found to exist between the GSPS cluster (via secA and secV) and the OxSR cluster (via sodB and oxyR). A strong PIN was also observed between the GSPS cluster (most of the candidate genes) and almost all SRGs of the DdSR cluster (except fnr).

2.6. SRGs and Pili Systems

Initially, all 35 SRGs (found in Table 2) and 26 pili system candidate genes (from our previous study) were used as input files for the String analysis. The output file (no. of edges expected: 38; found: 104) contained a total no. of 43 nodes representing genes with an average node degree of 4.84 (Figure 5).
Among these 43 nodes, 12 are represented by different pili candidate genes and are present as three separate clusters representing three different pili systems: type 4 pili or T4P (pilA, pilB, pilN, pilT, pilU, and gspE), tat pili (tatB and tatC), and tad pili (tadB, tadC, tadZ, and cpaF) (Figure 5). The SRGs were distributed into five clusters as before (Figure 5). A PIN was observed between T4P (via pilT) and DdSR (via gmk) and MSR (via recA), between Tat (via tatC and tatD) and MSR (via nlpD), DdSR (via recA), and NxSR (via aniA), and between Tad (via tadZ) and OSR (via envZ).

3. Discussion

3.1. MSR SRGs of Rs

In the present study, we observed that protein–protein interactions between σE and two-component systems such as Bae and Rcs were distinct in Rs. We report the interaction with DegQ (Table 2), a type 1 protease responsible for cleaving RseA and releasing σE/rpoE. The most crucial component of membrane stress, σE is extensively studied in human pathogens and is known to bind with the cytoplasmic membrane with the help of RseA [49]. The RecA- and LexA-guided SOS response is a well-studied phenomenon in several human enteric pathogens [19]. LexA is also reported to control the virulence gene expression in pathogenic bacteria [20]. Membrane permeability is affected by lexA and recA mutations in Escherichia coli K12 [50]. Genotoxic agents (e.g., mitomycin C) can significantly increase membrane vesicle production by promoting RecA-dependent stress responses [51]. Graphene oxide was reported to induce MSR in Rs through a similar mechanism [52].

3.2. OSR SRGs of Rs

Membrane stress is generally caused by a change in osmotic pressure, low pH, exposure to antimicrobial cationic peptides, and misfolding of envelope proteins such as secretin. Rs may survive in pure water and grow at a water activity (aW) near 1 [53]. Evidence suggests that regulating cytoplasmic composition and hydration is a key objective of cellular homeostasis in such conditions [54]. OmpR is a crucial factor for overcoming osmotic stress while entering the host and controlling the Bae two-component systems (Table 2). MacIntyre et al. showed that the disaccharide trehalose contributes to OSR of R. solanacearum strain GMI1000 [55]. Rs mutants that were unable to synthesize trehalose suffered significant growth inhibition in the presence of ionic (NaCl) or nonionic (high–molecular weight polyethylene glycol) osmotic stress [55].

3.3. OxSR SRGs of Rs

Disruption of the membrane releases components of the electron transport chain named quinones and flavoproteins, which generate reactive oxygen species (ROS) from molecular oxygen [49]. In the present investigation, along with molecular oxygen scavenging enzymes such as catalase and superoxide dismutase, OxyR was also reported (Table 2). Flores-Cruz and Allen previously reported the crucial role of OxyR as an oxidative stress response regulator in Rs [9]. The cysteine residues in the OxyR protein undergo redox conformational changes and can sense different components of stress factors [10]. Similar to E. coli, the expression of Dps in Rs was reported to be increased during the stationary phase and could be mediated through OxSR [7]. Vinegar residue substrate was reported to inhibit oxidative damage in R. solanacearum HB511 via the regulation of excess ROS production [56]. Recently, Yang et al. reported that OxSR in R. solanacearum CQPS-1 induced by daphnetin is mediated by SRGs, such as katGb, coxM, Rsp0993, Rsc2493, and dnaK [57].

3.4. NxSR SRGs of Rs

Disruption of the membrane releases components of the electron transport chain (quinones and flavoproteins). These components are known to generate reactive nitrogen species (RNS), such as peroxynitrite (ONOO) and nitric oxide (NO). RNS are very toxic to any kind of cell [58]. In this study, both ROS-scavenging enzyme nitric-oxide reductase (coded by NorB) and nitrosative stress-responsive factors (NsrR and AniA) were reported in the Rs genomes under investigation (Table 2). A recent study elucidated the nitrosative stress response in R. solanacearum GMI1000 and reported that NorA, HmpX, and NorB play a critical role in reducing nitrosative stress during plant pathogenesis [59].

3.5. DdSR SRGs of Rs

ROS, RNS, and low pH are known causal factors of DNA damage [60,61]. The σS factor and DksA are responsible for overcoming amino acid limitations by activating several amino acid biosynthesis genes [18,62]. Brown and Allen reported the involvement of three SRGs, uvrA1 (ipx64), dinF (ipx62), and dnaN (ipx61), in the SOS response of R. solanacearum K60 against plant-derived DNA-damaging substances [63].

3.6. Global Strategy in Stress Response of Rs

The versatility of the pathogenicity of Rs depends on the genetic diversity [64]. The diversity of Rs is broad, and thus, the genetic variation among different phylotypes is distinct. The cluster analysis of SRGs demonstrated that there are no significant differences among the Rs species complex (Figure S1). To explain these mutually incompatible results, it is important to understand the mosaic structure of Rs genomes. Due to remarkable genetic diversity, the genomes of Rs possess a number of similar regions scattered within variable regions [64]. Usually, genes encoding stress-responsive enzymes and factors are subject to strong levels of purifying selection, similar to most housekeeping genes [65]. The present findings suggest that the stress-responsive genes in Rs are under pressure of purifying selection. Our present finding is also supported by the experimental evidence provided by Puigvert et al. [66]. Transcriptomic analysis of R. solanacearum UY031 in their research found only two differentially expressed genes between resistant and susceptible plant accessions, indicating that the bacterial component plays a minor role in the establishment of disease [66].
Rs enters its plant host through root lesions, followed by multiplication and colonization in the root cortex and xylem, respectively. Finally, wilt symptoms develop, and then the plants die, releasing the pathogen back into the soil. From the five consistent clusters of SRGs in Rs, the stress responses can be classified into five categories (Figure 6a). Studying biological networks enables a deeper investigation of biological systems [67], and PINs have been constructed in many bacteria [68,69]. In the present study, our analysis indicates that SRGs work together interactively to fight against diverse stresses (Figure 6b). Currently, very few reports explore the relationship between clusters of SRGs and bacterial secretion systems. SRGs of the DdSR cluster were found to play a central role among all five clusters (Figure 6a). SRGs of all other clusters are linked to DdSR via PINs; however, PINs were completely absent between these other clusters. This finding strongly suggests that all other stress responses (i.e., MSR, OSR, OxSR, and NxSR) act through DdSR (Figure 6a). It also implies that the expression of SRGs of a single cluster cannot be treated as a signature of progressive stress responses.

3.7. Interactions among T3SS and SRGs of Rs

In the present investigation, PINs were observed between the T3SS cluster (representing eight genes: hrcC, hrcJ, hrcN, fliF, fliG, fliH, flhA, and flgJ) and the MSR cluster of SRGs (Figure 6b). Flores-Kim et al. previously reported the interaction between T3SS and MSR in Yersinia [66], which aligns with our findings. The secretion systems of Rs form a molecular arsenal that facilitates niche adaptation, host invasion, and evasion of plant defenses. As for many bacterial pathogens, the main virulence determinant in Rs is T3SS [70], which injects a number of effector proteins into plant cells, causing disease in hosts or a hypersensitive response in resistant plants. In this work, the PIN between the MSR cluster of SRGs and T3SS indicates the critical role of SRGs in infection or disease establishment (Figure 3).

3.8. Interactions among GSPS and SRGs of Rs

Similarly, PINs were also recorded among the GSPS cluster (representing nine genes: secA, secB, secD, secE, sceF, secG, secY, ftsY, and ffh), MSR cluster, and OxSR cluster. Interestingly, interactions between Sec and MSR in Escherichia coli and the Sec secretion system and OxSR in Pseudomonas aeruginosa were also reported [71,72]. As in most bacteria, the Sec pathway is the main route by which secretory proteins are exported across the cytoplasmic membrane in Rs [73]. Upon their exit from the membrane translocase complex, the proteins must acquire a properly folded conformation to gain full activity. Heat stress and overproduction of specific secretory proteins can result in the accumulation of misfolded proteins outside the cytoplasm [74]. The PIN between the DdSR cluster of SRGs and GSPS, along with the PIN between the MSR cluster of SRGs and GSPS, clearly indicates the importance of SRGs in the regulation of GSPS (Figure 4). The accumulation of misfolded proteins produces cellular stress responses leading to refolding or degradation of abnormally folded, non-functional proteins [75]. Rs may sense the secretion stress via the extracytoplasmic sigma factor σE.

3.9. Interactions among Different Pili Systems and SRGs of Rs

Furthermore, we explored the interactions among the T4P cluster, MSR, and DdSR, which was reported by a previous study with E. coli [76]. A strong PIN was observed between the T4P cluster, the MSR cluster, and NxSR. Rs binds to host cells by a multitude of pili. Therefore, the adhesion strength of a multipili-binding bacterium largely depends on the cooperativity of the attaching pili [73]. Interestingly, different types of pili of Rs were found to be under the control of different SRG clusters. For example, T4P pili are under the control of DdSR and MSR clusters, and Tat pili are controlled by MSR and NxSR clusters, while Tad pili are under the control of the OSR cluster (Figure 6b).

3.10. Network of SRGs and Bacterial Secretion Systems

In this work, we show the interaction of various SRG clusters with specific secretion components. For example, T3SS, GSPs, and T4P are related to MSR and DdSR, while Tad and Tad interact with OSR and NxSR, respectively (Figure 6b). A detailed transcriptome analysis of Rs provides the experimental proof that genes encoding T3SS (hrpY, hrpX, hrpK, and hrcT) and T4P (pilG, pilH, pilN, pilM, pilY, pilW, and fimV) are upregulated during host infection, along with peroxidases, catalases (katE and katG), and alkyl hydroperoxide reductases (ahpC1 and ahpF), suggesting adaptations to combat stress during plant infection [66]. Recently, Hendrich et al. also provided experimental proof that NxSR regulates T3SS in Rs [77]. Taken together, these findings indicate a close connection between bacterial stress management and secretion systems.

4. Materials and Methods

4.1. Data Mining

Completely assembled and annotated genomes (as of 31 March 2022) of Rs were downloaded either from the NCBI genome database (https://www.ncbi.nlm.nih.gov/genome) (accessed on 31 March 2022) or from the KEGG genome database (https://www.genome.jp/kegg/genome) (accessed on 31 March 2022), which includes 110 genomes (Table S4). Sequences of stress-responsive enzymes and stress regulators were retrieved from the UniPort database.

4.2. Mining of Bacterial Genes Coding for Stress-Responsive Enzymes and Stress Regulators

The present study examined a protocol for capturing the maximum number of candidate genes associated with bacterial stress responses (Figure 7).
To develop this protocol, information regarding bacterial stress-responsive enzymes and stress regulators was collected from published articles, followed by the preparation of a dataset. This dataset was used as a primary reference query in subsequent steps (Table 1).
Initially, the presence of these candidate genes within Rs genomes was searched using the KEGG platform, and orthologs were identified from the KO database of the KEGG platform (Table 2).
Nucleotide sequences of KO orthologs of the SRGs are presented in Table S1. The secondary query was carried out on Rs genomes using KO ortholog protein sequences. To examine the diversity of these genes across Rs genomes, microbial genome BLASTs were performed using protein sequences and the megaBLAST algorithm (Tables S1 and S2). Furthermore, the results obtained from the pBLAST search were again verified with the help of the KO orthology database in KEGG among the available Rs genomes.

4.3. Mining of Bacterial Secretion System Genes

For secretion system components, we used the KO numbers published in our previous study with Bradyrhizobium spp. [48] as reference data, and the results are provided in Table S4.

4.4. Protein–Protein Interaction Network (PIN) Study

PINs are an important clue for understanding the cellular physiological process, and STRING is one of the best platforms to construct and visualize such interaction networks [78]. The STRING v11.0 (http://string-db.org) (accessed on 22 May 2022) server is used for predicting structure–function aspects of genes and proteins [79]. Based on the approximate probability that a predicted link exists between two proteins in the same metabolic map in the KEGG database, different ranges of confidence scores of protein–protein interactions were assigned. The SRG sequences of Rs were aligned with Rs PS107 using BLASTX. The first aligned protein with an E-value below 1 × E−10 was considered a homologous protein. Then, these homologous proteins and their corresponding interactions were extracted from the whole interaction dataset of the related organism to compose the model organism-based protein–protein interaction sub-network. In order to obtain high-quality PINs, we considered the interactions with the highest confidence limits (0.9) from String. Here, we used this platform to explore the relation between different SRGs and generated an interrelation network. We also used this platform to study PINs between SRGs and bacterial secretion systems.

4.5. Statistical Analysis

All of the genes coding for stress-responsive enzymes and stress regulators were selected based on the BLAST score (Table S3). To access the distribution of SRGs within different phylotypes of Rs, the data obtained from BLAST results were arranged in a raw data matrix, and a heatmap was generated using Morpheus (https://software.broadinstitute.org/morpheus/) (accessed on 25 May 2022). The cluster analysis was performed with the help of Past 3.17 using binary matrix data (based on the presence or absence of an ortholog) [80]. The neighbor-joining method was used to construct a cluster dendrogram using 1000 as the bootstrap value.

5. Conclusions

This study provides a detailed report of the SRGs of Rs, their distribution, classification, and interaction with bacterial secretion systems. This is the first report of SRGs across Rs genomes. From this investigation, it is proposed that the SRGs of DdSR play a pivotal role among all other clusters. Multivariate analysis with SRGs provides evidence in support of a purifying selection of the stress management system in Rs. Interactions between SRGs and bacterial secretion systems in Rs were also elucidated for the first time.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/pathogens11070730/s1. Figure S1: Cluster dendrogram based on the distribution of stress-responsive genes within R. solanacerium species complex; Table S1: Secondary dataset for stress-responsive enzymes and factors of R. solanacearum using KO orthologs; Table S2: Secondary dataset for stress-responsive enzymes and factors of R. solanacearum using KO orthologs derived from STRING results; Table S3: Normalized (10%) pBLAST 2.2.26 score for the SRGs of R. solanacearum; Table S4: Bacterial secretion systems of Rs.

Author Contributions

Conceptualization, P.C.; methodology, P.C., P.B. and G.B.; formal analysis, P.C. and G.B.; resources, A.K.M., F.-S.Q. and S.S.; data source, A.K.M., F.-S.Q.; writing—original draft preparation, P.C., P.B., G.B.; software, P.C. and G.B.; writing—review and editing, G.B., S.S. and P.B.; supervision, P.C. and P.B.; funding acquisition, F.-S.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Research Foundation of Korea (NRF) (2018R1A6A1A03025124) and the Ministry of Health & Welfare, Republic of Korea (HV20C0085, HV20C0142).

Institutional Review Board Statement

This article does need any statement from Institutional Review Board.

Informed Consent Statement

This article does not contain any studies with human participants and/or animals performed by any of the authors. Formal consent was not required in this study.

Data Availability Statement

The genomes of Ralstonia solanacearum were downloaded from the NCBI RefSeq database (https://www.ncbi.nlm.nih.gov/genome/genomes/490) (accessed on 31 March 2022).

Acknowledgments

We are very thankful to M.U.C. Women’s College, University of Burdwan, Bardhaman, West Bengal, India; the University of Calcutta, West Bengal, India; and the University of Illinois Urbana-Champaign, Illinois, USA, for providing necessary support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schandry, N.; de Lange, O.; Prior, P.; Lahaye, T. TALE-like Effectors Are an Ancestral Feature of the Ralstonia solanacearum Species Complex and Converge in DNA Targeting Specificity. Front. Plant Sci. 2016, 7, 1225. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Salanoubat, M.; Genin, S.; Artiguenave, F.; Gouzy, J.; Mangenot, S.; Arlat, M.; Billault, A.; Brottier, P.; Camus, J.C.; Cattolico, L.; et al. Genome sequence of the plant pathogen Ralstonia solanacearum. Nature 2002, 415, 497–502. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Chattopadhyay, P.; Banerjee, G.; Handique, P.J. Use of an abscisic acid-producing Bradyrhizobium japonicum isolate as biocontrol agent against bacterial wilt disease caused by Ralstonia solanacearum. J. Plant Dis. Prot. 2022, 1–11. [Google Scholar] [CrossRef]
  4. Gutarra, L.; Herrera, J.; Fernandez, E.; Kreuze, J.; Lindqvist-Kreuze, H. Diversity, pathogenicity, and current occurrence of bacterial wilt bacterium Ralstonia solanacearum in Peru. Front. Plant Sci. 2017, 8, 1221. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Filloux, A. Bacterial Regulatory Networks; Caister Academic Press: London, UK, 2012. [Google Scholar]
  6. Hengge-Aronis, R. Signal Transduction and Regulatory Mechanisms Involved in Control of the σ S (RpoS) Subunit of RNA Polymerase. Microbiol. Mol. Biol. Rev. 2002, 66, 373–395. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Colburn-Clifford, J.M.; Scherf, J.M.; Allen, C. Ralstonia solanacearum Dps Contributes to Oxidative Stress Tolerance and to Colonization of and Virulence on Tomato Plants. Appl. Environ. Microbiol. 2010, 76, 7392–7399. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Mandal, S.; Kar, I.; Mukherjee, A.K.; Acharya, P. Elicitor-Induced Defense Responses in Solanum lycopersicumagainst against Ralstonia solanacearum. Sci. World J. 2013, 2013, 561056. [Google Scholar] [CrossRef] [Green Version]
  9. Flores-Cruz, Z.; Allen, C. Necessity of OxyR for the Hydrogen Peroxide Stress Response and Full Virulence in Ralstonia solanacearum. Appl. Environ. Microbiol. 2011, 77, 6426–6432. [Google Scholar] [CrossRef] [Green Version]
  10. Taylor, P.D.; Inchley, C.J.; Gallagher, M.P. The Salmonella typhimurium AhpC Polypeptide Is Not Essential for Virulence in BALB/c Mice but Is Recognized as an Antigen during Infection. Infect. Immun. 1998, 66, 3208–3217. [Google Scholar] [CrossRef] [Green Version]
  11. Vázquez-Torres, A. Redox Active Thiol Sensors of Oxidative and Nitrosative Stress. Antioxid. Redox Signal. 2012, 17, 1201–1214. [Google Scholar] [CrossRef] [Green Version]
  12. Karlinsey, J.E.; Bang, I.-S.; Becker, L.A.; Frawley, E.R.; Porwollik, S.; Robbins, H.F.; Thomas, V.C.; Urbano, R.; McClelland, M.; Fang, F.C. The NsrR regulon in nitrosative stress resistance of Salmonella enterica serovar Typhimurium. Mol. Microbiol. 2012, 85, 1179–1193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Crack, J.; Green, J.; Thomson, A.J. Mechanism of Oxygen Sensing by the Bacterial Transcription Factor Fumarate-Nitrate Reduction (FNR). J. Biol. Chem. 2004, 279, 9278–9286. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Lobato, L.; Bouzhir-Sima, L.; Yamashita, T.; Wilson, M.T.; Vos, M.H.; Liebl, U. Dynamics of the heme-binding bacterial gas-sensing dissimilative nitrate respiration regulator (DNR) and activation barriers for ligand binding and escape. J. Biol. Chem. 2014, 289, 26514–26524. [Google Scholar] [CrossRef] [Green Version]
  15. Raivio, T.L. MicroReview: Envelope stress responses and Gram-negative bacterial pathogenesis. Mol. Microbiol. 2005, 56, 1119–1128. [Google Scholar] [CrossRef] [PubMed]
  16. Humphreys, S.; Rowley, G.; Stevenson, A.; Anjum, M.F.; Woodward, M.J.; Gilbert, S.; Kormanec, J.; Roberts, M. Role of the Two-Component Regulator CpxAR in the Virulence of Salmonella enterica Serotype Typhimurium. Infect. Immun. 2004, 72, 4654–4661. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Humphreys, S.; Stevenson, A.; Bacon, A.; Weinhardt, A.B.; Roberts, M. The Alternative Sigma Factor, ς E, Is Critically Important for the Virulence of Salmonella typhimurium. Infect. Immun. 1999, 67, 1560–1568. [Google Scholar] [CrossRef]
  18. Fang, F.C.; Libby, S.J.; Buchmeier, N.A.; Loewen, P.C.; Switala, J.; Harwood, J.; Guiney, D.G. The alternative sigma factor katF (rpoS) regulates Salmonella virulence. Proc. Natl. Acad. Sci. USA 1992, 89, 11978–11982. [Google Scholar] [CrossRef] [Green Version]
  19. Shinagawa, H. SOS response as an adaptive response to DNA damage in prokaryotes. In Stress-Inducible Cellular Responses; Birkhäuser: Basel, Switzerland, 1996; pp. 221–235. [Google Scholar] [CrossRef]
  20. Waldor, M.K.; Friedman, D.I. Phage regulatory circuits and virulence gene expression. Curr. Opin. Microbiol. 2005, 8, 459–465. [Google Scholar] [CrossRef]
  21. Njoroge, J.W.; Nguyen, Y.; Curtis, M.M.; Moreira, C.G.; Sperandio, V. Virulence Meets Metabolism: Cra and KdpE Gene Regulation in Enterohemorrhagic Escherichia coli. mBio 2012, 3, e00280-12. [Google Scholar] [CrossRef] [Green Version]
  22. Yoon, H.; McDermott, J.E.; Porwollik, S.; McClelland, M.; Heffron, F. Coordinated Regulation of Virulence during Systemic Infection of Salmonella enterica Serovar Typhimurium. PLoS Pathog. 2009, 5, e1000306. [Google Scholar] [CrossRef] [Green Version]
  23. Troxell, B.; Sikes, M.L.; Fink, R.C.; Vazquez-Torres, A.; Jones-Carson, J.; Hassan, H.M. Fur Negatively Regulates hns and Is Required for the Expression of HilA and Virulence in Salmonella enterica Serovar Typhimurium. J. Bacteriol. 2011, 193, 497–505. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Patzer, S.I.; Hantke, K. The Zinc-responsive Regulator Zur and Its Control of theznu Gene Cluster Encoding the ZnuABC Zinc Uptake System in Escherichia coli. J. Biol. Chem. 2000, 275, 24321–24332. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Kinkel, T.L.; Roux, C.M.; Dunman, P.M.; Fang, F.C. The Staphylococcus aureus SrrAB Two-Component System Promotes Resistance to Nitrosative Stress and Hypoxia. mBio 2013, 4, e00696-13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Durand, S.; Braun, F.; Lioliou, E.; Romilly, C.; Helfer, A.-C.; Kuhn, L.; Quittot, N.; Nicolas, P.; Romby, P.; Condon, C. A Nitric Oxide Regulated Small RNA Controls Expression of Genes Involved in Redox Homeostasis in Bacillus subtilis. PLoS Genet. 2015, 11, e1004957. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Alvarez, A.F.; Rodriguez, C.; Georgellis, D. Ubiquinone and Menaquinone Electron Carriers Represent the Yin and Yang in the Redox Regulation of the ArcB Sensor Kinase. J. Bacteriol. 2013, 195, 3054–3061. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Sardiwal, S.; Kendall, S.L.; Movahedzadeh, F.; Rison, S.C.G.; Stoker, N.G.; Djordjevic, S. A GAF Domain in the Hypoxia/NO-inducible Mycobacterium tuberculosis DosS Protein Binds Haem. J. Mol. Biol. 2005, 353, 929–936. [Google Scholar] [CrossRef]
  29. Karlinsey, J.E.; Maguire, M.E.; Becker, L.A.; Crouch, M.-L.V.; Fang, F.C. The phage shock protein PspA facilitates divalent metal transport and is required for virulence of Salmonella enterica sv. Typhimurium. Mol. Microbiol. 2010, 78, 669–685. [Google Scholar] [CrossRef] [Green Version]
  30. Pagels, M.; Fuchs, S.; Pané-Farré, J.; Kohler, C.; Menschner, L.; Hecker, M.; McNamarra, P.J.; Bauer, M.C.; Von Wachenfeldt, C.; Liebeke, M.; et al. Redox sensing by a Rex-family repressor is involved in the regulation of anaerobic gene expression in Staphylococcus aureus. Mol. Microbiol. 2010, 76, 1142–1161. [Google Scholar] [CrossRef] [Green Version]
  31. Darwin, A.J.; Stewart, V. Nitrate and Nitrite Regulation of the Fnr-dependentaeg-46.5Promoter ofEscherichia coliK-12 is Mediated by Competition Between Homologous Response Regulators (NarL and NarP) for a Common DNA-binding Site. J. Mol. Biol. 1995, 251, 15–29. [Google Scholar] [CrossRef]
  32. Noriega, C.E.; Lin, H.-Y.; Chen, L.-L.; Williams, S.B.; Stewart, V. Asymmetric cross-regulation between the nitrate-responsive NarX-NarL and NarQ-NarP two-component regulatory systems from Escherichia coli K-12. Mol. Microbiol. 2010, 75, 394–412. [Google Scholar] [CrossRef] [Green Version]
  33. Kondoh, H.; Ball, C.B.; Adler, J. Identification of a methyl-accepting chemotaxis protein for the ribose and galactose chemoreceptors of Escherichia coli. Proc. Natl. Acad. Sci. USA 1979, 76, 260–264. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Greer-Phillips, S.E.; Alexandre, G.; Taylor, B.L.; Zhulin, I.B. Aer and Tsr guide Escherichia coli in spatial gradients of oxidizable substrates. Microbiology 2003, 149, 2661–2667. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Falconi, M.; Colonna, B.; Prosseda, G.; Micheli, G.; Gualerzi, C.O. Thermoregulation of Shigella and Escherichia coli EIEC pathogenicity. A temperature-dependent structural transition of DNA modulates accessibility of virF promoter to transcriptional repressor H-NS. EMBO J. 1998, 17, 7033–7043. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Ono, S.; Goldberg, M.D.; Olsson, T.; Esposito, D.; Hinton, J.C.D.; Ladbury, J.E. H-NS is a part of a thermally controlled mechanism for bacterial gene regulation. Biochem. J. 2005, 391, 203–213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Kamp, H.D.; Higgins, D.E. A Protein Thermometer Controls Temperature-Dependent Transcription of Flagellar Motility Genes in Listeria monocytogenes. PLoS Pathog. 2011, 7, e1002153. [Google Scholar] [CrossRef]
  38. Prugnola, A.; Arico, B.; Manetti, R.; Rappuoli, R.; Scarlato, V. Response of the bvg regulon of Bordetella pertussis to different temperatures and short-term temperature shifts. Microbiology 1995, 141, 2529–2534. [Google Scholar] [CrossRef] [Green Version]
  39. Herbst, K.; Bujara, M.; Heroven, A.K.; Opitz, W.; Weichert, M.; Zimmermann, A.; Dersch, P. Intrinsic Thermal Sensing Controls Proteolysis of Yersinia Virulence Regulator RovA. PLoS Pathog. 2009, 5, e1000435. [Google Scholar] [CrossRef]
  40. Quade, N.; Mendonca, C.; Herbst, K.; Heroven, A.K.; Ritter, C.; Heinz, D.W.; Dersch, P. Structural Basis for Intrinsic Thermosensing by the Master Virulence Regulator RovA of Yersinia. J. Biol. Chem. 2012, 287, 35796–35803. [Google Scholar] [CrossRef] [Green Version]
  41. Reines, A.E.; Deller, A.T. Parsec-scale radio emission from the low-luminosity active galactic nucleus in the dwarf starburst galaxy Henize 2–10. Astrophys. J. 2012, 750, L24. [Google Scholar] [CrossRef] [Green Version]
  42. Chakraborty, S.; Mizusaki, H.; Kenney, L.J. A FRET-Based DNA Biosensor Tracks OmpR-Dependent Acidification of Salmonella during Macrophage Infection. PLoS Biol. 2015, 13, e1002116. [Google Scholar] [CrossRef] [Green Version]
  43. Marcus, E.A.; Sachs, G.; Wen, Y.; Scott, D.R. Phosphorylation-dependent and Phosphorylation-independent Regulation ofHelicobacter pyloriAcid Acclimation by the ArsRS Two-component System. Helicobacter 2015, 21, 69–81. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Gore, A.L.; Payne, S.M. CsrA and Cra Influence Shigella flexneri Pathogenesis. Infect. Immun. 2010, 78, 4674–4682. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Converse, P.D.; DeShon, R.P. A tale of two tasks: Reversing the self-regulatory resource depletion effect. J. Appl. Psychol. 2009, 94, 1318–1324. [Google Scholar] [CrossRef] [PubMed]
  46. Boyd, J.; Oza, M.N.; Murphy, J.R. Molecular cloning and DNA sequence analysis of a diphtheria tox iron-dependent regulatory element (dtxR) from Corynebacterium diphtheriae. Proc. Natl. Acad. Sci. USA 1990, 87, 5968–5972. [Google Scholar] [CrossRef] [Green Version]
  47. Lieser, S.A.; Davis, T.C.; Helmann, J.D.; Cohen, S.M. DNA-Binding and Oligomerization Studies of the Manganese(II) Metalloregulatory Protein MntR from Bacillus subtilis. Biochemistry 2003, 42, 12634–12642. [Google Scholar] [CrossRef]
  48. Banerjee, G.; Basak, S.; Roy, T.; Chattopadhyay, P. Intrinsic role of bacterial secretion systems in phylogenetic niche conservation of Bradyrhizobium spp. FEMS Microbiol. Ecol. 2019, 95, fiz165. [Google Scholar] [CrossRef]
  49. Fang, F.C.; Frawley, E.R.; Tapscott, T.; Vázquez-Torres, A. Bacterial Stress Responses during Host Infection. Cell Host Microbe 2016, 20, 133–143. [Google Scholar] [CrossRef] [Green Version]
  50. Lage, C.; Leitao, A. Membrane permeability and sensitivity to lethal heat are affected by lexA and recA mutations in Escherichia coli K12. Braz. J. Med. Biol. Res. Rev. Bras. De Pesqui. Med. E Biol. 1994, 27, 2383–2389. [Google Scholar]
  51. Mozaheb, N.; Mingeot-Leclercq, M.-P. Membrane vesicle production as a bacterial defense against stress. Front. Microbiol. 2020, 11, 600221. [Google Scholar] [CrossRef]
  52. Chen, J.; Li, S.; Luo, J.; Zhang, Y.; Ding, W. Graphene oxide induces toxicity and alters energy metabolism and gene expression in Ralstonia solanacearum. J. Nanosci. Nanotechnol. 2017, 17, 186–195. [Google Scholar] [CrossRef]
  53. Tomlinson, D.L.; Elphinstone, J.G.; Soliman, M.Y.; Hanafy, M.S.; Shoala, T.M.; Abd El-Fatah, H.; Agag, S.H.; Kamal, M.; Abd El-Aliem, M.M.; Fawzi, F.G.; et al. Recovery of Ralstonia solanacearum from canal water in traditional potato-growing areas of Egypt but not from designated Pest-Free Areas (PFAs). Eur. J. Plant Pathol. 2009, 125, 589–601. [Google Scholar] [CrossRef]
  54. Wood, J.M. Osmotic stress. In Bacterial Stress Responses, 2nd ed.; Storz, G., Hengge, R., Eds.; Wiley: Hoboken, NJ, USA, 2010; pp. 133–156. [Google Scholar] [CrossRef]
  55. MacIntyre, A.M.; Barth, J.X.; Pellitteri Hahn, M.C.; Scarlett, C.O.; Genin, S.; Allen, C. Trehalose synthesis contributes to osmotic stress tolerance and virulence of the bacterial wilt pathogen Ralstonia solanacearum. Mol. Plant-Microbe Interact. 2020, 33, 462–473. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. He, M.; Shah Jahan, M.; Wang, Y.; Sun, J.; Shu, S.; Guo, S. Compost amendments based on vinegar residue promote tomato growth and suppress bacterial wilt caused by Ralstonia Solanacearum. Pathogens 2020, 9, 227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Yang, L.; Wei, Z.; Li, S.; Xiao, R.; Xu, Q.; Ran, Y.; Ding, W. Plant secondary metabolite, daphnetin reduces extracellular polysaccharides production and virulence factors of Ralstonia solanacearum. Pestic. Biochem. Physiol. 2021, 179, 104948. [Google Scholar] [CrossRef]
  58. Poole, R.K. Nitric oxide and nitrosative stress tolerance in bacteria. Biochem. Soc. Trans. 2005, 33, 176–180. [Google Scholar] [CrossRef]
  59. Truchon, A.N.; Hendrich, C.G.; Bigott, A.F.; Dalsing, B.L.; Allen, C. NorA, HmpX, and NorB cooperate to reduce NO toxicity during denitrification and plant pathogenesis in Ralstonia solanacearum. Microbiol. Spectr. 2022, 10, e0026422. [Google Scholar] [CrossRef]
  60. Jeong, K.C.; Hung, K.F.; Baumler, D.J.; Byrd, J.J.; Kaspar, C.W. Acid stress damage of DNA is prevented by Dps binding in Escherichia coliO157:H7. BMC Microbiol. 2008, 8, 181. [Google Scholar] [CrossRef] [Green Version]
  61. Richardson, A.E.; Barea, J.-M.; McNeill, A.M.; Prigent-Combaret, C. Acquisition of phosphorus and nitrogen in the rhizosphere and plant growth promotion by microorganisms. Plant Soil 2009, 321, 305–339. [Google Scholar] [CrossRef]
  62. Henard, C.A.; Tapscott, T.; Crawford, M.A.; Husain, M.; Doulias, P.-T.; Porwollik, S.; Liu, L.; McClelland, M.; Ischiropoulos, H.; Vázquez-Torres, A. The 4-cysteine zinc-finger motif of the RNA polymerase regulator DksA serves as a thiol switch for sensing oxidative and nitrosative stress. Mol. Microbiol. 2014, 91, 790–804. [Google Scholar] [CrossRef] [Green Version]
  63. Brown, D.G.; Allen, C. Ralstonia solanacearum genes induced during growth in tomato: An inside view of bacterial wilt. Mol. Microbiol. 2004, 53, 1641–1660. [Google Scholar] [CrossRef]
  64. Castillo, J.A.; Greenberg, J.T. Evolutionary Dynamics of Ralstonia solanacearum. Appl. Environ. Microbiol. 2007, 73, 1225–1238. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Feil, E.J.; Spratt, B.G. Recombination and the Population Structures of Bacterial Pathogens. Annu. Rev. Microbiol. 2001, 55, 561–590. [Google Scholar] [CrossRef] [PubMed]
  66. Puigvert, M.; Guarischi-Sousa, R.; Zuluaga, P.; Coll, N.S.; Macho, A.P.; Setubal, J.C.; Valls, M. Transcriptomes of Ralstonia solanacearum during root colonization of Solanum commersonii. Front. Plant Sci. 2017, 8, 370. [Google Scholar] [CrossRef] [Green Version]
  67. Wang, Z.; Li, S.; Li, F.; Yang, H.; Yang, F.; Xiang, J. Characterization of two types of vascular endothelial growth factor from Litopenaeus vannamei and their involvements during WSSV infection. Fish Shellfish. Immunol. 2015, 47, 824–832. [Google Scholar] [CrossRef] [PubMed]
  68. Kim, K.K.; Kim, H.B. Protein interaction network related to Helicobacter pylori infection response. World J. Gastroenterol. 2009, 15, 4518. [Google Scholar] [CrossRef] [Green Version]
  69. Marchadier, E.; Carballido-López, R.; Brinster, S.; Fabret, C.; Mervelet, P.; Bessières, P.; Noirot-Gros, M.-F.; Fromion, V.; Noirot, P. An expanded protein-protein interaction network in Bacillus subtilis reveals a group of hubs: Exploration by an integrative approach. Proteomics 2011, 11, 2981–2991. [Google Scholar] [CrossRef]
  70. Van Gijsegem, F.; Genin, S.; Boucher, C. Conservation of secretion pathways for pathogenicity determinants of plant and animal bacteria. Trends Microbiol. 1993, 1, 175–180. [Google Scholar] [CrossRef]
  71. Sikdar, R.; Simmons, A.R.; Doerrler, W.T. Multiple envelope stress response pathways are activated in an Escherichia coli strain with mutations in two members of the DedA membrane protein family. J. Bacteriol. 2013, 195, 12–24. [Google Scholar] [CrossRef] [Green Version]
  72. da Cruz Nizer, W.S.; Inkovskiy, V.; Versey, Z.; Strempel, N.; Cassol, E.; Overhage, J. Oxidative Stress Response in Pseudomonas aeruginosa. Pathogens 2021, 10, 1187. [Google Scholar] [CrossRef]
  73. González, E.T.; Brown, D.G.; Swanson, J.K.; Allen, C. Using the Ralstonia solanacearum Tat Secretome To Identify Bacterial Wilt Virulence Factors. Appl. Environ. Microbiol. 2007, 73, 3779–3786. [Google Scholar] [CrossRef] [Green Version]
  74. Sarvas, M. Post-translocational folding of secretory proteins in Gram-positive bacteria. Biochim. Biophys. Acta (BBA)-Mol. Cell Res. 2004, 1694, 311–327. [Google Scholar] [CrossRef] [PubMed]
  75. Wickner, S. Posttranslational Quality Control: Folding, Refolding, and Degrading Proteins. Science 1999, 286, 1888–1893. [Google Scholar] [CrossRef] [PubMed]
  76. Lieberman, J.A.; Frost, N.A.; Hoppert, M.; Fernandes, P.J.; Vogt, S.L.; Raivio, T.L.; Blanpied, T.A.; Donnenberg, M.S. Outer membrane targeting, ultrastructure, and single molecule localization of the enteropathogenic Escherichia coli type IV pilus secretin BfpB. J. Bacteriol. 2012, 194, 1646–1658. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Hendrich, C.G.; Truchon, A.N.; Dalsing, B.L.; Allen, C. Nitric Oxide Regulates the Ralstonia solanacearum Type 3 Secretion System. Microbiol. Spectr. 2022, 10, e00264-22. [Google Scholar]
  78. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2018, 47, D607–D613. [Google Scholar] [CrossRef] [Green Version]
  79. Pramanik, K.; Ghosh, P.K.; Ray, S.; Sarkar, A.; Mitra, S.; Maiti, T.K. An in silico structural, functional and phylogenetic analysis with three dimensional protein modeling of alkaline phosphatase enzyme of Pseudomonas aeruginosa. J. Genet. Eng. Biotechnol. 2017, 15, 527–537. [Google Scholar] [CrossRef]
  80. Hammer, Ø.; Harper, D.A.; Ryan, P. Paleontological statistics software package for education and data analysis. Paleontol. Electron. 2001, 4, 1–9. [Google Scholar]
Figure 1. STRING analysis results showing protein–protein interaction network in stress-responsive genes of R. solanacerium.
Figure 1. STRING analysis results showing protein–protein interaction network in stress-responsive genes of R. solanacerium.
Pathogens 11 00730 g001
Figure 2. Heatmap showing the distribution of stress-responsive genes among R. solanacerium genomes. Row min = 0; row max = 1.
Figure 2. Heatmap showing the distribution of stress-responsive genes among R. solanacerium genomes. Row min = 0; row max = 1.
Pathogens 11 00730 g002
Figure 3. STRING analysis results showing protein–protein interaction network within stress-responsive genes and type III secretion system (T3SS) of R. solanacerium.
Figure 3. STRING analysis results showing protein–protein interaction network within stress-responsive genes and type III secretion system (T3SS) of R. solanacerium.
Pathogens 11 00730 g003
Figure 4. STRING analysis results showing protein–protein interaction network within stress-responsive genes and general secretory protein system (GSPS or Sec system) of R. solanacerium.
Figure 4. STRING analysis results showing protein–protein interaction network within stress-responsive genes and general secretory protein system (GSPS or Sec system) of R. solanacerium.
Pathogens 11 00730 g004
Figure 5. STRING analysis results showing protein–protein interaction network within stress-responsive genes and multiple (T4P, Tat, and Tad) pili systems of R. solanacerium.
Figure 5. STRING analysis results showing protein–protein interaction network within stress-responsive genes and multiple (T4P, Tat, and Tad) pili systems of R. solanacerium.
Pathogens 11 00730 g005
Figure 6. (a) Diagram showing interactions among clusters of stress-responsive genes (SRGs) of R. solanacerium: MSR (membrane stress response), OSR (osmotic stress response), OxSR (oxidative stress response), NxSR (nitrosative stress response), and DdSR (DNA damage stress response); thin lines with “X” show the absence of interactions; thick lines show the presence of interactions. (b) Diagram showing interactions of clusters of stress-responsive genes (SRGs) with the secretion systems of R. solanacerium: Tad (tight adherence), Tat (twin arginine translocation), T4P (type-4 pili), T3SS (type-3 secretion system), and GPSS (general protein secretion system); thick lines show interactions among SRGs clusters; thin lines depict the interaction between SRGs and secretion system.
Figure 6. (a) Diagram showing interactions among clusters of stress-responsive genes (SRGs) of R. solanacerium: MSR (membrane stress response), OSR (osmotic stress response), OxSR (oxidative stress response), NxSR (nitrosative stress response), and DdSR (DNA damage stress response); thin lines with “X” show the absence of interactions; thick lines show the presence of interactions. (b) Diagram showing interactions of clusters of stress-responsive genes (SRGs) with the secretion systems of R. solanacerium: Tad (tight adherence), Tat (twin arginine translocation), T4P (type-4 pili), T3SS (type-3 secretion system), and GPSS (general protein secretion system); thick lines show interactions among SRGs clusters; thin lines depict the interaction between SRGs and secretion system.
Pathogens 11 00730 g006
Figure 7. Schematic diagram of the study protocol.
Figure 7. Schematic diagram of the study protocol.
Pathogens 11 00730 g007
Table 1. Reference dataset for bacterial stress-responsive enzymes and factors.
Table 1. Reference dataset for bacterial stress-responsive enzymes and factors.
Sl. No.RegulatorStress ResponseReported CandidatesReferences
1OxyROxidative stressSalmonella enterica[10]
2DksAOxidative stress-[11]
3NsrRNOS. enterica[12]
4FNROxygen limitationE. coli[13]
5DNRNitrosative stress-[14]
6BaeMembrane stress-[15]
7CpxMisfolded OMPs at acid pHS. enterica[16]
8RcsMembrane stress-[15]
9σEMembrane stressS. enterica[17]
10RseAMembrane stress-[18]
11DegSMembrane stress-[18]
12RecADNA damage-[19]
13LexADNA damageVibrio cholerae[20]
14FruRCarbohydrate starvationE. coli[21]
15SpoTAmino acid starvationS. enterica[22]
16FURIron starvationS. enterica[23]
17ZURZinc starvationS. enterica[24]
18SrrA/SrrBOxygen limitationS. enterica[25]
19ResD/ResEOxygen limitationS. enterica[26]
20ArcB/ArcAOxygen limitationE. coli[27]
21DosR/DosSOxygen limitationM. tuberculosis[28]
22PspEnvelop stressS. enterica[29]
23RexOxygen limitationS. aureus[30]
24NarQ/NarPOxygen limitationS. enterica[31]
25NarX/NarLOxygen limitationE. coli[32]
26TrgOxygen limitation-[33]
27TsrOxygen limitationE. coli[34]
28AerOxygen limitationE. coli[34]
29H-NSTemperatureS. flexneri[35]
30HhaTemperatureE. coli[36]
31GmaRTemperatureL. monocytogenes[37]
32MogRTemperatureL. monocytogenes[37]
33BvgA/BvgSTemperatureB. pertussis[38]
34RovATemperatureY. pestis[39]
35SlyATemperatureS. enterica[40]
36PhoP/PhoQTemperatureY. enterocolitica[41]
37PmrA/PmrBTemperatureY. enterocolitica[41]
38OmpREnvelopeS. enterica[42]
39DpsNutrient limitationR. solanacearum[7]
40ArsAAcid pHH. pylori[43]
41BvgATemperature regulatesB. pertussis[38]
42CsrAStarvationS. flexneri[44]
43DosR (NosR)Hypoxia/NO.M. tuberculosis[45]
44DtxRIron deprivationC. diphtheriae[46]
45MntRManganese starvationBacillus subtilis[47]
46GmaRTemperatureListeria monocytogenes[37]
47σSStarvationS. enterica[18]
Table 2. Functional classification of stress-responsive genes (SRG) of Rs.
Table 2. Functional classification of stress-responsive genes (SRG) of Rs.
Sl. No.SRGKO No.FunctionCluster
1nlpDK06194Lipoprotein NlpD | (GenBank) putative outer membrane metallopeptidase lipoprotein nlpD/MSR
2degQK04772Serine protease DegQ (EC:3.4.21.-) | (GenBank) putative serine protease do-like precursor (degP)
3rseBK03598Sigma-E factor negative regulatory protein RseB | (GenBank) rseB; Sigma-E factor regulatory protein
4rseAK03597Sigma-E factor negative regulatory protein RseA | (GenBank) rseA; Sigma-E factor negative regulatory protein
5rpoEK03088RNA polymerase sigma-70 factor, ECF subfamily (misfolded OMPs at acid pH; required for virulence) | (GenBank) rpoE
6lexAK01356Repressor LexA (EC:3.4.21.88) (Regulates toxin production) | (GenBank) lexA
7recXK03565Regulatory protein (Modulates recA activity) | (GenBank) recX
8recAK03553Catalyzes the hydrolysis of ATP in the presence of single-stranded DNA, RecA | (GenBank) recAn
9baeRK07664Two-component system, OmpR family, response regulator BaeR for mdtABCD and acrD | (GenBank) baeROSR
10baeSK07642Two-component system, OmpR family, sensor histidine kinase BaeS (EC:2.7.13.3) | (GenBank) baeS
11ompRK07659Two-component system, OmpR family, phosphate regulon response regulator OmpR | (GenBank) ompR
12ompRK07638Two-component system, OmpR family (EC:2.7.13.3) | (GenBank) putative sensory histidine kinase in Two-component regulatory system with OmpR
13envZK07638Two-component system, OmpR family, osmolarity sensor histidine kinase EnvZ (EC:2.7.13.3) | (GenBank) envZ
14rcsBK07687Two-component system, NarL family, captular synthesis response regulator RcsB | (GenBank) putative response regulator receiver
15KataK03781 Catalase (EC:1.11.1.6) | (GenBank) kataOxSR
16katBK03781Catalase [EC:1.11.1.6] (Paraquat-inducible catalase isozyme B) | (GenBank) katB
17katGK03782Catalase-peroxidase (EC:1.11.1.21) | (GenBank) katG
18oxyRK04761LysR family transcriptional regulator, hydrogen peroxide-inducible genes activator | (GenBank) oxyR
19sodBK04564 Superoxide dismutase, Fe-Mn family (EC:1.15.1.1) | (GenBank) sodB
20sodCK04565Superoxide dismutase, Cu-Zn family (EC:1.15.1.1) | (GenBank) sodC
21DpsK04047Starvation-inducible DNA-binding protein | (GenBank) putative DNA-binding protein
22nsrRK13771Rrf2 family transcriptional regulator, nitric oxide-sensitive transcriptional repressor (regulates genes required for virulence) | (GenBank) nsrRNxSR
23norBK04561Nitric oxide reductase subunit B (transmembrane protein) (EC:1.7.2.5) | (GenBank) norB
24aniAK00368Nitrite reductase (NO-forming) (EC:1.7.2.1) | (GenBank) ainA
25FnrK01420CRP/FNR family transcriptional regulator, anaerobic regulatory protein | (GenBank) fnrDdSR
26DnrK21563CRP/FNR family transcriptional regulator, dissimilatory nitrate respiration regulator | (GenBank) putative transcription regulator protein
27dskAK06204 DnaK suppressor protein | (GenBank) dskA
28rpoZK03060DNA-directed RNA polymerase subunit omega (EC:2.7.7.6) | (GenBank) rpoZ
29rpoSK03087RNA polymerase nonessential primary-like sigma factor (regulates expression of plasmid virulence genes) | (GenBank) rpoS
30rpoBK03043DNA-directed RNA polymerase subunit beta (EC:2.7.7.6) | (GenBank) RNA polymerase, beta subunit
31rpoAK03040DNA-directed RNA polymerase subunit alpha (EC:2.7.7.6) | (GenBank) RNA polymerase, alpha subunit
32polAK02335DNA polymerase I (EC:2.7.7.7) | (GenBank) polA
33h-nsK03746DNA-binding protein H-NS | (GenBank) Histone-like nucleoid-structuring protein H-NS
34FurK03711Fur family transcriptional regulator, ferric uptake regulator | (GenBank) fur
35GmkK00942Guanylate kinase | (GenBank) gmk
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Banerjee, G.; Quan, F.-S.; Mondal, A.K.; Sur, S.; Banerjee, P.; Chattopadhyay, P. Interrelation between Stress Management and Secretion Systems of Ralstonia solanacearum: An In Silico Assessment. Pathogens 2022, 11, 730. https://doi.org/10.3390/pathogens11070730

AMA Style

Banerjee G, Quan F-S, Mondal AK, Sur S, Banerjee P, Chattopadhyay P. Interrelation between Stress Management and Secretion Systems of Ralstonia solanacearum: An In Silico Assessment. Pathogens. 2022; 11(7):730. https://doi.org/10.3390/pathogens11070730

Chicago/Turabian Style

Banerjee, Goutam, Fu-Shi Quan, Amit Kumar Mondal, Shantanu Sur, Pratik Banerjee, and Pritam Chattopadhyay. 2022. "Interrelation between Stress Management and Secretion Systems of Ralstonia solanacearum: An In Silico Assessment" Pathogens 11, no. 7: 730. https://doi.org/10.3390/pathogens11070730

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