Journal of Bacteriology, January 2004, p . 556-565, Vol . 186, No . 2

Comparative Genomics of Rickettsia prowazekii Madrid E and Breinl Strains

Hong Ge,^1,2 Yao-Yu Eric Chuang,³ Shuping Zhao,³ Min Tong,^1,2 Mong-Hsun Tsai,³ Joseph J . Temenak,^1,4 Allen L . Richards,¹ and Wei-Mei Ching^1,2*

Rickettsial Diseases Department, Infectious Diseases Directorate, Naval Medical Research Center, Silver Spring, Maryland 20910,¹ Department of Preventive Medicine and Biometrics, Uniformed Services University of The Health Sciences, Bethesda, Maryland 20814,² Microarray Laboratory, Radiation Oncology Sciences Program, National Cancer Institute, National Institutes of Health, Gaithersburg, Maryland 20877,³ Division of Vaccines and Related Products Applications, Office of Vaccines Research and Review, Center for Biologics Evaluation and Research, Food and Drug Administration, Rockville, Maryland 20852⁴

Received 22 August 2003/ Accepted 10 October 2003

 
Rickettsia prowazekii, the causative agent of epidemic typhus, has been responsible for millions of human deaths . Madrid E is an attenuated strain of R . prowazekii, while Breinl is a virulent strain . The genomic DNA sequence of Madrid E has recently been published . To study the genomic variations between Madrid E (reference) and Breinl (test) DNAs, cohybridization experiments were performed on a DNA microarray containing all 834 protein-coding genes of Madrid E . Of the 834 genes assessed, 24 genes showed 1.5- to 2.0-fold increases in hybridization signals in Breinl DNA compared to Madrid E DNA, indicating the presence of genomic variations in 3% of the total genes . Eighteen of these 24 genes are predicted to be involved in different functions . Southern blot analysis of five genes, virB4, ftsK, rfbE, lpxA, and rpoH, suggested the presence of an additional paralog(s) in Breinl, which might be related to the observed increase in hybridization signals . Studies by real-time reverse transcription-PCR revealed an increase in expression of the above-mentioned five genes and five other genes . In addition to the elevated hybridization signals of 24 genes observed in the Breinl strain, one gene (rp084) showed only 1/10 the hybridization signal of Madrid E . Further analysis of this gene by PCR and sequencing revealed a large deletion flanking the whole rp084 gene and part of the rp083 gene in the virulent Breinl strain . The results of this first rickettsial DNA microarray may provide some important information for the elucidation of pathogenic mechanisms of R . prowazekii .

 
Rickettsiae are obligate intracellular gram-negative bacteria that belong to the alpha subdivision of Proteobacteria (48) . Rickettsial diseases are widely distributed throughout the world and are characterized by sudden onset of acute fever, skin eschars, headache, and, in cases of the fulminant disease, disseminated intravascular coagulation that leads to failure of the cardiovascular system . Diagnosis of rickettsial diseases is difficult due to the lack of specific signs and symptoms and the inability to grow the bacterium on axenic microbiological media for confirmation (32) . There are no licensed vaccines available for the prevention of rickettsial diseases in the United States . Moreover, three rickettsial agents, Rickettsia prowazekii, Rickettsia rickettsii, and Coxiella burnetii, have been listed as biothreat agents by the Centers for Disease Control and Prevention in the United States .

R . prowazekii is the etiologic agent of epidemic typhus in humans . Its presence is commonly associated with war, famine, and social disruption . Typhus killed >3 million people and infected in excess of 20 million people during World War I, and it killed several million more in World War II (20) . As recently as 1997, there was an epidemic of typhus in Burundi that infected well over 100,000 people (31) . Other outbreaks in Eastern Europe, Russia, Ethiopia, and South America further support the view that epidemic typhus is a reemerging infectious disease worldwide (20, 21) . Present-day threats posed by typhus bacteria include natural exposure to antibiotic-resistant strains and the use of laboratory-manipulated strains as biological warfare agents (28) . The World Health Organization has estimated that 104,000 casualties (19,000 dead and 85,000 incapacitated) would arise following a hypothetical air deployment of 50 kg of typhus agent (50) .

The pathogenic mechanism of rickettsiae is unclear but may involve potential virulence factors, such as phospholipases (37, 46, 49), antigenic surface proteins (10, 19, 24, 46), hemolysins (29), type IV secretion systems (2, 25, 26), and the recently described nudix hydrolase (18) . Significant biological differences exist between the virulent Breinl strain of R . prowazekii and the attenuated avirulent Madrid E strain (42) . The Breinl strain grows well in macrophage-like cells, while the Madrid E strain grows poorly . The Breinl strain has a greater cytotoxic effect on macrophage-like cells treated with alpha/beta interferon (IFN-/ß) or IFN- than the Madrid E strain, and the Breinl strain is less sensitive to IFN- treatment than the Madrid E strain in L929 cells (42) . However, the genes responsible for the variation in virulence between the two strains of R . prowazekii have not been well characterized . The published genomic sequence of the Madrid E strain (2) provides a foundation for comparative genomic studies of R . prowazekii. The circular genome consists of 1.1 Mbp containing 834 protein-coding genes, 12 pseudogenes, and 38 RNA genes . About 24% of the genome belongs to the noncoding region, which represents the highest proportion of noncoding DNA detected in any microbial genome sequenced so far (2) . The comparisons of functional profiles of Madrid E genes with those of mitochondrial genes revealed many similarities in ATP production and transport systems (2) .

The purpose of these studies was to compare the genomic content of the avirulent Madrid E strain with that of the virulent Breinl strain . We report the construction of an R . prowazekii DNA microarray with 834 PCR-amplified gene products from Madrid E and the comparative genomic studies with DNAs from Madrid E and Breinl using this microarray . The microarray results led to the discovery of a large deletion including the rp084 gene in the Breinl strain . Southern blot analysis of five dissimilar genes revealed the possible presence of additional paralogs, and the increased expression of 10 virulence-related genes was confirmed by real-time reverse transcription (RT)-PCR . Elucidation of the roles of these genes could lead to a better understanding of the mechanisms for the pathogenesis of R . prowazekii . In addition, insights gained from these studies may have implications for the development of new diagnostic strategies, treatment targets, and vaccine candidates for epidemic typhus .

 
Rickettsial strains and purifications. R . prowazekii strain Breinl was originally isolated from a patient in Warsaw, Poland, in 1922 and was brought to the United States in 1928 . This strain was passaged multiple times and used as a representative strain of epidemic typhus . The strain used in this study was in its 157th passage . Madrid E was the egg passage-attenuated strain of the Madrid strain, which was originally isolated from a typhus patient in Madrid (2) . The Madrid E strain used here was in its 281st passage . The exact passage histories of the two strains are unknown because of the lack of data collected early in the cultivation of the strains .

Six-day pathogen-free chicken embryos were inoculated with yolk sac suspensions from ampoules and incubated at 35°C . The infected yolk sacs were pooled, shaken with glass beads and further macerated by pipetting, diluted with brain heart infusion, and centrifuged at 10,000 x g for 30 min . The pellets were washed twice with buffer K36 (0.1 M KCl, 0.015 M NaCl, 0.05 M potassium phosphate, pH 7.0) by resuspension and sedimentation at 17,300 x g . The pellet was suspended in one-half volume of the original brain heart infusion, and 5 ml of this suspension was layered over 20 ml of 10% bovine albumin in buffer K36 . After centrifugation at 480 x g for 20 min, small sediments were discarded and the supernatants were diluted with 3 volumes of K36 and centrifuged at 10,000 x g for 30 min . The washed pellets were resuspended in K36 and subjected to two consecutive cycles of Renografin density gradient purification (13) . In both cases, 2.5 ml of rickettsial suspension was layered over the following gradients . (i) The first was a 30 to 45% linear gradient of Renografin for centrifugation at 25,000 rpm (SW27 rotor) for 1 h . After centrifugation, the areas of the gradient above and below the rickettsial band were removed with a Pasteur pipette . Then, the rickettsial band was drawn into a syringe through a 14-gauge cannula, diluted 10 times with K36, and centrifuged at 17,300 x g for 10 min . The pellets were suspended in K36 . (ii) The second was a 20 to 45% linear gradient of Renografin with centrifugation at 25,000 rpm (SW25.1 rotor) . The suspension from the first gradient centrifugation was processed through the second gradient centrifugation . The pure rickettsial band was drawn, and the final pellet was suspended in K36 . Live rickettsiae purified by this method were essentially free of host material .

Construction of microarray slides . (i) DNA extraction and amplification. Renografin-purified R . prowazekii Madrid E was resuspended in ATL lysis buffer (Qiagen, Valencia, Calif.) . After proteinase K digestion and RNase A treatment according to the instructions of the manufacturer, DNA was extracted three times with phenol-chloroform-isoamyl alcohol and precipitated with ethanol-sodium acetate . PCR primers for all 834 genes covering the entire open reading frame (ORF) of each gene from R . prowazekii Madrid E were purchased from Research Genetics (Huntsville, Ala.) . PCRs were carried out in 100-µl reaction mixtures using a 96-well PE-9700 thermocycler (Applied Biosystems, Foster City, Calif.) . The reaction mixture contained 15 ng of DNA, 0.1 mM deoxynucleoside triphosphate, 100 ng of each primer, 1x PCR buffer, 1.5 mM MgCl₂, and 2.5 U of Taq polymerase (Invitrogen-Life Technology, San Diego, Calif.) . The reaction conditions were as follows: initial denaturation at 94°C for 30 s and 33 cycles of denaturation (94°C for 45 s), annealing (53°C for 45 s), and extension (72°C for 45 s) . Final extension was performed at 72°C for 10 min . Most PCR products were 0.5 to 2 kb in size . For 147 genes that were not amplified or that yielded nonspecific products or faint bands under the conditions described above, the PCR conditions were adjusted or new primers were designed .

(ii) Preparation of glass slides for microarray. Microscope slides were selected manually in order to exclude those slides with dust or scratches . The slides were submerged in NaOH-ethanol for 2 h and then coated with polylysine after the residual NaOH was completely removed by vigorously washing the slides for 20 min in double-distilled H₂O . The slides were immediately spun at 550 rpm for 3 min and then stored in desiccators at room temperature until spotting was performed .

(iii) Spotting and posttreatment of slides. All PCR products were purified by ethanol-acetate precipitation and robotically printed onto polylysine-coated slides by the computerized OMNI Grid Arrayer (GeneMachines, San Carlos, Calif.) . The eight-pin model was used for spotting the DNAs of 834 genes plus a 96-gene duplicate as the internal control . The spotted DNA was cross-linked to polylysine by UV irradiation and treated with succinic anhydride to minimize background staining .

Hybridization of genomic DNA. Genomic DNAs from the Madrid E and Breinl strains were used as templates for direct incorporation of fluorescent nucleotide analogs (Cy3 and Cy5 dCTPs; Amersham, Piscataway, N.J.) by a random-priming polymerization reaction (6) . After denaturation of the DNA, the two 50-µl labeling reaction mixtures containing 2 µg of template DNA (one from Madrid E and another from Breinl); 1x random-primer reaction buffer (Invitrogen-Life Technology); 1.2 mM (each) dATP, dGTP, and dTTP; 0.6 mM dCTP; 3.75 nmol of Cy3-dCTP or Cy5-dCTP; and 40 U of Klenow fragment (Invitrogen-Life Technology) were incubated at 37°C for 2 h . By adding Cot-1 DNA, yeast tRNA, and Tris-EDTA buffer, the labeled DNA was then concentrated, denatured, and applied to microarray slides . The hybridization was carried out at 65°C overnight, and then the slides were washed for 5 min in 2x SSC buffer (1x SSC is 0.15 M NaCl plus 0.015 M sodium citrate) with 0.03% sodium dodecyl sulfate, followed by 5 min in 1x SSC and then 1 to 2 min in 0.2x SSC buffer . In a subsequent experiment, each DNA was reciprocally labeled in order to eliminate artifacts due to differences in labeling intensity .

Imaging and data analysis. The hybridized arrays were scanned at 10-µm resolution on a GenePix 4000 scanner (Axon Instruments, Inc., Foster City, Calif.) at variable photomultiplier tube voltage settings to obtain maximal signal intensities with <1% probe saturation . The Cy5- and Cy3-labeled DNA samples were scanned at 635 and 532 nm, respectively . The resulting TIFF images were analyzed by Gene Pix Pro version 3.0 . For each spot, the background intensity was subtracted from the spot signal intensity to produce a channel-specific value . The Cy3/Cy5 ratios for all targets were determined after background subtraction . Gene cluster analysis was performed using the National Cancer Institute-Center for Information Technology microArray database system . Hierarchical clustering was performed according to the method of Eisen et al . using an "uncentered" correlation algorithm (14) . The cluster data were filtered for signal intensities and background of >1.5-fold in both channels . The raw data were converted into log₂ data for the analysis of the genomic distribution versus log₂ ratio .

Southern blot analysis. Southern blot hybridization was carried out by using standard protocols (35) with the following modifications . Probes were prepared using a PCR DNA biotinylation kit (Kirkegaard & Perry Laboratories, Gaithersburg, Md.) that yielded products between 0.5 and 1.0 kb in size . MacVector version 7.0 software (Genetics Computer Group, Madison, Wis.) was used to select restriction enzymes for the digestion of genomic DNA to generate proper length distributions and to predict signal sizes in Southern blot analysis . Rickettsial DNA was digested with EcoRI/XbaI and HindIII for the detection of ftsK, rfbE, lpxA, and rpoH . HindIII was replaced by EcoRV/ClaI for the detection of virB4. None of these restriction enzymes cut Madrid E DNA sequences that were designed for the generation of probes . The DNA fragments were separated on 0.8% agarose gels by overnight electrophoresis and then transferred onto a positive nylon membrane (Kirkegaard & Perry Laboratories) overnight by alkaline capillary transfer . The labeled DNA probe was hybridized to the membrane, and signal was detected with a Random Primer DNA biotinylation kit (Kirkegaard & Perry Laboratories) according to the manufacturer's recommendations .

PCR and RT-PCR for rp084 gene. PCR and RT-PCR were performed with primers specific for the entire rp084 gene (forward, 5' ATG ACT ATA ACA GGC AAA TT 3'; reverse, 5' TAT TTT TCT TTT GCC TCA CT 3') . The reaction conditions were as follows: initial denaturation at 94°C for 30 s; 32 cycles of denaturation (94°C for 45 s), annealing (52°C for 45 s), and extension (72°C for 45 s); and a final extension step at 72°C for 10 min . Amplification of the 16 S rRNA gene was done as an internal control under the same conditions as for rp084 with the following primers: 5' GTT CGG AAT TAC TGG GCG TA 3' (forward) and 5' AAT TAA ACC GCA TGC TCC AC 3' (reverse) .

DNA sequencing. PCR primers that flanked the sequence of Madrid E from the end of rp082 to the beginning of rp085 were designed . A PCR product from Madrid E was calculated to be a 4.3-kb fragment . The PCR product from Breinl was 2.5 kb and was directly sequenced with the Big-Dye terminator cycle-sequencing kit using the automated sequencing machine (Applied Biosystems) . In addition, the full-length genes virB4, ftsK, lpxA, rfbE, and rpoH of Breinl and Madrid E were PCR amplified and cloned into the BlueScript vector (Stratagene, La Jolla, Calif.) and sequenced as described above . DNA sequence analyses and alignments were performed with MacVector version 7.0 and Blast software .

RNA extraction and RT. Total RNA was isolated from the Renografin-purified rickettsial strains Madrid E and Breinl with TRIZOL (Invitrogen-Life Technology) according to the manufacturer's instructions . The RNA samples were treated with DNase I (Invitrogen-Life Technology) at 37°C for 1 h and then extracted with phenol-chloroform-isoamyl alcohol two or three times . The quality and quantity of RNA were confirmed by visualization of the agarose gel and spectrophotometer analysis . Two micrograms of total RNA was used for cDNA synthesis at 42°C for 50 min with SuperScript II (Invitrogen-Life Technology), followed by RNase H treatment at 37°C for 20 min . A reaction with no reverse transcriptase was included as a control . All procedures were performed for both strains under identical conditions to allow comparison .

Real-time RT-PCR. The primer pairs and probes of selected genes for real-time PCR were designed using Primer Express software (ABI, Weiterstadt, Germany) as shown in Table 1 . The probes were labeled with 6-carboxyfluorescein and 6-carboxytetramethyl-rhodamine at the 5' and 3' ends, respectively . PCR cycling was conducted using the SmartCycler (Cepheid, Sunnyvale, Calif.) in reaction mixtures with total volumes of 25 µl containing 1x TaqMan Universal PCR Master Mix (Applied Biosystems), 0.2 µM each primer, 0.1 µM probe, and 1 µl of diluted cDNA template . After incubation at 50°C for 2 min and denaturation at 95°C for 10 min, the reaction went to 45 to 60 cycles of denaturation (94°C for 15 s) and extension (60°C for 1 min) . The 16S rRNA gene was included as an endogenous reference, and the comparative C_T (threshold cycles) method was applied by using arithmetic formulas (2^-C_T) (ABI system; AB, Piscataway, N.J.) . In this method, the amount of target was normalized to that of the reference gene relative to the calibrator . The efficiency of the target amplification and the efficiency of the reference amplification were tested by assessing the variations of C_T at different amounts of template from the reference and sample . No-cDNA templates in PCR and RT products without reverse transcriptase were used as controls . Real-time RT-PCR assays were performed in triplicate for each sample, and a mean value and standard deviation were calculated for the relative RNA expression levels .

 
Comparison of genomic DNAs from Madrid E and Breinl by microarray. The cohybridization of Cy3-labeled Madrid E DNA (green fluorescence signal) and Cy5-labeled Breinl DNA (red fluorescence signal) was applied on the glass microarray slides . The resulting signal intensities from the two fluorophores were compared for each gene in order to detect the genomic differences, e.g., DNA duplications or deletions, in the two strains . It should be noted that no plasmid has been found in R . prowazekii . The differences in signal should not be due to the presence, absence, or different copy numbers of a plasmid . Our results demonstrated that R . prowazekii genes are highly conserved between the two strains, because most of them appeared to have yellow fluorescence signals (Fig . 1A to D) . However, quite a few genes appeared to have varied signals . Figure 1E shows a distribution of log₂ ratio versus genomic location for all genes studied . Those with a fluorescence ratio of 1.5 or 0.667 were selected for further evaluation . Ratios between 0.67 and 1.5 are not reliable for predicting the presence of differences in gene content in microarray experiments . A total of 24 genes showed >1.5-fold- but <2.0-fold-higher hybridization signals, and one gene showed <0.667-fold-lower hybridization signals, with virulent Breinl DNA (Fig . 2) than with Madrid E DNA . Reciprocal-labeling experiments were conducted to rule out the possibility that the difference in intensity was due to a variation in Cy3 or Cy5 incorporation . The ratio of fluorescence for rp141 (rpoC) originally showed <0.667-fold intensity but was confirmed to be >1.5 (1.702 and 1.601) in two subsequent experiments (Fig. 2) . Thus, only one gene showed <0.667-fold lower in the ratio of the Breinl strain to the Madrid E strain . These 25 genes accounted for 3% (25 of 834) of all genes detected on the array . Eighteen out of 25 genes fell within functional categories based upon sequence homology to other bacterial genes (2) . The remaining seven genes are unmatched to any homologous function (Table 2) (2) .

 
Since the fluorescence ratios of 24 genes were >1.5, we evaluated the possibility of additional gene copies in Breinl compared to Madrid E in light of a report of different copy numbers of homologous genes . There are 54 paralogous gene families comprising 147 gene products reported in the genome of R . prowazekii Madrid E (2) . This observation may indicate specific adaptations of bacteria to their environments . The presence of four copies of SpoT in R . prowazekii and Rickettsia conorii might be related to the adaptation of Rickettsiae to environmental starvation in their vectors (25) . Other possible factors contributing to the increased fluorescence signals may be differences in the number of gene repeat sequences, gene insertions in Breinl, or deletion of genes in Madrid E .

DNA sequencing and Southern blot analysis. To rule out sequence variation among the 25 genes in Breinl and Madrid E, 5 of them (ftsk, rfbE, lpxA, rpoH, and virB4) were cloned and sequenced . The results demonstrated that the sequences of all five genes were identical . No insertions, duplications, deletions, or mutations were seen .

To evaluate the possibility of multiple copies or paralogs of the disparate genes and to explain why the hybridization signal ratios were between 1.5 and 2.0 with these genes, Southern hybridizations were performed for the above-mentioned five genes (Fig . 3) . Total DNAs were isolated from the R . prowazekii Madrid E and Breinl strains, respectively, and then digested with restriction enzymes, which were selected by MacVector version 7.0 software based on the published Madrid E sequence . The restriction enzymes were predicted to cut frequently within Madrid E DNA but not within the regions for probes . The Southern blots were hybridized with five DNA probes synthesized by PCR . The resulting blots displayed one strong band each for ftsK, rfbE, lpxA, and rpoH in all lanes of Madrid E and Breinl DNAs digested with EcoRI/XbaI and HindIII, which was consistent with the predicted single copies of these genes in the Madrid E genome (2) . However, the blots of Breinl DNA digested by EcoRI/XbaI appeared to have an additional weak band for all four genes . These additional bands were not seen in the Breinl strain after HindIII digestions . The blots hybridized with the rp103 probe for the detection of virB4 displayed two major bands in Madrid E . One was located at the expected position of 3,596 bp containing rp103. The second band resulted from the hybridization of the probe to another virB4 (rp784) and was located at a predicted position of 1,426 bp . However, the blots from Breinl DNA showed a third band at 1.3 kb, in addition to the two similar bands seen in Madrid E for virB4 . The genome sequence for Madrid E has confirmed two paralogs (rp103 and rp784) of virB4, which showed 30.9 and 20.3% DNA identities to Agrobacterium tumefactions virB4 (2) . Further sequencing alignments for rp103 and rp784 showed only 49 and 25% identities in DNA and protein, respectively, dispersed over the whole ORF, but both contain a consensus VirB4 domain . These findings with virB4 demonstrated that probe rp103 could hybridize to rp784 in addition to itself . Correspondingly, the third band in Breinl may represent another hybridization between rp103 and a paralog of virB4 . All these data support the notion that there may be an additional paralog(s) of virB4, ftsK, rfbE, lpxA, and rpoH in the Breinl strain that is not found in the Madrid E strain . Since paralog genes are not identical in their DNA sequences, despite belonging to the same functional family, it is conceivable that the signals due to paralogous hybridization may be not as strong as those due to homologous hybridization . We are not sure why there are no weak bands due to paralogs in Breinl DNA for ftsK, rfbE, lpxA, and rpoH after HindIII digestion, as well as for virB4 after EcoRV/ClaI digestion . It is possible that the paralogs have restriction cutting sites within the probe-complementary regions, or it may simply be that due to weak hybridization the signals were too faint to be seen . These data collectively suggested that an additional paralog(s) of the five genes might exist in the virulent Breinl strain, providing an explanation for the observation of gene hybridization ratios between 1.5 and 2.0 from microarray experiments .

 
Real-time RT-PCR analysis. To estimate the expression profiles of these genes in yolk sac-grown R . prowazekii, we conducted a comparative quantification of the gene expression profiles for 10 of the genes with detected differences in the microarray . Specific primers and probes were designed for each gene in order to perform real-time PCR . The target sequences were amplified from both Madrid E and Breinl . Control experiments without reverse transcriptase generated weak signals with C_T values near 60 cycles . Other controls without cDNA templates were negative in all experiments . The collected data were analyzed using the standard comparative 2^-C_T method . With this method, the expression profiles of the two strains were compared using the 16S rRNA gene as the internal control . The results showed that all real-time PCRs were positive for the corresponding target sequences with C_T values from 27 to 42 . The expression levels of all 10 genes were higher in the virulent Brienl strain than in the avirulent Madrid E strain, ranging from 1.3 to 9.9 times higher . The genes with twofold-higher or greater expression in Breinl were virB4, rfbE, lpxA, ampG2, pal, ileS, and uvrB . These results are summarized in Table 3 . The real-time PCR and microarray data for these 10 genes are illustrated in Fig. 4 . The expression ratios of lpxA, ampG2, pal, and uvrB were higher than those of other genes compared to the DNA ratios by microarray experiments, suggesting different transcriptional regulations in Breinl than in Madrid E .

 
Analysis of rp084. One of the intriguing observations in this study was that the rp084 gene exhibited only 1/10 (0.121, averaging from 0.158 and 0.084) of the intensity in Breinl DNA as in Madrid E DNA (Fig . 1 and 2) . To verify the microarray result, amplification of rp084 by PCR and RT-PCR were performed for the two strains . Neither the DNA nor the mRNA of this gene was detected in Breinl, but both were present in Madrid E (Fig . 5A and B) . PCR amplification of the region flanking rp083, rp084, and nearby noncoding regions produced a 4.3-kb fragment in Madrid E (expected size) and a 2.5-kb product in Breinl . Sequencing analysis of the amplicon from Breinl revealed that there was a deletion of 1,812 bp, including part of rp083, all of rp084, the noncoding region between rp083 and rp084, and part of the noncoding region after rp084 (Fig . 5C) . We are not sure of the significance of this deletion of genomic DNA for the virulent Breinl strain . The presence of rp084 was detected by PCR amplification in another virulent strain of R . prowazekii, MADII, but not in the Ananier strain (data not shown) . Both rp083 and rp084 belong to the category of unknown function in Madrid E (2) . There were no obvious similarities in either the primary DNA or protein sequence after an exhaustive database search . Following a BLAST database search with the predicted mass (18,232 Da) and pI (9.52) for rp084, three mitochondrial proteins, CTAG (cytochrome c oxidase assembly protein), CYCJ (cytochrome c-type biogenesis protein), and OSCP (ATP synthase oligomycin sensitivity conferral protein), with the same pI and very similar molecular masses of 18 kDa were identified . They show similarities in the predicted secondary structures (data not shown) and putative biochemical characteristics, such as overall hydrophilicity, a high content of basic amino acids (higher pI value), and a predicted -helical secondary structure, all of which are indicative of mitochondrion-related proteins (17) . CTAG found in the rickettsia Cowdria ruminantium was shown to affect cytochrome c oxidase synthesis at a terminal stage (45) . CYCJ might be required for the biogenesis of c-type cytochromes in Bradyrhizobium japonicum and Rhizobium etli (33, 41) . In Drosophila melanogaster and Saccharomyces cerevisiae, OSCP was implicated in proton conduction during ATP synthesis (4, 9, 34) . Cytochrome c is believed to be an important factor in the initiation of the apoptotic cascade . The cytochrome c released in response to stress could activate a caspase initiator, caspase 9, and trigger the subsequent activation of downstream executioner caspases in the mammalian host (22, 52) . Further functional studies of rp084 are necessary to demonstrate the specific role of this rickettsial gene .

 
Potential virulence genes identified in this study. The virB4, ftsK, lpxA, and rfbE genes are known to be virulence or virulence-related genes, and their overrepresentation in the virulent Breinl strain might indicate that they have roles in the virulence of R . prowazekii . Although there are limited reports of the study of these genes in Rickettsia, they have been investigated in Escherichia coli and other bacteria .

(i) virB4. The R . prowazekii Madrid E genome possesses two virB4 paralogs (rp103 and rp784) that are homologous to the virB4 genes encoding bacterial transport systems, such as T-DNA export in A . tumefaciens, with DNA sequence identities of 30.9 and 20.3%, respectively (2) . VirB4 is believed to be an essential virulence protein in A . tumefaciens as a cytoplasmic membrane ATPase (11, 16, 23) . VirB4 also functions in plasmid transport in E . coli, in Pt1 toxin export in Bordetella pertussis, in the Cag system in Helicobacter pylori, and in the icm/dot system in Legionella pneumophila (reference 12 and references therein) . However, the role of virB4 in rickettsiae has not been elucidated . In addition, other virB family members have been reported in R . prowazekii Madrid E, such as virB8-11 and virD4 (2) .

(ii) ftsK. FtsK is predicted to be involved in bacterial cell division . FtsK may possess two functions: the membrane domain (FtsK_N) targets the protein to the closing septum (47), and the cytoplasmic portion (FtsK_C) segregates chromosome dimers during cell division (39, 40, 51) . A mutation in ftsK44 appeared to block the late phase in the cell division pathway (5) . This information suggests that the FtsK protein plays a crucial role in cell division by triggering chromosome segregation . There has been no report of the study of ftsK in rickettsiae .

(iii) lpxA and rfbE. lpxA encodes UDP-N-acetylglucosamine (UDP-GlcNAc) acyltransferase, which is required for the first step of lipid A (endotoxin) biosynthesis (30) . LpxA catalyzes the transfer of an R-3-hydroxyacyl chain from R-3-hydroxyacyl-acyl carrier protein to the glucosamine 3-OH of UDP-GlcNAc (1) . The growth and virulence of most gram-negative bacteria require lipid A (27) . It is also a potent activator of innate immunity in animal systems (30, 43) . The rfbE gene encodes an enzyme necessary for the biosynthesis of the O-antigen side chain in E . coli O157:H7 (7) . Rickettsial lipopolysaccharide appears to be an important immunogen, as well as a virulence factor (36, 38) . Thus, rfbE has a potential role in the pathogenesis of R . prowazekii .

(iv) . RpoH. The RNA polymerase sigma-32 factor (RpoH) has been characterized as a primary component of the heat shock response in the survival of E . coli (8) . In response to a sudden heat shock or other stresses, the levels of RpoH rise transiently through enhanced synthesis and protein stabilization . In R . prowazekii, sigma-73 factor has been overexpressed in E . coli, and it can stimulate transcription by promoting RNA polymerase binding (3) . Rickettsial heat shock proteins, including DnaK, GroEL, and GroES, have been detected by Western blot analysis (15), but the regulatory mechanism for heat shock or stress is unknown .

In conclusion, we have successfully assembled the first genomewide DNA microarray for R . prowazekii and used this array to compare genomic variation between the virulent Breinl and avirulent Madrid E strains . Differences in the number of paralogs may be present for a few genes, such as virB4, ftsK, rfbE, lpxA, and rpoH, in Breinl and Madrid E . The absence of the whole rp084 and partial rp083 genes in Breinl is of unknown relevance to rickettsial virulence . Our results in this report may provide new insights into the molecular pathogenesis of R . prowazekii, such as the differences in the numbers of paralog genes, that may be useful for the development of improved diagnostic assays and vaccines . Further investigations are planned using this microarray to perform differential expression profile analysis of rickettsial mRNAs from infected mammalian cells and arthropod vectors .

 
This work was supported by the Naval Medical Research and Development Command, Work Unit 6000.ID.RD.A0239 .

The opinions and assertions contained here are those of the authors and are not to be construed as official .

The Renografin-purified rickettsiae were prepared under Gregory A . Dasch's direction . We thank Edwin Oaks for his critical comments during the preparation of the manuscript .

 
* Corresponding author . Mailing address: Rickettsial Diseases Department, Infectious Diseases Directorate, Naval Medical Research Center, Silver Spring, MD 20910 . Phone: (301) 319-7438 . Fax: (301) 319-7460 . E-mail: chingw@nmrc.navy.mil .

1. Anderson, M . S., H . G . Bull, S . M . Galloway, T . M . Kelly, S . Mohan, K . Radika, and C . R . Raetz. 1993 . UDP-N-acetylglucosamine acyltransferase of Escherichia coli. The first step of endotoxin biosynthesis is thermodynamically unfavorable . J . Biol . Chem . 268:19858-19865 .
2. Andersson, S . G., A . Zomorodipour, J . O . Andersson, T . Sicheritz-Ponten, U . C . M . Alsmark, R . M . Podowski, A . K . Naslung, A . S . Eriksson, H . H . Winkler, and C . G . Kurland. 1998 . The genome sequence of Rickettsia prowazekii and the origin of mitochondria . Nature 396:130-140.
3. Aniskovitch, L . P., and H . H . Winkler. 1996 . Rickettsia prowazekii sigma factor sigma 73 can be overexpressed in Escherichia coli and promotes RNA polymerase binding and transcription . Microbiology 142:901-906.
4. Bateson, M., R . J . Devenish, P . Nagley, and M . Prescott. 1999 . Single copies of subunits d, oligomycin-sensitivity conferring protein, and b are present in the Saccharomyces cerevisiae mitochondrial ATP synthase . J . Biol . Chem . 274:7462-7466 .
5. Begg, K . J., S . J . Dewar, and W . D . Donachie. 1995 . A new Escherichia coli cell division gene, ftsK. J . Bacteriol . 177:6211-6222.
6. Behr, M . A., M . A . Wilson, W . P . Gill, H . Salamon, G . K . Schoolnik, S . Rane, and P . M . Small. 1999 . Comparative genomics of BCG vaccines by whole-genome DNA microarray . Science 284:1520-1523 .
7. Bilge, S . S., J . C . Vary, Jr., S . F . Dowell, and P . I . Tarr. 1996 . Role of the Escherichia coli O157:H7 O side chain in adherence and analysis of an rfb locus . Infect . Immun . 64:4795-4801.
8. Bukau, B. 1993 . Regulation of the Escherichia coli heat-shock responses . Mol . Microbiol . 9:671-680.
9. Caggese, C., G . Ragone, B . Perrini, R . Moschetti, V . De Pinto, R . Caizzi, and P . Barsanti. 1999 . Identification of nuclear genes encoding mitochondrial proteins: isolation of a collection of D . melanogaster cDNAs homologous to sequences in the Human Gene Index database . Mol . Gen . Genet . 261:64-70.
10. Ching, W . M., G . A . Dasch, M . Carl, and M . E . Dobson. 1990 . Structural analyses of the 120-kDa serotype protein antigens of typhus group rickettsiae . Comparison with other S-layer proteins . Ann . N . Y . Acad . Sci . 590:334-351.
11. Dang, T . A., and P . J . Christie. 1997 . The VirB4 ATPase of Agrobacterium tumefaciens is a cytoplasmic membrane protein exposed at the periplasmic surface . J . Bacteriol . 179:453-462.
12. Dang, T . A., X . R . Zhou, B . Graf, and P . J . Christie. 1999 . Dimerization of the Agrobacterium tumefaciens VirB4 ATPase and the effect of ATP-binding cassette mutations on the assembly and function of the T-DNA transporter . Mol . Microbiol . 32:1239-1253.
13. Dasch, A . G., and E . Weiss. 1977 . Characterization of the Madrid E strain of Rickettsia prowazekii purified by renografin density gradient centrifugation . Infect . Immun . 15:280-286.
14. Eisen, M . B., P . T . Spellman, P . O . Brown, and D . Botstein. 1998 . Cluster analysis and display of genome-wide expression patterns . Proc . Natl . Acad . Sci . USA 95:14863-14865 .
15. Eremeeva, M . E., W.-M . Ching, Y . Wu, D . J . Silverman, and G . A . Dasch. 1998 . Western blotting analysis of heat shock proteins of Rickettsiales and other eubacteria . FEMS Lett . 167:229-237.
16. Errington, J., J . Bath, and L . J . Wu. 2001 . DNA transport in bacteria . Nat . Rev . Mol . Cell Biol . 2:538-545.
17. Forsha, D., C . Church, P . Wazny, and R . O . Poyton. 2001 . Structure and function of Pet100p, a molecular chaperone required for the assembly of cytochrome c oxidase in Saccharomyces cerevisiae. Biochem . Soc . Trans . 29:436-441.
18. Gaywee, J., S . Radulovic, J . A . Higgins, and A . F . Azad. 2002 . Transcriptional analysis of Rickettsia prowazekii invasion gene homolog (invA) during host cell infection . Infect . Immun . 70:6346-6354 .
19. Gilmore, R . D., N . Joste, and G . A . McDonald. 1989 . Cloning, expression and sequence analysis of the gene encoding the 120 kD surface-exposed protein of Rickettsia rickettsii. Mol . Microbiol . 3:1579-1586.
20. Gross, L. 1996 . How Charles Nicolle of the Pasteur Institute discovered that epidemic typhus is transmitted by lice: reminiscences from my years at the Pasteur Institute in Paris . Proc . Natl . Acad . Sci . USA 93:10539-10540 .
21. Kelly, D . J., A . L . Richards, J . Temenak, D . Strickman, and G . A . Dasch. 2002 . The past and present threat of rickettsial diseases to military medicine and international public health . Clin . Infect . Dis . 34:S145-S169.
22. Li, P., D . Nijhawan, I . Budihardjo, S . M . Srinivasula, M . Ahmad, E . S . Alnemri, and X . Wang. 1997 . Cytochrome c and d ATP-dependent formation of Apaf-1/caspase-9 complex initiates an apoptotic protease cascade . Cell 91:479-489.
23. Liu, G., G . C . Draper, and W . D . Donachie. 1998 . Dimerization of the Agrobacterium tumefaciens VirB4 ATPase and the effect of ATP-binding cassette mutations on the assembly and function of the T-DNA transporter . Mol . Microbiol . 29:893-903.
24. McDonald, G . A., R . L . Anacker, and K . Garjian. 1987 . Cloned gene of Rickettsia rickettsii surface antigen: candidate vaccine for Rocky Mountain spotted fever . Science 235:83-85.
25. Ogata, H., S . Audic, P . Renesto-Audiffren, P . E . Fournier, V . Barbe, D . Samson, V . Roux, P . Cossart, J . Weissenbach, J.-M . Claverie, and D . Raoult. 2001 . Mechanisms of evolution in Rickettsia conorii and R . prowazekii. Science 293:2093-2098 .
26. Ohashi . N., N . Zhi, Q . Lin, and Y . Rikihisa. 2002 . Characterization and transcriptional analysis of gene clusters for a type IV secretion machinery in human granulocytic and monocytic ehrlichiosis agents . Infect . Immun . 70:2128-2138 .
27. Onishi, H . R., B . A . Pelak, L . S . Gerckens, L . L . Silver, F . M . Kahan, M . H . Chen, A . A . Patchett, S . M . Galloway, S . A . Hyland, M . S . Anderson, and C . R . Raetz. 1996 . Antibacterial agents that inhibit lipid A biosynthesis . Science 274:980-982 .
28. Rachek, L . I., A . M . Tucker, H . H . Winkler, and D . O . Wood. 1998 . Transformation of Rickettsia prowazekii to rifampin resistance . J . Bacteriol . 180:2118-2124 .
29. Radulovic, S., J . M . Troyer, M . S . Beier, A . O . Lau, and A . F . Azad. 1999 . Identification and molecular analysis of the gene encoding Rickettsia typhi hemolysin . Infect . Immun . 67:6104-6108 .
30. Raetz, C . R . H. 1996 . Bacterial lipopolysaccharides: a remarkable family of bioactive macroamphiphiles, p . 1035-1063 . In F . C . Neidhardt et al . (ed.), Escherichia coli and Salmonella: cellular and molecular biology . American Society for Microbiology, Washington, D.C.
31. Raoult, D., V . Roux, J . B . Ndihokubwayo, G . Bise, D . Baudon, G . Martet, and R . Birtles. 1997 . Jail fever (epidemic typhus) outbreak in Burundi . Emerg . Infect . Dis . 3:357-360.
32. Raoult, D., R . J . Birtles, M . Montoya, E . Perez, H . Tissot-Dupont, V . Roux, and H . Guerra. 1999 . Survey of three bacterial louse-associated diseases among rural Andean communities in Peru: prevalence of epidemic typhus, trench fever, and relapsing fever . Clin . Infect . Dis . 29:434-436.
33. Ritz, D., L . Thony-Meyer, and H . Hennecke. 1995 . The cycHJKL gene cluster plays an essential role in the biogenesis of c-type cytochromes in Bradyrhizobium japonicum. Mol . Gen . Genet . 247:27-38.
34. Rubinstein, J., and J . Walker. 2002 . ATP synthase from Saccharomyces cerevisiae: location of the OSCP subunit in the peripheral stalk region . J . Mol . Biol . 321:613-619.
35. Sambrook, J., D . W . Russell, and J . Sambrook. 2000 . Molecular cloning: a laboratory manual, 3rd ed . Cold Spring Harbor Laboratory Press, Cold Spring Harbor, N.Y.
36. Schramek, S., R . Brezina, and J . Kazar. 1977 . Some biological properties of an endotoxic lipopolysaccharide from the typhus group rickettsiae . Acta Virol . 5:439-441.
37. Silverman, D . J., L . A . Santucci, N . Meyers, and Z . Sekeyova. 1992 . Penetration of host cells by Rickettsia rickettsii appears to be mediated by a phospholipase of rickettsial origin . Infect . Immun . 60:2733-2740.
38. Smith, D . T., N . F . Conant, and H . P . Willett. 1963 . The Rickettsiaceae and Chlamydiaceae . Epidemic typhus, p . 842-849 . In H . Zinsser (ed.), Zinsser microbiology . Appleton-Century Crofts, New York, N.Y.
39. Steiner, W., and P . L . Kuempel. 1998 . Cell division is required for resolution of dimmer chromosomes at the dif locus of Escherichia coli. Mol . Microbiol . 27:257-268.
40. Steiner, W., G . Jiu, W . D . Donachie, and P . Kuempel. 1999 . The cytoplasmic domain of FtsK protein is required for resolution of chromosome dimers . Mol . Microbiol . 31:579-583.
41. Tabche, M . L., E . G . Garcia, J . Miranda, J . E . Escamilla, and M . Soberon. 1998 . Rhizobium etli cycHJKL gene locus involved in c-type cytochrome biogenesis: sequence analysis and characterization of two cycH mutants . Gene 208:215-219.
42. Turco, J., and H . H . Winkler. 1991 . Comparison of properties of virulent, avirulent, and interferon-resistant strains . Infect . Immun . 59:1647-1655.
43. Ulevitch, R . J., and P . S . Tobias. 1995 . Receptor-dependent mechanisms of cell stimulation by bacterial endotoxin . Annu . Rev . Immunol . 13:437-457.
44. Vogel, J . P., H . L . Andrew, S . K . Wong, and R . R . Isberg. 1998 . Conjugative transfer by the virulence system of Legionella pneumophila. Science 279:873-876 .
45. Waghela, S . D., F . R . Rurangirwa, S . M . Mahan, C . E . Yunker, T . B . Crawford, A . F . Barbet, M . J . Burridge, and T . C . McGuire. 1991 . A cloned DNA probe identifies Cowdria ruminantium in Amblyomma variegatum ticks . J . Clin . Microbiol . 29:2571-2577.
46. Walker, D . H., H . M . Feng, and V . L . Popov. 2001 . Rickettsial phospholipase A2 as a pathogenic mechanism in a model of cell injury by typhus and spotted fever group rickettsiae . Am . J . Trop . Med . Hyg. 65:936-942 .
47. Wang, L., and J . Lutkenhaus. 1998 . FtsK is an essential cell division protein that is localized to the septum and induced as part of the SOS response . Mol . Microbiol . 29:731-740.
48. Weiss, E., and J . M . Moulder. 1984 . The rickettsias and chlamydias, p . 687-704 . In N . R . Krieg and J . G . Holt (ed.), Bergey's manual of systematic bacteriology, vol . 1 . Williams & Wilkins, Baltimore, Md.
49. Winkler, H . H., and R . M . Daugherty. 1989 . Phospholipase A activity associated with the growth of Rickettsia prowazekii in L929 cells . Infect . Immun . 57:36-40.
50. World Health Organization. 1970 . Health aspects of chemical and biological weapons . World Health Organization, Geneva, Switzerland.
51. Yu, X . C., E . K . Weihe, and W . Margolin. 1998 . FtsK is a bifunctional protein involved in cell division and chromosome localization in Escherichia coli. J . Bacteriol . 180:6424-6428 .
52. Zou, H., W . J . Henzel, X . Liu, A . Lutschg, and X . Wang. 1997 . Apaf-1, a human protein homologous to C . elegans CED-4, participates in cytochrome c-dependent activation of caspase-3 . Cell 90:405-413.

Free Online Full-text Article

What Is Protein?, What Is Functional Genomics?, What Is Biofilm?, What Is Bioreactor?, What Is Yeast?, r, Microbes, i, Microbe, a, Bacteria, s, Bacterium, o, Microorganism, e, Bacillus, c, Yeasts, i, Escherichia coli, c, Candida albicans, n, Escherichia coli, n, Streptococci, e, Escherichia coli, n, Anaerobes, e, Cell suspensions, i, Escherichia coli, r, Anaerobic bacteria, c, Microorganisms, e, Bacteriological, r, Escherichia coli, o, Bacillus subtilis, n, Streptococcal, i, Enterobacteriacea, s, Salmonella typhimurium, e, Pasteurella, c, Microbial, o, S. cerevisiae




 

© 2005 Transgalactic Ltd (manufacturer of Bioscreen C software) | Privacy Statement | P.O. Box 1393, 00101 Helsinki, Finland, phone: +358 9 85172920, fax: +358 9 8749481, e-mail: microbiology@bionewsonline.com
 

Last modified: May 25, 2005