Previous Article | Next Article ![]()
Journal of Bacteriology, July 2006, p. 5228-5239, Vol. 188, No. 14
0021-9193/06/$08.00+0 doi:10.1128/JB.00507-06
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Unité des Rickettsies, UMR 6020 CNRS, IFR48, Faculté de Médecine, 27, Boulevard Jean Moulin, Marseille, France,1 Institut de Pharmacologie Moléculaire et Cellulaire, UMR 6097 CNRS/UNSA, Sophia Antipolis, France2
Received 10 April 2006/ Accepted 24 April 2006
|
|
|---|
|
|
|---|
T. whipplei is a gram-positive bacterium which for a long time was thought to be resistant to cultivation, but it was isolated in eukaryotic cells and propagated in culture at 37°C in 2000 (45). This has enabled further characterization of this pathogen, including genome sequencing (6, 44). With a 0.92-Mb genome, T. whipplei, which is located phylogenetically in the high-G+C-content gram-positive bacterial group between the genus Cellulomonas and the actinomycete clade (36), is a reduced-genome species belonging to the Actinobacteria with lower G+C contents. Like other reduced-genome bacteria, T. whipplei does not have several genes that regulate transcription and are classified in the K functional category in the Clusters of Orthologous Groups database. Based on the intracellular nature of these microorganisms, it was hypothesized that this deficiency was due to the rather stable environment inside host cells, which makes extensive gene regulation useless (34). To date, very little is known about the natural habitats of T. whipplei and the way that it infects hosts. Therefore, while T. whipplei does not contain most regulatory elements, this microorganism is thought to be an environmental agent (37) that is able to adapt to a wide variety of stress conditions.
Temperature change is the most common stress that all living organisms encounter in natural habitats. To overcome critical situations which could be generated by extreme temperatures, bacteria have evolved complex and specific mechanisms that are called cold shock and heat shock responses (43, 57). To deal with heat stress, bacteria overexpress heat shock proteins (HSP), including chaperones encoded by the dnaK and groE operons and ATP-dependent proteases (Clp and Lon). Chaperones prevent misfolding and aggregation of partially denatured proteins, while proteases degrade these proteins (21). Expression of HSP is regulated at the transcriptional level by various positive and negative mechanisms, depending on the bacterium. Positive control is provided by alternative sigma factors that target RNA polymerase to the heat shock gene promoter (22). In gram-positive bacteria, HSP are regulated mainly by less sophisticated mechanisms, called repressor mechanisms, which arose early in evolution (40). Such regulatory elements correspond to specific DNA sequences located in promoter regions, providing a simple and economical regulation pathway. An inverted repeat having a conserved consensus sequence (TTGGCGCTC-N9-GAGTGCTAA) and termed CIRCE (controlling inverted repeat of chaperone expression) (63) has been found in several bacteria in association with the groESL or dnaK operon. Another DNA sequence, the HspR-associated inverted repeat (HAIR motif), was described as a regulon under the control of HspR (50). This regulon has been found mainly in high-G+C-content gram-positive bacteria, particularly actinomycetes (52). HSP are classified according to their regulation characteristics, which differ for each bacterial species. In gram-positive bacteria, four regulatory classes have been defined. The class I genes include the genes encoding the classical chaperones DnaK, GroES, and GroEL that are negatively controlled by the hrcA-encoded protein which binds to the CIRCE sequence (26). The class II genes include genes encoding more than 40 stress proteins regulated by a
B factor. Class III gene expression, under negative control of the class III stress gene repressor (ctrR), recognizes a repeat sequence in the promoter. Class IV comprises other genes whose regulation is not known (9). Overlap between class I and class III HSP has been observed for streptococci (9).
The molecular basis of the cold shock response is also well documented. Bacteria respond to low temperatures by induction of a set of proteins that include the cold shock protein (CSP) family, which is autoregulated by an unusual long 5' untranslated region in the mRNA transcript (43). The high level of conservation of CSP in different bacterial lineages suggests that such proteins are very ancient. Moreover, the presence of cold shock domains in both eubacteria and eukaryotes indicates that cold shock domain-like proteins were present before the divergence of bacteria and eukarya/archaea (19). While this ancient family of small proteins is universally conserved among bacteria, a few exceptions have been observed in archaea (Methanococcus jannaschii, Methanococcus thermoautotrophicum, and Archaeoglobus fulgidus) and some bacteria, such as Synechocystis sp. strain PCC6803 and Borrelia burgdorferi (19). T. whipplei also lacks CSP (6, 44). It has been hypothesized that in cyanobacteria CSP are replaced by small cold-inducible proteins (Rbps) homologous to the RNA-binding domain found in eukaryotic RNA-binding proteins (48), but no rbp-like sequence was found in T. whipplei.
Here, the adaptive responses and the regulation mechanisms of T. whipplei exposed to various sudden temperature shifts were investigated by using a specific microarray. This work was the first attempt to analyze the transcriptome regulation of this pathogenic bacterium whose origin is unknown.
|
|
|---|
The following criteria were used to identify optimal forward and reverse primers for generating PCR products specific for each of the selected ORFs: (i) the primers were highly specific for the T. whipplei genome; (ii) the amplicons were roughly 260 bp long; (iii) the amplicons were highly specific for the corresponding genes; (iv) each primer contained 20 to 25 bases, and the annealing temperature of the primers was
60°C; and (v) each primer contained an additional sequence in the 5' area specific for forward and reverse primers. The universal sequences added were used for a second round of amplification.
ORF-specific fragments were amplified using T. whipplei DNA as the template under following conditions: 40 cycles of 20 s of denaturation at 95°C, 30 s of annealing at 60°C, and 30 s of extension at 72°C, with an initial 5 min of denaturation at 95°C and final extension at 72°C for 5 min. Following amplification, PCR products were visualized by ethidium bromide staining after migration in a 2% agarose gel. Purified amplicons (QIAquick 96-well purification kit; QIAGEN, Courtaboeuf, France) were sequenced (3100 genetic analyzer; Applied Biosystems, Courtaboeuf, France) to check the identity of amplified samples and were processed by performing a second PCR with universal primers.
Amplified DNA fragments resuspended in a 50% dimethyl sulfoxide solution were printed in quadruplicate on Nexterion Slide A+ (Schott AG, Mainz, Germany) by using an SDDC-2 (Bio-Rad, Marnes-la-Coquette, France). A set of 23 artificial genes that serve as analytical controls to validate, filter, and normalize microarray data contained in a Universal ScoreCard kit (Amersham Biosciences, Orsay, France) was also printed on glass slides. The spotted slides were cross-linked by using a UV Stratalinker (Stratagene, La Jolla, CA) with a total energy of 250 mJ and then were stored in a dry environment at room temperature.
Strain, medium, and growth conditions. All experiments were performed with mid-log cultures of T. whipplei strain Twist (35) grown at 37°C in Dulbecco modified Eagle medium/F12 medium (Invitrogen Life Technologies, Carlsbad, CA) supplemented with 10% fetal calf serum, 1% L-glutamine, and 1% human nonessential amino acids (Invitrogen). Bacterial growth was monitored by flow cytometry counting using a Microcyte portable flow cytometer (Optoflow AS, Oslo, Norway) and by quantitative PCR as previously described (15, 46).
RNA purification. For the heat shock experiments, flasks containing bacteria were transferred into a water bath maintained at 43°C for the following times: 15, 30, and 60 min. For cold shock assays, a bacterial suspension was incubated at either 4°C or 28°C for 1 or 6 h. For all experimental conditions tested, the basal level of transcripts was estimated by using RNA extracted from different culture flasks collected separately before incubation under the temperature stress conditions. Considering both the culture phase and the low growth rate of this bacterium (28 h) (46), we assumed that the effects observed were due to the thermal stress. RNA was extracted by brief sonication of bacteria in 1 ml of Trizol reagent (Invitrogen Life Technologies) as described previously (11). Purified RNA samples resuspended in diethyl pyrocarbonate-treated water were treated with RNase-free DNase Set (QIAGEN) and cleaned on an RNeasy column (QIAGEN) to remove DNA contamination. The amount and quality of each RNA sample were checked by automated capillary gel electrophoresis using a Bioanalyzer 2100 with RNA Nano LabChips (Agilent, Palo Alto, CA). Four hundred milliliters of a T. whipplei axenic culture was required to obtain 10 µg of purified RNA.
RNA labeling. Fluorescently labeled cDNA was prepared using a CyScribe first-strand cDNA labeling kit (Amersham Biosciences). Briefly, 10 µg of total RNA supplemented with control RNAs (Universal ScoreCard kit; Amersham Biosciences) was annealed with random nonamer primers. Total RNA was directly reverse transcribed and labeled by incorporation of Cy3-dCTP or Cy5-dCTP (Amersham Biosciences) using CyScript reverse transcriptase. The remaining RNA was then degraded by alkaline hydrolysis treatment, and after neutralization, labeled first-strand cDNA was purified with a CyScribe GFX purification kit (Amersham Biosciences). Before hybridization, the levels of Cy3 and Cy5 incorporation were quantified by measuring the absorbance at 550 and 650 nm, respectively. Hybridizations were performed with incorporation levels ranging from 80 to 200 pmol of fluorochromes per µg of cDNA.
Hybridization and washing. Microarrays that were washed in 0.1% sodium dodecyl sulfate (SDS) (1 min) and rinsed with H2O (1 min) were boiled for 3 min in H2O in order to denature spotted double-strand DNA. Following prehybridization at 42°C for 45 min in 5x SSC-0.1% SDS-0.1% bovine serum albumin, the microarrays were rinsed in H2O (1 min) and dried with compressed nitrogen (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate). Hybridization was then carried out using two samples of cDNA (10 µg each) that were labeled with Cy3- or Cy5-dCTP, pooled, evaporated (SpeedVac concentrator), and resuspended in 6 µl of nuclease-free H2O. The mixture was heated at 95°C for at least 2 min and then cooled on ice for 30 s before addition of 7.5 µl of microarray hybridization buffer (supplied with the CyScribe kit) and 15 µl of 100% (vol/vol) formamide and applied to a microarray slide with a glass coverslip (24 by 60 mm). Following 18 h of hybridization at 42°C, the microarrays were washed in 2x SSC-0.2% SDS for 10 min, in 2x SSC, and then in 0.2x SSC. Finally, the microarrays were dried with compressed nitrogen and scanned at a resolution of 10 µm by using a ScanArray Express (Perkin-Elmer, Boston, MA).
Analysis of microarray data. All microarray results have been deposited in the GEO database (http://www.ncbi.nlm.nih.gov/geo/) (4) under GEO Series accession number GSE3693. Signal intensity and local background measurements were obtained for each spot by using the ScanArray Express microarray analysis system software (version 2.1.8; Perkin-Elmer). The data filtering and normalization results were then processed with the Microsoft Excel software. Spots with median background-corrected signal intensities in both channels that were less than twice the median background intensity were not included in further analyses. Data normalization was performed for the remaining spots by using total intensity normalization methods. The normalized test signal/reference signal log ratio for each spot was recorded. The data were processed by using the TMEV software (www.tigr.org/software/TM4/) (47). An analysis of variance test was applied to the data, and genes with a P value of <0.02 were considered to have significant differential expression. Significant changes in gene expression were identified with significance analysis microarrays (56) using two-class paired data and a 1.5-fold cutoff. All experiments were conducted three times, including dye swapping, which yielded 12 measurements per gene (representing four technical and three biological replicates). Gene expression was determined by determining the mean of the 12 values obtained.
Real-time RT-PCR. Reverse transcription (RT) from 2 µg of total RNA with 1 µg of random hexamer primers was carried out to generate cDNA. A real-time PCR was then performed for each cDNA preparation (1:20 dilution), using the LightCycler system (Roche) together with the SYBR Green master mixture. Eight genes exhibiting either up- or down-regulation (as shown by microarray data) were selected, as were invariant targets (namely, leuS, mgt, and TWT639). The levels of expression of these three genes, which were not sensitive to temperature changes, were used for normalization of real-time RT-PCR values. For each primer pair, a standard curve was constructed with genomic DNA of T. whipplei as the template. The relative expression ratios of target genes were calculated using the Pfaffl model (42). This model incorporates the amplification efficiencies of the target and reference (normalization) genes to correct for differences between the two assays.
Identification of potential DNA binding sites. Inverted-repeat elements, namely CIRCE (63) and HAIR (18) motifs, were searched using the Genome2D software (3, 63) with the whole T. whipplei genome (44).
|
|
|---|
![]() View larger version (36K): [in a new window] |
FIG. 1. Relative normalized fluorescence intensities of DNA microarrays. (A) Comparison of cDNAs derived from the same culture of bacteria grown at 37°C. (B) Comparison of cDNAs derived from two independent cultures of bacteria grown at 37°C and at 4°C. (C) Comparison of cDNAs derived from two independent cultures of bacteria grown at 37°C and at 28°C. (D) Comparison of cDNAs derived from two independent cultures of bacteria grown at 37°C and at 43°C. The upper and lower dashed lines indicate 1.5-fold changes in the signal intensities.
|
|
View this table: [in a new window] |
TABLE 1. Genes with significant differences (P < 0.02 and 1.5-fold change) between expression at 37°C and expression at 4°C, 28°C, and 43°C
|
![]() View larger version (21K): [in a new window] |
FIG. 2. Comparison of transcription measurements obtained by microarray and real-time RT-PCR assays. The relative transcriptional levels for 11 genes listed in Table 2 were determined by microarray analysis and real-time RT-PCR. The real-time RT-PCR log2 values were plotted against the microarray log2 values.
|
|
View this table: [in a new window] |
TABLE 2. Oligonucleotide primers used for real-time RT-PCR
|
![]() View larger version (27K): [in a new window] |
FIG. 3. Genes differentially expressed at 4°C and 43°C compared to 37°C grouped by functional classification according to The Institute for Genome Research T. whipplei genome database (http://www.tigr.org/). The percentages of genes were calculated from the number of genes belonging to each functional category, as follows: column 1, amino acid biosynthesis; column 2, biosynthesis of cofactors, prosthetic groups, and carriers; column 3, cell envelope; column 4, cellular processes; column 5, central intermediary metabolism; column 6, DNA metabolism; column 7, energy metabolism; column 8, fatty acid and phospholipid metabolism; column 9, protein fate; column 10, protein synthesis; column 11, purines, pyrimidines, nucleosides, and nucleotides; column 12, regulatory functions; column 13, transcription; column 14, transport and binding proteins; column 15, unknown function.
|
![]() View larger version (20K): [in a new window] |
FIG. 4. Kinetics of dnaK and groEL transcription for T. whipplei exposed to heat shock. The relative levels of transcription of groEL2 (TWT441) and dnaK (TWT750) were determined by real-time RT-PCR using primers listed in Table 2. The data obtained were expressed as the ratio of the values obtained with bacteria exposed to heat shock to the values obtained with untreated cells grown at 37°C. The histograms are representative of two distinct experiments.
|
![]() View larger version (27K): [in a new window] |
FIG. 5. Evidence for two HAIR motifs upstream of the dnaK operon in T. whipplei. (A) Schematic representation of the dnaK operon, including six genes up-regulated with a 15-min heat shock at 43°C. Two homologues of the HAIR motif (18), designated TW_HAIR1 (5'-CATGAGTCGATATGACTCAAT-3') and TW_HAIR2 (5'-CTTGAGTCATTACATGTCAAG-3'), were identified upstream of this region. (B) ClustalW alignment with the HAIR motif initially described for Mycobacterium tuberculosis H37Rv (53).
|
![]() View larger version (77K): [in a new window] |
FIG. 6. Circular representation of the T. whipplei Twist transcriptome. The outermost (first) circle indicates the nucleotide positions. The second and third circles indicate the ORF locations on the plus and minus strands, respectively. The fourth, fifth, and sixth circles indicate the microarray expression profiles for T. whipplei at 4°C versus 37°C, at 28°C versus 37°C, and at 43°C versus 37°C, respectively. The expression profiles for each gene, which are shown with a color-coded, base 2 logarithmic scale, were independently centered about zero. Green indicates decreased expression relative to the expression at 37°C. Red indicates increased expression relative to the expression at 37°C. Yellow indicates no change in expression relative to the expression at 37°C.
|
|
|
|---|
This analysis also highlighted the finding that under heat shock conditions, the ribC gene is overexpressed (12-fold increase). This gene is implicated in the biosynthesis pathway for riboflavin (vitamin B2), a vitamin not synthesized by humans, and is involved in colonization and persistence of Helicobater pylori (13). Bacterial entry into a host cell corresponds to significant environmental changes that are mimicked in part by heat shock, which has been associated with up-regulation of bacterial virulence factors (33). In this respect, we believe that ribC could play a role in T. whipplei pathogenicity. Other putative virulence factor could be encoded by ispDF (3.8-fold increase), which has been described as a potential drug target in several human pathogens, including H. pylori, Campylobacter jejuni, and Treponema pallidum (16). In prokaryotes, heat shock also mimics antibiotic stress. This fits well with enhanced expression of folE, a gene involved in folate biosynthesis. Indeed, sulfonamide-sensitive bacteria, like T. whipplei (7), must synthesize their own folates (12). Finally, the function of infC, which is mainly associated with a cold-sensitive phenotype (23), is unclear.
Evidence for a complex adaptive response of T. whipplei to cold shock. In experiments aimed at identifying other putative virulence factors of T. whipplei (33), we analyzed the transcriptional variations induced by 6 h of incubation of T. whipplei at 28°C. In order to obtain a better understanding of the effect of cold resistance of T. whipplei, we also incubated a mid-log bacterial suspension at 4°C for the same time. The results obtained showed that all genes regulated at 28°C were regulated at 4°C but that the magnitude of regulation was higher, indicating that 28°C is lower than the optimal temperature for growth of T. whipplei in vitro. Therefore, we did not identify virulence genes which could be down-regulated at this temperature.
Experiments carried out with severe cold conditions (i.e., 4°C) revealed that while T. whipplei lacks classical CSP, this microorganism has numerous specific adaptive mechanisms for responding to cold. The most striking feature is the paradoxical up-regulation of the GroEL2 transcript. This transcript encodes a member of the HSP family and is expected to be down-regulated at low temperatures (5, 25). While unconventional, an increase in transcription of HSP60 genes at low temperatures was previously described for the hyperthermoacidophilic archaeon Sulfolobus shibatae, which can adapt to and survive in extreme natural environments (30). Another chaperonin-encoding gene, clpP1, was also up-regulated, and its role in bacterial tolerance of low temperatures has been suggested previously (14). These HSP were not modified upon heat shock. Whether T. whipplei uses these HSP to compensate for a CSP defect is questionable. It could be hypothesized that in T. whipplei, such chaperonines ensure protein folding in thermal or acidic environments that are lethal to other microorganisms that have been characterized. Other major changes are related to the fatty acid pathway (bccA, fabD, kasA, cdsA), a feature in line with membrane adaptation described for other bacteria exposed to cold. As the physical barrier between living cells and their environment, the plasma membrane is indeed highly susceptible to changes in environmental temperatures. While the lipid bilayers of most organisms are mostly fluid at physiological temperatures, they undergo a reversible change of state at lower temperatures (38). T. whipplei probably circumvents alterations in the cell membrane by incorporating fatty acids that have lower melting points in order to restore membrane fluidity and function. Control of the membrane physical properties can also be afforded by observed changes in the transcriptional levels of several genes encoding the WISP proteins, which are specific T. whipplei membrane proteins (6).
As reviewed by Gualerzi et al. (23), the so-called "winter shopping list" is highly specific for each bacterial species based on its lifestyle. Thus, analysis of the T. whipplei transcriptome at 4°C revealed that there was 1.5-fold or greater up-regulation of several genes involved in energy metabolism, a situation theoretically encountered when bacteria replicate. Some genes encoding ribosomal proteins (rpmB, rplT, rpsN, and rplP) were also up-regulated, a result consistent with the fact that a temperature downshift induces stabilization of secondary structures of nucleic acids, leading to reduced efficiency of mRNA translation and transcription (28). Therefore, other ribosomal proteins encoding genes (rpsK, rpsL, and rplK) were down-regulated to the same extent; the reason for such antagonistic regulation is unclear. Cold shock also promotes down-regulation of rph and rng, which encode two RNA stabilization factors. At the same time, both rnc (RNase III) and pcnA [poly(A) polymerase] are up-regulated. RNase III degrades and processes selective RNAs during starvation (2), and poly(A) polymerase synthesizes a poly(A) tail on some specific mRNAs, promoting fast and selective degradation of mRNA-poly(A) (62). In bacteria, RNA degradation permits the cell to scavenge nucleotides for resynthesis of RNAs. We speculate that T. whipplei cold shock adaptation is accompanied by repression and activation of key enzymes which act in concert to induce specific changes in gene expression that lead to metabolic adaptation.
Up-regulation of the purM, purC, purF, and pyrB genes involved in purine/pyrimidine pathways (purF and purM form a stimulon) was also observed. Overexpression of these enzymes associated with down-regulation of prs leads to a decrease in the 5-phosphoribosyl 1-pyrophosphate pool in T. whipplei. In the same nucleotide metabolism pathway (http://www.genome.jp/kegg/pathway/map/map00230.html), deoD (EC 2.4.2.1), which encodes purine nucleoside phosphorylase, is repressed and might led to an increase in the amount of guanosine 3',5'-bispyrophosphate (ppGpp). The amount of this alarmone is decreased following temperature downshifts in Escherichia coli (29), but up-regulation of alarmone has been reported for Bacillus subtilis (57). Whether (p)ppGpp regulation plays a role in the cold shock response remains to be determined.
The T. whipplei genome has a large deficiency of genes encoding enzymes involved in amino acid synthesis (44, 46). This explains the finding that, in contrast to what was observed for most bacteria that have been studied (20), the levels of only two genes belonging to this category were increased after a cold shock. A whole up-regulated cold regulon composed of genes classified as genes that encode putative ABC transporters led us to speculate that T. whipplei keeps in storage nutrients essential for cold hibernation.
T. whipplei dnaK operon is controlled by HspR. We demonstrated that in T. whipplei, the dnaK operon and groEL2 are regulated in opposite ways, and we found two HAIR operator cassettes upstream of dnaK. These conserved structural motifs are probably recognized by the HspR repressor also present in the T. whipplei dnaK regulon. While a few examples of the HspR repressor-HAIR operator system have been described for gram-positive bacteria, such a regulation mechanism is widespread in actinomycetes (18, 52) and has also been found recently in Deinococcus radiodurans (49) and C. jejuni (1). Of note is the fact that HSP themselves can participate in the regulation process. Thus, autoregulation of DnaK acting in cooperation with HspR was demonstrated in Streptomyces coelicolor (8). No other HAIR sites were detected in the T. whipplei genome.
Thus, how T. whipplei deals with temperature variation stresses is rather confusing. The data obtained demonstrated that only minor transcriptional modifications regulated by an ancestral negatively regulated mechanism occurred at 43°C. Similar patterns have been reported for both M. genitalium (580 kb) (39) and B. aphidicola (641 kb) (60). This result is consistent with the fact that, like other reduced-genome bacteria, T. whipplei has lost most of the regulatory sequences and protein regulators (34, 44). In this context, the strong cold shock response of T. whipplei was striking. We previously observed that T. whipplei was highly resistant to freezing and remained viable for 3 years at 80°C without cryoprotective agents (unpublished data). Our results emphasize the fact that microorganisms can exhibit distinct adaptive pathways involving regulatory cascades which are still not known and which may play a role in the survival of reduced-genome bacteria. While gene deletion should be a complementary approach for exploring putative regulatory mechanisms, genetic manipulation of T. whipplei has not been successful so far. Therefore, transcriptional analysis of T. whipplei (10) is likely to reflect the adaptability of this organism having an environmental origin and could increase our understanding of this deadly and poorly known pathogen.
We acknowledge the excellent support of the Nice-Sophia Antipolis Transcriptome Platform of the Marseille-Nice Genopole, in which the microarrays were constructed.
Supplemental material for this article may be found at http://jb.asm.org/. ![]()
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»