Individual calculation of effective dose and risk of malignancy based on monte carlo simulations after whole body computed tomography
- Select a language for the TTS:
- UK English Female
- UK English Male
- US English Female
- US English Male
- Australian Female
- Australian Male
- Language selected: (auto detect) - EN

Play all audios:

ABSTRACT Detailed knowledge about radiation exposure is crucial for radiology professionals. The conventional calculation of effective dose (ED) for computed tomography (CT) is based on dose
length product (DLP) and population-based conversion factors (k). This is often imprecise and unable to consider individual patient characteristics. We sought to provide more precise and
individual radiation exposure calculation using image based Monte Carlo simulations (MC) in a heterogeneous patient collective and to compare it to phantom based MC provided from the
National Cancer Institute (NCI) as academic reference. Dose distributions were simulated for 22 patients after whole-body CT during Positron Emission Tomography-CT. Based on MC we calculated
individual Lifetime Attributable Risk (LAR) and Excess Relative Risk (ERR) of cancer mortality. EDMC was compared to EDDLP and EDNCI. EDDLP (13.2 ± 4.5 mSv) was higher compared to EDNCI
(9.8 ± 2.1 mSv) and EDMC (11.6 ± 1.5 mSv). Relative individual differences were up to −48% for EDMC and −44% for EDNCI compared to EDDLP. Matching pair analysis illustrates that young age
and gender are affecting LAR and ERR significantly. Because of these uncertainties in radiation dose assessment automated individual dose and risk estimation would be desirable for dose
monitoring in the future. SIMILAR CONTENT BEING VIEWED BY OTHERS VALIDATION OF SSDE CALCULATION IN A MODERN CT SCANNER AND CORRELATION WITH EFFECTIVE DOSE Article Open access 19 February
2025 CHARACTERIZING IMAGING RADIATION RISK IN A POPULATION OF 8918 PATIENTS WITH RECURRENT IMAGING FOR A BETTER EFFECTIVE DOSE Article Open access 14 March 2024 THE DOSE–RESPONSE
CHARACTERISTICS OF FOUR NTCP MODELS: USING A NOVEL CT-BASED RADIOMIC METHOD TO QUANTIFY RADIATION-INDUCED LUNG DENSITY CHANGES Article Open access 29 June 2020 INTRODUCTION The Euratom
council directive (_2013/59_/_Euratom_) emphasizes the need for patient radiation dose monitoring in clinical routine1. Computed tomography (CT) is indispensable for contemporary patient
care, but many studies suggest a relation between low dose protracted radiation exposure and an increased incidence of malignancy based on the linear Non-threshold Dose-Response Model2.
Therefore, detailed knowledge of radiation exposure from CT examinations is crucial for radiology healthcare professionals in clinical routine to detect and address increased risks of
malignancy. Widely used parameters for radiation exposure assessment like the volumetric CT dose index (CTDIvol) and the Dose Length Product (DLP) characterize scanner radiation output, but
are unable to take individual patient characteristics into account3. Conversion to effective dose (EDDLP) is feasible using population-based conversion factors (k) that take the averaged
radiosensitivity in defined anatomic regions into account4. A variety of different k-factors are recommended in the literature for different volumes (e.g. head, thorax, abdomen, pelvis) and
different CT-scanners5,6. These k-factors are mostly derived from phantom models that try to represent average patient anatomy in western populations, but are unable to respect individual
anatomy like missing organs due to aplasia or resection, organ hypo- or hypertrophy, skeletal deformations and metal implants. Nevertheless, these programs, like e.g. the National Cancer
Institute (NCI) dosimetry system for CT, can be considered as current academic reference (EDNCI)7. However, the routinely performed calculation of EDDLP is imprecise and the degree of
difference to the real ED in each individual patient is unknown. Risk estimates that are calculated on these imprecise data could then easily be misinterpreted. Therefore, conventionally
calculated EDDLP can be used for comparison between different CT-scanners or different examination protocols, but should not be used for comparison between different patient collectives or
for radiation risk estimation8,9. Stochastic estimation of radiation exposure from CT is feasible by Monte Carlo (MC) simulations. Romanyukha _et al_. used computational phantoms and
MC-simulations to calculate body size-specific conversion factors from regression curves10. Huda _et al_. were able to show that this phantom based technique is also feasible to derive organ
doses and to use these for cancer risk estimation11. Only few studies with small numbers of mostly pediatric patients used MC simulations for individual dose and risk estimation per
patient, probably due to the high computational power needed12,13. The main drawback of MC calculations is that they are limited to the reconstructed body volume. Therefore, overexposure in
z-axis and scattered radiation to the organs outside the directly exposed volume cannot be taken into account. Detailed information about biological effects of low-level ionizing radiation
is available from the National Academy of Sciences report number seven about _Health Risks from Exposure to Low Levels of Ionizing Radiation_ (BEIR VII). The authors of BEIR VII report
follow the linear Non-threshold Dose-Response Model and presume that the risk of cancer incidence and mortality due to medical examinations are mainly based on low-level ionizing radiation
dose, age and gender14. The aim of this study was to provide individual calculations of ED in whole-body exposure derived from MC calculations (EDMC) and to compare them with the
conventional method (EDDLP) and the academic reference (EDNCI) following an equivalence hypothesis in a heterogeneous study collective. Furthermore, cancer risk estimates that are based on
EDMC should allow to evaluate the impact of patient characteristics on individual risk of malignancy. Whole body CT from combined Positron Emission Tomography (PET)-CT examinations were
exclusively selected in order to minimize the indeterminate radiation dose to body parts outside the imaging volume that cannot be assessed by MC. METHODS STUDY DESIGN Twenty-two patients
from a collective of 34 consecutive patients, with clinical indication for whole body PET-CT, were retrospectively included (Fig. 1). Inclusion criteria were the complete and artifact-free
coverage of the scan volume, sufficient intra-venous contrast injection and full availability of individual patient data including age, height, weight and body mass index (BMI). Exclusion
criteria were severe artifacts (n = 0) and examinations without contrast injection (n = 12). Clinical indications were infectious diseases (n = 8), bronchial carcinoma (n = 6), head and neck
cancer (n = 3), cancer of unknown primary (n = 2), malignant melanoma (n = 2) and retroperitoneal fibrosis (n = 1). Written informed consent was obtained and archived for each patient. The
study was approved by the local ethics committee of the Friedrich-Alexander University Erlangen-Nuremberg and complied with the Declaration of Helsinki. IMAGE ACQUISITION All examinations
were performed with continuous volume acquisition using a Siemens Biograph True Point 64 PET-CT (Siemens Healthcare GmbH, Forchheim, Germany). Indication, supervision and image analysis were
performed by senior physicians with extensive experience in PET-CT (>6 years). CT was performed from head to mid-femur by a single spiral acquisition in a cranio-caudal direction using a
tube voltage of 120 kV, anatomic tube current modulation with a reference tube-current time product of 170 mAs, pitch 0.8, gantry rotation time 0.5 seconds and 64 × 0.6 mm slice acquisition
by 32 × 0.6 mm detector collimation and z-axis flying focal spot double sampling15. Unlike in diagnostic CT, where different body parts (e.g. head/neck and thorax/abdomen) are often
examined in multiple contrast phases and different positions, PET-CT uses single spiral acquisitions with only one single bolus injection allowing for complete coverage of all body organs
without overlapping of anatomic regions. Therefore, the effect of scattered radiation to all radiosensitive tissues can be taken into account by MC calculations and the effect of overranging
in the beginning and end of the spiral acquisition can be neglected16. All patients received body weight adapted intravenous injection of 18F-fluorodeoxyglucose (3 MBq/kg). After an average
resting period of 91 (±25.6) minutes iodine containing contrast agent (Ultravist 370®, Bayer-Schering Healthcare, Berlin, Germany) was mechanically administered into an antebrachial cannula
using a power injector (Accutron CT-D, Medtron AG, Saarbruecken, Germany), and the single spiral CT acquisition was started with a fixed delay of 70 s. The PET acquisitions were run
subsequently in 7–8 steps, depending on patient length and exam volumes. For this study full field of view (500 × 500 mm) images were reconstructed using a smooth filtered back projection
kernel (B31), slice thickness 5.0 mm and increment 5.0 mm. The resulting voxel size was 0.98 × 0.98 × 5.00 mm. CONVENTIONAL CALCULATION OF EFFECTIVE DOSE Tube current time product, CTDIvol
and DLP were recorded for each examination from the PET-CT scan protocol as provided by the scanner. Details about their definition and calculation are reported elsewhere in literature6. The
main drawback of these monitoring techniques is that they are unable to provide information about the biological impact of radiation exposure and should only be taken as an index of
radiation output by the CT system for comparative purposes9,17. The biologically relevant effective dose (EDDLP) was calculated for each examination by multiplication of DLP and the
region-specific conversion factor. For this study kBody (0.015 mSv/mGy∙cm) was used as recommended by the American Association of Physicists in Medicine6. Dose exposition related to
18F-fluorodeoxyglucose (FDG)-PET were not considered for ED calculation. PHANTOM BASED CALCULATION OF EFFECTIVE DOSE Academic reference ED was calculated using a dedicated CT dosimetry
software tool provided by the National Cancer Institute (NCI-CT)7, which combines reference phantoms provided by the International Commission on Radiological Protection (ICRP) and MC
simulations of a reference CT scanner (Somatom Sensation 16, Siemens Healthcare GmbH, Forchheim, Germany). Individual examination parameters, gender, age, height and body weight were used as
input for this mathematical phantom based calculation (EDNCI). Dedicated descriptions about computational methods of this software are available in the literature7. ORGAN SEGMENTATION All
radiosensitive organs or tissues including all remainders, which are mentioned in IRCP report 103 were segmented using a dedicated software package (ITK-SNAP 1.8, Penn Image Computing and
Science Laboratory, Philadelphia, USA). Semi-automatic threshold segmentation with manual corrections was used for the lungs, cortical bone, bone marrow, liver, spleen, brain, heart, muscle
and skin. All other organs or tissues were segmented manually in a slice per slice fashion. Therefore, missing organs due to aplasia or resection, organ hypo- or hypertrophy, skeletal
deformations and metal implants were considered. For hollow organs, only voxels contributing to the wall were considered to contain radiosensitive tissues. Segmentation of lymphatic tissues
was restricted to visible lymphatic nodes (Fig. 2). All segmentations were performed by a specially trained radiology resident (4 years experience) and reviewed from an experienced
board-certified radiologist (9 years experience). MONTE CARLO BASED CALCULATION OF EFFECTIVE DOSE Extensive mathematical models in simulation techniques can be utilized to calculate dose
distribution for each voxel of CT scan volumes using the attenuation values and geometry from DICOM datasets as input18,19. All MC dose simulations were carried out using the software
package ImpactMC (VAMP GmbH, Erlangen, Germany). Software details and information concerning the Monte Carlo calculation algorithms are reported elsewhere18,20. Multiplication of the
resulting relative dose values per voxel with air kerma provides the absorbed dose for each voxel. The air kerma is a scanner specific value, which describes the kinetic energy transferred
in air dependent on geometry and tube settings. It can be considered as measurement of x-ray beam intensity without object, which was 18.3 mGy for this study. The averaged absorbed dose over
an organ volume, defined from segmentations, provides the absorbed organ dose. Equivalent organ dose was calculated by multiplication with the radiation weighting factor for photons (WR =
1) and used for subsequent effective organ dose calculations (EDOrgan) by multiplication with the dedicated tissue weighting factor (WT) as recommended in ICRP report 10321. The sum of all
effective organ doses provides EDMC for each individual (Fig. 3). The relative contribution of the respective EDOrgan to the EDMC was calculated following Eq. 1:
$${{\rm{C}}}_{{\rm{Organ}}}={{\rm{ED}}}_{{\rm{Organ}}}/{{\rm{ED}}}_{{\rm{MC}}}$$ (1) Individual conversion factors (kMC) were calculated using Eq. 2 to illustrate the differences compared to
the conventional calculation method using DLP and kbody. $${{\rm{k}}}_{{\rm{MC}}}={{\rm{ED}}}_{{\rm{MC}}}/{\rm{DLP}}$$ (2) RADIATION RISK ASSESSMENT Lifetime Attributable Risk (LAR) is
defined as risk of disease of an exposed cohort in comparison to a non-exposed cohort. Calculations of LAR for cancer mortality were based on report VII about Biologic Effects of Ionizing
Radiation (BEIR VII). The BEIRV VII report refers to an individual radiation dose exposure of 100 mGy. Assuming a linear risk distribution between decennium age intervals LARMC was derived
from EDMC by linear interpolation between the younger (Ny) and older (NO) cohort as shown in equation 3 using patient age (Ap) and the age of the younger cohort (Ay) as input values22.
Equation 3: Calculation of individual LAR $$LA{R}_{MC}=\left({N}_{y}-\left(({N}_{y}-{{\rm{N}}}_{0})\ast \frac{{A}_{p}-{A}_{Y}}{10\,{\rm{years}}}\right)\right)\ast
\frac{E{D}_{Organ}}{100\,{\rm{mGy}}}$$ The Excess Relative Risk (ERRMC), as a measurement of the exceeding risk of an exposed person compared to a non-exposed person, was calculated using
the solid cancer mortality in the United States as baseline (female: 17500/100000; male: 22100/100000) using equation 423. Equation 4: Calculation of individual ERR
$$ER{R}_{MC}=LA{R}_{MC}/LA{R}_{baseline}$$ PAIRWISE PATIENT COMPARISON Because of the rather small patient collective in this study, mainly due to the extensive effort needed for
segmentation of all radiosensitive tissues, pairwise patient comparisons were selected to highlight the influence of individual patient characteristics on the radiation dose parameters and
radiation risk estimation. Matching patient pairs, each with two comparable values and one variable parameter, were found for the parameters age, sex and BMI. STATISTICAL ANALYSIS All
statistical analyses were performed using the software package SPSS Statistics Version 21 (IBM, Somers, NY, USA). Normal distribution of the data was tested by Kolmogorov-Smirnov and
Shapiro-Wilk test. Normally distributed data is presented as mean ± standard deviation. Median and range are provided if no normal distribution was assumed. Illustration is provided as
Bland-Altman plots. Spearman’s rank order test was used to test for correlations between tube current, BMI, EDDLP and EDMC. The significance level was defined as p < 0.05. RESULTS
DEMOGRAPHICS OF THE STUDY COLLECTIVE Five out of 22 patients (23%) were female and 17 male (77%). The mean age was 57.3 ± 14.3 years and the mean BMI was 26.0 ± 5.9 kg/m². Three patients
(13.6%) were younger than 40 years. The female patient group was younger (53.6 ± 17.6 years) compared to the male patient group (58.4 ± 13.6 years). Female patients had a lower mean BMI
(21.2 ± 2.5 kg/m²) compared to male patients (27.4 ± 5.88 kg/m2). From the male subgroup 9 patients (52.9%) suffered from overweight (BMI > 25), and three patients (17.6%) had severe
overweight (BMI > 30). One female and one male patient were considered as underweight (BMI < 19). No overweight patient was found in the female subgroup. Detailed patient
characteristics are shown in Table 1. ORGAN SPECIFIC RADIATION DOSE EXPOSURE Mean equivalent organ dose was 13.0 ± 3.5 mGy. Highest equivalent organ dose values were found in high
attenuating structures such as the cortical bone and the thyroid gland. Lowest values were found in profound organs like the extra-thoracic (ET) respiratory region and the uterus. Many
superficial organs like muscles (11.5 ± 2.0 mGy), breast tissue (10.9 ± 0.8 mGy) and salivary glands (14.3 ± 2.8 mGy) also had rather low equivalent organ dose values. Despite its profound
position well perfused kidneys had very high equivalent organ dose values (17.0 ± 2.1 mGy). Effective organ dose was mainly influenced by WT, nicely demonstrated by the thyroid gland. It has
the second highest mean equivalent organ dose (21.2 ± 5.4 mGy), but is only ranked 7th highest effective organ dose (8% of total organ dose distribution) because of the rather low WT
(0.04). In contrast, the lungs ranked only 7th in equivalent organ dose (13.4 ± 1.9 mGy), but received highest effective organ dose levels in female and male patients (12% and 14% of the
total EDMC) due to high WT (0.12). Colon, lungs and stomach contributed to more than 50% of the total EDMC in male patients and to 48% in female patients (Fig. 4). Differences between WT and
the relative contribution of the effective organ dose to whole body EDMC (COrgan) reflect the influence of individual radiation dose distribution on each organ dose. Highest positive
differences between WT and COrgan were found for the bone surface (♀: +102.1%; ♂: +134.8%), the thyroid gland (♀: +89.1%; ♂: +78.1%) and the kidneys (♀: + 49.4%; ♂: +47.1%). The male gonads
also had substantially higher COrgan than expected from its rather low WT ( + 38.2%), but the female gonads had the highest negative difference of all organs (−21.8%). Other organs with high
negative differences were the adrenal glands (♀: −18.3%; ♂: −11.8%). Detailed organ dose information is shown in Table 2 and organ-based LAR of cancer mortality is provided for all organs
listed in BEIR VII. Considering the similar effective organ doses of the lung (1.48 ± 0.15 mSv) and the breast (1.30 ± 0.09 mSv) the LAR for pulmonary malignant disease (13.25 ±
4.24/100.000) was remarkably higher compared to the LAR for breast cancer mortality (2.46 ± 1.62/100.000). EFFECTIVE DOSE AND INDIVIDUAL RADIATION RISK ASSESSMENT EDDLP (13.2 ± 4.52 mSv) was
higher than EDMC (11.6 ± 1.47 mSv) and EDNCI (9.8 ± 2.1 mSv). Particularly high differences between EDMC and EDDLP were found for patients with relatively high or low radiation dose
exposure. Mean difference between both methods was −1.7 mSv (Fig. 5a). The same tendency was found for the comparison of EDNCI and EDDLP but with higher mean negative differences −3.4 mSv
(Fig. 5b). A comparison between the three calculation methods is illustrated as boxplot (Fig. 6). The range of radiation dose was substantially smaller in both advanced methods (EDMC: 5.6
mSv, EDNCI: 10.2 mSv) compared to the conventional technique (EDDLP: 19.3 mSv). Effective mAs and DLP linearly increased with higher BMI in Spearman’s rank order test (rs = 0.961 and 0.949).
This effect should be mainly due to the anatomy-based tube current modulation algorithm. Therefore, also BMI and EDDLP had a high correlation coefficient (rs = 0.949). The correlation
between EDMC and BMI (rs = 0.644) was less and differences between EDDLP and the advanced techniques was especially high in over- and underweight patients. For example, the highest EDDLP of
26.3 mSv was found in a 77-year-old man suffering from adiposity grade 3 (BMI 40.0 kg/m²) while EDNCI (17.3 mSv) and EDMC (13.9 mSv) were substantially lower (34% and 47% less). Lowest EDDLP
of 7.0 mSv was calculated for a 30-year-old underweight man (BMI 16.8 kg/m²), which was close to EDNCI (7.1 mSv), but substantially lower than EDMC (8.5 mSv, 21% higher). Consequently, the
individually calculated conversion factors (kMC) have a range from 53% to 140% when referred to kBody from literature. The mean kMC (0.014 ± 0.004 mSv/mGy cm) approximately reflects the
established value (0.015 mSv/mGy cm), which therefore seems to be suitable for regular weight patients. However, only one fourth of the patients in this study (n = 6, 27.2%) had less than
10% difference between kMC (0.0135–0.0165 mSv/mGy cm) and kBody (mean BMI 22.95 ± 1.58, range 20.7–25.0 kg/m²). The values of LAR were especially high in young female patients. The highest
value (60/100.000) was calculated for a 30-year-old, normal weight woman (BMI 21.6 kg/m²). In contrast, low LAR values were calculated for older, predominantly male patients. The lowest
value (23/100.000) was found for a 76-year-old, overweight man (BMI 29.8 kg/m²). Extent of EDMC, EDDLP and LAR often differed considerably when compared between young and old or underweight
and obese patients, while the interaction between these individual parameters was rather difficult to predict. Detailed results of radiation dose and risk calculations are provided in Table
3. PAIRWISE PATIENT COMPARISON Two male patients with matching constitution (BMI 24 vs. 23 kg/m²) but different age (36 versus 67 years) had comparable DLP (783 vs. 724 mGy · cm, −7.5%),
EDDLP (11.8 vs. 10.9 mSv, −7.5%), EDNCI (9.4 vs. 9.0 mSv, −4.3%) and EDMC (10.9 vs 10.0 mSv, −8.4%). Nevertheless, calculations for LAR of cancer mortality (41.3 vs. 27.1/100.000; −34.4%)
and ERR (0.19 vs. 0.12%, −36.8%) differed considerably. The young age seems to be the decisive factor among all the considered factors, accounting for about 50% higher risk estimates (Fig.
7a). Two male patients with matching age (76 and 77 years) but different BMI (29.7 vs. 40.0 kg/m², +34.7%) were exposed to considerably different DLP (998 vs. 1751 mGy·cm, +75.4%) based on
anatomic tube current modulation. Therefore, EDDLP calculations provided very high values for the obese patient (15.0 vs. 26.3 mSv, +75.5%). Dose differences calculated with EDNCI (10.5 vs.
17.3 mSv, +64.5%) were also considerably high, while EDMC values ware nearly comparable (11.9 vs. 13.8 mSv; +15.9%). Differences in LARMC (22.8 vs. 25.1/100.000; +9.6%) and ERR (0.10 vs
0.11%; +10.0%) were even less (Fig. 7b). BMI seems to have only a minor influence on effective dose and risk estimation, while conventional methods would have resulted in an overestimation
of almost 100%. Two patients with matching age (56 and 54 years) and BMI (21.5 vs 21.0 kg/m²) but different sex had similar DLP (634 vs. 630 mGy·cm; +0.63%). ED was the same using the
conventional method in both patients (EDDLP = 9.5 mSv) and almost the same with the NCI method (EDNCI = 8.7 vs 9.0 mSv; −3.3%). EDMC calculation was slightly higher for the female patient
(12.3 vs. 10.3 mSv; +19.4%), but LAR (53.1 vs. 35.4/100000; +50.2%) and ERR (0.3% vs. 0.16%; +87.5%) were substantially higher. Sex seems to have a dominant impact on the risk estimation.
Moreover, the unfavorable combination of age and sex leads to the highest values in the entire study collective, despite average ED values (Fig. 7c). DISCUSSION Individual radiation dose
assessment and risk calculation is feasible by image based Monte Carlo simulations and organ segmentations in an adult clinical routine collective that underwent full body exposure in a
single spiral acquisition. This is in good agreement with the findings of Li _et al_. and Tian _et al_. who evaluated comparable techniques for radiation dose estimation in small pediatric
collectives (n = 2 and n = 42)16,24. A comparable approach for cancer risk estimation in adults based on anthropomorphic phantoms was described by Huda _et al_.11. In contrast to these prior
studies, we provide a combination of individual radiation dose analysis, voxel-based organ segmentation and cancer risk estimation in direct comparison to such phantom based estimation and
the DLP method. The need for such an advanced dose monitoring method, due to several uncertainties with the conventional techniques, has been raised elsewhere in literature20,25. In our
study high equivalent organ doses were strongly related to radiodensity (e.g. bone surface and high contrast medium uptake in the kidneys, the thyroid gland and the extra-thoracic
respiratory region), which seems to be a much stronger predictor for equivalent organ dose than anatomical organ position. This is in good agreement with the contrast media related increase
of DNA double-strand break foci after radiation exposure reported in a prior study for CT26. WT is the strongest predictor of effective organ dose. For example, the thyroid gland and the
kidneys have the second and third highest equivalent organ dose (21.2 ± 3.2 mSv and 17.0 ± 2.1 mSv). However, due to low tissue weighting factors (WT = 0.04 and 0.0092) effective organ dose
of the thyroid gland and the kidneys were only sixth and twelfth highest. Our results confirm that individual patient characteristics have a considerable impact on radiation dose
calculation, which is underrepresented by the conventional method. The calculated error in our adult collective (−48 to 39%) is comparable to the previously published study results for
children (−63 to 28%)19. Individually calculated kMC from our study reveals that the routinely used kbody from literature can be applied to larger clinical collectives, but are limited in
their value for individual risk assessment6. The DLP-method slightly overestimates the effective dose in general, but it seems to be appropriate in regular weight male patients. However,
especially in female and underweight patients, underestimation of the effective dose may be critical for risk perception in the clinical setting. Moreover, its overestimation in obese
patients may lead to restrained use of high exposure parameters, and then to poor or insufficient image quality. Individual anatomic characteristics (e.g. missing organs due to aplasia or
resection, organ hypo- or hypertrophy, skeletal deformations, metal implants) may lead to considerable changes of radiation dose distribution and consequently of individual risk. Even
phantom based estimation methods, like EDNCI in this study, are unable to take these individual properties into account7. Bland-Altman analysis demonstrated that these phantom based
estimates have the tendency to underestimate ED in general, while the effect of only moderate increase of ED in obese patients despite the very high energy exposure is confirmed. The
influence of personal risk profiles on cancer risk estimates are not yet represented in standard radiation dose reporting systems. It can further increase the error of perceived and true
risk of radiation dose from CT examinations. The combination of individual dose distributions and risk assessment has been presented for a few other examinations in literature, such as
absorptiometry and spine radiographs, but not for spiral CT examinations27. Pairwise patient comparisons provide a comprehensive overview of the extent of error with regard to age, BMI and
sex. Especially younger patients are prone to higher LAR and ERR since more active cell division and longer life expectancy after radiation exposure is presumed28. In our example the ERR of
a patient in his mid-thirties was 1.6-fold the ERR of a patient in his mid-sixties and the ERR for a female patient in her mid-fifties was almost doubled compared to a matching male patient.
The increasing tube current by anatomy-based modulation in obese patients seems to be of less importance to the patients’ cancer risk. CTDI and EDDLP were 1.8-fold higher for a high-grade
obese patient compared to an overweight patient with matching age and sex, while the ERR was only 10% higher. Therefore, future dose assessment strategies should not only focus on absolute
dose values. Furthermore, reporting LAR or ERR seems to be much more appropriate in clinical routine and a more comprehensible value for the patient-physician interaction. Some limitations
have to be considered while interpreting this study. First, the study population consisted of a heterogeneous, small and retrospectively selected patient collective with an imbalanced ratio
of female to male patients. The small number of patients is mainly due to the high computing power required for MC calculations and the work intense manual organ segmentation that we used
for this study. We estimate that the increasing computing power and automated organ segmentations by the application of artificial intelligence could overcome this limitation in the near
future29,30. Second, only full body dose exposure is reported in this study. This avoids indeterminate scattered radiation to radiosensitive structures and over-ranging effects, but limits
the findings to this clinically rather rare indication. Further evaluations for limited examination volumes could become feasible with retrospective simulations from these data in larger
collectives by future studies. Third, all examinations were conducted on the same scanner. The findings of this study can therefore not automatically be transferred to other CT systems5,24.
Forth, all tissue weighting factors that are provided in literature so far are averaged for sex and age, which may probably limit their applicability for patient specific risk estimation.
Fifth, the linear Non-threshold Dose-Response Model itself is discussed controversially among experts for diagnostic dose levels31,32. Sixth, LAR and ERR calculations in this study are only
related to low dose radiation exposure. They illustrate the individual radiation dose related risk to develop cancer. Obviously, the overall individual cancer risk for a primary or even a
secondary cancer and also the life expectancy is potentially much more influenced by age, genetic make-up, efficiency of DNA damage repair, therapy-related adverse effects and many other
influencing factors. The importance of patient specific dose surveillance is illustrated by this study. EDDLP can be used for radiation dose assessments in larger collectives, but individual
considerations require advanced techniques like EDMC. Conventional methods tend to underestimate radiation dose in underweight and female patients. Therefore, these patients have to be
evaluated with special care. The influence of young patient age and female gender on cancer risk estimates is very high. Thus, radiation dose assessment should not only provide whole body or
organ dose measurements but also individual risk calculations, which could be included in future surveillance programs. DATA AVAILABILITY The datasets generated during and analyzed during
this study are available from the corresponding author on reasonable request. REFERENCES * EC. COUNCIL DIRECTIVE 2013/59/EURATOM of 5 December 2013 laying down basic safety standards for
protection against the dangers arising from exposure to ionising radiation, a. r. D. * Laurier, D. _et al_. The International Nuclear Workers Study (Inworks): A Collaborative Epidemiological
Study to Improve Knowledge About Health Effects of Protracted Low-Dose Exposure. _Radiation protection dosimetry_ https://doi.org/10.1093/rpd/ncw314 (2016). Article Google Scholar * Huda,
W. & Mettler, F. A. Volume CT dose index and dose-length product displayed during CT: what good are they? _Radiology_ 258, 236–242, https://doi.org/10.1148/radiol.10100297 (2011).
Article PubMed Google Scholar * Blackwell, C. R. & McCullough, E. C. A chamber and electrometer calibration factor as determined by each of the five AAPM accredited dosimetry
calibration laboratories. _Medical physics_ 19, 207–208, https://doi.org/10.1118/1.596880 (1992). Article ADS CAS PubMed Google Scholar * Hill, K. D. _et al_. Radiation Safety in
Children With Congenital and Acquired Heart Disease: A Scientific Position Statement on Multimodality Dose Optimization From the Image Gently Alliance. _JACC. Cardiovascular imaging_ 10,
797–818, https://doi.org/10.1016/j.jcmg.2017.04.003 (2017). Article PubMed PubMed Central Google Scholar * AAPM. The Measurement, reporting, and management of radiation dose in CT,.
_AAPM Report No. 96_ (2008). * Lee, C., Kim, K. P., Bolch, W. E., Moroz, B. E. & Folio, L. NCICT: a computational solution to estimate organ doses for pediatric and adult patients
undergoing CT scans. _Journal of radiological protection: official journal of the Society for Radiological Protection_ 35, 891–909, https://doi.org/10.1088/0952-4746/35/4/891 (2015). Article
ADS CAS Google Scholar * Martin, C. J. Effective dose: how should it be applied to medical exposures? _The British journal of radiology_ 80, 639–647,
https://doi.org/10.1259/bjr/25922439 (2007). Article CAS PubMed Google Scholar * McCollough, C. H., Christner, J. A. & Kofler, J. M. How effective is effective dose as a predictor of
radiation risk? _AJR. American journal of roentgenology_ 194, 890–896, https://doi.org/10.2214/AJR.09.4179 (2010). Article PubMed Google Scholar * Romanyukha, A., Folio, L., Lamart, S.,
Simon, S. L. & Lee, C. Body Size-Specific Effective Dose Conversion Coefficients for Ct Scans. _Radiation protection dosimetry_ 172, 428–437, https://doi.org/10.1093/rpd/ncv511 (2016).
Article CAS PubMed Google Scholar * Huda, W. & He, W. Estimating cancer risks to adults undergoing body CT examinations. _Radiation protection dosimetry_ 150, 168–179,
https://doi.org/10.1093/rpd/ncr376 (2012). Article PubMed Google Scholar * Li, X. _et al_. Patient-specific dose estimation for pediatric chest CT. _Medical physics_ 35, 5821–5828,
https://doi.org/10.1118/1.3026593 (2008). Article ADS PubMed PubMed Central Google Scholar * Kost, S. D. _et al_. Patient-specific dose calculations for pediatric CT of the chest,
abdomen and pelvis. _Pediatric radiology_ 45, 1771–1780, https://doi.org/10.1007/s00247-015-3400-2 (2015). Article PubMed PubMed Central Google Scholar * In Health Effects of Exposure to
Low Levels of Ionizing Radiations: Time for Reassessment? (1998). * Kyriakou, Y., Kachelriess, M., Knaup, M., Krause, J. U. & Kalender, W. A. Impact of the z-flying focal spot on
resolution and artifact behavior for a 64-slice spiral CT scanner. _European radiology_ 16, 1206–1215, https://doi.org/10.1007/s00330-005-0118-9 (2006). Article PubMed Google Scholar *
Li, X. _et al_. Patient-specific radiation dose and cancer risk estimation in CT: part II. Application to patients. _Medical physics_ 38, 408–419, https://doi.org/10.1118/1.3515864 (2011).
Article ADS PubMed Google Scholar * McCollough, C. H. _et al_. CT dose index and patient dose: they are not the same thing. _Radiology_ 259, 311–316,
https://doi.org/10.1148/radiol.11101800 (2011). Article PubMed PubMed Central Google Scholar * Deak, P., van Straten, M., Shrimpton, P. C., Zankl, M. & Kalender, W. A. Validation of
a Monte Carlo tool for patient-specific dose simulations in multi-slice computed tomography. _European radiology_ 18, 759–772, https://doi.org/10.1007/s00330-007-0815-7 (2008). Article
PubMed Google Scholar * DeMarco, J. J. _et al_. A Monte Carlo based method to estimate radiation dose from multidetector CT (MDCT): cylindrical and anthropomorphic phantoms. _Physics in
medicine and biology_ 50, 3989–4004, https://doi.org/10.1088/0031-9155/50/17/005 (2005). Article ADS CAS PubMed Google Scholar * Jarry, G., DeMarco, J. J., Beifuss, U., Cagnon, C. H.
& McNitt-Gray, M. F. A Monte Carlo-based method to estimate radiation dose from spiral CT: from phantom testing to patient-specific models. _Physics in medicine and biology_ 48,
2645–2663 (2003). Article ADS CAS Google Scholar * The 2007 Recommendations of the International Commission on Radiological Protection. ICRP publication 103. _Annals of the ICRP_ 37,
1–332, https://doi.org/10.1016/j.icrp.2007.10.003 (2007). * Rampinelli, C. _et al_. Exposure to low dose computed tomography for lung cancer screening and risk of cancer: secondary analysis
of trial data and risk-benefit analysis. _Bmj_ 356, j347, https://doi.org/10.1136/bmj.j347 (2017). Article PubMed PubMed Central Google Scholar * Lee, W. C. Excess relative risk as an
effect measure in case-control studies of rare diseases. _PloS one_ 10, e0121141, https://doi.org/10.1371/journal.pone.0121141 (2014). Article CAS PubMed Google Scholar * Tian, X. _et
al_. Pediatric chest and abdominopelvic CT: organ dose estimation based on 42 patient models. _Radiology_ 270, 535–547, https://doi.org/10.1148/radiol.13122617 (2014). Article PubMed
PubMed Central Google Scholar * Jansen, J. T. & Shrimpton, P. C. Development of Monte Carlo simulations to provide scanner-specific organ dose coefficients for contemporary CT.
_Physics in medicine and biology_ 61, 5356–5377, https://doi.org/10.1088/0031-9155/61/14/5356 (2016). Article ADS CAS PubMed Google Scholar * Grudzenski, S., Kuefner, M. A., Heckmann,
M. B., Uder, M. & Lobrich, M. Contrast medium-enhanced radiation damage caused by CT examinations. _Radiology_ 253, 706–714, https://doi.org/10.1148/radiol.2533090468 (2009). Article
PubMed Google Scholar * Law, M. _et al_. Cumulative Effective Dose and Cancer Risk of Pediatric Population in Repetitive Whole-Body Scan Using Dual-Energy X-Ray Absorptiometry. _Journal of
clinical densitometry: the official journal of the International Society for Clinical Densitometry_, https://doi.org/10.1016/j.jocd.2017.09.005 (2017). * Kalra, M. K., Sodickson, A. D.
& Mayo-Smith, W. W. CT Radiation: Key Concepts for Gentle and Wise Use. _Radiographics: a review publication of the Radiological Society of North America, Inc_ 35, 1706–1721,
https://doi.org/10.1148/rg.2015150118 (2015). Article Google Scholar * McBee, M. P. _et al_. Deep Learning in Radiology. _Academic radiology_, https://doi.org/10.1016/j.acra.2018.02.018
(2018). * Lee, H. _et al_. Pixel-Level Deep Segmentation: Artificial Intelligence Quantifies Muscle on Computed Tomography for Body Morphometric Analysis. _Journal of digital imaging_ 30,
487–498, https://doi.org/10.1007/s10278-017-9988-z (2017). Article PubMed PubMed Central Google Scholar * Shore, R. _et al_. Implications of recent epidemiologic studies for the linear
nonthreshold model and radiation protection. _Journal of radiological protection: official journal of the Society for Radiological Protection_, https://doi.org/10.1088/1361-6498/aad348
(2018). * Kase, K. R. Radiation protection principles of NCRP. _Health physics_ 87, 251–257 (2004). Article CAS Google Scholar Download references AUTHOR INFORMATION AUTHORS AND
AFFILIATIONS * Department of Radiology, University Hospital Erlangen, Erlangen, Germany Markus Kopp, Wolfgang Wuest, Michael Brand, Matthias Wetzl, Wolfram Nitsch, Michael Uder &
Matthias May * Institute of Medical Microbiology and Hygiene, University of Technology Dresden, Dresden, Germany Tobias Loewe * Department of Nuclear Medicine, University Hospital Erlangen,
Erlangen, Germany Daniela Schmidt & Michael Beck * Siemens Healthineers GmbH, Forchheim, Germany Bernhard Schmidt Authors * Markus Kopp View author publications You can also search for
this author inPubMed Google Scholar * Tobias Loewe View author publications You can also search for this author inPubMed Google Scholar * Wolfgang Wuest View author publications You can also
search for this author inPubMed Google Scholar * Michael Brand View author publications You can also search for this author inPubMed Google Scholar * Matthias Wetzl View author publications
You can also search for this author inPubMed Google Scholar * Wolfram Nitsch View author publications You can also search for this author inPubMed Google Scholar * Daniela Schmidt View
author publications You can also search for this author inPubMed Google Scholar * Michael Beck View author publications You can also search for this author inPubMed Google Scholar * Bernhard
Schmidt View author publications You can also search for this author inPubMed Google Scholar * Michael Uder View author publications You can also search for this author inPubMed Google
Scholar * Matthias May View author publications You can also search for this author inPubMed Google Scholar CONTRIBUTIONS Contributed to the writing of the manuscript: M.K., M.M., T.L. and
W.W. Study design, measurements and statistical analysis: M.K., M.M., W.W., M.Br., M.W., T.L., W.N., B.S., D.S., M.Be. and M.U. All authors reviewed the manuscript. CORRESPONDING AUTHOR
Correspondence to Markus Kopp. ETHICS DECLARATIONS COMPETING INTERESTS M.K., M.M., W.W. and M.U. are members of the Siemens Speaker’s Bureau. B.S. is an employee of Siemens Healthineers
GmbH. M.W., T.L., M.Br., W.N., D.S. and M.Be. declare no competing financial interests. All authors declare no competing non-financial 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 licensed under a
Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are
included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and
your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this
license, visit http://creativecommons.org/licenses/by/4.0/. Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Kopp, M., Loewe, T., Wuest, W. _et al._ Individual Calculation of
Effective Dose and Risk of Malignancy Based on Monte Carlo Simulations after Whole Body Computed Tomography. _Sci Rep_ 10, 9475 (2020). https://doi.org/10.1038/s41598-020-66366-2 Download
citation * Received: 21 January 2019 * Accepted: 14 May 2020 * Published: 11 June 2020 * DOI: https://doi.org/10.1038/s41598-020-66366-2 SHARE THIS ARTICLE Anyone you share the following
link with will be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the Springer Nature
SharedIt content-sharing initiative