Abstract
Rats at various ages were observed to present with different radiosensitivity and bioavailability for radiotracers commonly used in preclinical research. We evaluated the effect of age-induced changes in body weight on radiation dose calculations. A series of rat models at different age periods were constructed based on the realistic four-dimensional digital rat whole-body (ROBY) computational model. Particle transport was simulated using the MCNPX Monte Carlo code. Absorbed fractions (AFs) and specific absorbed fraction (SAFs) of monoenergetic photons/electrons and S values of eight positron-emitting radionuclides were calculated. The SAFs and S values for most source-target pairs were inversely correlated with body weight. Differences between F-18 S values for most source-target pairs were between −1.5% and −2%/10 g difference in bodyweight for different computational models. For specific radiotracers, the radiation dose to organs presents a negative correlation with rat body weight. The SAFs for monoenergetic photons/electrons and S values for common positron-emitting radionuclides can be exploited in the assessment of radiation dose delivered to rats at different ages and weights. The absorbed dose to organs is significantly higher in the low-weight young rat model than in the adult model, which would result in steep secondary effects and might be a noteworthy issue in laboratory animal internal dosimetry.
THE RAT AS A MODEL of human disease is prized as a preeminent laboratory animal in biomedical research owing to the various advantages it offers in translational animal research. This includes the fact that the rodent genome is similar to the human genome and a large number of rodents can be reared in a short time; moreover, the breeding of rats in the laboratory is easy and less costly. At the beginning of the twentieth century, The Wistar Institute developed the first notable rat strain as a model organism for biological and medical research. 1 The Wistar rat is one of the most popular rat strains used for laboratory research, and more than half of currently used laboratory rat strains are descendents from the original colony established by this institute. 2 They are widely used in preclinical trials for tracer development, test of novel medical imaging methods, and many other biomedical research areas.3–5 Clear pictures of disease in small animals can be obtained using dedicated high-resolution preclinical positron emission tomography (PET) instrumentation, which enables researchers to track molecular events in vivo and probe the biological basis of human disease. However, in typical preclinical PET–computed tomography (CT) imaging studies, the animal is often exposed to nonnegligible radiation exposure, which may preclude the use of the animal in longitudinal studies, where it is used as its own control. In this context, the accurate assessment of radiation dose delivered to laboratory animals in experiments involving the use of stand-alone or hybrid imaging techniques (eg, PET, PET-CT) is of paramount importance.
For this purpose, various computational models have been developed to mimic the morphology and internal anatomic structures of small animals and used to assess the radiation dose delivered to laboratory animals in preclinical research studies. Hui and colleagues developed the first computational mouse model using simple equation-based mathematical functions to represent internal organs. 6 Further improvement led to finer mathematical small-animal models containing more internal organs and even realistic tumor models.7–11 With the popularity and widespread availability of Monte Carlo–based radiation transport packages and the growth of high-performance computing, including cloud computing, voxel-based small-animal models of different specimens employing voxel matrices derived from tomographic images, which best preserve anatomic authenticity, have been developed for imaging and radiation dosimetry applications.12–21 However, there is an intrinsic variability in organ mass and total body weight between these computational models owing to individual differences between studies specimens. Moreover, the body weight and length of laboratory animals used in longitudinal studies vary according to age periods. The rat at different ages can be linked to human beings in various growth phases22–24 and was observed to have different radiosensitivity and bioavailability for radiotracers.25–30 However, the intrinsic disadvantages of traditional mathematical and voxel-based models, such as the loss of most anatomic details and the difficulty in reshaping the contours of organs, obstructed the progress of radiation dosimetry research for small animals at different age periods. More recently, nonuniform rational B-spline surfaces (NURBS)-based modeling techniques have become popular for Monte Carlo–based radiation transport calculations.21,31,32 The NURBS-based model holds most anatomic characteristics of internal organs and tissues and offers the flexibility of deforming these structures more realistically. As such, it became a precious tool for the assessment of age-dependent small-animal radiation dosimetry.
In this work, we used the four-dimensional digital rat whole-body (ROBY) NURBS-based Wistar rat model to produce a series of models representing the laboratory rat at 10 different age periods ranging between 4 weeks (prepubescent period) and 45 weeks (adult phase). For the adult rat older than 300 days, weight growth was reported to be due mostly to fat deposition. 33 Monte Carlo–based particle transport simulations of monoenergy photons/electrons and spectra corresponding to decay schemes of positron-emitting radionuclides (C-11, N-13, O-15, F-18, Cu-64, Ga-68, Y-86, and I-124) were performed using MCNPX (Los Alamos National Laboratory, Los Alamos, NM) to calculate the absorbed fractions (AFs), specific absorbed fractions (SAFs), and S values for the considered series of models. The obtained radiation dosimetry database can be used for the assessment of radiation dose to small animals in experiments involving the use of specimens at different age periods. Understanding the variability of the absorbed dose with body weight is important to appreciate the uncertainties involved in small-animal internal radiation dosimetry calculations.
Materials and Methods
Design of Computational Rat Models
Given that the Wistar rat is one of the most popular rat strains used in small-animal research, we constructed a series of anatomic models using the realistic ROBY model 31 derived from magnetic resonance images of a Wistar specimen according to published weight and length values for this strain. The growth chart was obtained from 2-year study data on the Wistar Hannover GALAS rat 34 and used to calculate the length of rats at different age periods based on the weight-length formulas (no gender distinction) given by Donaldson 35 :
The body weight and length values summarized in Table 1 were used as the basis for modification of the original ROBY model. The models representing 4-, 6-, 8-, 10-, 12-, 15-, 21-, 27-, 35-, and 45-week-old Wistar rats were constructed and referred to as rat models 1 to 10, respectively. The resulting organ masses were also adjusted to follow reference data reported in the literature. 35 The resulting organ and body masses were calculated by multiplying each organ volume by its corresponding density. Traditional hollow organs (eg, heart, stomach, intestines, and bladder) were treated as uniform organs, with mass corresponding to the sum of wall and content. The skeleton was treated as a mixture of 70% bone, 25% red bone marrow, and 5% yellow bone marrow. 36 Table 2 summarizes the body size (lateral size, anteroposterior size, and length), model weight, and calculated mass for each organ of the constructed anatomic models. Figure 1 shows ventral-dorsal views of representative rat models.
Growth Chart and Length for the Wistar Rat at Different Ages
The growth chart was obtained from 2-year study data on the Wistar Hannover GALAS rat. 34
The length was calculated according to the weight-length formula given by Donaldson. 35
Summary of Organ Masses and Body Sizes of the 10 Computational Rat Models at Different Ages

Three-dimensional visualization of representative computational rat models showing (A) the 4-, (B) 8-, (C) 12-, (D) 21-, and (E) 45-week-old models.
Monte Carlo Simulations
The modified NURBS-based rat models were voxelized using the ROBY software 31 and used as input to MCNPX Monte Carlo code version 2.5.c 37 to simulate the transport and interaction of photons, electrons, and positrons in the computational rat models. The number of voxels for each identified region was calculated and multiplied by the voxel volume (0.5 × 0.5 × 0.5 mm3) and corresponding tissue density to yield the region mass. The density and chemical composition of each organ were assumed to be similar to those of human beings.38,39 Uniformly distributed monoenergetic photon and electron sources, with 20 discrete photon energies ranging from 0.01 to 3.0 MeV and 13 selected electron energies ranging from 0.1 to 3.0 MeV, and eight positron-emitting sources (C-11, N-13, O-15, F-18, Cu-64, Ga-68, Y-86, and I-124) were simulated in 14 chosen source organs. A total of 6 × 106 primary particle histories were generated, resulting in a statistical uncertainty, in terms of coefficient of variation, of less than 2% in most cases. Ten parallel MCNPX simulation processes were run on a 12-core workstation with 2.4 GHz Intel Xeon processors and 32 GB RAM running Microsoft Windows 7 OS. The average calculation time corresponding to the number of particles reported for each simulation process is 156 minutes. The decay data of the positron-emitting radionuclides were taken from the Health Physics Society electronic resource. 40
Absorbed Dose Calculation
The equation for mean absorbed dose calculation according to the Medical Internal Radiation Dose (MIRD) formalism is given by41,42
where D(rT,TD) = the mean absorbed dose to a target region rT over the dose integration period; TD; Ã (rS,t) = the time-dependent activity in source tissue rS; and S(rT←rS) = the S value describing the equivalent dose rate to the target organ rT per unit activity in the source organ rS.
For certain radionuclides, the S value is given by
where Ei is the individual energy of the ith nuclear transition, Yi is number of ith nuclear transitions per nuclear transformation, M(rT) is the mass of the target tissue rT, and ϕ (rT←rS,Ei) is the AF of radiation energy Ei in the target organ for the ith radiation type originating in the source organ, which describes the proportion of energy released in source organs deposited in the target organ. The ratio of AF and target mass is known as the SAF, which describes the mean AF in the target organ.
If Ã(rS,t) is normalized to a unit administered activity A0, we obtain the mean absorbed dose per administered activity (absorbed dose coefficient):
The fractional activity Ã(rS,t)/A0 at a given time point of activity measurement is calculated by dividing the multiplicative combination of measured uptake (%ID/g) and the tissue weight by 100. The time–fractional activity curve (TAC) for considered source organs rs is derived from graphs of the fractional activity Ã(rS,t)/A0 versus time (in seconds). The data are fitted to two-component exponential curves 43 using the nonlinear regression method. The activity as a function of time is given by:
where C1 and C2 are coefficients and k1 and k2 are time constants.
Results
The AFs and SAFs for monoenergetic photons and electrons for 16 source organs irradiating 17 target organs were calculated. Supplemental Tables S1 to S15 (available online only) provide detailed photon and electron SAF results and S values of the eight positron-emitting radionuclides for the considered source-target pairs in selected rat models. Figure 2 shows photon self-absorbed fractions for the liver and the heart for the 10 generated rat models. Comparisons with the results of Xie and colleagues, Stabin and colleagues, and Peixoto and colleagues are also made.8,12,18 A consistent trend can be observed for the different rat models. The self-absorbed fractions present a negative correlation with photon energy and a positive correlation with the rat model's weight because increasing photon energy facilitates the escape of recoil electrons and scattered photons from the source organ, whereas the increased weight is accompanied by a larger volume, which extends the distance of particle transport and causes more energy deposition in the source organ. Figure 3 shows photon and electron self-absorbed SAFs for the stomach and kidney. The self-absorbed SAFs decrease when the particle energy or the body weight increases; the latter outcome can be deduced from the definition of the SAF itself. In human internal radiation dosimetry, the electron cross-absorbed SAF for organs is commonly neglected 44 because the range of low-energy electrons is smaller than the dimension of most human organs. 45 However, given that the body dimension of rats is much smaller, it is necessary to evaluate the cross-absorbed SAF in small-animal dosimetry calculation. Figure 4 illustrates photon and electron cross-absorbed SAFs for the intestine irradiating the stomach and kidney. The increasing curves of electron cross-absorbed SAFs in Figure 4, B and D, indicate that more electrons escape from the source organ and stop in surrounding tissues when energy increases. The difference between the shapes of self-absorbed (see Figure 3, A and C) and cross-absorbed (see Figure 4, A and C) photon SAF curves is more pronounced in the lower-energy portion. The rising portion of cross-absorbed photon SAF curves is linked to the increase in the energy deposition in target organs from particles escaping from source organs as the initial energy increases. Above a certain initial photon energy, more and more particles escape from the rat body, consequently causing the drop of photon cross-absorbed SAFs. Figure 5 shows the SAFs for the total body irradiating the skeleton, which expresses high radiosensitivity and is considered one of the most important dose-limiting tissues. As can be seen in Figure 5A, the skeletal cross-absorbed photon SAFs are 170 to 530% higher in the 4-week-old model and 24 to 47% higher in the 15-week-old model compared to the 45-week-old model. On the other hand, the skeletal cross-absorbed electron SAFs are about 450% higher in the 4-week-old model and 43% higher in the 15-week-old model than in the 45-week-old model. Figure 6 shows photon and electron SAF ratios between various rat models (from the 4-week-old model to the 35-week-old model) to the rat model 10 (45 weeks old) for the total body irradiating itself. It was observed that when the emitted energy increases, the absolute difference in photon self-absorbed SAFs of total body between various models starts initially with a rapid decrease followed by a steady portion in the energy range 0.05 to 1.2 MeV and finally drops gradually. For electron SAFs of the total body irradiating itself, the absolute difference between various models slightly decreases while increasing electron energy.


Photon self-absorbed specific absorbed fractions (SAFs) for (A) the stomach and (C) the kidney and electron self-absorbed SAFs for (B) the stomach and (D) the kidney.

Photon cross-absorbed specific absorbed fractions (SAFs) for the intestine irradiating (A) the stomach and (C) the kidney and electron cross-absorbed SAFs for the intestine irradiating (B) the stomach and (D) the kidney.

(A) Photon and (B) electron specific absorbed fractions (SAFs) for the total body irradiating the skeleton.

Ratios of self-absorbed specific absorbed fractions (SAFs) for the total body between the nine computational rat models and the oldest rat model (46 weeks) for (A) photons and (B) electrons.
Figure 7 and Figure 8 show self-absorbed S values for the kidney, skeleton, heart, and lung and cross-absorbed S values of different source-target organ pairs for the constructed rat models and the eight considered positron-emitting radionuclides (C-11, N-13, O-15, F-18, Cu-64, Ga-68, Y-86, and I-124). For all positron-emitting radionuclides, the S values decrease, whereas the model weight increases. As can be seen in Figure 7, among the considered radionuclides, O-15 emits positrons with high energy and frequency and produces the highest self-absorbed S values in most organs. Interpreting the dosimetric results of cross-absorbed S values is more complicated than interpreting those of self-absorbed S values. Y-86 generates the highest cross-absorbed S values in most source-target organ pairs for different models, whereas Ga-68 induces the highest cross-absorbed S values in adjacent target organs for models with a lower weight. The relative differences in S values per 10 g difference in total body mass (%/10 g) for F-18 between the 4-week-old and 45-week-old models are summarized in Table 3. The relative difference ranges between − 1.5% and −2%/10 g for most source-target pairs. Figure 9 shows the contributions of positrons and the two annihilation photons to the total body self-absorbed S values of C-11, N-13, O-15, and F-18 in the 10 rat models. It can be seen that for these four radionuclides, the larger contribution (more than 70%) to the self-absorbed S values is due to positron interactions, whereas the contribution of the two annihilation photons increases gradually with the weight.

Representative self-absorbed S values for the (A) kidney, (B) skeleton, (C) heart, and (D) lung as a function of body mass for different positron-emitting radionuclides.

Representative cross-absorbed S values for (A) the liver irradiating the stomach, (B) the kidney irradiating the intestine, (C) the liver irradiating the spleen, and (D) the kidney irradiating the heart as a function of body mass for different positron-emitting radionuclides.

Proportions of contributions to self-absorbed S values of total body from positrons and the two annihilation photons for (A) C-11, (B) N-13, (C) O-15, and (D) F-18.
Percent Difference between F-18 S Values per 10 g difference (%/10g) in Total Body Mass from the 4-Week-Old to the 45-Week-Old Rat Model
We selected from the literature a 11C-labeled probe used in small-animal PET imaging to evaluate the effect of age-induced changes in rat weight on absorbed dose and illustrate the application of the produced S values to the assessment of the corresponding absorbed dose to small animals using the MIRD schema. 11C-labeled 3-amino-4-(2-dimethylaminomethyl-phenylsulfanyl)-benzonitrilev (11C-DASB) is a recently introduced probe for imaging serotonin transporters using PET. 46 According to the reported biodistribution data, internal organs are divided into two groups: organs with time-varying uptake given in Wilson and colleagues 46 (Figure 10A) and organs/tissues in the remaining body (Figure 10B), assuming a uniform uptake assigned according to the residual activity. The TACs of each organ were obtained by fitting the biodistribution data of 11C-DASB in the rat according to equation 5 and used to calculate the absorbed dose to organs in the 10 rat models. The absorbed dose results are illustrated in Figure 10, where the comparisons are made between the different models.

Absorbed doses from 11C-DASB to internal organs of the 10 computational rat models at different ages. 46
Discussion
Monte Carlo techniques are widely used in small-animal radiation dose calculations, in which trustworthy evaluation of radiation dosimetry is grounded on reliable computational models. However, a discrepancy between a computational model and the real anatomy of laboratory specimen commonly occurs, especially when the inanimate computational model is used to mimic a living and growing laboratory animal. Given that laboratory animals have substantially different physical characteristics and present with significantly different radiosensitivity and bioavailability for various radiopharmaceuticals at different age periods, it is recommended to evaluate the age-dependent absorbed dose to small animals.
In this work, we constructed a series of rat models at 10 different age periods and used Monte Carlo calculations to compare the AFs and SAFs of monoenergy photons and electrons and the S values of commonly used positron-emitting radionuclides in various source-target pairs. The self-absorbed AFs are positively correlated with the age and the mass because enlarging the source region results in more energy deposition of particles. The SAFs for all source-target pairs decrease when increasing body mass. Further significant changes are more apparent at low-photon energies, which indicate that the contribution to radiation dose from low-energy photon emitters is more sensitive to the change in weight and may result in a significantly higher absorbed dose in the young rat.
Changes in S values with age and weight were also evaluated for eight positron-emitting radionuclides commonly used in PET molecular imaging–based research. For all source-target pairs, S values decrease when increasing the age and weight. The contribution from the two annihilation photons to the total S values of positron emitters in the whole body is positively correlated with the age-dependent body weight. These results suggest that the activity administered to small animals for PET imaging produces a higher equivalent dose rate in the young rat than in the adult rat. Moreover, given that some tissues of young animals have higher sensitivity to ionizing radiation, this might result in higher radiation risks.30,47–49 Assessment of the internal radiation dose delivered to small animals in experimental research is important for determining the radiotoxicity and therapeutic effect of radiotracers. In this regard, assessment of the variability of S values with changes in weight at different ages is an important issue. Considering F-18 as an example, the difference in S values between the 4-week-old rat model to the 45-week-old rat model for most source-target pairs ranges between −1.5% and −2% per 10 g of difference in body weight.
Based on the calculated S values and biodistribution data gathered from small-animal PET studies, the absorbed doses of 11C-DASB in the 10 developed anatomic rat models were calculated and compared. For all evaluated organs, the calculated absorbed dose decreases with increasing body weight. The difference between absorbed dose in these organs ranges between −1.7% and −2.1% per 10 g difference in body weight from the 4-week-old rat model to the 45-week-old rat model. The above analysis suggests that for all organs, the absorbed dose in low-weight young specimens would be considerably higher than that in adult specimens. This dose-weight/dose-age correlation is also observed in humans, where the absorbed dose per unit activity administered in pediatric patients is significantly higher than that in adult patients. 50
Conclusion
A systematic study was performed to evaluate how dosimetric results vary with changes in age-related body weight. We investigated separately the effects of changes in small-animal weight on the calculation of AFs and SAFs of monoenergy photons/electrons and S values of positron-emitting radionuclides, as well as the organ-absorbed doses for a 11C-labeled probe. For most organs in the rat body, the self-absorbed AF is positively correlated to the weight, whereas the SAF and radionuclide-specific S value present an inverse correlation with body weight. In the rat body, the total S value for positron emitters is mostly contributed by positron interactions. For certain radiotracers, the absorbed dose to organs is negatively correlated with the rat body weight and would be significantly higher in low-weight young animals than in adult animals. Small-animal radiation dosimetry should be evaluated more carefully in multimodality molecular imaging experiments involving the use of young laboratory animals. The developed computational rat models can also be used for external radiation dose calculations. The produced SAFs of monoenergy photons/electrons and S values of positron-emitting radionuclides can be exploited for the assessment of radiation dose delivered by other radionuclides and various radiotracers used in small-animal trials.
Footnotes
Acknowledgment
Financial disclosure of authors: This work was supported by the Swiss National Science Foundation under grant SNSF 31003A-125246, the Geneva Cancer League, and the Indo-Swiss Joint Research Programme ISJRP 138866.
Financial disclosure of reviewers: None reported.
