Transcriptional profiling of Bordetella pertussis reveals requirement of RNA chaperone Hfq for Type III secretion system functionality.

Bordetella pertussis, the causative agent of human whooping cough (pertussis) produces a complex array of virulence factors in order to establish efficient infection in the host. The RNA chaperone Hfq and small regulatory RNAs are key players in posttranscriptional regulation in bacteria and have been shown to play an essential role in virulence of a broad spectrum of bacterial pathogens. This study represents the first attempt to characterize the Hfq regulon of the human pathogen B. pertussis under laboratory conditions as well as upon passage in the host and indicates that loss of Hfq has a profound effect on gene expression in B. pertussis. Comparative transcriptional profiling revealed that Hfq is required for expression of several virulence factors in B. pertussis cells including the Type III secretion system (T3SS). In striking contrast to the wt strain, T3SS did not become operational in the hfq mutant passaged either through mice or macrophages thereby proving that Hfq is required for the functionality of the B. pertussis T3SS. Likewise, expression of virulence factors vag8 and tcfA encoding autotransporter and tracheal colonization factor, respectively, was strongly reduced in the hfq mutant. Importantly, for the first time we demonstrate that B. pertussis T3SS can be activated upon contact with macrophage cells in vitro.

ATP-binding cassette protein

colony forming unit

outer membrane protein

P-value

quantitative reverse transcription polymerase chain reaction

Steiner-Scholte medium

Type III secretion system

wild-type

Introduction

Bordetella pertussis is a Gram-negative bacterial pathogen of the human respiratory tract and causative agent of whooping cough (pertussis). Infection by B. pertussis is especially severe in infants and despite vaccination, B. pertussis remains a major cause of children mortality worldwide, predominantly in developing countries. However, pertussis is now on the rise also in developed countries and is reemerging in some highly vaccinated populations. The increase in pertussis cases has been mostly attributed to short-lived immunity induced by the acellular vaccine and escape from immunity due to pathogen adaptation and antigen variation. Therefore, there is an urgent need for a better understanding of molecular mechanisms underlying the pathogenesis of B. pertussis infection that would lead to identification of novel factors as potential components of new vaccine formulations.

In order to successfully infect the host, B. pertussis produces a broad spectrum of virulence factors including filamentous hemagglutinin, pertactin, fimbriae, pertussis toxin and adenylate cyclase toxin. The type III secretion system (T3SS) allowing injection of proteins referred to as effectors directly into the cytosol of eukaryotic host cells has also been identified in Bordetella subspecies. The T3SS is a supramolecular structure consisting of an inner membrane export apparatus, a basal body that spans from inner to the outer membrane, an extracellular needle and the needle tip (translocon) complex. In the closely related B. bronchiseptica, this “injectisome” is composed of approximately 30 proteins including chaperones and effectors which are encoded by the bsc locus. Bordetella T3SS represents an important virulence factor as it subverts innate and adaptive immune responses of the host and is required for long term colonization of trachea of rats and mice. Although the T3SS of B. pertussis is highly homologous to that of B. bronchiseptica and bsc genes are actively transcribed, laboratory-adapted high-passage B. pertussis strains like Tohama I do not produce and secrete T3SS substrates when grown in vitro. However, the functionality of T3SS can be reverted upon contact with the host as evidenced by production of the T3SS tip complex protein Bsp22 in B. pertussis cells recovered from infected mice. Intriguingly, the T3SS activation is reversible as the production of Bsp22 protein is lost again after successive in vitro passages. Bsp22 is the most abundant T3SS substrate secreted by B. bronchiseptica and is considered as a hallmark of T3SS functionality. Deletion of bsp22 gene in B. bronchiseptica disrupts T3SS-mediated cytotoxicity and affects colonization levels and persistence in the lower respiratory tract of mice. Interestingly, production of Bsp22 was observed also in fresh clinical isolates of B. pertussis, thereby further confirming that contact with the host is crucial for switching on the expression of T3SS.

The expression of the majority of B. pertussis virulence factors is controlled at transcriptional level by a two-component system encoded by the bvg locus (for a review see ref.). This consists of the transmembrane sensor kinase BvgS and of the DNA-binding response regulator BvgA which in its phosphorylated form binds to promoter regions and activates transcription of dependent virulence genes. Expression of the bsc locus is regulated at both transcriptional and posttranscriptional levels by the BvgAS system and by factors encoded by the btr locus which is adjacent to bsc genes (). The BtrS protein whose production is tightly activated by the BvgAS system bears homology to the family of extracytoplasmic function sigma factors and is required for transcription of bsc and btr loci. BtrU, BtrV and BtrW factors carry high similarity to a partner switching module known in Gram-positive bacteria and are required for secretion and stability of T3SS substrates.

Figure 1.

Comparative transcriptomic analysis of laboratory-adapted wt and Δhfq strains of B. pertussis. (A) Genomic organization and open reading frames of the Type III secretion cluster in B. pertussis. The genes comprising both bsc and regulatory btr loci are presented as arrows. The drawing was created using CLC Sequence viewer version 6.6.2 (CLC bio A/S). (B) Functional categorization of genes significantly affected by loss of Hfq function. Data are expressed as numbers of genes within category which were significantly down- (green bars) or upregulated (red bars) in the hfq mutant cells grown to exponential (E) or stationary (S) phase of growth. Total numbers of B. pertussis genes associated with a specific category are shown in brackets. Significantly overrepresented functional categories are depicted with asterisks, P ˂ 0.05. (C) Changes in transcript abundance patterns of selected virulence factors in B. pertussis. The heat map presents comparison of gene expression profiles between Δhfq and wt cells (log2FC Δhfq/wt) grown to exponential (E) or stationary (S) phase as determined by microarray analysis. Data are averaged from 2 biological replicates, each gene is represented by 2 array elements. All expression profiles are in rows and are represented using the color scale on the top. (D) RT-qPCR analysis of bsp22, bteA, vag8 and brkA transcript levels in wt and Δhfq cells. Total RNA was isolated from cells grown to exponential (E; white columns) and stationary (S; gray columns) phases and RT-qPCR was performed using primers specific for bsp22, bteA, vag8 and brkA genes. Target expression levels in wt and Δhfq cells were normalized to the reference gene rpoB and the wt/Δhfq expression ratios were calculated. The values are means plus standard deviations from 3 experiments. (E) Western blot analysis of Vag8 and BrkA protein levels in Tohama I wt and Δhfq strain. Both strains were grown in SS medium and cells were harvested at exponential (E) and stationary (S) phase. Samples of cell lysates equivalent to 0.1 OD600 unit of initial culture were separated on 10% SDS-PAGE gel and analyzed by immunoblotting using anti-Vag8 and anti-BrkA antibodies. Relevant parts of the membranes are shown. The result is representative of 3 experiments.

In many pathogenic bacteria the RNA chaperone Hfq and small non-coding regulatory RNAs (sRNAs) emerged as key players in posttranscriptional regulation of virulence and physiological fitness. Posttranscriptional regulation of B. pertussis virulence has not yet been studied extensively, however, several sRNA were identified in B. pertussis and recently, we have shown that the RNA chaperone Hfq is required for virulence of B. pertussis. The Tohama I Δhfq strain produced decreased amounts of adenylate cyclase toxin and secreted reduced amounts of pertussis toxin. Consequently, the hfq mutant was clearly attenuated in the mouse respiratory model of infection as its lethality as well as its capacity to colonize mouse lungs was strongly reduced when compared to the wt strain. These data encouraged us to apply microarray analysis in order to evaluate the global impact of hfq deletion on gene expression with specific focus on virulence. Through differential analysis of expression profiles of laboratory-adapted and mouse-passaged wt and Δhfq strains we aimed at the identification of Hfq regulon in B. pertussis. In this work we show that loss of Hfq function has a deep impact on expression of several virulence factors including T3SS. Our data clearly indicate that Hfq is strongly required for functionality of T3SS in B. pertussis.

Results

Hfq affects global gene expression in B. pertussis

To evaluate the effect of hfq deletion on global gene expression profiles, DNA microarrays covering 3554 B. pertussis protein-coding genes were hybridized with total RNA purified from wt and Δhfq cells grown to exponential and stationary phase of growth. A cutoff of at least twofold variation in expression (log2FC >1; <−1; P < 0.05) was applied to identify genes with significant changes in transcript abundance between wt strain and hfq mutant. In the exponential phase, 368 genes (≈10.6%) fulfilled these criteria and displayed differences in transcript abundance ( and Table S1). Among them, 212 genes were downregulated and 156 upregulated in the Δhfq strain. In the stationary phase, the transcriptional profiles exhibited similar pattern with 351 (10.1%) significantly affected genes (Table S2). Among them, 186 genes were downregulated and 165 upregulated. Functional class categorization of filtered Hfq-dependent genes showed significant enrichment (P < 0.05) in some gene categories (). In the exponential phase, “Energy metabolism” was the only functional category significantly enriched for Hfq-dependent genes (20.9%; P < 0.002) since the “Pathogenicity and Adaptation” category did not reach the statistical significance (14.0%; P = 0.054). In the stationary phase, genes related to “Pathogenicity and Adaptation” (13.6%; P < 0.05) and “Information transfer” (15.1%; P < 0.015) categories were significantly overrepresented. Notably, in the exponential and stationary phase 80% (28/35) and 88% (30/34) of genes related to “Pathogenicity and adaptation” category were downregulated in the hfq mutant, respectively, being fully in accordance with reduced virulence of the mutant. Interestingly, a large portion of affected genes (17/25 and 20/36, respectively) within “Information transfer” category belonged to genes encoding 30S and 50S ribosomal proteins. While not significantly enriched, the highest numbers (125 and 110, respectively) of affected genes were related to largest functional category of “Surface proteins.” Besides large number of genes predicted to encode putative membrane proteins this group also included several genes encoding ATP binding cassette (ABC) transport proteins which expression in the exponential phase was largely upregulated in the mutant (Tables S1 and S2). Interestingly, 16 of 19 genes (84.2%) associated with regulatory function were downregulated in the Δhfq strain in the exponential phase.

Lack of Hfq results in reduced expression of several virulence genes including T3SS

Within the set of genes related to pathogenicity and displaying significantly lower transcript levels in the mutant regardless of the growth phase, we identified the T3SS bsc gene cluster (). Among the significantly downregulated genes were bsp22, bopD and bopN encoding the T3SS tip complex protein, the pore-forming translocator protein and the T3SS effector, respectively. In contrast, the expression of the btr locus which comprises T3SS regulatory genes was not significantly affected. Notably, expression of bteA and btcA genes which are not part of the bsc locus and encode another T3SS effector protein BteA and its chaperone BtcA, respectively, was also significantly downregulated in the hfq mutant. Similarly, expression of other virulence genes including vag8 autotransporter, serum resistance gene brkA (only in the exponential phase) and tracheal colonization factor tcfA was also significantly decreased in the mutant strain, suggesting that Hfq has an impact on a broad spectrum of virulence genes (). In contrast to the majority of virulence genes and in agreement with our previous data, the ptx/ptl locus was significantly upregulated in the hfq mutant (Tables S1 and S2).

In order to verify our microarray data, we have analyzed transcript levels of several virulence genes including the T3SS locus. We have quantified the levels of the bsp22, bteA, vag8 and brkA transcripts in both wt and Δhfq strains of B. pertussis by RT-qPCR. As shown in , the relative expression levels of bsp22, bteA, vag8 and brkA genes were decreased in the hfq mutant thereby confirming our microarray data. Furthermore, Western blot analysis of wt and Δhfq cell lysates showed that levels of Vag8 and BrkA proteins mirrored the relative transcript levels and thereby provided another layer of evidence for credibility of our microarray data ().

Hfq is required for regained production of Bsp22 protein upon passage in mice

Transcriptional profiling of wt and Δhfq strains adapted to laboratory culture conditions suggested that Hfq is required for expression of several virulence factors including T3SS. Therefore, we asked whether the differential expression of the T3SS locus would have an impact on production of T3SS components. Laboratory-adapted B. pertussis strains like Tohama I do not produce T3SS substrates under standard in vitro conditions, however, Gaillard and coworkers have shown that after passage through mice the recovered Tohama I cells produced Bsp22 into SS medium and kept the production for 3 additional passages in vitro. Therefore, we decided to passage both strains in mice and monitor their ability to produce and secrete the T3SS tip complex protein Bsp22 (marker of T3SS functionality).Toward this aim, CD-1 mice were infected with sublethal doses of both wt and Δhfq strain (1 × 105 CFU/mouse) and colonization profiles of both strains were examined over 23 days. As illustrated by and in line with our previous observations, the wt strain was able to proliferate in lungs and its CFU counts increased by more than one order of magnitude within 8 days of infection whereas the hfq mutant was unable to multiply and was cleared from mouse lungs with a significantly higher efficacy. In the course of the infection, the viable wt and Δhfq bacteria were recovered from mouse lungs at each experimental time point, subcultivated in SS medium and assayed for Bsp22 accumulation in whole-cell lysates and in culture supernatants. As shown in , Bsp22 protein accumulated in whole cell lysates as well as in culture supernatants of mouse-passaged wt bacteria recovered at days 8, 12 and 23 post infection. Strikingly, Bsp22 protein was absent in cell lysates and culture supernatants of passaged Δhfq bacteria. These results indicated that Hfq is required for both expression and production of Bsp22 and consequently, for T3SS function in B. pertussis.

Figure 2.

Hfq is required for Bsp22 and Vag8 production in B. pertussis. (A) Colonization of mouse lungs after intranasal challenge with the wt or Δhfq strain. 15 mice/strain were infected intranasally with 1 × 105 CFU/mouse and groups of 3 animals were sacrificed 2 h (day 0) and 5, 8, 12 and 23 days after infection. Lung homogenates were plated in duplicates on BGA plates and CFU of the wt (•) and Δhfq (○) strains were counted. The values are means plus standard deviations. *, P˂ 0.05; **, P˂ 0.01 (Δhfq versus wt). The dashed line represents the limit of detection. (B) Westernblot analysis of Bsp22 and Vag8 protein levels in the wt or Δhfq strains recovered from lung homogenates in the course of colonization experiment (see panel A). Samples equivalent to 0.1 OD600 unit (pellets, P) or 1 OD600 unit (supernatants, S) of initial culture were analyzed by immunoblotting using anti-Bsp22 and anti-Vag8 antibodies. Relevant parts of the membranes are shown.

Transcriptional profiling of passaged strains confirms requirement of Hfq for virulence

Next, we assumed that transcriptional profiling of wt and Δhfq bacteria passaged in mice would not only allow us to further corroborate the effect of hfq deletion on production of T3SS components but could also enlighten the global contribution of Hfq to B. pertussis adaptation and survival during infection. Therefore, in order to elaborate this assumption, total RNA was purified from wt and Δhfq bacteria recovered at days 0 and 12 post infection and used for microarray analysis to determine changes in gene expression occurring upon passage in mice. A cutoff of at least twofold variation in expression (log2FC >1; <−1; P < 0.05) between wt and hfq mutant strains (Δhfq versus wt at day 0 and at day 12) or in the course of infection (day 12 vs. day 0 for both strains) was applied to identify genes with significant changes in transcript abundance. Comparison of expression profiles of both strains at day 0 revealed 346 (10.1%) significantly affected genes (, left histogram) showing a distribution into functional categories similar to our initial experiment (). In the hfq mutant, 198 genes were downregulated and 148 genes were upregulated when compared to the wt strain (Table S3). None of the functional categories was enriched significantly, however, categories “Energy metabolism” (14.7%; P = 0.062), “Information transfer” (13%; P = 0.074) and “Pathogenicity and Adaptation” (13.2%; P = 0.058) were highly represented. Within the latter category, expression of majority of the bsc locus and bteA and vag8 genes was significantly diminished in the hfq mutant. The difference in expression of these virulence factors became even more apparent when we compared the transcriptomic profiles of wt and Δhfq bacteria recovered 12 days post infection (Table S4). As indicated in (right histogram), the total number of affected genes (359, 10.4%) did not surprisingly differ from day 0, however, the “Pathogenicity and Adaptation” category became the only significantly overrepresented functional category (18.4%; P < 1 × 10−4) with 40 out of 46 genes downregulated in the hfq mutant. Likewise, 12 days post infection the difference in expression of vag8 and T3SS-specific genes between wt and Δhfq strains further increased (). Interestingly, expression of the btrU gene which plays a role in secretion of T3SS substrates was also significantly downregulated in the hfq mutant upon infection. In addition, genes involved in iron metabolism such as bhuU and BP0546 encoding hemin permease and bacterioferritin, respectively, also exhibited lower transcript levels in the Δhfq strain at day 12 (Table S4). Within the set of genes observed to have significantly higher transcript levels in the hfq mutant 12 days post infection were bvgR and fim2 encoding BvgR regulatory factor and serotype 2 fimbrial subunit, respectively.

Figure 3.

Comparative transcriptomic analysis of mouse-passaged wt and Δhfq strains. (A) Functional categorization of genes significantly affected by loss of Hfq function in B. pertussis cells recovered at day 0 (D0) and at day 12 (D12) post infection. Data are expressed as numbers of genes within category which were significantly down- (green bars) or upregulated (red bars) in the hfq mutant. Total numbers of B. pertussis genes associated with a specific category are shown in brackets. Significantly overrepresented functional categories are depicted with asterisks, P ˂ 0.05. (B) Changes in transcript abundance patterns of selected virulence factors in B. pertussis. The heat map presents comparison of gene expression profiles between Δhfq and wt cells recovered (log2FC Δhfq/wt) at day 0 and day 12 post infection as determined by microarray analysis. All expression profiles are in rows and are represented using the color scale on the top, each gene is represented by 2 array elements. (C) RT-qPCR analysis of bsp22, bteA, bopN and vag8 transcript levels in wt and Δhfq bacteria upon infection. Differential expression analysis was performed using total RNA isolated from wt and Δhfq bacteria recovered from infected mice at day 0 (white columns) and day 12 (gray columns). The values are means plus standard deviations from 3 replicates.

When compared to wt strain, relative abundance of bsp22, bteA, bopN and vag8 transcripts as determined by RT-qPCR () was highly decreased in the Δhfq strain and this trend was further intensified in samples collected at day 12 post infection, thereby confirming the output of microarray data. In addition, Vag8 levels were determined in cell lysates of both passaged strains by Westernblot analysis. As seen in and in agreement with microarray and RT-qPCR data, the hfq mutant produced significantly reduced amounts of Vag8 protein while the Vag8 production in wt strain was further increased upon infection.

In order to characterize the infection-associated changes in gene expression in individual strains we have also compared transcriptional profiles in RNA samples collected at day 12 and day 0 in both wt and Δhfq strains by microarray analysis using the same threshold parameters (Fig. S1A). Comparison of transcript levels determined at day 12 with those obtained at day 0 revealed 115 and 78 genes to be affected in the wt and Δhfq strain, respectively (Fig. S1A, Tables S5 and S6). This analysis showed that upon contact with the host the T3SS-specific genes as well as the vag8 gene are significantly upregulated in the wt strain (Fig. S1B). This pattern of expression was observed also in the hfq mutant although to a much lower extent. Notably, expression of the btrS gene encoding T3SS regulatory factor was upregulated in both strains while transcript levels of btrU were increased only in the wt strain. Relative transcript levels of selected genes as determined by RT-qPCR were in agreement with microarray data. In contrast to Δhfq strain, the expression levels of bsp22, bteA, bopN and vag8 genes assessed in wt bacteria recovered at day 12 were increased when compared to those recovered at day 0 (Fig. S1C).

Hfq is required for expression of T3SS in B. pertussis cells recovered from infected macrophages

Our data indicated that Hfq is required for T3SS activation in mouse-passaged bacteria. During infection, resident and infiltrating macrophages are one of the first cells of the innate immune system to sense B. pertussis cells and thus play a very important role in the clearance of the pathogen. Therefore, we asked if the contact with macrophages in vitro would be sufficient for activation of T3SS in B. pertussis cells. Thus, in order to corroborate this hypothesis, we infected RAW 264.7 macrophages with wt and Δhfq bacteria and viable bacteria were recovered 24 and 48 hours post infection. Control bacteria which were not in contact with macrophages (time 0) together with recovered bacteria were further subcultivated in SS medium and culture supernatants were tested for presence of the Bsp22 protein. As shown in , wt cells recovered at both time points produced Bsp22 while the Δhfq cells did not. This result indicated that upon contact with macrophage cells the functionality of T3SS (assessed as production of Bsp22 protein) can be regained only in wt bacteria and that similarly to mouse-passaged bacteria, the production of Bsp22 can be observed in subsequent in vitro cultures of recovered bacteria. In order to further refine our observation, we have repeated macrophage infection experiment only with the wt strain and this time the viable cells were recovered already 3, 6 and 12 hours post infection. Control bacteria (time 0) and recovered bacteria were cultivated in SS medium and samples for RNA as well as for Westernblot analysis were harvested at each time point. As presented in , wt cells which were recovered from macrophages produced and secreted Bsp22 protein already 3 hours upon infection indicating that relatively short-time interaction with immune cells was sufficient for activation of T3SS. In addition, quantitative PCR analysis of RNA samples () revealed that upon infection of macrophages, expression of bopN, bsp22 and bteA genes was increased and thus copied the expression patterns observed in mouse-passaged strains.

Figure 4.

Expression of T3SS is activated only in wt B. pertussis cells upon infection of macrophages. (A) Westernblot analysis of Bsp22 protein levels in culture supernatants of wt or Δhfq bacteria recovered from RAW 264.7 macrophages. Bacteria recovered from macrophages 24 and 48 hours post infection were used to inoculate liquid cultures. Samples of supernatants equivalent to 1 OD600 unit of initial culture were analyzed by immunoblotting using anti-Bsp22 antibody. Relevant part of the membrane is shown. The result is representative of 3 experiments. (B) Westernblot analysis of Bsp22 protein levels in culture supernatants of wt bacteria recovered from RAW 264.7 macrophages. Viable bacteria, recovered from macrophages 3, 6 and 12 hours post infection, were used to inoculate liquid cultures. Samples of supernatants equivalent to 1 OD600 unit of initial culture were analyzed by immunoblotting using anti-Bsp22 antibodies. Relevant part of the membrane is shown. The result is representative of 3 experiments. (C) RT-qPCR analysis of bopN, bsp22 and bteA transcripts levels in wt bacteria upon infection of macrophages. Differential expression analysis was performed using total RNA isolated from viable wt bacteria recovered from macrophages 3 (white columns), 6 (gray columns) and 12 hours (black columns) post infection. The values are means plus standard deviations from 3 replicates.

Discussion

Recently we demonstrated that the hfq mutant of B. pertussis is significantly attenuated in virulence and explained this observation by reduced production and/or secretion of the key virulence factors adenylate cyclase and pertussis toxins by the mutant strain. Nevertheless, we realized that additional factors may be responsible for the reduced virulence of the mutant. In the present study, we therefore aimed at full characterization of the Hfq regulon in B. pertussis with focus on factors involved in pathogenicity. Transcriptional profiling has shown that Hfq is a global regulator in B. pertussis as the deletion of hfq gene had a large effect on gene expression with more than 10% of significantly affected genes in all experiments. Among the most overrepresented functional categories were “Information transfer,” “Energy metabolism” and “Pathogenicity and Adaptation” categories. Within the last category, comparative transcriptomic analysis identified additional virulence factors that play important roles in pathogenicity of B. pertussis and whose expression was significantly downregulated in the Δhfq strain. Besides tracheal colonization factor tcfA, serum resistance factor brkA and autotransporter vag8, a major part of the T3SS locus and T3SS effector gene bteA were the most prominent virulence factors downregulated in the hfq mutant. Quantitative PCR analysis confirmed our transcriptomic data and provided further evidence for Hfq involvement in T3SS expression. T3SS is one of the most remarkable results of long-standing association of bacteria with eukaryotic hosts and Hfq was reported to play specific roles in T3SS activity in several pathogens. Therefore, in our next experiments we concentrated on elucidation of the Hfq requirement for T3SS expression. Since the expression and production of T3SS components such as Bsp22 in laboratory-adapted strains is switched on only when the cells come into contact with the host, both wt and Δhfq strains were passaged through mice. Comparative analysis of transcriptomic profiles of Tohama I wt and Δhfq cells recovered from mouse lung revealed that expression of several T3SS-specific genes including the bteA effector was significantly induced in the wt strain, while this effect was much less obvious in the mutant strain. Hfq is known to act at the posttranscriptional level and therefore, we speculated that Hfq may exert its effect rather indirectly via T3SS-specific regulator BtrS encoded within the btr locus. However, we did not observe any substantial difference in btrS gene expression between wt and Δhfq strains. In fact, when we compared gene expression profiles in bacteria recovered at day 12 and day 0, the btrS expression was similarly increased in both wt (log2FC = 1.3; P < 1 × 10−6) and Δhfq strains (log2FC = 0.9, P < 2× 10−5, data not shown) suggesting that Hfq does not affect btrS expression and may employ another mechanism for T3SS activation. Remarkably, our microarray data indicated that several putative transcriptional factors are expressed at significantly reduced levels in the hfq mutant. Although none of them belongs to the AraC family of transcriptional activators which were shown to activate T3SS in Yersinia, it would be of interest to determine their possible role in regulation of T3SS activity in B. pertussis. In this study we confirmed the observation that a) in the laboratory-adapted Tohama I strain, T3SS genes are actively transcribed in vitro and b) while T3SS substrates such as Bsp22 protein are not produced and secreted under in vitro growing conditions, this ability can be regained upon passage in the host. Furthermore, infection experiment showed that differential expression of the bsp22 gene translated into strikingly altered levels of Bsp22 protein, a hallmark of T3SS functionality in B. pertussis. Indeed, Westernblot analysis of bacteria recovered from infected mice showed that in remarkable contrast to the wt strain, Bsp22 protein did not accumulate in both cell lysates and culture supernatants of the hfq mutant. The absence of Bsp22 protein in spite of reduced but detectable bsp22 transcription indicated that Hfq affects T3SS functionality also at posttranscriptional level and consequently that Hfq is required for synthesis of T3SS substrates. Recently, Hfq was reported to activate synthesis of T3SS proteins via a small RNA in Y. pestis and it would be of high importance to search for an analogous sRNA in B. pertussis. It should be mentioned that in contrast to the Δhfq strain, the btrU gene was found upregulated in the wt strain upon infection, however, this factor was shown to be required for secretion but not production of T3SS substrates. Therefore, its reduced expression cannot explain undetectable accumulation of Bsp22 protein inside the Δhfq cells. The precise mechanism of B. pertussis T3SS activation upon contact with the host remains unclear, however, the rapid activation of Bsp22 production upon infection and persistence of this effect for several in vitro passages suggests an epigenetic regulation.

Among other virulence factors which were significantly affected upon passage in mice, expression of tracheal colonization factor tcfA was significantly reduced in the Δhfq strain when compared to the wt strain. Recently, we reported that the hfq mutant of B. pertussis exhibits a striking colonization defect when used in mixed infection experiment with the wt strain. We understood this result primarily as a consequence of reduced production of adenylate cyclase toxin by the hfq mutant, however, herein presented data indicate that impaired expression of T3SS gene cluster and tracheal colonization factor tcfA may as well significantly contribute to reduced capacity of the hfq mutant to colonize mouse lungs. In support, T3SS-deficient strains of B. bronchiseptica and B. pertussis were cleared much faster from the respiratory tract of infected animals. Remarkably, we could recapitulate our data obtained in in vivo infection experiments also in in vitro infection model of murine macrophages. In contrast to Δhfq cells, we observed increased expression of T3SS genes and regained ability to produce Bsp22 protein in wt cells recovered from RAW 264.7 macrophages. These data represent the first evidence that T3SS can be activated also upon contact with immune cells in vitro and therefore, we consider macrophage cells as a very versatile and cost-effective infection model for gene expression studies in B. pertussis.

Notably, the expression of autotransporter vag8 gene as determined by microarray and RT-qPCR analyses was also strongly decreased in the hfq mutant and almost mirrored T3SS expression patterns under all tested conditions. Likewise, as evidenced by immunoblotting, the production of Vag8 protein was significantly reduced in the hfq mutant under standard laboratory conditions. Furthermore, in contrast to Δhfq cells, the expression and production of Vag8 was further increased in wt cells upon infection and thereby again copied the T3SS expression pattern. These results suggested that besides being regulated by the BvgAS system the expression of non-adjacent bsc and vag8 genes is coordinated by some yet unknown trans-encoded factor. Vag8 plays an important role during infection as it interacts with C1 esterase inhibitor and confers to serum resistance. In addition, Vag8 is presumed but not shown to be an accessory factor of T3SS and our data would support this hypothesis. In contrast to other virulence factors, the transcription of the ptx/ptl locus, encoding pertussis toxin and its secretion apparatus, was significantly upregulated in the mutant. This observation is fully in line with our previous study where we have shown that in spite of increased transcription of ptx/ptl locus the amount of toxin secreted by the hfq mutant was lower compared to the wt.

While focused on virulence genes, microarray analysis revealed a great number of Hfq-dependent genes which classified into several categories. Within the set of genes with known function and which were discovered to have higher transcript abundance in the hfq mutant emerged genes encoding the components of ABC transporter systems and amino acid transporters (“Surface proteins” category) and 30S and 50S ribosomal proteins (“Information transfer” category). Overexpression of both ABC transporter and ribosomal genes seems to be a typical phenomenon observed in hfq mutants of other bacteria. Among genes displaying deregulated transcript levels in the hfq mutant were also genes encoding several periplasmic and outer membrane proteins. This was not surprising as we have previously observed some differences in protein composition of the membrane fractions between the wt and Δhfq strain. Several studies indicated that lack of Hfq results in irregular production of OMPs which affects the membrane integrity and cell envelope stability. Collectively, misregulation of ABC transporter genes, unbalanced expression of ribosomal genes and aberrant expression/production of OMPs might result in pleiotropic defects of hfq mutants such as slower growth and therefore, may as well explain the observed growth deficiency of the Tohama I hfq mutant.

Materials and Methods

Bacterial strains and growth conditions

The Bordetella pertussis Tohama I strain and its isogenic hfq deletion mutant were grown on Bordet-Gengou agar (BGA) plates supplemented with 15% defibrinated sheep blood for 3 to 4 days at 37°C. For liquid cultures, bacteria were grown in modified Stainer-Scholte (SS) medium supplemented with 0.1% cyclodextrin (Sigma) and 0.5% casamino acids (Difco) at 37°C. Samples for RNA isolation and protein analyses were taken at exponential (OD600 ≈1.3 and 0.7 for wt and Δhfq strains, respectively) and early stationary (OD600 ≈3.0 and 1.3 resp.) phases of growth.

Immunoblotting

B. pertussis wt and Δhfq strains were cultivated in SS medium at 37°C and cell pellets from 1-ml aliquots were lysed and boiled in sample buffer. For analysis of supernatant fractions, 20-ml aliquots of cell-free supernatants were filtered through 0.22 μm filters, precipitates obtained with 10% (w/v) trichloracetic acid were washed with 80% acetone, dissolved in sample buffer and boiled. Samples equivalent to 1 OD600 unit (supernatants) or 0.1 OD600 unit (whole-cell lysates) were separated on SDS-PAGE gels and transferred on nitrocellulose membrane. Membranes were probed with polyclonal antibodies raised against Bsp22 (laboratory stock), Vag8 (kindly provided by Dimitri Diavatopoulos) and BrkA (kindly provided by Rachel Fernandez) followed by incubation with corresponding IgG antibodies conjugated with horse radish peroxidase. The antibody-antigen complexes were visualized using SuperSignal West Femto chemiluminescent substrate (Thermo) according to standard protocol.

Preparation of the microarrays

Oligonucleotide probes were designed for 3554 annotated coding sequences of Bordetella pertussis Tohama I strain (BX470248.1, GI:33591069) using the OligoArray v2.1 software. One oligonucleotide for each of the annotated genes was designed except for the non-coding genes (rRNA, tRNA and smRNA(ssrA)) and for the transposase genes of IS elements. Out of the total 3554 probes 3462 were 60-nucleotides long with a Tm between 80 and 99°C, 60 probes were 50-nucleotides long with a Tm between 70 and 99°C and 32 probes were 45-nucleotides long with a Tm between 70 and 99°C. Aminated oligonucleotides were obtained from Sigma Aldrich and were resuspended in 20 μL of spotting solution (sciSPOT-cDNA, Scenion) at a final concentration of 20 μM. 96 oligonucleotide probes (randomly selected among the 3554 probes to be produced in duplicate) were mixed together, spotted as one probe and used as positive control of hybridization (intrinsic control providing an average intensity value for each study). All designed probes including the control mixture of 96 probes and human gene oligonucleotide-probes (negative controls) were spotted in duplicate (non-adjacent) on aldehydsilane coated glass slides (Nexterion Slide AL, Schott) using a QArray II spotter (Genetix) equipped with 12 ArrayHit 946MP3 pins. Spotting was done at 20-22°C and in 50–55% relative humidity. After spotting, glass slides were incubated at 120°C for 1 h and then treated just before usage with 132 mM NaBH4 in PBS containing 25% of ethanol for 30 minutes. The slides were then washed twice in 0.2% SDS solution and twice in distilled water before being dried by centrifugation.

RNA preparation and microarray analysis

B. pertussis wt and Δhfq cells were grown in 2 independent biological replicates in SS medium and harvested by centrifugation (6000 rpm, 10 min, 4°C) at exponential and stationary phases of growth. Total RNA was extracted from pellets using Trizol (Sigma) as recommended by the manufacturer and further treated with DNase I (Promega) to remove any contaminating genomic DNA. Purified RNA was checked for concentration and integrity using the Bioanalyzer 2100 device (Agilent Technologies). Labeling with Cyanine-5 (Cy5) and Cyanine-3 (Cy3) was performed by random primed reverse transcription starting from 10 μg of purified total RNA and using SuperscriptIII (Invitrogen) reverse transcriptase according to manufacturer protocol and amplified cDNA was purified on Qiaquick PCR purification kit (Qiagen). A mix of 1 μg of Cy5-cDNA + 1 μg of Cy3-cDNA was then prepared in 490 μL of hybridization buffer (5× SSC, 40% deionised formamide, 5× Denhart's buffer, 0.1% SDS, 1 mM sodium pyrophosphate). The labeled mixture was used for hybridization using a gasket slide (Agilent Technologies) for 14 to 16 h at 52°C under agitation. For each sample a Cy5/Cy3 dye-swap technical replicate starting from the same RNA preparation was prepared and hybridized. After hybridization, slides were subsequently washed at room temperature in 2× SSC, 0.2%SDS for 5 minutes, in 0.5× SSC for 10 minutes, in 0.05× SSC for 5 minutes, and 0.01× SSC for 2 minutes before being dried and scanned using InnopScan700 (Innopsys). Background correction as well as within- and between-array normalization was performed from raw data using LIMMA (Linear Models for Microarray Data) package. After normalization, identification of genes with significantly modulated expression was performed using moderated Student's t-test with empirical Bayes shrinkage of standard errors. Statistics were corrected for multiple testing using false-discovery rate approach. In order to select genes significantly affected by hfq deletion, absolute log fold change values (log2FC) of transcript abundance was set to at least 1 (2-fold variation) and threshold for adjusted P-value (P) of less than 0.05 was considered as indicative of significance. Significantly affected genes were classified into functional categories as assigned by Park and coworkers with following modifications. The categories “Degradation of small molecules” and “Degradation of large molecules” were combined into category “Degradation,” category “Pathogenicity, Adaptation and Chaperones” was abbreviated to “Pathogenicity and Adaptation.” Genes comprising T3SS bsc and btr loci were moved from “Surface protein” category into “Pathogenicity and Adaptation” category. Statistical significance of overrepresentation of functional categories within the sets of genes with modified transcript abundance was assessed using the Fisher exact test and corrected for multiple testing using Benjamini and Hochberg method. Clustering of microarray data was performed using MeV software from TIGR.

Quantitative PCR (RT-qPCR)

To study the quantitative effects on gene expression, 1 μg of total isolated RNA was reverse transcribed in duplicates into cDNA following the manufacturer's instructions in a 25-μl reaction using the Reverse Transcription System (Promega). All primers (Table S7) were designed to anneal at 60°C and were analyzed for secondary structures and for cross-dimers in primer pairs. RT-qPCR was performed on Bio-Rad CFX96 instrument using SYBR® Green JumpStart™ Taq ReadyMix™ (Sigma). Briefly, 200 nmol/l of each primer together with 40 ng of reverse transcribed RNA were used in a 20-μl qPCR reaction volume, with an initial step at 95°C for 2 min, followed by 40 cycles of 95°C for 15 s, 60°C for 30 s, 72°C for 30 s, followed by recording of the melting curve. The rpoB gene was used as reference gene and relative gene expression was quantified using amplification efficiency values.

Infection of macrophages

B. pertussis cells were suspended in Dulbecco's Modified Eagle's medium (DMEM) supplemented with 10% bovine serum albumin and the RAW 264.7 murine macrophages were infected at a multiplicity of infection of 10 bacteria per cell (MOI 10). After addition of bacteria, tissue culture plates were centrifuged for 5 min at 300 g to facilitate bacterial interaction with macrophage cells. After 60 min of incubation at 37°C with 5% CO2, non-adherent bacteria were removed by 3 washing steps in DMEM. Then, medium containing 100 μg/ml of polymyxin B sulfate was added for 60 min at 37°C with 5% CO2 to kill extracellular bacteria. Next, infected macrophages were washed with fresh medium containing 30 μg/ml of polymyxin B sulfate and further incubated at 37°C with 5% CO2. In order to recover viable intracellular bacteria, infected macrophages were lysed with sterile water at different times post infection as specified in legends of . Lysates as well as bacteria which were not in contact with macrophages (time point 0; control sample) were plated onto BGA plates and obtained bacteria were subsequently cultivated in SS medium for RNA isolation and Westernblot analysis.

In vivo infection of mice

Four-week-old CD-1 Swiss outbred mice were obtained from Charles River Laboratories and maintained under specific pathogen free conditions. Overnight cultures of B. pertussis strains were diluted in sterile phosphate buffered saline (PBS, pH 7.4) to adjust appropriate cell concentrations. Sublethal doses (1 × 105 CFU/mouse) of both wt and Δhfq bacteria were administered intranasally in suspensions of 50 μl and in parallel, the number of delivered bacteria was determined by plating serial dilutions of challenge suspensions on BGA plates. Cohorts of 15 animals were challenged per strain and 3 mice per group were sacrificed 2 hours and 5, 8, 12 and 23 days post infection. The lungs were removed aseptically and homogenized in 2 ml of PBS using a cell grinder (Heidolph RZR2020). The resulting homogenates were serially diluted with PBS and plated in duplicates on BGA plates and colonies were counted to calculate the CFU values per mouse. Student's t test (SigmaPlot 1, Systat Software Inc., California, USA) was used to compare the mean CFU values between wt and Δhfq strains at each time point. Differences were considered significant at P values of ≤ 0.05. Recovered bacteria were subcultured in SS medium and samples for RNA isolation and westernblot analysis were taken at exponential phase of growth.

All animal experiments were approved by the Animal Welfare Committee of the Institute of Microbiology of the ASCR, v.v.i., in Prague, Czech Republic. Handling of animals was performed according to the Guidelines for the Care and Use of Laboratory Animals, the Act of the Czech National Assembly, Collection of Laws No. 149/2004, inclusive of the amendments on the Protection of Animals against Cruelty and Public Notice of the Ministry of Agriculture of the Czech Republic, Collection of Laws No. 207/2004, on care and use of experimental animals.