Previous Article | Next Article ![]()
Infection and Immunity, February 2006, p. 1339-1351, Vol. 74, No. 2
0019-9567/06/$08.00+0 doi:10.1128/IAI.74.2.1339-1351.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Center for Human Bacterial Pathogenesis Research, Department of Pathology, Baylor College of Medicine, Houston, Texas 77030,1 Center for Molecular and Translational Human Infectious Diseases Research, The Methodist Hospital Research Institute, and Department of Pathology, The Methodist Hospital, Houston, Texas 770302
Received 1 November 2005/ Accepted 18 November 2005
|
|
|---|
M5005_Spy_0680 mutant strain produced significantly larger abscesses after subcutaneous inoculation into mice, consistent with a hypervirulence phenotype. The mutant strain had significantly higher in vitro expression of several proven and putative virulence genes, including scpA, encoding a peptidase that inactivates complement protein C5a. Together, the data provide new information about previously uncharacterized GAS TCSs. |
|
|---|
To sense and respond to changes in external conditions, bacteria often use a specialized gene pair referred to as a two-component gene regulatory system (TCS). Changes in the environment are detected by a membrane-localized sensor kinase that transduces a signal to a receiver protein (response regulator), a process resulting in altered gene expression. TCSs have been well studied in many pathogenic bacteria, and they are known to regulate expression of many key virulence factors, such as extracellular toxins.
The eight GAS genomes that have been sequenced have 13 TCSs, including 12 that are well conserved in all strains (2, 3, 6, 9, 24, 32, 33, 36) (Fig. 1A). Only 3 of these 12 TCSs (fasBCA, covRS, and ihk/irr) have been studied in detail, and information is available about their involvement in host-pathogen interactions (7, 8, 15, 40). Genetic analysis of fasBCA mutant strains revealed that this TCS assists in tissue invasion and destruction by regulating expression of genes encoding secreted virulence factors, including streptokinase, streptolysin S fibrin, and fibrinogen binding proteins. Expression of the fasBCA genes in exponential phase is affected by environmental signals, such as a temperature change from 29°C to 37°C (34). The CovR/S TCS has been studied extensively and is known to be a major negative regulator of virulence genes (5, 7, 16). For example, inactivation of the CovR/S TCS significantly increases virulence for mice and resistance to complement-mediated opsonophagocytic killing in vitro (7, 16) and accordingly increases expression of key virulence factors, such as hasABC, skaA, sagA, speB, and speMF (5, 7). Expression microarray analysis found that CovR/S directly or indirectly influences expression of 15% of genes in the GAS genome (7). Expression of this TCS is altered by the external concentration of Mg2+ ions (11). It has been suggested that a low concentration of Mg2+ in extracellular body fluids prevents derepression of CovR-regulated genes and allows for their full expression in host tissues (11).
![]() View larger version (24K): [in a new window] |
FIG. 1. Two-component gene regulatory systems in group A Streptococcus. (A) Chromosomal location of genes encoding 13 TCSs in sequenced GAS genomes. Red color denotes the four TCSs analyzed in this study. Triangles denote the direction of transcription. (B) Chromosomal regions flanking the four TCSs studied. The open rectangles and triangles below M5005_Spy_0680, M5005_Spy_0784, M5005_Spy_0830, and M5005_Spy_1280 indicate the region of the gene deleted in each mutant. Narrow red lines with arrows denote putative regulatory interactions between response regulator, TCS promoter, and other binding sites. Red circles with + and symbols denote putative activation and repression by a response regulator, respectively. Numbers refer to open reading frame designations in serotype M1 strain M5005. The putative function of inferred proteins is based on homology searches. 0678, 5' nucleotidase; 0679, GTP pyrophosphokinase; 0680, TCS response regulator; 0681, TCS sensor kinase; 0682, mevalonate kinase; 0683, diphosphomevalonate decarboxylase; 0684, phosphomevalonate kinase; 0685, isopentenyl-diphosphate delta-isomerase; 0780 to 0783, PTS system, mannose/fructose family components IIA-D; 0784, TCS sensor kinase; 0785, TCS response regulator; 0786, iron(III)-binding protein; 0787, hypothetical protein; 0788, acetyltransferase; 0789, hydrolase (HAD superfamily); 0826 to 0829, spermidine/putrescine transport system PotABC and spermidine/putrescine binding protein; 0830, TCS response regulator; 0831, TCS sensor kinase; 0832, malate-sodium symport; 0833, NAD-dependent malic enzyme; 1277, arginine repressor ArgR; 1278, hypothetical protein; 1279, hypothetical cytosolic protein; 1280, TCS sensor kinase; 1281, TCS response regulator; 1282, peptide methionine sulfoxide reductase MsrA/MsrB; 1283, thiol:disulfide interchange protein TlpA; 1284, cytochrome c-type biogenesis protein CcdA.
|
In the present work, we used isogenic strains, expression microarray analysis, and mouse infection models to investigate four TCSs of unknown function in GAS. We focused our studies on TCSs that have only weak homology with TCS genes of known or inferred function in other bacteria.
|
|
|---|
Strain M5005 and its isogenic derivatives were grown in Todd-Hewitt broth (Difco Laboratories) supplemented with 0.2% yeast extract (THY medium) at 37°C in an atmosphere of 5% CO2-20% O2. THY agar or tryptose agar with 5% sheep blood (Becton Dickinson) was used as solid medium. THY agar supplemented with spectinomycin (150 µg/ml) was used for selection of antibiotic-resistant mutant derivative strains. Cloning experiments were performed with Escherichia coli Nova Blue (Novagen). Ampicillin (100 µg/ml) or spectinomycin (150 µg/ml) was used for selection of E. coli clones when required. The oligonucleotide primers used to construct mutant strains are presented in Table 1.
|
View this table: [in a new window] |
TABLE 1. Primers used in this study
|
Transformation of GAS.
Strain M5005 was grown overnight at 37°C in THY broth supplemented with 40 mM L-threonine and 250 mM sucrose, subcultured to fresh medium, and grown to an optical density at 600 nm (OD600) of 0.2 to 0.25. Bacteria were collected by centrifugation, washed three times with ice-cold 0.5 M sucrose, and resuspended in 0.5 M sucrose with 20% glycerol. Competent cells (100 µl) were mixed with dialyzed DNA and transferred to a cooled cuvette with a 0.1-cm electrode gap. After 30 min of incubation, the bacteria-DNA mixture was subjected to a pulse of 1.8 kV, 25 µF, 400
. The bacteria were mixed immediately with 1 ml THY containing 250 mM sucrose, incubated for 2 to 3 h at 37°C, and plated onto THY medium containing spectinomycin.
Construction of isogenic mutant strains with inactivated TCS response regulator genes. Four isogenic mutant strains were constructed, each having an inactivated response regulator gene (M5005_Spy_0680, M5005_Spy_0785, M5005_Spy_0830, and M5005_Spy_1281) of a TCS (Fig. 1B). These genes are designated SPy0874, SPy1062, SPy1106, and SPy1556, respectively, in the genome of strain SF370 (Fig. 1A). The gene numbers differ between strains SF370 and M5005 because the genome of the latter strain contains several large blocs of prophage DNA not present in the genome of strain SF370 (6, 36). To construct each mutant strain, the coding sequence of the response regulator gene was replaced with a spectinomycin resistance cassette by a double crossover strategy. DNA fragments encoding the upstream region of four response regulator genes were amplified by PCR with primer pairs (Table 1) 0680LF-0680LR (530 bp), 0785LF-0785LR (557 bp), 0830LF-0830LR (562 bp), and 1281LF-1281LR (465 bp). The purified gene fragments were digested with NdeI/XmaI and cloned into pUC19 (Fermentas) digested with appropriate enzymes to yield plasmids p0680-5', p0785-5', p0830-5', and p1281-5'. The downstream region of each gene encoding the target response regulators was amplified with primers 0680RF-0680RR (490 bp), 0785RF-0785RR (436 bp), 0830RF-0830RR (469 bp), and 1281RF-1281RR (466 bp), cloned into pSTblue-1 vector (Novagen), and digested with HincII/AvrII to generate plasmids p0680-3', p0785-3', p0830-3', and p1281-3'. The spectinomycin resistance cassette (spc) containing the add9 gene, which was shown to have no polar effects on downstream genes (19), was excised from plasmid pSL60-2 with SmaI and cloned into the PmlI site of each plasmid containing the downstream fragments of the response regulator genes. The orientation of the spc cassette was determined by PCR, and plasmid clones containing the spc cassette and downstream insert present in the same orientation were used for subsequent cloning steps. Fragments with the spc cassette fused with GAS DNA fragments were excised with Ecl136II/SnaBI and subcloned into the SmaI site of plasmids p0680-5', p0785-5', p0830-5', and p1281-5' containing the upstream region of each gene to be inactivated. The resulting plasmids (pM0680, pM0785, pM0830, and pM1281) have the upstream and downstream regions flanking each response regulator with the spc cassette cloned between them rather than the gene encoding the response regulator. The final recombinant DNA constructs were linearized with ScaI and transformed into electrocompetent strain M5005. GAS transformants were selected on solid medium containing spectinomycin. Chromosomal DNA isolated from spectinomycin-resistant colonies was used in PCR screening with pairs of primers 0680L-0680R, 0785L-0785R, 0830L-0830R, and 1281L-1281R. Primer pairs were designed to complement the chromosomal DNA sequence located outside of the targeted integration site. Putative recombinants were confirmed to contain the correct mutation by PCR analysis, DNA sequencing, and Southern hybridization using the spectinomycin cassette as a probe (Fig. 2). Successful inactivation of the target gene also was confirmed by the results of the microarray experiment, which showed the absence of transcripts of the inactivated response regulator gene (see below).
![]() View larger version (22K): [in a new window] |
FIG. 2. Confirmation of correct genetic construction of mutant strains. (A) PCR analysis of fragments amplified from chromosomes of mutant strains. Fragments were amplified with primers designed against target sequences outside of the chromosomal integration site. The size of the spc cassette differs from the size of the inactivated gene, which means that the sizes of fragments amplified from the mutant and parental strains differ, as shown. DNA fragments were amplified with primer pairs 0680L/R (lanes 1 and 2), 0785L/R (lanes 3 and 4), 0830L/R (lanes 5 and 6), and 1281L/R (lanes 7 and 8). Chromosomal DNA from mutant strains was used in PCRs visualized in the following lanes: 1 ( M5005_Spy_0680), 3 ( M5005_Spy_0785), 5 ( M5005_Spy_0830), and 7 ( M5005_Spy_1281). Lanes 2, 4, 6 and 8 contain the PCR products generated using wild-type strain M5005 chromosomal DNA. (B) Southern hybridization analysis of mutant strains. Chromosomal DNA isolated from wild-type strain M5005 (lane 1) and mutant strains M5005_Spy_0680 (lane 2), M5005_Spy_0785 (lane 3), M5005_Spy_0830 (lane 4), and M5005_Spy_1281 (lane 5) was digested with HincII and probed with a DNA fragment encoding the spc cassette. Each mutant strain contains only one insertion of the spc cassette.
|
3 h after the bacteria entered stationary phase. The total number of RNA samples analyzed was 100 (five strains times five biologic replicates times four time points). The bacteria were mixed with RNA Protect reagent (QIAGEN), and GAS cells were collected by centrifugation and stored at 80°C until processing. Labeled cDNA (targets) used for hybridization with the custom chip were prepared according to the procedure described by Affymetrix for use with Pseudomonas aeruginosa (1) with modifications described recently (8, 39). To minimize nonbiological variability caused by the multistep sample preparation and hybridization procedure, the targets were prepared according to a complicated experimental protocol designed by Partek Inc. and used in other studies (8, 39). All bacterial samples prepared on one day (5 strains x 4 time points = 20) were processed at the same time using 96-well plates for all steps of the protocol, from RNA isolation to labeling. All samples were placed near the center of the 96-well plate to minimize possible spurious variation induced by a plate-edge effect. Similarly, to avoid potential "right-to-left," "top-to-bottom," or "center-to-edge" spurious effects, the samples were arranged in a computer-generated matrix that took into account experimental parameters (strain and growth phase). This strategy minimizes the possibility that trends identified for changes in gene transcripts are confounded by unrelated parameters, such as culture batch or steps in target preparation. Each group of 20 samples prepared simultaneously was hybridized simultaneously with the custom GeneChips. Other conditions or parameters that in principle can detrimentally affect the quality of the resulting data (chip lot, wash batch, position, or scan order) were randomly assigned.
Chip hybridization data were acquired and normalized using Affymetrix GeneChip Operating Software (GCOS). Hybridization intensity values were normalized to the mean intensity of all GAS genes present on the chip using GCOS version 1.0, and only genes with a "present" signal were included for further analysis. Tables generated by GCOS containing normalized hybridization values obtained from all 100 chips were exported to Partek Pro 6.0 package (Partek Inc.) for statistical validation of data and subsequent analysis. Hybridization data were analyzed using a three-way analysis of variance (ANOVA) and linear contrasts to identify genes differentially expressed between each knockout strain and wild-type strain M5005. The blocking factors used for the ANOVA were strain, growth phase, and culture batch. Using this approach, factors such as wash batch, chip lot, position on plate, etc., can become covariate in the ANOVA to reduce noise. A PartekPro principal component analysis (PCA) and visualization system was used to assess chip quality and chip-to-chip variability. Input information combined hybridization intensity values and information about sample preparation and hybridization.
All expression microarray data have been deposited in the Gene Expression Omnibus (GEO) database at NCBI at http://www.ncbi.nlm.nih.gov/geo/.
TaqMan confirmation of microarray analysis.
To confirm that obtained microarray data represent true changes in gene expression in mutant strains, we validated the data with quantitative real-time PCR (TaqMan). From transcripts detected by microarray analysis as differentially expressed, we selected at random several (8 to 11 transcripts per mutant strain) of them for analysis. cDNA for TaqMan confirmation was generated from the RNA samples used for microarray analysis (two of five biological replicates) using the same procedure as for microarrays. TaqMan reactions for each gene tested were prepared in quadruplicate for each of two biological replicates and were amplified using Platinum Quantitative PCR SuperMix-UDG (Invitrogen). The transcript amounts in each mutant were standardized to an internal control gene (proS) and compared with standardized expression in the wild-type strain (
CT method).
Mouse infection experiments.
The two mouse models used for infection experiments have been described extensively in previous studies (18, 20, 21). GAS strains used for mouse infection studies were grown in THY medium to exponential phase (OD600
0.5), harvested, washed twice with ice-cold phosphate-buffered saline, and frozen at 70°C in aliquots. Thawed aliquots of strain M5005 and the isogenic mutant strains were adjusted to the same CFU/ml by diluting with phosphate-buffered saline. The number of CFU used for mouse inoculation was determined by plating on sheep blood agar. Male immunocompetent hairless mice (4 to 6 weeks old) of strain Crl:SKH1-hrBR (Charles River) were used for subcutaneous inoculation studies. The animals were anesthetized with isoflurane and inoculated with 2 x 107 CFU of GAS in 100 µl. To monitor abscess formation and infection progression, animal weight and abscess size were measured every day for the first week and twice a week thereafter for a total of 28 days. Abscess length (L) and width (W) values were used to calculate abscess area (A =
[L/2] x [W/2]) and volume (V = [4/3]
[L/2]2 x [W/2]) using the equation for a spherical ellipsoid (20). Outbred CD-1 Swiss male mice (4 to 6 weeks old) (Harlan) were used for intraperitoneal inoculation with a dose of 5 x 107 CFU in 250 µl. Near-mortality was monitored every 2 h for the first 48 h after infection and then every 6 h for the next 5 days.
|
|
|---|
There was no observable difference in colony phenotype or in vitro grow rate between the wild-type strain and the four mutant strains on blood agar or tryptose agar (data not shown). For example, there was no difference in colony size, zone of hemolysis, or mucoid appearance of the colonies. Similarly, there was no difference between the wild-type and mutant strains in bacterial chain length (data not shown). Consistent with these observations, the growth rate in THY medium of each mutant strain was identical to the wild-type parental strain (Fig. 3).
![]() View larger version (20K): [in a new window] |
FIG. 3. Growth profiles of isogenic mutant strains in THY medium. Arrows denote time points when bacteria were collected for expression microarray experiments. WT, wild-type strain M5005. Mutant strain numbers are shown in the insert for each figure.
|
Chip quality and chip-to-chip variability were assessed using a PartekPro PCA and visualization system. Importantly, the PCA discriminated very well between the data from the chips representing early logarithmic, mid-logarithmic, late logarithmic, and early stationary growth phases (Fig. 4). These results indicated that the transcriptome profile data are of high quality and demonstrated extensive remodeling of the transcript profile over time during in vitro culture.
![]() View larger version (63K): [in a new window] |
FIG. 4. PCA of data obtained from expression microarray analysis of the isogenic mutant strains. The PCA plot captures the variance in a data set in terms of principal components and displays the most significant of these on the x and y axes. The plot indicates that the transcriptome data are of high quality, as the samples cluster together according to growth phase: early logarithmic, red; mid-logarithmic, blue; late logarithmic, green; and stationary, purple.
|
![]() View larger version (21K): [in a new window] |
FIG. 5. Correlation of microarray and TaqMan assays. Several genes from each mutant strain were randomly chosen for data confirmation. Values obtained using both methods were plotted on the x and y axes separately for genes tested in each mutant. Each dot represents one gene tested by both methods.
|
M5005_Spy_0680,
M5005_Spy_0785, and
M5005_Spy_0830), expression of about 15% of all genes was significantly altered compared to the wild-type strain at one or more time points. In contrast, the transcriptome of mutant strain
M5005_Spy_1281 was far more altered compared to the wild-type strain, with 41% of transcripts detected after hybridization affected at one or more time points (Fig. 6). Changes in transcript levels correlated with growth phase in three of the four mutant strains (Fig. 7, left panel for each category; see Fig. S1 in the supplemental material). For example, in the
M5005_Spy_0680 mutant strain, the majority of affected transcripts occurred in the late logarithmic and stationary phases. Similarly, a massive change in transcript profile was observed in the
M5005_Spy_1281 mutant strain in the stationary phase, with over 70% of all changes occurring at that time point. Transcripts in the
M5005_Spy_0785 mutant strain were mostly affected in the early exponential and stationary phases. Finally, mutant strain
M5005_Spy_0830 differed prominently from the other three mutant strains, in that no correlation was observed between growth phase and altered transcript profile (Fig. 7).
![]() View larger version (79K): [in a new window] |
FIG. 6. Pairwise comparison of the number of genes affected in the four TCS mutant strains. The numbers and percentages presented are based on the number of genes detected after hybridization (n = 1,632).
|
![]() View larger version (53K): [in a new window] |
FIG. 7. Sixteen functional categories of genes affected in the four tested mutant strains. In each category panel, the four left bars of the graph represent the cumulative number of changed transcripts during the experiment in four TCS mutants. Each bar represents genes changed during early logarithmic (red), mid-logarithmic (blue), late logarithmic (green), and stationary (purple) growth phase. The fifth (black) bar represents the number of genes in each category annotated in the M5005 genome. Four yellow panels represent the percentage of changed genes in a particular category in relation to number of affected genes in the mutant strain; the last bar (black) represents the percentage of annotated genes in the M5005 genome.
|
M5005_Spy_0680 strain revealed that more genes were down-regulated (n = 206) than up-regulated (n = 136) at one or more of the time points analyzed. The transcript representing the M5005_Spy_0681 kinase was up-regulated in the mutant strain, suggesting that its cognate response regulator acts as a repressor of its own operon (Fig. 1B).
Similar to mutant strain
M5005_Spy_0680, more genes were down-regulated (n = 181) than up-regulated (n = 106) in the
M5005_Spy_0785 mutant strain. However, the transcript encoding the cognate sensor kinase was down-regulated, suggesting that it acts as an activator of its own operon (Fig. 1B). In addition, the transcripts representing an apparent operon (M5005_Spy_0780-0783) located directly upstream of this TCS were greatly down-regulated (
200-fold), suggesting that M5005_Spy_0785 also positively regulates expression of this operon encoding carbohydrate metabolism proteins (Fig. 1B).
Inactivation of the M5005_Spy_0830 response regulator affected the smallest number of GAS transcripts, with no bias toward down-regulated (n = 115) or up-regulated (n = 103) genes. In general, similar numbers of transcript changes occurred in all growth phases. The level of down-regulation was quite high in the strain with mutations of two contiguous genes (M5005_Spy_0832-0833) in an operon encoding a putative NAD-dependent malic enzyme (40-fold decrease) and malate-sodium symport (27-fold decrease). The transcript level of M5005_Spy_0831, encoding the sensor kinase component of this system, was decreased in the mutant strain, consistent with the idea that the response regulator is an activator of this TCS (Fig. 1B).
Transcript factor changes for genes in the
M5005_Spy_1281 strain varied from 8.1-fold (M5005_Spy_0149/SPy0175, encoding component of the phosphotransferase system [PTS]) to +5.4-fold for araD/SPy0179. The studied response regulator presumably acts as repressor of its own transcription, as the transcript encoding the sensor kinase was up-regulated (Fig. 1B).
Gene categories affected in mutant strains.
The annotated genes in the genome of strain M5005 were assigned to 16 functional categories. The percentage of annotated genes in each functional category was compared with the percentage of transcripts affected during at least one time point in the four analyzed mutant strains, regardless of growth phase (Fig. 7). This type of comparison can be useful because it can suggest a role for regulators of unknown function on the basis of over- or underrepresented gene categories. In general, for most gene categories, the percentage of affected genes in the mutant strains did not differ dramatically from the percentage of annotated genes in the wild-type strain. However, some differences were observed (Fig. 7, yellow panels in each functional category). For example, a prominent feature of the transcriptome of strain
M5005_Spy_0680 was the percentage of affected genes involved in carbohydrate and nucleotide metabolism. The number of affected transcripts was especially high in the transition from logarithmic to stationary phase (Fig. 7, left panels).
Mouse virulence of the four TCS mutant strains.
Inactivation of many TCSs in pathogenic bacteria alters virulence due to changes in gene expression, including but not limited to genes encoding virulence factors such as extracellular or cell wall-anchored proteins (5, 7, 15, 41). Inactivation of bacterial TCSs usually decreases virulence, but exceptions to this general observation have been reported (7, 28). Two well-described mouse models of infection were used to test the hypothesis that inactivation of these four GAS TCSs alters virulence. We first assessed the mouse virulence of the mutant strains by intraperitoneal inoculation. No significant virulence difference was observed between the wild-type parental strain and the four isogenic mutant strains as assessed by Kaplan-Meier survival curve analysis (data not shown). Next, the virulence of the four mutant strains relative to the wild-type strain was analyzed in a soft-tissue model of GAS infection. These experiments were conducted because it is possible that genes required for maximum virulence in this model differ from the intraperitoneal injection model. No significant difference in virulence was observed for three of the four mutant strains (
M5005_Spy_0785,
M5005_Spy_0830, and
M5005_Spy_1281). In contrast, mutant strain
M5005_Spy_0680 caused significantly larger abscesses than wild-type parental strain M5005, consistent with a hypervirulence phenotype in this model (Fig. 8).
![]() View larger version (50K): [in a new window] |
FIG. 8. Virulence difference between the wild type and the isogenic M5005_Spy_0680 mutant strain. (A) The abscess areas at day 2, day 6, and day 12 after subcutaneous infection of mice with wild-type strain M5005 and isogenic mutant derivative strain M5005_Spy_0680. Each dot represents the abscess area for an individual animal at the time point noted. Horizontal lines denote mean abscess area. Statistical significance was calculated with a t test. WT, wild-type strain M5005. (B) Abscesses developed by mice injected with the wild type (left) and the M5005_Spy_0680 mutant strain (right).
|
M5005_Spy_0680 hypervirulent mutant strain.
To gain initial insight into the possible molecular basis of the large-abscess phenotype of strain
M5005_Spy_0680, we compared the transcripts affected in this mutant strain with the wild-type parental organism, with special focus on genes potentially involved in virulence. The transcript level of several genes known to enhance virulence was significantly increased in the mutant strain, including scpA, sagGI, silD, and sic. In addition, the gene (prtS) encoding a recently described secreted protease that cleaves and inactivates interleukin-8 (4) was up-regulated in the early logarithmic phase of growth. The largest category of genes with altered transcripts was genes of unknown function, with the majority of changes observed in the stationary and/or late logarithmic growth phases. Although it is not known whether any of these genes encode virulence factors, we note that four of them (M5005_Spy_0354, M5005_Spy_0355, M5005_Spy_0666, and M5005_Spy_0721) have typical or probable gram-positive secretion signal sequences. Importantly, these four transcripts were all increased, consistent with an enhanced-virulence phenotype. As noted above, inactivation of the M5005_Spy_0680 response regulator also affected the expression of multiple genes apparently involved in bacterial metabolism. For example, decreased expression of M5005_Spy_0025-0031 genes (purMNDEK) and M5005_Spy_0347-0349 genes (nrdFIE), involved in nucleotide metabolism, and genes involved in arginine utilization (M5005_Spy_1237-1238/SP1506-1507) were observed. In addition, carbohydrate utilization genes were prominently affected, including several PTS transport systems and genes apparently involved in maltose (M5005_Spy_1055-1067/SP1291-1304) and tagatose (M5005_Spy_1636/SP1921) utilization pathways. These and other findings are summarized in Fig. 9.
![]() View larger version (51K): [in a new window] |
FIG. 9. Inactivation of the M5005_Spy_0680 gene affects levels of multiple transcripts in GAS. M5005_Spy_0680/0681 TCS responds to unknown signals and controls (directly or indirectly) the transcription of genes involved in virulence and multiple cellular and metabolic processes.
|
|
|
|---|
Analysis of the in vitro growth phenotype of mutant strains can yield useful data, but studies of this type do not address the role of specific genes in pathogen interaction with an intact host. Thus, we used two well-characterized mouse models to test the hypothesis that inactivation of each of these four TCS genes altered virulence. The
M5005_Spy_0680 mutant strain caused significantly larger skin lesions in mice than the wild-type parental strain M5005, whereas the other mutants were not altered in mouse virulence. Although we did not address the molecular basis of increased virulence of the
M5005_Spy_0680 mutant strain, we identified clues that may help to explain this phenotype. Expression of 310 genes was directly or indirectly influenced by this TCS during at least one of the time points analyzed, including several (scpA, sagGI, silD, and sic) known to participate in GAS pathogenesis. In addition, we note that the prtS gene encoding a protease that inactivates interleukin-8 (4) was up-regulated in the early logarithmic phase of growth, a process that could contribute to immune evasion in vivo. Expression also was increased for 24 hypothetical proteins with no homologs in other bacteria. None of these 24 up-regulated genes encoded proteins with an LPXTG cell wall-anchoring motif, a common feature of many virulence factors made by gram-positive bacteria. However, 5 of these 24 proteins have putative membrane-spanning motifs, raising the possibility that they are involved in transport of molecules (including virulence factors). In addition, four up-regulated proteins have gram-positive secretion signal sequences and thus may be located extracellularly. We speculate that study of the involvement of these proteins in host-pathogen interaction might be fruitful.
Global analysis of in vitro transcription.
Inactivation of each of the four response regulators affected expression of a very large number of genes, consistent with previous transcriptome studies of CovR/CovS and Ihk/Irr (7, 40). For example, the percentage of affected genes during at least one experimental time point varied from 12% to 41% of genes, as observed for the
M5005_Spy_1281 mutant strain. In addition to the large number of affected genes, changes in transcription profiles were clearly growth phase dependent. Most genes were affected at only one or two of the time points analyzed, suggesting that each TCS acts primarily at a distinct phase of growth. Virtually all transcription changes were relatively modest, rarely exceeding a twofold difference (see Fig. S1 in the supplemental material). Taken together, these two observations suggest that most of the changes in gene expression were indirect effects mediated by secondary effectors rather than primary effects caused by direct interaction of the response regulator with many different promoter regions. It is also possible that some of the changes in transcript level were due to cross talk (a sensor kinase phosphorylates a noncognate response regulator) or cross-regulation (a sensor kinase phosphorylates more than one response regulator). Although cross talk and cross-regulation have been reported for E. coli and Pseudomonas spp., they have not been studied in GAS (10, 30, 38, 42).
The transcriptome data permit us to identify common themes regarding genes influenced by each of the four TCSs studied, thereby providing potential insight about functions of each TCS. The M5005_Spy_0680/0681 TCS apparently affects nucleotide and complex carbohydrate metabolism, especially during transition from logarithmic to stationary phase. Although speculative, it is possible that the observed down-regulation of the ccpA homolog in late exponential phase participates. The CcpA (carbon catabolite repression) protein is a major regulator of sugar transport and metabolism in gram-positive bacteria (43). Finally, we note that Virtaneva et al. (39) recently reported a significant correlation between expression of these genes and clinical pharyngitis in experimentally infected monkeys. This observation supports the idea that the M5005_Spy_0680/0681 TCS plays an important role in host-pathogen interaction.
No distinct global theme was revealed by analysis of the transcripts influenced by the M5005_Spy_0784/0785 TCS. However, one apparent transcriptional unit may be directly regulated by this TCS. Four genes (M5005_Spy_0780-0783) encoding a putative mannose/fructose-PTS system were down-regulated over 200-fold in the late logarithmic and stationary growth phases in the mutant strain. These genes are located directly upstream of the two genes encoding this TCS (Fig. 1B). This fact, taken together with the extremely high level of transcript down-regulation observed in the mutant strain, suggests that it is a direct activator of this operon.
Two observations suggest that the M5005_Spy_0830/0831 TCS is involved in malate metabolism. First, two genes (M5005_Spy_0832-0833) located directly downstream of this TCS were highly down-regulated in the mutant strain (see Table S1 in the supplemental material). These genes encode homologs of a malate-sodium symporter (M5005_Spy_0832) and an NAD-dependent malic enzyme (M5005_Spy_0833). Second, consistent with these findings, a distant Bacillus subtilis homolog of this TCS is involved in malate transport, inducible by external malate (37).
The most prominent finding related to the M5005_Spy_1280/1281 TCS was the very large number of altered gene transcripts observed in the stationary phase. Only 16 transcripts were altered in the early logarithmic phase, whereas 473 were changed in the stationary phase, suggesting that this TCS functions as a major "switch" to change cell physiology late in growth. We find it remarkable that such a massive transcriptome change did not produce a significant difference in mouse virulence, an observation stressing the complexity of GAS pathogenesis.
Most expression microarray studies have concentrated on analysis of single TCS mutants at one growth point rather than multiple strains and different growth phases (8, 22, 25, 29, 41). Our work adds to a growing list of reports describing expression array analysis of multiple TCSs in bacteria (14, 26, 27). One theme emerging from these investigations is that great variance exists in the number of genes directly or indirectly regulated by each TCS. For example, microarray analysis of mutations in all 32 E. coli TCSs identified three classes of transcriptional profiles based on the number of affected genes (minimal, moderate, or dramatic changes) (26). Our results mirror these findings. We found that three of the four GAS TCS mutants had moderate changes in transcriptome, whereas the transcript profile of
M5005_Spy_1281 was dramatically changed. Given the very large number of transcripts altered in the
M5005_Spy_1281 mutant strain, it is possible that it is a central actor in interfacing with multiple downstream regulators. Thus, inactivation of M5005_Spy_1281 would result in significant dysregulation of a large number of downstream genes. Consistent with this, expression of a very large number of transcription regulator genes was altered in the mutant strain (see Table S1 in the supplemental material).
Supplemental material for this article may be found at http://iai.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»