Mini Review ARTICLE
Bacteriophage Production Models: An Overview
- 1Laboratorio de Microbiología, Instituto de Biología, Pontificia Universidad Católica de Valparaíso, Valparaíso, Chile
- 2Laboratorio de Biotecnología, Instituto de Nutrición y Tecnología de los Alimentos, Universidad de Chile, Santiago, Chile
- 3Instituto de Ciencias Biomédicas, Facultad de Ciencias de la Salud, Universidad Autónoma de Chile, Santiago, Chile
The use of bacteriophages has been proposed as an alternative method to control pathogenic bacteria. During recent years several reports have been published about the successful use of bacteriophages in different fields such as food safety, agriculture, aquaculture, and even human health. Several companies are now commercializing bacteriophages or bacteriophage-based products for therapeutic purposes. However, this technology is still in development and there are challenges to overcome before bacteriophages can be widely used to control pathogenic bacteria. One big hurdle is the development of efficient methods for bacteriophage production. To date, several models for bacteriophage production have been reported, some of them evaluated experimentally. This mini-review offers an overview of different models and methods for bacteriophage production, contrasting their principal differences.
Bacteriophages have come into focus of scientific research, as they play an important role in almost every microbial community. As viral predators of bacteria, they have a substantial influence on microbial populations and dynamics in different environments. There are several reviews dealing with the role of bacteriophages in different habitats such as the seas or the human body (Clokie and Mann, 2006; Wahida et al., 2016; Łusiak-Szelachowska et al., 2017). Since their discovery more than 100 years ago, separately by Frederick Twort and then by Felix D’Herelle (Salmond and Fineran, 2015), bacteriophages have been used in eastern European countries for medical treatment of bacterial infections, whereas in the rest of the world antibiotics were the protagonists (Myelnikov, 2018). Nowadays, as infections with multi-resistant bacteria have become a worldwide threat (Zaman et al., 2017), patients from all over the world are treated at the Eliava Institute of Bacteriophages, Microbiology, and Virology, in Tbilisi, Georgia, which has perhaps the longest experience in bacteriophage therapy (Kutateladze and Adamia, 2008), and also in the Bacteriophage Therapy Unit of the Ludwik Hirszfeld Institute of Immunology and Experimental Therapy in Wrocław, Poland (Międzybrodzki et al., 2012). The application of bacteriophages could be not only a valuable solution in the medical sector, but as well as in other fields where bacteria can have a negative impact.
Some companies in the United States such as OmniLytics Inc. (Sandy, UT, United States) and Intralytix Inc. (Baltimore, MD, United States) have developed different bacteriophage products for the application as disinfectants in the food industry that can be used against Salmonella, Escherichia coli and Listeria monocytogenes. In Europe, a Dutch company Micreos BV (Wageningen, Netherlands) also marketed bacteriophage products against Salmonella and E. coli and a German company, Fink Tec (Hamm, Germany), targeting E. coli (Moye et al., 2018). Wider application of bacteriophages is expected in the food value chain including agriculture and aquaculture, where a broad spectrum of diverse plant and fish pathogens causes significant economic losses (Buttimer et al., 2017; Doss et al., 2017).
Although some bacteriophage products are already being commercialized, an effective, constant and controllable process for bacteriophage production has yet to be achieved. The production of phages in laboratories can be considered a routine process, and protocols are well defined; however, these processes are not easily scaled up. Industrial entities have the main interest to obtain reliable methods for phage production that allow scale-up of the process. However, the solution is not easy, due to the biological nature of the system and the diverse types of interactions that occur between phages and bacteria.
There have been several attempts to generate reliable methods for bacteriophage production. Some researchers have used a theoretical approach with simulation models, while others took the practical approach through experiments. This mini-review examines selected examples of both approaches, contrasting their principal differences.
Generalities in Bacteriophage Production
The biological nature of bacteriophages forces their reproduction in the host cell. Therefore, a method for bacteriophage production requires a production process involving at least two operating units, growth of the host bacteria and bacteriophage propagation (or infection). It is important to consider basic parameters for bacterial growth and phage infection, such as the selected substrates for the bacterium and the optimal temperature, both for bacterial growth and phage infection, since these factors may influence the infectivity of phages (Tokman et al., 2016). Similarly, is important to know the biology of the phage to be produced, including the different infection parameters such as adsorption rate, burst size and latent period; however, as will be discussed later, these parameters can change depending on the infection conditions (Santos et al., 2014). More importantly, is recommended to have a deep understanding of the specific interactions that may occur between the bacterial host and the selected phage, such as the presence of a CRISPR-cas system in the bacterium, because these factors may have a strong influence on the phage infection process (Levin et al., 2013). It is also recommended to select a non-virulent bacterial strain as host. Bacteriophage production at the industrial level will require large quantities of the host bacterium, so avoiding the use of virulent drug-resistant, and especially multi-resistant pathogens should be mandatory in a phage production process (Torres-Barceló, 2018). The same applies for bacteria carrying prophages, because they could be induced during the process, altering the final outcome (Stewart and Levin, 1984).
A reliable process for large-scale bacteriophage production can be very elusive, as data obtained in a laboratory are not always useful for scaling up biological processes (Kwok, 2010). Researchers have tried to fill this gap mainly through studies on bacteriophage production based on computer simulations, some of them validated experimentally. Here, we will analyze first theoretical studies focused on phage production models and then selected studies that have been validated experimentally. All the cases agree with the assay criteria for further purification and validation of a bacteriophage-based product, and some of them are included in both sections (Santos et al., 2014; Nabergoj et al., 2018a).
Theoretical Models for Bacteriophage Production
To describe a process of phage production through a mathematical model it is important to define the kinetic parameters to include in the model. The three basic parameters for phage production are the populations of susceptible uninfected bacteria, phage-infected bacteria, and free phages (Krysiak-Baltyn et al., 2016). Starting from this, different models have included additional variables such as resistant uninfected bacteria (Santos et al., 2014; Chaudhry et al., 2018) or multiple bacterial species (Levin et al., 1977). All these populations interact controlled by kinetic parameters associated with bacterial growth and phage infection. It is considered well known which constants are important for bacteria; however, this is still in discussion for bacteriophages. There is consensus that adsorption constant, latency period and burst size are important variables to consider; however, their relevance in the model varies between different studies. Moreover, different authors use different nomenclature to define the kinetic parameters, which is one of the main difficulties to establish comparisons between different models and to unify the general knowledge on this topic. For instance, the adsorption rate of phages (indicator of phage particles adsorbed to bacteria) is commonly referred to by the symbol “δ”; however, Beretta and Kuang (1998) used the symbol “K,” which can be also the symbol for Monod’s constant of substrate specificity “Ks.” Other examples of different nomenclature can be found in Table 1. As in other biological process, it is expected that authors working in the field of phage-bacteria growth models agree on a specific algebraic vocabulary or include a clear explanation of terms and units in their articles and a clear nomenclature, as recently stated by Krysiak-Baltyn et al. (2018). Based on the nomenclature used by other authors (Table 1) we propose the use of Greek characters to name the different kinetic parameters in phage reproduction. Burst size can be symbolized by β, adsorption rate by δ, eclipse time by ε and phage decay rate by λ. The only exceptions would be phage concentration, which is commonly known as “P,” and latency time, known as “L.” Uniformity in this mathematical language will facilitate the understanding and data mining for future academic or industry reviewers.
Beginning with Campbell (1961), many efforts have been made to describe models of bacteriophage production, describing phage population behavior under several conditions and methods. Table 1 summarizes different phage production models, given as differential or integral equations (depending on each author’s decision), mentioning specific considerations for each model.
Phage production models are generally consistent in describing phage population change over time. This may be represented as a kinetic change in phage particles or plaque forming units (PFU) per unit of time, final concentrations obtained after a batch process, or during a period of time in a continuous process. In spite of the general consensus, these models differ in several statements. Models proposed by Campbell (1961) and Beretta and Kuang (1998) are consistent in balancing phage particles with generation terms (liberation of bacteriophage particles per unit of time) and loss of free bacteriophage due to adsorption or decay rates; these models are useful due to their simplicity and the use of standard phage growth parameters like adsorption rate, burst size and latency time, and are a quick way to simulate batch production processes, but they could not fit processes like resistant bacteria populations or phage evolution over time. These models also tend to underestimate the influence of parameters such as burst size and latent period, while more recent models have shown the importance of these parameters and how they can vary depending on other factors (Santos et al., 2014; Nabergoj et al., 2018b).
An interesting model recently proposed by Santos et al. (2014) considers the influence of bacterial growth rate on the phage adsorption constant and a normal distribution equation that rules the values of latent period, taking into account the variability in these parameters. This model has proven to be very useful because it provides the opportunity to evaluate the influence of the substrate in phage production, and including the bacterial growth rate in the model offers an indirect tool to consider the physiological state of the bacteria during the process. The dependence of bacteriophage infection parameters on bacterial growth rate was later also explored by other authors (Krysiak-Baltyn et al., 2018; Nabergoj et al., 2018b); Nabergoj and colleagues found that burst size increased linearly with bacterial growth rate, while adsorption constant and latent period decreased.
Other models have explored the influence of multiple bacterial species, and the occurrence of bacterial resistance (Levin et al., 1977; Santos et al., 2014; Chaudhry et al., 2018). Although the aim of these studies was not always to develop methods for phage production, they are useful to describe potential situations than can occur during the process. These models include variables associated with bacterial resistance selection and reversion rates as a function of the bacterial population (availability of less or more susceptible bacteria over time), defining conditions in which susceptible and resistant bacteria can co-exist, such as a strong selective disadvantage in resistant bacteria (for instance lower growth rate), and/or the existence of a spatial refuge (or density refuge) where (below which) the phage is not able to infect the bacteria. Chaudhry et al. (2018) offered an interesting explanation of how phages can persist in populations dominated by resistant bacteria, suggesting that the latter could produce susceptible bacteria at frequencies that would allow phage replication. Interestingly, this phenomenon has been suggested before (Bastías et al., 2010). The generation of phage-resistant strains in phage production systems could be a cause of concern and should therefore be included in the development of new methods, to minimize this possibility. Several authors have suggested that this problem can be avoided with the phage production setup, which will be discussed in the next section.
Another interesting study is that of Krysiak-Baltyn et al. (2018), which also incorporates variable infection parameters as a function of bacterial growth rate, and estimates operational cost and productivity in a simulated two-stage process system. One of the important conclusions of this theoretical study is that the optimal substrate concentration for bacterial growth should not be necessarily the same for bacteriophage production, and according their analysis the cost per mL of phages at a concentration of 4 × 1010 phages/mL could be as low as $ 1.78 × 10–2. It would be interesting to have an experimental validation of this estimation, and to determine how it is adapted to different economies or countries.
Finally, bacteriophage evolution must be considered in a production process as well, since the phages might increase or decrease their efficiency to infect bacteria over time (Lenski and Levin, 1985). This concept could be represented as infection rates in host-range experiments, where even methods for host-range expansion can be achieved for phage therapy applications (Mapes et al., 2016). This situation has been simulated in batch cultures, showing that appearance of phage mutants is strongly dependent on the genetic flexibility of phages (rates of mutations) (Levin and Bull, 2004). The ability to predict phage evolution during production would be helpful to set up a production process minimizing the probability of altering the lytic properties of phages. The reviewed articles show that bacteriophage production models are an important approach that can help to find the best strategies, however, they need to be validated experimentally.
Experimental Experiences in Bacteriophage Production
There are several practical studies involving phage production. Some are focused on phage production in bioreactors, while others focus on the evaluation and optimization of the process. As expected, these experiences also consider a step of bacterial growth and phage infection/propagation in flasks and bioreactors (Table 2). This data is useful to give insights on how certain models of host-bacteriophage can be used for propagation and increase levels of phage production. The most common host-phage systems used are E. coli strains and their phages, probably due to the amount of information regarding these bacteria-phage systems (E. coli phages T3, T4, and T7) and the lack of information on other bacteria-phage systems.
According to one reported case, the titers obtained can be as high as 1.2 × 1016 PFU mL–1 in a batch bioreactor (5 L) (Sochocka et al., 2015). This level of production agrees with the production needed for therapeutic purposes (>1 1010 PFU mL–1), considering purification steps, the decay rate of phages and stability or shelf life (Naghizadeh et al., 2018). Other authors have also reported promising levels of production of 5 × 1012 PFU mL–1 in 1.2 L (Warner et al., 2014), and 2.4 × 1013 PFU day–1 in 1 L (Nabergoj et al., 2018a; Table 2).
It is difficult to establish comparisons about which method could be more efficient, since they use different culture procedures and different host-bacteriophage systems. Batch culture is the cheapest (not simplest) way to produce bacteriophages, but it is highly limited by the maximum volume of the equipment available, total operation times and substrate availability (higher concentrations can be inhibitory for bacterial growth). Continuous culture has higher scalability when optimizing bacterial dilution rate via inlet and outlet flux modification. Besides, regulating the dilution rate will allow direct control over the bacterial growth rate, which has a direct influence on infection parameters like burst size, adsorption constant and latent period (Mancuso et al., 2018; Nabergoj et al., 2018b). Dilution rate can also be used to increase the productivity of the system, as was shown by Nabergoj et al. (2018a) where a maximum phage productivity of 109 phages mL–1 h–1 was achieved with a low dilution rate of 2 h–1 in a 1 L cellstat system. A continuous operating system can be permanently operative, and is therefore the most convenient way to produce a biotechnological product for a company. However, they are hard and expensive to implement, requiring constant monitoring to maintain the steady state. A totally continuous process for bacterial growth and bacteriophage production could increase the probability of bacteriophage resistance occurrence if specific countermeasures are not adopted (Middelboe et al., 2001).
Some authors have suggested the implementation of two-stage processes, one exclusively for bacteria production and a second for propagation of phages (Schwienhorst et al., 1996; Sauvageau and Cooper, 2010; Nabergoj et al., 2018a). This can be achieved with a cellstat system, where two bioreactors are connected in series with a constant flow through the system. In this case, the flow rate between the reactors and the volume in each reactor (and dilution rate and bacterial growth rate by addition) can be controlled in order to reach maximum productivity (Nabergoj et al., 2018a). Another interesting setup proposed by Sauvageau and Cooper (2010) consists of a semi-continuous system of a two stage, self-cycling process. In this case, each stage functions similarly to a batch culture, where the bacteria are first grown separately from the phage, and then introduced to the phage propagation stage when an appropriate concentration is reached, thus allowing the initiation of the infection process using a desired multiplicity of infection (Kasman et al., 2002). This setup also has the advantage of not requiring permanent monitoring to maintain the steady state of continuous systems, and has been used to obtain productivity of 7.59 × 1014 PFU mol CO2–1 (Sauvageau and Cooper, 2010). Both examples, the cellstat system and the two stage self-cycling process, have the great advantage that bacteria are grown in the absence of phages, therefore bacteriophage resistance is not favored during the process.
Finally, is important to note that there are some parameters that are not always reported in studies about bacteriophage production. For instance, parameters such as aeration proportion or air inflow into the bioreactor are mentioned only in two reports (Sauvageau and Cooper, 2010; Santos et al., 2014), even though this is one of the most important parameters in the production of bacteria at the industrial level. Information about other parameters such as energy transfer, different substrate utilization, bioreactor design, agitation, propellers and materials of construction in bacteriophage production is scarce or inexistent.
The rediscovery of the potential use of phages in a broad spectrum of applications is very exciting and promising. The evidence suggests that systems for bacteriophage production that reduce the probability of bacteriophage resistance occurring should be preferred, such as a cellstat or a two-stage self-cycling process. These options would also allow controlling variables to increase the productivity of the process. Nevertheless, models for production of bacteriophages are far from being established and can be improved in several ways. There are still many challenges to overcome. Further studies on optimized large-scale bacteriophage production, infrastructure and equipment costs, different safety concerns and application dosage are needed, and experience suggests that these challenges should be faced with collaborative efforts of academic and industrial partners.
Finally, it is important to note that most of the models for bacteriophage production can be applied within a specific range of values for parameters of phage infection and bacterial growth. Therefore, regardless of the important advances in phage production models and setups, the deep knowledge of the specific phage-bacteria system will be always the first requirement in order to establish an efficient phage production system.
RG, SL, and RB conceived the work and wrote the manuscript. KG, GH, and JR wrote the sections of the manuscript. All authors contributed to the bibliographic revision, manuscript revision, read, and approved the submitted version.
This work was funded by the CONICYT PFCHA DOCTORADO/2016 21161133 and Postdoctorado PUCV 2018.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
RG thanks CONICYT PFCHA DOCTORADO/2016 21161133. SL thanks Postdoctorado PUCV 2018.
Bastías, R., Higuera, G., Sierralta, W., and Espejo, R. T. (2010). A new group of cosmopolitan bacteriophages induce a carrier state in the pandemic strain of Vibrio parahaemolyticus. Environ. Microbiol. 12, 990–1000. doi: 10.1111/j.1462-2920.2010.02143.x
Bull, J. J., Millstein, J., Orcutt, J., and Wichman, H. A. (2006). Evolutionary feedback mediated through population density, illustrated with viruses in chemostats. Am. Nat. 167, E39–E51. doi: 10.1086/499374
Chaudhry, W. N., Pleška, M., Shah, N. N., Weiss, H., McCall, I. C., Meyer, J. R., et al. (2018). Leaky resistance and the conditions for the existence of lytic bacteriophage. PLoS Biol. 16:e2005971. doi: 10.1371/journal.pbio.2005971
Doss, J., Culbertson, K., Hahn, D., Camacho, J., and Barekzi, N. (2017). A review of phage therapy against bacterial pathogens of aquatic and terrestrial organisms. Viruses 9:50. doi: 10.3390/v9030050
González-Menéndez, E., Arroyo-López, F. N., Martínez, B., García, P., Garrido-Fernández, A., and Rodríguez, A. (2018). Optimizing propagation of Staphylococcus aureus infecting bacteriophage vB_sauM-phiIPLA-RODI on Staphylococcus xylosus using response surface methodology. Viruses 10:153. doi: 10.3390/v10040153
Kasman, L. M., Kasman, A., Westwater, C., Dolan, J., Schmidt, M. G., and Norris, J. S. (2002). Overcoming the phage replication threshold: a mathematical model with implications for phage therapy. J. Virol. 76, 5557–5564. doi: 10.1128/JVI.76.11.5557-5564.2002
Krysiak-Baltyn, K., Martin, G. J. O., Stickland, A. D., Scales, P. J., and Gras, S. L. (2016). Computational models of populations of bacteria and lytic phage. Crit. Rev. Microbiol. 42, 942–968. doi: 10.3109/1040841X.2015.1114466
Lenski, R. E., and Levin, B. R. (1985). Constraints on the coevolution of bacteria and virulent phage: a model, some experiments, and predictions for natural communities. Am. Nat. 125, 585–602. doi: 10.1086/284364
Levin, B. R., Moineau, S., Bushman, M., and Barrangou, R. (2013). The population and evolutionary dynamics of phage and bacteria with CRISPR-mediated immunity. PLoS Genet. 9:e1003312. doi: 10.1371/journal.pgen.1003312
Łusiak-Szelachowska, M., Weber-Da̧browska, B., Jończyk-Matysiak, E., Wojciechowska, R., and Górski, A. (2017). Bacteriophages in the gastrointestinal tract and their implications. Gut Pathog. 9:44. doi: 10.1186/s13099-017-0196-7
Mancuso, F., Shi, J., and Malik, D. J. (2018). High throughput manufacturing of bacteriophages using continuous stirred tank bioreactors connected in series to ensure optimum host bacteria physiology for phage production. Viruses 10:537. doi: 10.3390/v10100537
Mapes, A. C., Trautner, B. W., Liao, K. S., and Ramig, R. F. (2016). Development of expanded host range phage active on biofilms of multi-drug resistant Pseudomonas aeruginosa. Bacteriophage 6:e1096995. doi: 10.1080/21597081.2015.1096995
Middelboe, M., Hagström, A., Blackburn, N., Sinn, B., Fischer, U., Borch, N. H., et al. (2001). Effects of bacteriophages on the population dynamics of four strains of pelagic marine bacteria. Microb. Ecol. 42, 395–406. doi: 10.1007/s00248-001-0012-1
Międzybrodzki, R., Borysowski, J., Weber-Da̧browska, B., Fortuna, W., Letkiewicz, S., Szufnarowski, K., et al. (2012). Clinical aspects of phage therapy. Adv. Virus Res. 83, 73–121. doi: 10.1016/B978-0-12-394438-2.00003-7
Nabergoj, D., Kuzmić, N., Drakslar, B., and Podgornik, A. (2018a). Effect of dilution rate on productivity of continuous bacteriophage production in cellstat. Appl. Microbiol. Biotechnol. 102, 3649–3661. doi: 10.1007/s00253-018-8893-9
Naghizadeh, M., Karimi Torshizi, M. A., Rahimi, S., and Dalgaard, T. S. (2018). Synergistic effect of phage therapy using a cocktail rather than a single phage in the control of severe colibacillosis in quails. Poult. Sci. 98, 653–663. doi: 10.3382/ps/pey414
Santos, S. B., Carvalho, C., Azeredo, J., and Ferreira, E. C. (2014). Population dynamics of a Salmonella lytic phage and its host: implications of the host bacterial growth rate in modelling. PLoS One 9:e102507. doi: 10.1371/journal.pone.0102507
Smrekar, F., Ciringer, M., Jančar, J., Raspor, P., Štrancar, A., and Podgornik, A. (2011). Optimization of lytic phage manufacturing in bioreactor using monolithic supports. J. Sep. Sci. 34, 2152–2158. doi: 10.1002/jssc.201100182
Sochocka, M., Tomczyk, T., Sobczyński, M., Szermer-Olearnik, B., and Boratyński, J. (2015). The kinetics of Escherichia coli B growth and bacteriophage T4 multiplication in SM-1 novel minimal culture medium. J. Gen. Appl. Microbiol. 61, 75–81. doi: 10.2323/jgam.61.75
Tokman, J. I., Kent, D. J., Wiedmann, M., and Denes, T. (2016). Temperature significantly affects the plaquing and adsorption efficiencies of listeria phages. Front. Microbiol. 7:631. doi: 10.3389/fmicb.2016.00631
Keywords: bacteriophage, bacteriophage production, bacteriophage therapy, phage therapy, phage production models
Citation: García R, Latz S, Romero J, Higuera G, García K and Bastías R (2019) Bacteriophage Production Models: An Overview. Front. Microbiol. 10:1187. doi: 10.3389/fmicb.2019.01187
Received: 16 January 2019; Accepted: 09 May 2019;
Published: 04 June 2019.
Edited by:Robert Czajkowski, University of Gdańsk, Poland
Reviewed by:Konstantin Anatolievich Miroshnikov, Institute of Bioorganic Chemistry (RAS), Russia
Diana Gutiérrez, Institute of Dairy Products of Asturias (IPLA), Spain
Malgorzata Barbara Lobocka, Institute of Biochemistry and Biophysics (PAN), Poland
Copyright © 2019 García, Latz, Romero, Higuera, García and Bastías. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Roberto Bastías, firstname.lastname@example.org
†These authors have contributed equally to this work