Differential transcriptional profile of Corynebacterium pseudotuberculosis in response to abiotic stresses
- Anne Cybelle Pinto1,
- Pablo Henrique Caracciolo Gomes de Sá2,
- Rommel T J Ramos2,
- Silvanira Barbosa2,
- Hivana P Melo Barbosa2,
- Adriana Carneiro Ribeiro2,
- Wanderson Marques Silva1,
- Flávia Souza Rocha1,
- Mariana Passos Santana1,
- Thiago Luiz de Paula Castro1,
- Anderson Miyoshi1,
- Maria P C Schneider2,
- Artur Silva2 and
- Vasco Azevedo1Email author
© Pinto et al.; licensee BioMed Central Ltd. 2014
Received: 4 July 2013
Accepted: 13 December 2013
Published: 9 January 2014
The completion of whole-genome sequencing for Corynebacterium pseudotuberculosis strain 1002 has contributed to major advances in research aimed at understanding the biology of this microorganism. This bacterium causes significant loss to goat and sheep farmers because it is the causal agent of the infectious disease caseous lymphadenitis, which may lead to outcomes ranging from skin injury to animal death. In the current study, we simulated the conditions experienced by the bacteria during host infection. By sequencing transcripts using the SOLiDTM 3 Plus platform, we identified new targets expected to potentiate the survival and replication of the pathogen in adverse environments. These results may also identify possible candidates useful for the development of vaccines, diagnostic kits or therapies aimed at the reduction of losses in agribusiness.
Under the 3 simulated conditions (acid, osmotic and thermal shock stresses), 474 differentially expressed genes exhibiting at least a 2-fold change in expression levels were identified. Important genes to the infection process were induced, such as those involved in virulence, defence against oxidative stress, adhesion and regulation, and many genes encoded hypothetical proteins, indicating that further investigation of the bacterium is necessary. The data will contribute to a better understanding of the biology of C. pseudotuberculosis and to studies investigating strategies to control the disease.
Despite the veterinary importance of C. pseudotuberculosis, the bacterium is poorly characterised; therefore, effective treatments for caseous lymphadenitis have been difficult to establish. Through the use of RNAseq, these results provide a better biological understanding of this bacterium, shed light on the most likely survival mechanisms used by this microorganism in adverse environments and identify candidates that may help reduce or even eradicate the problems caused by this disease.
KeywordsDifferential gene expression Transcripts RNAseq SOLID™ Stress C. pseudotuberculosis
Corynebacterium pseudotuberculosis is a Gram-positive pathogenic bacterium belonging to the class Actinobacteria, which is a member of the Corynebacterium, Mycobacterium, Nocardia and Rhodococcus genera (the CMNR group). The CMNR group shares several common characteristics, including (i) the organisation of the cell wall, which is mainly composed of peptidoglycan, arabinogalactan and mycolic acids and (ii) the high G + C content of the genome (47-74%) .
The bacterium causes caseous lymphadenitis disease, which affects small ruminants and large animals (such as horses and cattle) worldwide and can infect humans . Therefore, there is an urgent need to control the disease through the development of effective vaccines, therapies and diagnostic kits.
Because C. pseudotuberculosis is a facultative intracellular microorganism found preferentially in macrophages in the host, during the infection process, the bacterium is exposed to a number of environmental changes that are far from ideal .
After phagocytosis, the phagosome quickly becomes acidic (pH ~ 4.6-5.0)  negatively affecting the metabolism and damaging macromolecules in the invading cell. In addition, other intracellular stresses negatively affect the microorganism, including oxidative thermal shock and nitrosative, surface, osmotic and starvation stresses; however, the bacterium manages to escape and persist in the environment .
To survive in this environment, the pathogen must mount an immediate and adequate, protective response that is reflected initially by transcriptional changes in specific sets of genes . In this context, sigma factors, which coordinate the expression of these genes under different types of stresses are important ; these factors include sigS in E.coli, which is involved in trehalose synthesis during osmotic stress, with trehalose serving as an important osmoprotectant in this type of stress . In Mycobacterium tuberculosis, RT-PCR (real-time polymerase chain reaction) was used to demonstrate the transcriptional profile of 10 sigma factors during the exponential growth phase. The role of these sigma factors was analysed under different conditions of stress, and a number of the factors demonstrated increased expression in response to one kind of stimulus, whereas others responded to more than one stimulus. The resistance to different environmental stresses is associated with the ability of the pathogenic bacteria to survive in the host, and different sigma factors play fundamental roles in the survival of the pathogen. For example, in M. smegmatis, sigB is involved in the response to oxidative stress, and sigE was shown to play a role in the regulation of genes involved in the response to acid stress, thermal shock and sodium dodecyl sulphate (SDS) exposure [5, 6].
The functions of sigma factors include involvement in the adaptation to stress, the interaction of the bacterium with the extracellular medium and in a number of cases, with bacterial virulence .
For C. pseudotuberculosis 1002, there is still no information regarding the role of sigma factors in the regulation of genes involved in bacterial survival throughout the infection process. In addition, few virulence determinants contributing to bacterial survival have been identified. Therefore, investigations related to the control of caseous lymphadenitis have been difficult . To date, the virulence determinants most studied in C. pseudotuberculosis infection include the following: the PLD (phospholipase D) protein, an exoprotein that is also considered leukotoxic, contributing to the formation of lesions and the destruction of caprine macrophages during infection ;the fagABC operon and the fagD gene, which play a role in the virulence of the bacterium and have been identified as genes involved in iron acquisition ; the high concentration of cell wall lipids, which renders the microorganism resistance to digestion by cellular enzymes and allows it to persist as a facultative intracellular parasite  and CP40, identified as an immunogenic protein that exhibits proteolytic activity as a serine protease .
The availability of the C. pseudotuberculosis strain 1002 genome (access number CP001809) has allowed the further investigation and characterisation of the microorganism, which is poorly characterised despite its importance to agribusiness. Therefore, to generate additional information, the transcriptional profile of C. pseudotuberculosis was analysed using cDNA sequencing with SOLiD™ 3 plus next-generation technology (Life TechnologiesTM, CA). Using this technology, we investigated the molecular characteristics of the microorganism that allow it to persist in the host and the mechanisms the bacteria use to escape the host immune response. Additionally, in an attempt to contribute to the elimination of caseous lymphadenitis in caprine and ovine populations, we investigated a large number of targets related to C. pseudotuberculosis virulence by simulating the specific conditions tolerated by the bacterium when invading the host (i.e., acidic, osmotic and high temperature stresses).
Results and discussion
The sequencing of cDNA is an attractive technology for the investigation of gene expression in prokaryotic organisms because it provides a high level of coverage and high sensitivity for the detection of transcripts at considerably lower costs compared to traditional methods . Currently, sequencing technology is considered the gold standard for the analysis of gene expression levels . Therefore, because the microorganism can survive environmental changes in the host during infection, we analysed the transcriptional profile of C. pseudotuberculosis strain 1002 using RNAseq technology. Acid, osmotic and thermal shock stresses were used to identify genes involved in bacterial tolerance of these unfavourable environments.
Number of total readings obtained during sequencing
Filtered Ribosomal Transcripts
Number of unique readings mapped in the genome of C. pseudotuberculosis 1002 and coverage of the transcripts in the genome
Uniquely mapped readings
The automatic annotation of the genome using the Fgenes software (http://www.softberry.com), followed by manual curation of the C. pseudotuberculosis strain 1002 genome, identified 2,090 coding regions. From the cDNA sequencing performed using SOLiDTM at the beginning of the exponential phase, 2,055 transcripts active in control were identified (equivalent to 98.32% of the transcribed genome) and 35 genes (1.67%) were considered non-transcripts, exhibiting a RPKM value of 0. Under the different stress conditions, 2,065 (98.80%) transcripts were produced under osmotic stress, 2,063 (98.70%) transcripts under thermal stress and 2,064 (97.76%) under acid stress.
Among the differentially expressed genes, the genes exhibiting a 2-fold change in expression (at least 2x relative to the control) were selected for analysis. The fold-change values were calculated based on the RPKM value between the stress and the control, in which a value greater than one indicated induced gene expression and a value less than one indicated repressed gene expression.
In a previous study reported by our group , sequencing of transcripts from C. pseudotuberculosis strain Cp31 (biovar equi) was performed using an Ion Torrent platform (Life Technologies). In this study, a comparative analysis between 2 rRNA depletion methodologies was performed, and the transcripts were submitted to ab initio assembly. Both transcriptomes were then submitted to gene ontology analysis according to biological processes and molecular functions. Data were obtained only under physiological conditions using brain heart infusion (BHI) media. The authors observed that few transcripts represented genes involved in pathogenicity or the cell adhesion process and concluded that contact with the host may influence the induction of these transcripts.
The present study simulated some of the environmental conditions that the pathogen faces in the host. The cell adhesion process was among those most represented in the conditions described below, and these levels showed greater than a 2-fold change. This result indicated that there were a greater number of transcripts representing genes that participate in the cell adhesion process when compared to the physiological condition. Together, these findings demonstrate that contact with the host most likely influences the transcription of genes essential for bacterial survival.
Genes induced in the biological processes of the osmotic stress stimulon
The majority of the genes induced under osmotic stress were part of the oxidoreduction process, represented by 4 genes (see Additional file 1: Figure S1). Normal aerobic metabolism induces the production of active oxygen molecules, which are increased following exposure to certain environments . The results demonstrated that the growth of the bacterium was reduced but was not interrupted under these conditions; therefore, the bacterium survived in the environment. This observation was confirmed by the biosynthetic process, which was comprised of 3 genes (see Additional file 1: Figure S1), indicating that the bacterium was able to invest energy into the replication process and survived in the unfavourable environment.
The adhesion process, comprising 3 genes, was prominent under osmotic stress. Adhesion is essential for the initiation of the infectious process because the bacterium-host interaction establishes a pathogenic relationship. The Cp1002_0988 gene, encoding a hypothetical protein, is located within 1 of the pathogenicity islands of C. pseudotuberculosis 1002, indicating its importance in the development of the disease. Additionally, this gene exhibited a 6.8-fold change in expression levels compared with the control (see Additional file 2: Table S1). The Cp1002_1764 and Cp1002_1765 genes, also identified as encoding hypothetical proteins, exhibited 2.7-fold and 4-fold changes in expression, respectively, compared with the control.
Pathogenic bacteria have developed highly sophisticated signal transduction systems that control the coordinated expression of a number of virulence determinants in response to environmental stresses, and changes in osmolarity contribute to the expression of the genes ; therefore, further studies to identify the proteins encoded by these genes and to evaluate the true contribution of these proteins will be necessary.
Genes induced in the biological processes of the acid stress stimulon
The Cp1002_2043 gene exhibited a 7.8-fold change in expression (see Additional file 4: Table S2). This gene, which may encode the Dps protein, was linked to the processes of stress response and cellular oxidoreduction (see Additional file 3: Figure S2). This protein protects the organism against oxidative stress because it stores iron in a bioavailable form, reducing the possibility of the production of reactive oxygen molecules. Under acidic conditions, the production of molecules that produce reactive oxygen species is increased. Reports have demonstrated that in cells in which the pH is decreased, the ratio of HOO– to O2– increases , increasing the chance of producing hydrogen peroxide (H2O2). Therefore, the increased number of genes constituting the cellular oxidoreduction process is justified. Another important gene in the oxidoreduction process was Cp1002_0173, which exhibited a 4-fold change in expression and is presumed to encode a catalase that plays a role in reducing the concentration of H2O2 in the cell.
The mechanisms used by the bacterium to resist the damage caused by reactive oxygen species are essential for survival within the macrophage . Therefore, these proteins might promote the survival of the bacterium in the acid medium starting at the early exponential growth phase.
The Cp1002_1192 msrB (methionine sulphoxide peptide reductase) gene, which was 1 of the genes present in the oxidoreduction process, exhibited a high (16-fold) change in expression. Because the MsrA protein may be required to maintain the role of adhesins, the high increase in msrB expression suggested that the protein contributes to the survival of the pathogen in the host, the resistance to oxidative stress in vitro and the adhesion capability of eukaryotic cells . The methionine (Met) residues in proteins exposed at the cell surface are thought to be involved in capturing reactive oxygen species, and a complex comprised of MsrA and MsrB reduces the oxidised Met residues, removing the reactive oxygen species . Met is the amino acid most sensitive to reactive oxygen species, and the oxidation of the Met residue in a protein alters the protein structure (or prevents translation), drastically affecting the function . There is evidence that only the MsrB domain is present in C. pseudotuberculosis 1002. However, the loss of the msrA-msrB domain or msrB alone in Helicobacter pylori resulted in a reduction in virulence in mouse models, likely because of the oxidation of important proteins . Therefore, the process of cellular oxidoreduction may involve genes that contribute predominantly to the maintenance and persistence of the bacterium in media harmful to the cell.
Genes induced in the biological processes of the thermal shock stimulon
The adhesion process comprises the Cp1002_1765 gene, which may encode a secreted protein. Cp1002_1765 exhibited a 2-fold change in expression compared with the control, and it is an important candidate for studies related to controlling the disease caseous lymphadenitis.
Among the genes involved in the process of the response to stress, the Cp1002_1895 gene, which may encode a heat shock regulatory protein (HspR), exhibited the largest fold-change in expression (4x relative to the control). This protein acts as a negative regulator of the expression of genes encoding chaperones and proteases in different bacteria under physiological conditions. The heat shock proteins play a key role in cellular metabolism under all growth conditions, monitoring the folding, assembly and translocation of cellular proteins . In M. tuberculosis, HspR represses the operon formed by the dnaK-grpE-dnaJ-hspR genes through interaction at the HAIR (HspR-associated inverted repeats) region located in the 5′UTR region of the genes. In C. pseudotuberculosis, the hspR gene is located in the reverse strand, below the dnaJ, GrpE and dnaK genes, indicative of their regulation by HspR.
It has been suggested that the HspR protein acts as a repressor of other genes linked to virulence/pathogenicity . A study in M. tuberculosis demonstrated that the partial disruption of heat-shock regulation influences virulence because the bacterium loses the ability to establish a chronic infection . Mutation of hspR in the organism produced increased expression of the DnaK chaperone, which is highly antigenic, resulting in an enhanced immune response in the host. Furthermore, it was demonstrated that DnaK acts as co-repressor for HspR; the activity of HspR is dependent on DnaK. A proposed mechanism for the regulation of hspR expression is that under conditions of heat shock, the operon is induced, leading to the increased synthesis of DnaK and HspR. When the concentration of these proteins reaches a critical level, they bind to the promoter of the operon, resulting in repression . To investigate whether this scenario occurred in C. pseudotuberculosis, we analysed the gene encoding DnaK and demonstrated that the expression of the gene was induced 2.5-fold relative to the control, which is consistent with the proposed mechanism (see Additional file 6: Table S3).
Distribution of genes among the simulated conditions
Among the shared repressed genes were those encoding proteins related mostly to energy metabolism, sugar transport, amino acids that contribute greatly to the maintenance and replication of the organism in the environment , such as argJ, nanK, opp (see Additional file 8: Table S5) and hypothetical proteins. These data are consistent with the colony-forming unit experiments, which demonstrated a decrease in replication in the stressful environments (Figure 1). A reduction in growth is a survival strategy and occurs in essentially all stressful situations because the organism’s ability to perceive the environment and to modulate the response-controlling mechanisms of resistance, metabolism and other processes  which are increased, specific and suitable for the new conditions, is essential.
Expression of genes encoding sigma factors under simulated conditions
In their natural environment, or in the host during the infection process, the bacterium is exposed to disturbances that require a fast and adaptive response to ensure the survival of the pathogen. Therefore, it is necessary for the bacterium to change the pattern of expression of the genes encoding proteins that directly combat the deleterious nature of the stress . Certain sigma factors play an important and fundamental role in regulating the expression of virulence genes ; therefore, it is important to analyse these genes individually.
In the C. pseudotuberculosis strain 1002 genome, 8 genes encoding sigma factors were identified, which included the essential sigma factor SigA, non-essential and alternative SigB and 6 alternative factors belonging to the group of extracytoplasmic factors SigC, SigD, SigE, SigH, SigK and SigM, which are dispensable and induced frequently in response to specific conditions, for example, in response to stress.
It is possible that sigA encodes sigma factor RpoD, alternatively named sigma 70, which promotes the binding of RNA polymerase to specific sites by activating the transcription of most genes essential for exponential growth in Escherichia coli (E. coli) . In C. pseudotuberculosis 1002, the gene contains the 4 domains conserved in the sigma 70 family. Sporadically, σA can act as an alternative sigma factor that is specifically required for the expression of virulence genes. An investigation of the role of sigA in E. coli in which the rpoD gene was induced (from a lack of amino acids and thermal shock), suggested that the protein was involved in the mechanism of recovery from stress . Streptomyces spp. contains several homologues of the principal factor that are not essential for replication under normal conditions but appear to play a role under certain growth conditions .
Because sigA demonstrated changes in expression between the control and the stress conditions in C. pseudotuberculosis 1002, it may play a role as an alternative sigma factor, regula  ting genes involved in the maintenance of the bacterium.
The gene encoding sigma M was considered differentially expressed only under osmotic stress. This gene may not be expressed preferentially in the early exponential phase. Microarray analysis of C. glutamicum in the exponential phase (OD610nm = 0.2 to 0.3) demonstrated that the disruption of sigM did not affect the level of transcription of genes induced by thermal shock, which implies that this sigma factor is not involved in the regulation of gene expression in response to this stress .
In M. tuberculosis, expression of sigM occurred in the stationary phase and only under conditions of thermal stress . In C. glutamicum, experiments revealed that the deletion of sigM caused a reduction in the number of viable cells under conditions of heat shock, cold shock and disulphide (an oxidative stress subtype) stress in the exponential growth phase (OD600nm = 0.7). Furthermore, experiments using real-time PCR demonstrated that the transcription of sigM increased significantly after the application of the stresses. These results suggest that sigma factor M is involved in the stress response , albeit with a strong indication that the involvement occurs at the late growth phase.
Because there was an increase in the expression of sigM under osmotic stress in C. pseudotuberculosis 1002, a study of the regulon will be required to identify the genes encoding proteins that assist in the persistence of the bacterium in this hostile environment.
Sigma B is a very common sigma factor in the stress response. The role of the protein in the response to acidity tolerance was established in Bacillus subtilis, Brevibacterium flavum, Listeria monocytogenes and Staphylococcus aureus[43–45]. In B. subtilis, SigB regulates the majority of stress responses, thereby contributing to the transcription of more than100 genes . Another study demonstrated a reformulation of transcription during the infection process and the relevance of sigma B in controlling the expression of virulence genes important for adaptation to the intestinal environment in L. monocytogenes. In L. monocytogenes, σB contributed to cell survival under different stress conditions , and the absence of the factor reduced the ability of the species to invade epithelial cells . It is known that σB contributes to virulence in several Gram-positive pathogens , and these observations underscore the need to study the factor separately in C. pseudotuberculosis and to identify the regulon that contributes to the survival and escape of the bacterium from the host immune system.
The sigE gene, encoding the presumed factor SigE, a member of the σ70 subfamily exhibiting extracytoplasmic function, regulates functions related to perception and response to changes in the periplasm and in the extracytoplasmic environment. In Haemophilus influenzae, the expression of rpoE increased 102-fold after phagocytosis by macrophages, and the survival of Haemophilus influenzae containing a mutated rpoE was reduced compared to wild type . According to a study in Vibrio cholerae, rpoE mutants attenuate virulence, and the ability of the bacteria to colonise the mouse intestines is reduced.
According to a study in C. pseudotuberculosis 1002 , compared with the wild type strain, the mutant 1002 strain is more sensitive to agents that generate nitrosative, acid (pH 5.5) and surface stresses, which indicates the important role SigE plays in the persistence of the bacterium in hostile environments.
The investigation of a mutant sigE in M. tuberculosis H37Rv  demonstrated that the mutant was more sensitive to various environmental stresses, such as thermal shock, SDS and oxidative agents but not to acidity. In a report on M. Smegmatis, the mutant sigE was more sensitive to acidity, hydrogen peroxide, SDS and thermal shock than the wild type. In the M. tuberculosis H37Rv study , the authors determined that sigE is required for the stress response and for the ability to grow and survive inside macrophages. Furthermore, it was suggested that sigE influences the level of sigB in the cell; however, the expression of sigB is not completely influenced by sigE because the level of sigB was reduced in the mutant strain, whereas the mRNA levels of other sigma factors were not affected.
In C. pseudotuberculosis 1002 under acidic conditions, sigE was induced and the expression of both sigE and sigB was significant, whereas under conditions where the mRNA level of sigE was not considered significant, the sigB expression was also not significant. These sigma factors may influence each other at this exponential phase; however, studies that are more specific will be necessary to verify this hypothesis.
The sigH gene, which encodes the extracytoplasmic sigma factor H, was induced (within the established cutoff value) only under the acidic conditions. In M. tuberculosis, the sigH factor is involved in the response to different stresses and it has been suggested that it regulates the expression of genes involved in the intracellular survival of the microorganism. Additionally, proteins that are part of the sigH facto regulon may interact with the host’s immune system, modulating its response . This observation suggests that at the beginning of the C. pseudotuberculosis 1002 replication process, the expression of the sigH gene and its regulon are required for the bacterium to persist in an acidic environment.
In C. glutamicum, a mutation in the sigH gene blocked the transcription of sigM, and the identification of the sigH promoter upstream of sigM implied that the factor is under the direct transcriptional control of sigH. It was determined that in C. pseudotuberculosis 1002, when sigH expression increased under conditions of stress, sigM expression also increased, although the sigM transcripts did not exhibit significant induction during this growth phase. These results suggest that the sigH gene responded to an environment similar to that encountered by the bacterium in the host, highlighting the importance of identifying the regulon and the genes responsible for the persistence of the bacterium in the environment.
Identification of non-coding RNA
Identification in silico of ncRNA in the genome of Corynebacterium pseudotuberculosis strain 1002 using the RFAM programme
In most cases, the ncRNAs in bacteria play a role in regulating cellular response during environmental change, which is beneficial to the organism, which needs to adapt quickly and efficiently to these changes and may be impaired by the cell because of the energy cost involved in the expression of a large number of genes. Compared to the synthesis of regulatory proteins, less energy is needed for the synthesis of a small RNA (sRNA) . Among the ncRNAs detected in the genome, riboswitches (elements that control expression in response to various metabolites and can function as sensors or binding sites for a receptor protein that senses a change in the cellular environment) were identified through the analysis of active transcripts. Furthermore, other elements that may be responsible for the synthesis of peptidoglycan, which promotes the expression of genes important under stress conditions [56, 57], were identified.
Next-generation technology allowed the identification of a large number of genes presumed to be required by C. pseudotuberculosis 1002 for survival in unfavourable environments, such as acidity, thermal shock and osmotic stress. A number of these genes encode hypothetical proteins, which highlight the need for further investigation of the microorganism. Among the most relevant biological processes identified under all simulated conditions were the processes of adhesion, stress response and oxidoreduction. In these processes, genes involved in the virulence of the organism were affected and should be investigated in more detail to identify the roles they play in the cell, especially inside the host. Furthermore, it is believed that the identification of these genes may contribute to the development of vaccines that are more effective, diagnostic kits and therapies for caseous lymphadenitis.
The expression of sigma factors varied under the different conditions and at the beginning of the exponential phase, these important factors involved in the regulation of genes required for the maintenance of microorganisms under different environmental conditions were observed. Understanding the regulon of the genes would clarify the biology of the organism.
The predicted ncRNAs were identified by the coverage of the transcripts and these factors presumably contribute to the regulation of genes related to the persistence of the bacterium in harmful environments. Further studies will be required to confirm the role of the ncRNAs in escaping the host immune system and contributing to the survival of the bacterium in hostile environments.
The data regarding the expression of induced or repressed genes do not necessarily indicate protein translation; therefore, future experiments will focus on understanding the biology of the transcripts and their products.
Culture conditions: obtaining bacterial cells
Corynebacterium pseudotuberculosis strain 1002 was grown in petri dishes containing BHI media (broth composed of (g/L): calf-brain infusion 200.00, beef-heart infusion 250.00, proteose peptone 10.00, dextrose 2.00, sodium chloride 5.00, di-sodium phosphate 2.50 pH 7.4 ± 0.2 at 25°C) at room temperature (RT). One colony was used to prepare the pre-inoculum in 20 mL of BHI media supplemented with 0.05% Tween 80. The culture was grown overnight at 37°C in a shaker at 160 rpm. One millilitre of this pre-inoculum was used to prepare the inoculum in an Erlenmeyer flask containing 100 mL fresh BHI, and this culture was incubated at 37°C at 160 rpm. This preparation was monitored until the beginning of the exponential growth phase (A600 = 0.2), which was reached approximately 2.5 hours after the initial inoculation (see Additional file 10: Figure S4).
Application of stresses
After the culture reached the beginning of the exponential growth phase, the inoculum was divided into 4 50-mL Falcon tubes (1 for each condition), each containing a final volume of 20 mL, and these tubes were then centrifuged for 3 minutes at 8,000 rpm at RT. The pellet was resuspended in fresh BHI specific to each condition. For the acid stress condition, the media was supplemented with hydrochloric acid (which the pH changed to 5). Osmotic stress was achieved with 2 M NaCl, and thermal stress was induced by resuspending the pellet in BHI medium pre-heated to 50°C. In the control condition, bacterial pellets were resuspended in BHI medium at a physiological condition. After the addition of culture media, the tubes were kept in a shaker at 37°C and 160 rpm for 15 minutes, with the exception of the thermal stress sample that was subjected to a temperature of 50°C. An aliquot of each condition was used for decimal dilutions from 10-1 to 10-6, from which 10-4 to 10-6 bacteria were seeded in BHI agar, and petri dishes were kept at 37ºC for 48 hours for viability analysis and colony counting (this step was performed in duplicate). The remaining sample was subjected to centrifugation at RT for 3 minutes at 8,000 rpm, and the pellet was resuspended in 2 ml of RNAlater, according to the manufacturer’s instructions.
The bacteria suspended in RNAlater® buffer were subjected to total RNA extraction using the ChargeSwitch® total RNA cell kit (Invitrogen, USA) in accordance with the manufacturer’s recommendations, including the following adaptations: after the addition of the lysis buffer (Invitrogen), the material was transferred to 2-mL tubes partially filled with 1-mm diameter glass microbeads (Bertin Technologies). The cells were lysed mechanically using a Prescellys 24 homogeniser, set at 6,500 rpm, for 2 cycles (15 seconds per cycle) with an interval of 30 seconds between the cycles. The samples were centrifuged for 1 minute, and the supernatant transferred to fresh 2-ml tubes and incubated in a dry bath at 60°C for 15 minutes (represents the complete original protocol). DNase was added to eliminate the residual genomic DNA. The elution of the total RNA from the magnetic beads was performed using 100 μL of milli-Q RNase-free water. The amount of total RNA was assessed using a Qubit® 2.0 fluorometer (Invitrogen).
mRNA enrichment through rRNA depletion
To enrich the mRNA, rRNA from each total RNA sample was removed using the Ribominus™ Transcriptome Isolation kit for yeast and bacteria (Invitrogen, USA), in accordance with the manufacturer’s recommendations. The rRNA-depleted RNA was used for cDNA synthesis using the SOLiD™ Total RNA-Seq kit in accordance with the standard protocol recommended by the manufacturer, and the material was quantified in a Qubit® 2.0 fluorometer (Invitrogen).
Sequencing in SOLiD™
The depleted RNA was fragmented using RNase III in preparation for amplification of the cDNA library, which was produced by reverse transcription from adapters attached to the ends of the RNA molecules, in accordance with the SOLiD™ Total RNA-Seq kit protocol (Life Technologies™, CA). Next, 6% denaturing polyacrylamide gel electrophoresis was performed and fragments of appropriate sizes (150 to 250 bases) were cut from the gel for cDNA amplification using PCR. Following recommended protocols, the cDNA was purified and the sizes were confirmed using 2% agarose electrophoresis. The PCR amplification in emulsion was performed using primers complementary to the adapters, in accordance with the Applied Biosystems SOLiD™ 3 Plus System Templated Bead Preparation Guide. After amplification, the microspheres were deposited onto slides for sequencing in accordance with the manufacturer’s recommendations. The SOLiD™ 3 Plus system was used to sequence the 50-nucleotide RNA reads.
Analysis in silico
After obtaining the reading files using the SOLiD™ technique, the data were loaded into Bioscope version 1.2.1-5 programme (Life Technologies™, CA), and the gene expression data were obtained and quantified as reads per kilobase of coding sequence per million reads (RPKM) . The DEGseq programme  was used to identify differentially expressed genes. The programme uses the output file from Bioscope and the RPKM values as the input file. For differentially expressed genes, a cutoff value of p <0.001 was defined. The data from each stimulon were analysed using Blast2GO (http://www.blast2go.com), and the results were exported to CoreStImulon programme  to search more rapidly for the genes comprising each biological process defined by Blast2GO.
To predict the non-coding RNA, a similarity search was performed in the RFAM database  using the script rfam_scan-1.0.3.pl (http://rfam.sanger.ac.uk/), and sRNAs smaller than 70 bp were discarded to minimise the number of false positives.
Transcriptomic coverage for all sRNAs annotated by RFAM was confirmed by manual curation in the program Artemis using the file BAM.
high-throughput sequencing of cDNA libraries
Reads per kilobase of coding sequence per million mapped
peptide methionine sulphoxide reductase
complementary DNA synthesised from RNA
- Sig and σ:
Brain heart infusion broth.
The authors thank the Genome and Proteome Network of the State of Pará [Rede Paraense de Genômica e Proteômica] for supporting the development of this work and the Research Support Foundation of the State of Pará [Fundação de Amparo a Pesquisa do Estado do Pará] for funding the study. The authors also thank the National Council for Scientific and Technological Development [Conselho Nacional de Desenvolvimento Científico e Tecnológico – CNPq] and the Coordination for Improvement of Higher Education Personnel [Coordenação de Aperfeiçoamento de Pessoal de Nível Superior, CAPES] for the financial support provided to the authors.
- Dorella FA, et al: Review article Corynebacterium pseudotuberculosis : microbiology , biochemical properties , pathogenesis and molecular studies of virulence. Vet Res. 2006, 37: 201-218. 10.1051/vetres:2005056.PubMedView ArticleGoogle Scholar
- Kazmierczak MJ, Wiedmann M, Boor KJ: Alternative Sigma Factors and Their Roles in Bacterial Virulence. Society. 2005, 69: 527-543.Google Scholar
- Jin Y, Tian Y, Zhang W, Jang SH, Jen AK, Meldrum DR: Tracking bacterial infection of macrophages using a novel red-emission pH sensor. Anal Bioanal Chem. 2010, 398 (3): 1375-1384. 10.1007/s00216-010-4060-6.PubMed CentralPubMedView ArticleGoogle Scholar
- Ramos JL, Gallegos MT, Marqués S, Ramos-González MI, Espinosa-Urgel M, Segura a: Responses of Gram-negative bacteria to certain environmental stressors. Curr Opin Microbiol. 2001, 4: 166-171. 10.1016/S1369-5274(00)00183-1.PubMedView ArticleGoogle Scholar
- Manganelli R, Dubnau E, Tyagi S, Kramer FR, Smith I: Differential expression of 10 sigma factor genes in Mycobacterium tuberculosis. Mol Microbiol. 1999, 31 (2): 715-724. 10.1046/j.1365-2958.1999.01212.x.PubMedView ArticleGoogle Scholar
- Wu Q-L, Kong D, Lam K, Husson RN: A mycobacterial extracytoplasmic function sigma factor involved in survival following stress. J Bacteriol. 1997, 179 (9): 2922-2929.PubMed CentralPubMedGoogle Scholar
- Missiakas D, Raina S: The extracytoplasmic function sigma factors: role and regulation. Mol Microbiol. 1998, 28 (6): 1059-1066. 10.1046/j.1365-2958.1998.00865.x.PubMedView ArticleGoogle Scholar
- Tashjian J, Campbell S: Interaction between caprine macrophages and Corynebacterium pseudotuberculosis: An electron microscopic study. Am J Vet Res. 1983, 44 (4): 690-693.PubMedGoogle Scholar
- Billington SJ, Esmay PA, Songer JG, Jost BH: Identification and role in virulence of putative iron acquisition genes from Corynebacterium pseudotuberculosis. FEMS Microbiol Lett. 2002, 208: 41-45. 10.1111/j.1574-6968.2002.tb11058.x.PubMedView ArticleGoogle Scholar
- Ca M, Gyles CL: Characterization of strains of corynebacterium pseudotuberculosis. Can J Comp Med. 1982, 46: 206-208.Google Scholar
- Wilson MJ, Brandon MR, Walker J: Molecular and biochemical characterization of a protective 40-kilodalton antigen from Corynebacterium pseudotuberculosis. Infect Immun. 1995, 63: 206-211.PubMed CentralPubMedGoogle Scholar
- Sorek R, Cossart P: Prokaryotic transcriptomics: a new view on regulation, physiology and pathogenicity. Nat Rev Genet. 2010, 11: 9-16.PubMedView ArticleGoogle Scholar
- Sendler E, Johnson GD, Krawetz SA: Local and global factors affecting RNA sequencing analysis. Anal Biochem. 2011, 419 (2): 317-322. 10.1016/j.ab.2011.08.013.PubMedView ArticleGoogle Scholar
- Jozefczuk S, Klie S, Catchpole G, Szymanski J, Cuadros-Inostroza A, Steinhauser D, Selbig J, Willmitzer L: Metabolomic and transcriptomic stress response of Escherichia coli. Mol Syst Biol. 2010, 6: 364-PubMed CentralPubMedView ArticleGoogle Scholar
- Mortazavi A, Williams BA, Mccue K, Schaeffer L, Wold B: Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008, 5: 7-8.View ArticleGoogle Scholar
- Wang L, Feng Z, Wang X, Wang X, Zhang X: DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010, 26 (1): 136-138. 10.1093/bioinformatics/btp612.PubMedView ArticleGoogle Scholar
- Castro TL: Ion Torrent-based transcriptional assessment of a Corynebacterium pseudotuberculosis equi strain reveals denaturing high-performance liquid chromatography a promising rRNA depletion method. Microb Biotechnol. 2013, 6 (2): 1-10.Google Scholar
- Pinto AC, Ramos RT, Silva WM, Rocha FS, Barbosa S, Miyoshi A, Schneider MP, Silva A, Azevedo V: The core stimulon of Corynebacterium pseudotuberculosis strain 1002 identified using ab initio methodologies. Integr Biol. 2012, 4 (7): 789-794. 10.1039/c2ib00092j.View ArticleGoogle Scholar
- Farr S, Kogoma T: Oxidative stress responses in Escherichia coli and Salmonella typhimurium. Microbiol Rev. 1991, 55 (4): 561-585.PubMed CentralPubMedGoogle Scholar
- Sleator RD, Hill C: Bacterial osmoadaptation: the role of osmolytes in bacterial stress and virulence. FEMS Microbiol Rev. 2002, 26 (1): 49-71. 10.1111/j.1574-6976.2002.tb00598.x.PubMedView ArticleGoogle Scholar
- Ratledge C, Dover LG: Iron metabolism in pathogenic bacteria. Annu Rev Microbiol. 2000, 54: 881-941. 10.1146/annurev.micro.54.1.881.PubMedView ArticleGoogle Scholar
- Sasindran SJ, Saikolappan S, Dhandayuthapani S: Methionine sulfoxide reductases and virulence of bacterial pathogens. Future Microbiol. 2007, 2: 619-630. 10.2217/174609220.127.116.119.PubMedView ArticleGoogle Scholar
- Delaye L, Becerra A, Orgel L, Lazcano A: Molecular Evolution of Peptide Methionine Sulfoxide Reductases (MsrA and MsrB): On the Early Development of a Mechanism That Protects Against Oxidative Damage. J Mol Evol. 2007, 64 (1): 15-32. 10.1007/s00239-005-0281-2.PubMedView ArticleGoogle Scholar
- Dhandayuthapani S, Jagannath C, Nino C, Saikolappan S, Sasindran SJ: Methionine sulfoxide reductase B (MsrB) of Mycobacterium smegmatis plays a limited role in resisting oxidative stress. Tuberculosis (Edinb). 2009, 89 (Suppl 1): S26-32.View ArticleGoogle Scholar
- Alamuri P, Maier RJ: Methionine sulphoxide reductase is an important antioxidant enzyme in the gastric pathogen Helicobacter pylori. Mol Microbiol. 2004, 53: 1397-1406. 10.1111/j.1365-2958.2004.04190.x.PubMedView ArticleGoogle Scholar
- Rajan LA, Joseph TC, Thampuran N, James R: Functional Characterization and Sequence Analysis of Choline Dehydrogenase from Escherichia coli. Genet Eng Biotechnol J. 2010Google Scholar
- Das Gupta T, Bandyopadhyay B, Das Gupta SK: Modulation of DNA-binding activity of Mycobacterium tuberculosis HspR by chaperones. Microbiology. 2008, 154: 484-490. 10.1099/mic.0.2007/012294-0.PubMedView ArticleGoogle Scholar
- Stewart GR, Wernisch L, Stabler R, Mangan JA, Hinds J, Laing KG, Young DB, Butcher PD: Dissection of the heat-shock response in Mycobacterium tuberculosis using mutants and microarrays. Microbiology. 2002, 148 (10): 3129-3138.PubMedView ArticleGoogle Scholar
- Bandyopadhyay B, Gupta TD, Roy D, Gupta SKD: DnaK Dependence of the Mycobacterial Stress-Responsive Regulator HspR Is Mediated through Its Hydrophobic C-Terminal Tail. J Bacteriol. 2012, 194 (17): 4688-4697. 10.1128/JB.00415-12.PubMed CentralPubMedView ArticleGoogle Scholar
- Isabella VM, Clark VL: Deep sequencing-based analysis of the anaerobic stimulon in Neisseria gonorrhoeae. BMC Genomics. 2011, 12: 51-10.1186/1471-2164-12-51.PubMed CentralPubMedView ArticleGoogle Scholar
- Weerasinghe JP, Dong T, Schertzberg MR, Kirchhof MG, Sun Y, Schellhorn HE: Stationary phase expression of the arginine biosynthetic operon argCBH in Escherichia coli. BMC Microbiol. 2006, 6: 14-10.1186/1471-2180-6-14.PubMed CentralPubMedView ArticleGoogle Scholar
- Weihofen W, Berger M, Chen H, Saenger W, Hinderlich S: Structures of human N-Acetylglucosamine kinase in two complexes with N-Acetylglucosamine and with ADP/glucose: insights into substrate specificity and regulation. J Mol Biol. 2006, 364: 388-399. 10.1016/j.jmb.2006.08.085.PubMedView ArticleGoogle Scholar
- Goodell EW, Higgins CF: Uptake of cell wall peptides by Salmonella typhimurium and Escherichia coli. J Bacteriol. 1987, 169 (8): 3861-3865.PubMed CentralPubMedGoogle Scholar
- Rohde KH, Veiga DF, Caldwell S, Balázsi G, Russell DG: Linking the transcriptional profiles and the physiological states of Mycobacterium tuberculosis during an extended intracellular infection. PLoS Pathog. 2012, 8 (6): e1002769-10.1371/journal.ppat.1002769.PubMed CentralPubMedView ArticleGoogle Scholar
- Hecker M, Völker U: General stress response of Bacillus subtilis and other bacteria. Adv Microb Physiol. 2001, 44: 35-91.PubMedView ArticleGoogle Scholar
- Gomez M, Doukhan L, Nair G, Smith I: sigA is an essential gene in Mycobacterium smegmatis. Mol Microbiol. 1998, 29: 617-628. 10.1046/j.1365-2958.1998.00960.x.PubMedView ArticleGoogle Scholar
- Helmann JD, Chamberlin MJ: Structure and function of bacterial sigma factors. Annu Rev Biochem. 1988, 57 (1): 839-872. 10.1146/annurev.bi.57.070188.004203.PubMedView ArticleGoogle Scholar
- Taylor WE, Straus DB, Grossman AD, Burton ZF, Gross CA, Burgess RR: Transcription from a heat-inducible promoter causes heat shock regulation of the sigma subunit of E. coli RNA polymerase. Cell. 1984, 38 (2): 371-381. 10.1016/0092-8674(84)90492-6.PubMedView ArticleGoogle Scholar
- Buttner M, Chater K, Bibb M: Cloning, disruption, and transcriptional analysis of three RNA polymerase sigma factor genes of Streptomyces coelicolor A3 (2). J Bacteriol. 1990, 172 (6): 3367-3378.PubMed CentralPubMedGoogle Scholar
- Ehira S, Teramoto H, Inui M, Yukawa H: Regulation of Corynebacterium glutamicum heat shock response by the extracytoplasmic-function sigma factor SigH and transcriptional regulators HspR and HrcA. J Bacteriol. 2009, 191: 2964-2972. 10.1128/JB.00112-09.PubMed CentralPubMedView ArticleGoogle Scholar
- Agarwal N, Woolwine SC, Tyagi S, Bishai WR: Characterization of the Mycobacterium tuberculosis sigma factor SigM by assessment of virulence and identification of SigM-dependent genes. Infect Immun. 2007, 75 (1): 452-461. 10.1128/IAI.01395-06.PubMed CentralPubMedView ArticleGoogle Scholar
- Nakunst D, Larisch C, Hüser AT, Tauch A, Pühler A, Kalinowski J: The extracytoplasmic function-type sigma factor SigM of Corynebacterium glutamicum ATCC 13032 is involved in transcription of disulfide stress-related genes. J Bacteriol. 2007, 189: 4696-4707. 10.1128/JB.00382-07.PubMed CentralPubMedView ArticleGoogle Scholar
- Kovács T, Hargitai A, Kovács KL, Mécs I: pH‐dependent activation of the alternative transcriptional factor σB in Bacillus subtilis. FEMS Microbiol Lett. 1998, 165 (2): 323-328.PubMedGoogle Scholar
- Bischoff M, Entenza J, Giachino P: Influence of a Functional sigB Operon on the Global Regulators sar and agr inStaphylococcus aureus. J Bacteriol. 2001, 183 (17): 5171-5179. 10.1128/JB.183.17.5171-5179.2001.PubMed CentralPubMedView ArticleGoogle Scholar
- Ferreira A, Sue D, O'Byrne CP, Boor KJ: Role of Listeria monocytogenes σB in survival of lethal acidic conditions and in the acquired acid tolerance response. Appl Environ Microbiol. 2003, 69 (5): 2692-2698. 10.1128/AEM.69.5.2692-2698.2003.PubMed CentralPubMedView ArticleGoogle Scholar
- Vijay K, Brody MS, Fredlund E, Price CW: A PP2C phosphatase containing a PAS domain is required to convey signals of energy stress to the σB transcription factor of Bacillus subtilis. Mol Microbiol. 2000, 35 (1): 180-188. 10.1046/j.1365-2958.2000.01697.x.PubMedView ArticleGoogle Scholar
- Toledo-Arana A, Dussurget O, Nikitas G, Sesto N, Guet-Revillet H, Balestrino D, Loh E, Gripenland J, Tiensuu T, Vaitkevicius K, et al: The Listeria transcriptional landscape from saprophytism to virulence. Nature. 2009, 459: 950-956. 10.1038/nature08080.PubMedView ArticleGoogle Scholar
- Ferreira A, O'Byrne CP, Boor KJ: Role of sigB in Heat, Ethanol, Acid, and Oxidative Stress Resistance and during Carbon Starvation in Listeria monocytogenes. Appl Environ Microbiol. 2001, 67 (10): 4454-4457. 10.1128/AEM.67.10.4454-4457.2001.PubMed CentralPubMedView ArticleGoogle Scholar
- Garner M, Njaa B, Wiedmann M, Boor K: Sigma B contributes to Listeria monocytogenes gastrointestinal infection but not to systemic spread in the guinea pig infection model. Infect Immun. 2006, 74 (2): 876-886. 10.1128/IAI.74.2.876-886.2006.PubMed CentralPubMedView ArticleGoogle Scholar
- Craig JE, Nobbs A, High NJ: The Extracytoplasmic Sigma Factor, sigE, Is Required for Intracellular Survival of Nontypeable Haemophilus influenzae in J774 Macrophages. Infect Immun. 2002, 70 (2): 708-715. 10.1128/IAI.70.2.708-715.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Kovacikova G, Skorupski K: The alternative sigma factor σE plays an important role in intestinal survival and virulence in Vibrio cholerae. Infect Immun. 2002, 70 (10): 5355-5362. 10.1128/IAI.70.10.5355-5362.2002.PubMed CentralPubMedView ArticleGoogle Scholar
- Pacheco LG, Castro TL, Carvalho RD, Moraes PM, Dorella FA, Carvalho NB, Slade SE, Scrivens JH, Feelisch M, Meyer R, et al: A Role for Sigma Factor sigma(E) in Corynebacterium pseudotuberculosis Resistance to Nitric Oxide/Peroxide Stress. Front Microbiol. 2012, 3: 126-PubMed CentralPubMedView ArticleGoogle Scholar
- Manganelli R, Voskuil MI, Schoolnik GK, Smith I: The Mycobacterium tuberculosis ECF sigma factor sigmaE: role in global gene expression and survival in macrophages. Mol Microbiol. 2001, 41: 423-437. 10.1046/j.1365-2958.2001.02525.x.PubMedView ArticleGoogle Scholar
- Dutta NK, Mehra S, Martinez AN, Alvarez X, Renner NA, Morici LA, Pahar B, Maclean AG, Lackner AA, Kaushal D: The stress-response factor SigH modulates the interaction between Mycobacterium tuberculosis and host phagocytes. PloS One. 2012, 7 (1): e28958-10.1371/journal.pone.0028958.PubMed CentralPubMedView ArticleGoogle Scholar
- Mikulík K, Palečková P, Felsberg J, Bobek J, Zídková J, Halada P: SsrA genes of streptomycetes and association of proteins to the tmRNA during development and cellular differentiation. Proteomics. 2008, 8 (7): 1429-1441. 10.1002/pmic.200700560.PubMedView ArticleGoogle Scholar
- Mehta NB, et al: Riboswitches: classification, function and insilico approach. Int J. 2010, 1: 409-420.Google Scholar
- Weinberg Z, Wang JX, Bogue J, Yang J, Corbino K, Moy RH, Breaker RR: Comparative genomics reveals 104 candidate structured RNAs from bacteria, archaea, and their metagenomes. Genome Biol. 2010, 11 (3): R31-10.1186/gb-2010-11-3-r31.PubMed CentralPubMedView ArticleGoogle Scholar
- Bessho Y, Shibata R, Sekine S-i, Murayama K, Higashijima K, Hori-Takemoto C, Shirouzu M, Kuramitsu S, Yokoyama S: Structural basis for functional mimicry of long-variable-arm tRNA by transfer-messenger RNA. Proc Natl Acad Sci. 2007, 104 (20): 8293-8298. 10.1073/pnas.0700402104.PubMed CentralPubMedView ArticleGoogle Scholar
- Wassarman KM: Small RNAs in bacteria: diverse regulators of gene expression in response to environmental changes. Cell. 2002, 109 (2): 141-144. 10.1016/S0092-8674(02)00717-1.PubMedView ArticleGoogle Scholar
- Gardner PP, Daub J, Tate J, Moore BL, Osuch IH, Griffiths-Jones S, Finn RD, Nawrocki EP, Kolbe DL, Eddy SR: Rfam: Wikipedia, clans and the “decimal” release. Nucleic Acids Res. 2011, 39 (suppl 1): D141-D145.PubMed CentralPubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.