Reviewer #2 (Public Review):
This paper by Patel and Matange focuses on understanding the evolutionary response of E. coli cells to the antibiotic trimethoprim (TMP). Mutations in the gene folA, which encodes the dihydrofolate reductase (DHFR) enzyme - an established target of TMP, are known to mediate intrinsic resistance to TMP. This work shows that de-repression of the PhoQ/PhoP pathway via inactivation of MgrB, which is a negative feedback inhibitor of PhoQ, leads to TMP tolerance. Further, this response to TMP is due to the upregulation of DHFR expression. Similarly, inactivation of MgrB and a corresponding increase in the PhoQ/PhoP-regulated gene expression is a prominent mechanism for acquired colistin resistance in clinical isolates of Klebsiella pneumoniae. The authors performed adaptive laboratory evolution of E. coli under TMP selection and identified factors contributing to the transition of TMP-tolerant bacterial cells to TMP-resistant cells. At high TMP concentrations, the cells become resistant via mutations in DHFR. Cells evolved under low (sub-MIC) concentrations of TMP develop mutations inactivating the RpoS sigma factor to offset the cost of PhoQ de-repression. The data presented in the paper are clear and the conclusions are mostly valid. However, the authors need to modify parts of the main text and figure presentations for clarity and a better/more straightforward interpretation of the results by the reader. In general, the results obtained in this study explain well the evolutionary consequences of these mutations and the pathway for the acquisition of the mutations.
1) The authors find that mutations in the mgrB locus precede mutations in folA during E. coli's response to TMP. Why only sequence 5 of the 10 TMPR mutants? Was this subset chosen for sequencing based on any specific criteria? Below are some follow-up comments.
a. Do any of the mutations cause growth defects relative to the wild-type strain?
b. Line 103: What are the mutations in folA promoter region? Only mutations in the coding sequence are listed in table 1 and figure 1A.
c. Line 109: The authors speculate that IS-element insertions in the mgrB promoter region reduce its expression, maybe they can provide a reference here from previous studies that have analyzed such mutations. Also, including details of the length/size of these insertion elements within table 1 would be helpful.
d. Line 111: the phrase "stop-codon readthrough" is misleading. The authors should rephrase to clarify that the single nucleotide deletion leads to a shift in the reading frame leading to an altered protein sequence at the C-terminal end.
2) Based on growth assays including competitions, and measurements of folA gene expression in mgrB-deficient E. coli cells, the authors conclude that tolerance to TMP is caused by PhoP-dependent upregulation of DHFR.
a. The authors should rewrite the text (lines 143-155) to make the experimental design of the competitions more obvious to the reader. Indicating either within the figure legend or main text what ∆mgrB/total means would definitely make analysis of the figure and results easier for the reader The reader needs to go to the materials section to get a full understanding how exactly this experiment was performed.
b. In Figure 1C, the IC50 value for ∆phoP is similar to that of wild type. If PhoP-dependent expression of folA important for TMP tolerance/resistance, shouldn't we expect to see a lower IC50, similar to that of ∆mgrB∆phoP? Intriguingly, the data for wild type in Figure 1C appears to be in conflict with the data in Figure 3B, please clarify.
c. In Figure 1D, it is hard to figure out the exact strains and conditions of each competition. For instance, the ratios 10:1, 100:1 and 1000:1 needs to be clearly labeled, "wild type: mgrB" or "wild type: specific mutant" as applicable, the label on the X-axis is misplaced. Does "WmgrB" refer to ∆mgrB? If yes, change to ∆mgrB. Fitness values need a label or put into a table.
d. Line 172: incorrect figure citation, replace Figure 2B with 2A.
e. Lines 180-181: Only 5 out of the 10 TMPR isolates were sequenced and found to have mutations in the mgrB locus. In the absence of sequencing data confirming such mutations in TMPR 6-10 isolates, the increased levels of DHFR cannot be attributed to loss of mgrB.
f. In Figure 2C, it would be helpful to show the GFP fluorescence data for the single deletions, ΔphoP and ΔrpoS, to further support the claim that TMP tolerance via DHFR upregulation is PhoP dependent. In addition, the X-axis should specify the promoter reporter that was used.
g. Lines 181-183: reference for the previous work on W30G folA is missing.
h. In Figure 2, there is a discrepancy in the level of DHFR observed for both TMPR2 and 3 isolates in panels D and E - the DHFR protein levels are much higher in panel E. Can the authors explain this discrepancy, especially given the W30G mutation in TMPR3 (expected to show reduced levels of DHFR)? Is the same amount of protein loaded in both experiments? If so, why are the levels of protein different (and vastly different for TMPR3)? Better quantification of the western blots depicting the signal for the replicates would be helpful.
3) The data presented here also show that mgrB and folA mutations act in synergy in TMP resistant E. coli.
a. It would be useful to the reader to include a table listing the MIC values in Figure 3. The plate images showing the E-tests are difficult to read and less helpful in interpreting the MICs and can be moved to the supplement.
b. In Figure 3E (and lines 234-238), what was the strain background used for DHFR overexpression? The details are missing from the paper.
4) To follow the adaptive pathway for TMP resistance, the authors sequenced genomes of TMP-resistant isolates.
a. Line 283: How many strains were sequenced at each time point? "3 to 5" is confusing.
b. In Figure 4, the data points/symbols and lines are hard to read in both panels A and B. These graphs can be replotted with open symbols or different colors to help the reader analyze the figure much more easily.
c. Overall, it is still unclear how folA expression is regulated by PhoP regulation. An alternate hypothesis is that loss of MgrB may influence folA gene expression in a PhoP independent manner. Have the authors ruled out this possibility?