- 1Centro de Investigación en Estudios Avanzados del Maule (CIEAM), Vicerrectoría de Investigación y Postgrado, Universidad Católica del Maule, Talca, Chile
- 2Laboratorio de Microbiología Costera, Facultad de Ciencias del Mar y Recursos Biológicos, Universidad de Antofagasta, Antofagasta, Chile
- 3Departamento de Matemáticas, Facultad de Ciencias Básicas, Universidad de Antofagasta, Antofagasta, Chile
- 4Centro de Biotecnología de los Recursos Naturales, Facultad de Agronomía y Ciencias Forestales, Universidad Católica del Maule, Talca, Chile
- 5Departamento de Matemática, Física y Estadística, Facultad de Ciencias Básicas, Universidad Católica del Maule, Talca, Chile
- 6Instituto de Ciencias Naturales Alexander von Humboldt, Facultad de Ciencias del Mar y Recursos Biológicos, Universidad de Antofagasta, Antofagasta, Chile
A mathematical first-order difference equation was designed to predict the dynamics of the phage-bacterium adsorption process in aquatic environments, under laboratory conditions. Our model requires knowledge of bacteria and bacteriophage concentrations and the measurements of bacterial size and velocity to predict both the number of bacteriophages adsorbed onto their bacterial host and the number of infected bacteria in a given specific time. It does not require data from previously performed adhesion experiments. The predictions generated by our model were validated in laboratory. Our model was initially conceived as an estimator for the effectiveness of the inoculation of phages as antibacterial therapy for aquaculture, is also suitable for a wide range of potential applications.
Introduction
Viruses infecting bacteria, known as bacteriophages or phages, are the most common and diverse entities in the biosphere, especially in aquatic environments (Breitbart and Rohwer, 2005; Suttle, 2005). In the ocean, its numbers are generally near 10 millions of particles per mL of water, a value that decreases with water depth, coastal distance depth, and coastal distance (Cochlan et al., 1993; Paul et al., 1993). Its major abundance has been observed in coastal sediments where their numbers can reach up to one billion per mL (Bergh et al., 1989; Paul et al., 1993; Hewson et al., 2001). Besides this, they are important for mobilizing and transferring genetic information (Abedon, 2008), and promoting the evolution of organisms. Phage-host interactions are models of synergic evolution with a vast population size and a short duplication time. This fact makes these systems ideal for the study of evolution and community-dynamics and for the quick observation of the ecological principles (Lenski and Levin, 1985; Weitz et al., 2005; Brockhurst et al., 2007).
Félix d’Herelle discovered bacteriophages in 1917 and immediately visualized their potential to control bacterial diseases (D’Herelle, 1917). Nowadays, due to the excessive use of antibiotics (Blazquez et al., 2002; Harbottle et al., 2006; Madhusudana Rao and Lalitha, 2015), the strategy of using lytic bacteriophages as antimicrobial agents, has become an interesting alternative for controlling bacterial infections (Matsuzaki et al., 2005; Almeida et al., 2009; Carvalho et al., 2012; Verbeken et al., 2014; Keen and Adhya, 2015; Kutter et al., 2015; Sarhan and Azzazy, 2015). The main advantages of bacteriophage usage are their high specificity (even at strain level), infective ability, and their exponential growth curve in the presence of a target bacteria, with a higher efficiency than antibiotics (Levin and Bull, 1996). According to Drake et al. the mutation rates of bacteriophages are of three orders of magnitude higher than bacteria (Drake et al., 1998), which means, they have the potential to overcome the bacterial resistance mechanism (Buckling and Rainey, 2002; Brockhurst et al., 2007). However, there are diverse strategies of bacteriophage infection and bacterial resistance, allowing for their natural coexistence and equilibrium both at ecological (Schrag and Mittler, 1996; Rabinovitch et al., 2003; Heilmann et al., 2012) and genomic levels (Buckling and Rainey, 2002; Barrangou et al., 2007; Stern and Sorek, 2011; Samson et al., 2013). However, if we want to apply methodologies based on bacteriophages as biological controllers, it is important to understand the dynamics that guide the phage-host interaction.
The law of mass action, particle collision, and diffusion of the virion are the determinant factors for the collision and adsorption of bacteriophages onto bacterial surfaces (Koch, 1960; Stent, 1963; Payne and Jansen, 2000; Buckling and Rainey, 2002; Moldovan et al., 2007; Cairns et al., 2009) are used in most mathematical approaches. The more accepted models define the initial adsorption as an irreversible process, which depend on the time (Gani, 1965; Chang and Chang, 1969; Murray and Jackson, 1992; Carletti, 2007; Stopar, 2008; Cairns et al., 2009) and are based on the Poisson distribution (Stent, 1963; Stopar, 2008). The determining factors for these models are bacteriophage and bacterial density (Beretta and Kuang, 1998; Payne and Jansen, 2001; Cairns et al., 2009), temperature (Moldovan et al., 2007), and the presence of cellular receptors (Schwartz, 1976; Moldovan et al., 2007) among others. However, some of these models minimize the surrounding space so that the bacterium is unable to reach and influence due to its swimming movements and the size of the area (Yassky, 1962; Gani, 1965; Lowe et al., 1987; Moldovan et al., 2007). Finally all of these deviations were commonly fixed through the generation of a phage adsorption rate constant (k) a posteriori, performed in a set of in-vitro experiments (Hyman and Abedon, 2009). On the other hand, the probability of an encounter between a bacteriophage and its host has been attributed to factors such as abundance, movement, shape, and size (Murray and Jackson, 1992; Payne and Jansen, 2001). The size and number of the host defines the quantity of the occupied space in a determined volume, which logically establishes the probability of a bacteriophage-host encounter. Besides, if one considers, that some bacteria possess high motility in comparison to the diffusion coefficients of virion particles, it is imperative to include this factor in the calculation of the initial adsorption. It is estimated that in seawater, it is possible to find a mean value of 1 × 106 bacteria per mL (Cochlan et al., 1993). This means that the available space for each one of them is around 0.1 μL (1 × 10–6 mL).
Previous research has show that bacteriophages are present in an average concentration of one magnitude order over the total amount of bacteria (Wommack and Colwell, 2000), and up to three magnitudes of order as the maximum in highly populated areas (Bergh et al., 1989; Paul et al., 1993; Hewson et al., 2001). This means that a planktonic bacterium could theoretically find one phage virion every 10 μm in its traveling path when moving in any direction. Some marine bacteria can exceed the speed of 300 μm/s (Mitchell et al., 1995), in that case it would find a total of 100,000 bacteriophages per hour. On the other hand, the movement of a bacteriophage only depends on its diffusion coefficient that, in general, has a value close to 5 × 10–8 cm–2s–1 (Murray and Jackson, 1992), which, in most cases, is negligible compared to the speed of motile bacteria (Joiner et al., 2019).
In this work, we have designed and evaluated a model based on the interaction dynamics involving the strictly lytic bacteriophage vB_VorS-PVo5 (PvB) and its host bacteria Vibrio ordalii ATCC-33509 (VO). VO is an important pathogen that causes great losses for the aquaculture industry (Schiewe et al., 1981). Their specific phage (PvB) is important because of its potential to be used as an antibacterial agent in phage therapy in aquaculture (Echeverría-Vega et al., 2016). It was isolated from mollusk inhabitants of the Antofagasta coast (Chile). Our model was developed after analyzing the variation and recalculation of the probability of encounters between phages and bacteria in discrete time periods. It is possible to model the adsorption process by using the measurements of the initial parameters such as the total amount of initial bacteriophages, velocity, size, and number of bacteria and thus, we can predict both the number of adsorbed phages and infected bacteria at the same time. Unlike the classical adsorption theory (Hyman and Abedon, 2009) which determinates an adsorption constant a posteriori through a battery of previous tests, our model uses only the measured physical parameters of the particles involved (phages and bacteria) to model the adsorption a priori. As a result, it also reduces the risk of experimental errors. The proposed model fits with the data measured and its construction is based on the probability of a collision between particles, therefore it could be applied to other phage-bacterium systems and can be extrapolated to predict any other phenomena involving particle collisions.
Materials and Methods
Bacterial Strain and Growth Conditions
A Chilean wild-type of the flagellated Gram-negative Vibrio ordalii (VO) was used for this study. VO is a marine pathogenic bacterium originally isolated from infected Argopecten purpuratus (Lamarck 1819) larvae (Riquelme et al., 1995) and have a high impact on the salmon farming industry. VO was routinely grown and maintained in a liquid tryptic soy broth (TSB) and trypticase soy agar (TSA) supplemented with a 50 % total volume of sterile seawater at 25 °C (Oxoid, United States).
Culture of Phage and Bacteria
The phage vB_VorS-PVo5 (PvB) is described as a lytic strain belonging to the Siphoviridae family (Echeverría-Vega et al., 2016). It was obtained from soft tissues of Perumytilus purpuratus (Lamarck 1819) recovered from the coast of Antofagasta, Chile using a previously described methodology (Echeverría-Vega et al., 2019).
Transmission Electron Microscope Analysis (TEM)
All the preparations was performed following the protocol described by Carlson (Carlson, 2005) with some slights modifications. Stabilized carbon formvar grids (EMS, United States) were used to mount the phage samples for TEM. A drop (10 μL) of the sample containing a viral concentration greater than 1 × 108 viral particles per milliliter (VP/mL) was placed on the grid. After 15 min the liquid was removed with an absorbent paper. Subsequently, 10 μL of a staining reagent solution (1 % Uranyl Acetate in water) was placed on the grid for 15 min. Excess water was removed with blotting paper and then, allowed to air dry for 15 min. The samples were observed using a Phillips 100 CM TEM.
Setting of Initial Parameters
Counting the bacterial number and size was completed by the DAPI-staining method (Porter and Feig, 1980). For this, 400 cells for a sample were counted in triplicate. The bacterial size was obtained from the mean of several measurements of 100 cells, using images obtained from the Zeiss Axiolab fluorescence microscope (Zeiss, Germany). The bacterial volume (VB) was calculated considering the cells as cylindrical by using the formula V = πr2l, where r is the radius and l is the length. The cell including the flagellum has a length of 13 μm (Schiewe et al., 1981) and Vibrio oscillates with a wavelength of 3,2 μm and can turn in angles up to 90° when tumbling (Xie et al., 2011). Bacteria can tumble every 1 s in average and furthermore generate turbulences in the medium when rotating (Joiner et al., 2019). Then, bacteria are likely to interact with phages over a space beyond their place of positioning. All these factors were taken into account while defining the “Vital space” (VS) which is the entire potential coverage zone where the bacterium can interact with the phage or the volume that the bacterium is capable of reaching when it is in a determined space in a specific period of time. Our hypothesis is that VS is one spherical volume with a radius of 5 times the size of ΔlB surrounding each bacterium. Bacteriophage number counting was achieved using the method of double-layer agar (Carlson, 2005; Kropinski et al., 2009). For that purpose, a triplicate of four dilutions were performed in TSB medium (1/10, 1/100, 1/1000, and 1/10000) containing the bacteriophage. The plates were incubated overnight to perform the counting of plaques by visual inspection the next day. The measurements of the bacteriophage size were completed by transmission electron microscopy (TEM) images. The bacteriophage volume (PB) was calculated considering the head as a sphere and the tail as a cylinder. For the measurements of bacterial movement, one sample was taken from the liquid culture in the exponential phase (16 h after the renovation) and then it was observed under the Zeiss Axiolab microscope (Zeiss, Germany) equipped with an objective grid with 100 squares of 100 μm2 each. One hundred measurements of the time taken by one bacterium to cross a zone of 20 μm was carried out to calculate bacterial velocity. For the bacteriophage movement, a diffusion constant of 5 × 10–8 cm–2s–1 was considered (Murray and Jackson, 1992).
Adsorption Curves and Determination of Growth Parameters
A bacterial growth curve was generated with the data obtained after measuring the turbidity of the liquid cultures by UV-visible spectrophotometry in triplicate. To determine the effect of the phage on bacterial growth, VO was inoculated in 500 mL of TSB medium with a final concentration of 5 × 106 colony-forming units per milliliter (CFU/mL) and distributed in 72 glass tubes. An inoculum of 5 × 107 plaque-forming units per milliliter (PFU/mL) of bacteriophage was added to 36 tubes to obtain a multiplicity of infection (MOI) = 10. During the first hour of cultivation, the absorbance was measured every 15 min. Then, from the second hour to the fifth and a half hour, it was measured every 30 min. A final measurement was taken after 24 h. To determine the latency period and the burst-size of the phage under study, the growth curve was achieved using the one-step methodology (Carlson, 2005; Hyman and Abedon, 2009). For this, 1 × 107 CFU/mL of VO with a total of 1 × 102 PFU/mL of PvB was used as the inoculum. A 100 μL sample was collected after every 15 min to grow and later count the phage plaques in double agar plates. The latency time is defined as the period in which there is no significant variations in the number of phages (slope = 0). The burst-size was calculated by dividing the number of total phages obtained in the beginning of the logarithmic phase with the initially added phages. In order to empirically determine the rate of the adsorption of a phage, the previously described established protocols (Hyman and Abedon, 2009; Kropinski, 2009) were performed with some modifications. A VO liquid culture was grown at 25°C in TSB medium supplemented with 50 % sterile marine water, an inoculum of 1 × 108 CFU/mL was added after 24 h in the exponential phase. Triplicate assays using two different dilutions (1 × 106 and 1 × 105 CFU/mL) of this bacterial culture in the exponential phase were performed in 50 ml plastic tubes. Two different concentrations of bacteriophages (8000 and 800 PFU/mL) were added. Then, we took 1 mL samples once every minute for a total time of 10 min and transferred them to a sterile tube containing chloroform to lyse the bacteria. Samples were agitated for 10 s in a vortex, and the supernatant was immediately extracted to count the free phage virions according to the methodology described above.
Modeling of the Adsorption Process
The simulation and graphics were completed using the Maple software (Maplesoft, Cybernet Systems Co., Ltd., Japan) and Microsoft Excel software (Microsoft Office, United States). The simulations consider:
– Infected bacteria, increasing the number of bacteria and phages keeping the MOI = 1.
– Infected bacteria over time, maintaining the bacterial number at 106 CFU and the MOI > 1.
– Infected bacteria by varying its speed, maintaining the bacterial concentration in 1 × 106 CFU and MOI = 1.
– Adsorbed phages, maintaining constant the number of phages at 1 × 106 UFP/mL and by varying the bacterial number.
We performed adsorption experiments using the previously described protocol and default parameters with both 8000 and 800 phages, respectively. We plotted the adsorption curves comparing with both the classical adsorption model (Hyman and Abedon, 2009) and a cubic polynomial regression model which best fit the experimental data. We evaluated the predictability of our model comparing the experimental curves with those generated by simulations. This model was deposited into BioModels (Chelliah et al., 2014) and assigned the identifier MODEL2007190001.
Results
The phage PvB has been described in previous studies (Echeverría-Vega et al., 2016) as a strain belonging to the Siphoviridae family (Order Caudovirales), which groups viruses containing double-stranded DNA chains. It produced lysis plaques of ∼ 2 mm of diameter in the agar plate when infected with VO. The observation of the viral particles in electronic transmission microscopy revealed an icosahedral capsid of ∼ 85 nm of diameter and a tale of ∼ 150 nm long, which gives a volume of 2,1 × 10–17 mL (Figure 1).
The bacterial length for VO is approximately 1.8 μm with 0.6 μm of width and its calculated volume is 5 × 10-13 mL (VB) which is equivalent to a sphere of 1 uL diameter (ΔlB). Vital space VS was calculated as 5 × 10-10 mL. The bacterial mean velocity was 30 μm/s (SE = 0.49). The growth test showed a strong diminishing of VO growth when PvB was inoculated to the culture, which indicated a clear inhibitory effect of the phage on the bacterial population after 120 min, VO lysis still increased until 1,450 min (Figure 2).
The corresponding measurements were taken on the growth curve of PvB by the one-step method. The lapsed time between the adsorption and the bacterial lysis (latency time) was estimated as 100 min and the burst-size was around 152 phages/infected-bacteria. These parameters were considered in the typical range for a Siphoviridae family member (Yu et al., 2013) and allowed us to use it for further tests oriented to the bacterial control (Figure 3).
Our proposed model of adsorption considered the established theory according to which the encounters should randomly occur as per the Poisson equations (Ellis and Delbrück, 1939). If adsorption is considered as a dependent function of the space, in which phage and bacteria are found, the total volume of the liquid medium in which they cohabit, will be understood as an integration of spaces of the size of one bacterial vital space (VS). In this way, the probability of the encounter of a phage with the bacteria is computed according to the total number of spaces.
We define “Total Slot” (ET) as the theoretical number of bacterial vital spaces that fits in the total volume, in other words, the total volume used in the cell culture divided by the volume of one vital space (Equation 1).
ET: total slots.VT: total volume.VS: bacterial vital space.
For the case of VO, the VS is 5 × 10–10 mL, which means in 1 mL, there is around 2 × 109 ET available.
If an amount of phage (insignificant volume of 2x10–17 mL) is added (FT) to a known volume of a liquid medium (VT), it will be distributed using probability as defined according to the Poisson equation. Thus, the event X to find n phages in one slot E is defined by the probability:
where λ is the average of phages in a slot (E).
If n=0 and λ = FT(1/ET), then P0=P(X = 0) is the probability of not finding phages in one slot and it is defined by:
Then, the probability of finding at least one phage in one slot (or simply adsorption probability) is:
P0: probability that no phage is in one place.
Pa: probability that at least one phage fits in a punctual slot = (1−P0).
FT: total phage added initially.
ET: total slots.
Bacterial velocity is a critical parameter in this process. If we know the bacterial velocity (υB), we can define the interval of time used by one bacterium to move a segment of length ΔlB, equivalent to the diameter of one sphere of the same volume of the bacterium and it is defined as:
The averages for VO are: ΔlB = 1 μm and υB = 30 μm/s. Therefore Δt = 1/30 s and the time taken by VO to move from E to E′ (or E(t)) to Et +Δt) is 0.033s approximately.
Considering initial conditions in t0, we can define the steps of our discrete model as follows:
where k = 1, 2, 3,…, and therefore we define the variables as:
Etk: number of slots with phages at time tk.
Ftk: number of free phages at time tk.
Btk: number of bacteria with adsorbed phages measured at time tk.
Using these variables, we can write the model as:
Equations 7, 8, and 9 show that at the initial time t0, adsorption probability determines the number of slots occupied with phages. Its initial number and the proportion of empty slots (without bacteria) determine the number of virions (free phages). On the other hand, the number of bacteria multiplied by the probability of adsorption generates the number of infected bacteria (with at least one phage adsorbed).
To explain the model in a very simplistic way, in Figure 4 we can see a total volume VT divided in cubes representing the total space (ET = 300). If we add FT = 100 phages, the probability that one slot contains phages (Pa) is 28% and its distribution is shown in Figure 4A. The probability of a phage not fitting in an empty slot (P0) is 72% (Eq. 4). If BT = 10, the probability that a bacteria fits in the infected place is given by P. We get the initial number of infected bacteria Bt0 = 3 by multiplying Pa with the bacterial number (Eq. 9). We also get free phages in t0 (Ft0 = 67) from Eq. 8 (Fig. 4B). For each time interval “Δt”, the bacteria will move 1ΔlB. Figure 4C shows 10 time intervals. This new positioning of bacteria generates a new adsorption probability (Figure 4D), that is:
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/550979/fmicb-11-550979-HTML/image_m/fmicb-11-550979-g004.jpg)
Figure 4. Schematic 3D representation of the phage-bacteria adsorption process. A grid of 100x100x3 squares represent the total space (ET). The sides of squares are 10ΔlB. One vital space (VS) is the spherical volume that fits inside one cube (represented in A). Red dots on yellow squares represent slots with phages. (A) Total phages (FT = 100) are distributed randomly across the space (sizes are not in scale), green ellipses are uninfected bacteria with their vital space represented by a green mesh sphereVS. (B) When one bacterium arrives in a yellow slot (with at least one phage inside), it is infected and turns into red. The phages in these vital spaces are adsorbed. (C) Bacteria move 10ΔlB in random directions, which is indicated by black arrows. (D) Previously infected bacteria (red colored) are able to gain more adsorbed phages when moving to a new infected slot. All the uninfected bacteria that move to an infected slot become infected. A dark halo surrounding bacteria and phages, indicates they are located on the lower level and a white halo indicates they are on the upper level.
and
Once one bacterium moves from one slot to another, it carries the adsorbed phages along with it and thus leaves the slot empty where the phages were previously. This diminishes the probability that one bacterium is infected the next time (tk+1). Besides, it is feasible that this bacterium would continue adsorbing phages when it moves. Then, the number of infected bacteria in each interval depends on the number of slots containing phages in each time-lapse (Δt). From this, the followings points can be inferred:
– The total infected bacteria in a time lapse Δt will be the sum of the bacteria that was placed in an infected slot in time tk–1, and
– The total amount of infected bacteria in tk will be the total of infected bacteria in time tk–1 plus the ones that were infected between time tk–1 and tk. This new number of infected bacteria is calculated by multiplying the number of not infected (free) bacteria with the ratio between total slots with phages in tk (Etk) and total slots (ET).
– The total amount of free bacteria is calculated as the difference between total bacteria and those infected in tk–1.
This can be summarized in the following equation:
This shows that there is a clear dependency on infected bacteria within the space. The total slots with phages which remain for the next time-interval (Etk) are the slots with phages in tk−1, minus the probability that one bacterium moves to a slot with phages (Eq. 7). No limit in the number of adsorbed phages on each bacteria was considered. Therefore we have the following equations:
And for phages:
where correspond to the per-capita availability of slots E and free phages F in a period of time Δt (i.e., for N=E or N=F). From equations 7, 8, 9, 12, 13, and 14 we can see the following equations system:
Furthermore, from equation 13 and 14 we calculated their explicit solutions:
and
By combining, we can define slots in term of phages:
With that, Eq. System 1 can be reduced to a first order system with two states (B and F) given by:
Phage-Bacteria System Analysis
Adsorption Time
From Eq. 6, the system times are given by tk = t0 + kΔt, so it is possible to determine the time when there were no more free phages (or when they were all absorbed by bacteria), that is, we can calculate the value of K as the number of approximate time periods at which the system stabilizes. From the Eq. 16, we look for Ftk < 1, and thus we obtain the inequality:
and if:
then, K is the lowest integer greater than h. As K gives the number of time intervals that must occur until 0.99 phages remains, the number of phages will be approximated to the nearest greater integer, i.e., > 0.5, considered as one unit. Then the new value of K is obtained from Ftk < 1/2, resulting in:
Incorporating MOI
If the multiplicity of infection (MOI = FT/BT) is incorporated into the model, we get a reinterpretation of Eq. System 2, as:
Stationary Solutions
As Eq. system 3 is a two-dimensional discrete model, the study of its dynamics is subject to determining its equilibrium solutions or equilibrium points and its associated stability. We understand that for the equilibrium solution of any pair B∗, F∗ that Btk = B∗ and Ftk = F∗, for all K≥0 so that stationary and constant system solutions are available. Then from the Eq. System 3, we propose the following option:
and
after solving it we get
and
From the resolution of Eq 23 and Eq 24, the system has infinite stationary or equilibrium solutions of the form e = (B∗, 0), with B∗ ≤ BT, that is, solutions showing that all or part of the total bacteria could be infected after the kth time step in which the phages are depleted. For the stability analysis of the equilibria, we will use the classic theory of determination of the eigenvalues of the Jacobian matrix J, where J is obtained by the partial derivatives with respect to phages and bacteria on the right side of Eq. system 3, that is:
When evaluating the equilibria e = (B∗, 0) on J we get:
where it can be seen that J(e) is a triangular matrix, so its eigenvalues are on the diagonal as: λ1 = 1 y . Both of the eigenvalues have a positive real part, establishing that the obtained equilibria e will have unstable nodes, that is, the long-term dynamics of the system will tend to favor the total infection of bacteria.
Simulations and Real Data
We used the software Matlab (The MathWorks, R2016b) to simulate different scenarios for phage adsorption and compared them with the real data, considering the previously measured parameters for the phage PvB and VO.
Simulation Varying Parameters
Default parameters were VB = 5.24 × 10–10 mL; VT = 1 mL; ΔlB = 1 μm; υB = 30 μm/s; BT = 106 CFU; FT = 106 PFU. Using Eq. System 3, the simulations were performed considering default parameters and modifying BT, FT, MOI, and υB as shown in Figure 5.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/550979/fmicb-11-550979-HTML/image_m/fmicb-11-550979-g005.jpg)
Figure 5. The simulations generated from the Eq. System for Btk and Ftk. (A) Increasing the number of bacteria and phages without varying the MOI = 1. (B) Keeping the bacterial number constant at 106 CFU and varying the MOI. (C) Keeping BT = 106 CFU and MOI = 1 and varying υB. (D) Varying BT and keeping FT = 106 PFU/mL.
The simulations in Figure 5 were performed with Eq. system 3 and they showed the number of infected bacteria and the number of free phages with respect to time. Besides, the free phages reach zero in all simulations. In Figure 5A, the adsorption was shown to keep the MOI = 1 but changes both the total number of phages and bacteria present in the medium. It was observed that the adsorption curves depend on the amount/number of phages and bacteria present in the community. Therefore, as the total number of individuals present (phages and bacteria) increased, the adsorption slope increased. This establishes that as the dynamics are controlled by the time required for all phages to be adsorbed, then not all bacteria will become infected when MOI = 1. We performed a second simulation varying the MOI, but kept the number of bacteria constant at 106 CFU (Figure 5B).
The model predicted that an infection of the entire bacterial population would occur only with an MOI of ten or more. This is in agreement with our stability result that indicates that the total infection is reached in the equilibrium e = (BT, 0). On the other hand, increasing this value does not imply a significant increase in adsorption time and dynamics. The effect of the variation of the bacterial velocity was also simulated (Figure 5C), in which case, the adsorption time was reduced with higher bacterial velocity. Finally, the effect of reducing the number of bacteria present in the medium was simulated, keeping the number of phages constant at 106 PFU/mL (Figure 5D). In this case, it was observed that by decreasing the concentration of bacteria in the medium, it was not possible to reach the equilibrium e = (BT, 0). In summary, the MOI could explain the total number of infected bacteria, that is, only if MOI = > 10, all bacteria were infected, and the equilibrium e = (BT, 0) were achieved, which would be the theoretical prediction of stability.
Model Comparison With Real Data
To evaluate the predictability offered by the model, we compared in-vitro data from the initial adsorption of PvB phage on VO with simulations of the same, including, a cubic polynomial regression model of the data as an external modeling source. We considered the default parameters with two different initial phage concentrations, 800 and 8000 PFU.
The cubic adjustment made with the real data gave us a regression model (with ai determined parameters) that explains 90% of the variability of the data, measured with the coefficient of determination R2. The construction of the fit model provided us with a reliable tool for comparison with our proposed phage-bacteria absorption model (Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/550979/fmicb-11-550979-HTML/image_m/fmicb-11-550979-g006.jpg)
Figure 6. Simulation for adsorption dynamics of the phage PvB on Vibrio ordalii versus experimental data using default parameters with (A) 8000 phages and (B) 800 phages. Cubic adjustment (with R = 0.9560 for (A) and R = 0.9356 for (B) and classical model curves (using the methodology proposed by Hyman and Abedon, 2009) were included for comparison.
Discussion
Bacteriophages are the most abundant organisms in the biosphere (Clokie et al., 2011; Aziz et al., 2015). Phages cohabit with their host bacteria in an equilibrium state, modulated by both phage infection rate and bacterial resistance mechanisms (Bohannan and Lenski, 2000; Wei et al., 2010; Bull et al., 2014). This is comparable to predator-prey dynamics, but in this case the encounters are governed by probability. In the present work, we built a mathematical model using the major factors determining the initial encounter between phages and bacteria in the function of time. This was done to predict the number of infected bacteria and remnant free phages. At first, we included the amounts of bacteria and phages as determining factors in concordance with previous reports (Yassky, 1962; Murray and Jackson, 1992; Moldovan et al., 2007). However, taking into account the environmental dynamics, we included the living space as a function of cell size and the bacterial velocity, in a probabilistic function that allowed for predicting such encounters more accurately. The data required to feed the model was easily obtainable from simple laboratory experiments (Echeverría-Vega et al., 2019) or public databases, (i.e., BioNumbers database) (Milo et al., 2010). The model assumes that the phage movement occurs by simple diffusion (Yassky, 1962; Gani, 1965; Schwartz, 1976) at approximately 5 × 108 cm2s–1 (Murray and Jackson, 1992), which is negligible compared to the VO velocity of 30 μm s-1. Additionally, it considers the bacterial size and its vital space or slot, which we defined as all the volume that the bacterium is capable of reaching when it is in a determined space in a specific period of time. Nevertheless, we did not consider the phage size, since it is three orders of magnitude lower, so it does not affect the vital space of the bacteria in normal conditions. In this way, it is possible to generate infection curves for phages and bacteria and therefore predict the number of adsorbed phages and infected bacteria over time. Some existing models consider phage infection as a simplified process based on the mass-action kinetics theory by means of an empirically determined factor called adsorption constant k (Krueger, 1931; Stent, 1963) without taking into account temporal dynamics. A novel model proposed by Joiner et al. (2019) considers bacterial movement and medium effect in the collision frequency. However, our design differs from the previously mentioned research as our model is capable of predicting the sustained decrease of free phages and the number of infected bacteria over time considering temporary changes in the probabilistic functions that modulate the curves. The results obtained from our model show that the principal factors in the initial adsorption in the phage-bacteria binary interaction, are the number and speed of host bacteria. Other environmental factors, like community complexity and habitat components were not included in our model.
Figure 5B shows that an MOI < 10 is not sufficient for the infection of all the bacteria present due to the occurrence of multiple infections, which eliminates the phages from the system before the time required to infect all bacteria has been reached. Nonetheless, one can see in Figure 5A that the concentration of bacteria and phages, while maintaining a constant MOI, causes a considerable effect on the curves, which means the MOI must not be taken as the only factor that modulates phage infections. On the other hand, the increase in speed and/or bacterial concentration significantly decreases the infection time (Figures 5C,D) which means that when the bacteria are highly active they are more likely to be infected by phages in less time (a kill the winner strategy), favoring the establishment of equilibrium states in complex environmental systems (Thingstad, 2000). The experimental validations performed based on the temporal dynamic model were shown to be highly concordant with those predicted. The importance of the present study lies in the generation of a model capable of describing the dynamics of the initial phage infection stage, which is crucial information for designing treatment strategies through phage therapy. Our work is based on statistical predictions using measurable physical parameters of phages and bacteria and hence does not require prior experimental adsorption tests like the classical model (Hyman and Abedon, 2009). Our model has the potential to be used as a basis for modeling both phage-bacteria interactions as well as other natural processes that involve random particle collisions.
Conclusion
Our model uses statistics to predict the dynamics of the phage-bacterium adsorption process in aquatic environments. It does not require data from previous adhesion experiments and can make predictions using simple measurements of speed, volume, and number of particles. The model is useful in predicting both the number of bacteriophages adsorbed on their bacterial host and the number of infected bacteria in a given specific time. The comparison with experimentally obtained curves demonstrates that the model is applicable in explaining these dynamics. The model was conceived to predict the effectiveness of inoculation in phage therapy, but is also suitable for other physical processes based on particle collision.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.
Author Contributions
AE-V and RA conceived the case study. PM-V, CS-S, and AE-V performed the experimental research. AE-V, MA, RD-V, and MS performed and tested the mathematical model. AE-V, FG, and RA checked and discussed the experiment results and wrote the manuscript. All the authors read and approved the final manuscript and critically reviewed the document. All authors contributed to the article and approved the submitted version.
Conflict of Interest
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.
Funding
This research work was carried out and financed under the following two projects:
1. Project CORFO FIC-R 13IDL2-18530 (RA).
2. Doctorate scholarship CONICYT FIC-R N°21092008 (AE-V).
References
Abedon, S. (2008). Bacteriophage Ecology. Population Growth, Evolution, and Impact of Bacterial Viruses. Advances in Molecular and Cellular Microbiology 15. New York, NY: Cambridge University Press. doi: 10.1017/CBO9780511541483
Almeida, A., Cunha, Á, Gomes, N. C., Alves, E., Costa, L., and Faustino, M. A. (2009). Phage therapy and photodynamic therapy: low environmental impact approaches to inactivate microorganisms in fish farming plants. Mar. Drugs 7, 268–313. doi: 10.3390/md7030268
Aziz, R. K., Dwivedi, B., Akhter, S., Breitbart, M., and Edwards, R. A. (2015). Multidimensional metrics for estimating phage abundance, distribution, gene density, and sequence coverage in metagenomes. Front. Microbiol. 6:381. doi: 10.3389/fmicb.2015.00381
Barrangou, R., Fremaux, C., Deveau, H., Richards, M., Boyaval, P., Moineau, S., et al. (2007). CRISPR provides acquired resistance against viruses in prokaryotes. Science 315, 1709–1712. doi: 10.1126/science.1138140
Beretta, E., and Kuang, Y. (1998). Modeling and analysis of a marine bacteriophage infection. Math. Biosci. 149, 57–76. doi: 10.1016/S0025-5564(97)10015-3
Bergh, O., Borsheim, K. Y., Bratbak, G., and Heldal, M. (1989). High abundance of viruses found in aquatic environments. Nature 340, 467–468. doi: 10.1038/340467a0
Blazquez, J., Oliver, A., and Gomez-Gomez, J. M. (2002). Mutation and evolution of antibiotic resistance: antibiotics as promoters of antibiotic resistance? Curr. Drug Targets 3, 345–349. doi: 10.2174/1389450023347579
Bohannan, B. J. M., and Lenski, R. E. (2000). Linking genetic change to community evolution: insights from studies of bacteria and bacteriophage. Ecol. Lett. 3, 362–377. doi: 10.1046/j.1461-0248.2000.00161.x
Breitbart, M., and Rohwer, F. (2005). Here a virus, there a virus, everywhere the same virus? Trends Microbiol. 13, 278–284. doi: 10.1016/j.tim.2005.04.003
Brockhurst, M. A., Morgan, A. D., Fenton, A., and Buckling, A. (2007). Experimental coevolution with bacteria and phage: the Pseudomonas fluorescens φ2 model system. Infect. Genet. Evol. 7, 547–552. doi: 10.1016/j.meegid.2007.01.005
Buckling, A., and Rainey, P. B. (2002). Antagonistic coevolution between a bacterium and a bacteriophage. Proc. Royal Soc. Lond. B Biol. Sci. 269, 931–936. doi: 10.1098/rspb.2001.1945
Bull, J. J., Vegge, C. S., Schmerer, M., Chaudhry, W. N., and Levin, B. R. (2014). Phenotypic Resistance and the Dynamics of Bacterial Escape from Phage Control. PLoS One 9:e94690. doi: 10.1371/journal.pone.0094690
Cairns, B. J., Timms, A. R., Jansen, V. A. A., Connerton, I. F., and Payne, R. J. H. (2009). Quantitative models of in vitro bacteriophage-host dynamics and their application to phage therapy. PLoS Pathog 5:e1000253. doi: 10.1371/journal.ppat.1000253
Carletti, M. (2007). Mean-square stability of a stochastic model for bacteriophage infection with time delays. Math. Biosci. 210, 395–414. doi: 10.1016/j.mbs.2007.05.009
Carlson, K. (2005). “Working with bacteriophages: common techniques and methodological approaches,” in Bacteriophages: Biology and use, ed. E.K.y.A Sulakvelidze (Boca Raton, FL: CRC Press), 437–494. doi: 10.1201/9780203491751.ax1
Carvalho, E. D., Gianmarco Silva, D., and Silva, R. J. (2012). Health and Environment in Aquaculture. Rijeka: InTech. doi: 10.5772/2462
Chang, M. L., and Chang, T. S. (1969). Direct solution of Markovian phage attachment to bacteria in suspension. Math. Biosci. 5, 9–18. doi: 10.1016/0025-5564(69)90032-7
Chelliah, V., Juty, N., Ajmera, I., Ali, R., Dumousseau, M., Glont, M., et al. (2014). BioModels: ten-year anniversary. Nucl. Acids Res. 43, D542–D548. doi: 10.1093/nar/gku1181
Clokie, M. R. J., Millard, A. D., Letarov, A. V., and Heaphy, S. (2011). Phages in nature. Bacteriophage 1, 31–45. doi: 10.4161/bact.1.1.14942
Cochlan, W. P., Wikner, J., Steward, G. F., Smith, D. C., and Azam, F. (1993). Spatial distribution of viruses, bacteria and chlorophyll a in neritic, oceanic and estuarine environments. Mar. Ecol. Prog. Ser. 92, 77–87. doi: 10.3354/meps092077
D’Herelle, F. (1917). Sur un microbe invisible an- tagoniste des bacilles dysentkriques. Comp. Rend. Acad. Sci. 165, 373–375.
Drake, J. W., Charlesworth, B., Charlesworth, D., and Crow, J. F. (1998). Rates of spontaneous mutation. Genetics 148, 1667–1686.
Echeverría-Vega, A., Morales-Vicencio, P., Saez-Saavedra, C., Ceh, J., and Araya, R. (2016). The complete genome sequence and analysis of vB_VorS-PVo5, a Vibrio phage infectious to the pathogenic bacterium Vibrio ordalii ATCC-33509. Stand. Genomic Sci. 11, 1–8. doi: 10.1186/s40793-016-0166-6
Echeverría-Vega, A., Morales-Vicencio, P., Saez-Saavedra, C., Gordillo-Fuenzalida, F., and Araya, R. (2019). A rapid and simple protocol for the isolation of bacteriophages from coastal organisms. MethodsX 6, 2614–2619. doi: 10.1016/j.mex.2019.11.003
Ellis, E. L., and Delbrück, M. (1939). The growth of bacteriophage. J. Gen. Physiol. 22, 365–384. doi: 10.1085/jgp.22.3.365
Gani, J. (1965). Stochastic phage attachment to Bacteria. Biometrics 21, 134–139. doi: 10.2307/2528357
Harbottle, H., Thakur, S., Zhao, S., and White, D. G. (2006). Genetics of Antimicrobial Resistance. Anim. Biotechnol. 17, 111–124. doi: 10.1080/10495390600957092
Heilmann, S., Sneppen, K., and Krishna, S. (2012). Coexistence of phage and bacteria on the boundary of self-organized refuges. Proc. Natl. Acad. Sci. U.S.A. 109, 12828–12833. doi: 10.1073/pnas.1200771109
Hewson, I., O’Neil, J. M., Fuhrman, J. A., and Dennison, W. C. (2001). Virus-Like particle distribution and abundance in sediments and overlying waters along eutrophication gradients in two subtropical estuaries. Limnol. Oceanogr. 46, 1734–1746. doi: 10.4319/lo.2001.46.7.1734
Hyman, P., and Abedon, S. T. (2009). “Practical Methods for Determining Phage Growth Parameters,” in Bacteriophages: Methods and Protocols, Volume 1: Isolation, Characterization, and Interactions, eds M. R. J. Clokie and A. M. Kropinski (Totowa, NJ: Humana Press), 175–202. doi: 10.1007/978-1-60327-164-6_18
Joiner, K. L., Baljon, A., Barr, J., Rohwer, F., and Luque, A. (2019). Impact of bacteria motility in the encounter rates with bacteriophage in mucus. Sci. Rep. 9:16427. doi: 10.1038/s41598-019-52794-2
Keen, E. C., and Adhya, S. L. (2015). Review of phage therapy: current research and applications. Clin. Infect. Dis. 61, 141-142. doi: 10.1093/cid/civ257
Koch, A. L. (1960). Encounter efficiency of coliphage-bacterium interaction. Biochim. Biophys. Acta 39, 311–318. doi: 10.1016/0006-3002(60)90168-2
Kropinski, A. M. (2009). “Measurement of the Rate of Attachment of Bacteriophage to Cells,” in Bacteriophages: Methods and Protocols, Volume 1: Isolation, Characterization, and Interactions, eds M. R. J. Clokie and A. M. Kropinski (Totowa, NJ: Humana Press), 151–155. doi: 10.1007/978-1-60327-164-6_15
Kropinski, A. M., Mazzocco, A., Waddell, T. E., Lingohr, E., and Johnson, R. P. (2009). “Enumeration of Bacteriophages by Double Agar Overlay Plaque Assay,” in Bacteriophages: Methods and Protocols, Volume 1: Isolation, Characterization, and Interactions, eds M. R. J. Clokie and A. M. Kropinski (Totowa, NJ: Humana Press), 69–76. doi: 10.1007/978-1-60327-164-6_7
Krueger, A. P. (1931). The sorption of bacteriophage by living and dead susceptible bacteria : I. Equilibrium conditions. J. Gen. Physiol. 14, 493–516. doi: 10.1085/jgp.14.4.493
Kutter, E. M., Kuhl, S. J., and Abedon, S. T. (2015). Re-establishing a place for phage therapy in western medicine. Fut. Microbiol. 10, 685–688. doi: 10.2217/fmb.15.28
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. Natural. 125, 585–602. doi: 10.1086/284364
Levin, B. R., and Bull, J. J. (1996). Phage therapy revisited: the population biology of a bacterial infection and its treatment with bacteriophage and antibiotics. Am. Natural. 147, 881–898. doi: 10.1086/285884
Lowe, G., Meister, M., and Berg, H. (1987). Rapid rotation of flagellar bundles in swimming bacteria. Nature 325, 637–640. doi: 10.1038/325637a0
Madhusudana Rao, B., and Lalitha, K. V. (2015). Bacteriophages for aquaculture: are they beneficial or inimical. Aquaculture 437, 146–154. doi: 10.1016/j.aquaculture.2014.11.039
Matsuzaki, S., Rashel, M., Uchiyama, J., Sakurai, S., Ujihara, T., Kuroda, M., et al. (2005). Bacteriophage therapy: a revitalized therapy against bacterial infectious diseases. J. Infect. Chemother. 11, 211–219. doi: 10.1007/s10156-005-0408-9
Milo, R., Jorgensen, P., Moran, U., Weber, G., and Springer, M. (2010). BioNumbers—the database of key numbers in molecular and cell biology. Nucl. Acids Res. 38, (Suppl. 1), D750–D753. doi: 10.1093/nar/gkp889
Mitchell, J. G., Pearson, L., Dillon, S., and Kantalis, K. (1995). Natural assemblages of marine bacteria exhibiting high-speed motility and large accelerations. Appl. Environ. Microbiol. 61, 4436–4440. doi: 10.1128/AEM.61.12.4436-4440.1995
Moldovan, R., Chapman-McQuiston, E., and Wu, X. L. (2007). On kinetics of phage adsorption. Biophys. J. 93, 303–315. doi: 10.1529/biophysj.106.102962
Murray, A. G., and Jackson, G. A. (1992). Viral Dynamics: a model of the effects of size, shape, motion and abundance of single-celled planktonic organisms and other particles. Mar. Ecol. Prog. Ser. 89, 103–116. doi: 10.3354/meps089103
Paul, J. H., Rose, J. B., Jiang, S. C., Kellogg, C. A., and Dickson, L. (1993). Distribution of viral abundance in the reef environment of Key Largo, Florida. Appl. Environ. Microbiol. 59, 718–724. doi: 10.1128/AEM.59.3.718-724.1993
Payne, R. J. H., and Jansen, V. A. A. (2000). Phage therapy: the peculiar kinetics of self-replicating pharmaceuticals. Clin. Pharmacol. Therap. 68, 225–230. doi: 10.1067/mcp.2000.109520
Payne, R. J. H., and Jansen, V. A. A. (2001). Understanding bacteriophage therapy as a density-dependent kinetic process. J. Theor. Biol. 208, 37–48. doi: 10.1006/jtbi.2000.2198
Porter, K. G., and Feig, Y. S. (1980). The use of DAPI for identifying and counting aquatic microflora. Limnol. Oceanogr. 25, 943–948. doi: 10.4319/lo.1980.25.5.0943
Rabinovitch, A., Aviram, I., and Zaritsky, A. (2003). Bacterial debris—an ecological mechanism for coexistence of bacteria and their viruses. J. Theor. Biol. 224, 377–383. doi: 10.1016/S0022-5193(03)00174-7
Riquelme, C., Hayashida, G., Toranzo, A., Vilches, J., and Chavez, P. (1995). Pathogenicity studies on a Vibrio anguillarum-related (VAR) strain causing an epizootic in Argopecten purpuratus larvae cultured in Chile. Dis. Aqua. Organ. 22, 135–141. doi: 10.3354/dao022135
Samson, J. E., Magadan, A. H., Sabri, M., and Moineau, S. (2013). Revenge of the phages: defeating bacterial defences. Nat. Rev. Micro. 11, 675–687. doi: 10.1038/nrmicro3096
Sarhan, W. A., and Azzazy, H. M. (2015). Phage approved in food, why not as a therapeutic? Exp. Rev. Anti Infect. Ther. 13, 91–101. doi: 10.1586/14787210.2015.990383
Schiewe, M. H., Trust, T. J., and Crosa, J. H. (1981). Vibrio ordalii sp. nov.: a causative agent of vibriosis in fish. Curr. Microbiol. 6, 343–348. doi: 10.1007/BF01567009
Schrag, S. J., and Mittler, J. E. (1996). Host-parasite coexistence: the role of spatial refuges in stabilizing bacteria-phage interactions. Am. Nat. 148, 348–377. doi: 10.1086/285929
Schwartz, M. (1976). The adsorption of coliphage lambda to its host: effect of variations in the surface density of receptor and in phage-receptor affinity. J. Mol. Biol. 103, 521–536. doi: 10.1016/0022-2836(76)90215-1
Stern, A., and Sorek, R. (2011). The phage-host arms-race: Shaping the evolution of microbes. Bioessays 33, 43–51. doi: 10.1002/bies.201000071
Stopar, D. (2008). Modeling Bacteriophage Population Growth. Bacteriophage Ecology. Cambridge: Cambridge University Press.
Thingstad, T. F. (2000). Elements of a theory for the mechanisms controlling abundance, diversity, and biogeochemical role of lytic bacterial viruses in aquatic systems. Limnol. Oceanogr. 45, 1320–1328. doi: 10.4319/lo.2000.45.6.1320
Verbeken, G., Huys, I., Pirnay, J.-P., Jennes, S., Chanishvili, N., Scheres, J., et al. (2014). Taking bacteriophage therapy seriously: a moral argument. Biomed Res. Int. 2014:8. doi: 10.1155/2014/621316
Wei, Y., Ocampo, P., and Levin, B. R. (2010). An experimental study of the population and evolutionary dynamics of Vibrio cholerae O1 and the bacteriophage JSF4. Proc. Royal Soc. B Biol. Sci. 277, 3247–3254. doi: 10.1098/rspb.2010.0651
Weitz, J. S., Hartman, H., and Levin, S. A. (2005). Coevolutionary arms races between bacteria and bacteriophage. Proc. Natl. Acad. Sci. U.S.A. 102, 9535–9540. doi: 10.1073/pnas.0504062102
Wommack, K. E., and Colwell, R. R. (2000). Virioplankton: viruses in aquatic ecosystems. Microbiol. Mol. Biol. Rev. 64, 69–114. doi: 10.1128/MMBR.64.1.69-114.2000
Xie, L., Altindal, T., Chattopadhyay, S., and Wu, X.-L. (2011). Bacterial flagellum as a propeller and as a rudder for efficient chemotaxis. Proc. Natl. Acad. Sci. U.S.A. 108:2246. doi: 10.1073/pnas.1011953108
Yassky, D. (1962). A Model for the Kinetics of Phage Attachment to Bacteria in Suspension. Biometrics 18, 185–191. doi: 10.2307/2527456
Keywords: bacteriophage adsorption dynamics, model for bacteriophage infection, phage-host interaction, particle collision, phage therapy
Citation: Echeverría-Vega A, Morales-Vicencio P, Saez-Saavedra C, Alvarez MA, Gordillo F, Del-Valle R, Solís ME and Araya R (2020) Characterization of the Bacteriophage vB_VorS-PVo5 Infection on Vibrio ordalii: A Model for Phage-Bacteria Adsorption in Aquatic Environments. Front. Microbiol. 11:550979. doi: 10.3389/fmicb.2020.550979
Received: 11 April 2020; Accepted: 11 September 2020;
Published: 30 October 2020.
Edited by:
Robert Czajkowski, University of Gdańsk, PolandReviewed by:
Jeremy J. Barr, Monash University, AustraliaPanos G. Kalatzis, University of Copenhagen, Denmark
Copyright © 2020 Echeverría-Vega, Morales-Vicencio, Saez-Saavedra, Alvarez, Gordillo, Del-Valle, Solís and Araya. 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: Alex Echeverría-Vega, aecheverria@ucm.cl; Rubén Araya, ruben.araya@uantof.cl