 Commentary
 Open Access
 Published:
The prospective mathematical idea satisfying both radiation hormesis under low radiation doses and linear nonthreshold theory under high radiation doses
Genes and Environment volume 42, Article number: 4 (2020)
Abstract
It has yet to be determined whether or not the probability of developing cancer due to radiation exposure levels of low doses is proportional to the dose. Herein, for radiation hormesis occurring at low doses, mathematical models using functions that take a mountainlike shape having two inflection points are considered. The following perspectives were obtained: (i) When the probability of developing cancer decreases at radiation levels above the natural background dose, the radiation hormesis effect occurs up to ~ 12.4 mSv. (ii) When there is a proportional relationship at ≥750 mSv, the radiation hormesis effect occurs up to ~ 225 mSv. Thus, by performing studies at the molecular and cellular levels for radiation doses at ≤16.8 or 307 mSv, it is possible to investigate carcinogenesis resulting from low radiation doses.
Background
As radiation has the ability to ionize substances, the idea that, as far as is possible, it is best to avoid any exposure is generally accepted. As the basis for this, the International Commission on Radiological Protection (ICRP) employs a model whereby the probability of developing cancer is proportional to the radiation exposure dose (Linear nonthreshold theory: LNT) [1]. Certainly, a proportional relationship has been described at ≥100 mSv, but it has not been confirmed whether there is a proportional relationship for cases of < 100 mSv [1]. In addition, the radiation hormesis theory is thought to actually have beneficial effects on health [2,3,4,5,6,7,8,9,10]. Therefore, I wondered if there was an idea that LNT and hormesis could hold at the same time. In this paper, such a mathematical idea is considered.
Probability of developing cancer D(x) and the inhibition effect R(x)
Hereafter, the radiation dose is defined as x, and the probability of developing cancer as D(x). Accepting the fact that at high doses, the probability of developing cancer is proportional to the dose, and taking the constant of proportionality to be k, Eq. 1 holds true. Here, for a hormesis region to be present at low doses, an inhibitor factor that reduces D(x) becomes necessary. The inhibition effect relating to the inhibitor factor is described by R(x). Taking this inhibition effect into account, D(x) may be defined as in Eq. 2 [11]. Here, it is assumed that both D(x) and R(x) are continuous functions. It should be noted that D(x) ≥ 0 is always satisfied, since there is no possibility that the probability of developing cancer becomes negative.
For actual in vivo cases, there exist many inhibitor factors including DNA repair, removal of active oxygen, and apoptosis. However, compared to the inhibitor factor that has the greatest effect on D(x), the effects of other inhibitor factors on D(x) are small or nonexistent. Thus, in the following, for the sake of simplicity, a single inhibitor factor to be present is considered.
Next, the form of R(x) is considered. Assuming the inhibitor factor is triggered only by radiation, when x is 0 mSv the value of R must be 0 (R(0) = 0). As the radiation dose increases, R(x) increases, but if the radiation increases excessively, the inhibitor factor itself becomes inactivated by the radiation, and thus R(x) begins to decrease at a certain dose. As a result, as x approaches infinity, the value of R becomes 0. Two forms of graphs for R(x) having these characteristics can be imagined, and are shown in Fig. 1. The radiation dose at which R(x) reaches a maximum is defined as x_{1}. The graph in Fig. 1a has a single point of inflection with x > x_{1}. On the other hand, the graph in Fig. 1b has a single point of inflection in each of the regions 0 < x < x_{1} and x > x_{1}.
When D(x) has a hormesis region for x > 0, the graph of D(x) has the form shown in Fig. 2a. For the case where R(x) takes the form in Fig. 1a, D(x) takes the form in Fig. 2b or c. Fig. 2b clearly has no hormesis region. Furthermore, although Fig. 2c has a hormesis region, it does not fulfill the condition that D(x) ≥ 0. Thus, in order for D(x) to have the form shown in Fig. 2a, it is necessary for R(x) to have the form shown in Fig. 1b.
In this paper, for R(x) to have the form shown in Fig. 1b, i chose the simple function given in Eq. 3 among several functions. The constant a is understandably positive. The value of x_{1} in Eq. 3 is 2/a, so x_{1} depends solely on the constant a. By varying a, the radiation dose x_{1} where the inhibitor amount reaches a maximum can be freely changed, and R(x) can be freely adjusted with respect to D(x).
Substituting Eq. 3 into Eq. 2, Eq. 4 is obtained. D(x) is a continuous function, and when it has the maximum hormesis region based on fulfillment of the condition D(x) ≥ 0, it has the form seen in Fig. 2d. The value of x for the local minimum of D(x) is defined as x_{2}, for which Eqs. 5 and 6 must be satisfied. Then, x_{2} = 1/a and ka = 1/e are obtained.
The value of x for the local maximum of D(x) is defined as x_{3}, and the other value of x where D has the same value as the local maximum D(x_{3}) is defined as x_{4}. Previously, x_{4} has been defined as being the zero equivalent point (ZEP) [3]. As x_{3} and x_{4} cannot be solved analytically, approximate values were obtained through numerical calculations using the graphing calculator “Grapher 2.5” as x_{3} = ~ 0.285/a and x_{4} = ~ 1.469/a. It is clear that x_{2}, x_{3}, and x_{4} do not depend on k, and x_{4} is ~ 5.15 times greater than x_{3}.
The radiation hormesis effect posits that the probability of developing cancer decreases at radiation levels above the natural background dose. Therefore, up to a certain dose above the natural background, D(x) should decrease. Thus, the natural background radiation dose can be taken to be between x_{3} and x_{2}. In order to maximize the hormesis region, x_{3} is set to the natural background dose, and then ZEP extends to up to ~ 5.15 times the value of the natural background dose. In addition, it is noted that the multiple, 5.15, is independent of the constant a. When the worldwide average dose of the natural background radiation is taken to be 2.4 mSv [12], the maximum ZEP is ~ 12.4 mSv.
Furthermore, as 2.4 mSv is x_{3} = ~ 0.285/a, the value of x corresponding to 100 mSv (x_{5}) becomes ~ 11.875/a. Here, it is considered whether Eq. 4 can be approximated to Eq. 1 at 100 mSv. Changing the form of Eq. 4 into that of Eq. 7 and substituting ka = 1/e, Eq. 8 is obtained. Furthermore, by substituting x_{5} = ~ 11.875/a, Eq. 9 is obtained. At 100 mSv and above, Eq. 4 satisfied by ka = 1/e can be approximated to Eq. 1. Thus, taking Eq. 4 as the model, linearity is satisfied above 100 mSv.
Lastly, when ZEP is ~ 12.4 mSv, x_{1} for the maximum of R(x) becomes ~ 16.8 mSv. That is, the inhibitor factor of which the amount is a maximum at ~ 16.8 mSv may show the maximum hormesis effect.
Maximum hormesis region when the x ^{2} term in Eq. 3 is replaced with an x ^{n} term
If the x^{2} term in Eq. 3 is replaced with an x^{n} term, a graph of the form of Fig. 1b can still be achieved. That is, Eq. 10 is expressed in place of Eq. 4. Taking the case of the maximum hormesis region, Eqs. 5 and 6 must be satisfied. Therefore, x_{2} = (n–1)/a and ka^{n–1} = {(n–1)/e}^{n–1} are obtained.
Since the case when n = 2 has already been considered in the second section, the cases where n ≥ 3 will be considered in sequence. As these cases are also impossible to solve analytically, the solutions are determined using numerical calculations. For n = 3, 4, 5, and 10, the maximum ZEP is ~ 10.1, ~ 8.9, ~ 8.1, and ~ 6.3 mSv, respectively. That is to say, as n increases, the maximum ZEP becomes smaller. Thus, Eq. 4 (Eq. 10 when n = 2) is best suited to consider the maximum hormesis region.
Condition for a hormesis region to be present for Eq. 4
The condition for ka giving the maximum hormesis region was considered in the second section. In contrast, in this section, the condition for ka under which the hormesis region begins to appear (Fig. 3) is determined.
The condition for which the hormesis region begins to appear is given by Eqs. 11 and 12, and then \( {x}_2=\left(2\sqrt{2}\right)/\mathrm{a} \) and \( ka=2\left(\sqrt{2}1\right){e}^{\sqrt{2}2} \) were determined. Therefore, when combined with the conclusions of the second section, the condition for which the hormesis region appears is 1/e ≤ ka < \( 2\left(\sqrt{2}1\right){e}^{\sqrt{2}2} \). Expressed to three significant digits, this corresponds to ~ 0.368 ≤ ka < ~ 0.461. Thus, for Eq. 4 to have a hormesis region, the restrictive condition must be fulfilled.
Reconsidering the second section when there is a proportional relationship at ≥750 mSv
Siegel et al. asserted that threshold is 0.75 Gy [13]. Reconsidering the assumption that x_{3} is the natural background radiation dose in the second section 2, the condition by which a proportional relationship should approximately hold at ≥750 mSv and accepting a ≤ 10% error is imposed. Using Eq. 8, it is necessary for Eq. 13 to be satisfied. Thus, x_{6} corresponding to 750 mSv was determined to be ~ 4.890/a. From x_{6}, it was determined that x_{3} = ~ 0.285/a and x_{4} = ~ 1.469/a correspond respectively to ~ 43.7 and ~ 225 mSv. Summarizing the results of the above calculations, when satisfying the proportional relationship with an error within 10% at ≥750 mSv, the maximum hormesis region becomes 43.7–225 mSv. In addition, x_{1} for the maximum of R(x) becomes ~ 307 mSv.
Conclusion and implication
When the probability of developing cancer decreases at radiation levels above the natural background dose, the maximum ZEP becomes ~ 12.4 mSv, and at the same time, a proportional relationship is approximately obtained at ≥100 mSv. At ~ 16.8 mSv, R(x) reaches a maximum. Additionally, for Eq. 4, a hormesis region appears when ~ 0.368 ≤ ka < ~ 0.461.
When there is a proportional relationship at ≥750 mSv, the maximum ZEP becomes ~ 225 mSv. At ~ 307 mSv, R(x) reaches a maximum.
Since statistically measuring D(x) at low doses is effectively not possible, analyzing the following three points would help clarify the radiation hormesis effect, perhaps making it possible to determine the probability of developing cancer at low doses.
 (i)
Finding a factor which expressed inhibition effect versus dose has the approximate form of Fig. 1b.
 (ii)
Analyzing the variations of the inhibitor factor in the region up to ~ 16.8 or 307 mSv.
 (iii)
Determining k, which indicates the correlation between D(x) and R(x).
Although preliminary, it is felt that the results and discussions presented in this paper may be of potential use to other researchers. Furthermore, if such inhibition factors are identified, it might possibly lead to a method of effectively reducing the cancer rates.
Availability of data and materials
All data generated during this study are included in this published article.
Abbreviations
 Eq:

Equation
 ICRP:

International Commission on Radiological Protection
 LNT:

Linear nonthreshold theory
 ZEP:

Zero equivalent point
References
 1.
ICRP. Lowdose extrapolation of radiationrelated cancer risk. ICRP publication 99. International commission on radiological protection. Oxford: Elsevier; 2006.
 2.
Socol Y, Dobrzynski L, Doss M, Feinendegen LE, Janiak MK, Miller ML, Sanders CL, Scott BR, Ulsh B, Vaiserman A. Ethical issues of current healthprotection policies on lowdose ionizing radiation. Dose Response. 2014;12:342–8.
 3.
Luckey TD. Radiation hormesis. Boca Raton: CRC Press, Inc.; 1991.
 4.
Sanders CL. Radiation hormesis and the linearnothreshold assumption. Berlin: Springer; 2010.
 5.
Feinendegen LE, Bond VP, Sondhaus CA. The dual response to lowdose irradiation: induction vs. prevention of DNA damage. In: Yamada T, Mothersill C, Michael BD, Potten CS. Biological effects of low dose radiation. Excerpta Medica. International Congress Series 1211. Amsterdam: Elsevier; 2000, p. 3–17.
 6.
Sutou S. A message to Fukushima: nothing to fear but fear itself. Genes Environ. 2016;38:12.
 7.
Sutou S. The 10th anniversary of the publication of genes and environment: memoir of establishing the Japanese environmental mutagen society and a proposal for a new collaborative study on mutagenic hormesis. Genes Environ. 2017;39:9.
 8.
Sutou S. Lowdose radiation from Abombs elongated lifespan and reduced cancer mortality relative to unirradiated individuals. Genes Environ. 2018;40:26.
 9.
Sutou S, Koeda A, Komatsu K, Shiragiku T, Seki H, Yamakage K, Niitsuma T, Kudo T, Wakata A. Collaborative study of thresholds for mutagens: proposal of a typical protocol for detection of hormetic responses in cytotoxicity tests. Genes Environ. 2018;40:20.
 10.
Rattan SIS, Bourg EL. Hormesis in health and disease. Boca Raton: CRC Press, Inc.; 2014.
 11.
Kiefer J. Biological radiation effect. Berlin Heidelberg: SpringerVerlag; 1990.
 12.
UNSCEAR. Sources and effects of ionizing radiation. UNSCEAR 2008 Report. New York: United Nations Publications; 2008.
 13.
Siegel JA, Sacks B, Socol Y. The LSS cohort of atomic bomb survivors and LNT. Comments on “Solid Cancer incidence among the life span study of atomic bomb survivors: 19582009” (Radiat Res 2017; 187:51337) and “Reply to the comments by Mortazavi and Doss” (Radiat Res 2017; 188:36971). Radiat Res. 2017;188:463–4.
Acknowledgements
Special thanks to Prof. H. Miyazawa (Tokushima Bunri University).
Funding
This work was supported by my school fund from Tokushima Bunri University. In addition, my idea in this manuscript was gotten on the basis of my other studies supported by research grants from Radiation Effects Association, from the Nakatomi Foundation, from KAKENHI and from the Japan Prize Foundation, and then I am grateful for the foundations.
Author information
Affiliations
Contributions
The author read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The author declares that he has no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Kino, K. The prospective mathematical idea satisfying both radiation hormesis under low radiation doses and linear nonthreshold theory under high radiation doses. Genes and Environ 42, 4 (2020). https://doi.org/10.1186/s4102102001454
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s4102102001454
Keywords
 Radiation hormesis
 LNT theory
 Inhibition effect