Mycobacterium tuberculosis H37Rv was grown in Dubos medium (Difco Dubos Broth Base, USA) and supplemented with either: (1) 0.2% dextrose (D), or (2) long chain fatty acids (F) (oleic acid, stearic acid and palmitic acid at a final concentration of 0.001% each) as main carbon sources at 37°C. Exponential phase cultures at an OD₆₀₀ = 0.4 were obtained at day 7 (with dextrose) and day 8 (with FA), according to the growth curve published previously by our group (Rodríguez et al., 2014). Exponential cultures were submitted to hypoxic conditions by leaving a ratio of air volume to liquid medium of 0.5 in flasks with tightly sealed caps. Cultures were agitated with a 50 mm teflon-coated magnetic stirrer bar at 70 rpm, as previously reported (Wayne and Hayes, 1996). NRP1 and NRP2 states were defined, according to Wayne and Hayes (1996), with a parallel culture supplemented with methylene blue (1.5 μg/mL), which indicates approximate oxygen depletion. In this way, the NRP1 phase (fading of methylene blue) was reached at day four in both, fatty acids and dextrose medium (FNRP1 and DNRP1); and NRP2 (complete decolorization of methylene blue) was reached at day eight (FNRP2 and DNRP2). RNA isolation and CFU/ml quantification were carried out 24 h after reaching the corresponding NRP phase.

RNA from each growing condition (exponential, NRP1 and NRP2) was isolated as previously described (Rodríguez et al., 2014). Briefly, cells were harvested and resuspended in guanidium chloride buffer (6 M guanidinium chloride, 0.1% Tween 80, 1 mM 2-mercaptoethanol, 10 mM EDTA) and lysed mechanically in a FastPrep (Thermo Scientific). Nucleic acids were purified using phenol- chloroform- isoamyl alcohol and RNA was selectively precipitated with absolute ethanol and washed three times with Trizol reagent (Invitrogen). RNA integrity was estimated with bioanalyzer (Agilent Technologies) and quantified by Nano-drop ND 1000 (Thermo Scientific). RNASeq libraries were prepared as previously describe (Rodríguez et al., 2014). Briefly, total RNA was fragmented, submitted to poly A tailing and end repairing (NEB reagents, USA). Strand-specific labeling was carried out by ligation of a 5′ hybrid DNA-RNA primer. Reactions for the synthesis of the first cDNA strand was carried out with Super Script II reverse transcriptase (Invitrogen, USA), Illumina's poly T primer and dNTPs (20 mM). Illumina adaptors and barcodes were ligated by PCR following manufacturer instructions and pair end sequencing was performed on Illumina HiSeq sequencer.

After trimming, those sequences with quality scores over 30 (Illumina 1.3+) and a minimal length of 50 bases were used for the analysis. Quality assessment of the reads was done using the FASTX toolkit v. 0.0.13 (http://hannonlab.cshl.edu/fastx_toolkit/index.html). The reads were mapped against the Mtb H37Rv reference genome (accession number NC 000962.2) using Bowtie v. 0.12.7 (Langmead et al., 2009). Resulting SAM files were used with the Tophat package (Langmead et al., 2009) to compute gene expression, differential gene expression and statistical significance. Gene expression was computed as RPKM (Reads Per Kilobase Million of mapped Reads). Statistical significance was determined on a binary basis between conditions by False Discovery Rate (FDR). We defined genes with a FDR of <0.05 as over-expressed.

The total transcripts of selected genes were measured by real-time qRT-PCR in a LightCycler 480 Instrument (Roche). Quantification was performed with gene-specific primers (see Supplementary Table 1) and SYBR green (Molecular Probes, Inc.). Samples were subjected to 40 cycles of amplification (denaturation at 95°C for 30 s, specific annealing temperature for 20 s, and extension at 72°C for 30 s) with a final extension at 72°C for 5 min. To ensure that the fluorescence levels detected were due to the amplification of a specific product, a melting curve analysis was performed. Absolute quantification was carried out by obtaining a standard curve for each set of primers according to 10-fold dilutions of by known amounts of Mtb H37Rv chromosomal DNA (1,000, 10,000, 100,000, and 1,000,000 theoretical copies). Crossing point values were interpolated to standard curve to obtain gene expression (number of gene copies per μg of RNA). Normalization of these data was performed by using 16S rRNA expression levels.

The RNAseq data sets have been deposited in NCBI's Gene Expression Omnibus (Edgar et al., 2002) under accession number GSE119225.

To gain insight into the adaptation of Mtb to hypoxia in a lipid environment we performed a transcriptomic analysis of Mtb following the Wayne and Hayes model (Wayne and Hayes, 1996) with growing conditions that included dextrose or LC-FA as main carbon sources (Rodríguez et al., 2014). NRP1 phase was reached at day four after exponential phase in both fatty acids and dextrose medium (FNRP1 and DNRP1); NRP2 phase (complete fading of methylene blue) was reached at day eight (FNRP2 and DNRP2) (see Materials and Methods). Table 1 shows the global results obtained for each condition. The sequences retained after cleaning were between 10.34 and 18.31 million reads, which cover the Mtb genome as indicated by the plateau reached in the saturation curves when mapping the reads against the Mtb genome (data not shown). To facilitate the analysis, and allow the comparison with previous studies, the data were normalized as reads per kilobase per million reads (RPKMs). As expected, we observed a reduction in about 20 to 35% in reads mapping to coding sequences (CDS) in both, the dextrose hypoxia model (DNRP1 and DNRP2) and the fatty acid hypoxia model (FNRP1 and FNRP2) (Table 1). This reduction was accompanied with an increase of reads of intergenic regions (IGRs). The high mapping reads against IGRs in the presence of fatty acids, especially in FNRP2, is similar to the IGRs expression observed in the stationary phase of bacilli growing in LC-FA (Rodríguez et al., 2014) when bacteria develop the dormant phenotype.

We next searched for parallels between our data and that of other hypoxia models, including the DosR and EHR (Voskuil et al., 2004; Rustad et al., 2008) as well as a persistence condition developed in the presence of antibiotics (Keren et al., 2011) (Figure 1). For these comparisons, we used data from early hypoxic stages (DNRP1 and FNRP1). Genes found to be shared between conditions are indicated in the Supplementary Table 2. As expected, data from DNRP1 had a higher number of genes in common with those from the dosR regulon (Voskuil et al., 2004) and EHR (Rustad et al., 2008), with 17 genes and 98 genes in common, respectively; Only two genes were found to be common to all these three conditions. Up to 32 genes were commonly over-expressed among DNRP1, persistence and EHR; and near 40 genes were common among persisters and either EHR or DNPR1 (Figure 1 and Supplementary Table 2).

Interestingly, we found that 5 out of 10 genes that were common among DNRP1/EHR/FNRP1 were also over-expressed in the persister cells (Supplementary Table 1). Those five genes are Rv0251c, Rv1221, Rv2050, Rv2694c, and Rv2745c encoding for the heat shock protein hsp, the alternative sigma factor SigE, the RNA polymerase-binding protein RbpA, the conserved protein Rv2694c and the transcriptional regulatory protein ClgR, respectively. These results suggest common pathways of adaptation to stresses conditions represented by hypoxia and persistence in antibiotics. Surprisingly, we could not find genes in common between FNPR1 and dosR regulon, while 11 and 21 genes were shared respectively between FNRP1 and the EHR and DNRP1 datasets. Of those, the gene Rv2137c, encoding a conserved hypothetical protein, was unique between FNRP1 and EHR dataset.

Taking into consideration that in our previous study 27 genes belonging to the dosR regulon were found to be over-expressed in the stationary phase cultures in the presence of LC-FAs (Rodríguez et al., 2014) these results suggest that the dosR regulon response was further unchanged when Mtb enter into hypoxia under our model of LC-FA lipid environment. Aguilar-Ayala and co-workers found similar results using the same LC-FA condition plus cholesterol as carbon source (Aguilar-Ayala et al., 2017) reinforcing our findings. Notably, no common genes were found between the five conditions compared (Figure 1 and Supplementary Table 2).

To determine changes in the metabolic function of Mtb during its adaptation to hypoxia, we compared functional categories of differentially expressed genes by Mtb grown in the three conditions including dextrose (namely, DE, DNRP1, and DNRP2) as well as the three conditions including LC-FA (FE, FNRP1, and FNRP2), as corresponding carbon sources (Figure 2 and Supplementary Table 3).

According to our results, about 10% of the PE/PPE genes showed over-expression during the entrance to early hypoxia independently of the carbon source (Figures 2A,B). Sixteen PE/PPE genes were over-expressed in DNRP1 vs. DE, meanwhile only three were over-expressed in FNRP1 vs. FE and PPE31 (Rv1807) common to both conditions. PPE31 has been identified as an essential gene for Mtb in vivo (Fishbein et al., 2015) and, according to our data, could also play some role in the adaptation to early hypoxia of the tubercle bacilli.

The genome of Mtb has 214 TFs (Cole et al., 1998), a relative high number considering the size of its genome. Recent work using ChIP-seq and over-expression of 206 TFs (Rustad et al., 2014; Minch et al., 2015) provided a comprehensive transcriptional map of Mtb. Using the data reported by these authors, we identified all the significant over-expressed TFs in the DNRP1 vs. DE and FNRP1 vs. FE (Supplementary Table 3). Then, we compared all the target genes of those TF identified using ChIP-seq by Minch's data (Minch et al., 2015) with our over-expressed genes in each of the conditions studied (Supplementary Table 4). This analysis shows that while in hypoxia in the presence of dextrose there were 37 over-expressed TFs, in the fatty acid environment only 2 TFs were over-expressed (Figures 3, 4 and Supplementary Table 4). The Rv0081 gene, a central regulatory gene in response to hypoxia in the presence of dextrose (Galagan et al., 2013), was also the regulatory master in our DNRP1 model, confirming its pivotal role in controlling Mtb adaptation to low oxygen in dextrose. However, this regulator was not over-expressed in hypoxia when LC-FA was present (See Supplementary Table 3). This result could indicate that the adaptation machinery of the bacilli is already prepared to support a hypoxic stress when Mtb uses lipids as carbon source. The high number of TFs that were upregulated in the DNRP1 vs. DE condition (Figure 3 and Supplementary Table 4) indicates the strong adaptive changes required for the tubercle bacilli upon entrance to hypoxia when dextrose is the carbon source.

Only two TFs were upregulated comparing FNRP1 vs. FE: Rv3765c (tcrX), a two component transcriptional regulator (Bhattacharya and Das, 2011) and Rv2745c (clgR), which encodes a transcriptional protein involved in the regulation of proteases and chaperons (McGillivray et al., 2014). These data demonstrate that the adaptation to hypoxia in the presence of LC-FA supposes a minor regulatory change to the tubercle bacilli. Our data suggest that the relation of ClgR with hypoxia is maintained also in the presence of lipids. A complete set of data concerning TFs and its regulated genes are showed in Supplementary Table 4.

In our study, we detected several sigma factors over-expressed during the early hypoxic stage NRP1 in cultures including both carbon sources (DNRP1 and FNRP1) relative to exponential cultures (DE and FE) (Table 2). Up to four sigma factors were over-expressed in DNRP1 compared to DE. Two of them (sigB and sigE) were also over-expressed in FNRP1 relative to FE. It can be inferred that those two sigma factors are key in the early adaptation to hypoxia of Mtb independently of the carbon source. The fact that sigE was the sigma factor with the higher level of over-expression in both hypoxia models also supports its pivotal role in the early response to low oxygen by Mtb (Table 2). In addition, no significant changes were detected in the level of expression of sigma factors when comparing early (NRP1) with late (NRP2) hypoxia conditions in both carbon sources (see Supplementary Table 3) which suggests that sigma factors participate in the initial adaptation to hypoxia of Mtb, while remaining unchanged when bacteria face a lower level of oxygen.

The high level of IGRs expression detected when the tubercle bacilli enter into hypoxia (Table 1) prompted us to analyse the contribution of small RNAs (sRNAs), including non-coding RNAs (ncRNAs) in the process. The sRNAs participating in the adaptation of Mtb to hypoxia in the two cultured conditions studied are summarized in Figure 5. Nine and 12 sRNAs showed over-expression in the Mtb adaptation to hypoxia in the presence of dextrose and LC-FA, respectively (Figure 5 and Supplementary Table 5). Independently of the carbon source, the 4.5S RNA was over-expressed in all hypoxia conditions. Two more ncRNAs showed a significant higher expression only in the presence of LC-FA: MTS1338 and MTS0194 (Figure 5). Several ncRNAs seemed to have an important role in early and late hypoxia in the presence of dextrose (DNRP1 vs. DE and DNRP2 vs. DE). These are mcr3 and MTS2822 (Figure 5). Of note, MTS2822 was also over-expressed in conditions with LC-FA as carbon source (FNRP1 vs. FE and FNRP2 vs. FE) (Figure 5 and Supplementary Table 5).

Finally, significant higher expression of some sRNAs was also related to the adaptation to late hypoxia (Figure 5). MTS2823 and the stable 10S RNA were detected with significantly high expression in NRP2 in both carbon sources. In addition, MTS0997 was highly expressed in late hypoxia only when LC-FA was the carbon source. MTS2823 is a well-known ncRNA, whose levels of expression highly increased when the bacilli enter stationary phase in LC-FA (Rodríguez et al., 2014) and has also been linked to the slow-down of the Mtb growth (Arnvig and Young, 2012).

Only two toxins were detected with significant higher expression during the adaptation of Mtb to hypoxia, with LC-FA as carbon source (FNRP1 and FNRP2; Table 3). Conversely, up to four complete TA systems, three single toxins and one single antitoxin were over-expressed when hypoxic cultures included dextrose (DNRP1 and DNRP2; Table 3). These results suggest that TA modules play a major role in the adaptation of the tubercle bacilli to hypoxia when dextrose was the carbon source but not in the presence of LC-FA. Members of the vapBC family were upregulated in dextrose, namely vapB10 and vapC37 in DNRP1, as well as the toxin vapC20 in DNRP2 (Table 3). Of note, this system has been previously associated to hypoxia in Mtb (Fernandez-Garcia et al., 2016) and might have a relevant function during late hypoxia in dextrose. The ParDE2 system was suggested to be a main participant in the inhibition of bacterial growth in the related bacteria Mycolicibacterium smegmatis (Oren and Garrity, 2018), indicating a putative role of this TA system in dormancy and stress survival (Gupta, 2009). Our results showed that ParDE1 instead or ParDE2 could play this role in the adaptation of Mtb to hypoxia in dextrose (Table 3). Two of the three members of the HigBA family (High inhibition of growth) were also over-expressed in DNRP1 compared to DE (Table 3). The higBA1 gene was detected with significant higher expression in DNRP1 compared to DE. HigBA1 is a member of the tripartite system TAC. This module contains three components, toxin, antitoxin, and a chaperone (HigBA1+Chaperone) (Sala et al., 2014). TAC is considered one of the main contributors to the Mtb survival under stress situations (Fernandez-Garcia et al., 2016), including hypoxia (Ramage et al., 2009). On the other hand, higBA2 was also upregulated in DNRP1, suggesting for the first time that adaptation to hypoxia could be a putative function for the HigBA2 system. We only detected the involvement of a single toxin member of the MazEF family in DNRP1 (Table 3) in agreement with the previous data that described the connection of these TA family with persistence in the presence of antibiotics (Singh and Barry, 2010). Contrary to the adaptation to hypoxia in the presence of dextrose, only two unclassified toxins were upregulated in LC-FA during hypoxia, it can be inferred, therefore, that TA modules do not participate in the survival of Mtb to hypoxia in this LC-FA environment.

To gain insight into the metabolic route taken for the tubercle bacilli to survive to a hypoxic stress, we analyzed the genes with higher expression shared between the four hypoxic conditions tested: NRP1 vs. exponential and NRP2 vs. NRP1 of Mtb growing in dextrose as well as growing in LC-FA (Table 4 and Supplementary Table 6). We found 47 genes with higher differential expression under hypoxia in all conditions, those genes were all over-expressed in DNRP1 in relation to DE with the exception of the gene Rv2554c; only eight genes were over-expressed in DNRP2 in relation to DNRP1 (Table 4). On the other hand, 21 and 22 genes were over-expressed in conditions comparing FNRP1 with FE and FNRP2 with FNRP1, respectively. This result again suggests that Mtb undertake a strong metabolic change upon entrance to hypoxia when dextrose is the carbon source, and a more gradual change in their adaptation to hypoxia in the presence of LC-FA. Most of the genes with significant high expression in DNPR1 were shared with any of the other conditions tested. Five genes were only detected in hypoxia when dextrose was the carbon source (DNRP1 and DNRP2) including tatB (traslocase), PE29 and PE20 (PE-PPE family proteins), alaS (Alanyl-tRNA synthetase) and the hypothetical protein coded by Rv3205 (Table 4). Interestingly, about half of the genes over-expressed in early hypoxia in dextrose (DNRP1 vs. DE) were also over-expressed in late hypoxia in the presence of LC-FA (FNRP2). One of those genes was iciA, member of the LysR family, involved in regulation of the chromosome replication. That result suggests that the level of stress found by the bacteria in dextrose at low level of hypoxia was only reached at high level of hypoxia when the bacilli grew in the lipid environment tested, in agreement with the results found in the analysis of functional categories (Figures 2A,B).

A complete set of RPKM data concerning genes included in Table 4 are showed in Supplementary Tables 3, 6.

From a global perspective we found similarities and differences in the metabolic pathways undertaken by the bacilli when comparing hypoxia in the presence of dextrose vs. conditions including LC-FA as main carbon source. In agreement with our previous results (Rodríguez et al., 2014), we observed an increased level of reads mapping to IGRs in LC-FA (Table 1), suggesting a relevant role for small RNAs and other IGRs when LC-FA are the culture media. We also observed striking differences when analyzing the functional categories of the genes of interest during the adaptation to hypoxia comparing both carbon sources. The abundance of genes belonging to the functional category of information pathways showed opposite trends in both conditions, being reduced in DNRP1 vs. DE but increased in FNRP1 vs. FE (Figure 2). This could indicate that the basic metabolic activity, important for the bacilli under hypoxia in LC-FA, was not that relevant when dextrose was the carbon source. A similar result was previously described when cholesterol was added to LC-FA media (Aguilar-Ayala et al., 2017), which therefore could represent a general condition associated to the Mtb adaptation to early hypoxia in lipid environments. According to these authors, the addition of cholesterol makes the bacteria more metabolically active in NRP1. In fact, the addition of cholesterol increased the percentage of reads mapped to CDS in NRP1 comparing to those found in the absence of cholesterol (Aguilar-Ayala et al., 2017) (Table 1). In relation to this, it has been found that a high-cholesterol diet was related to higher risk of active TB (Soh et al., 2016).

We observed that the number of over-expressed genes encoding regulatory proteins was different between hypoxic cultures supplemented with both carbon sources (Figure 2, Supplementary Table 3). The participation of those genes was relevant in early hypoxia in dextrose (Figure 3) and less relevant in LC-FA (Figure 4). Among sigma factors, sigE is one of the most studied of Mtb and it has been proposed as a central regulator of the stress response of the pathogen (Datta et al., 2011). Recently it was demonstrated that SigE has a major role in determining the level of basal tolerance of Mtb to antitubercular drugs (Pisu et al., 2017). We detected that sigE and sigB increased upon entrance to hypoxia regardless the nature of the carbon source, reinforcing the role described for sigE as a central regulator of hypoxia in Mtb (Du et al., 2016) a condition associated to drug tolerance. Our results indicate that sigB could be an additional component of this central system (Table 2). Supporting this notion is the fact that both sigE and sigB were found over-expressed during persistence (Du et al., 2016; Flentie et al., 2016) another stress condition for the bacilli.

We found that the central hypoxia regulator in dextrose, namely Rv0081 (Galagan et al., 2013) although over-expressed in NRP1 in both carbon sources, play a main role in DNRP1 (Figure 3) compared to FNRP1 (Figure 4) where other two transcriptional factors, namely TcrX and ClgR, seemed to have a main participation (Supplementary Figure 1 and Table 3).

Aiming to define a hypoxic signature (Table 4), we identified five transcriptional regulators with significantly high expression in DNRP1. Four of those have also significant high expression in FNRP2 and the remaining one, namely clgR (Rv2745), is also highly expressed in FNRP1. According to our data, this regulator may play a main role in the early adaptation to hypoxia in the LC-FA environment (Figure 4 and Supplementary Figure 1). ClgR is a transcriptional regulatory protein that has been involved in the maintenance of the membrane integrity of Mtb during response to stress (Veatch and Kaushal, 2018). In fact, clgR is induced in Mtb upon several stressful conditions including redox stress by diamide, SDS, hypoxia and low pH (McGillivray et al., 2014). In addition, a main role in the reactivation of bacilli after dormancy-induced by hypoxia was also assigned to this regulator (McGillivray et al., 2015). Our data in DNRP1 (Supplementary Table 3) were in accordance with the description of ClgR as negative regulator of Clp proteases in dextrose (McGillivray et al., 2015). This result could be related with the increased level of toxin-antitoxin activity found in dextrose (Table 3) through inactivation of antitoxins by Clps. On the contrary, the expression of Clp proteases did not change upon entrance to FNRP1 (Supplementary Table 3) suggesting a different role for ClgR in LC-FA. We found that the tubercle bacillus is in a reductive stress under LC-FA lipid environment (Rodríguez et al., 2014); therefore, the increased expression found for this regulator in FNRP1 (Figure 4 and Supplementary Table 3) could indicate that its role in lipid media was related with their participation in redox stresses. In agreement with this is the previous finding that the induction of clgR during redox stress did not result in the induction of clp genes (McGillivray et al., 2014). Our results highlight, for the first time, the influence of the environment in the complex regulatory activity of ClgR (Veatch and Kaushal, 2018).

One of the common regulatory factors between DNRP1 and FNRP2 was iciA (Rv1985c) (Figure 3 and Table 4). IciA is an in vitro replication initiation inhibitory protein (Kumar et al., 2009; Zhou et al., 2010; Marcus et al., 2016) probably involved in the entrance into periods of quiescence through the control of the chromosomal replication. Induction of iciA has been previously described during nutrient starvation (Zhou et al., 2010). Interestingly, iciA appears to be required at early hypoxia in dextrose, but it is not required until late hypoxia in the presence of LC-FA, what suggests, again, that that lipid environment is a less stressful condition for the tubercle bacilli. Notably, a single gene (Rv2554c) was only over-expressed in NRP2 regardless the carbon source, besides showing significant high expression in DNRP2 compared to FNRP2. This gene encodes a holliday junction resolvase-like protein, suggesting the implication of the reparation of double-strand breaks process in late hypoxia.

Other interesting regulators involved in the adaptation to hypoxia in LC-FA are sRNAs. These regulatory molecules were more relevant in the adaptation when LC-FA was the carbon source (Figure 5 and Supplementary Table 5). It is worth to mention that several ncRNAs with high expression in the adaptive pathway to hypoxia had been previously associated to slow down growth by Mtb (Arnvig and Young, 2012). This was the case of MTS1338, MTS0997, MTS0194, and MTS2823 (Figure 5) showing this activity as an apparently common feature of ncRNAs in this bacillus. The over-expression of either MTS1338 or MTS0194 also causes slow-down of the Mtb growth rate when the bacteria is cultured in the low-potassium in vitro dormancy model (Haning et al., 2014; Ignatov et al., 2014). On the other hand, the over-expression of MTS2823 causes down regulation of many genes involved in the energy metabolism (Arnvig and Young, 2012).

The MTS0194 is induced upon hydrogen peroxide and acid stress, two conditions that Mtb faces inside macrophages during active infection (Arnvig and Young, 2012). Interestingly, the cooper-inducible regulatory protein CsoR binds MTS0194 (Minch et al., 2015). CsoR loses its repressor activity under high levels of Cu (Marcus et al., 2016). The high level of expression found for MTS0194 in FNRP1 and FNRP2 suggests the loss of that repressor activity under hypoxia in LC-FA (see Supplementary Table 3) with increasing level of expression of MTS0194 as a consequence (see Supplementary Table 5 and Figure 1). This result could indicate that hypoxia in LC-FA represents a condition closer to that described inside the mycobacterial phagosome, where a high level of Cu is present (Marcus et al., 2016).

Regardless the culture media, some stable RNAs were detected with high expression in early and late hypoxia. The 4.5S RNA showed very high level of expression both in early and late hypoxia, being higher when Mtb is in dextrose compared to LC-FA (Supplementary Table 5). On the other hand, the 10S RNA also showed high level of expression in late hypoxia. These two RNAs are involved in key metabolic activities that help to the adequate function of the protein synthesis by the cell. The 4.5S RNA has been identified in Mtb as part of the signal recognition particle (SRP) together to the Ffh protein, and it is involved in the recognition of signal peptides emerging from the ribosomes, assisting the transport from the cell toward their functional compartments (Arnvig and Young, 2009; Palaniyandi et al., 2012). The 10S RNA, encoded by the gene ssrA, has been identified as stable tmRNA participating in the trans-translation in Mtb (Personne and Parish, 2014). The tmRNAs are responsible for recycling stalled ribosomes and thus ensuring their availability for protein synthesis, avoiding the accumulation of abnormal proteins by the bacteria. Different from other bacteria that have several systems to recover stalled ribosomes, Mtb has only the system represented by tmRNA and the accessory protein SmpB (Rv3100c). Contrary to smpB, the ssrA gene is essential in Mtb (Personne and Parish, 2014). Our data suggest for the first time the key role that both small RNAs may play in the survival of the tubercle bacillus during hypoxia. Our data also confirm the main contribution of small stable RNAs in the adaptation of Mtb to hypoxia, being that role particularly relevant in the dormancy-like lipid environment represented by LC-FA.

TA systems are considered relevant for persisting bacteria, helping their survival under stress conditions through the selected downregulation of targeted genes (Gerdes and Maisonneuve, 2012). In addition, the expression of these TA systems is regulated by environmental cues, which are relevant in the context of persistent infections, such as those caused by Salmonella enterica, Helicobacter pylori, and Mtb (De la Cruz et al., 2013; Cárdenas-Mondragón et al., 2016; Slayden et al., 2018). Our results suggest that TA modules play a role in the adaptation of Mtb to hypoxia in dextrose (Table 3) similarly to the survival mechanism switched on by the bacteria in a persistence-related stress condition (Butt and Titball, 2016). On the contrary, TA modules do not have a relevant participation in hypoxia in LC-FA (Table 3). It can be inferred that hypoxia does not represent a relevant stress for the tubercle bacilli when Mtb uses LC-FA as a carbon source. Interestingly, Aguilar-Ayala and co-workers (Aguilar-Ayala et al., 2017) detected participation of complete TA systems in NRP1 condition when cholesterol was added to the LC-FA medium. This result again suggested that cholesterol, a metabolic source of propionate, makes this last condition more toxic to the bacilli (Galagan, 2014) and could represent a stress environment that requires a decreased cell activity mediated by TA modules for survival.

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.