
Citation: Celina Gómez, Luigi Gennaro Izzo. Increasing efficiency of crop production with LEDs[J]. AIMS Agriculture and Food, 2018, 3(2): 135-153. doi: 10.3934/agrfood.2018.2.135
[1] | Cristina Anton, Alan Yong . Stochastic dynamics and survival analysis of a cell population model with random perturbations. Mathematical Biosciences and Engineering, 2018, 15(5): 1077-1098. doi: 10.3934/mbe.2018048 |
[2] | Huan Kong, Guohong Zhang, Kaifa Wang . Stability and Hopf bifurcation in a virus model with self-proliferation and delayed activation of immune cells. Mathematical Biosciences and Engineering, 2020, 17(5): 4384-4405. doi: 10.3934/mbe.2020242 |
[3] | Cameron Browne . Immune response in virus model structured by cell infection-age. Mathematical Biosciences and Engineering, 2016, 13(5): 887-909. doi: 10.3934/mbe.2016022 |
[4] | Allen L. Nazareno, Raymond Paul L. Eclarin, Eduardo R. Mendoza, Angelyn R. Lao . Linear conjugacy of chemical kinetic systems. Mathematical Biosciences and Engineering, 2019, 16(6): 8322-8355. doi: 10.3934/mbe.2019421 |
[5] | Amer Hassan Albargi, Miled El Hajji . Mathematical analysis of a two-tiered microbial food-web model for the anaerobic digestion process. Mathematical Biosciences and Engineering, 2023, 20(4): 6591-6611. doi: 10.3934/mbe.2023283 |
[6] | Giulia Bertaglia, Lorenzo Pareschi, Giuseppe Toscani . Modelling contagious viral dynamics: a kinetic approach based on mutual utility. Mathematical Biosciences and Engineering, 2024, 21(3): 4241-4268. doi: 10.3934/mbe.2024187 |
[7] | Hsiu-Chuan Wei . Mathematical modeling of tumor growth: the MCF-7 breast cancer cell line. Mathematical Biosciences and Engineering, 2019, 16(6): 6512-6535. doi: 10.3934/mbe.2019325 |
[8] | Yan Fu, Tian Lu, Meng Zhou, Dongwei Liu, Qihang Gan, Guowei Wang . Effect of color cross-correlated noise on the growth characteristics of tumor cells under immune surveillance. Mathematical Biosciences and Engineering, 2023, 20(12): 21626-21642. doi: 10.3934/mbe.2023957 |
[9] | Giulio Caravagna, Alex Graudenzi, Alberto d’Onofrio . Distributed delays in a hybrid model of tumor-Immune system interplay. Mathematical Biosciences and Engineering, 2013, 10(1): 37-57. doi: 10.3934/mbe.2013.10.37 |
[10] | Svetlana Bunimovich-Mendrazitsky, Yakov Goltser . Use of quasi-normal form to examine stability of tumor-free equilibrium in a mathematical model of bcg treatment of bladder cancer. Mathematical Biosciences and Engineering, 2011, 8(2): 529-547. doi: 10.3934/mbe.2011.8.529 |
A large variety of models of interactions between agents have been developed within the framework of thermostatted kinetic theory [1,2] and notably applied to traffic [3], swarming [4] and cell competition [5,6,7,8]. The interactions between immune system and cancer are known to be complex. The competition between immune system cells and cancer cells can be compared to the interactions between immune system cells and viruses. In both cases, pathogenic cells induce immune system responses [9,10,11,12]. Different cell types are involved in the defense mechanism. Three types of cells are commonly recognized as being involved in the immune response, antigen-presenting cells, T cells and B cells. Dentritic cells ingest cancer cells, isolate antigens and present them to T cells, which leads to their activation and proliferation [13]. Activated T cells mutate into killer T cells, which destroy cancer cells. Like innate dendritic cells, B cells are antigen-presenting cells developed by the adaptive immune system. B cells produce antibodies, specific to each pathogen, which bind to the surface of an invading cell and mark it for destruction by T cells. When a cancer cell is killed, additional antigens are released, which stimulates further activation, i.e. learning of new T and B cells. In addition, regulatory T cells are produced to control the response of immune system. They adjust the level of antibody secretion, i.e. the regulation of B and T cell production and activation. While T cells are activated, cancer cells mutate and acquire the ability to avoid detection and destruction by the immune system. Specifically, cancer cells express proteins on their surface that induce the inactivation of immune cells, thus promoting the proliferation of cancer cells [14,15]. Through cellular interactions, cancer cells learn not to be detected, misleading the regulatory T cells, which in turn limits the production of immune system cells and let the cancer cells proliferate [16,17]. In the model, the learning processes and activation of both immune system cells and cancer cells is reproduced by the increase of a real quantity called activity and carried by each cell [2].
Stimulating the immune system of a patient can lead to disappointing results and conversely, chemotherapy may boost immune system response [18]. For some cancers, immunotherapy leads to a long remission followed by a fatal outcome, reflecting the complex relationship between the immune system and the cancer. This phenomenon, known as the three Es characterizes the following three phases commonly observed in a patient, elimination, equilibrium and escape of cancer, introduced in the theory of immunosurveillance as early as the mid-20th century [9,19]. The initial elimination of cancer cells is only apparent and corresponds to the attack by killer T cells, before the cancer cells have learned to mislead the immune system. During the pseudo-equilibrium phase, a small number of cancer cells remain and learn from the immune system cells. Eventually, the cancer cells are able to proliferate and the cancer escapes the control of the immune system. C-Reactive Protein (CRP) levels have been used as markers of inflammation in response to antigen exposure in relation to the progression of many cancers, such as melanoma and lung cancer [20]. Oscillations of CRP concentration in advanced cancer patients have been observed and interpreted in terms of immune regulatory cycles [21,22]. The oscillations of CRP concentration reveal the antagonistic effects of immune system stimulation on a tumor. Many macroscopic models of interactions between cancer and immune system leading to oscillating solutions for concentrations after a Hopf bifurcation can be found in the recent literature [23,24,25,26], in particular after introduction of one or more delays in the equations [27,28], but without direct connection with a description at cell scale.
We recently designed a minimal model of interactions between cancer and immune system cells in a homogeneous system, ignoring the spatial description of a tumor [7,8]. The spatio-temporal description of the immune system [29] or tumors [30] often involves partial differential equations for macroscopic quantities, with nevertheless an emphasis on the importance of multiscale modeling [31]. In this work we introduce a space-velocity description of tumor growth at cell scale in the framework of thermostatted kinetic theory. The homogeneous model includes autocatalytic cell production and increased cell activity during favorable interactions. In chemical kinetics, autocatalysis refers to the formation of a product which is also a reactant. Hence, the production rate of the autocatalytically formed species depends on its own concentration, leading to long induction times, when the species is present in low concentration at the beginning of the reaction, and then to explosive behaviors, when the concentration exceeds a threshold. Autocatalysis is a crucial ingredient in the destabilization of a steady state and the emergence of more complex behaviors such as oscillations [32,33]. In the model, the fluctuations of cell activity are regulated by a so-called thermostat [2]. The introduced interaction processes only consider the death of a cell when its activity is not favorable. However, all cells [34], including memory T cells involved in the long-lasting protection acquired after vaccination [35,36,37,38], have a finite life span, whether or not they are very active. Thus, the spontaneous death of educated cells leads to a loss of information. In order to take this phenomenon into account, the model introduces the regulation of activity fluctuations, by imposing that its second moment remains constant. By analogy with the dissipation of energy in a system maintained at constant temperature [1], the activity regulator is also called a thermostat.
We design a kinetic Monte Carlo algorithm inspired by the direct simulation Monte Carlo (DSMC) method initially devoted to dilute gases [39,40] and simulate the evolution of cell position, velocity, nature, and activity. Contrary to the standard DSMC method in which the velocities of the dilute gas particles are updated due to binary collisions, jump processes are introduced to randomize the direction of cell velocities in the present inhomogeneous model involving dense soft matter. Specifically, the algorithm that we develop harnesses the discretization of space into boxes and the implementation of interactions only between cells randomly chosen inside a same spatial box, as suggested by standard DSMC. In spite of this powerful procedure, which earns the name of Monte Carlo to the algorithm, the simulations remain expensive in computation time, because of the large number of cancer cells often created in the case of cancer escape. Our aim is to check if the inhomogeneous model contains enough biological ingredients to reproduce the complex antagonistic effects observed in cancer evolution, including the "three Es" (elimination, equilibrium, and escape of cancer from immunosurveillance) [19,26] and spontaneous oscillations in the number of immune system cells observed during tumor growth [21,22] but not included in the homogeneous model [7,8]. In addition, the simulations of an inhomogeneous system reported in the present paper give access to the evolution of the geometry of the tumor, its possible deformation, shifting, and splitting that were clearly not affordable in the framework of a homogeneous description.
The paper is organized as follows. In section 2, we present the inhomogeneous model and the kinetic equations governing the evolution of the probability distributions of cancer, immune system and normal cells. The kinetic Monte Carlo algorithm is made precise. Simulation results for an initially localized tumor are presented in section 3. Specifically, we examine if, for a poorly effective thermalization, the number of cancer cells rapidly drops, then slightly increases during a long quasi-steady regime, and eventually explodes, although the autocatalytic production of immune system cells is much more favored than cancer cell formation. The existence of oscillations for the total number of immune system cells is also discussed. A sensitivity analysis is performed and the different behaviors observed as the field that regulates activity and the rate constant related to the production of immune system cells vary are shown. The variations of the pseudo-period of oscillations with the field and the rate constant are determined. Section 4 contains conclusions.
The competition between immune system cells and cancer cells can be compared to the interactions between immune system cells and viruses. In both cases, the foreign cells induce responses from the immune system [10,11,12]. Different types of cells are involved in the defense mechanism [41,42]. Dentritic cells ingest cancer cells, isolate antigens and present them to T cells, leading to their activation and proliferation [13]. Killer T cells are able to destroy cancer cells and regulatory T cells control the number of lymphocytes and prevent their overproduction [18,43]. In parallel with the learning of the immune system cells, the learning of the cancer cells takes place, which may lead to cancer escape. Through cell interactions, the cancer cells learn not to be detected, misleading the regulatory T cells which in turn limit the production of immune system cells and let cancer cells proliferate [16,17].
Recently, we have introduced a homogeneous model of competition between cancer cells and immune system cells which reproduces cell interactions, activation, and learning processes [7,8] and accounts for the three Es of immunotherapy. We propose to develop a space-velocity description of tumor growth, which requires to model the updating of cell position x and velocity v in addition to their nature j and activity u. As in the homogeneous model, we introduce normal cells n, cancer cells c, and immune system cells i as a whole with the aim of introducing as few variables and parameters as possible and providing qualitative interpretations of clinical observations.
Inhomogeneous initial conditions with tumor cells located in the middle of the system are preferably considered. A Gaussian distribution Pu(u) of mean μ and variance σ is used to generate the initial activities of the cells, regardless of their nature. Possible heterogeneity in initial activity distribution may be considered. All the cells are supposed to have the same speed which remains constant. The direction of cell velocity v is defined by an angle θv∈]0,2π] randomly chosen using a uniform distribution Pv(θv).
Cell activity u and nature j=c,i,n, where c stands for cancer, i for immune, and n for normal, may change during local interactions that are assumed to leave position x and velocity v of cancer cells and immune system cells unchanged. We consider only three processes which include interaction, activation, proliferation or death, in a similar manner as autocatalytic chemical processes
{c(v,u)+n(v′,u′)κcn(v,u,v′,u′)→c(v,u+ε)+c(v′,u′)R⟶n(v′′,u′′) | (2.1) |
i(v,u)+c(v′,u′)κic(v,u,v′,u′)→i(v,u+ε)+i(v′,u′) | (2.2) |
c(v,u)+i(v′,u′)κci(v,u,v′,u′)→c(v,u+ε)+c(v′,u′) | (2.3) |
where the rate constants κcn(v,u,v′,u′), κic(v,u,v′,u′), and κci(v,u,v′,u′) have the following nontrivial dependence on the activities and the velocities of the interacting couple:
κcj(v,u,v′,u′)=kcj(u−u′)H(u−u′)H(θinter−θ(v,v′))forj=n,i | (2.4) |
κic(v,u,v′,u′)=kic(u−u′)H(u−u′)H(θinter−θ(v,v′)) | (2.5) |
where kcn, kic and kci are constant, θinter∈]0,π] is the maximum angle between the velocity directions permitting interaction and θ(v,v′) is defined as follows:
θ(v,v′)={∣θv−θv′∣if∣θv−θv′∣≤π∣∣θv−θv′∣−2π∣if∣θv−θv′∣>π | (2.6) |
where θv∈]0,2π] is the angle defining the direction of v. The condition on the directions of the velocities of the interacting pair imposed by the Heaviside function H(θinter−θ(v,v′)) ensures that the two cells remain close during a sufficient time for activation and learning to take place.
The learning process is reproduced by the small increase ε of the activity of the cell with the already larger activity before the interaction. The process described in Eq. (2.1) has the form of an autocatalytic production of cancer cells, the presence of a cancer cell being a pre-requisite for the formation of another cancer cell, as in a chemical process catalysed by a product of the reaction.
Equation (2.1) accounts for both the mutation of a normal cell into a cancer cell and the division of a cancer cell. The specific rate constant expression ensures that proliferation of cancer cells is more favorable when they have already learned to blend into their environment, i.e. have increased their activity u. The reservoir R of normal cells maintains their number constant and plays also the role of reservoir of activity. The notion of reservoir is useful in far-from-equilibrium systems able to self-organize. The formation of complex structures requires that the system is maintained far from equilibrium, which occurs through the consumption of energy or matter. A reservoir is an ideal device which is supposed to be able to maintain the concentration of a given species constant by instantaneously absorbing incoming fluxes or injecting the species to compensate for losses [32]. In the present model, the condition on the activities of the interacting cells imposed by the Heaviside function H(u−u′) introduces a selection. The normal cell which disappears has a smaller activity than the interacting cancer cell, which could induce a non physical increase in the mean activity of the remaining normal cells. The model provides for reinjecting a normal cell with a randomly chosen activity u″ using the initial Gaussian distribution Pu(u) of fixed mean value μ, so that the mean activity of normal cells remains constant, consistent with biology. The velocity v″ of the reinjected normal cell is randomly chosen using the initial uniform distribution Pv(θv).
The process given in Eq. (2.2) corresponds to the autocatalytic production of immune system cells and accounts for both the death of a cancer cell attacked by a killer T cell and the proliferation of immune system cells with increased activity, i.e. the knowledge of the presence of cancer cells in their environment. As in Eq. (2.1), the rate constant depends on a Heaviside step function H(u−u′) implying that the increase of the number of immune cells is locally more probable when they have a larger activity than the neighboring cancer cells.
The regulation of the number of immune cells by T-regulatory cells and the learning process of cancer cells are taken into account in the third process given in Eq. (2.3). This process accounts for both the death of an immune system cell and the division of a cancer cell. The rate constant being proportional to H(u−u′), the proliferation of cancer cells locally occurs if their activity is larger than the one of the surrounding immune system cells. Hence, the model accounts for the ability of trained cancer cells to mislead the immune system and trigger the detrimental response of T-regulatory cells.
In order to introduce a regulation of the explosive cell production due to autocatalytic processes, we impose that the rate of each process is proportional to the difference of activities of the interacting pair. Specifically, according to Eq. (2.3), further growth of a cluster of cancer cells with a smaller activity than the surrounding immune system cells is not likely to be observed. Following Eq. (2.2), an analogous phenomenon is expected for immune system cells: A local increase in the number of immune system cells will not be sustained if their activity is smaller that the one of the neighboring cancer cells. According to Eqs. (2.4) and (2.5), the expression of the rate constants, proportional to the relative activity of the interacting cells introduces a remote analogy to the collision term of the Boltzmann equation, which depends on the relative velocity of the encounters [44].
In order to mimic the preferential increase in the number of cancer cells due to the deleterious response of T-regulatory cells rather than by mutation of normal cells, we assign values of kci in Eq. (2.3) larger than the values of kcn in Eq. (2.1).
The processes given in Eqs. (2.1-2.3) associate learning process with systematic increased cell activity. However, cell interactions that are not explicitly taken into account in the scheme make learning less effective. In the model the regulation of the activity fluctuations is ensured by a thermostat mimicking information loss due to dissipation. Thermalization, characterized by a field E and a friction coefficient α, acts on the activity uk of each cell k according to [1,5]
dukdt=E−αuk | (2.7) |
Choosing
α=⟨u⟩E⟨u2⟩ | (2.8) |
where ⟨u⟩ is the mean activity of the entire system, ensures that the second moment ⟨u2⟩ remains constant.
All cells are supposed to move at constant speed in all directions. We introduce a stochastic turning operator that modifies the direction of cell velocities inside a cone of apex angle θjump with 0<θjump≤π/2. Hence, the rate constant associated with the velocity jumps obeys
κv(v,v′)=kvH(θjump−∣θv−θv′∣) | (2.9) |
where kv is constant.
The distribution function fj for cells of type j=c,i,n obeys the following kinetic equation
(∂t+v⋅∇x)fj(t,x,v,u)+∂u((E−αu)fj(t,x,v,u))=Ij+Vj | (2.10) |
where v⋅∇x is the standard advection operator, ∂u((E−αu)fj) is the thermalization operator, Ij is the interaction operator and Vj is the velocity randomization operator. Following Eqs. (2.1-2.3), the interaction operator Ic associated with cancer cells reads
Ic=∫R3∫R+κcn(v,u−ε,v′,u′)fc(t,x,v,u−ε)fn(t,x,v′,u′)dv′du′+∫R3∫R+κcn(v′,u′,v,u)fc(t,x,v′,u′)fn(t,x,v,u)dv′du′−∫R3∫R+κcn(v,u,v′,u′)fc(t,x,v,u)fn(t,x,v′,u′)dv′du′−∫R3∫R+κic(v′,u′,v,u)fi(t,x,v′,u′)fc(t,x,v,u)dv′du′+∫R3∫R+κci(v,u−ε,v′,u′)fc(t,x,v,u−ε)fi(t,x,v′,u′)dv′du′+∫R3∫R+κci(v′,u′,v,u)fc(t,x,v′,u′)fi(t,x,v,u)dv′du′−∫R3∫R+κci(v,u,v′,u′)fc(t,x,v,u)fi(t,x,v′,u′)dv′du′ | (2.11) |
with κcj(v,u,v′,u′) and κic(v,u,v′,u′) given in Eqs. (2.4) and (2.5). The two first terms of the right-hand side account for the positive contributions to the distribution function of cancer cells of velocity v and activity u after the interaction of a cancer cell and a normal cell through the process given in Eq. (2.1). Specifically, the first term corresponds to the formation of a cancer cell with the desired properties from a cancer cell of initial activity u−ε after update by the increment ε. All contributions are obtained by integrating over all acceptable velocities v′ and activities u′ of the interacting normal cell. The second term is obtained when assigning the velocity v and activity u of the disappearing normal cell to a newly formed cancer cell and integrating over all acceptable velocities v′ and activities u′ of the interacting cancer cell. The fourth term is related to the process given in Eq. (2.2) and the fifth to seventh terms are associated with Eq. (2.3). The interaction operator Ii for immune system cells is formed following analogous rules as for cancer cells:
Ii=∫R3∫R+κic(v,u−ε,v′,u′)fi(t,x,v,u−ε)fc(t,x,v′,u′)dv′du′+∫R3∫R+κic(v′,u′,v,u)fi(t,x,v′,u′)fc(t,x,v,u)dv′du′−∫R3∫R+κic(v,u,v′,u′)fi(t,x,v,u)fc(t,x,v′,u′)dv′du′−∫R3∫R+κci(v′,u′,v,u)fc(t,x,v′,u′)fi(t,x,v,u)dv′du′ | (2.12) |
For normal cells, the interaction operator In reads
In=−∫R3∫R+κcn(v′,u′,v,u)fc(t,x,v′,u′)fn(t,x,v,u)dv′du′+Pu(u)Pv(θv)∬R3∬R+κcn(v′,u′,v″,u″)fc(t,x,v′,u′)fn(t,x,v″,u″)dv′dv″du′du″ | (2.13) |
The second term of the right-hand side is related to the action of the reservoir, which injects a normal cell of velocity v and activity u exactly at the same rate as a normal cell of velocity v″ and activity u″ just disappeared through Eq. (2.1). The contribution to the evolution of fn(t,x,v,u) includes integration over all velocities v″ and weighting by Pv(θv), integration over all activities u″ and weighting by Pu(u), as well as integration over all acceptable velocities v′ and activities u′ of the interacting cancer cell.
Finally, the velocity randomization operator Vj which models the velocity-jump process inside a cone reads
Vj=∫R3κv(v,v′)(fj(t,x,v′,u)−fj(t,x,v,u))dv′ | (2.14) |
where the rate constant κv(v,v′) of clockwise or anticlockwise jump of velocity direction is given in Eq. (2.9).
Graeme Bird [39,40] introduced an efficient algorithm, the direct simulation Monte Carlo (DSMC) method, to simulate the Boltzmann equations associated with a dilute gas [45] and successfully adapted to concentrated solutions [46]. The performance of DSMC with respect to molecular dynamics is related to considering collisions only for particles belonging to a same spatial box and following an acceptance-rejection technique to treat the collisions. Our aim is to show that these hypotheses can be satisfactorily extended to cellular dynamics in order to directly solve the three kinetic equations given in Eq. (2.10) for the different cell types j=c,i,n. We propose the following algorithm. During a time step Δt, cell natures, activities, positions and velocities are updated according to the model of interactions, thermalization and velocity jumps described in the previous section. Specifically, for given values of the constants kcn, kic, and kci, the maximum number of interactions between pairs is determined in each spatial box and the corresponding interactions between randomly chosen pairs are actually accepted provided that the conditions on activities and velocities expressed in Eqs. (2.4) and (2.5) are met. The nature and the activity of the interacting cells are then updated. After the interactions in all spatial boxes have been considered, the activity of each cell is thermalized according to Eqs. (2.7). Cell positions are modified according to their velocities, which may result in changing of spatial box. Finally, the directions of cell velocities are updated following the stochastic deflection process of rate constant given in Eq. (2.9).
In the next section we present the simulation results obtained for a two-dimensional inhomogeneous system.
Two-dimensional simulations are performed in a square of side lx=ly=30 divided in boxes of side Δx=Δy=1. Periodic boundary conditions are applied. Initially, 20 normal cells and 20 immune system cells per box are homogeneously spread in the entire system. In order to study the growth of an initially localized tumor, we start from 20 cancer cells per box in a central discretized disc of diameter D=5 and no cancer cells elsewhere. The effects of the rate constant values on the system behavior have been studied in detail in a homogeneous system [7,8] and we have checked that they are not sensitively different in the inhomogeneous case. The absolute values of the rate constants are arbitrary. We choose to present the results obtained in a boosted immune system such that kic>kci. This choice results in the more favorable autocatalytic formation of immune system cells than cancer cells for populations with similar activities. We impose kci>kcn in order to comply with the much more probable division of cancer cells under the pressure exerted by regulatory T cells through Eq. (2.3) than the mutation of normal cells into cancer cells reproduced by Eq. (2.1). We impose kci/kcn=103. We refer to the sensitivity analysis performed in a homogeneous system for a detailed discussion of the effect of the parameters introduced in the model [8]. In order to speed up the simulation, we accept the interaction between two cells without condition on the direction of their velocities, which amounts to choosing θinter=π. The rate constant for velocity deflection kv and the maximum angle of velocity jump θjump are chosen sufficiently large for the randomization of velocity direction to be efficient. We have kv=0.1 and θjump=π/6.
Before performing a sensitivity analysis and showing the different behaviors encountered as essential parameters of the model vary, we focus on two nontrivial behaviors reproducing clinical observations, the three Es of immunotherapy [19] and oscillations of the number of immune system cells [21,22]. The same parameter values are chosen to illustrate the two typical behaviors, except the ratio kic/kci, equal to 20 for the 3Es and 10 for the oscillations.
Figure 1 shows the evolution of the total number Nc of cancer cells and the total number Ni of immune system cells in the entire system. Due to the action of the reservoir, the number of normal cells is constant by construction in agreement with biological requirements. In the simulations, the total number of normal cells is set at 18000 in a system of 30×30 spatial boxes. A phenomenon analogous to the three Es (elimination, equilibrium, and escape of cancer) observed in immunotherapy [19] is reproduced by the inhomogeneous model for a critical value of the field E depending on the values of the other parameters. The initial elimination phase observed in a homogeneous system [7,8] is not observed and the number Nc of cancer cells first slightly increases. Then, Nc remains nearly constant during a long pseudo-equilibrium phase. The elimination phase would be recovered for an initial number of cancer cells larger than the value of the plateau reached during the equilibrium phase. Finally, Nc explodes, becomes greater than the number of cells of the immune system, Ni, which in turn fluctuates with a large amplitude. Shortly after a very large fluctuation, Ni cancels out and the cancer escapes from immunosurveillance. In the inhomogeneous system, we rather observe the 2Es, for equilibrium and escape. If one excepts the very beginning of the evolution, the behaviors of the homogeneous and inhomogeneous systems are similar and reproduce a very long induction period before explosion, obtained for sufficiently inefficient thermalization in a small interval of field values.
Cancer cells invade the entire system despite kci<kic: Inefficient control of activity fluctuations severely limits the chances of success of immunotherapy. Following the expression of the rate constant κci given in Eq. (2.4), the local important production of cancer cells is related to the learning of these cells, i.e. to the increase of their activity with respect to the neighboring immune system cells. The learning process implemented in the model consists of a small activity increase through cell interactions and a rate dependence on the relative activity of the interacting pairs. This simple procedure is refined enough to reproduce the complex learning process of cancer cells, which can mislead the immune system and induce a harmful regulatory T cell response that wrongly orders a decrease in the production of immune cells.
The spatial description gives access to the evolution of the clusters of cancer cells and immune system cells. A cluster of cancer cell is defined as a continuous set of spatial boxes with neighboring boxes containing at least 20 cancer cells, having a common side and not only a common corner. The value 20 corresponds to the initial number of cancer cells per box in the small, localized tumor. An analogous definition is used for a cluster of immune system cells. The size of a cluster is given by the number of boxes it contains. We use the Hoshen-Kopelman algorithm to build the clusters [47,48]. Figure 2 gives the evolutions of the averaged number of clusters and the size of the biggest cluster for each cell type for the same conditions as in Figure 1.
During the short transient initial period preceding the pseudo-equilibrium phase, the number of clusters of immune system cells rapidly increases while the size of the biggest cluster decreases from the total system size 30×30 to a small value around 10, revealing the rapid fragmentation of the initial big cluster. During the same period, the number of clusters of cancer cells slightly fluctuates between 1 and 0 and the size of the biggest cluster instantaneously drops from 21 to 0, meaning that the number of cancer cells per box rapidly becomes smaller than 20. During the long pseudo-equilibrium phase, the number of clusters of immune system cells linearly decreases due to two phenomena, the coalescence of growing clusters and the decrease of the number of immune system cells, in agreement with the slight increase of the size of the biggest cluster. During the same phase, the number of clusters and the size of the biggest cluster of cancer cells vanish, as if cancer has been eliminated. Then, the number of clusters of cancer cells exponentially grows implying a global, nonlocal escape of cancer while the number of clusters of immune system cells drops. Eventually, the clusters of cancer cells coalesce and a single cluster invades the entire system while the immune system disappears: Cancer becomes out of control. System size has an impact on the time at which an initially localized tumor invades the entire system but does not affect the qualitative behavior of the system.
The main contribution of the inhomogeneous description is shown in Figure 3. The total numbers Ni of immune system cells and the mean activities, uc, and ui, of both types of cells display pseudo-oscillations of quite well-defined period and increasing amplitude. The total number Nc of cancer cells rapidly becomes much larger than the total number Ni of immune system cells and the chosen scale prevents from showing the end of the evolution of Nc. The extrema of the mean activities uc and ui are associated with typical spatial configurations as shown in the different snapshots of Figures 3-5. A minimum of uc is rapidly followed by a maximum of ui and then a maximum of uc.
As shown in Figure 3, a minimum of uc corresponds to a nearly symmetric disc of cancer cells attacked by a small cluster of immune system cells due to a local favorable fluctuation of ui. Rapidly, the cluster of immune system cells grows and, after a small delay, the mean activity ui reaches a maximum, as displayed in Figure 6. The cluster of immune system cells has an exotic shape which corresponds to a domain of the central tumor in which the number of cancer cells is depleted. According to Figure 5, the symmetry of the entire system is rapidly restored: due to their larger activities in average, the cancer cells deplete the immune system cells in the central part of the tumor, where they are more abundant. As a result, the immune system cells are pushed outwards the disc, forming an external ring. Over a longer period of time, the ring of immune system cells will dissolve, leading to a minimum of Ni, a nearly homogeneous distribution of immune system cells, and a nearly perfect disc of cancer cells. A new local fluctuation of ui will appear and trigger the cycle again.
As exemplified by the Brusselator model [32], a chemical model presenting time oscillations contains the autocatalytic production of one species X to the detriment of another species Y as well as a feedback mechanism which regenerates Y and consumes X. In the model of cancer/immune system interaction, the conditions on the relative activity of an interacting pair of cells given in Eqs. (2.2) and (2.3) introduce such a feedback mechanism. In a small, localized cluster containing a high number of immune system cells with high activity, the autocatalytic production of immune cells by the process given in Eq. (2.2) is first favored. Then, the total number of cancer cells Nc decreases. The effect on the average activity ui is non-trivial, the activity of the already present immune system cell being increased by ε and the newly formed immune system cell inheriting the lower activity of the interacting cancer cell. Simultaneously, the interaction tends to deplete the cancer cell population of its lowest activity elements. As a result, the number of cancer cells decreases but their mean activity ui increases. The autocatalytic production of cancer cells is then promoted by the process given in Eq. (2.3), which in turn tends to decrease the total number of immune cells Ni but increase their average activity ui. The system is again in the situation where the process given in Eq. (2.2) favors the autocatalytic production of immune cells and a new cycle begins. Hence, activity-weighted rate constants for the autocatalytic productions of immune system cells and cancer cells is able to induce an oscillating behavior and reproduce immune regulatory cycles observed in advanced cancer patients and revealed by oscillations of C-Reactive Protein concentration [20,21].
The interaction processes do not only select cells depending on their activity, they also introduce criteria on their velocity direction. Small values of the angle of interaction θinter coupled to inefficient velocity randomization, i.e. small values of both the rate constant kv and the maximum angle of velocity jump θjump, could induce a bias in the results for large values of v/Δt with respect to box size Δx. Specifically, a cluster of autocatalytically formed immune system cells with increased activity would move in a similar direction during a certain time, inducing a persistent anisotropy in cluster shape, which would induce a persistent deformation of tumor shape. Adapted parameter values could be chosen to simulate specific tumor shapes. Conversely, efficient velocity randomization is essential to restore the initial tumor symmetry in the presence of velocity selection during cell interactions. The simulation results shown in Figures 3-5 are given for θinter=π, i.e. in the absence of selection of the angle of the interacting pair, which prevents the formation of clusters with a selected velocity. Cell speed and time step obey v=10−3 and Δt=1, respectively. The choice of kv=0.1 for the rate constant of velocity deflection leads to about 100 velocity jumps before a cell leaves a box, which consists of an efficient velocity randomization. Hence, a randomly occurring cluster of immune system cells surrounded by cancer cells gives rise to the isotropic propagation of a radial front of immune system cells. Then, the cluster depletes from the center. Indeed, the cancer cells in the boxes occupied by the initial cluster of immune system cells which have not been consumed have an increased activity and are thus able to inhibit the immune system cells, leading to the formation of an external ring of immune system cells with radial velocity, similar to the one observed in Figure 4.
The evolution of the averaged number of clusters and size of the biggest cluster are given in Figure 6 for each cell type and for the parameter values of Figures 3-5 leading to oscillations. The number of clusters of cancer cells slightly increases during the entire duration of the simulation, the biggest cluster presenting a first fluctuating growing period followed by a faster increase accompanied by more regular oscillations. The behavior of the size of the biggest cluster is close to the one of the mean number Nc of cancer cells, which shows that the tumor swarms little, remains localized and grows around the initial cluster of cancer cells. The size of the biggest cluster of cancer cells does not exceed 600 boxes and remains smaller than system size equal to 30×30 cells, making explicit that the oscillating behavior shown in Figures 3-5 is not affected by small size effects and boundary conditions. The dependence of the pseudo-period of oscillations on system size is checked in the next subsection.
The first increase of the number of clusters of immune system cells reveals the fragmentation of the big cluster initially occupying the entire system. Then the number of clusters of immune system cells decreases while the size of the biggest cluster remains constant in average, proving that clusters disappear due to immune system cells destruction and not to coalescence. As shown in the snapshot of Figure 4, the maxima of Ni correspond to the periodic appearance of a small cluster containing a large number of immune system cells in the central disc occupied by the large killing population of cancer cells, in which the immune system cells will be soon consumed before a new cycle.
Qualitatively different behaviors are observed as some relevant parameters of the model vary. In a homogeneous system, the field E controlling activity fluctuations and one of the rate constants, kic or kci, related to immune system cell or cancer cell production proved to have a significant impact on the phenomena [8]. In the case of an inhomogeneous system, we perform systematic analyses for increasing values of either E or kic with the other parameters set to the values shown in the caption of Figure 1.
Figures 7 and 8 illustrate typical evolutions of the numbers of cells and mean activities, respectively, as the field E varies for kic=0.05, chosen smaller than kci. The same vertical scale is chosen for all subfigures of a given figure. For a sufficiently small value of the field, typically E≲10−5, corresponding to an inefficient control of activity fluctuations, a fast monotonous convergence toward elimination of the immune system and cancer growth is obtained. As E increases, the evolutions of both cell numbers and activities become oscillating, in a chaotic manner for E≃5×10−5 and a more regular way illustrated in Figures 3-5 for E≃10−4. For a critical value of the field, E≃1.4×10−4, associated with a bifurcation, irregular fluctuations appear after a long induction period and the behavior illustrates the 3Es, with the final escape of cancer from immunosurveillance. For a slightly larger value of the field, E≃1.5×10−4, the behavior is completely different and the system monotonically converges toward a steady state in which the activities of the two cell types are equal: Cancer is controlled, exactly as in the case of a sufficiently well thermalized homogeneous system [8]. The choice kic>kci implies that the final steady number of immune system cells is larger than the final steady number of cancer cells.
The results shown in Figures 7 and 8 correspond to simulations stopped at different times such that the immune system cells have disappeared, the number of cells in a spatial box exceeds a threshold set to 3000 as obtained for the critical value E≃1.4×10−4, or when the mean activities of the cancer cells or the immune system cells are equal. The two first cases obtained for E≤1.4×10−4 correspond to cancer escape. The last case observed for E≥1.5×10−4 is associated with cancer control and the simulation time is then supposed to mimic the life expectancy of a patient dying of a disease other than cancer. The characteristic times of the evolution are completely different depending on the field values, in particular because the 3Es are associated with a bifurcation, i.e. a divergence of the induction time before explosion [32]. Sufficiently far from the critical field value, the relaxation toward the stationary state is much faster.
Figure 9 represents the variation of the mean pseudo-period of oscillations of the total number of immune system cells Ni versus E in the field range where they are observed. The uncertainty about their determination is minimum for E≃10−4 associated with the most regular oscillations of both Ni, ui, and uc and increases as the boundaries of the domain of existence get closer. Even for E≃10−4 the time between two maxima of Ni decreases over the evolution, leading to rather large error bars for the pseudo-period. The averaged pseudo-period ⟨T⟩ slightly decreases as E increases. The value of ⟨T⟩ in a system of 30×30 boxes is satisfactorily compared to the result obtained in a system of 20×20 boxes. The difference between the two results is much smaller than the variability of the pseudo-period during time, proving that the simulation results deduced from a 30×30 system or even a 20×20 system are not sensitive to boundary effects and reliable.
The different typical behaviors encountered as the rate constant kic controlling the autocatalytic formation of immune system cells increases are represented in Figures 10 and 11, all other parameters being fixed.
As shown in Figure 11 for sufficiently small values of kic obeying kic≲0.007, the system rapidly and monotonously converges to a steady solution with the same final value uendc=uendi of the mean activities of the two cell types c and i. The final cell numbers obey Nendc>Nendi for kic<kci, they coincide for kic=kci, and obey Nendc<Nendi for kci<kic≲0.007, as shown in Figure 10. The symmetric roles played by the constants kic and kci clearly appear. In all these cases, cancer is controlled. On the contrary, for sufficiently large kic, such that kic≳0.01, cancer escapes from the surveillance of the immune system. The case obtained for kic=0.01, already shown in Figure 1, illustrates the 3Es, with a long induction period and the eventual explosion of cancer cell number and both activities. For kic≃0.05, rather regular oscillations of cell numbers and activities are obtained. Before disappearing for larger values of kic, the oscillations become erratic for kic≃0.25. The value kic=0.01 can be considered as critical, smaller values leading to monotonous cancer control and similar mean activities of both cell types, larger values to oscillating cancer escape and larger mean activity for cancer cells. The results presented in Figures 10 and 11 are nonintuitive since increasing kic, which should reinforce the production of immune system cells, induces the exit from the domain where cancer is controlled. This result is explained by the feedback effect associated with the activity-dependence of the rate constant κic given in Eq. (2.5). At constant field E, increasing kic locally depletes the population of cancer cells from its most active representatives, which induces a local increase of cancer cell activity. This fluctuation of cancer activity is then responsible for cancer cell production through the autocatalytic process given in Eq. (2.3), which explains the eventual explosion of the number of cancer cells and cancer escape from immunosurveillance. In other words, a larger value of kic leads to a larger critical value Ec of the field associated with the 3Es and delimiting the domains between cancer escape and cancer control: As shown in Figures 7 and 8, for kic=0.05 we find Ec≃1.4×10−4 and in Figures 10 and 11, for kic=0.01 the bifurcation arises for a smaller critical field value Ec≃10−4. A stronger stimulation of the immune system requires a more efficient regulation of activity fluctuations for cancer to be controlled.
Figure 12 shows the variation of the mean pseudo-period of oscillations ⟨T⟩ versus the rate constant kic in the range in which they exist. The uncertainty is much smaller than the variation of the averaged pseudo-period and we conclude without ambiguity that ⟨T⟩ increases as kic increases. It was not easy to anticipate such a result, due, once again, to the antagonistic effects included in the model. At first glance, increasing kic should lead to the faster production of immune system cells according to the autocatalytic process given in Eq. (2.2), which should decrease the time between two maxima of Ni or ui. However, as shown in Figures 4 and 5, the formation of a cluster of immune system cells associated with a maximum of mean activity ui is rapidly followed by a maximum for the mean activity uc of the cancer cells, due to the depletion of small activity cancer cells induced by the process given in Eq. (2.2). Then, the very active cancer cells consume the immune system cells according to the autocatalytic process given in Eq. (2.3). As a result, the feedback effect on the activities is stronger than the increase of kic and the time between two maxima of Ni or ui is longer for larger kic.
We end the sensitivity analysis of the system to the field E and the rate constant kic with a few remarks about the effect of the other parameters. As explicitly shown in the case of a homogeneous system [8], the role of kci is symmetrical to the one of kic. The results displayed in Figures 10 and 11 for kic<kci, kic=kci, and kic>kci illustrate this property. A small change of cell speed v modifies the critical values of the field and the constant kic for which the 3Es are observed but does not qualitatively alter the behavior of the system. However, another issue, far beyond the scope of this article, could be addressed in a different domain of parameter space. The model could be extended to the study of the anisotropic growth of a tumor by choosing smaller interaction angles θinter, larger cell speeds v, smaller rate constants kv of velocity randomisation and smaller deflection angles θjump in order to favor the formation of clusters of cancer cells with a high activity moving in the same direction during a sufficiently long time to distort the tumor.
The interactions between cancer and the immune system are complex and nontrivial antagonistic effects have been observed in patients with long-standing tumors and revealed by oscillations of some markers of inflammation [21]. We extend a model of cancer and immune system competition [7,8] to an inhomogeneous system in order to follow the dynamics of tumor growth at cell scale within the framework of thermostatted kinetic theory. The model includes cell-cell interactions mimicking mutation, death, division, regulation, and learning associated with increased activity. Activity fluctuations are thermostatted to account for dissipation of information through nonspecific cell interactions. We implement a kinetic Monte Carlo algorithm to directly simulate the kinetic equations in the spirit of DSMC method [39]. The impact of thermalization on the behavior of the system is important and a bifurcation occurs as the field controlling activity fluctuations increases. For small field values, cancer invades the entire system whereas, for large field values, the system reaches a steady state with non vanishing numbers of cancer cells and immune system cells. Cancer is then controlled and the tumor reaches a stationary size. For a critical value of the field for which large activity fluctuations are eventually generated, the model reproduces the three Es of immunoediting [19], leading to the unexpected escape of cancer from immune system control after the elimination and equilibrium phases. For a field value slightly smaller than the critical value, the model is sufficiently complex to reproduce an observed feature of tumor growth associated with rather regular oscillations of the immune response.
We did not observe time oscillations in simulations of a homogeneous system [7,8]. Their stimulation requires non trivial initial conditions that are spontaneously generated in the simulations of the inhomogeneous system. In particular oscillations may begin when a cluster of immune system cells have a sufficiently higher activity than the neighboring cancer cells. The simulations have been performed for an efficient cell velocity randomization which ensures the rapid recovery of the initial symmetry of the tumor. Following the evolution of the mean number of clusters as well as the size of the biggest cluster reveals that the tumor remains a connected space apart from the loss of very small clusters on its periphery. The model involves two autocatalytic reactions, the one producing cancer cells, the other, immune system cells, with a regulation through activity-dependent rate constants. In the present case, large-activity immune system cells kill low-activity cancer cells, leading to increased activity for the remaining cancer cells. These educated cancer cells are then able to kill low-activity immune system cells. Consequently, the mean activity of the remaining immune system cells increases and a new cycle may begin.
It is worth noting that decreasing the rate constant kic controlling the autocatalytic production of immune system cells has a similar effect on the behavior of the system as increasing the field E. The same succession of bifurcations from cancer proliferation associated with chaotic oscillations, periodic oscillations, long induction related to the 3Es, and cancer control is observed as shown when comparing Figures 7 and 10. Counterintuitively, larger rates of autocatalytic production of immune system cells require more efficient thermalization for cancer to be controlled, due to the feedback effect induced by activity-dependent rate constants.
To our knowledge, the results of the model we developed are the first to account for the 3Es of immunotherapy and oscillations of immune system response within the framework of thermostatted kinetic theory and, in particular, at cell scale. As a perspective, macroscopic equations for concentrations and activities could be derived from the kinetic equations in the limit of appropriate space and time scaling [5,6]. This derivation could give a microscopic interpretation at the scale of cell interactions to mean-field models [23,24,25,26] and to the delays [27,28] introduced in the equations for macroscopic quantities such as the concentrations of different cell types.
The results suggest improved treatment protocols, combining vaccination and chemotherapy at proper times. Reducing the dominance of regulatory T cells over killer T cells may be crucial in immunotherapy: Whereas killer T cells destroy cancer cells, regulatory T cells are susceptible to improperly reduce the number of immune system cells, due to misleading information conveyed by cancer cells. Although the model introduces a single type of immune system cells, it takes into account the role of killer T cells through the process (2.2) and the role of regulatory T cells through the process (2.3). Vaccination can be used to maximize killer T cell response and chemotherapy or radiotherapy, to deplete, inter alia, regulatory T cell population. Hence, the choice of the timing of vaccination, on the one hand, and chemotherapy, on the other hand, with respect to the immune regulatory cycle revealed for example by C-Reactive Protein oscillations, could be crucial in the treatment of some tumors. The model, which reproduces intrinsic oscillations could be used to minimize the averaged number of cancer cells after periodic series of immune system boosting followed by immune system cell depletion at different times of the cycle. The results could give some hints on how a tumor reacts to periodic double disturbances mimicking vaccination and chemotherapy.
All authors declare no conflicts of interest in this paper.
[1] | Barta DJ, Tibbitts TW, Bula RJ, et al. (1992) Evaluation of light-emitting diode characteristics for a space-based plant irradiation source. Adv Space Res 12: 141–149. |
[2] | Bula RJ, Morrow RC, Tibbitts TW, et al. (1991) Light-emitting diodes as a radiation source for plants. HortScience 26: 203–205. |
[3] | Emmerich JC, Morrow RC, Clavette T, et al. (2004) Plant research unit lighting system development. SAE Technical Paper Series: 2004-01-2454. |
[4] | Morrow RC, Duffie NA, Tibbitts TW, et al. (1995) Plant response in the ASTROCULTURE flight experiment unit. SAE Technical Paper Series: 951624. |
[5] | Massa GD, Kim HH, Wheeler RM, et al. (2008) Plant productivity in response to LED lighting. HortScience 43: 1951–1956. |
[6] | Inada K (1976) Action spectra for photosynthesis in higher plants. Plant Cell Physiol 17: 355–365. |
[7] | McCree KJ (1972) The action spectrum absorptance and quantum yield of photosynthesis in crop plants. Agric Meteorol 9: 191–216. |
[8] | Brown CS, Schuerger AC, Sager JC (1995) Growth and photomorphogenesis of pepper plants under red light-emitting diodes with supplemental blue or far-red lighting. J Am Soc Hortic Sci 120: 808–813. |
[9] | Goins GD, Yorio NC, Sanwo-Lewandowski MM, et al. (1998) Life cycle experiments with Arabidopsis grown under red light-emitting diodes (LEDs). Life Support Biosph Sci 5: 143–149. |
[10] | Hoenecke ME, Bula RJ, Tibbitts TW (1992) Importance of blue photon levels for lettuce seedlings grown under red-light-emitting diodes. HortScience 27: 427–430. |
[11] |
Tripathy BC, Brown CS (1995) Root-shoot interaction in the greening of wheat seedlings grown under red light. Plant Physiol 107: 407–411. doi: 10.1104/pp.107.2.407
![]() |
[12] | Yorio NC, Wheeler RM, Goins GD, et al. (1998) Blue light requirements for crop plants used in bioregenerative life support systems. Life Support Biosph Sci 5: 119–128. |
[13] |
Cope KR, Snowden MC, Bugbee B (2014) Photobiological interactions of blue light and photosynthetic photon flux: effects of monochromatic and broad-spectrum light sources. Photochem Photobiol 90: 574–584. doi: 10.1111/php.12233
![]() |
[14] |
Goins GD, Yorio NC, Sanwo MM, et al. (1997) Photomorphogenesis, photosynthesis, and seed yield of wheat plants grown under red light-emitting diodes (LEDs) with and without supplemental blue lighting. J Exp Bot 48: 1407–1413. doi: 10.1093/jxb/48.7.1407
![]() |
[15] |
Hogewoning SW, Douwstra P, Trouwborst G, et al. (2010) An artificial solar spectrum substantially alters plant development compared with usual climate room irradiance spectra. J Exp Bot 61: 1267–1276. doi: 10.1093/jxb/erq005
![]() |
[16] | Bugbee B (2016) Toward an optimal spectral quality for plant growth and development: the importance of radiation capture. Acta Hort 1134: 1–12. |
[17] | Assmann SM, Shimazaki K-I (1999) The multisensory guard cell. Stomatal responses to blue light and abscisic acid. Plant Physiol 119: 809–815. |
[18] |
Kaufman LS (1993) Transduction of blue-light signals. Plant Physiol 102: 333–337. doi: 10.1104/pp.102.2.333
![]() |
[19] |
Lin C, Ahmad M, Cashmore AR (1996) Arabidopsis cryptochrome 1 is a soluble protein mediating blue light-dependent regulation of plant growth and development. Plant J 10: 893–902. doi: 10.1046/j.1365-313X.1996.10050893.x
![]() |
[20] | Cope KR, Bugbee B (2013) Spectral effects of three types of white light-emitting diodes on plant growth and development: absolute versus relative amounts of blue light. HortScience 48: 504–509. |
[21] |
Poulet L, Massa GD, Morrow RC, et al. (2014) Significant reduction in energy for plant-growth lighting in space using targeted LED lighting and spectral manipulation. Life Sciences in Space Research 2: 43–53. doi: 10.1016/j.lssr.2014.06.002
![]() |
[22] |
Smith H, Whitelam GC (1997) The shade avoidance syndrome: multiple responses mediated by multiple phytochromes. Plant Cell Environ 20: 840–844. doi: 10.1046/j.1365-3040.1997.d01-104.x
![]() |
[23] |
Snowden MC, Cope KR, Bugbee B (2016) Sensitivity of seven diverse species to blue and green light: interactions with photon flux. Plos One 11: e0163121. doi: 10.1371/journal.pone.0163121
![]() |
[24] |
Wang XY, Xu XM, Cui J (2015) The importance of blue light for leaf area expansion, development of photosynthetic apparatus, and chloroplast ultrastructure of Cucumis sativus grown under weak light. Photosynthetica 53: 213–222. doi: 10.1007/s11099-015-0083-8
![]() |
[25] | Ouzounis T, Rosenqvist E, Ottosen CO (2015) Spectral effects of artificial light on plant physiology and secondary metabolism: a review. HortScience 50: 1128–1135. |
[26] | Mitchell CA, Stutte GW (2015) Sole-source lighting for controlled-environment agriculture. NASA Technical Reports. |
[27] |
Chen L, Lin CC, Yeh CW, et al. (2010) Light converting inorganic phosphors for white light-emitting diodes. Materials 3: 2172–2195. doi: 10.3390/ma3032172
![]() |
[28] |
Brodersen CR, Vogelmann TC (2010) Do changes in light direction affect absorption profiles in leaves? Funct Plant Biol 37: 403–412. doi: 10.1071/FP09262
![]() |
[29] |
Sun JD, Nishio JN, Vogelmann TC (1998) Green light drives CO2 fixation deep within leaves. Plant Cell Physiol 39: 1020–1026. doi: 10.1093/oxfordjournals.pcp.a029298
![]() |
[30] |
Terashima I, Fujita T, Inoue T, et al. (2009) Green light drives leaf photosynthesis more efficiently than red light in strong white light: revisiting the enigmatic question of why leaves are green. Plant Cell Physiol 50: 684–697. doi: 10.1093/pcp/pcp034
![]() |
[31] | Frantz JM, Joly RJ, Mitchell CA (2000) Intracanopy lighting influences radiation capture, productivity, and leaf senescence in cowpea canopies. J Am Soc Hortic Sci 125: 694–701. |
[32] | Kim HH, Goins GD, Wheeler RM, et al. (2004) Green-light supplementation for enhanced lettuce growth under red- and blue-light-emitting diodes. HortScience 39: 1617–1622. |
[33] |
Lu N, Maruo T, Johkan M, et al. (2012) Effects of supplemental lighting with light-emitting diodes (LEDs) on tomato yield and quality of single-truss tomato plants grown at high planting density. Environmental Control in Biology 50: 63–74. doi: 10.2525/ecb.50.63
![]() |
[34] |
Park Y, Runkle ES (2017) Far-red radiation promotes growth of seedlings by increasing leaf expansion and whole-plant net assimilation. Environ Exp Bot 136: 41–49. doi: 10.1016/j.envexpbot.2016.12.013
![]() |
[35] |
Zhen S, van Iersel MW (2017) Far-red light is needed for efficient photochemistry and photosynthesis. J Plant Physiol 209: 115–122. doi: 10.1016/j.jplph.2016.12.004
![]() |
[36] | Blom TJ, Tsujita MJ, Roberts GL (1995) Far-red at end of day and reduced irradiance affect plant height of easter and asiatic hybrid lilies. HortScience 30: 1009–1012. |
[37] | Chia PL, Kubota C (2010) End-of-day far-red light quality and dose requirements for tomato rootstock hypocotyl elongation. HortScience 45: 1501–1506. |
[38] | Decoteau DR, Friend HH (1991) Growth and subsequent yield of tomatoes following end-of-day light treatment of transplants. HortScience 26: 1528–1530. |
[39] |
Decoteau DR, Kasperbauer MJ, Daniels DD, et al. (1988) Plastic mulch color effects on reflected light and tomato plant-growth. Sci Hort 34: 169–175. doi: 10.1016/0304-4238(88)90089-1
![]() |
[40] | Ilias IF, Rajapakse N (2005) The effects of end-of-the-day red and far-red light on growth and flowering of Petunia × hybrida 'countdown burgundy' grown under photoselective films. HortScience 40: 131–133. |
[41] |
Kasperbauer MJ, Peaslee DE (1973) Morphology and photosynthetic efficiency of tobacco leaves that received end-of-day red or far-red light during development. Plant Physiol 52: 440–442. doi: 10.1104/pp.52.5.440
![]() |
[42] |
Yang ZC, Kubota C, Chia PL, et al. (2012) Effect of end-of-day far-red light from a movable LED fixture on squash rootstock hypocotyl elongation. Sci Hort 136: 81–86. doi: 10.1016/j.scienta.2011.12.023
![]() |
[43] | Fraszczak B (2013) Effect of short-term exposure to red and blue light on dill plants growth. HortScience 40: 177–185. |
[44] |
Chinchilla S, Izzo LG, van Santen E, et al. (2018) Growth and physiological responses of lettuce grown under pre-dawn or end-of-day sole-source light-quality treatments. horticulturae 4: 8. doi: 10.3390/horticulturae4020008
![]() |
[45] |
Jishi T, Kimura K, Matsuda R, et al. (2016) Effects of temporally shifted irradiation of blue and red LED light on cos lettuce growth and morphology. Sci Hort 198: 227–232. doi: 10.1016/j.scienta.2015.12.005
![]() |
[46] |
Sung IK, Takano T (1997) Effects of supplemental blue-and red-lights in the morning twilight on the growth and physiological responses of cucumber seedlings. Environmental Control in Biology 35: 261–265. doi: 10.2525/ecb1963.35.261
![]() |
[47] |
Goto E (2003) Effects of light quality on growth of crop plants under artificial lighting. Environmental Control in Biology 41: 121–132. doi: 10.2525/ecb1963.41.121
![]() |
[48] |
Tallman G, Zeiger E (1988) Light quality and osmoregulation in Vicia guard cells: evidence for involvement of three metabolic pathways. Plant Physiol 88: 887–895. doi: 10.1104/pp.88.3.887
![]() |
[49] | Gómez C, Mitchell CA (2015) Growth responses of tomato seedlings to different spectra of supplemental lighting. Hortscience 50: 112–118. |
[50] | Thomas B, Vince-Prue D (1996) Daylength perception in short-day plants. In: Photoperiodism in Plants, 2 Eds. London: Academic Press, 118–142. |
[51] | Hamaker CK (1998) Influence of photoperiod and temperature on flowering of Asclepias tuberosa, Campanula carpatica 'Blue Clips', Coreopsis grandiflora 'Early Sunrise', Coreopsis verticillata 'Moonbeam', Lavandula angustifolia 'Munstead', and Physostegia virginiana 'Alba'. MS Thesis, Michigan State University. |
[52] |
Runkle ES, Heins RD, Cameron AC, et al. (2001) Photocontrol of flowering and stem extension of the intermediate-day plant Echinacea purpurea. Physiol Plant 112: 433–441. doi: 10.1034/j.1399-3054.2001.1120318.x
![]() |
[53] | Craig DS (2012) Determining effective ratios of red and far-red light from light-emitting diodes that control flowering of photoperiodic ornamental crops. MS Thesis, Michigan State University. |
[54] |
Meng Q, Runkle ES (2015) Low-intensity blue light in night-interruption lighting does not influence flowering of herbaceous ornamentals. Sci Hort 186: 230–238. doi: 10.1016/j.scienta.2015.01.038
![]() |
[55] |
Meng Q, Runkle ES (2017) Moderate-intensity blue radiation can regulate flowering, but not extension growth, of several photoperiodic ornamental crops. Environ Exp Bot 134: 12–20. doi: 10.1016/j.envexpbot.2016.10.006
![]() |
[56] | Whitman CM, Heins RD, Cameron AC, et al. (1998) Lamp type and irradiance level for daylength extensions influence flowering of Campanula carpatica 'Blue clips', Coreopsis grandiflora 'Early Sunrise', and Coreopsis verticillata 'Moonbeam'. J Am Soc Hortic Sci 123: 802–807. |
[57] | Nadarajah N (2011) Is solid state lighting ready for the incandescent lamp phase-out? SPIE Optical Engineering + Applications Conference 8123: 812302. |
[58] |
Runkle ES, Padhye SR, Oh W, et al. (2012) Replacing incandescent lamps with compact fluorescent lamps may delay flowering. Sci Hort 143: 56–61. doi: 10.1016/j.scienta.2012.05.028
![]() |
[59] | Meng QW, Runkle ES (2014) Controlling flowering of photoperiodic ornamental crops with light-emitting diode lamps: a coordinated grower trial. HortTechnology 24: 702–711. |
[60] |
Craig DS, Runkle ES (2016) An intermediate phytochrome photoequilibria from night interruption lighting optimally promotes flowering of several long-day plants. Environ Exp Bot 121: 132–138. doi: 10.1016/j.envexpbot.2015.04.004
![]() |
[61] | Craig DS, Runkle ES (2013) A moderate to high red to far-red light ratio from light-emitting diodes controls flowering of short-day plants. J Am Soc Hortic Sci 138: 167–172. |
[62] | Meng QW, Runkle ES (2016) Control of flowering using night-interruption and day-extension LED lighting, In: Kozai T, Fujiwara K and Runkle ES, LED Lighting for Urban Agriculture. Springer, Singapore, 191–201. |
[63] |
Meng QW, Runkle ES (2017) Moderate-intensity blue radiation can regulate flowering, but not extension growth, of several photoperiodic ornamental crops. Environ Exp Bot 134: 12–20. doi: 10.1016/j.envexpbot.2016.10.006
![]() |
[64] | Hamamoto H, Yamazaki K (2009) Reproductive response of okra and native rosella to long-day treatment with red, blue, and green light-emitting diode lights. HortScience 5: 1494–1497. |
[65] |
Jeong SW, Park S, Jin JS, et al. (2012) Influences of four different light-emitting diode lights on flowering and polyphenol variations in the leaves of Chrysanthemum (Chrysanthemum morifolium). J Agric Food Chem 60: 9793–9800. doi: 10.1021/jf302272x
![]() |
[66] | Hamamoto H, Shimaji H, Higashide T (2003) Budding and bolting responses of horticultural plants to night-break treatments with LEDs of various colors. J Agric Meteorol 59: 103–110. |
[67] | Kohyama F, Whitman C, Runkle ES (2014) Comparing flowering responses of long-day plants under incandescent and two commercial light-emitting diode lamps. HortTechnology 24: 490–495. |
[68] | Lang SP, Struckmeyer BE, Tibbitts TW (1983) Morphology of intumescence development on tomato plants. J Am Soc Hortic Sci 108: 266–271. |
[69] |
Morrow RC, Tibbitts TW (1988) Evidence for involvement of phytochrome in tumor-development on plants. Plant Physiol 88: 1110–1114. doi: 10.1104/pp.88.4.1110
![]() |
[70] | Craver JC, Miller CT, Williams KA, et al. (2014) Ultraviolet radiation affects intumescence development in ornamental sweetpotato (Ipomoea batatas). HortScience 49: 1277–1283. |
[71] |
Kubota C, Eguchi T, Kroggel M (2017) UV-B radiation dose requirement for suppressing intumescence injury on tomato plants. Sci Hort 226: 366–371. doi: 10.1016/j.scienta.2017.09.006
![]() |
[72] | Rud NA (2009) Environmental factors influencing the physiological disorders of edema on ivy geranium (Pelargonium Peltatum) and intumescences on tomato (Solanum Lycopersicum). MS Thesis, Kansas State University. |
[73] | Eguchi T, Hernández R, Kubota C (2016) Far-red and blue light synergistically mitigate intumescence injury of tomato plants grown under UV-deficit light environment. HortScience 51: 712–719. |
[74] | Wollaeger HM, Runkle ES (2014) Growth of impatiens, petunia, salvia, and tomato seedlings under blue, green, and red light-emitting diodes. HortScience 49: 734–740. |
[75] | Wollaeger HM, Runkle ES (2015) Growth and acclimation of impatiens, salvia, petunia, and tomato seedlings to blue and red light. HortScience 50: 522–529. |
[76] | Wheeler RM, Morrow RC (2010) Physiological disorders in closed, controlled environment crops. NASA Technical Reports. |
[77] |
Trouwborst G, Hogewoning SW, van Kooten O, et al. (2016) Plasticity of photosynthesis after the 'red light syndrome' in cucumber. Environ Exp Bot 121: 75–82. doi: 10.1016/j.envexpbot.2015.05.002
![]() |
[78] |
Ohashi-Kaneko K, Takase M, Kon N, et al. (2007) Effect of light quality on growth and vegetable quality in leaf lettuce, spinach and komatsuna. Environmental Control in Biology 45: 189–198. doi: 10.2525/ecb.45.189
![]() |
[79] | Son KH, Oh MM (2013) Leaf shape, growth, and antioxidant phenolic compounds of two lettuce cultivars grown under various combinations of blue and red light-emitting diodes. HortScience 48: 988–995. |
[80] | Johansen NS, Vanninen I, Pinto DM, et al. (2011) In the light of new greenhouse technologies: 2. Direct effects of artificial lighting on arthropods and integrated pest management in greenhouse crops. Ann Appl Biol 159: 1–27. |
[81] | Schuerger AC, Brown CS (1997) Spectral quality affects disease development of three pathogens on hydroponically grown plants. HortScience 32: 96–100. |
[82] | Vanninen I, Pinto DM, Nissinen AI, et al. (2010) In the light of new greenhouse technologies: 1. Plant-mediated effects of artificial lighting on arthropods and tritrophic interactions. Ann Appl Biol 157: 393–414. |
[83] |
Carvalho SD, Folta KM (2014) Environmentally modified organisms – Expanding genetic potential with light. Crit Rev Plant Sci 33: 486–508. doi: 10.1080/07352689.2014.929929
![]() |
[84] | Kanto T, Matsuura K, Yamada M, et al. (2009) UV-B radiation for control of strawberry powdery mildew. VI International Strawberry Symposium 842: 359–362. |
[85] |
Kobayashi M, Kanto T, Fujikawa T, et al. (2014) Supplemental UV radiation controls rose powdery mildew disease under the greenhouse conditions. Environmental Control in Biology 51: 157–163. doi: 10.2525/ecb.51.157
![]() |
[86] |
Imada K, Tanaka S, Ibaraki Y, et al. (2014) Antifungal effect of 405-nm light on Botrytis cinerea. Lett Appl Microbiol 59: 670–676. doi: 10.1111/lam.12330
![]() |
[87] |
Tokuno A, Ibaraki Y, Ito S-i, et al. (2012) Disease suppression in greenhouse tomato by supplementary lighting with 405 nm LED. Environmental Control in Biology 50: 19–29. doi: 10.2525/ecb.50.19
![]() |
[88] | Kudo R, Ishida Y, Yamamoto K (2011) Effects of green light irradiation on induction of disease resistance in plants. VI International Symposium on Light in Horticulture 907: 251–254. |
[89] | Kudo R, Yamamoto K (2013) Effects of green light irradiation on Corynespora leaf spot disease in Perilla. Hortic Res 12: 13–157. |
[90] | Kudo R, Yamamoto K, Suekane A, et al. (2009) Development of green light pest control systems in plants. I. Studies on effects of green light irradiation on induction of disease resistance. SRI Res Rep 93: 31–35. |
[91] |
Suthaparan A, Torre S, Stensvand A, et al. (2010) Specific light-emitting diodes can suppress sporulation of Podosphaera pannosa on greenhouse roses. Plant Dis 94: 1105–1110. doi: 10.1094/PDIS-94-9-1105
![]() |
[92] | Wang H, Jiang YP, Yu HJ, et al. (2010) Light quality affects incidence of powdery mildew, expression of defence-related genes and associated metabolism in cucumber plants. Eur J Plant Pathol 127: 125–135. |
[93] |
Patel JS, Zhang SA, McGrath MT (2016) Red light increases suppression of downy mildew in basil by chemical and organic Products. J Phytopathol 164: 1022–1029. doi: 10.1111/jph.12523
![]() |
[94] | Bishop AL, Bellis GA, McKenzie HJ, et al. (2006) Light trapping of biting midges Culicoides spp. (Diptera: Ceratopogonidae) with green light-emitting diodes. Australian Journal of Entomology 45: 202–205. |
[95] | Chen TY, Chu CC, Henneberry TJ, et al. (2004) Monitoring and trapping insects on poinsettia with yellow sticky card traps equipped with light-emitting diodes. HortTechnology 14: 337–341. |
[96] |
Duehl AL, Cohnstaedt AR, Teal P (2011) Evaluating light attraction to increase trap efficiency for Tribolium castaneum (Coleoptera: Tenebrionidae). J Econ Entomol 104: 1430–1435. doi: 10.1603/EC10458
![]() |
[97] |
Katsuki M, Omae Y, Okada K, et al. (2012) Ultraviolet light-emitting diode (UV LED) trap for the West Indian sweet potato weevil, Euscepes postfasciatus (Coleoptera: Curculionidae). Appl Entomol Zool 47: 285–290. doi: 10.1007/s13355-012-0113-y
![]() |
[98] | McQuate GT (2014) Green light synergistally enhances male sweetpotato weevil response to sex pheromone. Sci Rep 4: 4499. |
[99] |
Sonoda S, Kataoka Y, Kohara Y, et al. (2014) Trap catches of dipteran insects using ultraviolet LED (light emitting diode) and water-pan trap. Jpn J Appl Entomol 58: 32–35. doi: 10.1303/jjaez.2014.32
![]() |
[100] |
Stukenberg N, Ahrens N, Poehling HM (2018) Visual orientation of the black fungus gnat, Bradysia difformis, explored using LEDs. Entomol Exp Appl 166: 113–123. doi: 10.1111/eea.12652
![]() |
[101] |
Stukenberg N, Gebauer K, Poehling HM (2015) Light emitting diode(LED)-based trapping of the greenhouse whitefly (Trialeurodes vaporariorum). J Appl Entomol 139: 268–279. doi: 10.1111/jen.12172
![]() |
[102] |
Yoon J-b, Nomura M, Ishikura S (2012) Analysis of the flight activity of the cotton bollworm Helicoverpa armigera (Hubner) (Lepidoptera: Noctuidae) under yellow LED lighting. Jpn J Appl Entomol 56: 103–110. doi: 10.1303/jjaez.2012.103
![]() |
[103] |
Liu H, Gao Z, Deng SZ, et al. (2018) The photokinesis of oriental fruit flies, Bactrocera dorsalis, to LED lights of various wavelengths. Entomol Exp Appl 166: 102–112. doi: 10.1111/eea.12648
![]() |
[104] | Nelson JA, Bugbee B (2014) Economic analysis of greenhouse lighting: light emitting diodes vs. high intensity discharge fixtures. Plos One 9: e99010. |
[105] | Wallace C, Both AJ (2016) Evaluating operating characteristics of light sources for horticultural applications. VIII International Symposium on Light in Horticulture 1134: 435–443. |
[106] | Frantz JM, Chun C, Joly RJ, et al. (1998) Intracanopy lighting of cowpea canopies in controlled environments. Life Support Biosph Sci 5: 183–190. |
[107] | Massa GD, Mitchell CA, Emmerich JC, et al. (2005) Development of a reconfigurable LED plant-growth lighting system for equivalent system mass reduction in an ALS. SAE Technical Paper Series: 2005–01–2955. |
[108] | Dueck TA, Janse J, Eveleens BA, et al. (2012) Growth of tomatoes under hybrid LED and HPS lighting systems. Acta Hort 925: 335–342. |
[109] | Gómez C, Mitchell CA (2016) Physiological and productivity responses of high-wire tomato as affected by supplemental light source and distribution within the canopy. J Am Soc Hortic Sci 141: 196–208. |
[110] |
Pettersen RI, Torre S, Gislerod HR (2010) Effects of intracanopy lighting on photosynthetic characteristics in cucumber. Sci Hort 125: 77–81. doi: 10.1016/j.scienta.2010.02.006
![]() |
[111] |
Trouwborst G, Schapendonk A, Rappoldt K, et al. (2011) The effect of intracanopy lighting on cucumber fruit yield-Model analysis. Sci Hort 129: 273–278. doi: 10.1016/j.scienta.2011.03.042
![]() |
[112] | Gómez C, Mitchell CA (2014) Supplemental lighting for greenhouse-grown tomatoes: intracanopy LED towers vs. overhead HPS lamps. Acta Hort 1037: 855–862. |
[113] | Brault D, Gueymard C, Boily R, et al. (1989) Contribution of HPS lighting to the heating requirements of a greenhouse. Am Soc Agric Engr 89: 4039. |
[114] | van Iersel MW (2017) Optimizing LED lighting in controlled environment agriculture, In: Gupta SD, Light Emitting Diodes for Agriculture: Smart Lighting. Springer Nature Singapore Pte Ltd. 59–80. |
[115] | Pinho P, Hytönen T, Rantanen M, et al. (2012) Dynamic control of supplemental lighting intensity in a greenhouse environment. Lighting Res Technol 45: 295–304. |
[116] | Clausen A, Maersk-Moeller HM, Corfixen Soerensen J, et al. (2015) Integrating commercial greenhouses in the smart grid with demand response based control of supplemental lighting. International Conference on Industrial Technology Management Science: 199–213. |
[117] |
Schwend T, Beck M, Prucker D, et al. (2016) Test of a PAR sensor-based, dynamic regulation of LED lighting in greenhouse cultivation of Helianthus annuus. Eur J Hortic Sci 81: 152–156. doi: 10.17660/eJHS.2016/81.3.3
![]() |
[118] |
van Iersel MW, Gianino D (2017) An adaptive control approach for light-emitting diode lights can reduce the energy costs of supplemental lighting in greenhouses. HortScience 52: 72–77. doi: 10.21273/HORTSCI11385-16
![]() |
[119] | van Iersel MW, Dove S (2016) Maintaining minimum light levels with LEDs results in more energy-efficient growth stimulation of begonia. 2016 Conference program, American Society for Horticultural Science. |
[120] | van Iersel MW, Mattos E, Weavers G, et al. (2016) Using chlorophyll fluorescence to control lighting in controlled environment agriculture. VIII International Symposium on Light in Horticulture 1134: 427–433. |
[121] | van Iersel MW, Weaver G, Martin MT, et al. (2016) A Chlorophyll fluorescence-based biofeedback system to control photosynthetic lighting in controlled environment agriculture. J Am Soc Hortic Sci 141: 169–176. |
[122] |
Carstensen AM, Pocock T, Bankestad D, et al. (2016) Remote detection of light tolerance in basil through frequency and transient analysis of light induced fluorescence. Comput Electron Agric 127: 289–301. doi: 10.1016/j.compag.2016.06.002
![]() |
[123] |
Kozai T (2013) Resource use efficiency of closed plant production system with artificial light: concept, estimation and application to plant factory. Proc Jpn Acad Ser B 89: 447–461. doi: 10.2183/pjab.89.447
![]() |
[124] |
Al-Kodmany K (2018) The vertical farm: a review of developments and implications for the vertical city. Buildings 8: 24. doi: 10.3390/buildings8020024
![]() |
[125] | Kozai T, Fujiwara K, Runkle ES (2016) LED lighting for urban agriculture. Springer Nature Singapore Pte Ltd. |
[126] | Akiyama T, Kozai T (2016) Light environment in the cultivation space of plant factory with LEDs, In: Kozai T, Fujiwara K and Runkle ES, LED Lighting for Urban Agriculture. Springer Nature Singapore Pte Ltd. 91–109. |
[127] | Ibaraki Y (2016) Lighting efficiency in plant production under artificial lighting and plant growth modeling for evaluating the lighting efficiency, In: Kozai T, Fujiwara K and Runkle ES , LED Lighting for Urban Agriculture. Springer Nature Singapore Pte Ltd. 151–161. |
[128] | Hayashi E (2016) Current status of commercial plant factories with LED lighting market in Asia, Europe, and other regions, In: Kozai T, Fujiwara K and Runkle ES, LED Lighting for Urban Agriculture. Springer Nature Singapore Pte Ltd. 295–308. |
[129] | Higgins C (2016) Current status of commercial vertical farms with LED lighting market in North America, In: Kozai T, Fujiwara K and Runkle ES, LED Lighting for Urban Agriculture. Springer Nature Singapore Pte Ltd. 309–315. |
[130] | Hayashi E, Higgins C (2016) Global LED lighting players, economic analysis, and market creation for PFALs, In: Kozai T, Fujiwara K and Runkle ES, LED Lighting for Urban Agriculture. Springer Nature Singapore Pte Ltd. 317–345. |
[131] | Coyle BD, Ellison B (2017) Will consumers find vertically farmed produce "out of reach"? Choices 32: 1–8. |
[132] | Yano Y, Nakamura T, Maruyama A (2016) Consumer perception and understanding of vegetables produced at plant factories with artificial lighting, In: Kozai T, Fujiwara K and Runkle ES, LED Lighting for Urban Agriculture. Springer Nature Singapore Pte Ltd. 347–363. |
[133] | Wheeler RM (2010) Plants for human life support in space: From Myers to Mars. Grav Space Res 23: 25–35. |
[134] | Massa GD, Emmerich JC, Morrow RC, et al. (2006) Plant-growth lighting for space life support: a review. Gravitational and Space Biology 19: 19–30. |
[135] | Duffie NA, Zhou W, Morrow RC, et al. (1995) Control and monitoring of environmental parameters in the ASTROCULTURE flight experiment. SAE Technical Paper Series: 951627. |
[136] | Massa GD, Dufour NF, Craver JA, et al. (2017) VEG-01: Veggie hardware validation testing on the International Space Station. Open Agriculture 2: 33–41. |
[137] | Massa GD, Wheeler RM, Morrow RC, et al. (2016) Growth chambers on the International Space Station for large plants. Acta Hort 1134: 215–222. |
[138] | Mitchell CA, Dougher TAO, Nielsen SS, et al. (1996) Costs of providing edible biomass for a balanced vegetarian diet in a controlled ecological life support system, In: Suge H, Plants in Space Biology. Tohoku University Press, Sendai, Japan, 245–254. |
[139] | Wheeler RM, Mackowiak CL, Stutte GW, et al. (1996) NASA's biomass production chamber: a testbed for bioregenerative life support studies. Adv Space Res 18: 215–224. |
[140] | Cuello JL (2002) Latest developments in artificial lighting technologies for bioregenerative space life support. Proceedings of the Fourth International ISHS Symposium on Artificial Lighting, 49–56. |
[141] |
Drysdale AE, Ewert MK, Hanford AJ (2003) Life support approaches for Mars missions. Adv Space Res 31: 51–61. doi: 10.1016/S0273-1177(02)00658-0
![]() |
[142] | Bourget CM (2008) An introduction to light-emitting diodes. HortScience 43: 1944–1946. |
[143] |
De Micco V, Aronne G (2008) Biometric anatomy of seedlings developed onboard of Foton-M2 in an automatic system supporting growth. Acta Astronaut 62: 505–513. doi: 10.1016/j.actaastro.2008.01.019
![]() |
[144] |
Kitaya Y, Kawai M, Tsuruyama J, et al. (2001) The effect of gravity on surface temperature and net photosynthetic rate of plant leaves. Adv Space Res 28: 659–664. doi: 10.1016/S0273-1177(01)00375-1
![]() |
[145] |
Kwon T, Sparks JA, Nakashima J, et al. (2015) Transcriptional response of Arabidopsis seedlings during spaceflight reveals peroxidase and cell wall remodeling genes associated with root hair development. Am J Bot 102: 21–35. doi: 10.3732/ajb.1400458
![]() |
[146] |
Levinskikh MA, Sychev VN, Derendyaeva TA, et al. (2000) Analysis of the spaceflight effects on growth and development of super dwarf wheat grown on the Space Station Mir. J Plant Physiol 156: 522–529. doi: 10.1016/S0176-1617(00)80168-6
![]() |
[147] |
Manzano AI, Matia I, Gonzalez-Camacho F, et al. (2009) Germination of Arabidopsis seed in Space and in simulated microgravity: alterations in root cell growth and proliferation. Microgravity Sci Tec 21: 293–297. doi: 10.1007/s12217-008-9099-z
![]() |
[148] |
Johnson CM, Subramanian A, Pattathil S, et al. (2017) Comparative transcriptomics indicate changes in cell wall organization and stress response in seedlings during spaceflight. Am J Bot 104: 1219–1231. doi: 10.3732/ajb.1700079
![]() |
[149] |
Millar KD, Kumar P, Correll MJ, et al. (2010) A novel phototropic response to red light is revealed in microgravity. New Phytol 186: 648–656. doi: 10.1111/j.1469-8137.2010.03211.x
![]() |
[150] |
Vandenbrink JP, Herranz R, Medina FJ, et al. (2016) A novel blue-light phototropic response is revealed in roots of Arabidopsis thaliana in microgravity. Planta 244: 1201–1215. doi: 10.1007/s00425-016-2581-8
![]() |
[151] |
Adkins E, Eapen S, Kaluwile F, et al. (2010) Off-grid energy services for the poor: introducing LED lighting in the Millennium Villages Project in Malawi. Energy Policy 38: 1087–1097. doi: 10.1016/j.enpol.2009.10.061
![]() |
[152] |
Pode R (2010) Solution to enhance the acceptability of solar-powered LED lighting technology. Renew Sust Energ Rev 14: 1096–1103. doi: 10.1016/j.rser.2009.10.006
![]() |
1. | Carlo Bianca, On the Modeling of Energy-Multisource Networks by the Thermostatted Kinetic Theory Approach: A Review with Research Perspectives, 2022, 15, 1996-1073, 7825, 10.3390/en15217825 | |
2. | G. Morgado, L. Masurel, A. Lemarchand, C. Bianca, 2022, Chapter 12, 978-3-031-12514-0, 225, 10.1007/978-3-031-12515-7_12 | |
3. | Carlo Bianca, A decade of thermostatted kinetic theory models for complex active matter living systems, 2024, 50, 15710645, 72, 10.1016/j.plrev.2024.06.015 | |
4. | Carlo Bianca, Thermostatted Kinetic Theory Structures in Biophysics: Generalizations and Perspectives, 2024, 4, 2673-9909, 1278, 10.3390/appliedmath4040069 | |
5. | Gabriel Morgado, Annie Lemarchand, Carlo Bianca, From Cell–Cell Interaction to Stochastic and Deterministic Descriptions of a Cancer–Immune System Competition Model, 2023, 11, 2227-7390, 2188, 10.3390/math11092188 |