Single START-domain protein Mtsp17 is involved in transcriptional regulation in Mycobacterium smegmatis

About the Authors:

Ying Zhou

Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Software, Visualization, Writing – original draft, Writing – review & editing

Affiliation: Department of Bone and Joint Surgery, The First Affiliated Hospital of Jinan University, Guangzhou, China

ORCID logo https://orcid.org/0000-0002-6347-2502

Tianying Zhong

Roles Writing – review & editing

Affiliation: Guangdong Province Green and High Performance Novel Materials Engineering Research Center, Jiangmen Polytechnic, Jiangmen, China

Wenjing Wei

Roles Funding acquisition, Resources, Writing – original draft

Affiliation: Center for Tuberculosis Control of Guangdong Province, Guangzhou, China

Zhuhua Wu

Roles Data curation, Methodology, Resources

Affiliation: Center for Tuberculosis Control of Guangdong Province, Guangzhou, China

Anping Yang

Roles Data curation, Investigation, Writing – review & editing

Affiliation: School of Medicine, Foshan University, Foshan, Guangdong, China

Ning Liu

Roles Funding acquisition, Writing – original draft, Writing – review & editing

* E-mail: [email protected] (NL); [email protected] (MW); [email protected] (XZ)

Affiliation: Department of Bone and Joint Surgery, The First Affiliated Hospital of Jinan University, Guangzhou, China

Ming Wang

Roles Methodology, Software, Validation, Writing – original draft, Writing – review & editing

* E-mail: [email protected] (NL); [email protected] (MW); [email protected] (XZ)

Affiliation: Key Laboratory of RNA Biology, Institute of Biophysics, Chinese Academy of Sciences, Beijing, China

Xiaoli Zhang

Roles Conceptualization, Funding acquisition, Methodology, Resources, Writing – original draft, Writing – review & editing

* E-mail: [email protected] (NL); [email protected] (MW); [email protected] (XZ)

Affiliation: School of Medicine, Foshan University, Foshan, Guangdong, China

Introduction

Mtsp17 (Rv0164) is a 17-kDa protein isolated from Mycobacterium tuberculosis (MTB) culture filtrates that has been shown to be a T-cell and B-cell antigen [1,2]. It is highly conserved across the Mycobacteria. MTB Mtsp17 shares less than 30% homology with its homologues in other genera, making it not only an interesting anti-TB drug target but also an immunogen that may be worth investigating in vaccine development. Although mtsp17 has been shown to be an essential gene in both Mycobacterium smegmatis (MSM) and MTB [3,4], its physiological function is presently unclear. Considering its essentiality in the genus Mycobacterium, further investigation of Mtsp17 will likely provide new insights into the physiology of Mycobacteria.

Mtsp17 folds in a similar manner to polyketide cyclases/aromatases (CYC/ARO), even though they share less than 20% sequence identity [3]. CYC/ARO proteins are monofunctional enzymes in Type II polyketide synthase (PKS) complexes and transform linear poly-β-ketone intermediates into aromatic polyketides [5]. Twenty-one Type I and three Type III PKSs have been annotated in Mycobacteria, but Type II PKS proteins have not yet been identified [6]. Structural comparisons and docking studies have suggested that Mtsp17 is likely not an enzyme that catalyzes polyketide cyclization and aromatization [3]. The essential functions of Mtsp17 in mycobacteria are therefore likely to operate by mechanisms that are different from CYC/ARO proteins.

Mtsp17 is a single-domain steroidogenic acute regulatory protein-related lipid transfer (START) domain protein. START domains are lipid-binding domains that function as polyketide cyclases/aromatases in lipid metabolism. START domains also bind ligands during non-vesicular traffic between intracellular compartments and fuse with DNA binding domains to modulate gene transcription [7–9]. As mycobacteria lack internal organelles such as the endoplasmic reticulum and Golgi apparatus, the START domain in Mtsp17 may regulate transcription through a DNA contact-independent mechanism. Comprehensive profiling of the genes and pathways regulated by Mtsp17 should provide insights into its physiological functions.

In this study, we performed RNA-sequencing (RNA-seq) based transcriptomics analysis to obtain a global profile of genes regulated by Mtsp17 in MSM. We compared the transcriptomes of wildtype and the mtsp17-complemented strains (mtsp17 induced and not induced), and discovered that Mtsp17 has transcriptional regulatory properties. This transcriptional regulatory role was further validated by qPCR and its implications are discussed.

Materials and methods

Bacterial strains and culture conditions

As knockout of mtsp17 was lethal in MSM, a conditional Mtsp17 knockout strain (M0129C) was created by a specialized transduction procedure described previously [3] in which plasmid-encoded Mtsp17 was induced by tetracycline while the genomic mtsp17 gene was deleted (M0129C was also the mtsp17-complemented strain). Over-expression of Mtsp17 was accomplished with a pMV261 plasmid that constitutively over-expressed Mtsp17 using a Hsp60 promoter [10]. Mycobacterium smegmatis strains were cultured in Middlebrook 7H9 medium (Difco, Baltimore, MD, USA) supplemented with 0.5% glycerol and 0.05% Tween 80 at 37°C. 25 μg.mL-1 kanamycin and 20 ng.mL-1 tetracycline (Sangon Biotech, Shanghai, China) were required for culturing the mtsp17-complemented strain in which the genomic mtsp17 gene was replaced by a sacB-hyg cassette and an extrachromosomal mtsp17 gene was introduced using a pMind-derived plasmid. After cultures of the complemented strain reached exponential phase in the presence of tetracycline (OD600 = 0.4), they were transferred into fresh culture medium and divided into two halves. Tetracycline (20 ng.mL-1) was added to one half of the culture (M0129C_T20, induced) and incubated for 4 hours, while the other half was incubated without tetracycline (M0129C_T0, not induced) for 4 more hours. Total RNA was then isolated from these two cultures of the complemented strain for sequencing. An exponential phase culture of Mycobacterium smegmatis mc2155 was used as a control (Wt).

RNA isolation and transcriptome analysis

10 ml of exponentially growing cultures was harvested by centrifugation (4°C, 3000 rpm, 10 min). Total RNA was isolated from three biological replicates of each sample using a FastPrep Instrument and FastRNA Pro Blue Kits (MP Biomedicals, Solon, OH, USA) according to the manufacturer’s instructions. After treatment with DNase I, the RNA was precipitated with ammonium acetate/isopropanol. The quality and quantity of each RNA sample were assessed using an Agilent 2100 Bioanalyzer (Agilent Technologies, Beijing, China) and a NanoDrop 1000 spectrophotometer (Thermofisher, Shanghai, China). A Ribo-Zero rRNA Removal Kit (Epicentre, Madison, WI, USA) was used to remove bacterial rRNA from total RNA preparations.

cDNA libraries were constructed using a NEBNextUltraTM RNA library Prep Kit (New England Biolabs, Beijing, China) and sequenced on a HiSeq 4000 sequencer (Illumina). The processed reads were mapped to the Mycobacterium smegmatis MC2 155 genome (NC_008596) using Bowtie2 v2.2.6 with default parameters. FPKM values for each gene were compared and differentially expressed genes (DEGs) were determined with HTSeq v0.6.0 and DESeq v1.22.1. Genes with a multiple hypothesis-adjusted p-value below 0.05 (adjusted p-value < 0.05) were considered as DEGs. Functional enrichment analysis was carried out using GOseq v1.22 based on the Gene Ontology database.

Real-time fluorescence quantification PCR (qPCR) analysis

qPCR was used to measure the expression of selected genes in a CFX96 instrument (Bio-Rad). The GoScript Reverse Transcription System (Promega, Beijing, China) was used for cDNA preparation, and a 2x SYBR Green mix (Vazyme, Nanjing, China) and gene-specific primers (S1 Table) were used to amplify the genes (PCR cycle: 95°C, 10 s; 60°C, 30 s; 40 cycles). Primers were optimized to give 90–110% amplification efficiency and had a single melting temperature. sigA (MSMEG_2758) was used as an internal normalization standard. Relative quantification was performed using the ΔΔCT method. Fold changes presented are means ± standard deviations of three technical replicates.

Results

RNA-sequencing results uncover the importance of Mtsp17 in transcription regulation

To investigate if Mtsp17 functions in transcription regulation, we looked for DEGs in the transcriptomes of the conditional Mtsp17 knockout strain (M0129C) and the wild-type strain (Wt) (S2 Table). Although mtsp17 was down regulated in the three comparisons (M0129C_T0 versus T20; Wt versus M0129C_T20; Wt versus M0129C_T0), the number of DEGs varied dramatically in each comparison. Overall, 4% of protein-coding genes were differentially-expressed (fold change > 2 and adjusted p-value < 0.05) in the M0129C_T0 versus T20 comparison compared to 24–25% in other 2 comparisons (Table 1). 50 DEGs (fold change > 2) were common to the three comparisons (Fig 1A, S3 Table), and of these, 24 genes (including mtsp17) showed the same trends in fold change. We noted that 7 of the 24 genes had homologues in MTB (Fig 1B, Table 2). Similar trends in the expression of DEGs in the 3 comparisons suggest that Mtsp17 plays an important role in regulating gene expression. Fifteen of the 24 genes were up-regulated and 8 were down-regulated, implying that Mtsp17 can act as both an activator and a repressor.

[Figure omitted. See PDF.]

Fig 1. Differentially expressed genes in pair-wise comparisons.

(A) Venn diagram showing the numbers of significantly differentially expressed genes (fold change > 2, adjusted p-value < 0.05). (B) Differentially-expressed genes that show similar fold change trends. (C) Principal component analysis of the global expression patterns of the 3 samples. (D) Hierarchical clustering of the 2146 differentially expressed genes (adjusted p-value < 0.05).

https://doi.org/10.1371/journal.pone.0249379.g001

[Figure omitted. See PDF.]

Table 1. Numbers of differentially expressed genes.

https://doi.org/10.1371/journal.pone.0249379.t001

[Figure omitted. See PDF.]

Table 2. DEGs common to all three of the pair-wise comparisons.

https://doi.org/10.1371/journal.pone.0249379.t002

Expression patterns in the two complemented strains were moderately different

Principal component analysis of the global expression patterns of the 3 samples clustered Wt separately from the two complemented strain samples. M0129C_T0 clustered apart from M0129C_T20, with the majority of the variance occurring along PC2 (Fig 1C). Hierarchical clustering of the DEGs (adjusted p-value < 0.05) showed that the gene expression pattern of the M0129C_T20 sample was similar to M0129C_T0 and was different from Wt (Fig 1D). These results indicate that the difference in expression patterns between M0129C_T20 and M0129C_T0 was more moderate than that between the complemented strain and the Wt strain. Comparing the DEGs from the M0129C_T0 and T20 samples is thus more likely to reveal the mechanism by which Mtsp17 regulates gene expression.

qPCR analysis is consistent with transcriptome data

Down regulation of Mtsp17 affected the expression of 1645 genes, inducing the expression of 808 genes and repressing that of 837 genes (adjusted p-value < 0.05) (Fig 2A). We randomly selected 13 DEGs for qPCR analyses to verify the gene expression changes observed in our transcriptome data. Trends in fold changes of the qPCR data matched our RNA-seq results (Fig 2B), validating the down-regulated Mtsp17 transcriptome profile generated by RNA sequencing.

[Figure omitted. See PDF.]

Fig 2. Differentially expressed genes in the M0129C_T0 versus M0129C_T20 comparison.

(A) Volcano plot showing differentially expressed genes. Green, down-regulated genes. Red, up-regulated genes. Blue, genes not differentially expressed. Padj, adjusted p-value. (B) qPCR analyses of 13 randomly selected DEGs. (C) Significantly enriched gene ontology (GO) terms (adjusted p-value < 0.05).

https://doi.org/10.1371/journal.pone.0249379.g002

Mtsp17 acts as a transcriptional activator of the desA1 gene

Gene ontology (GO) analysis of DEGs using GOseq v1.22 indicted that 7 GO terms, including “structural constituent of ribosome” (GO:0003735) and “ribonucleoprotein complex” (GO:0030529), were significantly enriched (adjusted p-value < 0.05) (Fig 2C). Twenty-seven of the 58 ribosomal proteins for which quantitative data was obtained by RNA-seq were differentially-expressed (adjusted p-value < 0.05), however the fold changes of all 27 differentially expressed ribosomal proteins (12 30S ribosomal proteins and 15 50S ribosomal proteins) were less than 2.

275 genes were significantly differentially expressed (adjusted p-value < 0.05 and fold change > 2, S4 Table) and 4 of the 275 genes were essential for mycobacterial growth, but only desA1 (MSMEG_5773) was down-regulated (Log2 fold change = -1.62). DesA1 was also positively correlated with Mtsp17 in qPCR data from the conditional Mtsp17 knockout and over-expression strains (Fig 3).

[Figure omitted. See PDF.]

Fig 3. Fold changes of DesA1 and Mtsp17 mRNA.

Mtsp17 KD, a conditional Mtsp17 knockout strain, incubated in the absence of the inducer tetracycline for 4 additional hours. Mtsp17 OE, the Mtsp17 over-expression strain.

https://doi.org/10.1371/journal.pone.0249379.g003

DesA1 is a desaturase involved in mycolic acid biosynthesis in mycobacteria [11]. Changes in mycolic acids are central to the adaptation of MTB to various environments [12,13]. Transcription of desA1 is known to be regulated by the transcription factor MadR, and is reported to be not affected by other transcription factors [13,14]. MadR-mediated desA1 repression is evolutionarily conserved and leads to a decrease in mycolic acid biosynthesis and ensuing loss of mycobacterial viability [12]. By analogy, the decline in bacterial growth in Mtsp17 insufficient strains may be due to a decrease in desA1 transcription.

More transcriptome data related to Mtsp17-desA1 regulation were searched in the Gene Expression Omnibus (GEO) database and are summarized in Table 3. Overexpression of MadR did not significantly alter levels of mtsp17 mRNA in either MTB or MSM [13,14]. When 206 transcription factors were overexpressed, mtsp17 mRNA also showed little change [14]. These phenomena suggest that Mtsp17 may be an upstream factor of transcription factors, and support 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 is likely a transcriptional activator of the desA1 gene.

[Figure omitted. See PDF.]

Table 3. Log2 fold change values of genes in different pair-wise comparisons.

https://doi.org/10.1371/journal.pone.0249379.t003

Mtsp17 is involved in regulation of the SigF regulon via an anti-SigF antagonist

Transcription in MSM is carried out by a multi-subunit RNA polymerase and one of 28 sigma subunits [15]. The sigma subunit ensures the transcription machinery to initiate transcription of particular genes [16,17]. Sigma factors are post-translationally regulated by a partner switching system (PSS) in which anti-sigma factors sequester the sigma factor from the transcription machinery and anti-sigma factor antagonists lift the repression of sigma factors [18–20].

MSMEG_0586, an anti-SigF antagonist, was one of the top 20 differentially expressed genes, showing an 8.9-fold decrease (log2 fold change = -3.16, adjusted p-value = 7.3E-40) in the M0129C_T0 versus M0129C_T20 comparison, while other sigma factors or anti-sigma factors or anti-anti-sigma factors showed log2 fold changes of no more than 1.5. MSMEG_0586 is involved in the PSS of the SigF regulatory pathway [20]. In light of significant down-regulation of the anti-SigF antagonist MSMEG_0586 and no obvious change in other SigF PSS factors in the conditional Mtsp17 knockout strain, we hypothesized that transcriptomic changes partially overlapped between the conditional Mtsp17 knockout and the SigF knockout strains. When DEGs in the M0129C_T0 versus M0129C_T20 comparison were compared with DEGs in the SigF knockout versus wildtype comparison [21], 70 of 79 DEGs with log2 fold changes > 1.5 and adjusted p-values of < 0.05 showed the same fold change trends (S5 Table). Furthermore, 10 of the 79 DEGs were involved in transcriptional regulatory mechanisms and they all had the same varying tendency in the SigF knockout versus wildtype comparison and M0129C_T0 versus M0129C_T20 comparison, except for MSMEG_0529 (Fig 4A). The above results show that the majority of significantly differentially expressed genes regulated by Mtsp17 belong to the SigF regulon.

[Figure omitted. See PDF.]

Fig 4. Comparisons of gene expression in Mtsp17 and SigF mutants.

(A) 10 differentially expressed genes with log2 fold change > 1.5 that are involved in transcriptional regulatory mechanisms and showed the same varying tendency in the SigF knockout (SigF KO) and conditional Mtsp17 knockout (Mtsp17 KD) strains. (B) Fold changes of MSMEG_0586 (anti-anti-SigF) and MSMEG_5590 mRNA in the conditional Mtsp17 knockout (Mtsp17 KD), Mtsp17 over-expression (Mtsp17 OE) and SigF knockout (SigF KO) strains. Data of the SigF knockout strain are from GSE19774.

https://doi.org/10.1371/journal.pone.0249379.g004

To further understand the relationship between Mtsp17 and the SigF regulon, we measured levels of MSMEG_0586 and MSMEG_5590 mRNA in the Mtsp17 over-expression strain by qPCR. MSMEG_5590 was chosen as it belongs to the SigF regulon and was significantly correlated with MSMEG_0586 in the SigF knockout strain (Fig 4B SigF KO). Both MSMEG_0586 and MSMEG_5590 were down-regulated in the conditional Mtsp17 knockout mutant and up-regulated in the Mtsp17 over-expression strain (Fig 4B). The above results indicate that Mtsp17 regulates transcription of MSMEG_0586 and alterations in MSMEG_0586 mRNA levels lead to transcriptional changes in the SigF regulon.

Discussion

MTB genes that are essential are good targets for drug discovery [22]. Elucidating the mechanism of essential genes is of benefit to the development of more effective vaccines and better drugs. Mtsp17 is an essential protein in both MTB and MSM, but understanding of its function is limited. In the transcriptomic analysis presented here, we have identified genes regulated by Mtsp17 and demonstrated for the first time that a START family protein with a mono domain can regulate transcription in MSM.

The anti-SigF antagonist MSMEG_0586 was one of the top 20 differentially expressed genes in the conditional Mtsp17 knockout mutant and showed an 8.9-fold decrease. Since MSMEG_0586 is as an activator of SigF by post-translational regulation and 70 of 79 DEGs with log2 fold changes > 1.5 belong to the SigF regulon (S5 Table), transcriptional changes in these SigF regulon genes can be explained by alterations in MSMEG_0586 expression by Mtsp17 and lead to post-translational regulation of the SigF regulon by the PSS of the SigF regulatory pathway. As both SigF and Mtsp17 are conserved in the Mycobacteria and are transcriptionally regulated in response to nutrient depletion and oxidative stress [10,20], we propose that the regulation of the SigF regulon by Mtsp17 may operate in conjunction with the SigF regulatory pathway during stress responses.

Expression of PknK (MSMEG_0529) showed different trends in terms of fold changes in the Mtsp17 and SigF mutants (Fig 4A). As pknK mRNA increases 1.9-fold on SigF over-expression [14] and SigF knockout has no effect on pknK mRNA levels [21], we propose the presence of SigF dependent (when SigF is abundant) and SigF independent (when SigF is inadequate) PknK regulatory pathways. Here, down-regulation of mtsp17 led to a decrease in the anti-SigF antagonist MSMEG_0586 (8.9-fold) and upregulation of pknK (2.9-fold). These results correspond to the presence of SigF independent regulation of PknK when SigF is functionally inadequate due to down regulation of the anti-SigF antagonist.

The START protein PCTP interacts with Pax3, a mammalian transcription factor, to regulate gene transcription [23]. In the case of MSM, START protein Mtsp17 interacts with the global transcription factor CarD [10], and the percentage of down-regulated proteins detected here for the Mtsp17 mutant was similar to that reported for the CarDR47E mutant (3% versus 3.1% [24]). However, it is reported that there are no obvious transcriptional changes in DesA1 in CarD mutants, and the anti-SigF antagonist increases in both the loss-of-function (with weakened affinity for RNAP or DNA) and gain-of-function (with increased affinity to RNAP) mutants of CarD [24]. The discrepancy in the transcriptional changes of DesA1 and the anti-SigF antagonist between the Mtsp17 and CarD mutants suggests it is unclear whether CarD plays a key role in this Mtsp17-induced regulatory pathway. Furthermore, we cannot rule out the possibility that other DNA binding proteins are involved in this pathway. Further investigation of Mtsp17 interacting proteins will be necessary to provide a complete understanding of the transcriptional regulatory network of Mtsp17.

This work shows the landscape of Mtsp17 transcriptional regulation and has uncovered important roles of Mtsp17 in the regulatory network of the DesA1 and SigF regulon. As Mtsp17, DesA1 and the SigF regulon are all associated with stress [10,13,25], transcriptional regulation by Mtsp17 likely acts together with other regulatory mechanisms to help coordinate signal transfer under stresses such as nutritional deficiency. Elucidating the detailed mechanism of the Mtsp17 transcriptional regulatory network identified here will require further study.

Supporting information

S1 Table. Primers used in this study.

https://doi.org/10.1371/journal.pone.0249379.s001

(XLSX)

S2 Table. Gene expression profiles of M. smegmatis mc2155 strain (Wt) and mtsp17-complemented strains.

https://doi.org/10.1371/journal.pone.0249379.s002

(XLSX)

S3 Table. 50 DEGs (fold change > 2, adjusted p-value < 0.05) overlapping between the three comparisons.

https://doi.org/10.1371/journal.pone.0249379.s003

(XLSX)

S4 Table. 275 significantly differentially expressed genes (fold change > 2, adjusted p-value < 0.05) in the M0129C_T0 versus T20 comparison.

https://doi.org/10.1371/journal.pone.0249379.s004

(XLSX)

S5 Table. 70 DEGs (log2 fold change > 1.5, adjusted p-value < 0.05) showing the same change tendency in the M0129C_T0 versus M0129C_T20 comparison and SigF knockout versus wildtype comparison.

https://doi.org/10.1371/journal.pone.0249379.s005

(XLSX)

Acknowledgments

The authors would like to express their thanks to Beijing Allwegene BioTech for technical assistance with RNA-seq, and to Dr Joy Fleming (Institute of Biophysics, Chinese Academy of Sciences, Beijing) for help with manuscript preparation.

Citation: Zhou Y, Zhong T, Wei W, Wu Z, Yang A, Liu N, et al. (2021) Single START-domain protein Mtsp17 is involved in transcriptional regulation in Mycobacterium smegmatis. PLoS ONE 16(4): e0249379. https://doi.org/10.1371/journal.pone.0249379

References

1. Lim JH, Kim HJ, Lee KS, Jo EK, Song CH, Jung SB, et al. Identification of the new T-cell-stimulating antigens from Mycobacterium tuberculosis culture filtrate. FEMS microbiology letters. 2004;232(1):51–9. pmid:15019734.

2. Ireton GC, Greenwald R, Liang H, Esfandiari J, Lyashchenko KP, Reed SG. Identification of Mycobacterium tuberculosis antigens of high serodiagnostic value. Clinical and vaccine immunology: CVI. 2010;17(10):1539–47. pmid:20719988; PubMed Central PMCID: PMC2952993.

3. Zheng S, Zhou Y, Fleming J, Zhou Y, Zhang M, Li S, et al. Structural and genetic analysis of START superfamily protein MSMEG_0129 from Mycobacterium smegmatis. FEBS letters. 2018;592(8):1445–57. pmid:29512898.

4. DeJesus MA, Gerrick ER, Xu W, Park SW, Long JE, Boutte CC, et al. Comprehensive Essentiality Analysis of the Mycobacterium tuberculosis Genome via Saturating Transposon Mutagenesis. mBio. 2017;8(1). pmid:28096490; PubMed Central PMCID: PMC5241402.

5. Zhang Z, Pan HX, Tang GL. New insights into bacterial type II polyketide biosynthesis. F1000Research. 2017;6:172. pmid:28299197; PubMed Central PMCID: PMC5321127.

6. Cole ST, Brosch R, Parkhill J, Garnier T, Churcher C, Harris D, et al. Deciphering the biology of Mycobacterium tuberculosis from the complete genome sequence. Nature. 1998;393(6685):537–44. pmid:9634230.

7. Li S, Chen N, Li F, Mei F, Wang Z, Cheng X, et al. Characterization of wheat homeodomain-leucine zipper family genes and functional analysis of TaHDZ5-6A in drought tolerance in transgenic Arabidopsis. BMC plant biology. 2020;20(1):50. pmid:32005165; PubMed Central PMCID: PMC6993422.

8. Schrick K, Bruno M, Khosla A, Cox PN, Marlatt SA, Roque RA, et al. Shared functions of plant and mammalian StAR-related lipid transfer (START) domains in modulating transcription factor activity. BMC biology. 2014;12:70. pmid:25159688; PubMed Central PMCID: PMC4169639.

9. Horenkamp FA, Valverde DP, Nunnari J, Reinisch KM. Molecular basis for sterol transport by StART-like lipid transfer domains. The EMBO journal. 2018;37(6). pmid:29467216; PubMed Central PMCID: PMC5852651.

10. Zhou Y, Wei W, Fleming J, Ye C, Zheng S, Liu F, et al. Mycobacterium smegmatis MSMEG_0129 is a nutrition-associated regulator that interacts with CarD and ClpP2. The international journal of biochemistry & cell biology. 2020;124:105763. pmid:32389745.

11. Singh A, Varela C, Bhatt K, Veerapen N, Lee OY, Wu HH, et al. Identification of a Desaturase Involved in Mycolic Acid Biosynthesis in Mycobacterium smegmatis. PloS one. 2016;11(10):e0164253. pmid:27741286; PubMed Central PMCID: PMC5065219 data and materials.

12. Ehrt S, Schnappinger D, Rhee KY. Metabolic principles of persistence and pathogenicity in Mycobacterium tuberculosis. Nature reviews Microbiology. 2018;16(8):496–507. pmid:29691481; PubMed Central PMCID: PMC6045436.

13. Peterson EJ, Bailo R, Rothchild AC, Arrieta-Ortiz ML, Kaur A, Pan M, et al. Path-seq identifies an essential mycolate remodeling program for mycobacterial host adaptation. Molecular systems biology. 2019;15(3):e8584. pmid:30833303; PubMed Central PMCID: PMC6398593.

14. Rustad TR, Minch KJ, Ma S, Winkler JK, Hobbs S, Hickey M, et al. Mapping and manipulating the Mycobacterium tuberculosis transcriptome using a transcription factor overexpression-derived regulatory network. Genome biology. 2014;15(11):502. pmid:25380655; PubMed Central PMCID: PMC4249609.

15. Hurst-Hess K, Biswas R, Yang Y, Rudra P, Lasek-Nesselquist E, Ghosh P. Mycobacterial SigA and SigB Cotranscribe Essential Housekeeping Genes during Exponential Growth. mBio. 2019;10(3). pmid:31113892; PubMed Central PMCID: PMC6529629.

16. Chen J, Boyaci H, Campbell EA. Diverse and unified mechanisms of transcription initiation in bacteria. Nature reviews Microbiology. 2021;19(2):95–109. pmid:33122819.

17. Vishwakarma RK, Brodolin K. The sigma Subunit-Remodeling Factors: An Emerging Paradigms of Transcription Regulation. Frontiers in microbiology. 2020;11:1798. pmid:32849409; PubMed Central PMCID: PMC7403470.

18. Chauhan R, Ravi J, Datta P, Chen T, Schnappinger D, Bassler KE, et al. Reconstruction and topological characterization of the sigma factor regulatory network of Mycobacterium tuberculosis. Nature communications. 2016;7:11062. pmid:27029515; PubMed Central PMCID: PMC4821874.

19. Oh Y, Song SY, Kim HJ, Han G, Hwang J, Kang HY, et al. The Partner Switching System of the SigF Sigma Factor in Mycobacterium smegmatis and Induction of the SigF Regulon Under Respiration-Inhibitory Conditions. Frontiers in microbiology. 2020;11:588487. pmid:33304334; PubMed Central PMCID: PMC7693655.

20. Andrews ESV, Rzoska-Smith E, Arcus VL. Post-transcriptional modulation of the SigF regulon in Mycobacterium smegmatis by the PhoH2 toxin-antitoxin. PloS one. 2020;15(7):e0236551. pmid:32726339; PubMed Central PMCID: PMC7390352.

21. Singh AK, Dutta D, Singh V, Srivastava V, Biswas RK, Singh BN. Characterization of Mycobacterium smegmatis sigF mutant and its regulon: overexpression of SigF antagonist (MSMEG_1803) in M. smegmatis mimics sigF mutant phenotype, loss of pigmentation, and sensitivity to oxidative stress. MicrobiologyOpen. 2015;4(6):896–916. pmid:26434659; PubMed Central PMCID: PMC4694148.

22. Johnson EO, LaVerriere E, Office E, Stanley M, Meyer E, Kawate T, et al. Large-scale chemical-genetics yields new M. tuberculosis inhibitor classes. Nature. 2019;571(7763):72–8. pmid:31217586.

23. Kanno K, Wu MK, Agate DS, Fanelli BJ, Wagle N, Scapa EF, et al. Interacting proteins dictate function of the minimal START domain phosphatidylcholine transfer protein/StarD2. The Journal of biological chemistry. 2007;282(42):30728–36. pmid:17704541.

24. Zhu DX, Garner AL, Galburt EA, Stallings CL. CarD contributes to diverse gene expression outcomes throughout the genome of Mycobacterium tuberculosis. Proceedings of the National Academy of Sciences of the United States of America. 2019;116(27):13573–81. pmid:31217290; PubMed Central PMCID: PMC6613185.

25. Humpel A, Gebhard S, Cook GM, Berney M. The SigF regulon in Mycobacterium smegmatis reveals roles in adaptation to stationary phase, heat, and oxidative stress. Journal of bacteriology. 2010;192(10):2491–502. pmid:20233930; PubMed Central PMCID: PMC2863567.