Sequential phage delivery can outperform cocktails by delaying cross-resistance evolution
Data files
Apr 01, 2026 version files 595.17 KB
-
Figure_1A.csv
4.86 KB
-
Figure_1A.R
461 B
-
Figure_1B.csv
4.36 KB
-
Figure_1B.R
459 B
-
Figure_2.csv
5.61 KB
-
Figure_2.R
465 B
-
Figure_3.csv
7.40 KB
-
Figure_3.R
417 B
-
Figure_4.csv
30.46 KB
-
Figure_4.R
435 B
-
Figure_5.csv
639 B
-
Figure_6A.csv
134.52 KB
-
Figure_6B.csv
117.43 KB
-
Figure_6C.csv
134.56 KB
-
Figure_S1.csv
1.30 KB
-
Figure_S10.csv
729 B
-
Figure_S11.csv
749 B
-
Figure_S12.csv
746 B
-
Figure_S2.csv
2.71 KB
-
Figure_S3.csv
2.92 KB
-
Figure_S4.csv
1.01 KB
-
Figure_S5.csv
502 B
-
Figure_S6.csv
467 B
-
Figure_S7.csv
878 B
-
Figure_S8.csv
859 B
-
Figure_S9.csv
636 B
-
figure5.R
358 B
-
Figure6_Table_S16.R
1.55 KB
-
README.md
32.50 KB
-
Table_S10.csv
4.96 KB
-
Table_S10.R
484 B
-
Table_S11.csv
56.72 KB
-
Table_S11.R
2.66 KB
-
Table_S14.csv
689 B
-
Table_S15.csv
551 B
-
Table_S16.csv
2.85 KB
-
Table_S17.csv
19.37 KB
-
Table_S2.csv
2.08 KB
-
Table_S2.R
1.16 KB
-
Table_S3.csv
1.49 KB
-
Table_S3.R
1.10 KB
-
Table_S4.csv
450 B
-
Table_S5.csv
793 B
-
Table_S6.csv
1.48 KB
-
Table_S8.csv
4.80 KB
-
Table_S9.csv
3.55 KB
Abstract
Antimicrobial resistance has renewed interest in bacteriophage therapy, yet bacterial evolution frequently undermines treatment efficacy. Combination phage therapy is commonly implemented as simultaneous phage cocktails, but whether this is optimal remains in question. Here, we experimentally compare simultaneous versus sequential administration of two phages, an evolved λ called ‘λtrn’ and T2, on Escherichia coli K-12 under controlled laboratory conditions. Across replicated experiments, treatment outcome depended strongly on delivery strategy, dosing order, and timing. Contrary to expectations, sequential delivery consistently achieved greater and more sustained bacterial suppression than simultaneous cocktails, although only when T2 initiated the sequence. Phenotypic assays revealed that treatment differences are driven by the accessibility and timing of cross-resistance evolution. λtrn-First treatments rapidly selected for cross-resistant bacteria prior to exposure to the second phage, rendering subsequent treatment ineffective. In contrast, T2-First sequential treatments delayed or limited cross-resistance and frequently produced single-phage resistance or collateral sensitivity. Cocktail treatments showed intermediate dynamics, with cross-resistance evolving more slowly but consistently. Whole-genome sequencing identified distinct genetic routes to cross-resistance, including regulatory mutations in envZ affecting expression of the phage receptor OmpF, as well as envelope-modifying, mucoidy-associated mutations. Engineering envZ mutations into unevolved backgrounds confirmed the mutations sufficiency to confer low-cost cross-resistance. Together, these results demonstrate that phage therapy efficacy depends not only on phage composition but on how selection pressures are ordered in time, highlighting evolutionary steering as a powerful principle for multi-phage therapy design.
Dataset DOI: 10.5061/dryad.stqjq2chf
Description of the data and file structure
Here, we experimentally compare simultaneous versus sequential administration of two phages, an evolved λ called ‘λtrn’ and T2, on Escherichia coli K-12 under controlled laboratory conditions. Bacteria and phage were cocultured for 14 days, with phages administered simultaneously (cocktail) or sequentially.
Files and variables
File: Figure6_Table_S16.R
Description: R code used to generate figures 6A-C. To import data, save as an Excel file and import using R Studio's Excel import function, or import .csv file using the "read.csv" function.
- Figure 6A: which shows Growth curves showing bacterial density (OD600) measured over 24 hours for individual λ-First isolates plotted alongside the wild-type (WT; black).
- Figure 6B: Growth curves showing bacterial density (OD600) measured over 24 hours for individual T2-First isolates plotted alongside WT (black).
- Figure 6C: Growth curves showing bacterial density (OD600) measured over 24 hours for individual cocktail isolates plotted alongside WT (black).
- This code was also used to generate the maximum exponential growth rates shown in Table S16.
File: Table_S2.R
Description: R code used to perform the global multiple linear regression modeling on the initial experiment data. Modeling shows phage delivery strategy, dose timing, and starting phage shape treatment efficacy for Table S2. To import data, save the corresponding data file (name corresponds to R file name) as an Excel file and import using R Studio's Excel import function, or import .csv file using the "read.csv" function.
File: Table_S3.R
Description: R code used to perform the multiple linear regression modeling focused on multiphage treatment shows phage delivery strategy, dose timing, and starting phage shape treatment efficacy shown in Table S3. To import data, save the corresponding data file (name corresponds to R file name) as an Excel file and import using R Studio's Excel import function, or import .csv file using the "read.csv" function.
File: Table_S11.R
Description: R code used to perform the multiple linear regression modeling showing day of isolation, starting phage, replicate population, phage administration method, and phage impact on the evolution of resistance, shown in Table S11. To import data, save the corresponding data file (name corresponds to R file name) as an Excel file and import using R Studio's Excel import function, or import .csv file using the "read.csv" function.
File: Table_S10.R
Description: R code used to perform Multiple linear regression modeling on:
- The 14-day follow-up study, dataset.
- The 5-day follow-up data set.
Code notations "#14 day follow-up study" and "#5 day follow-up study" indicate which code pertains to each data set. To import data, save the corresponding data file (name corresponds to R file name) as an Excel file and import using R Studio's Excel import function, or import .csv file using the "read.csv" function.
File: Figure_2.csv
Description: Subset of bacterial population dynamics from the initial suppression experiment illustrating variation across phage treatments. Bacterial density (OD600) was measured once daily for two weeks across multiple phage suppression assays and a no-phage control. λ-First and T2-First represent opposite sequential treatment orders, while the Cocktail treatment included both phages at each dosing event. All data shown correspond to the Late dosing regime, in which phages were added at the start of the experiment and again on day 3.
Variables
- “Treatment” data type: string, indicates which treatment data is from.
- “Days” data type: integer, indicates the day of the experiment the measurement was taken on.
- “Replicate” data type: integer, indicates which replicate population the measurement was taken from.
- “OD600” data type: float, this column has the bacterial density measurement in OD600.
File: Figure_3.csv
Description: Population dynamics for treatments T2-First Late, and Phage Cocktail Late across a 14-day follow-up experiment. Bacterial density is shown in OD600, measured in duplicate via Tecan plate reader. OD values used are an average of the duplicate plate readings.
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Days: see variable definition in "Figure_2.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- OD600: see variable definition in "Figure_2.csv" above
File: Figure_5.csv
Description: Resistance profiles for three strains independently edited with an EnvZ mutation (Q45P), the strain the mutation was discovered in “R3.2”, and wild type “WT”. Resistance was quantified using the resistance coefficient (RC) protocol as described in Methods & Materials.
Variables
- Isolate: data type: string, indicates the bacterial isolate tested for resistance.
- Phage: data type: string, this column indicates the phage used to test for phage-resistance.
- RC: data type: float, this column shows the resistance coefficient (RC).
- SD: data type: float, this column shows the standard deviation (SD) of the resistance coefficient data.
- CI: data type: float, this column should the 95% confidence interval (CI) for the resistance coefficient data.
- Tukey's: data type: float, this column shows the Tukey’s HSD post hoc p value.
- Significance: data type: string, this column has asterisks indicating if the resistance estimates are significantly different than the sensitive ancestor, with p < 0.001 denoted as **, & p<0.0001 as ***.
File: Figure_4.csv
Description: Cross-resistance profiles of bacterial isolates sampled across treatments and time points. Each point represents the resistance coefficients of a single isolate. Pale symbols denote isolates not significantly different from the ancestor, while opaque symbols indicate significance in at least one resistance dimension. Color coldness reflects the day on which the isolate was sampled.
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Population: “Population” data type: string, indicates which replicate population the measurement was taken from.
- T2RC: data type: float, shows the resistance coefficient (RC) of the bacterial isolate, indicating how susceptible or resistant the isolate is to T2 phage.
- T2RC Tukey's: data type: float, shows the Tukey’s HSD post hoc analysis p-value generated after ANOVA analysis.
- T2RC Significance: data type: integer, this column is a binary with a 1 indicating a significant “T2RC Tukey’s” value, and a 0 indicated a non-significant (>0.05) p value.
- LRC: data type: float, this shows the resistance coefficient (RC) of the bacterial isolate, indicating how susceptible or resistant the isolate is to lambda-trn phage.
- LRC Tukey's: data type: float, shows the Tukey’s HSD post hoc analysis p-value generated after ANOVA analysis.
- LRC Significance: data type: integer, this column is a binary with a 1 indicating a significant “LRC Tukey’s” value, and a 0 indicated a non-significant (>0.05) p value.
- Day: data type: integer, shows the day of the experiment the bacteria was isolated from.
- Isolate: data type: string, shows the abbreviated isolate identifier, indicating the treatment, replicate population, colony, and day of isolation, with “L” for Lambda-First, “P” for phage cocktail, and “T” for T2-First (Example: L1.1D1 is Lambda-First replicate population 1, colony 1 isolated on day 1).
- Both Significance: data type: string, this column indicates whether or not “T2RC Significance” or “LRC Significance” were significant, if either “T2RC Significance” or LRC Significance” has a value of “1” then the column shows “TRUE”, if both have a value of “0” then the column shows “FALSE”/
- AlphaLT: data type: integer/float, this column has a value used to code the alpha/opacity of data points, if the “Both Significance” column has “TRUE” then the “AlphaLT” column has a value of 1, if the “Both Significance” column has “FALSE” then the “AlphaLT” column has a value of 0.25.
File: Figure_6A.csv
Description: Growth curves showing bacterial density (OD600) measured over 24 hours for individual λ-First isolates plotted alongside the wild-type (WT; black), and Lambda-first isolates in red. (L1.3 is Lambda-First Late Day 1 R1 colony #3, L3.2 is Lambda-First Late Day 3 R3 colony #2, & L5.2 is Lambda-First Late Day 3 R5 colony #2).
Variables
- Strain: data type: string, this column shows the bacterial isolate being studied,
- Replicate: data type: integer, this column shows the growth rate experiment replicate number.
- Time: data type: integer, this column shows the time in seconds that the OD600 measurement was taken. OD600 was measured every 5 minutes.
- OD: data type: float, this column shows the OD600 measurement.
File: Figure_6C.csv
Description: Growth curves showing bacterial density (OD600) measured over 24 hours for individual cocktail isolates plotted alongside WT (black). (P1.2 is Cocktail Late Day 1 R1 colony #2, P2.1 is Cocktail Late Day 2 R2 colony #1, & P3.2 is Cocktail Late Day 3 R3 colony #2).
Variables
- Strain: See variable definition in "Figure_6A.csv" above
- replicate: see variable definition in "Figure_6A.csv" above
- time: See variable definition in "Figure_6A.csv" above
- OD: See variable definition in "Figure_6A.csv" above
File: Figure_6B.csv
Description: Growth curves showing bacterial density (OD600) measured over 24 hours for individual T2-First isolates plotted alongside WT (black). (T1.1 is T2-First Late Day 2 R1 colony #1, T3.3 is T2-First Late Day 3 R3 colony #3, & T5.3 is T2-First Late Day 3 R5 colony #3).
Variables
- Strain: See variable definition in "Figure_6A.csv" above
- replicate: See variable definition in "Figure_6A.csv" above
- time: See variable definition in "Figure_6A.csv" above
- OD: See variable definition in "Figure_6A.csv" above
File: Figure_S1.csv
Description: Phage titer in Phage Cocktail Late in the first 3 days of the experiment by the host. Samples of frozen coculture from each replicate population were defrosted, and spot assays were done on 3 hosts: WT K12 E. coli, LamB- OmpF- E. coli, and fadL- OmpF- E. coli.
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Host: data type: string, this column indicates which bacterial strain the phage were spotted onto in the spot assay.
- Phage: data type: string, this column indicates which phage titer is being shown
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- LogPFU: data type: float, this column shows the log-transformed titer data.
- pfu: data type: float, this column shows the phage titer in plaque-forming units (pfu) per ml.
File: Figure_S2.csv
Description: Phage titer in T2-First Late in the first 3 days of the experiment by host. Samples of frozen coculture from each replicate population were defrosted, and spot assays were done on 3 hosts: WT K12 E. coli, LamB- OmpF- E. coli, and fadL- OmpF- E. coli.
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Host: see variable definition in "Figure_S1.csv" above
- Phage: see variable definition in "Figure_S1.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- LogPFU: see variable definition in "Figure_S1.csv" above
- pfu: see variable definition in "Figure_S1.csv" above
File: Figure_S3.csv
Description: Phage titer in Lambda-First Late in the first 3 days of the experiment by host. Samples of frozen coculture from each replicate population were defrosted, and spot assays were done on 3 hosts: WT K12 E. coli, LamB- OmpF- E. coli, and fadL- OmpF- E. coli.
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Host: see variable definition in "Figure_S1.csv" above
- Phage: see variable definition in "Figure_S1.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- LogPFU: see variable definition in "Figure_S1.csv" above
- pfu: see variable definition in "Figure_S1.csv" above
File: Figure_S4.csv
Description: . Phage titers in the 𝜆-First Late treatment on day 14 the experiment, measured on different host genotypes. T2 titer shown in blue, 𝜆trn in red, and both phages in green.
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Host: see variable definition in "Figure_S1.csv" above
- Phage: see variable definition in "Figure_S1.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- titer: data type: float, shows the phage titer in pfu/mL
- LogPFU: see variable definition in "Figure_S1.csv" above
File: Figure_S5.csv
Description: Phage titers in the Cocktail treatment on day 14 the experiment, measured on different host genotypes. T2 titer shown in blue, 𝜆trn in red, and both phages in green.
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Host: see variable definition in "Figure_S1.csv" above
- Phage: see variable definition in "Figure_S1.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- titer: see variable definition in "Figure_S4.csv" above
- LogPFU: see variable definition in "Figure_S1.csv" above
File: Figure_S6.csv
Description: Phage titers in the T2-First Late treatment on day 14 the experiment, measured on different host genotypes. T2 titer shown in blue, 𝜆trn in red, and both phages in green
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Host: see variable definition in "Figure_S1.csv" above
- Phage: see variable definition in "Figure_S1.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- ttiter: see variable definition in "Figure_S4.csv" above
- LogPFU: see variable definition in "Figure_S1.csv" above
File: Table_S2.csv
Description: Multiple linear regression modeling shows phage delivery strategy, dose timing, and starting phage shape treatment efficacy.
Variables
- Type: data type: string, this column indicates the type of treatment, “Sequential”, “Cocktail”, or “Monophage”.
- Timing: data type: string, this column indicates the dosing schedule for the treatment, “Early”, “Mid”, or “Late”.
- Starting Phage: data type: string, this column shows which phage was given as the first dose, “T2” Lambda”, or “Both”.
- OD: data type: float, shows the mean OD per replicate, calculated by taking the average of the daily OD values of per replicate across the entire duration of the experiment
File: Table_S3.csv
Description: Multiple linear regression modeling focused on multiphage treatment shows phage delivery strategy, dose timing, and starting phage shape treatment efficacy.
Variables
- Type: see variable definition in "Table_S2.csv" above
- Timing: see variable definition in "Table_S2.csv" above
- SP: Data type: string, Starting phage (SP) indicates which phage was given as the first dose (T2, Lambda, or both)
- OD: see variable definition in "Table_S2.csv" above
File: Table_S4.csv
Description: Pairwise comparison between the resistance coefficient WT benchmark , ln(WTP /WTNP), of each phage shows that T2 is a more suppressive phage under our experimental conditions.
Variables
- T2 WT Standards: T2 Killing standard used to calculate the resistance coefficient. Calculated by ln(WTP /WTNP), where WTP is the cell density after 30-minute exposure to T2 phage and *WTNP is the cell density of WT E. coli cultured without phage for the same duration.
- Lambda WT Standards: Lambda Killing standard used to calculate the resistance coefficient. Calculated by ln(WTP /WTNP), where WTP is the cell density after 30 minute exposure to lambda-trn phage and *WTNP is the cell density of WT E. coli cultured without phage for the same duration.
File: Table_S5.csv
Description: Pairwise comparisons show T2-First Sequential treatments outperform cocktail treatments. Comparisons were made between bacterial density (OD) 24 hours after the second dose was administered.
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Day: data type: integer, indicates which day the data being compared are from
- Replicate: see variable definition in "Figure_2.csv" above
- OD: see variable definition in "Figure_2.csv" above
File: Table_S6.csv
Description: Pairwise comparisons show 𝜆-First Sequential treatments are outperformed by all other treatments. Comparisons made between bacterial density (OD) 24 hours after the second dose was administered.
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Day: data type: integer, indicates which day the data being compared are from
- Replicate: see variable definition in "Figure_2.csv" above
- OD: see variable definition in "Figure_2.csv" above
File: Table_S8.csv
Description: ANOVA comparisons among all three experimental runs show no significant difference. ANOVA comparisons among the initial experiment, the 5-day follow-up, and the 14-day follow-up show no significant difference in OD across experimental runs. The mean bacterial density for each replicate population averaged across the full 14-day experiment was used for this analysis
Variables
- Treatment ID: data type: string, Abbreviated treatment name
- Treatment: see variable definition in "Figure_2.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- Mean OD: data type: float, Bacterial density for each replicate population averaged across the full 14-day experiment
- Experiment: data type: string, indicates which experiment the data being compared is from (Initial Experiment, 5-day follow-up, or 14-day follow-up)
File: Table_S9.csv
Description: Pairwise comparisons show T2-First Sequential treatments outperform cocktail treatments in follow-up studies
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- OD600: data type: float, Average OD for a given replicate across the entire duration of the experiment
- Experiment: see variable definition in "Table_S8.csv" above
File: Table_S10.csv
Description: Multiple linear regression modeling shows phage delivery strategy, and dose timing shape treatment efficacy in follow-up studies
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- OD600: data type: float, Average OD for a given replicate across the entire duration of the experiment
- Experiment: see variable definition in "Table_S8.csv" above
- Type: see variable definition in "Table_S2.csv" above
- Timing: see variable definition in "Table_S2.csv" above
File: Table_S14.csv
Description: EnvZ mutants are significantly more resistant to both phages than the WT is.
Variables
- Isolate: see variable definition in "Figure_5.csv" above
- Isolate ID: data type: string, Abbreviated isolate identifier
- Replicate: see variable definition in "Figure_2.csv" above
- T2 RC: data type: float Resistance coefficient (RC) indicating the resistance/susceptibility of the isolate to phage T2
- Lambda RC: float Resistance coefficient (RC) indicating the resistance/susceptibility of the isolate to phage lambda
File: Table_S11.csv
Description: Multiple linear regression modeling shows day of isolation, starting phage, replicate population, phage administration method, and phage impact the evolution of resistance.
Variables
- Full: data type: string, Full treatment name
- Treatment: see variable definition in "Figure_2.csv" above
- Tube: data type: string, Unique replicate tube identifier, each replicate population has a unique identifier used throughout the experiment (Cocktail replicate population 1 = C1)
- Day: data type: integer, Day of isolation. This indicates the day of the experiment that the colony was isolated from
- Colony: Unique colony identifier. Each individual colony has a unique, non-repeated colony identifier. No two colonies share an identifier as colonies were randomly isolated.
- Phage: see variable definition in "Figure_5.csv" above
- RC: see variable definition in "Figure_5.csv" above
- SP: see variable definition in "Table_S3.csv" above
File: Table_S15.csv
Description: . Cross resistance of EnvZ mutants is not significantly different from 𝜆-First Late isolate L3.2. The table shows the ANOVA and Tukey's HSD results.
Variables
- Isolate: see variable definition in "Figure_5.csv" above
- Isolate ID: see variable definition in "Table_S14.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- T2 RC: see variable definition in "Table_S14.csv" above
- Lambda RC: see variable definition in "Table_S14.csv" above
File: Table_S16.csv
Description: Shows average maximum growth rates (𝝁Max) for each bacterial isolate, the results of ANOVA analysis comparing the maximum growth rates, and the Tukey’s HSD post hoc analysis results
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Day of Isolation: data type: integer. This indicates the day of the experiment that the colony was isolated from
- Replicate Population: data type: integer, indicates which replicate population the measurement was taken from.
- Colony: data type: integer, identifies the colony number; 3 colonies were isolated from each replicate population
- Isolate ID: see variable definition in "Table_S14.csv" above
- replicate: data type: integer, refers to the technical replicates of the bacterial isolates in the growth rate experiment
- mumax: data type: float, shows the maximum growth rate
- Block: data type: integer, indicates which experimental block the data is from
File: Table_S17.csv
Description:
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- OD: see variable definition in "Figure_2.csv" above
- LogCFU: data type: float Log10 transformed bacterial density from colony counts
- CFU: data type: float and integer (rows 470-479 are integers) Colony Forming Units (CFU), bacterial density given in CFU/mL, calculated from colony counts
File: Figure_1A.csv
Description: Subset of bacterial population dynamics from the initial suppression experiment illustrating variation across phage treatments.**** Bacterial density (OD600) was measured once daily for two weeks across multiple phage suppression assays and a no-phage control. 1A shows the λ-First sequential treatment, λ-Only monophage treatment, and the cocktail treatment.
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- OD: see variable definition in "Figure_2.csv" above
File: Figure_1B.csv
Description: Shows the T2-First sequential treatment, the T2-Only monophage treatment and the cocktail treatment. All data shown correspond to the Late dosing regime, in which phages were added at the start of the experiment and again on day 3.
Variables
- Treatment: see variable definition in "Figure_2.csv" above
- Days: see variable definition in "Figure_2.csv" above
- Replicate: see variable definition in "Figure_2.csv" above
- OD600: see variable definition of "OD" in "Figure_2.csv" above
File: Figure_S7.csv
Description: Timeseries plot showing the evolution of resistance to T2 in the T2-First bacterial populations in the first 6 days of the experiment. Each data point represents the average resistance coefficient (RC) of the three bacterial isolated from each replicate population. Error bars represent the 95% confidence interval. Missing data are from replicates where the phage went extinct.
Variables
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- RC: see variable definition in "Figure_5.csv" above
- CI: data type: float, 95% Confidence Interval (CI)
- 95% CI Upper Bound: data type: float, Resistance coefficient (RC) 95% Confidence Interval (CI) upper bound ( RC + C1)
- 95% CI Lower Bound : data type: float, Resistance coefficient (RC) 95% Confidence Interval (CI) lower bound ( RC - C1)
File: Figure_S8.csv
Description: Timeseries plot showing the evolution of resistance to λtrn in the T2-First bacterial populations in the first six days of the experiment. Each data point represents the average resistance coefficient (RC) of the three bacterial isolated from each replicate population. Error bars represent the 95% confidence interval. Missing data are from populations where the phage went extinct
Variables
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- RC: see variable definition in "Figure_5.csv" above
- CI: see variable definition in "Figure_7.csv" above
- 95% CI Upper Bound: see variable definition in "Figure_7.csv" above
- 95% CI Lower Bound: see variable definition in "Figure_7.csv" above
File: Figure_S10.csv
Description: Timeseries plot showing the evolution of resistance to λtrn in the Cocktail treatment bacterial populations in the first 6 days of the experiment. Each data point represents the average resistance coefficient (RC) of the 3 bacterial isolated from each replicate population. Error bars represent the 95% confidence interval.
Variables
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- RC: see variable definition in "Figure_5.csv" above
- CI: see variable definition in "Figure_7.csv" above
- Phage: see variable definition in "Figure_5.csv" above
File: Figure_S9.csv
Description: Time series plot showing the evolution of resistance to T2 in the Cocktail treatment in the bacterial populations in the first six days of the experiment. Each data point represents the average resistance coefficient (RC) of the three bacterial isolated from each replicate population. Error bars represent the 95% confidence interval.
Variables
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- RC: see variable definition in "Figure_5.csv" above
- CI: see variable definition in "Figure_7.csv" above
- Phage: see variable definition in "Figure_5.csv" above
File: Figure_S11.csv
Description: Timeseries plot showing the evolution of resistance to T2 in the λ-First treatment bacterial populations in the first 3 days of the experiment. Each data point represents the average resistance coefficient (RC) of the 3 bacterial isolated from each replicate population. Error bars represent the 95% confidence interval.
Variables
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- RC: see variable definition in "Figure_5.csv" above
- CI: see variable definition in "Figure_7.csv" above
File: Figure_S12.csv
Description: . Timeseries plot showing the evolution of resistance to λtrn in the λ-First treatment bacterial populations in the first 3 days of the experiment. Each data point represents the average resistance coefficient (RC) of the 3 bacterial isolated from each replicate population. Error bars represent the 95% confidence interval.
Variables
- Replicate: see variable definition in "Figure_2.csv" above
- Day: see variable definition of "Days" in "Figure_2.csv" above
- RC: see variable definition in "Figure_5.csv" above
- CI: see variable definition in "Figure_7.csv" above
File: figure5.R
Description: R code using ggplot2 to plot Figure 5; a bar plot showing the resistance coefficient scores of L3.2 and the engineered envZ mutants (edited to have the envZ mutation observed in L3.2). To import data, save the corresponding data file (name corresponds to R file name) as an Excel file and import using R Studio's Excel import function, or import .csv file using the "read.csv" function.
File: Figure_1A.R
Description: R code using ggplot2 to plot Figure 1A, which shows a subset of bacterial population dynamics from the initial suppression experiment, illustrating variation across phage treatments. 1A shows the λ-First sequential treatment, λ-Only monophage treatment, and the cocktail treatment. To import data, save as an Excel file and import using R Studio's Excel import function, or import .csv file using the "read.csv" function.
File: Figure_1B.R
Description: R code using ggplot2 to plot Figure 1B, which shows the T2-First sequential treatment, the T2-Only monophage treatment and the cocktail treatment. All data shown correspond to the Late dosing regime, in which phages were added at the start of the experiment and again on day 3. To import data, save the corresponding data file (name corresponds to R file name) as an Excel file and import using R Studio's Excel import function, or import .csv file using the "read.csv" function.
File: Figure_2.R
Description: R code used to generate Figure 2, which shows population dynamics for treatments T2-First Late, Lambda-First Late, and Phage Cocktail Late across the initial 14-day experiment. Bacterial density is shown in OD600, measured in duplicate via Tecan plate reader. OD values used are an average of the duplicate plate readings. To import data, save the corresponding data file (name corresponds to R file name) as an Excel file and import using R Studio's Excel import function, or import .csv file using the "read.csv" function.
File: Figure_3.R
Description: R code used to generate Figure 3 which shows population dynamics for treatments T2-First Late, and Phage Cocktail Late across a 14-day follow-up experiment. To import data, save the corresponding data file (name corresponds to R file name) as an Excel file and import using R Studio's Excel import function, or import .csv file using the "read.csv" function.
File: Figure_4.R
Description: R code used to generate Figure 4, which shows resistance profiles of bacterial isolates sampled across treatments and time points. Each point represents the resistance coefficients of a single isolate. Pale symbols denote isolates not significantly different from the ancestor, while opaque symbols indicate significance in at least one resistance dimension. Color coldness reflects the day on which the isolate was sampled. To import data, save the corresponding data file (name corresponds to R file name) as an Excel file and import using R Studio's Excel import function, or import .csv file using the "read.csv" function.
Code/software
R Studio Version 2026.01.0+392 (2026.01.0+392) and R version 4.5.2 (2025-10-31) were used to prepare plots, statistical analysis, and modeling.
ggplot2, and model_performance r packages were used to prepare plots (ggplot2) and perform model selection (model_performance)
