Abstract
Microdosimetry is a tool for the investigation of microscopic energy deposition of ionizing radiation. This work used
Introduction
The biological effect of ionizing radiation generally depends on the actual energy deposition of rays in organisms. 1 -4 Due to the randomness of interaction between rays and substances, the macroscopic dosimetry such as absorbed dose and dose rate can only reflect the average level of the energy deposition of the irradiated substance. When the scale is reduced to micron dimension or even smaller, a remarkable deviation occurs between the practical radiation energy deposition and the absorbed dose. 5 -7 In order to study the microscopic distribution of radiation energy, the basic concept of micro dosimetry was put forward. 5 Dose deposition distribution and fluctuation of ionizing radiation in microscopic levels such as cell and subcell were surveyed. Microdosimetry has been widely utilized in the study of radiation biological effect, solid radiation effect, radiation protection, radiation therapy and so on. 8 -10
The combination of microdosimetry and biology is of significance. Biological effect of radiation is induced by radiation energy deposition. For the systematic comprehension of the relationship between radiation dose and biological effect, microdosimetry and biological effect should be combined organically. Cell or animal models used previously in biological study were not practicable in physical modeling owing to excessive cells. Moreover, there are various kinds of cells in common animal models, which are lack of uniform evaluation standard. Therefore, an animal model with simple structure is required to establish a combined model of microdosimetry and biology.
As one of the classical model animals,
In the study of biological effect of ionizing radiation, the absorbed dose of
Materials and Methods
Culture and Maintenance of Nematodes
The N2 strain was obtained from the group of Huimin Zhang, College of Basic Medicine and Life Sciences, Soochow University. It was used as the wild-type
The petri dish with sufficient adult worms was selected and embryos were washed off the cell-strainer with M9 buffer. After being digested with bleach solution and washed with M9 buffer for thrice, embryos were incubated on non-seeded NGM plates without OP50 for 24 h at 20°C. The following day, synchronous L1 nematodes were transferred onto a prepared agar pad on a glass slide.
Irradiation
The external gamma radiation exposure was performed at the State Key Laboratory of Radiation Medicine and Protection, School of Radiation Medicine and Protection, Soochow University. The entire body of L1 nematodes was exposed to 25 Gy/min of gamma radiation (total doses ∼100 Gy). The worms were transferred to a fresh NGM agar medium containing OP50 at 20°C immediately after irradiation.
Effects on Somatic Growth and Reproduction
L1 nematodes were photographed with a microscope (Normaski). For the establish of nematode model, geometric parameters of the whole body and gonad of 30 nematodes were measured with imageJ software. The length (size) of nematodes was determined at 48 h of exposure onward from L1 stage.
The young adult nematodes were then placed on NGM plates and immersed with K-medium solution. A head thrash was defined when the length of body was bent to half. Times of head thrash per minute was randomly recorded with 15 nematodes in each group.
After exposure, the young adult nematodes were isolated and placed in NGM medium with one in each dish, subsequently transferred to fresh petri dish every 24 hours. The egg-laying amount for 3 consecutive days was recorded with 30 nematodes in each group using stereomicrograph. Nematodes died during the experiment were excluded in this study.
Model Establishment and Monte Carlo Simulation
In this work, a model was constructed and Monte Carlo simulation was performed on the irradiated nematodes to obtain the absorbed dose distribution of nematodes and gonads. Irradiation source, culture dish, culture medium, body and gonad of nematodes were taken into account in the model. Gamma photons with energy of 1.17 and 1.33 MeV emitted by 60Co decay were the irradiation source. The petri dish made up from PMMA was used with a diameter of 6 cm. For simplification, the medium was simulated as liquid water. According to the measured result of the size of whole body and gonads, cuboid with 384.0 µm in length, 25.8 µm in both height and width were simulated for nematode size, while cuboid with 20.3 µm in length, 10.0 µm in both height and width were simulated for gonadal size.
In this study, 200 nematodes were simulated to distribute randomly as a first approximation. That is, nematodes were cultured in a random distribution in the medium and they were assumed to have no intersections. The recent version (version 3.20) of the Particle and Heavy-Ion Transport code System (PHITS) was used in Monte Carlo simulation, 19,20 which has been widely used for various applications, such as low-energy neutron interaction, 21 beam transport functions, 22 and a microdosimetric tally function. 23 The cut-off energy of photon and electron were both set at 1 keV. To improve efficiency, variance of forced collision in the nematodes was declined. The simulated photon number were 20E8 with the maximum relative error less than 4.8%, and the error was calculated by more than 5 times. Simulations were executed on a system with 44 Intel(R) Xeon(R) E5-2696 and 256 GB RAMs.
Statistical Analysis
SPSS 20.0 software was used for statistical analysis. Non-paired
Results
Establishment of Microdosimetry Model and Dose Estimation for C . Elegans
To reveal the potential reason for the changes of potential biological effects after irradiation, nematodes model upon 60Co irradiation exposure was established by Monte Carlo software PHITS (Figure 1.a, b, c). Nematodes were presumed to be in water layer. When the average absorbed dose of the water layer (nematodes were not included) reached 100 Gy, the simulated absorbed doses for the whole body and gonad of 200 nematodes were listed as follows: the absorbed dose of body was 99.9 ± 57.8 Gy (19.6-332.2 Gy), that of gonads were 129.4 ± 558.8 Gy (9.5-6597 Gy). To reflect the general level, the extremely high doses of gonads were excluded, finally the gonad absorbed dose was determined to be 19.8 ± 5.1 Gy. Meanwhile, the frequency distribution of the estimated dose suggested that the dose of 50% nematodes body ranged from 50 to 110 Gy, and that of 11 cases (5.5%) exceeded the extreme value 210 Gy (Figure 1.d).After excluding extremly high doses, the dose of 50% gonads distributed between 16.0 and 22.7 Gy (Figure 1.e).

Microdosimetric model Establishment and dose estimation of C. elegans. a. Schematic diagram of body (red) and gonad (black) of C. elegans (cuboids). b. Schematic diagram of the relative position of body and gonad. The gonad was located in the ventral body. c. Schematic diagram of irradiation conditions. Nematodes were distributed relatively uniformly on the medium (blue). d. Frequency distribution diagram for dose of nematodes body. e. Frequency distribution diagram for nematodes gonads (with outliers excluded).
Fitting of Growth and Changes of Movement Behavior With the Dose of Nematodes
The effects of gamma radiation on

Effects of ionizing radiation on growth and movement of C. elegans. a. Size of nematode in the control group. b. Size of nematode in irradiation group. c. Effects of ionizing radiation on the length of C. elegans body. d. Effects of ionizing radiation on the head thrash frequency of C. elegans. e. Frequency distribution for dose of irradiated nematodes. f. Frequency distribution of head thrash frequency in the irradiation group. g. Frequency distribution profile of worm body dose and head thrash frequency. Asterisks indicate significant difference from control treatment (p-value < 0.05). Scale bar: 200 μm.
The relationship between the changes of head thrash after irradiation and the estimated microdose of C. elegans was investigated. Nematodes receiving the radiation dose higher than 99.9 Gy was assigned to one group based on the microdose frequency distribution. According to the frequency distribution of the head swing, nematodes with the frequency less than 84 times/min were allocated to one group. Finally, the contour analysis of the 2 processed frequency distribution curves demonstrated that there was no statistical difference between them (p > 0.05) (Figure 2.g), implying that there was connection between the frequency distribution of head swing and the estimated worm dose.
Fitting of the Egg Laying Amount Variation of Nematode With the Gonad dose
Injury of germ cells induced by ionizing radiation directly resulted in the decrease of the egg laying amount (Figure 3a, b). The total egg laying amounts of

Effects of ionizing radiation on C. elegans reproductive system and analysis. a. Fecundity of nematodes in the control group. b. Fecundity of nematodes in the irradiation group. c. Effects of ionizing radiation on egg laying amount. d. Frequency distribution of egg laying amount in irradiation group. e. Frequency distribution for dose of gonads after irradiation. f. Frequency distribution profile of gonadal dose and egg laying amount. Asterisks indicate significant difference from control treatment (p-value < 0.05).
The connection between alterations of spawning amounts after irradiation and gonad estimated result by microdosimetry was studied. According to the frequency distribution of dose, gonads exposed to high doses (≥19.8 Gy) of radiation were assigned to one group (Figure 3.d). Based on the frequency distribution of spawning amount from the results of reproductive capacity, nematodes with the spawning amount less than 14 were classified to one group (Figure 3.e). The contour analysis of the 2 processed frequency distribution curves showed there was no statistical difference between these 2 groups (
Discussion
In this study, a microdosimetry model of Caenorhabditis elegans irradiated by 60Co irradiation device was established by using Monte Carlo software PHITS for the first time. The difference of dose micro-distribution was verified by the determination of head thrash and brood size. Head thrash is an indicator of swimming behavior of nematodes, which is mainly regulated by cholinergic neurons. Among the cholinergic related genes, ric-3 is located in pharynx and body wall, 24 while lev-10 and cam-1 are located in muscle cells. 25,26 From the expression pattern of each gene, it can be seen that cholinergic neurons are relatively evenly distributed in the whole body. It has been reported that ionizing radiation inhibits movement through systemic mechanisms, which may involve motor neurons and / or body wall muscle cells. 27 Therefore, head thrash was used as an indicator to measure the dose of the worm body. The brood size is a recognized biological end point, 28 which is often used in toxicological studies and can reflect the damage of the reproductive system to a certain extent. 29,30 The brood size is mainly related to the development of germ cells of nematodes. It was reported that ionizing radiation can induce apoptosis of damaged germ cells through DNA damage checkpoints, resulting in reduction of the brood size. 16 Therefore, the brood size was selected as an indicator to reflect the gonadal dose in this study. In the study of radiobiology, it was noted that the bystander effect should not be ignored. The bystander effect refers to the phenomenon that the biological effects of ionizing radiation can be observed in unexposed neighboring cells to an irradiated cell or group of cells, 31 -33 which had been confirmed in C.elegans. Although the irradiation condition of this study was performed on whole body, which is different from the local irradiation used to study the bystander effect, the side effect may still affect the experimental results. Moreover, with the deepening of the later study, the irradiation conditions may change from whole body irradiation to local irradiation, so it is necessary to consider the bystander effect and distinguish the radiation biological effect of the irradiated tissues from the bystander effect.
In the establishment of microdosimetric model, the overall dose was concerned, so the overall geometric shape of worm body and gonad was mainly considered. However, the internal structure of nematode and gonad was not distinguished, which are regarded as homogeneous tissues. Due to the low simulated volume and the difficulty of tangent calculation, particles in surface transport may easily “leak” particles (“particle lost”). Thus, this study adopted a rectangle for simulation instead of slender ellipsoid and the cylinder.
The average dose of the worm was estimated to be 99.9 Gy by the microdosemetry model, which was consistent with the set radiation dose 100 Gy. However, there were remarkable differences between individual worms in the dose as for the microscopic dose distribution. The distinguish was more visible for gonads (may be one thousand times), implying that the difference of dose distribution may increase with the decrease of size of the objective. This result was consistent with the conclusion of Kellerer
At present, cell model was adopted in the study of microdosimetry. Nevertheless, except for weak maneuverability, cells can only be cultured in
The possibility of
Conclusion
C.elegans is known to be sensitivity of to ionizing gamma radiation. On the microscopic scale, the energy deposition of ionizing radiation is significant to study the biological effects of ionizing radiation.
Footnotes
Authors’ Note
Na Chen and Liang Sun conceived and designed the experiment. Wentao Yu and Huiqiang Long performed the biological experiments. Yidi Wang and Liang Sun performed the Microdosimetric experiments. Jin Gao analyzed the data. Wentao Yu wrote the article. Wentao Yu and Huiqiang Long and Jin Gao contributed equally to this work.
Acknowledgments
We are grateful to the group of Huimin Zhang, College of Basic Medicine and Life Sciences, Soochow University for providing the C. elegans strain.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (grant numbers 82003391, U186720), the Natural Science Research Projects of Colleges and Universities in Jiangsu Province (grant numbers 20KJB310007), the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD), the Collaborative Innovation Center of Radiological Medicine of Jiangsu Higher Education Institutions, China and the Nuclear Energy Development Project, China (No. 2016-1295).
