- Research article
- Open Access
Defining the nitrogen regulated transcriptome of Mycobacterium smegmatis using continuous culture
BMC Genomics volume 16, Article number: 821 (2015)
Nitrogen is essential for microbial growth and its importance is demonstrated by the complex regulatory systems used to control the transport, assimilation and metabolism of nitrogen. Recent studies are beginning to shed light on how mycobacteria respond to nitrogen limitation and several regulators (e.g., GlnR, PII) have been characterized at a molecular level. However, despite this progress, our knowledge of the transcriptional response of mycobacteria to nitrogen limitation and its regulation is confined to batch culture.
To gain further insight into the response of mycobacteria to nitrogen limitation, we developed a nitrogen-limited chemostat. We compared the transcriptional response of nitrogen-limited cells to carbon-limited cells using RNA-seq analysis in a continuous culture model at a constant growth rate.
Our findings revealed significant changes in the expression of 357 genes (208 upregulated, 149 downregulated; >2-fold change, false discovery rate <5 %) in response to nitrogen limitation in continuous culture. The vast majority of the GlnR regulon (68 %) was differentially expressed under nitrogen limitation in continuous culture and approximately 52 % of the 357 genes overlapped with a previously published study investigating the response of M. smegmatis to nitrogen limitation in batch culture, while expression of only 17 % of the genes identified in batch culture were affected in our chemostat model. Moreover, we identified a unique set of 45 genes involved in the uptake and metabolism of nitrogen that were exclusive to our chemostat model. We observed strong downregulation of pathways for amino acid catabolism (i.e., alanine, aspartate, valine, proline and lysine), suggesting preservation of these amino acids for critical cellular function. We found 16 novel transcriptional regulators that were directly or indirectly involved in the global transcriptomic response of M. smegmatis to nitrogen limitation and identified several non-coding RNAs that might be involved in the transcriptional or post-transcriptional regulation of nitrogen-regulated gene expression.
Using nitrogen-limited continuous culture we identified the nitrogen-responsive transcriptome of M. smegmatis, including a number of small non-coding RNAs implicated in controlling nitrogen-regulated gene expression.
Carbon and nitrogen are key components of organic material and their availability in the environment is necessary for growth and survival of microorganisms. The major mechanisms of carbon source utilization and its regulation are well defined for many groups of microorganisms, but our knowledge of nitrogen metabolism and its regulation is largely confined to enteric bacteria (for review see [1–3]). The importance of nitrogen for microbial growth is demonstrated by the complex regulatory systems used to control nitrogen assimilation in response to internal and external nitrogen levels. In some bacteria, nitrogen availability is sensed by the intracellular ratio of 2-oxoglutarate:glutamine . This signal ensures that the uptake of nitrogen sources are commensurate with the metabolic requirements of the organism i.e. carbon or nitrogen. High 2-oxoglutarate levels signal nitrogen limitation, while high glutamine levels signal nitrogen excess in bacteria [1, 4, 5].
In enteric bacteria, nitrogen metabolism is regulated by the two-component regulatory system NtrBC and the signal transduction protein PII . These regulatory systems coordinate nitrogen metabolism by regulating genes involved in ammonium assimilation, amino acid transport and nitrate metabolism through protein-protein interactions [1, 4, 6, 7]. An integral component of this regulation is GlnD, an uridylyl transferase/uridylyl-removing enzyme that senses and transmits the nitrogen status of the cell to the PII protein . The PII protein in its uridylylated state stimulates phosphorylation of the response regulator NtrC by the NtrB kinase, which results in expression of the NtrBC-regulated genes under nitrogen limitation . Homologues of the PII protein are found in many bacterial genomes, where they are mainly involved in the regulation of nitrogen metabolism. However, the PII protein may have additional roles, as recently described in cyanobacteria, where it regulates the metabolism of inorganic carbon .
Actinomycetes are a group of Gram-positive bacteria that inhabit a wide range of aquatic, terrestrial and human habitats. In actinobacteria, two different regulatory systems have been identified that regulate the transcriptional response to nitrogen limitation, AmtR, a TetR-type transcriptional regulator and GlnR, an OmpR-type transcriptional regulator [11–13]. AmtR was identified in Corynebacterium glutamicum to regulate more than 30 genes involved in nitrogen metabolism, including the PII protein, glutamine synthetase, glutamate synthase and urease . In Streptomyces coelicolor, GlnR mediates the transcriptional response to nitrogen limitation by elevating expression of the amtB-glnK-glnD operon and the glutamine synthetase glnA, while repressing transcription of the glutamate dehydrogenase gdhA [13, 14].
Slow-growing mycobacteria such as Mycobacterium tuberculosis appear to harbor only GlnR that was shown to regulate the expression of at least 33 genes in response to nitrogen limitation . M. tuberculosis can exploit different nitrogen sources, however, organic nitrogen sources (e.g. asparagine and aspartate) are preferred for extracellular and intracellular growth [16, 17]. Amino acids can function as both carbon and nitrogen sources, and although M. tuberculosis can metabolize a variety of amino acids in vitro, only a small number (aspartate, glutamate, asparagine and glutamine) can support growth of M. tuberculosis at the acidic pH that prevails in host macrophages . Inorganic nitrogen sources such as ammonium and nitrate are less efficient for growth than amino acids , but the reasons for this remain unknown.
The genome of the fast growing saprophytic actinobacterium Mycobacterium smegmatis contains copies of both global nitrogen regulators, GlnR and AmtR . It has recently been demonstrated that GlnR regulates the expression of more than one hundred genes in response to nitrogen limitation in M. smegmatis , while the AmtR regulon remains unknown. Previous microarray studies using batch culture nitrogen run out experiments have identified 1090 genes in M. smegmatis that are differentially expressed in response to nitrogen limitation . However, only a small subset of these genes is under control of GlnR, while the regulatory mechanisms for the majority have yet to be identified [20, 21]. During nitrogen run out studies, the growth rate changes significantly as nitrogen becomes depleted and therefore these types of experiments may not unequivocally uncover genes that respond solely to nitrogen limitation.
The aim of the current study was to deliver a molecular framework for how a mycobacterial cell responds to nitrogen-depleted and -replete conditions at a defined growth rate. To address this aim, we developed a nitrogen-limited chemostat for M. smegmatis at a constant growth rate of 0.12 h−1 (td = 5.7 h) followed by RNA-seq analysis to identify genes responding to nitrogen limitation.
Results and discussion
Development and validation of nitrogen-limited continuous culture to understand the global transcriptomic response to nitrogen limitation
The first step of this study was to establish a defined minimal medium limited for either nitrogen or carbon to identify the molecular response of M. smegmatis to nitrogen limitation. We modified HdB medium  by replacing 0.05 % (w/v) Tween-80 with 0.05 % (w/v) Tyloxapol, which cannot be metabolized as a carbon source, and (NH4)2SO4 with NH4Cl and K2SO4 (11.4 mM) to avoid simultaneous limitation for sulphur and nitrogen. Inorganic nitrogen sources (e.g. ammonium) are less efficient for growth than amino acids, however, we chose an inorganic nitrogen source to allow tighter control of the nitrogen to carbon ratio in the medium. Using batch culture, we identified glycerol as a suitable carbon source, with a high cell yield and no change in external pH during the course of the experiment compared to other carbon sources such as glucose, acetate or succinate (Additional file 1, Figure S1A). When cells were grown in batch culture on glycerol (25 mM) under nitrogen-replete conditions, growth ceased at an optical density (OD600) of 3.5 (Additional file 1, Figure S1B). When the culture was limited for nitrogen (carbon excess), the final OD600 was 2.1 (Additional file 1, Figure S1B). Based on these experiments we used 30 mM glycerol and 10 mM NH4Cl (nitrogen-replete)/1.25 mM NH4Cl (nitrogen-depleted) for continuous culture studies.
In continuous culture, M. smegmatis was grown at 50 % air saturation and at a dilution rate of 0.12 h−1 until steady-state conditions were achieved (typically after 3–4 resident times). Nitrogen-limited cultures (1.25 mM NH4Cl/30 mM glycerol) reached steady-state at an OD600 of 0.8 and no residual ammonium could be detected (Fig. 1a, top). The residual concentration of glycerol in the medium was approximately 18 mM indicating cells were not carbon limited (Fig. 1a, middle). At this stage, addition of 25 mM NH4Cl into the culture vessel led to an immediate increase in optical density, demonstrating that the cultures were indeed limited for nitrogen (Fig. 1a, top). Addition of 10 mM glycerol to a nitrogen-depleted culture had no effect (Fig. 1a, bottom). In the carbon-limited culture (10 mM NH4Cl/30 mM glycerol), the steady-state-phase was reached after 56 h at an OD600 of 4.3 with the residual ammonium concentration in the medium remaining constant at 0.5-1 mM (Fig. 1b, middle). A pulse of 1 mM NH4Cl into the culture vessel had no effect on growth (Fig. 1b, bottom), while the addition of glycerol led to a significant increase in OD600, demonstrating nitrogen excess (replete) under these conditions (Fig. 1b, top).
For each condition, three independent biological replicates were grown until steady-state-phase and cells were harvested after at least 3 volume changes in the chemostat vessel. The analysis of the transcriptome using RNA-sequencing revealed a total of 357 differentially expressed genes under nitrogen limitation [>2-fold change, false discovery rate (FDR) was <5 %] i.e. approximately 5 % of the M. smegmatis genome (Additional file 2, Table S1). Moreover, the usage of different thresholds (i.e. 2-fold, <1 % FDR; 2-fold, <0.1 % FDR; 1.5-fold, <0.1 % FDR) had no major effect on the number of differentially expressed genes. In total, 208 genes were upregulated and 149 were downregulated. Classification of these genes into functional categories showed major changes in transport proteins, genes that are associated with nitrogen and amino acid metabolism and genes that are functionally assigned as regulatory proteins. A large number of genes encoding for hypothetical proteins were differentially expressed (36 upregulated; 38 downregulated), however, their function in response to nitrogen depletion needs to be investigated (Additional file 3, Figure S2).
Nitrogen limitation studies in continuous culture versus batch culture
Previous work published by Williams et al. focused on the response of M. smegmatis to nitrogen stress in batch culture nitrogen run out experiments and showed differential expression of 1090 genes (574 upregulated, 516 downregulated) (Fig. 2) . Surprisingly, 903 of these genes were not differentially expressed in our continuous culture. In fact, only 17 % of the genes reported by Williams et al. responded to nitrogen limitation in continuous culture, including 70 genes that were predicted to be under control of GlnR (Fig. 2a) . Despite the differences in methodology between the two studies, we were able to identify a significant overlap in nitrogen metabolism related genes. For example, expression of a similar set of ammonium and nucleotide uptake systems (Table 1) and metabolic pathways (Table 2) were elevated in both studies. In addition, we identified a unique set of genes involved in the metabolism of amino acids, nucleotides and urea that were nitrogen-responsive in continuous culture (Table 2). We identified 36 additional genes, including those encoding for regulatory enzymes, which showed an inversed expression profile in continuous culture compared to batch culture (Additional file 4, Table S2). In batch culture it is often difficult to assign transcriptional changes to a single stimulus due to changes in growth rate, nutrient depletion and end-product buildup. It is striking that there was very little overlap in the downregulated genes in response to nitrogen limitation between batch culture and continuous culture (Fig. 2c). A major portion of the 497 downregulated genes in batch culture was associated with the general reduction in cellular metabolism due to a reduced growth rate. Development of a nitrogen-limited continuous culture enabled us to define transcriptional changes purely in response to nitrogen limitation excluding growth rate effects. We showed that using continuous culture we can identify the set of genes involved in nitrogen uptake and metabolism that were reported by Williams et al., however, this system allowed us to further extend our knowledge towards the transcriptional response to nitrogen limitation and reduce the number of genes responding to other environmental factors .
Nitrogen limitation activates the expression of genes involved in scavenging nitrogen sources in the environment
In this study, a total of 53 genes encoding for transporters or corresponding binding proteins were differentially expressed in response to nitrogen limitation (Additional file 2, Table S1). Of the upregulated genes, 18 encoded for transporters that are generally involved in uptake of both organic and inorganic nitrogen-containing compounds like amino acids (e.g. msmeg_2525), nucleotides (e.g. msmeg_5730) and ammonium (e.g. msmeg_4635) (Table 1). Previously published data by Berney and Cook showed downregulation of nine of these transporters in a carbon-limited continuous culture comparing slow growth (td 69 h) versus fast growth (td 4.6 h) . All three ammonium transporters and several amino acid (3 of 8) and nucleotide transporters (3 of 7) showed a significant decrease in transcript level, however transcription of all of these genes was enhanced upon nitrogen limitation. The GlnR regulon comprises 15 of these nitrogen compound transporters including the three ammonium transporters and the majority of nucleotide and amino acid transporters. The most prominent amino acid permease (msmeg_2525) was not under direct control of GlnR. In continuous culture, we identified 4.8-fold upregulation of a polar amino acid ABC transport system that has not been reported to be affected by nitrogen limitation. Further work is required, to identify the regulatory mechanisms that mediate the uptake of amino acids under nitrogen depletion. We observed a decrease in gene expression for the catabolism of five amino acids (Table 2), suggesting M. smegmatis attempts to balance the availability of ammonium and particular amino acids under nitrogen depletion. This finding is supported by a recent report of a proteasome-mediated amino acid recycling in M. smegmatis under nutrient limitation . The transcriptomic response of M. smegmatis to increase its ability for the uptake (scavenging) of nitrogen compounds is remarkably different from other actinobacteria like C. glutamicum, where no amino acid permeases were upregulated under nitrogen limitation .
Peptides appear to play an important role in replenishing the intracellular ammonium pool necessary for the different anabolic pathways in M. smegmatis (Table 2), (Fig. 3). Amongst the genes induced in response to nitrogen limitation we identified a peptidase, an aminopeptidase and an aliphatic amidase catalyzing the successive degradation of peptides to amino acids and the subsequent recovery of ammonium (Table 2). Furthermore, expression of a gene encoding the dipeptide transporter DppB (msmeg_1085) was induced 7.7-fold (FDR < 0.1 %). The dipeptide transporter facilitates the import of dipeptides into the cytoplasm, where they can be hydrolyzed in order to replenish the intracellular amino acid pool or serve for incorporation into proteins . The importance of (di-)peptides as nitrogen donors in mycobacteria remain to be explored, whereas previous work suggested that peptides cannot be used as carbon source .
Strong induction of nucleotide catabolism and recovery of ammonium
In our analysis we identified a novel set of genes that were significantly upregulated upon nitrogen limitation and appeared to be involved in the breakdown of nucleotides and inorganic nitrogen sources (Table 2). Nucleotides contain a high level of nitrogen and their degradation plays an important role in nitrogen metabolism in a large number of different microorganisms (e.g. Bacillus subtilis), where purine catabolic pathways are described as an alternative pathway of nitrogen utilization once primary nitrogen sources are exhausted . Genes encoding for proteins involved in the degradation of purine nucleotides were previously reported to be upregulated in M. smegmatis during nitrogen limitation . The enzyme guanine deaminase is upregulated 3.7-fold (FDR < 0.1 %) catalyzing the reaction from guanine to xanthine and releasing 1 mol of ammonium per 1 mol guanine. The three enzymes, uricase (up 2.8-fold, FDR < 0.1 %), transthyretin (up 3.5-fold, FDR < 0.1 %) and 2-oxo-4-hydroxy-4-carboxy-5-ureidoimidazoline decarboxylase (up 3.9-fold, FDR < 0.1 %) catalyze the subsequent steps of urate degradation resulting in the production of allantoin. This metabolite is then further degraded with the concomitant production of ureidoglycolate and urea, which is broken down in either a one-step reaction via urease or a two-step reaction via urea carboxylase/allophanate hydrolase enzyme complex to ammonium (Fig. 3). Transcription of the urease-encoding genes was elevated also in the batch culture, while genes encoding for an allophanate hydrolase (msmeg_0435-msmeg_0436) were upregulated only under nitrogen depletion in continuous culture .
We observed an upregulation of the pyrimidine nucleotide degradation gene cluster and detected differential expression of a third copy of a putative phenylhydantoinase (msmeg_3553; up 3-fold, FDR < 0.1 %) that plays an important role in the degradation of pyrimidine nucleotides. Catabolism of pyrimidine nucleotides can occur through three different pathways; a reductive pathway in which uracil is reduced to two molecules of CO2 and one molecule of ammonium and β-alanine, respectively, an oxidative pathway and a recently identified Rut pathway . M. smegmatis harbors a reductive uracil degradation pathway , which is induced upon nitrogen limitation. However, not all of the genes in this pathway have been identified in M. smegmatis. After the import of uracil by several nucleotide permeases, the first step of this reductive pathway is catalyzed by a dihydropyrimidine dehydrogenase that is able to reduce uracil and thymine to 5,6-dihydrouracil and 5,6-dihydrothymine (Fig. 3), followed by a reaction catalyzed by a phenylhydantoinase. This enzyme hydrolyzes the opening of the heterocyclic ring and β-alanine synthase (msmeg_3555; up 1.9-fold, FDR < 1 %) is catalyzing the last step of this pathway (Fig. 3).
Amino acid catabolism is strongly downregulated under nitrogen limitation
A large cohort of genes, particularly in the catabolism of the amino acids alanine, aspartate, valine, proline and lysine were downregulated 5-40-fold in response to nitrogen depletion (Table 2), (Fig. 3). Interestingly, most of these genes showed an elevated transcription in batch culture (Additional file 5, Table S3). The amino acid alanine has a substantial importance as a central metabolite in bacterial metabolism and also its role in the synthesis of peptidoglycan as D-/L-alanine is essential for bacteria. Several pathways have been described contributing to the biosynthesis of L-alanine and β-alanine, e.g. catabolism of amino acids (i.e. valine and cysteine) and transamination of pyruvate. In M. smegmatis, genes encoding for all of the aforementioned pathways were differentially expressed in response to nitrogen limitation (Table 2). The conformation of the amino acid β-alanine does not allow its incorporation into proteins, but it serves together with pantoate as precursor of coenzyme A (CoA) biosynthesis, which is essential for a functional TCA cycle, as well as fatty acid and cholesterol biosynthesis. Degradation of purine nucleotides via an L-aspartate-alpha-decarboxylase (PanD) results in production of β-alanine and PanD was identified as the predominant pathway of β-alanine synthesis in the closely related C. glutamicum, where a panD mutant exhibited β-alanine auxotrophy . However, expression of panD was strongly downregulated in M. smegmatis under nitrogen limitation. Valine degradation via the intermediates 3-methyl-2-oxobutanoate and 2-dehydropantoate to (R)-pantoate was repressed the same time, suggesting the demand to prevent unnecessary consumption of amino acids for CoA biosynthesis (Fig. 3).
A second pathway of L-aspartate (Asp) catabolism was differentially expressed in M. smegmatis under nitrogen limitation, which is linked to the concomitant biosynthesis of lysine. This pathway is a nine-step reaction including important metabolites such as L-aspartate semialdehyde (homoserine biosynthesis) and meso-2,6-diaminopimelate (constituent of bacterial cell walls). Interestingly, the initial steps of aspartate catabolism were repressed, while the degradation of meso-2,6-diaminopimelate to lysine was upregulated (Fig. 3). Lysine can act as donor of an amino group by transferring an ammonium group to 2-oxoglutarate to form glutamate under nitrogen excess, however, this pathway of lysine catabolism via a lysine aminotransferase is downregulated 7.5-fold (FDR < 1 %), indicating an intracellular accumulation of aspartate and lysine and suggesting a secondary function of these amino acids [32, 33]. Previous studies discussed the importance of intracellular lysine to control growth rate in mycobacteria and suggested a link between lysine accumulation and fatty acid metabolism . Proline has been shown to serve as mechanism for methylglyoxal detoxification, when anabolic and catabolic processes were imbalanced and we show that proline degradation was repressed under nitrogen depletion.
A large number of transcriptional regulatory systems are differentially expressed in response to nitrogen limitation
We identified 26 differentially expressed transcriptional regulators that are either directly or indirectly responding to nitrogen limitation in M. smegmatis (Table 3). Only a handful of these regulators have been characterized, including the nitrogen regulatory protein PII (msmeg_2426) and the PII adenylyl transferase (msmeg_2427). The mycobacterial copy of the PII protein is not required for the regulation of the glutamine synthetase activity and does not act as regulator of the transcriptional response to nitrogen limitation . This is in contrast to C. glutamicum, where the PII protein was identified as the sole signal transduction protein, binding to AmtR and releasing this repressor from its target DNA, in order to allow transcription of genes involved in nitrogen uptake, assimilation and metabolism .
Another well-described transcriptional regulator is the OmpR-type response regulator GlnR, which has been identified as a mediator of the transcriptomic response to nitrogen limitation in M. smegmatis . Determination of the GlnR regulon, by combining expression profiling of M. smegmatis wild type and a ∆glnR deletion mutant under nitrogen-starvation conditions in batch culture revealed a total of 103 genes directly controlled by GlnR . A large portion of these genes (72 %) were also differentially expressed under nitrogen limitation in our continuous culture model, however, GlnR itself was not among them, indicating a posttranslational regulatory mechanism for GlnR like in S. coelicolor . Comparison of the differentially expressed transcriptional regulatory proteins in batch culture versus continuous culture showed an overlap of six genes with four of these under control of GlnR (Table 3). We further observed 16 transcriptional regulatory proteins that were differentially expressed in continuous culture (Table 3). The TetR-like transcriptional regulator AmtR showed a 1.7-fold upregulation under nitrogen replete conditions and was therefore outside the selected cut-off. However, in our RNA-seq analysis we identified 127-fold downregulation of an antisense transcript (asRNA) of almost the entire msmeg_4300 gene under nitrogen-depleted conditions (Additional file 6, Figure S3). These data suggest a mechanism of post-transcriptional regulation by an asRNA where transcription of msmeg_4300 is enhanced upon nitrogen excess and modulates the translation efficiency of the AmtR encoding gene.
We identified several other cases of potential non-coding small RNAs that might regulate expression of genes involved in fatty acid and central carbon metabolism under nitrogen limitation (Additional file 7, Table S4). In our RNA-seq analysis, we observed a 102-fold downregulation of an asRNA of msmeg_4299 gene under nitrogen depletion, while msmeg_4299 showed a 2.7-fold upregulation under nitrogen limitation (Additional file 6, Figure S3). The intergenic region msmeg_3131-msmeg_3132 comprises the promotor regions of the genes msmeg_3131 (long-chain acyl-CoA synthetase) and msmeg_3132 (DNA-binding protein) and revealed a 4.7-fold increase in transcription under nitrogen excess (Additional file 6, Figure S3). Interactions between sRNAs and the 5’UTR of mRNAs affect translation efficiency, while an interaction with the 3’UTR of mRNAs usually do not affect translation, but mRNA stability . These findings are supported by a strong upregulation of genes involved in the degradation of fatty acids as well as the downregulation of genes involved in the TCA cycle as both metabolic pathways are CoA-dependent (Additional file 2, Table S1).
Interestingly, three regulatory proteins (one putative two-component system and two putative orphan response regulators) showed a different expression profile in continuous culture compared to batch culture. The operon msmeg_6236-msmeg_6238 encodes for a two-component system with unknown function and the gene msmeg_3297 encodes for a CadC-like transcriptional regulator that has been linked to pH homeostasis in E. coli . No change in external pH was observed in our experiments, suggesting the CadC-like protein in M. smegmatis was performing a different function.
Herein we report the transcriptomic response of M. smegmatis to nitrogen limitation in a continuous culture model at a defined growth rate (Fig. 4). We show that amino acid metabolism plays an important role in the adaptation of M. smegmatis to nitrogen depletion and identified 16 novel transcriptional regulators that were either directly or indirectly involved in the global transcriptomic response of M. smegmatis to nitrogen limitation. Several non-coding RNAs were differentially expressed suggesting transcriptional or post-transcriptional regulation of gene expression and we propose a regulatory mechanism involving a trans-acting asRNA for the AmtR protein in M. smegmatis. Comparison of our chemostat data to the global transcriptomic response in batch culture nitrogen run out experiments revealed that only 17 % of the previously described nitrogen-regulated genes overlapped between batch and continuous culture. Our data highlight the versatile metabolic capability of M. smegmatis and provide a molecular framework for understanding how environmental mycobacteria respond to nitrogen-depleted environments.
Bacterial strains, media and growth conditions
M. smegmatis mc2155 was grown in LBT or adapted HdB minimal medium supplemented with 0.2 % (w/v) glycerol (unless otherwise stated) as sole carbon source, NH4Cl in various concentrations as sole nitrogen source, 11.74 mM K2SO4 as sulphur source and 0.05 % (w/v) Tyloxapol at 37 °C with agitation (200 rpm). Aerobic starter cultures were inoculated to an initial optical density (OD600) of 0.05 and grown with agitation (200 rpm) at 37 °C. Chemostat bioreactors were prepared as previously described . Samples to measure residual glycerol concentration and ammonium concentration were taken in intervals of four to six hours. Glycerol concentration was measured according to Garland and Randle  and ammonium concentration was measured according to Weatherburn . After entering steady-state, the chemostat was left running for at least three volume changes before cell harvest.
RNA extraction and Reverse transcriptase PCR
Cells were lysed by four cycles at 4800 rpm for 30s in a mini-Beadbeater (Biospec Products), with 30s on ice between each of the cycles. Total RNA was extracted using TRIzol® reagent (Ambion) according to the manufacturers instructions. DNA was removed by treatment with 3 U RNase-free DNase using the TURBO DNA-free kit (Ambion) according to the manufacturers instructions. The quality of the RNA was checked with the Bioanalyzer (RIN >9) and the concentration was determined using a NanoDrop ND-100 spectrophotometer.
Next, depletion of ribosomal RNA and strand-specific library preparation were performed, using the EpiCentre ScriptSEq™ Complete Kit for Bacteria according to manufacturers instructions. Quantification of nucleic acids was performed using a Qubit® Fluorometer to ensure DNA contamination of the samples was less than 10 % and libraries were created. Prepared libraries have undergone a quality control using an Agilent 2100 Bioanalyzer: DNA 1000 Labchip, Quant-iT dSDNA HS Assay for quantification and Quant-iT RNA Assay and Quant-iT Protein Assay for percentage contamination check, using an Invitrogen Qubit® Fluorometer. Libraries were run on a Illumina MiSeq 300 cycle Kit_v2 with a paired-end (PE) read length of 2x150 and a PhiX control library was also loaded and used as control for the run.
Analysis of RNA sequencing data
Adapter sequences were removed from raw fastq files using Flexbar  and reads shorter than 50 bp were discarded. Sequences were then mapped against the M. smegmatis genome (GenBank NC_008596.1) using Bowtie2 with the”very sensitive” option. Counts for each gene, intergenic regions and counts on reverse strand were calculated with featureCounts , taking into account the reads’ strand-specificity. Statistical analysis and principal component analysis were performed using DESeq (Additional file 8, Figure S4 and Additional file 9, Figure S5) . We considered differentially expressed genes based on their fold-change values as well as p values adjusted for multiple testing with the Benjamini-Hochberg procedure, which controls false discovery rate (FDR, referred to as padj in DESeq results). Venn diagrams were generated using the BioVenn program . Gene functions were assigned using public databases (NCBI , UniProt , KEGG , Biocyc Database Collection ).
All RNA-sequencing data have been deposited in ArrayExpress and can be accessed through the accession number E-MTAB-3918. All other supporting data are included as additional files.
Merrick MJ, Edwards RA. Nitrogen control in bacteria. Microbiol Rev. 1995;59:604–22.
Reitzer L. Nitrogen assimilation and global regulation in Escherichia coli. Annu Rev Microbiol. 2003;57:155–76.
Leigh JA, Dodsworth JA. Nitrogen regulation in bacteria and archaea. Annu Rev Microbiol. 2007;61:349–77.
Ninfa AJ, Jiang P. PII signal transduction proteins: sensors of alpha-ketoglutarate that regulate nitrogen metabolism. Curr Opin Microbiol. 2005;8:168–73.
Engleman EG, Francis SH. Cascade control of Escherichia coli glutamine synthetase. II. Metabolite regulation of the enzymes in the cascade. Arch Biochem Biophys. 1978;191:602–12.
Arcondeguy T, Jack R, Merrick M. PII signal transduction proteins, pivotal players in microbial nitrogen control. Microbiol Mol Biol Rev. 2001;65:80–105.
Sant’Anna FH, Trentini DB, de Souto WS, Cecagno R, da Silva SC, Schrank IS. The PII superfamily revised: a novel group and evolutionary insights. J Mol Evol. 2009;68:322–36.
Jiang P, Peliska JA, Ninfa AJ. Enzymological characterization of the signal-transducing uridylyltransferase/uridylyl-removing enzyme (EC 184.108.40.206) of Escherichia coli and its interaction with the PII protein. Biochemistry. 1998;37:12782–94.
Keener J, Kustu S. Protein kinase and phosphoprotein phosphatase activities of nitrogen regulatory proteins NTRB and NTRC of enteric bacteria: roles of the conserved amino-terminal domain of NTRC. Proc Natl Acad Sci U S A. 1988;85:4976–80.
Hisbergues M, Jeanjean R, Joset F, de Marsac NT, Bédu S. Protein PII regulates both inorganic carbon and nitrate uptake and is modified by a redox signal in Synechocystis PCC 6803. FEBS Lett. 1998;463:216–20.
Burkovski A. Ammonium assimilation and nitrogen control in Corynebacterium glutamicum and its relatives: an example for new regulatory mechanisms in actinomycetes. FEMS Microbiol Rev. 2003;27:617–28.
Beckers G, Strösser J, Hildebrandt U, Kalinowski J, Farwick M, Krämer R, et al. Regulation of AmtR-controlled gene expression in Corynebacterium glutamicum: mechanism and characterization of the AmtR regulon. Mol Microbiol. 2005;58:580–95.
Tiffert Y, Supra P, Wurm R, Wohlleben W, Wagner R, Reuther J. The Streptomyces coelicolor GlnR regulon: Identification of new GlnR targets and evidence for a central role of GlnR in nitrogen metabolism in actinomycetes. Mol Microbiol. 2008;67:861–80.
Fink D, Weissschuh N, Reuther J, Wohlleben W, Engels A. Two transcriptional regulators GlnR and GlnRII are involved in regulation of nitrogen metabolism in Streptomyces coelicolor A3(2). Mol Microbiol. 2002;46:331–47.
Williams KJ, Jenkins VA, Barton GR, Bryant WA, Krishnan N, Robertson BD. Deciphering the metabolic response of Mycobacterium tuberculosis to nitrogen stress. Mol Microbiol. 2015. doi:10.1111/mmi.13091.
Gouzy A, Larrouy-Maumus G, Wu T-D, Peixoto A, Levillain F, Lugo-Villarino G, et al. Mycobacterium tuberculosis nitrogen assimilation and host colonization require aspartate. Nat Chem Biol. 2013;9:674–6.
Gouzy A, Larrouy Maumus G, Bottai D, Levillain F, Dumas A, Wallach JB, et al. Mycobacterium tuberculosis exploits asparagine to assimilate nitrogen and resist acid stress during infection. PLoS Pathog. 2014;10:e1003928.
Gouzy A, Poquet Y, Neyrolles O. Nitrogen metabolism in Mycobacterium tuberculosis physiology and virulence. Nat Rev Microbiol. 2014;12:729–37.
Amon J, Titgemeyer F, Burkovski A. A genomic view on nitrogen metabolism and nitrogen control in mycobacteria. J Mol Microbiol Biotechnol. 2009;17:20–9.
Jenkins VA, Barton GR, Robertson BD, Williams KJ. Genome wide analysis of the complete GlnR nitrogen-response regulon in Mycobacterium smegmatis. BMC Genomics. 2013;14:301.
Williams KJ, Bryant WA, Jenkins VA, Barton GR, Witney AA, Pinney JW, et al. Deciphering the response of Mycobacterium smegmatis to nitrogen stress using bipartite active modules. BMC Genomics. 2013;14:436.
Hartmans S, Bont JAM, Stackebrandt E. The genus mycobacterium-nonmedical. In: The prokaryotes third edition: archaea bacteria: firmicutes, actinomycetes. 3rd ed. New York, NY: Springer New York; 2006. p. 889–918.
Berney M, Cook GM. Unique flexibility in energy metabolism allows mycobacteria to combat starvation and hypoxia. PLoS One. 2010;5:e8614.
Elharar Y, Roth Z, Hermelin I, Moon A, Peretz G, Shenkerman Y, et al. Survival of mycobacteria depends on proteasome-mediated amino acid recycling under nutrient limitation. EMBO J. 2014;33:1802–14.
Silberbach M, Hüser A, Kalinowski J, Pühler A, Walter B, Krämer R, et al. DNA microarray analysis of the nitrogen starvation response of Corynebacterium glutamicum. J Biotechnol. 2005;119:357–67.
Zimmer DP, Soupene E, Lee HL, Wendisch VF, Khodursky AB, Peter BJ, et al. Nitrogen regulatory protein C-controlled genes of Escherichia coli: Scavenging as a defense against nitrogen limitation. Proc Natl Acad Sci USA Vol. 2000;97:14674–9.
Green RM, Seth A, Connell ND. A peptide permease mutant of Mycobacterium bovis BCG resistant to the toxic peptides glutathione and S-nitrosoglutathione. Infect Immun. 2000;68:429–36.
Schultz AC, Nygaard P, Saxild HH. Functional analysis of 14 genes that constitute the purine catabolic pathway in Bacillus subtilis and evidence for a novel regulon controlled by the PucR transcription activator. J Bacteriol. 2001;183:3293–302.
Loh KD, Gyaneshwar P, Markenscoff Papadimitriou E, Fong R, Kim K-S, Parales R, et al. A previously undescribed pathway for pyrimidine catabolism. Proc Natl Acad Sci U S A. 2006;103:5114–9.
Vogels GD, Van der Drift C. Degradation of purines and pyrimidines by microorganisms. Bacteriol Rev. 1976;40:403–68.
Dusch N, Pühler A, Kalinowski J. Expression of the Corynebacterium glutamicum panD gene encoding L-aspartate-alpha-decarboxylase leads to pantothenate overproduction in Escherichia coli. Appl Environ Microbiol. 1999;65:1530–9.
Csonka LN. Physiological and genetic responses of bacteria to osmotic stress. Microbiol Rev. 1989;53:121–47.
Berney M, Weimar MR, Heikal A, Cook GM. Regulation of proline metabolism in mycobacteria and its role in carbon metabolism under hypoxia. Mol Microbiol. 2012;84:664–81.
Griffin JE, Pandey AK, Gilmore SA, Mizrahi V, McKinney JD, Bertozzi CR, et al. Cholesterol catabolism by Mycobacterium tuberculosis requires transcriptional and metabolic adaptations. Chem Biol. 2012;19:218–27.
Williams KJ, Bennett MH, Barton GR, Jenkins VA, Robertson BD. Adenylylation of mycobacterial Glnk (PII) protein is induced by nitrogen limitation. Tuberculosis. 2013;93:198–206.
Repoila F, Darfeuille F. Small regulatory non-coding RNAs in bacteria: physiology and mechanistic aspects. Biol Cell. 2009;101:117–31.
Haneburger I, Eichinger A, Skerra A, Jung K. New insights into the signaling mechanism of the pH-responsive, membrane-integrated transcriptional activator CadC of Escherichia coli. J Biol Chem. 2011;286:10681–9.
Garland PB, Randle PJ. A rapid enzymatic assay for glycerol. Nature. 1962;196:987–8.
Weatherburn MW. Phenol-Hypochlorite reaction for determination of ammonia. Anal Chem. 1967;39:971–4.
Dodt M, Roehr JT, Ahmed R, Dieterich C. FLEXBAR-flexible barcode and adapter processing for Next-Generation Sequencing platforms. Biol (Basel). 2012;1:895–905.
Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014;30:923–30.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.
Hulsen T, de Vlieg J, Alkema W. BioVenn - a web application for the comparison and visualization of biological lists using area-proportional Venn diagrams. BMC Genomics. 2008;9:488.
NCBI. Resource coordinators: database resources of the national center for biotechnology information. Nucleic Acids Res. 2015;43(Database issue):D6–D17.
UniProt Consortium. UniProt: a hub for protein information. Nucleic Acids Res. 2015;43(Database issue):D204–12.
Kanehisa M, Goto S, Sato Y, Kawashima M, Furumichi M, Tanabe M. Data, information, knowledge and principle: back to metabolism in KEGG. Nucleic Acids Res. 2014;42(Database issue):D199–205.
Caspi R, Altman T, Billington R, Dreher K, Foerster H, Fulcher CA, et al. The MetaCyc database of metabolic pathways and enzymes and the BioCyc collection of pathway/genome databases. Nucleic Acids Res. 2014;42(Database issue):D459–71.
This work was supported by the Maurice Wilkins Centre for Molecular Biodiscovery and a Marsden Grant from the Royal Society of New Zealand. G.M.C. was supported by a James Cook Fellowship from the Royal Society of New Zealand. MP was supported by a University of Otago Doctoral Scholarship, the Webster Center for Infectious Diseases and the Otago School of Medical Sciences. We thank Claudia Sala and Stewart Cole for critical reading of the manuscript.
The authors declare that they have no competing interests.
MP designed and performed experiments, analyzed and interpreted the data and wrote the manuscript. AB performed the bioinformatics analysis of the RNA-sequencing data and helped with the data analysis. GMC supervised the project, performed data interpretation and wrote the manuscript. All authors have approved the final manuscript for publication.
Growth of M. smegmatis in batch culture with different carbon sources under nitrogen-depleted and nitrogen-replete conditions. (PDF 264 kb)
A complete list of genes that are differentially expressed in response to nitrogen limitation in M. smegmatis. (XLSX 2413 kb)
Functional characterization of differentially expressed genes. (PDF 272 kb)
Comparison of the gene expression profile between nitrogen-depleted continuous culture versus batch culture of M. smegmatis. (XLSX 88 kb)
A list of genes that are differentially expressed in continuous culture, but not in batch culture. (XLSX 89 kb)
Images of IGV files from differentially expressed intergenic regions and antisense-strand reads mentioned in the text. (PPTX 495 kb)
A list of differentially expressed intergenic and antisense-strand regions. (XLSX 12 kb)
Data quality assessment by sample clustering and visualization. (PPTX 147 kb)
Principal component plot of the samples as calculated from the variance stabilising transformation of the count data using DESeq and heatmap showing the Euclidean distances between the samples from the variance stabilising transformation of the count data using DESeq. (PPTX 72 kb)
About this article
Cite this article
Petridis, M., Benjak, A. & Cook, G.M. Defining the nitrogen regulated transcriptome of Mycobacterium smegmatis using continuous culture. BMC Genomics 16, 821 (2015). https://doi.org/10.1186/s12864-015-2051-x
- Nitrogen depletion
- Continuous culture
- Gene expression profiling