Refined analysis of the Campylobacter jejuni iron-dependent/independent Fur- and PerR-transcriptomes
© Butcher et al. 2015
Received: 27 October 2014
Accepted: 29 May 2015
Published: 4 July 2015
The genome of Campylobacter jejuni contains two iron activated Fur-family transcriptional regulators, CjFur and CjPerR, which are primarily responsible for regulating iron homeostasis and oxidative stress respectively. Both transcriptional regulators have been previously implicated in regulating diverse functions beyond their primary roles in C. jejuni. To further characterize their regulatory networks, RNA-seq was used to define the transcriptional profiles of C. jejuni NCTC11168 wild type, Δfur, ΔperR and ΔfurΔperR isogenic deletion mutants under both iron-replete and iron-limited conditions.
It was found that 202 genes were differentially expressed in at least one mutant under iron-replete conditions and 331 genes were differentially expressed in at least one mutant under iron-limited conditions. The CjFur and CjPerR transcriptomes characterized in this study were compared to those previously identified using microarray profiling and found to be more extensive than previously understood. Interestingly, our results indicate that CjFur/CjPerR appear to co-regulate the expression of flagellar biogenesis genes in an opposing and iron-independent fashion. Moreover the ΔfurΔperR isogenic deletion mutant revealed that CjFur and CjPerR can compensate for each other in certain cases, suggesting that both regulators may compete for binding to specific promoters.
The CjFur and CjPerR transcriptomes are larger than previously reported. In particular, deletion of perR results in the differential expression of a large group of genes in the absence of iron, suggesting that CjPerR may also regulate genes in an iron-independent manner, similar to what has already been demonstrated with CjFur. Moreover, subsets of genes were found which are only differentially expressed when both CjFur and CjPerR are deleted and includes genes that appear to be simultaneously activated by CjFur and repressed by CjPerR. In particular the iron-independent co-regulation of flagellar biogenesis by CjFur/CjPerR represents a potentially novel regulatory function for these proteins. These findings represent additional modes of co-regulation by these two transcriptional regulators in C. jejuni.
The microaerophilic bacterial pathogen Campylobacter jejuni can colonize and thrive in the gut of many different mammalian and avian hosts (humans, bovines, avians, porcines etc.) [1, 2], and C. jejuni is also known to persist in the environment during its transit between different hosts . In all these environments, C. jejuni must regulate several functions to be able to cope with these variables and often stressful conditions. Previous work has demonstrated that genes relating to iron metabolism and oxidative stress defense are key factors for C. jejuni’s survival in vivo [2, 4, 5]. The classical iron-regulatory protein in most Gram-negative bacteria, including C. jejuni, is the ferric uptake regulator (CjFur) protein, which senses intracellular iron and represses genes involved in iron acquisition accordingly [2, 6–8]. Interestingly C. jejuni also possesses a second iron responsive transcriptional regulator, the peroxide responsive regulator PerR that is responsible for regulating C. jejuni’s response to peroxide stress [4, 9]. This is in contrast with other gram negative bacteria whose peroxide sensors (e.g. OxyR, OhR) do not utilize iron as a cofactor to sense oxidative stress .
It has long been known that the members of the CjFur and CjPerR regulons overlap . The best example of this is the katA gene which encodes for catalase. Both CjFur and CjPerR must be deleted for this gene to be fully derepressed and unresponsive to iron repression . In contrast, there are other genes which appear to be solely regulated by either CjFur (e.g. cfrA, chuA) or CjPerR (e.g. ahpC) [2, 4, 8]. At times it has been difficult to clearly define whether a gene is regulated by CjFur and/or CjPerR due to the functional similarity of these two transcriptional regulators. Specifically as Fur family proteins, such as CjPerR, have been shown to use iron or manganese as co-factor in vitro, deletion of fur may indirectly influence CjPerR activity by altering iron levels in the cell. Moreover, the CjFur and CjPerR binding sequences are likely to be extremely similar based on studies of other organisms harboring multiple Fur family members .
There have been several studies that characterized the CjFur and CjPerR transcriptomes using microarray transcriptomic profiling of either Δfur or ΔperR isogenic mutant strains [2, 4, 8]. However directly comparing the profiles obtained in these studies can be difficult due to confounding factors that may alter gene expression independently from transcriptional regulator deletion. In C. jejuni these can include differences in growth conditions including the gas content used (e.g. gas packs vs. defined gas mixes) and the composition of the growth media [12, 13]. In addition, given that CjFur and CjPerR are known to co-regulate genes, there may be members of the CjFur regulon whose differences in expression are being masked by CjPerR and vice versa. We have therefore characterized the transcriptomes of our C. jejuni NCTC11168 wild-type strain and its corresponding Δfur, ΔperR and ΔfurΔperR isogenic deletion mutants using RNA-seq. These strains were grown following a standardized protocol to improve the reproducibility of the transcriptomic results. While there is increasing evidence that CjFur also regulates genes in vivo in the absence of iron (iron-independent regulation) , there has been little investigation as to whether CjPerR may also be able to regulate gene transcription in an iron-independent fashion. Thus the transcriptomes of C. jejuni NCTC11168 wild-type strain and its corresponding Δfur, ΔperR and ΔfurΔperR were characterized under iron-replete and iron-limited conditions to define the both iron-dependent /iron-independent transcriptomes of CjFur and CjPerR.
Our results indicate that the CjFur and CjPerR transcriptomes are larger than those previously reported in C. jejuni using microarray profiling. This is especially true with regards to iron-independent regulated genes for both CjFur and CjPerR. In addition we have found subsets of genes which are only differentially expressed when both CjFur and CjPerR are deleted and genes which appear to be simultaneously activated by CjFur and repressed by CjPerR.
Bacterial Strains and growth
The C. jejuni NCTC11168 wild-type, Δfur , ΔperR  and ΔfurΔperR  strains were routinely cultured on Mueller-Hinton (MH) agar plates under microaerophilic conditions (83 % N2, 4 % H2, 8 % O2 and 5 % CO2) at 37 °C in a MACS-VA500 workstation (Don Whitley, West Yorkshire, England). C. jejuni strains were grown overnight in MH agar overlaid with MH broth (biphasic cultures). Prior to inoculation, overnight cultures were washed in minimal essential media alpha (MEMα) (Invitrogen) to remove excess iron sources. Glass flasks used for growing strains for RNA extraction were soaked with 1 M HCl to solubilize bound metal ions and rinsed twice with ddH2O prior to sterilization. Strains were grown in a final volume of 50 mL MEMα supplemented with 10 mM pyruvate with freshly prepared FeSO4 added at a final concentration of 40 μM as needed to generate iron-replete growth conditions. In agreement with previous work , the various strains tested did not show any appreciable differences in their growth under the conditions tested (data not shown). All experiments were done in accordance with the University of Ottawa’s policies for working with biohazardous materials, permit #B05-BMI-1213.
Total RNA Isolation
C. jejuni was grown in MEMα under either iron-replete (+40 μM FeSO4) or iron-limited (no added iron) conditions until midlog phase was reached (OD600 ~ 0.2). Total RNA was then isolated as described previously [2, 7]. Briefly, a 10 % cold RNA stop solution used to preserve RNA integrity (10 % buffer saturated phenol pH 4.3 in absolute ethanol) was added to the media and cells were pelleted by centrifugation at 6000xg. Pellets were resuspended in TE buffer pH 8.0 and RNA was extracted using a hot-phenol chloroform method . RNA was precipitated overnight at −80 °C and the RNA pellet was washed five times in 80 % ethanol. RNA pellets were resuspended in TE buffer and treated with DNase I (Epicentre) to remove contaminating genomic DNA. Treated samples were purified using the RNeasy kit (Qiagen) and PCR was performed to confirm the absence of genomic DNA. RNA purity and quality was assessed using a Nanodrop and the Experion RNA StdSens analysis kit (Bio-Rad). Samples with high quality RNA were selected for Illumina sequencing library construction. Total RNA was isolated from three independent biological replicates of wild-type, Δfur, ΔperR and ΔfurΔperR strains under either iron-replete or iron-limited conditions (24 RNA extractions in total).
Illumina library construction and sequencing
Library construction, cluster generation and sequencing were performed at McGill University and Génome Québec Innovation (Montreal, Canada). Ribosomal RNAs were depleted using the RiboZero Meta-Bacteria kit (Epicentre) and strand-specific sequencing libraries were constructed using TruSeq kits (Epicentre). Cluster generation was done using standard Illumina protocols with sequencing performed on an Illumina HiSeq 2500 sequencer. The libraries were subjected to paired-end sequencing with a target read length of 100 nucleotides and the samples for iron-replete and iron-limited conditions were run on separate lanes (12 samples/lane). At least 12 million reads were obtained per library (Additional file 1: Table S1). The Illumina reads have been submitted to the Sequence Read Archive at NCBI with the accession number SRP044881.
Read alignment to the C. jejuni genome and bioinformatic analysis
Summary of differentially expressed genes in each strain under different growth conditions
Increased under iron limitation
Decreased under iron limitation
Increased under iron replete conditions
Decreased under iron replete conditions
Bioinformatic analysis of potential Fur and PerR consensus binding sites
The promoter regions of Fur/PerR regulated genes under either iron-limited or replete conditions were analyzed for the existence of potential binding sites. Genes were analyzed as a whole and also as subgroups based on their potential mode of regulation (iron dependent/independent and activation/repression). Regions consisting of 300 bp upstream the predicted start codon of each gene were retrieved from the RSAT server (http://www.rsat.eu/)  and analyzed for potential regulatory motifs by the MEME suite of analysis tools (http://meme-suite.org/) . The genes containing the iron-dependent Fur repression motif consisted of ceuB, cfbpA, cfrA, chaN, chuA, chuZ, cj0176c, cj0818, cj1384c, cj1587c, cj1658, dsbI, exbB2, katA and tonB3.
RT-qPCR confirmation of selected genes
RT-qPCR was used to confirm the differential expression under iron limited conditions of several genes including ahpC, chuA, cj0414, cj0415, cj1386, flgE1, flgE2, flgG, flgH, and fliK. These genes were selected based on the fact that they show inverse regulation by CjFur/CjPerR (flgE1, flgE2, flgG, flgH, and fliK) or where only differentially expressed in the ΔperR mutant (cj0414, cj0415, cj1386). AhpC and chuA were also included as controls for potential metal dependent regulation by CjPerR or CjFur respectively. RT-qPCR reactions were performed as previously described using slyD as an endogenous control . Genes were considered significantly expressed at a fold change >1.5 with a p value <0.05.
Confirming sample consistency via PCA analysis
Principal component analysis (PCA) was performed to show the overall structure of the transcriptomic data and identify sample outliers based on the growth media (iron-replete vs. deplete) or genetic background (wild-type vs. Δfur vs. ΔperR vs. ΔfurΔperR). The PCA showed distinct sample separation based on the iron content of the growth medium, highlighting the key influence of iron on the C. jejuni transcriptome (Fig. 1a). The ΔfurΔperR samples clustered closer to the Δfur samples as compared to the ΔperR samples in both iron replete and iron limited conditions. Moreover the samples for each strain tended to cluster together with two noticeable exceptions (Fig. 1a, b, arrows). One of the ΔperR samples (perr_plusIron_3) grown in the presence of iron did not group with other iron-replete samples and instead clustered with iron-limited samples (Fig. 1a). In addition one of the Δfur samples (fur_minusIron_2) grown in iron-limited conditions shows a distinctly different clustering pattern from its mates and the ΔfurΔperR samples in the PC1vsPC3 plot (Fig. 1b). Therefore both of these outlier samples were removed from downstream differential expression analysis, to avoid introducing artifacts or confounding the identification of differentially expressed genes. To note, including these outlier samples resulted in a significantly reduced number of differentially expressed genes with a high number of false negatives due to increased variation in the gene expression profiles (data not shown).
Previous work in our laboratory has also used RNA-seq to identify iron responsive genes in wild-type C. jejuni NCTC11168 . These previous sequencing results were compared to the current wild-type results to determine the similarity of the runs and whether the results from both studies could be merged for meta-analysis. PCA clearly separated the transcriptomes based on the sequencing runs (Additional file 11: Figure S4). In fact, the sequencing results for these two studies could be separated by two axes corresponding to either growth condition (iron-replete vs iron-limited) or study (current vs previous). Therefore to avoid biased results, the iron regulon was reanalyzed using only the data from this study and compared to the iron regulon that was determined previously.
Transcriptome signatures of the Δfur, ΔperR and ΔfurΔperR mutants under iron-replete conditions
The major objective of this work was to define the CjFur and CjPerR regulatory networks. Because CjFur and CjPerR regulated genes are known to be iron-repressed we first determined the CjFur and CjPerR regulon under iron-replete conditions as these genes would be differentially expressed between the wild-type and mutant strains. Globally, 86 genes were found to be differentially expressed in the Δfur mutant (73 up-regulated and 13 down-regulated), 111 genes in the ΔperR mutant (85 up-regulated and 26 down-regulated) and 95 genes in the ΔfurΔperR mutant (72 up-regulated and 23 down-regulated) (Additional file 6: Table S5). Genes found to be differentially expressed in at least one strain were subjected to hierarchical clustering analysis (using average linkage Euclidean distance). The analysis generated 8 main clusters, with clusters A-C representing down-regulated genes and clusters D-H representing up-regulated genes (Additional file 6: Table S5, Fig. 2). Clusters A, B and C contain genes that were primarily down-regulated in the Δfur, ΔfurΔperR or ΔperR strains respectively. The most down-regulated Δfur genes (Cluster A) consisted of the rrc-cj0011c operon (29.2 and 5.8 fold activation). Cluster B primarily consists of genes highly down-regulated (~1.5-8 fold) in the ΔfurΔperR strain, including cj1295-1298 and several flagellar genes such as flgM, fliS and flaA. Cluster C contains genes primarily down-regulated in the ΔperR strain (~2 fold) and includes several motility associated genes (cj1340c-cj1345c, maf6, maf7). Interestingly, no genes were found to be down-regulated in all three strains and only cj1340c (Δfur, ΔperR), flaA (Δfur, ΔfurΔperR) and cj1442c (ΔperR, ΔfurΔperR) were down-regulated in two of the three strains (Additional file 12: Figure S6).
Clusters D-F consist of genes that were slightly up regulated in at least one of the three strains (~2 fold). In contrast with the down regulated genes, many up regulated genes were up-regulated in multiple strains. Clusters G and h contain genes that were highly up-regulated in at least two strains (~6-120 fold). Genes in cluster g are greatly up-regulated in the Δfur and ΔfurΔperR strains. Most of these genes are directly involved in iron acquisition (e.g. hemin: chuABCDZ; enterobactin: cfrA, ceuBCDE; lactoferrin/transferrin: cfbpABC, ctuA, chaN; rhodotorulic acid: cj1658-p19, cj1660-1665; ferrous ions: feoB,cj1397) or supplying the energy required for iron uptake across the outer membrane (e.g. tonB1/B2/B3, exbB1/B2, exbD1/D2). Many of these genes have been previously shown to be directly regulated by Fur [8, 27, 28]. Cluster H consists of four genes that are highly up-regulated in all three strains and are involved in oxidative stress defense (katA, cj1386, ahpC and trxB). The C. jejuni alkyl hydroperoxidase (ahpC) has been previously shown to be directly regulated by PerR , while the catalase (katA), thioredoxin trxB and the heme trafficking protein cj1386 have been suggested to be PerR regulated [4, 9]. In addition katA, trxB and cj1386 have been proposed to be Fur regulated [2, 8], although Fur has only been shown to bind the katA promoter  (note that katA and cj1386 are not co-transcribed ). Most of the genes in both clusters G and H have also been previously found to be iron responsive in C. jejuni [2, 7, 8].
The promoter regions of differentially expressed genes in each strain (Δfur, ΔperR and ΔfurΔperR) under each condition (iron-replete and iron-limited) were analyzed by MEME to identify potential binding motifs. However we were only able to identify a consensus binding motif for Fur iron-dependent repressed targets (Additional file 13: Figure S5) that was essentially identical to those already published [2, 6].
Transcriptome signatures of the Δfur, ΔperR and ΔfurΔperR mutants under iron-limited conditions
In order to identify genes regulated by CjFur and CjPerR in an iron-independent fashion we also characterized the CjFur and CjPerR transcriptomes under iron-limited conditions. Globally, 31 genes were found to be differentially expressed in the Δfur mutant (11 up-regulated and 20 down-regulated genes), 373 genes in the ΔperR mutant (160 up-regulated and 113 down-regulated genes) and 72 genes in the ΔfurΔperR mutant (49 up-regulated and 23 down-regulated genes) (Additional file 7: Table S6). Hierarchical clustering analysis of the differentially expressed genes revealed 11 main clusters (A-K) of expression profiles with clusters I-J being primarily composed of genes that were down-regulated in at least one strain and clusters C, D, G and H containing genes that were up-regulated in at least one strain (Fig. 3; ~2-4 fold). Cluster A is composed of genes that were strongly up-regulated in the Δfur and ΔfurΔperR strains such as cj1356c, fdxA, rrc and cj0948c-0949c (~4-8 fold). Cluster I contains genes that were strongly down-regulated in the Δfur and ΔfurΔperR strains such as flaA, and the flagellar modification island genes cj1296-97 (~3-4 fold).
Cluster B contains genes up-regulated in the ΔperR mutant while being down-regulated in the Δfur and ΔfurΔperR mutants, while cluster F consists of genes down-regulated in the ΔperR mutant while being up-regulated in the Δfur and ΔfurΔperR mutants. Interestingly Cluster E is primarily composed of genes that were strongly down-regulated in the Δfur mutant (~3-4 fold down-regulation) and strongly up-regulated in the ΔperR mutant (~4-6 fold up-regulation) (also see Additional file 12: Figure S6). These include many different flagellar biogenesis genes such as flgE1, flgE2, flgG, flgH, flgI, flgJ, and fliK (to note, most of these genes have been previously shown to be required for C. jejuni motility [5, 31]). In addition, none of the genes present in Cluster E were found to be differentially expressed in the ΔfurΔperR mutant (Additional file 7: Table S6).
Gene-set enrichment analysis of Δfur, ΔperR and ΔfurΔperR under both iron-rich and iron-limited conditions
Cross platform correlation of the Δfur and ΔperR transcriptomes
The commonly identified up-regulated ΔperR genes under iron-replete conditions include the antioxidant genes katA, trxB and rrc as well as a variety of genes involved in flagellar biogenesis (e.g. flaD, flgD, fgI). In contrast with the microarray profiling, the RNA-seq analysis also identified the antioxidant genes ahpC and cj1386 as being CjPerR regulated. There were no CjPerR activated genes in common between the RNA-seq and microarray profiling under iron-replete conditions. The RNA-seq analysis for the ΔperR mutant under iron-limited conditions identified many differentially expressed genes that were not present in the microarray profiling. Many of the commonly identified up-regulated genes are also differentially expressed under iron-replete conditions (e.g. katA, cj1386, flgI, flgD). The ΔperR mutant under iron-limited conditions was also the only condition where there was a disagreement between the RNA-seq and microarray data for whether a gene was significantly up-regulated or down-regulated upon perR deletion (Fig. 6).
Comparison of the iron regulon from two independent RNA-seq experiments and microarray profiling
Overall, the RNA-seq analysis identified more differentially expressed genes as compared to the microarray analysis. Approximately 26 % of the genes identified in the RNA-seq analysis were also found to be differentially expressed in the microarray profiling experiments (180/686). In addition, there were 75 genes identified in the microarray analysis that were not found to be differentially expressed via RNA-seq (Additional file 10: Table S7).
In contrast, all of the genes identified as being iron regulated in a previous RNA-seq study were also found to be differentially expressed in this work . However, the fold changes determined in this previous RNA-seq study appear to be greater as compared to this work. In particular, the Log2(FC)’s in the previous study exceeded >25 in some cases as compared to a maximum of Log2(FC) of ~10 in this work (Fig. 7). Indeed, most of these extremely differentially expressed genes (e.g. cj0243c, cj0878, cj0466, cj0730, cj1722) were not identified as significant in the previous RNA-seq study due to high variability in the number of reads across biological replicates . However, all these genes were found to be significantly differentially expressed in this dataset.
Although the importance of the CjFur and CjPerR transcriptional regulators on C. jejuni’s in vivo fitness is now well established, there are still gaps in our understanding of their regulatory networks. Importantly, and in contrast to many other bacteria, CjFur regulates gene expression in an iron independent manner. However it is unknown whether CjPerR is also able to regulate genes independently from iron. As well, the degree to which the CjFur/CjPerR transcriptomes overlap is unclear as the transcriptome of a double deletion of fur and perR has not been described. In this study, we have used RNA-seq to profile the transcriptomes of Δfur, ΔperR and ΔfurΔperR isogenic deletion mutants in order to characterize the CjFur and CjPerR transcriptomes under both iron-replete and iron-limited conditions. In contrast to previously reported results, we have found that the transcriptomes for both of these metalloregulators are more extensive than anticipated. We have also found that these transcriptomes appear to reveal roles for both proteins in iron-dependant and iron-independent regulatory functions. Our work also highlights the importance of having multiple biological replicates when conducting genome wide analyses and to identify and remove outlier samples that may mask expression differences.
Unsurprisingly, the most up-regulated genes under iron-replete conditions in the isogenic deletion mutants are those previously characterized as being involved in iron acquisition (Fig. 2 [Cluster G]) or in oxidative stress defense (Fig. 2 [Cluster H]). These two functions were clearly separated with the iron acquisition genes up-regulated in the Δfur mutant and the oxidative stress defense genes up-regulated in the ΔperR mutant. The Δfur mutant had relatively few genes down-regulated under iron-replete conditions with only rrc and cj0011c forming a distinct CjFur activated cluster. In contrast, there were many more genes down-regulated in the ΔperR mutant and in the ΔfurΔperR mutant (Fig. 2 [Clusters BC]), suggesting that PerR may also activate gene expression in an iron-dependent form.
Interestingly there were a large number of genes differentially expressed in the isogenic deletion mutants under iron-limited conditions, which was unexpected as both CjFur and CjPerR are metalloregulators. In particular, there were 31 and 272 genes differentially expressed in the Δfur and ΔperR mutants respectively (Fig. 3). While CjFur has been previously shown to regulate gene expression in an iron-independent manner, this result suggests that CjPerR may also regulate genes independently from iron. However, additional studies using purified CjPerR protein without a metal cofactor would be required to demonstrate this specific mode of regulation as transcriptomic studies cannot distinguish between direct and indirect regulation.
Under iron-limited conditions the Δfur mutant also showed relatively few down-regulated genes as compared to up-regulated genes. The oxidative stress defense genes cj1386, fdxA and rrc were found to be regulated by CjFur independently from iron. Rrc has been previously shown to be differentially expressed in a Δfur mutant under iron-limited conditions by RT-qPCR and directly bound by metal free CjFur . Previous work has also shown cj1345c to be differentially expressed in Δfur mutant under iron-limited conditions using RT-qPCR and also shown to be directly bound by metal free CjFur via EMSA . However in this study, cj1345c’s change in expression level was not statistically significant. There was also a cluster of genes that appear to be up-regulated in the Δfur mutant and down-regulated in the ΔperR mutant under iron-limited conditions (Fig. 3 [Cluster E], Additional file 12: Figure S6). Interestingly these genes are primarily composed of flagellar biogenesis genes (flgD, flgE1, flgE2, flgG2, flgH, flgI, flgJ, fliK). Indeed, GSEA found that the KEGG category for flagellar assembly was enriched in both the Δfur and ΔperR iron-limited conditions, with most of the flagellar assembly apparatus under opposing regulation in these two mutants (Additional file 14: Table S8, Additional file 15: Table S9, Fig. 5). Further supporting the fact that these genes appear to be differentially regulated by CjFur/CjPerR is the fact that these genes were not differentially expressed in the ΔfurΔperR mutant (Additional file 12: Figure S6), suggesting that each protein exerts its influence independently and GSEA did not identify flagellar assembly as being significantly enriched in the ΔfurΔperR mutant (Additional file 14: Table S8, Additional file 15: Table S9). This represents another manner by which the CjFur and CjPerR transcriptomes may overlap with each other. In the other cases where CjFur and CjPerR are known to co-regulate the same gene (e.g. katA repression), both transcriptional regulators have the same effect on transcription. These interactions between the CjFur and CjPerR transcriptomes further contrast with other bacteria, where their corresponding Fur and PerR transcriptomes are separate from each other . It should also be noted that C. jejuni perR has also been shown to be autoregulated and thus this may account for some of the differences seen in the ΔperR mutant [4, 29].
When comparing the transcriptomes identified in this study using RNA-seq (CjFur, CjPerR and iron) to those identified using microarray profiling the most obvious observation is that the RNA-seq identified substantially more differentially expressed genes as compared to the microarray profiling (Fig. 6). This is likely due to several factors. Firstly, the greater dynamic range offered by RNA-seq allows for the assignment of genes being differentially expressed at a lower fold change as compared to microarrays (1.5 fold as compared to 2 fold) which naturally leads to more genes being assigned as differentially expressed. Indeed, most of the genes identified using both techniques were those that showed the greatest level of differential expression. In addition, microarrays often have difficulties identifying differentially expressed genes which are expressed at high levels in both the test and control conditions due to signal saturation. For example, previous microarray profiling did not detect any significant change in ceuBC/ctuA in the Δfur mutant and ahpC in the ΔperR mutant [2, 4]. This is despite the presence of Fur/PerR boxes in each promoter region and gel shift assays showing that each promoter was directly bound by either CjFur/CjPerR [8, 27, 29]. As RNA-seq works by directly counting the reads aligning to each gene it does not suffer from this same limitation and correctly identified these genes as being differentially expressed in the CjFur/CjPerR mutants respectively.
In this study, PCA analysis was used to determine the relative similarity of each strain’s biological replicates. This analysis identified two samples that appeared to be substantially different from their mates and that were subsequently discarded in the downstream analyses. Including these samples in the differential expression analysis leads to markedly different results depending on each strain. Including the erroneous ΔperR iron-replete sample in the differential analysis leads to an iron-limited gene profile (i.e. all iron acquisition genes become significantly up-regulated) whereas including the erroneous Δfur iron-limited sample leads to essentially no genes being identified as differentially expressed (data not shown). The potential for variability between samples influencing results is highlighted when comparing the iron regulon determined in this study with a previously published dataset . In the previous study, the iron-replete samples showed a greater level of variability (Additional file 11: Figure S4) that resulted in many genes not being identified as significantly differentially expressed. However, it should be noted that all the genes actually identified as being iron regulated were confirmed in this work. Moreover, while some of the variability present between these aforementioned samples could be deduced by plotting sample RPKMs and determining the R2 value for each pair of replicates, some samples identified as outliers in the PCA analysis could not be confidently distinguished from their mates based on R2 values alone (data not shown). As such, higher order analyses for determining sample similarity (PCA vs R2 values etc.) and multiple (>2) replicates for different samples are highly recommended for future studies employing RNA-seq for differential gene expression analysis. Another factor that could influence the effectiveness of RNA-seq analyses is the number of reads aligning to coding regions (i.e. sequencing depth). While sequencing depth is a difficult concept to define in transcriptomic studies (due to variations in the complexity of each unique transcriptome), increasing the number of reads aligning to each gene will likely increase the confidence in the downstream analyses . The number of reads aligning to genes is 10x greater in this current study as compared to the previous RNA-seq study identifying iron regulated genes in C. jejuni . This greater sequencing depth appears to have more successfully captured the distribution of transcripts expressed at low levels in either iron-replete/iron-limited conditions and resulted in attenuated fold changes as compared to previous work (Fig. 7, Panel b).
In conclusion, our results indicate that both the CjFur and CjPerR transcriptomes are larger than previously believed. In particular, deletion of perR appears to result in the differential expression of a large group of genes in the absence of iron, suggesting that CjPerR may also regulate genes in an iron-independent manner; similar to what has already been demonstrated with CjFur. We have also found subsets of genes which are only differentially expressed when both CjFur and CjPerR are deleted. Finally, we have genes which appear to be simultaneously activated by CjFur and repressed by CjPerR, which represents another manner in which these two transcriptomes interact with each other.
Availability of supporting data
The raw sequencing reads for this study have been submitted to the Sequence Read Archive at NCBI (http://www.ncbi.nlm.nih.gov/sra) with the accession number SRP044881. Other supporting information is available in the supplemental files associated with this manuscript.
This work was supported CIHR (MOP#84224) to AS and an Ontario Graduate Scholarship to JB. This work was also supported in part by a Norwich Research Park PhD studentship to RH, and the Biotechnology and Biological Sciences Research Council Institute Strategic Programme Grant BB/J004529/1. The authors also thank Annika Flint for her constructive comments on this manuscript.
- Young KT, Davis LM, Dirita VJ. Campylobacter jejuni: molecular biology and pathogenesis. Nat Rev Microbiol. 2007;5(9):665–79.PubMedView ArticleGoogle Scholar
- Palyada K, Threadgill D, Stintzi A. Iron acquisition and regulation in Campylobacter jejuni. J Bacteriol. 2004;186(14):4714–29.PubMed CentralPubMedView ArticleGoogle Scholar
- Bronowski C, James CE, Winstanley C. Role of environmental survival in transmission of Campylobacter jejuni. FEMS Microbiol Lett. 2014;356(1):8–19.PubMedView ArticleGoogle Scholar
- Palyada K, Sun YQ, Flint A, Butcher J, Naikare H, Stintzi A. Characterization of the oxidative stress stimulon and PerR regulon of Campylobacter jejuni. BMC Genomics. 2009;10:481.PubMed CentralPubMedView ArticleGoogle Scholar
- Flint A, Sun YQ, Butcher J, Stahl M, Huang H, Stintzi A. Phenotypic screening of a targeted mutant library reveals Campylobacter jejuni defenses against oxidative stress. Infect Immun. 2014. doi:10.1128/IAI.01528-13.PubMed CentralPubMedGoogle Scholar
- Butcher J, Sarvan S, Brunzelle JS, Couture JF, Stintzi A. Structure and regulon of Campylobacter jejuni ferric uptake regulator Fur define apo-Fur regulation. PNAS. 2012;109(25):10047–52.PubMed CentralPubMedView ArticleGoogle Scholar
- Butcher J, Stintzi A. The transcriptional landscape of Campylobacter jejuni under iron replete and iron limited growth conditions. PloS one. 2013;8(11), e79475.PubMed CentralPubMedView ArticleGoogle Scholar
- Holmes K, Mulholland F, Pearson BM, Pin C, McNicholl-Kennedy J, Ketley JM, et al. Campylobacter jejuni gene expression in response to iron limitation and the role of Fur. Microbiology. 2005;151(Pt 1):243–57.PubMedView ArticleGoogle Scholar
- van Vliet AH, Baillon ML, Penn CW, Ketley JM. Campylobacter jejuni contains two Fur homologs: characterization of iron-responsive regulation of peroxide stress defense genes by the PerR repressor. J Bacteriol. 1999;181(20):6371–6.PubMed CentralPubMedGoogle Scholar
- Dubbs JM, Mongkolsuk S. Peroxide-sensing transcriptional regulators in bacteria. J Bacteriol. 2012;194(20):5495–503.PubMed CentralPubMedView ArticleGoogle Scholar
- Fuangthong M, Helmann JD. Recognition of DNA by three ferric uptake regulator (Fur) homologs in Bacillus subtilis. J Bacteriol. 2003;185(21):6348–57.PubMed CentralPubMedView ArticleGoogle Scholar
- John A, Connerton PL, Cummings N, Connerton IF. Profound differences in the transcriptome of Campylobacter jejuni grown in two different, widely used, microaerobic atmospheres. Res Microbiol. 2011;162(4):410–8.PubMedView ArticleGoogle Scholar
- Verhoeff-Bakkenes L, Arends AP, Snoep JL, Zwietering MH, de Jonge R. Pyruvate relieves the necessity of high induction levels of catalase and enables Campylobacter jejuni to grow under fully aerobic conditions. Lett Appl Microbiol. 2008;46(3):377–82.PubMedView ArticleGoogle Scholar
- van Vliet AH, Wooldridge KG, Ketley JM. Iron-responsive gene regulation in a Campylobacter jejuni fur mutant. J Bacteriol. 1998;180(20):5291–8.PubMed CentralPubMedGoogle Scholar
- Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10(3):R25.PubMed CentralPubMedView ArticleGoogle Scholar
- Field D, Tiwari B, Booth T, Houten S, Swan D, Bertrand N, et al. Open software for biologists: from famine to feast. Nature Biotechnol. 2006;24(7):801–3.View ArticleGoogle Scholar
- Carver T, Harris SR, Berriman M, Parkhill J, McQuillan JA. Artemis: an integrated platform for visualization and analysis of high-throughput sequence-based experimental data. Bioinformatics. 2012;28(4):464–9.PubMed CentralPubMedView ArticleGoogle Scholar
- R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing; 2014.Google Scholar
- Robinson MD, McCarthy DJ, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.PubMed CentralPubMedView ArticleGoogle Scholar
- Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26(6):841–2.PubMed CentralPubMedView ArticleGoogle Scholar
- Sturn A, Quackenbush J, Trajanoski Z. Genesis: cluster analysis of microarray data. Bioinformatics. 2002;18(1):207–8.PubMedView ArticleGoogle Scholar
- Luo W, Friedman MS, Shedden K, Hankenson KD, Woolf PJ. GAGE: generally applicable gene set enrichment for pathway analysis. BMC Bioinformatics. 2009;10:161.PubMed CentralPubMedView ArticleGoogle Scholar
- Luo W, Brouwer C. Pathview: an R/Bioconductor package for pathway-based data integration and visualization. Bioinformatics. 2013;29(14):1830–1.PubMed CentralPubMedView ArticleGoogle Scholar
- Krzywinski M, Birol I, Jones SJ, Marra MA. Hive plots–rational approach to visualizing networks. Briefings Bioinformatics. 2012;13(5):627–44.View ArticleGoogle Scholar
- Thomas-Chollier M, Defrance M, Medina-Rivera A, Sand O, Herrmann C, Thieffry D, et al. RSAT 2011: regulatory sequence analysis tools. Nucleic acids research. 2011;39(Web Server issue):W86–91.PubMed CentralPubMedView ArticleGoogle Scholar
- Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, et al. MEME SUITE: tools for motif discovery and searching. Nucleic Acids Res. 2009;37(Web Server issue):W202–208.PubMed CentralPubMedView ArticleGoogle Scholar
- Miller CE, Rock JD, Ridley KA, Williams PH, Ketley JM. Utilization of lactoferrin-bound and transferrin-bound iron by Campylobacter jejuni. J Bacteriol. 2008;190(6):1900–11.PubMed CentralPubMedView ArticleGoogle Scholar
- Ridley KA, Rock JD, Li Y, Ketley JM. Heme utilization in Campylobacter jejuni. J Bacteriol. 2006;188(22):7862–75.PubMed CentralPubMedView ArticleGoogle Scholar
- Kim M, Hwang S, Ryu S, Jeon B. Regulation of perR expression by iron and PerR in Campylobacter jejuni. J Bacteriol. 2011;193(22):6171–8.PubMed CentralPubMedView ArticleGoogle Scholar
- Flint A, Sun YQ, Stintzi A. Cj1386 is an ankyrin-containing protein involved in heme trafficking to catalase in Campylobacter jejuni. J Bacteriol. 2012;194(2):334–45.PubMed CentralPubMedView ArticleGoogle Scholar
- Kamal N, Dorrell N, Jagannathan A, Turner SM, Constantinidou C, Studholme DJ, et al. Deletion of a previously uncharacterized flagellar-hook-length control gene fliK modulates the σ54-dependent regulon in Campylobacter jejuni. Microbiology. 2007;153(Pt 9):3099–111.PubMedView ArticleGoogle Scholar
- Haas BJ, Chin M, Nusbaum C, Birren BW, Livny J. How deep is deep enough for RNA-Seq profiling of bacterial transcriptomes? BMC Genomics. 2012;13:734.PubMed CentralPubMedView ArticleGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.