Differentially expressed genes varied dramatically among the WT and complemented strains
The sample M0129C_T20 was induced with tetracycline while M0129C_T0 was not. M. smegmatis mc2155 was applied as a control (Wt). To identify differentially expressed genes (DEGs), pair-wise comparisons of read counts between the complemented strain and wildtype were performed using DESeq (Supplementary Table 2). Although mtsp17 was down regulated in the 3 groups (M0129C_T0 versus T20; Wt versus M0129C_T20; Wt versus M0129C_T0), DEGs varied dramatically in quantity. Overall 4% protein-coding genes were DEGs (fold change > 2 and adjusted p-value < 0.05) in the M0129C_T0 versus T20 group while 24–25% in other 2 groups (Table 1). 50 DEGs (fold change > 2) overlapped between the three groups (Fig. 1A, Supplementary Table 3), among which 24 genes (including mtsp17) showed the same trends of fold change and 7 of these 24 genes have homolog in MTB (Fig. 1B, Table 2). The same change trends of DEGs among the 3 groups indicate that Mtsp17 has an important role of gene expression regulation. Furthermore, 15 genes were up regulated while 8 of the 23 DEGs were down regulated, implying that Mtsp17 could act as both an activator and a repressor.
Difference between the expression patterns of complemented strains was more moderate
Principal component analysis of the global expression patterns of the 3 samples demonstrated that wildtype clustered apart from the two complemented strain samples, and M0129C_T0 clustered apart from M0129C_T20 with the majority of variance occurring along PC2 (Fig. 1C). In addition, hierarchical clustering with the DEGs (adjusted p-value < 0.05) showed that gene expression pattern of the M0129C_T20 sample was more similar to M0129C_T0 (Fig. 1D). These results indicate that the difference between the expression patterns of M0129C_T20 and M0129C_T0 is more moderate than between the complemented strain and wildtype. Therefore focusing on the DEGs of the M0129C_T0 versus T20 group is more apt to unveil the mechanism that regulates the gene expression profile by Mtsp17 depression.
Table 1
Numbers of genes differentially expressed
Group | | M0129C_T0 /T20 | Wt/ M0129C_T20 | Wt/ M0129C_T0 |
mtsp17 ratio (log2 FC) | | -2.4 | -5.4 | -3 |
Differentially expressed genes (Padj < 0.05, log2 FC) | Down > 1 | 3% (207) | 11% (797) | 12%(861) |
Up > 1 | 1% (68) | 13%(896) | 13%(916) |
FC < 1 | 20%(1370) | 17% (1150) | 18%(1273) |
Others | | 76% | 59% | 57% |
Data were from RNA-seq with 3 replicates of each sample. FC, fold change. Padj, adjusted p-value. |
Table 2
Partial DEGs overlapped between the pair-wise comparisons
MSM | MTB | Essential[4] | Description | Functional category |
MSMEG_0129 | Rv0164 | Yes | Mtsp17 | Unknown |
MSMEG_3763 | Rv1686c | No | ABC transporter | Cell wall and cell processes |
MSMEG_0787 | Rv0411c | No | extracellular solute-binding protein | Cell wall and cell processes |
MSMEG_0788 | Rv0412c | No | hypothetical protein | Cell wall and cell processes |
MSMEG_2850 | Rv3178 | No | cell entry related family protein | Unknown |
MSMEG_1769 | Rv3288c | No | UsfY protein | Unknown |
MSMEG_4757 | Rv2524c | Yes | fatty acid synthase | Lipid metabolism |
Qpcr Analyses Verified The Rna Sequencing Data
Down regulation of Mtsp17 resulted in induction of 808 genes and repression of 837 genes (adjusted p-value < 0.05), representing 1645 transcriptional targets that related to Mtsp17 (Fig. 2A). 13 DEGs were randomly selected for qPCR analyses to verify the gene expression changes observed via RNA sequencing. The trend in fold changes of qPCR data matched RNA-seq results (Fig. 2B), thereby validating the down-regulated Mtsp17 transcriptome profile generated by RNA sequencing.
Mtsp17 acted as a transcriptional activator of desA1 gene
Gene ontology (GO) terms of the DEGs were assessed by GOseq v1.22 and 7 GO terms were significantly enriched (adjusted p-value < 0.05), including “structural constituent of ribosome” (GO:0003735) and “ribonucleoprotein complex” (GO:0030529) (Fig. 2C). In addition, quantity information of 58 ribosomal proteins was obtained in the RNA-seq data and 27 were DEGs (adjusted p-value < 0.05). But fold changes of all the 27 differentially expressed ribosomal proteins (twelve 30S ribosomal proteins and fifteen 50S ribosomal proteins) were less than 2.
Four of the 275 significantly differentially expressed genes (adjusted p-value < 0.05 and fold change > 2, Supplementary Table 4) were essential for mycobacterial growth, but only desA1 (MSMEG_5773) was down-regulated by depression of Mtsp17. DesA1 is a desaturase involved in mycolic acid biosynthesis in mycobacteria [10]. Changes of mycolic acids are central to adaptation of MTB to various environments [11, 12]. Transcription of MTB desA1 is regulated by the transcription factor MadR while overexpression of other transcription factors leads no change of desA1 mRNA [13]. The MSM MadR (encoded by MSMEG_0916) and DesA1 (encoded by MSMEG_5773) proteins show high similarity to the homologs in MTB (85% and 80% identity), respectively. The evolutionarily conserved MadR-mediated desA1 repression leads to decrease of mycolic acid biosynthesis and ensues loss of mycobacterial viability [12]. Analogously, the decline of bacterial growth in insufficient Mstp17 strains is likely due to depression of desA1 by down regulating mtsp17.
More transcriptome data related to Mtsp17-desA1 regulation were searched and summarized in Table 3. MadR overexpression made no significant change of mtsp17 mRNA in both MTB and MSM. The same phenotype was also found in MTB when 206 transcription factors were overexpressed respectively [13]. These phenomena indicate that Mtsp17 is probably an upstream factor of transcription factors which supports our hypothesis that Mtsp17 may regulate transcription through a DNA contact-independent mechanism. Under hypoxia, transcriptional change patterns of madR and desA1 were not in accordance with MadR repressing desA1 and could be explained by Mtsp17 activating desA1 (Table 3). The above results indicate that Mtsp17 plays as a transcriptional activator of desA1 gene.
Table 3
Log2 fold change value of genes in different pair-wise comparisons
Genes | MadR overexpression | MadR overexpression | Mtsp17 downregulation | Hypoxia |
madR | 1.48 | 1.82 | 0.13 | -0.16 |
desA1 | -1.46 | -1.75 | -1.62 | -3.25 |
mtsp17 | -0.18 | 0.06 | -2.41 | -1.53 |
Species | MSM | MTB | MSM | MTB |
GEO* | GSE116027 | GSE59086 | This study | GSE116353 |
*the Gene Expression Omnibus (GEO) series accession number. |
Characterization of SigF regulon in the mtsp17 complemented strain
Transcription in MSM is carried out by a multi-subunit RNA polymerase and one of the 28 sigma subunits [14]. The sigma subunit ensures the transcription machinery to initiate transcription of particular genes. Sigma factor is post-translationally regulated by its cognate anti-sigma factors, which sequester them from the transcription machinery and could be neutralized by anti-sigma factor antagonists [15].
There was one anti-sigma factor or its antagonist, which is an anti-sigma factor antagonist of SigF (MSMEG_0586, 8.9-fold decrease, adjusted p-value = 7.3E-40) among the top 20 differentially expressed genes in the M0129C_T0 versus M0129C_T20 group. SigF controls the expression of a particular subset of genes (SigF regulon) by altering RNA polymerase specificity [16]. The DEGs in the M0129C_T0 versus M0129C_T20 group (M0129C_T0/T20) were compared with fold changes in a SigF mutant [17]. 70 of seventy nine DEGs with log2 fold change more than 1.5 and adjusted p-value less than 0.05 showed the same trend of fold changes (Supplementary Table 5). Furthermore, 10 of the seventy nine DEGs are involved in transcriptional regulatory mechanism and they all showed the same trend in SigF knockout versus wildtype group and M0129C_T0 versus M0129C_T20 group except MSMEG_0529 (Fig. 3A). The above results indicate that the majority of significantly differentially expressed genes regulated by Mtsp17 belong to SigF regulon.
MSMEG_0529 encodes the Serine/threonine-protein kinase PknK and participates in regulatory pathways that slow growth of M. tuberculosis under stresses [18]. SigF overexpression produced moderate upregulation of pknK (log2 fold change = 0.91) [13]. Here, down regulation of mtsp17 led to decrease of the anti-SigF antagonist gene (8.9-fold) and upregulation of pknK (2.9-fold). These results indicate that Mtsp17 rather than SigF is the key factor of the pknK regulatory network. More details are needed to be explored to full fill the regulatory relationships between pknK gene and Mtsp17.