INTRODUCTION
Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), is refractory to treatment with a single antibiotic and, despite combination therapy, widespread multidrug resistance is a growing concern (
1,
2). New therapeutic approaches are required to subvert both tolerance and outright drug resistance in mycobacteria. In addition to establishing cell integrity and facilitating nutrient import, the mycobacterial cell wall and its outer membrane serve as critical determinants of mycobacterial drug tolerance and susceptibility (
3–5). Given the relationship between the cell wall and intrinsic antibiotic susceptibility, a more complete understanding of the genetic networks that affect its construction will facilitate the development of new therapeutics.
The synthesis and remodeling of the mycobacterial cell wall lipid/peptidoglycan components are highly regulated (
6), yet little is understood about how the membrane-associated proteome is properly synthesized. In bacteria, recent advances in techniques such as quantitative proteomics, ribosome profiling, and cryo-electron microscopy (cryo-EM) have revealed that an additional layer of regulation during protein synthesis serves to maintain cellular homeostasis (
7,
8) and enable appropriate synthesis of membrane proteins (
9). With this increased resolution, it is clear that ribosome function can be altered by proteins and RNA factors that associate to them in a message- and environment-specific manner (
10,
11). These associating factors help to control translation rate, protein folding, and localization, ultimately contributing to a cell with a spatially and temporally regulated proteome (
12–14).
In mycobacteria, it is plausible that regulation at the ribosome may contribute to pathogenicity and drug susceptibility (
15). Studies of ribosome-associated proteins such as HflX, Mpy, and LepA have linked these factors to drug tolerance in both
M. tuberculosis and
M. smegmatis (
16–18). While some of these enzymes have clear roles in ribosome stability and hibernation, less is known about how LepA associates with the mycobacterial ribosome and affects drug tolerance (
19). Specifically, the link between LepA and antibiotic tolerance in mycobacteria was uncovered in a screen for genetic determinants of single-cell heterogeneity and drug susceptibility. Mutations in several genes, including
lepA, altered the rate of antibiotic-mediated killing (
18).
LepA is a ribosome-dependent GTPase found in almost all organisms, from bacteria to human mitochondria (
19,
20). It uses four classical elongation factor protein domains to contact the ribosome and hydrolyze GTP, occupying the same position on the 70S ribosome as elongation factor G (EF-G) (
21). It is nonessential in most bacteria, but it has been speculated that LepA may confer a fitness benefit in certain growth conditions, such as altered cation concentrations or low pH (
22–27). Despite its conservation, the physiological role of LepA remains unclear, as underlined by two different proposed roles for this GTPase. In
Escherichia coli, loss of LepA results in decreased polysome formation, suggesting a role in initiation and ribosome assembly (
27–29). Alternatively, structural studies on multiple bacterial LepA homologs indicate that its C-terminal domain makes contact with the A/P-site tRNA of the 70S ribosome and may alter the conformation of the ribosome-tRNA complex (
30,
31), perhaps participating in translational quality control (
29,
32,
33).
Here, we investigate the mechanistic basis of altered drug susceptibility of a mycobacterial lepA deletion mutant. We find that LepA augments protein levels for certain members of the mycobacterial porin (Msp) family during translation, as well as a number of other membrane-associated proteins. LepA deficiency results in decreased synthesis of MspA, the major porin in M. smegmatis, and a reduction in cell permeability as measured by dye accumulation and killing by certain antibiotics. Thus, we find that LepA acts as a translational aid in mycobacteria, providing evidence for its role in maintenance of prominent members of the mycobacterial membrane proteome, and demonstrating ribosome-based control of membrane homeostasis.
DISCUSSION
LepA is highly conserved across the kingdoms of life, yet its cellular role remains unclear. While
lepA mutants are viable, they have a very specific permeability defect. Considering our data and that of others suggesting that porins are important for drug accessibility (
18,
46), it is unsurprising that both
lepA and
mspA mutants are less susceptible to rifampin and vancomycin. This change in uptake could be due to either indirect changes to outer membrane permeability resulting from loss of porin, or to direct decreases in drug transport through the porin itself. In fact, changes in MspA function alter the transport of a large number of nutrients across the cell wall (
41,
47–49). Our data also suggest that LepA aids the translation of other membrane-associated proteins, namely, AmiB, a cell wall amidase, and a member of a putative taurine transport operon, MSMEG_0114 (see Fig. S2 in the supplemental material). Furthermore, we found a number of membrane proteins whose levels increased in the absence of LepA (Data Set S1). We hypothesize that these changes are a compensatory response to the loss of the major porin family in
M. smegmatis, either through transcriptional changes or some unknown regulatory mechanism. Decreased nutrient uptake due to loss of porins may trigger the upregulation of a number of alternative processes to compensate for nutrient import. Alternatively, if LepA works in concert with the Sec translocon, loss of LepA could skew the balance of ribosome activity toward some other secretion machinery.
A number of observations from genome-wide transposon screens in mycobacteria suggest that LepA’s function centers around membrane processes. In
M. tuberculosis,
lepA becomes essential for growth in the absence of
ponA1, a prominent cell wall biosynthetic enzyme in
M. tuberculosis (
50). Loss of cell wall enzymes like PonA1 may result in collateral perturbation to the membrane proteome, causing LepA to become indispensable. Further, while
lepA is nonessential for
in vitro growth of
M. tuberculosis,
lepA transposon mutants grow poorly during murine infection (
51).
M. tuberculosis requires a number of different membrane complexes during infection (
52,
53), and we speculate that LepA-dependent synthesis of membrane processes is required for survival
in vivo. Finally, there is evidence both in mycobacteria and
E. coli that the synthesis of outer membrane porins and cell wall amidases is tightly regulated (
34,
54). LepA function in mycobacteria appears to be critical for maintenance of the mycomembrane, the complex network of lipids and proteins that coordinate cellular processes ranging from cell division to nutrient transport (
Fig. 4) (
5).
How does LepA control the abundance of MspA? Unlike in
E. coli, LepA does not alter the distribution of ribosomal subunits and assembled ribosomes in mycobacteria (
28,
29).
In vitro, LepA increased translation, indicating that it functions with assembled ribosomes. While we cannot rule out specific ribosomal protein defects from altered biogenesis (
29), LepA’s role in porin synthesis and our observations
in vitro lead us to conclude that, in a mycobacterial cell, LepA acts during translational elongation. LepA has been shown to have 70S-dependent GTPase activity, and its binding affinity to the ribosome is increased by its C-terminal domain (
55). The GTPase mutant’s lack of activity
in vitro and inability to complement the Δ
lepA phenotype or increase porin abundance supports the model that LepA acts at the 70S ribosome to augment porin levels (
Fig. 4D). As structural studies in
E. coli suggest there is no direct interaction between LepA and mRNA, some other mechanism, such as sensing ribosomal conformation or interaction with secretion machinery, likely underlies LepA function. Alternatively, mycobacterial LepA may have a distinct interaction with the mycobacterial ribosome relative to the well-studied LepA homologs.
We propose two possible roles for LepA as a translational elongation factor in
M. smegmatis: (i) in quality control for abundant mRNAs and (ii) in membrane protein synthesis. In
M. smegmatis,
mspA is one of the most abundant transcripts in the cell (Table S2) (
56). Often, highly transcribed mRNAs are also highly translated in bacteria (
57). While we do not find evidence of a relationship between RNA levels and LepA influence, we think it is appropriate to discuss the implications of translational control of an abundant message. As LepA competes with EF-G for ribosome binding (from work in
E. coli), it is possible that, at least for abundant messages, LepA acts as a quality control elongation factor that facilitates the continued and efficient synthesis of messages that demand a large pool of ribosomes. Loss of this factor could impact the translation rate and, thus, cotranslational folding for highly abundant messages. In the case of porins, this delay might exacerbate a nascent polypeptide’s ability to be recognized by the secretory system or other posttranslational factors. While LepA may have roles beyond the translation of membrane proteins, studies of translation in other organism suggest a necessity for strict translational control of membrane proteins (
58–61). LepA may also therefore be involved directly at the Sec translocon. Given the complexity of the mycobacterial cell wall, we hypothesize that LepA is one of a number of regulatory mechanisms that interact to properly synthesize cell wall-associated proteins. Our findings also align with observations from eukaryotic organelles, whereby LepA appears to influence the levels of critical membrane-bound respiration and photosystem complexes (
22,
23,
27,
62).
Pathogenic mycobacteria are intrinsically resistant to antibiotics and require combination therapies. Additionally, acquired drug resistance is a growing concern that demands urgent efforts to develop new treatments. For mycobacteria, defining the mechanisms by which the cell builds its membrane proteome is critical to identifying pathways that encode innate antibiotic resistance. We expect that mycobacteria utilize other forms of posttranscriptional control to build their membrane proteome. Thus, better understanding the interplay of cellular regulatory networks will inform antibiotic discovery efforts.
MATERIALS AND METHODS
Bacterial strains and growth conditions.
M. smegmatis strains were inoculated from frozen stocks into Middlebrook 7H9 medium supplemented with 0.2% glycerol, 0.05% Tween 80, and ADC (5 g/liter bovine serum albumin, 2 g/liter dextrose, 3 µg/ml catalase) and grown at 37°C. Appropriate antibiotics or inducing agents were used at the following concentrations in M. smegmatis: nourseothricin (Nat; 20 μg/ml), zeocin (Zeo; 20 μg/ml), kanamycin (Kan; 25 μg/ml), hygromycin B (Hyg; 50 μg/ml), and anhydrous tetracycline (aTc; 100 ng/ml). Transformations, performed for the construction of M. smegmatis strains, were plated onto LB agar plates supplemented with the appropriate antibiotic. Unless otherwise specified for an experiment, strains were grown to log phase (optical density at 600 nm [OD600] of 0.3 to 0.8) without antibiotics. For cloning purposes, E. coli strains were grown in LB broth or on LB agar with antibiotics as follows: Nat (40 μg/ml), Zeo (50 μg/ml), Kan (25 μg/ml), and Hyg (100 μg/ml).
Bacterial strain construction.
All bacterial strains constructed in this study can be found in Table S2 in the supplemental material. Description of the plasmids, primers, and recombinant DNA used to construct the strains can be found in Tables S3, S4, and S5, respectively. Generally, all plasmids were constructed by restriction digestion of the parental vector (with the desired antibiotic resistance gene and phage integration gene for
M. smegmatis propagation) and all inserts were prepared by amplifying gene fragments with 18- to 25-bp Gibson assembly overhangs. Vector and insert combinations were ligated together by Gibson isothermal assembly (
63). Plasmids were isolated from
E. coli and insert orientation was sequenced via Sanger sequencing. Sequencing reactions were carried out with an ABI3730xl DNA analyzer at the DNA Resource Core of Dana-Farber/Harvard Cancer Center.
Deletion mutants. The
lepA mutant, the Δ
lepA::
zeo strain (HR334), and the
mspA mutant, the Δ
mspA::
zeo strain (HR329), were previously constructed (
18). The Δ
lepA::
zeo Δ
mspA::
hyg strain (SF789) was constructed with HR334 as the parental strain, using double-stranded recombineering. For recombination, a linear double-stranded DNA (dsDNA) fragment was generated by amplifying the following fragments: the 500-bp upstream region of
mspA, the 500-bp downstream region of
mspA, and a
lox-
hyg-
lox fragment. The three fragments were ligated together using Gibson assembly. The deletion cassette was transformed into a Δ
lepA recombineering strain as previously described (
64,
65) and plated on Hyg to select for double mutants.
lepA and mspA alleles. Plasmid pSF121, used for lepA complementation, was generated using a parental vector (pCT94) that integrates into the L5 phage site and is marked with a Kan resistance (kan) gene. The vector was digested with XbaI and HindIII (New England BioLabs, Ipswitch, MA). lepA and its 5′ untranslated region (UTR; 300 bp upstream) were amplified via PCR. The vector and insert were gel extracted and ligated using Gibson assembly. The complemented strain (SF178 ΔlepA::zeo L5::lepA-kan) and the marked mutant ΔlepA strain (SF181 ΔlepA::zeo L5::empty-kan) were derivatives of HR334. Plasmid pSF417, used for lepA complementation in the CRISPRi experiments, was generated using a parental vector (pCT204) that integrates at the Tweety (Tw) phage integration site and is marked with a Hyg resistance (hyg) gene. The lepA insert was amplified from the M. smegmatis genome as above, and the parental vector was linearized by digesting with SspI and NdeI. The resulting vector was assembled as above.
For microscopy, mspA was fused to monomeric red fluorescent protein (mRFP), with no linker, and constitutively expressed from PrpsA, in a Tw-integrating vector with a nourseothricin resistance (nat) gene. The mspA-mRFP vector was built using Gibson assembly and transformed into ΔlepA and the complemented strain for localization of MspA-mRFP.
Candidate-luciferase fusions. All luciferase reporters were generated from the parental vector CT250, a Tw-integrating vector. The vector was linearized using NdeI and HindIII, to preserve the upstream promoter (P
rpsA). For each reporter, the candidate gene was amplified from the
M. smegmatis genome and the luciferase gene was amplified from a plasmid (pJR976) (
44). Using Gibson overhangs that contained glycine-serine-glycine (GSG) linkers, each reporter vector was constructed using Gibson assembly. Each reporter vector was transformed into Δ
lepA and the complemented strain.
Porin knockdown constructs. Knockdown of each porin was accomplished using mycobacterial CRISPRi, with knockdown systems constructed as previously described (
44). Porin knockdown vectors were created by annealing oligonucleotides for each porin and ligating these fragments into a linearized vector (pJR965, digested by BsmBI), containing the mycobacterial CRISPRi system. Knockdown vectors were cotransformed into HR334 with pSF417 or pSF418.
Calcein acetoxymethyl ester (AM) staining.
Strains were grown to log phase and stained with 0.5 µg/ml of calcein AM (Invitrogen, Carlsbad, CA) for 1 h. Strains were analyzed by flow cytometry on a MACSQuant (VYB excitation: 488 nm; emission filter: 525/50) in the same manner as previously described (
18). Median fluorescence was used from each replicate to compute an overall mean fluorescence intensity.
Kill curves.
Strains were grown to mid-log phase, diluted to an OD
600 of 0.05 and treated with 10× MICs of the following drugs: rifampin (20 µg/ml), isoniazid (40 µg/ml), vancomycin (4 µg/ml), and linezolid (500 ng/ml). Survival was assessed over time as described previously (
18).
Drug susceptibility assays.
Drug susceptibility was determined using a MIC assay, as described previously (
66). In 96-well plates, strains were diluted to 0.005 and tested in biological triplicate in serial dilutions of tetracycline (Sigma-Aldrich, St. Louis, MO), clarithromycin (Sigma-Aldrich, St. Louis, MO), chloramphenicol (Sigma-Aldrich, St. Louis, MO), amikacin (Sigma-Aldrich, St. Louis, MO), and erythromycin (Santa Cruz Biotechnology, Santa Cruz, CA). The highest concentrations of each drug tested were 4 µg/ml for tetracycline, 4 µg/ml for clarithromycin, 320 µg/ml for chloramphenicol, 3.2 µg/ml for amikacin, and 16 µg/ml for erythromycin. Plates were agitated at 37°C for 21 h. To determine MICs for each condition, 0.0002% resazurin was added to each well and plates were agitated at 37°C for 3 h. The first well with no growth (blue) in each concentration gradient was considered the MIC. A biological replicate, in this case, is considered a single row in a 96-well plate of drug and bacterial incubation, using bacteria from the same culture.
Purification of mycobacterial LepA.
M. smegmatis LepA and LepA H109A were each cloned with an N-terminal 6× His tag using Gibson assembly and expressed from pET28a in BL21
E. coli, as previously described for
E. coli LepA (
28). Briefly, 200 ml of log-phase culture was induced with 1 mM isopropyl-β-
d-1-thiogalactopyranoside (IPTG) for 4 h at room temperature. Cells were harvested at 5,000 ×
g for 10 min, and pellets were frozen at −80°C overnight. The pellet was lysed at room temperature (RT) for 30 min in BugBuster 10× protein extraction reagent (Millipore Sigma, St. Louis, MO). Cell lysates were clarified by centrifugation at 15,000 ×
g for 30 min at 4°C. Lysate was brought up to 30 mM imidazole, pH 7.6, and His-tagged LepA was extracted via an Ni-nitrilotriacetic acid (Ni-NTA) column (New England BioLabs, Ipswitch, MA) purification. Beads were collected in plastic columns with a 10-ml bed volume and washed with 4 × 10 ml wash buffer (50 mM Tris-HCl [pH 7.6], 300 mM KCl, 5% glycerol, 6 mM β-mercaptoethanol [BME], 30 mM imidazole). One-milliliter elution fractions were collected using elution buffer (50 mM Tris-HCl [pH 7.6], 40 mM KCl, 5% glycerol, 6 mM BME, 200 mM imidazole) and analyzed via SDS-PAGE. The cleanest elution fractions were pooled and dialyzed into 6 liters (3 × 2 liters) of storage buffer (50 mM Tris-HCl [pH 7.6], 50 mM KCl, 5% glycerol, 6 mM BME) using dialysis cassettes with a 10-kDa molecular weight cutoff (MWCO). Aliquots (10 μl) were flash-frozen with liquid nitrogen and stored at –80°C. LepA protein concentration was calculated using the Qubit protein assay kit (Thermo Fisher, Berkeley, MO).
In vitro translation.
To assess the effect of LepA on translation, in vitro translation reactions were prepared with purified mRNA. Plasmid pSF741 was used in a HiScribe T7 in vitro transcription kit (New England BioLabs, Ipswich, MA) to generate Venus mRNA. A master mix of purified Venus mRNA (500 ng per reaction) and PURExpress (New England BioLabs, Ipswich, MA) components was prepared in duplicate reactions with purified M. smegmatis LepA or LepA H109A (300 ng per reaction). When no LepA was added to the reaction, an equal volume of storage buffer was added in place of protein. Reactions were carried out in 20 μl in a 384-well plate for 4 h at 37°C, and fluorescence (measured at an excitation of 505 nm and an emission of 540 nm) was collected on a SpectraMax M2 microplate reader.
Ribosome analysis.
Preparation of mycobacterial ribosomes was performed as previously described for
E. coli (
67), yet optimized for
M. smegmatis. A culture of 500 ml of cells was grown to mid-log phase, filtered over 0.22-μm, 90-mm membranes (Millipore Sigma, St. Louis, MO) on a fritted glass microfiltration apparatus (Kimball-Chase, Rockwood, TN), and scraped into liquid nitrogen. A 500-μl aliquot of lysis buffer (20 mM Tris [pH 8], 10 mM MgCl
2, 100 mM NH
4Cl, 5 mM CaCl
2, 0.4% Triton X-100, 0.1% NP-40, 34 mg/ml chloramphenicol, 100 U/ml RNase-free DNase I) was added to the cell scrapes. Frozen cells and lysis buffer were ground in a Retsch 400 mixer mill using 10-ml grinding jars and 12-mm grinding balls at 15 Hz for 5 × 3 min. Cell lysates were thawed and clarified at 15,000 ×
g for 15 min at 4°C. Aliquots of 250 μl of lysate were layered onto a 10 to 40% linear sucrose gradient. The sucrose gradients were spun in a Beckman ultracentrifuge at 150,000 ×
g for 2.5 h at 4°C. The gradients were fractionated and analyzed using a gradient fractionator (BioComp Instruments, Inc., NB, Canada).
Proteomics and RNA sequencing.
Cultures of 60 ml of each strain were grown to log phase (OD
600 ∼0.4) and split into two parts to extract protein and RNA separately. Both aliquots were spun at 5,000 ×
g for 10 min. For proteomics, cells were resuspended in 500 μl of urea lysis buffer (8 M urea in 50 mM Tris [pH 8.2], 75 mM NaCl, Roche complete EDTA-free protease inhibitor cocktail tablet) and subjected to bead beating for 4 × 45 s with 3 min on ice in between. Cell lysates were spun down at 20,000 ×
g for 10 min at 4°C and the supernatant was isolated for proteomics sample preparation. For RNA sequencing, RNA was isolated as described previously (
68), depleted for rRNA using RiboZero (Epicenter, Madison, WI), and prepared for sequencing using KAPA stranded RNA-Seq library preparation kit (Millipore Sigma, St. Louis, MO).
Quantitative proteomics. Samples for quantitative proteomics experiments were processed as described previously (
69). Briefly, three biological replicates of each strain’s lysates were reduced with 5 mM dithiothreitol (DTT), alkylated with 10 mM iodoacetamide (IAA), and digested with endoproteinase Lys-C (Wako Laboratories) for 2 h at a 1:50 enzyme to substrate ratio at 30°C, followed by an overnight digestion with trypsin (Promega, Madison, WI) at a 1:50 enzyme to substrate ratio at 37°C. Reactions were quenched with neat formic acid (FA) to a final concentration of 1%. Digests were desalted using tC18 SepPak reversed phase cartridges (Waters, Milford, MA) following the manufacturer’s protocol. A tandem mass tag (TMT) isobaric labeling strategy was used for this experiment. An aliquot (50 μg) of each of the 9 samples were labeled by TMT10plex reagent (Thermo Fisher Scientific, Waltham, MA) following the manufacturer’s protocol. A pooled reference standard was generated by mixing equal amounts of each of the nine samples and included in the tenth channel of the TMT10plex. Labeling efficiency was assessed prior to quenching the reactions. Once sufficient (>99%) labeling efficiency was achieved, reactions were quenched and samples were mixed together. Combined sample was desalted using tC18 Sep-Pac reversed-phase cartridges, and the eluate was dried down completely. Sample was reconstituted and fractionated on a Zorbax 300 Extend-C
18 4.6- by 250-mm column (Agilent Technologies, Santa Clara, CA), as described previously (
69). Fractions were collected every minute during the gradient and further concatenated into a total of 24 fractions that were analyzed on a Q Exactive Plus mass spectrometer (MS) coupled to an EASY-nLC 1200 ultra-high-performance liquid chromatography (UHPLC) system (Thermo Fisher Scientific, Waltham, MA). One microgram of each of the fractions was injected on a 75-μm ID Picofrit column (New Objective, Woburn, MA) packed with Reprosil-Pur C
18-AQ 1.9 μm beads (Maisch, GmbH) in-house to a length of 22 cm. Sample was eluted at a 200-nl/min flow rate with solvent A of 0.1% FA–3% acetonitrile (ACN), solvent B of 0.1% FA–90% ACN and a gradient of 2 to 6% B in 1 min, 6 to 30% B in 84 min, 30 to 60% B in 9 min, 60 to 90% B in 1 min, and a hold at 90% B for 5 min. MS data were acquired in data-dependent mode with MS1 resolution of 70,000 and automatic gain control (AGC) of 3e6. MS/MS was performed on the most intense 12 ions with a resolution of 35,000, AGC of 5e4, isolation width of 1.6 amu, and normalized collision energy of 29. Data were extracted and searched against an
M. smegmatis database using Spectrum Mill MS proteomics workbench (Agilent Technologies, Santa Clara, CA). Extracted spectra were searched using carbamidomethylation of cysteines and TMT labeling of N termini and lysine residues as fixed modifications and methionine oxidation, asparagine deamidation, and protein N-terminal acetylation as variable modifications. Spectrum to database matching was controlled with peptide level false discovery rate (FDR) of less than 1%. Peptides were rolled into protein groups and subgroups in Spectrum Mill with a protein level FDR of 0%. Protein summary export consisting of a list of quantified proteins with a reporter ion ratio of every TMT channel to the pooled reference channel was generated for quantitation of proteins. TMT10 reporter ion intensities were corrected for isotopic impurities in the Spectrum Mill protein/peptide summary module using the afRICA correction method, which implements determinant calculations according to Cramer's rule and correction factors obtained from the reagent manufacturer’s certificate of analysis (
https://www.thermofisher.com/order/catalog/product/90406) for lot number SE240163 (
70). Proteins identified with 2 or more peptides were used for further statistical analysis. Comparisons of protein levels in each strain were assessed for significance using a two-sample moderated
t test with an adjusted
P value threshold of less than 0.05 for assessing significantly altered proteins. For visualization purposes in
Fig. 2, the protein and RNA ratios associated with LepA were excluded from the volcano plot.
RNA sequencing. Samples were sequenced on an Illumina HiSeq 2500 in paired-end mode with a read length of 125 bp. Approximately 4 million reads were collected for each sample. Reads were mapped to the genome sequence of
M. smegmatis mc
2 155 as a reference genome using Burrows-Wheeler aligner (BWA) (
71). A Python script was used to separate reads in .sam files that mapped to the positive strand and negative strand of the chromosome. Then reads mapping to each open reading frame (ORF) (in a strand-specific manner) were tabulated. The raw read counts were converted to fragments per kilobase per million reads (FPKMs) by dividing by gene length (in base pairs) and total reads in the sample and scaling up by 10
9. For analyses of differential gene expression, DESeq2 (
72) was used to estimate log fold changes according to a hierarchical model based on the negative binomial distribution, and
P values were calculated via a Wald test as a measure of significance.
P values were adjusted for a false-discovery rate (FDR) of 5% over all genes by the Benjamini-Hochberg procedure.
Luciferase assays. Strains were grown to log phase and luciferase assays were conducted using Nanoglo luciferase assay system (Promega, Madison, WI). Briefly, 100 μl of cells was mixed with 100 μl of Nanoglo reagent (prepared as the kit protocol described). Within 2 min, luminescence measurements were taken in a TECAN Spark 10M plate reader with an integration time of 1,000 ms. The OD600 was also measured in each well, and luminescent values were normalized by OD600 to obtain relative luminescence values.
Porin knockdown and contribution to LepA phenotype.
Strains were grown to log phase, diluted back into medium with or without aTc, and allowed to grow for 15 h to reach log phase. Cells were stained with calcein AM and analyzed by flow cytometry in biological triplicate, as described above.
Ethidium bromide uptake assay.
Strains were grown to log phase, washed with phosphate-buffered saline (PBS)-0.4% glycerol (PBS-G) and prepared at an OD600 of 0.8. An aliquot of 100 μl of each strain was mixed with 100 μl of 4 μg/ml of ethidium bromide (prepared in PBS-G) in a 96-well plate. Fluorescence was measured in a TECAN Spark 10M plate reader, using an excitation of 520 nm and emission of 600 nm.
Fluorescence microscopy and image analysis.
GroEL-Dendra2 and MmpL3-Dendra2 strains were a gift from the Mycobacterial Systems Resource (principal investigator, K. Derbshire). Live still imaging of MspA-mRFP and Dendra2 fusion strains was performed using a Nikon TI-E inverted, wide-field microscope equipped with a Plan Apo 100×, 1.45-numerical-aperture (NA) objective, Spectra X LED light source, and an Andor Zyla sCMOS camera. M. smegmatis samples were spotted on 2% agarose pads. NIS-Elements version 4.5 software was used for data acquisition and ImageJ software was used for processing.
Experimental replicates.
Unless otherwise noted all experiments were conducted at least twice, in biological triplicate.
Data analysis.
Protein functional enrichment analysis was performed using “Functional Annotation” software within the DAVID platform (
https://david.ncifcrf.gov/home.jsp). Specifically, InterPro terms considered biologically significant, and therefore visualized, refer to at least 6 proteins in the list of 80 proteins significantly altered by LepA. Statistical significance was determined using a modified Fisher’s exact test (
35). All other statistical measurements and tests are specified in the figure legends.
mRNA quantification.
mRNA was quantified as described previously (
44). Briefly, purified RNA (DNase treated) was used as the template for cDNA synthesis, following the manufacturer’s instructions with Superscript V (Life Technologies, Carlsbad, CA). RNA was removed from the reaction using alkaline hydrolysis and the cDNA was cleaned using column purification (Zymo Research, Irvine, CA). Quantitative PCR (qPCR) was performed on purified cDNA using iTaq Universal SYBR green Supermix (Bio-Rad, Hercules, CA). mRNA fold change was calculated using the threshold cycle (ΔΔ
CT) method, where porin transcript level was normalized by
sigA level in each genetic background.
Data availability.
The original mass spectra and sequence database have been deposited in the public proteomics repository MassIVE and are accessible at
ftp://[email protected] when providing the data set password “mycobacteria.”
The RNA-seq raw sequence files are deposited at BioProject accession number
PRJNA518044, and the gene expression levels (FPKMs) are deposited in GEO under accession number
GSE126130.
The further data that support these findings are available from the corresponding author upon reasonable request.
ACKNOWLEDGMENTS
We thank Allison Carey and Thibault Barbier for advice on the manuscript. We also thank Allen Buskirk for his reading of the manuscript. We thank Sarah Fortune and the Fortune lab for insightful discussion on the work. We thank Bill Neidermeyer and members of the Whelan lab for allowing us to use ribosome fractionation machinery. We are grateful to the Michael Welsh and the Walker lab for allowing us to use their protein purification equipment.
This research was supported by funding from the NIH/NIAID (U19 AI107774 and P01 AI095208 to E.J.R.). The DNA Resource Core of Dana-Farber/Harvard Cancer Center, where sequencing reactions were carried out, was funded in part by NCI Cancer Center support grant 2P30CA006516-48.
Author contributions are as follows. Conceptualization: S.R.S.F, F.G.T., E.H.R, C.L.D, and E.J.R.; Methodology: S.R.S.F, F.G.T, E.H.R, C.L.D, E.J.R, H.K., T.R.I, and S.A.C.; Investigation: S.R.S.F, F.G.T, E.H.R, I.D.W, A.W, L.W., H.K., and T.R.I.; Data Curation: S.R.S.F, H.K., and T.R.I.; Writing—Original Draft: S.R.S.F, F.G.T, E.H.R, and E.J.R.; Writing—Review & Editing: S.R.S.F, F.G.T, E.H.R, E.J.R., H.K., T.R.I., S.A.C., C.L.D, and I.D.W.
We have no conflicts of interest to declare.