This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Mutagenesis and Background Neutron Radiation

Augusto González Instituto de Cibernética, Matemática y Física, La Habana
Abstract

We suggest a possible correlation between the ionization events caused by the background neutron radiation and the experimental data on mutations with damage in the DNA repair mechanism, coming from the Long Term Evolution Experiment in E. Coli populations.

pacs:
61.80.Hg, 87.53.-j, 87.23.Kg

I Introduction

In microelectronics, single failure events sporadically occur which, in some areas, like plane and space navigation, could have catastrophic consequences. Preliminary estimations IBM and more recent experiments exp indicate a correlation between these events and the Background Neutron Radiation (BNR) RFN . The mechanism of failure is the collision of a neutron from the BNR with an atomic nucleus in the chip, leading to a shower of electrons and ions that locally changes the conductivity and shortcuts the device.

In the present paper, we suggest the BNR as a cause of genetic “fails” in living cells, that is one of the possible origins of the so called spontaneous mutations. Cells exposed to the shower of electrons and ions, caused by the collision of a neutron and a proton of water, could be anihilated or experience a permanent damage, in particular, a damage in the DNA. The frequency of such events is similar to the rate of appearance of mutations with damage in the DNA repair mechanism Lenski , as measured in the Long Term Evolution Experiment (LTEE), where E. Coli populations evolve under controlled conditions EELP .

II The LTEE in E. Coli populations

The LTEE is an experiment conduced by Prof. R. Lenski and his group at the Michigan State University EELP . Each day, the bacteria undergo 6 - 7 generations of binary evolution. In a year, around 3400 generations occur. This means that, since the experiment started in 1988, it passed 60000 generations.

In the experiment, 12 populations of bacteria, with a common ancestor, independently evolve. Every day, 0.1 ml of the bacterial culture is serially transferred to 9.9 ml of a glucose solution, and mantained under controlled temperature until the next day. The number of bacteria varies approximately as shown in Fig. 1. That is, grows according to the law N02t/t0N_{0}~2^{t/t_{0}} in the first 8h, until the glucose is depleted, and then reach a stationary state. In the last 16h there is no appreciable mortality. The dependence 2t/t02^{t/t_{0}} is due to the way of reproduction, by cellular division.

Refer to caption
Figure 1: Daily evolution of the number of bacteria in a culture in the LTEE.

The experiment shows a set of very interesting results Lenski . We shall stress only two of them. First, in a given population, the total number of single point mutations in the DNA, after 20000 generations of evolution, is estimated as 3×1083\times 10^{8}. That is, the rate of point mutuations is:

fSPM1s1.f_{SPM}\approx 1~s^{-1}. (1)

On the other hand, in 2 of the 12 cultures, after 2500 - 3000 generations, mutations with a damaged DNA repair and edit mechanism appeared and became numerically dominant. A third line evolved the mutator phenotype after 8500 generations, and a fourth after 15000 generations. According to Prof. Lenski, the mechanism through which the mutator becomes numerically dominant is roughly the following. Once the mutation appears, the rate of spontaneous mutations increases 100 times, as compared with cells in which the DNA repair mechanism is not damaged. Thus, the mutator has a higher chance to generate the next winner and become dominant in a relatively short time scale, around 250 generations.

A second aspect, stressed by Prof. Lenski, is that mutations in which the DNA repair mechanism is damaged are “deleterious”, in the sense that a segment of the DNA is removed.

III Ionization events caused by the BNR

With regard to the BNR, we may assume that the cells live in pure water. Indeed, water is the main component of the solution, and the pH should be close to 7 in order to preserve life pH . In these conditions, the important processes are the collisions between neutrons, from the BNR, and the Hydrogen nuclei (protons) of water. The ejected proton gives rise to a shower of ions and electrons that is extended approximately 0.1 mm along the proton trajectory.

In Fig. 2, we show the flux per unit energy of neutrons in the BNR RFN , FF, in units of Neutrons/(MeV s cm2); the total cross section for neutron-proton dispersion sigma , σtotal\sigma_{total}, in units of 10-24 cm2; and the product FσtotalF~\sigma_{total}, in units of 10-24 s-1 MeV-1, as functions of the energy of incident neutrons. This last magnitude is proportional to the probability that a neutron with a given energy collides with a proton in water. It can be noticed that only neutrons with energy lower than a few MeV have a significant effect.

Refer to caption
Figure 2: Neutron flux per unit of energy in the BNR, FF, total cross section of the n+p collision, σtotal\sigma_{total}, and the product FσtotalF~\sigma_{total}, as functions of the energy of the incident neutron.

From these data, we may estimate the probability of neutron-proton collisions:

Probn+p\displaystyle Prob_{n+p} =\displaystyle= Np01000MeVdEFσtotal\displaystyle N_{p}\int_{0}^{1000~MeV}{\rm d}E~F~\sigma_{total} (2)
=\displaystyle= 8×103s1,\displaystyle 8\times 10^{-3}s^{-1},

where Np=6.6×1023N_{p}=6.6\times 10^{23} is the number of protons in 10 ml of water.

In order to compute the energy transfer from the proton to water, we use the data in Fig. 2 and the so-called stopping power of protons in water, tabulated in Ref. sp, . By using a Monte Carlo algorithm, we arrive to the results shown in Fig. 3. According to Ref. sp, , energy losses are mainly due to the interaction of the proton with the electrons in the medium, leading to the ionization of water molecules. The basic process of ionization is: H2Oe+H2O+{\rm H_{2}O\to e+H_{2}O^{+}}, which requires an energy of 12.62 eV Ip . The ejected electron and H2O+{\rm H_{2}O^{+}} could lead to secondary ionization processes. Dividing the yy axis of Fig. 3 by 10, we obtain a rough estimate of the number of ions produced in each 100 nm step of the proton motion, that is around 300 ions at distances close to the n-p collision point, and 30 ions when distance is of the order of 0.1 mm.

Refer to caption
Figure 3: Energy transfer from the proton to the medium, in 100 nm steps, along the proton trajectory.

IV BNR effects on the LTEE

As mentioned above, n+p dispersion events in the glucose solution take place every 125 s. We already know that the shower of ions and electrons, created by the proton, is more intense in the first 0.1 mm along the proton trajectory. The bacteria touched by this ion shower could be destroyed or experience a permanent damage, especially in their DNA, which can be later inherited by the descendants. We shall stress that, for the DNA changes to be transmitted, the ionization event should take place in the first 8h of daily evolution, according to Fig. 1. Otherwise, there is practically no cellular division in the day it occurred, and the probability to pass to the next day is only 1/100.

The mean number of bacteria in the first 8h is:

N¯=N0808dt2t/t0=21.5N0,\bar{N}=\frac{N_{0}}{8}\int_{0}^{8}{\rm d}t~2^{t/t_{0}}=21.5~N_{0}, (3)

where N0=5×106N_{0}=5\times 10^{6} bacteria, and 28h/t0=1002^{8h/t_{0}}=100. Each bacterium occupies a mean volume of around 10 cm/3(21.5N0){}^{3}/(21.5~N_{0}), that is, a cube with sides 45 μ\mum long. In the first 0.1 mm=100 μ\mum of the ion shower, only 2 such cubes could be allocated. The probability that the shower touches a bacterium is, thus:

2ShowerVolumeCubeVolume=2l2×45μm(45μm)3,2~\frac{\rm Shower~Volume}{\rm Cube~Volume}=2~\frac{l^{2}\times 45~\mu m}{(45~\mu m)^{3}}, (4)

where ll is the lateral dimension of the ion shower. ll could be estimated from the Debye screening length of pure water:

λD=(kBTϵϵ02nq2)1/2,\lambda_{D}=\left(\frac{k_{B}T\epsilon\epsilon_{0}}{2nq^{2}}\right)^{1/2}, (5)

where kbk_{b} denotes the Boltzman constant, ϵ80\epsilon\approx 80 is the relative dielectric constant of water epsilon , q=1q=1 is the charge of the ions H+ y OH- in water, and nn their concentration:

n\displaystyle n =\displaystyle= 107×(3×1022molecules/cm3)\displaystyle 10^{-7}\times(3\times 10^{22}~{\rm molecules/cm^{3}}) (6)
=\displaystyle= 3×1015ions/cm3.\displaystyle 3\times 10^{15}~{\rm ions/cm^{3}}.

Taking all these numbers together, we get λD=500nm=0.5μm\lambda_{D}=500~nm=0.5~\mu m. And putting l=λDl=\lambda_{D} in Eq. (4), we get a probability of 2.5×1042.5\times 10^{-4}. Notice that ll is a magnitude of the same order of the E. Coli dimensions, thus the ion shower may cause strong effects on a bacterium.

We may compute the rate in which bacteria from a single population are touched by the BNR ionization events:

fBNR\displaystyle f_{BNR} \displaystyle\approx (2.5×104)×(8×103s1)\displaystyle(2.5\times 10^{-4})\times(8\times 10^{-3}~s^{-1}) (7)
\displaystyle\approx 2×106s1.\displaystyle 2\times 10^{-6}s^{-1}.

This number is very small, as compared with fSPMf_{SPM}, Eq. (1). However, it is consistent with the frequency of deleterious mutations, with damage in the DNA repair mechanism, mentioned in section II. Indeed, in Δt2400\Delta t\sim 2400 generations 1\sim 1 year 3×107\sim 3\times 10^{7} s, the BNR had a direct incidence on fBNR×Δt60f_{BNR}\times\Delta t\sim 60 bacteria. Some of them could have experienced damages in the DNA repair mechanism. The 100 times increase in the mutation rate could have given this subpopulation, after 100 - 600 generations, the possibility to generate beneficial mutations that would be fixed, allowing them to become numerically dominant. The fact that only 4 of 12 populations evolved in this way could be related to the probability 1/3\sim 1/3 that the BNR events take place in the first 8h of daily evolution.

Let us notice that we are assuming very fast BNR ionization events, as compared with the bacterial motion. Only those bacteria placed along the ion shower are affected by it. We may estimate the duration of such a event from:

τ1ϵϵ0σ,\tau_{1}\approx\frac{\epsilon\epsilon_{0}}{\sigma}, (8)

where σ=5.5×106Coul/(Vsm)\sigma=5.5\times 10^{-6}~{\rm Coul/(V~s~m)} is the conductivity of pure water conduct . That is, τ1104\tau_{1}\approx 10^{-4} s.

A second estimate for the duration comes from the difussion constants of ions in water D , D103μm2/sD\approx 10^{3}~\mu m^{2}/s. Taking λD=0.5μm\lambda_{D}=0.5~\mu m as a characteristic dimension, results in:

τ2λD2D2.5×104s.\tau_{2}\approx\frac{\lambda_{D}^{2}}{D}\approx 2.5\times 10^{-4}s. (9)

In both cases, the times are of the order of 10410^{-4} s. Taking into account that, at ambient temperatures, the typical speeds of bacterial motion are around 2 mm/s, only bacteria in contact with the ion shower, or very close to it, will be affected.

The fact that mutations with damage in the DNA repair mechanism are deleterious Lenski is also consistent with the nature of BNR ionization processes. Indeed, the electron and ion shower is highly energetic and may produce such damages in the DNA, especially in the first steps after the n+p collision.

We shall compare the concentration of produced ions with the concentration of spontaneous ions in water, Eq. (6). In each of the first 100 nm steps, the ejected proton creates around 300 ions. The induced concentration is, thus:

nind=3000.52×0.1×μm3=1.2×1016ions/cm3,n_{ind}=\frac{300}{0.5^{2}\times 0.1\times{\rm\mu m^{3}}}=1.2\times 10^{16}{\rm ions/cm^{3}}, (10)

that is, 4 times higher than nn given in Ec. (6). The presence of ions in such high concentrations is also a strong mutagenic factor.

V Concluding Remarks

In the present paper, we indicate a possible correlation between BNR ionization events and the LTEE observed rates of deleterious mutations with damages in the DNA repair and edit mechanism. In this way, we are indicating the probable origin of a class of “spontaneous” mutations.

The experimental confirmation of this possible correlation is plausible: restart the experiment by using fossils, and shield some of the evolving populations against the BNR. The shielded cultures should exhibit much lower rates for deleterious mutations with damages in the DNA repair mechanism. In around 1 - 2 years (2500 - 5000 generations), changes in mutation rates should be manifest.

On the other hand, a comment by Prof. Lenski Lenski that some cancer cells also exhibit damages in the DNA repair mechanism, motivates us to rise the hypothesis about the BNR as one os the processes triggering cancer. Other events, like inhalation of radioactive Radon contained in air through breathing, are recognized carcinogens Radon . Defficient feeding, infectious proceeses, etc could be considered as conditions creating an evolutive pressure over the expossed cells, similar to the limited amount of glucose in the LTEE. Under these conditions, the BNR induced deleterious mutations, with damages in the DNA repair and edit mechanism, and the subsequent rise in the rate of spontaneous mutations, could allow the mutators to generate well adapted individuals that could become numerically dominant. In order to check this hypothesis, a controlled experiment in animals could be designed, for example in mouses, which are widely used as models of cancer in humans mice .


Acknowledgements.
The author acknowledge C. Ceballos for the information on the role of BNR in microelectronics, and E. Altshuler, A. Cabo, C. Cruz, G. Martín and E. Moreno for their comments and criticism.

References

  • (1) J.F. Ziegler, H.W. Curtis, H.P. Muhlfeld et. al., IBM experiments in soft fails in computer electronics (1978-1994), IBM Journal of Research and Development, Vol. 40, No. 1, page 3, 1996.
  • (2) T. Nakamura, M. Baba, E. Ibe, Y. Yahagi, H. Kameyama, Terrestrial neutron-induced soft errors in advanced memory devices, Singapore, World Scientific, 2008; M. Nicolaidis, Ed., Soft errors in modern electronic systems, Springer, New York, 2010.
  • (3) M.S. Gordon, P. Goldhagen, K.P. Rodbell, et. al., IEEE Transactions on Nuclear Science 51 (6), 3427 (2004).
  • (4) R.E. Lenski, Phenotypic and genomic evolution during a 20000 generation experiment with the bacterium E. Coli, in J. Janick, Ed., Plant Breeding Reviews, Vol. 24, Part 2, page 225, 2004.
  • (5) R. Lenski, Summary data from the long-term evolution experiment, http://myxo.css.msu.edu/ecoli/summdata.html
  • (6) W. Boron, E.L. Boulpaep, Eds., Medical Physiology: A Cellular And Molecular Approach, Elsevier, 2009.
  • (7) J.W. Norbury, Nucleon-Nucleon Total Cross Section, NASA/TP-2008-215116.
  • (8) PSTAR : Stopping Power and Range Tables for Protons, http://physics.nist.gov/PhysRefData/Star/Text/PSTAR.html
  • (9) NIST data on the ionization potential of water, http://webbook.nist.gov/cgi/cbook.cgi?ID=C7732185&Mask=20
  • (10) M. Uematsu and E.U. Franck, J. Phys. Chem. Ref. Data, Vol. 9, No. 4, page 1291, 1980.
  • (11) R.H. Shreiner and K.W. Pratt, Primary Standards and Standard Reference Materials for Electrolytic Conductivity, NIST Special Publication 260-142, 2004 Ed.
  • (12) E.L. Cussler, Diffusion: Mass Transfer in Fluid Systems, New York, Cambridge University Press, 1997.
  • (13) National Research Council. Committee on Health Risks of Exposure to Radon: BEIR VI. Health Effects of Exposure to Radon. Washington, National Academy Press, 1999.
  • (14) http://www.nih.gov/science/models/mouse/resources/hcc.html