Purpose:

Radiotherapy for head and neck cancer may result in serious side effects, such as hyposalivation, impairing the patient's quality of life. Modern radiotherapy techniques attempt to reduce the dose to salivary glands, which, however, results in low-dose irradiation of the tissue stem cells. Here we assess the low-dose sensitivity of tissue stem cells and the consequences for tissue function.

Experimental Design:

Postirradiation rat salivary gland secretory function was determined after pilocarpine induction. Murine and patient-derived salivary gland and thyroid gland organoids were irradiated and clonogenic survival was assessed. The DNA damage response (DDR) was analyzed in organoids and modulated using different radiation modalities, chemical inhibition, and genetic modification.

Results:

Relative low-dose irradiation to the high-density stem cell region of rat salivary gland disproportionally impaired function. Hyper-radiosensitivity at doses <1 Gy, followed by relative radioresistance at doses ≥1 Gy, was observed in salivary gland and thyroid gland organoid cultures. DDR modulation resulted in diminished, or even abrogated, relative radioresistance. Furthermore, inhibition of the DDR protein ATM impaired DNA repair after 1 Gy, but not 0.25 Gy. Irradiation of patient-derived salivary gland organoid cells showed similar responses, whereas a single 1 Gy dose to salivary gland–derived stem cells resulted in greater survival than clinically relevant fractionated doses of 4 × 0.25 Gy.

Conclusions:

We show that murine and human glandular tissue stem cells exhibit a dose threshold in DDR activation, resulting in low-dose hyper-radiosensitivity, with clinical implications in radiotherapy treatment planning. Furthermore, our results from patient-derived organoids highlight the potential of organoids to study normal tissue responses to radiation.

Translational Relevance

Radiotherapy for cancer treatment may result in serious side effects due to coirradiation of normal tissue. In the case of head and neck cancers, a major complication is xerostomia caused by damage to the salivary glands. Optimizing radiotherapy planning using modern radiotherapy techniques results in increasing the volume of the salivary glands receiving relatively low doses per fraction. We show that salivary glands over respond to low doses with a strong reduction in secretory potential. Using a 3D organoid model we show that this is due to low-dose radiation sensitivity of murine and patient stem cells. Moreover, the data suggest that fractionated doses of less than 1 Gy may induce clinically relevant normal tissue effects. This study indicates the importance of sparing glandular tissue stem cells during radiotherapy even from relatively low radiation doses.

Currently, more than half of all patients with cancer are treated with radiotherapy. Unfortunately, radiotherapy is often accompanied by severe side effects induced by unavoidable coirradiation of surrounding normal tissues (1, 2). Even modern radiotherapy of patients with head and neck cancer (HNC) still results in functional impairment of coirradiated salivary glands, consequential xerostomia (dry mouth syndrome), alterations in speech and taste, difficulties with mastication and deglutition, dental caries, and oral infections (3, 4), which have a major impact on patient-reported quality of life (5, 6). Currently, no effective treatment exists for radiation-induced side effects in patients with HNC. Radiation-induced normal tissue damage critically depends on the stem cell response, which is essential for the preservation of homeostasis and function of normal tissue after irradiation (7–9). High-precision techniques, such as intensity-modulated radiotherapy (IMRT) and particle therapy (2), could be used to deliver the radiation dose to the tumor target while sparing the high-density stem cell region (8). However, these technologies differ greatly in the level to which dose to the high-density stem cell region can be reduced. Although IMRT can significantly reduce the mean dose to the gland, large volumes still receive a dose below 1 Gy per fraction.

Interestingly, there are several indications that cells can be disproportionately sensitive to radiation doses below 1 Gy, and that such low doses to specific regions of tissues may have a large impact on function and survival. Precisely targeted proton irradiation experiments on the rat cervical spinal cord indicated that besides the size of the irradiated volume, also the dose distribution determines the chance of paralysis of the fore and hind limbs (10–12). Importantly, very low doses to a region adjacent to a high-dose volume were shown to markedly reduce the tolerance of the spinal cord (10–12). Furthermore, a dose of only 1 Gy to the cranial region of the rat parotid gland had a greater effect on secretory function than higher doses of 2–5 Gy to the same region (13), indicating a nonmonotonous response with increasing radiation doses. This deviates from the expected monotonously decreasing dose response based on classic cell survival responses (14, 15).

Increased sensitivity to relatively low doses of radiation has been observed in in vitro survival assays of several, but not all, tumor cell lines, whereas the dose at which this disproportional sensitivity has been observed can vary depending on cell type, culture conditions, and endpoint studied (16, 17). Although the exact mechanism is not clear, some studies suggested the presence of a dose threshold in initiation of the DDR (18). So far, comparable observations for normal tissue cells are primarily limited to 2D in vitro survival curves of keratinocytes and fibroblasts (19). However, these cells are not stem cells responsible for long-term homeostasis and regeneration after radiation damage (7–9, 20). Knowledge of the effects of low-dose irradiation on normal tissue stem cells is presently restricted to the hematopoietic system, where Rodrigues-Moreira and colleagues showed that low-dose irradiation induces persistent oxidative stress and results in long-term detrimental effects in hematopoietic stem cells (21). However, whether other normal tissue stem cells show high sensitivity to low doses is currently unknown.

Recent developments in the culture of organoids, three-dimensional tissue resembling structures derived from adult stem cells (9, 22), have opened unique opportunities to study normal tissue development and pathologic responses (23). These structures contain highly potent self-renewing stem cells that differentiate to all tissue lineages and are able to regenerate many tissues (9, 24, 25). Moreover, they enable the in vitro study of normal tissue stem cell responses to radiation (26). In this study, we use organoids to investigate the effects of low-dose irradiation on normal tissue stem cells. First, we show that the impact of a low dose of radiation on rat parotid gland saliva secretion may be greater than that of higher doses. We demonstrate that this response is mimicked by, but not limited to, murine salivary gland stem cells in 3D organoid cultures (9, 26–28) and is due to insufficient DDR activation at lower doses. Next, we show that human salivary gland stem cells also isolated from patient biopsies show this response. Moreover, we demonstrate that low-dose radiation sensitivity is not limited to a single fraction, but accumulates after several fractions resulting in a pronounced additive response compared with that of an equal single fraction of a higher dose. Our data suggest that it may be important to take low-dose sensitivity of stem cells into account in the selection of radiotherapy technology and treatment plan optimization.

Isolation of murine salivary and thyroid gland cells

Murine salivary glands were dissected from 8- to 12-week-old female C57BL/6 mice. For experiments with cells from Atm/ mice, salivary glands dissected from 12- to 16-week-old Atm/ (129/SV × C57BL/6) mice were used. Generation and phenotyping of Atm/ mice were described previously (29, 30). Salivary gland cells were isolated and cultured to form organoids as described previously (27, 28, 31).

Murine thyroid glands were dissected from 8- to 12-week-old female C57BL/6 mice. Thyroid gland tissue derived from 3 mice was collected in one tube in Hank balanced salt solution (HBSS) containing 1% BSA and was mechanically digested using the GentleMACS dissociator (Miltenyi Biotec) followed by enzymatic digestion in HBSS/1% BSA buffer containing collagenase I (100 U/mL; Gibco) and dispase (1.5 U/mL; Gibco). Cells were seeded into one well of a 12-well plate in DMEM-F12 medium (Gibco) containing 1% penicillin/streptomycin (Gibco), glutamax (2 mmol/L, Gibco), EGF (20 ng/mL, Sigma Aldrich), fibroblast growth factor-2 (20 ng/mL, PeproTech), and 0.5% B27 supplement (Gibco). After 1 day, primary organoids (passage 0) were disrupted using 0.05% trypsin-EDTA (Gibco) into single cells and reseeded at a density of 2 × 104 cells/well in Matrigel (BD Biosciences) and medium supplemented with ROCK1-inhibitor (Y-27632; Sigma-Aldrich) was added to the gels.

All experiments were approved by the Ethical Committee on animal testing of the University of Groningen (Groningen, the Netherlands).

Human salivary gland isolation and culturing

This study was conducted in accordance with the Declaration of Helsinki, Ethical Principles for Medical Research Involving Human Subjects (adopted by the 18th WMA General Assembly, Helsinki, Finland, June 1964, latest amendment 64th WMA General Assembly, Fortaleza, Brazil, October 2013). Methods for the collection of nonmalignant human submandibular salivary gland biopsies, isolation of salivary gland cells, and subsequent passaging have recently been developed (32).

In short, nonmalignant submandibular tissue was collected (after informed consent and Institutional review board approval) from donors undergoing elective tissue head and neck dissection for squamous cell carcinoma of the oral cavity. Biopsies were collected after surgery in HBSS containing 1% BSA. Biopsies were simultaneously mechanically digested using the GentleMACS dissociator (Miltenyi Biotec) and enzymatically digested in HBSS/1% BSA buffer containing collagenase II and hyaluronidase. Cells were plated into 12-well plates at a concentration of 4 × 105 cells per well in DMEM-F12 medium (Gibco) containing 1% penicillin/streptomycin (Gibco), glutamax (2 mmol/L, Gibco), EGF (20 ng/mL, Sigma Aldrich), fibroblast growth factor-2 (20 ng/mL, PeproTech), 1% N2 supplement (Gibco), insulin (10 μg/mL, Sigma), and dexamethasone (1 μmol/L, Sigma). After 3–5 days, primary spheres (passage 0) were disrupted using 0.05% trypsin-EDTA (Gibco) into single cells, counted, and resuspended at 4 × 105 cells/mL. Twenty-five microliters of this solution was mixed with 50 μL Basement Membrane Matrigel (BD Biosciences) and reseeded at the center of 12-well plates. Twenty minutes after seeding, medium supplemented with ROCK1 inhibitor (Y-27632; Sigma-Aldrich) was added to the gels (passage 1) and subsequently incubated at 37°C, 5% CO2. One week after seeding, the media were replaced with dispase (1 mg/mL) and incubated for 30–60 minutes at 37°C to dissolve the Matrigel and release the spheres. Spheres were processed to single cells using 0.05% trypsin-EDTA, sphere and cell numbers were noted, and cells were reseeded in Matrigel.

Human salivary gland organoid irradiations were performed with organoids from 4 different patients at all doses, except 4 Gy (for which N = 3).

Irradiation treatment

Irradiation with X-rays was performed using the Xstrahl 200 X-ray Therapy System at a voltage of 220 kV with a dose rate of 0.52 Gy/minute (26). SOBP C-ion irradiations were performed with (high LET) 90 MeV/u carbon ions in a 3.5-mm long SOBP with a diameter of 30 mm, allowing accurate irradiation of the cells. This resulted in a dose-averaged LET of 150 ± 11 keV/μm at the center of the Matrigel samples, whereas the ± represents the LET at the anterior and posterior of the gel (26). Plateau C-ion irradiations were performed with an LET of 49 ± 2 keV/μm (26).

Irradiation in an intense neutron field with a dose rate of 0.1 Gy/minute were performed at the Physikalisch-Technische Bundesanstalt (Braunschweig, Germany). The neutrons were produced by the 9Be+d reaction on a thick Be− target within a collimator at a deuteron energy of 13 MeV. The energy distribution is broad and extends from about 0.1 to 10 MeV (33). The 0 degree, “free-in-air,” tissue-kerma–averaged mean neutron energy is <En> = 5.8 MeV (34) and the dose-averaged LET in water is 56 keV/μm (35). Absorbed dose to tissue was determined using a calibrated tissue-equivalent ionization chamber according to ICRU Report No. 45 (36). The relative standard uncertainty for the total dose determination was 6%.

The fraction of absorbed dose due to photon radiation was about 2.5%. Specifics regarding cell densities and time of irradiation postseeding are described in the specific assays.

Rat irradiations

Rat irradiations have been previously described in van Luijk and colleagues (13). In brief, 9- to 10-week-old male albino Wistar rats, weighing 230–250 g, were anesthetized and placed in a holder on a positioning rod by their incisors. Using the shoot-through technique, parotid glands were irradiated with 150 MeV protons at the Kernfysisch Versneller Instituut (Groningen, the Netherlands). Forty-five millimeter thick brass collimators were used to give 100%, 50% cranial, or 50% caudal irradiation. Stimulated saliva flow rate was determined 14 days prior to irradiation and up to 180 days postirradiation. Relative flow was calculated over 180 days and normalized to the area under the curve per dose to cranial (stem cell–containing) region. All experiments were carried out in agreement with the Netherlands Experiments on Animals Act (1977) and the European Convention for the Protection of Vertebrates used for Experimental Purpose (Strasbourg, 18.III.1986).

3D clonogenic survival assay

To determine sensitivity of salivary gland and thyroid gland stem/progenitor cells following irradiation, a 3D organoid–based survival assay was performed as described previously (26). Cells were seeded as single cells in Matrigel, as described above, 2 hours preirradiation with 0–4 Gy. Murine salivary gland stem cells were irradiated and replated at a density of 1 × 104 cells per well (0–2 Gy) or 3 × 104 cells per well (4 Gy). For neutron irradiation, cells were irradiated in suspension and subsequently seeded in Matrigel due to the radiation set-up. One-week postirradiation, spheres and cells were counted. Survival was calculated as follows:

For combination treatments with ATMi (2.5 μmol/L; KU-55933; Axon Medchem) or PARPi (0.05 μmol/L; KU-0058948; Axon Medchem), media were changed 24 hours prior to passaging for irradiation and replaced with media supplemented with inhibitor (or DMSO). After 24 hours, organoids were released from Matrigel, processed, and seeded as single cells (as above), and cultured with media supplemented with inhibitor (or DMSO control). Inhibitors were resuspended in DMSO and added to media at 1:1,000.

Wild-type murine salivary gland irradiations were performed at the start of passage 2. All other irradiations (murine thyroid gland, murine Atm/ and WT salivary gland mutant, and human salivary gland) were performed at the start of passage 3.

For fractionation experiments, spheres were counted and surviving fraction were calculated at both 7 days postseeding and 7 days post final irradiation of the fractionation schedule, similar as in previous experiments when seeding and single (and therefore final) irradiation were performed on the same day.

Because primary tissue–derived organoid cultures do not respond as reproducibly as cell lines, inhibition/irradiation experiments were performed using organoids isolated from the tissue of the same mouse with the same genetic background within 1 experiment. This may result in slightly different control curves between sets of experiments.

Immunofluorescent microscopy

For immunofluorescent microscopy, 4 × 104 single cells were seeded in Matrigel 2 hours preirradiation. Dispase was added 30 minutes prior to stated time points. Cells were centrifuged and washed with PBS. Cells were then resuspended in 2% paraformaldehyde, placed on Adhesion slides (Marienfeld), fixed for 15 minutes, and permeabilized with 0.2% Triton for 10 minutes.

All samples were incubated overnight at 4°C with primary antibody [anti–phospho-Histone H2A.X (Ser139); 1:500; Cell Signaling Technology, #2577] followed by incubation at room temperature with secondary antibody (Alexa-Fluor 488; 1:800; Life Technologies, goat anti-rabbit, A11008) for 90 minutes. Nuclear staining was performed using Hoechst 33342 (Molecular Probes, Life Technologies). Imaging was performed using TissueFAxs (Tissuegnostics), foci were analyzed and counted using ImageJ (NIH, Bethesda, MD).

For ADP-ribosylation microscopy, mouse salivary gland cells were seeded on growth factor–reduced Matrigel-coated immunofluorescent microscopy slides following 24 hours treatment with PARPi (or control). After 3 hours, cells were treated with 150 μmol/L H2O2 for 10 minutes, subsequently washed with PBS and fixed by incubating for 10 minutes in 4% PFA. After fixation cells were washed with PBS and then permeabilized by incubating with 0.3% Triton/2% BSA in PBS. Samples were incubated at room temperature for 1 hour with primary antibody (anti-PAR; 1:500; Trevigen, #4336-BPC-100) followed by incubation at room temperature with secondary antibody (Alexa-Fluor 488; 1:1,000; Life Technologies, goat anti-rabbit, A11008) for 60 minutes. Nuclear staining was performed using Hoechst 33342 (Molecular Probes, Life Technologies).

Statistical analysis

All values are represented as mean ± SEM (*, P ≤ 0.05) of at least three independent experiments. Student t test was used for testing statistical significance using GraphPad Prism (GraphPad software).

Glandular stem cells show a nonlinear response at low radiation doses

Previously, it was shown that irradiation of specific regions of the rat parotid gland with a low dose had a disproportionately high impact on the secretory function when compared with higher doses (13). We hypothesize that this is caused by the irradiation of the region that contains the glandular stem cells (8). Therefore, with this recently acquired knowledge, we reanalyzed previously obtained data to determine the in vivo effects of low doses of radiation on normal tissue stem cells. An ablative dose of 30 Gy of protons to the caudal 50% region of the rat salivary gland was used to activate compensatory stem cell regeneration in the cranial 50% region (Fig. 1A; ref. 13). Here, it was shown that subsequent irradiation of the cranial salivary gland region containing a higher density of stem cells (8) with doses from 0 to 10 Gy immediately after the 30 Gy ablative dose to the caudal region, results in a further decrease in pilocarpine-induced salivary flow compared with ablation of the caudal region alone (Fig. 1B). Interestingly, this was not dose dependent, as a dose of 1 Gy led to an unexpectedly strong reduction in salivary gland function, when compared with the 2–10 Gy doses (Supplementary Fig. S1), indicating that the in vivo function of salivary glands may indeed be excessively affected by low radiation doses to the regions with high stem cell densities.

To further determine whether stem cells may be affected disproportionally by low doses, we X-ray–irradiated mouse salivary gland organoid–derived cells, known to contain stem cells, using a recently described methodology (26). These cells have been shown to be able to extensively self-renew and differentiate to all cell types of the salivary gland, both in vitro and in vivo (9, 27, 28). In this culture system, only stem-like cells can expand and form secondary in vitro 3D organoids (Fig. 2A; refs. 9, 27, 28). The capability of a cell to form an organoid upon irradiation can, therefore, be used as a read-out of the radiation response of stem/progenitor cells (26). Indeed, irradiation of organoid-derived cells resulted in general in reduced levels of secondary 3D organoid formation indicating loss of stem/progenitor cells, resulting in a conventional dose-dependent radiation response with doses above 1 Gy (Fig. 2B and C). However, doses below 1 Gy showed a disproportionate reduction in 3D organoid–forming potential with a minimum salivary gland stem cell survival at 0.25 Gy and an increasing salivary gland stem cell survival for doses up to 1 Gy.

To show that this response was not specific for salivary gland stem cells only, we irradiated thyroid gland–derived stem cells, cultured as organoids in a similar manner to salivary gland stem cells. Indeed, the reduction in organoid-forming potential of thyroid gland stem/progenitor cells was very similar, albeit not as pronounced, to that of salivary gland stem cells with a minimum surviving fraction of thyroid gland stem cells also around 0.25 Gy (Fig. 2D). The observed response resembles previously observed low-dose hypersensitivity seen in cell lines (16, 18, 37), albeit with a much lower general radiation sensitivity at higher doses as observed in 3D culture systems of cell lines (26).

DDR after low-dose IR

A critical event determining the radiosensitivity of a cell is its ability to cope with DNA double-strand breaks (DSBs) after irradiation (38). It has been hypothesized that low-dose hypersensitivity originates from lack of DNA repair activation at lower doses (39). To test this, salivary gland organoid–derived cells were X-ray irradiated with 0.25 or 1 Gy. The resulting DSBs were quantified by determining the number of γH2AX foci (40) at 30 minutes and 24 hours postirradiation.

As expected, the initial (30 minutes) number of γH2AX foci per cell increased with increasing doses (Fig. 2E and F). However, similar levels of residual damage were observed at 24 hours following both 0.25 and 1 Gy (Fig. 2E and F). Quantification of cells with 4 or more foci (to reflect an estimate of induced damage based on the majority of cells at 0 Gy having 0–3 foci) showed similar results (Supplementary Fig. S2), indicating that a much greater proportion of DSBs are cleared 24 hours after 1 Gy irradiation than after 0.25 Gy. Taken together these data suggest a less efficient activation of the DDR at 0.25 Gy than 1 Gy.

To determine whether increased sensitivity below 0.5 Gy relates to an impairment in activation of the DDR (18, 39), the complexity of the DNA damage was modulated using high linear energy transfer (LET) particle radiation. Clustered damage induced by such densely ionizing particles is more difficult and slower to repair than damage induced by lower LET X-rays (41). Therefore, salivary gland- or thyroid gland–derived cells were irradiated with carbon ions of either 49 ± 2 keV/μm (Intermediate LET) or 150 ± 11 keV/μm (High LET) or neutrons (mixed LET) with a dose-averaged LET of 56 keV/μm (35). Irradiation with intermediate LET, high LET, and mixed LET indeed resulted in a LET-dependent exponential survival curve for both salivary gland stem cells (Fig. 3A and B; Supplementary Fig. S3A) and thyroid gland stem cells (Fig. 3C and D; Supplementary Fig. S3B) without clear signs of an increased resistance at doses ≥1 Gy in salivary gland stem cells or thyroid gland stem cells. Lower survival of stem cells observed in the ≥1 Gy dose range in response to irradiation with increasing LET may be due to an incapability to repair more complex breaks (41), therefore, it is possible that the DDR is more efficiently activated at doses ≥1 Gy X-rays in contrast to doses <1 Gy in salivary gland and thyroid gland stem cells.

The DDR includes signal transduction pathways that initiate DNA repair, but also cell-cycle checkpoint arrest and cell death (42). The absence of DDR induction has previously been suggested to be the cause of increased sensitivity at lower doses of radiation in 2D-cultured cell lines (43). To further substantiate this in normal tissue stem cells, we next investigated the effect of inhibition of key DDR proteins, such as PARP-1 and ATM. PARP-1 is a key component of single-strand DNA break (SSB) repair, but is also known to be involved in DSB alternative end joining (44). PARP-1 has been shown to play a role in the induction of DDR, possibly due to its delayed activation in response to DNA SSBs induced by ionizing radiation (45). Furthermore, PARP-1 is rapidly recruited to single strand breaks (45), and thus due to the large number of SSBs induced by ionizing radiation could potentially play a role in increased resistance in normal tissue stem cells.

Mouse salivary gland organoid–derived cells were irradiated with X-rays in the presence or absence of a PARP-1 inhibitor (PARPi; KU-0058948) at a concentration which we show to significantly decrease the activity of PARP-1 (Supplementary Fig. S4A and S4B) without being too toxic to the cells (Supplementary Fig. S4C). PARP-1 inhibition resulted in a delayed initiation of DDR visible only at 0.5 Gy (P = 0.0061; Fig. 4A). These data indicate that PARP-1 may play a role, although not essential, in the increased survival of normal tissue stem cells detected at doses ≥0.5 Gy and confirm the involvement of the DDR in salivary gland stem cells.

To verify that the induction of the DDR is essential for induced resistance in normal tissue stem cells, the DDR was impeded by modulating the protein kinase ATM, an apical activator of the DDR to DSBs (46). Failure to activate the ATM-dependent G2–M checkpoint has been suggested to explain low-dose hypersensitivity in established cell lines (37, 43, 46). Therefore, due to its early and pivotal role in many DDR processes, ATM could potentially also contribute to the increased resistance of normal tissue stem cells. Indeed, treatment with ATM inhibitor (ATMi), KU-55933, at a sublethal concentration of the inhibitor (Supplementary Fig. S4D), did not influence the low-dose response of the salivary gland stem cells (Fig. 4B); however, it did abolish the recovery resulting in a lower survival at doses > 0.5 Gy, similar to that observed after particle irradiation. To further confirm the role of ATM in the increased radioresistance in normal tissue stem cells, we irradiated salivary gland stem cells isolated from Atm/ mice. Here also the low-dose response was observed in both wild-type and Atm/ salivary gland stem cells; however, the induction of DDR appears to be completely absent in Atm/ salivary gland stem cells already at 0.5 Gy (Fig. 4C).

To investigate whether this was truly due to an impaired DDR, we next irradiated organoid-derived salivary gland cells in the presence or absence of the ATMi and determined the number of γH2AX foci at 30 minutes and 24 hours postirradiation. No significant differences were observed between ATMi-treated cells and DMSO controls at 30 minutes after both 0.25 and 1 Gy, indicating that the damage induction was not changed (Fig. 4D and E). After 24 hours, as expected, a pronounced reduction in γH2AX foci was observed in DMSO-treated salivary gland cells exposed to 1 Gy. However, no change in γH2AX foci was seen after 0.25 Gy, indicating that DSB repair may not have been initiated at this lower dose. As expected, a significant increase in γH2AX foci was observed at 24 hours after irradiation with 1 Gy in the presence of ATMi. Similar observations were made when the percentage of cells positive for γH2AX (4 or more foci) was scored (Supplementary Fig. S5).

Clinical relevance

To test whether the observed response in normal murine glandular tissue stem cells also occurs in human stem cells obtained from patient salivary gland biopsies (32), we established a stem cell survival dose–response curve. Indeed, low-dose hypersensitivity and increased radioresistance were also observed in human salivary gland stem cells (Fig. 5A and B), indicative of a potential clinical relevance.

For radiotherapy, the dose is mostly delivered at 1.8 to 2 Gy daily fractions to the tumor, resulting in a lower dose per fraction to large parts of the normal surrounding tissue. It could be possible that the low-dose responses are only seen after the first fraction and that the DDR is activated in a second or following fraction, rendering the effect to be of minimal clinical relevance. Therefore, we next irradiated mouse salivary gland organoid–derived cells with 4 daily fractionated doses of 0.25 Gy to mimic clinically relevant doses of radiation to the normal tissue and compared this with a single dose of 1 Gy. Organoids were counted and normalized at day 10 postseeding, being 7 days following the final fraction (Fig. 5C). In contrast to what would be expected in the case of DDR activation, we found that a single dose of 1 Gy consistently resulted in higher survival (92.1% ± 3.0%) than 4 fractions of 0.25 Gy with a total dose of 1 Gy (74.6% ± 0.5%; Fig. 5D and E). This was also true if survival was calculated at 7 days postseeding (Supplementary Fig. S6), as per previous experiments when seeding and single (and therefore final) irradiation were performed on the same day. These data indicate that multiple fractions of low radiation doses, received by normal tissue during fractionated radiotherapy, may be more detrimental to normal salivary gland stem cell survival than a single equivalent higher dose.

Taken together, our data show that both murine and human normal tissue stem cells exhibit a low-dose hypersensitivity and delayed initiation of DDR following higher X-ray doses.

In this study, we demonstrate a major impact of low radiation doses on in vivo salivary gland function, which may be explained by a lack of DDR at low doses resulting in hypersensitivity of normal tissue stem cells. Using an in vitro radiation survival assay of murine salivary gland and thyroid gland organoid–derived stem cells, we show similar effects after 0.25 Gy and 1 Gy of X-rays resulting from a lack or presence of DDR activation, respectively. Increasing the complexity of DSBs using high LET irradiation and inhibition of key DDR proteins did not affect the hypersensitivity at low doses, but resulted in delayed or abolished resistance at higher doses. Importantly, we show the presence of this phenomenon also in patient biopsy–derived submandibular salivary gland organoids. Finally, we demonstrate that multiple clinically relevant fractions of low doses had a far greater impact on survival of salivary gland stem cells than a single fraction of an equivalent dose. Taken together, our findings show the potential clinical relevance of low-dose hypersensitivity to X-rays in adult stem cells that may result in an unexpected functional decline of normal tissues receiving low radiation doses during radiotherapy treatment.

The recent discovery of a region within the salivary gland containing a high density of stem cells that should and can be spared from irradiation (8) further stresses the importance of stem cells in the radiation response of tissues. Indeed, sparing of the high-density stem cell region was suggested to result in an improved normal tissue functioning after irradiation; however, this did not account for a possible low-dose hypersensitivity parameter. Our in vivo rat data show that indeed normal tissue stem cells are hypersensitive to low doses of radiation (1 Gy), albeit this may only be a temporary increased dysfunction compared with 2 or 3 Gy as they displayed the same level of impairment after 180 days. It could be that in the long run a lower number of stem cells could eventually regenerate the tissue. However, the impact on quality of life of earlier salivary gland dysfunction and related xerostomia seems to remain persistent (5).

Using a recently developed technique that allows the culture of tissue-specific organoids from single cells (9, 28, 32), we were able to assess responses of tissue-derived stem cells in vitro. These organoids closely resemble organ tissues, containing tissue-specific stem cells and differentiated cell types, and are more physiologically relevant than traditional 2D monolayer culture systems (23). Moreover, these organoid 3D culture systems allow the expansion of genetically and phenotypically stable adult, organ-specific stem cells (9, 47). Passaging of organoid-derived single cells results in the formation of secondary organoids from stem cells only, as shown by their ability to form all salivary gland lineages and rescue in vivo salivary gland functionality after radiation damage in transplantation studies (9, 28). As such, the radiation response of stem cells may be reflected by the organoid-forming potential of the cells seeded as single cells (26). So far, only 2D cultures of fibroblast and keratinocytes have been studied as representative models of low-dose radiation effects in normal tissue (19). However, normal tissue organoid responses, in combination with the observation that 3D culture systems can give rise to responses more similar to patients than 2D cultures (23), provided us with the unique opportunity to assess stem cell responses in vitro.

The dependency for the occurrence of induced resistance on efficient activation of the DDR is indicated by our DDR inhibition experiments. Our data indicate that PARP-1 may play a nonessential role in the induced resistance, as PARP-1 inhibition simply increased the dose, at which increased resistance was observed but did not remove this effect. The role of PARP-1 in normal tissue stem cells has not been previously investigated at lower doses, but Chalmers and colleagues (45) suggested a role for PARP-1 in increased resistance of tumor cells. If true, this could have consequences for the therapeutic ratio of such a treatment. However, our data suggest that ATM plays an important role in induced resistance in normal salivary gland stem cells because inhibition or knock-out of ATM abrogates this response, through impairment of the DDR. This might be due to a G2-phase checkpoint, as suggested by Krueger and colleagues to occur in 2D-cultured fibroblast and tumor cell lines (39). Rodrigues-Moreira and colleagues demonstrated that low-radiation doses induce persistent oxidative stress, which activates Nrf2 in stem cells resulting in long-term side effects (21). We cannot exclude the possibility of Nrf2 and oxidative stress activation in our system; however, it should be noted that the work of Rodrigues-Moreira and colleagues investigated much lower doses (20 mGy; ref. 21). Cells of the hematopoietic system are known to be more radiosensitive than cells originating from other tissue (48); therefore, it is possible that the radiation response of hematopoietic stem cells differs from solid tissue stem cells. Furthermore, differences in culturing conditions (2D vs. 3D) may affect radiosensitivity (26, 49), and therefore could potentially influence the doses at which increased sensitivity is observed.

Currently, radiation treatment planning considers the consequence of radiation to increase monotonously with dose, whereas stem cell–sparing techniques are being implemented to reduce the dose to normal tissue stem cells (8). These techniques do not take into consideration that normal tissue stem cells may be more sensitive to lower doses. However, based on our data one could conclude that low-dose hypersensitivity may have a significant impact on the normal tissue stem cells and, consequently, on the results of treatment planning. Honore and Bentzen (50) used a modeling study to investigate the role of low-dose hypersensitivity in treatment planning and suggested that much of the benefit of treatment planning may be offset if low-dose hypersensitivity is indeed present. Next to this, it has recently been shown that reducing dose to the high-density stem cell region of the salivary gland during radiotherapy may improve post-therapy salivary gland function following treatment of HNC (8). However, together with the present data, the low-dose response may be an important factor that needs to be considered in treatment planning. Minimizing the dose to normal tissue stem cells may only be optimal if the low-dose hypersensitivity dose range can be avoided in structures within an organ/tissue containing a high stem cell density. Proton therapy could potentially reduce the dose in low-dose regions even further, which may help to avoid the low-dose response of normal tissue stem cells (51).

One proposed mechanism behind late normal tissue complications of radiotherapy is that the cells remaining in the tissue have been compromised with a potential accumulation of sublethal damage. However, it has previously been shown that organ dysfunction and lack of regeneration of the salivary gland is accentuated by ablation of the high density stem cell region (8). We cannot rule out that the effect that we see is due to sublethal damage to the stem cells. Our in vitro clonogenic survival assay is based on the ability to form an organoid, therefore we cannot directly determine whether the stem cells have undergone cell death or have survived but are unable to form an organoid.

Although we find a potentially important hypersensitive low-dose response in normal tissue stem cells, it should be noted that our model has some drawbacks. Though organoids resemble the tissue of origin, the radiation response of normal tissue is more complex than solely the stem cell response. Normal tissue effects, next to the response of stem cells, encompass cytokine, inflammatory, vasculature, and connective tissue responses (1). Interestingly, healthy normal tissue stem cells may modulate the irradiated surrounding by inducing expression of compensatory anti-inflammatory, antifibrotic, and pro-proliferation genes (33). Furthermore, we are aware that in organoid models, stem cells are pushed to proliferate (28) potentially increasing their intrinsic radiosensitivity. Under normal in vivo conditions, stem cells are largely quiescent and remain so until they are stimulated to divide by external factors, such as tissue damage (52). Despite this distinction, our in vivo rat data closely resemble the in vitro observed low-dose hypersensitivity, albeit at a slightly higher dose. Furthermore, while organoids have previously been used to model diseases and for drug screening (23), we show the potential of this system to study in vivo normal tissue responses to radiation, potentially offering a means of patient response prediction in the future.

In conclusion, low-dose hypersensitivity might play a role in the normal tissue response away from the tumor where large volumes are irradiated with low fraction doses upon, for example, IMRT. Reducing the volume and the dose to these regions using particle therapy, such as proton therapy, may more accurately spare designated regions and result in further optimization of radiotherapy for cancer treatment.

J.A. Langendijk reports receiving speakers bureau honoraria and is a consultant/advisory board member for IBA. No potential conflicts of interest were disclosed by the other authors.

Conception and design: P.W. Nagle, N.A. Hosper, L. Barazzuol, J.A. Langendijk, P. van Luijk, R.P. Coppes

Development of methodology: P.W. Nagle, N.A. Hosper, L. Barazzuol, A.L. Jellema, M. Baanstra, S. Brandenburg, J.A. Langendijk

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): P.W. Nagle, N.A. Hosper, L. Barazzuol, A.L. Jellema, M. Baanstra, M.-J. van Goethem, S. Brandenburg, U. Giesen, P. van Luijk

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): P.W. Nagle, N.A. Hosper, L. Barazzuol, P. van Luijk, R.P. Coppes

Writing, review, and/or revision of the manuscript: P.W. Nagle, N.A. Hosper, L. Barazzuol, S. Brandenburg, U. Giesen, J.A. Langendijk, P. van Luijk, R.P. Coppes

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): P.W. Nagle, N.A. Hosper

Study supervision: J.A. Langendijk, R.P. Coppes

We would like to thank Arnout van der Borden, Peter van der Hulst, and Jan Wierenga of the Department of Radiation Oncology, UMCG for their assistance with X-ray irradiation, and the operators and technical staff of the PTB accelerator facility in Braunschweig for their assistance with neutron irradiation. Salivary glands from Atm/ mice were kindly provided by Prof. PA Jeggo of the University of Sussex (Sussex, United Kingdom). The research leading to these results received funding from the European Atomic Energy Community Seventh Framework Programme (FP7/2007-2013) under grant agreement no. 295975 (ANDANTE).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
Bentzen
SM
. 
Preventing or reducing late side effects of radiation therapy: radiobiology meets molecular pathology
.
Nat Rev Cancer
2006
;
6
:
702
13
.
2.
Baumann
M
,
Krause
M
,
Overgaard
J
,
Debus
J
,
Bentzen
SM
,
Daartz
J
, et al
Radiation oncology in the era of precision medicine
.
Nat Rev Cancer
2016
;
16
:
234
49
.
3.
Dijkema
T
,
Terhaard
CH
,
Roesink
JM
,
Braam
PM
,
van Gils
CH
,
Moerland
MA
, et al
Large cohort dose-volume response analysis of parotid gland function after radiotherapy: intensity-modulated versus conventional radiotherapy
.
Int J Radiat Oncol Biol Phys
2008
;
72
:
1101
9
.
4.
Jellema
AP
,
Doornaert
P
,
Slotman
BJ
,
Leemans
CR
,
Langendijk
JA
. 
Does radiation dose to the salivary glands and oral cavity predict patient-rated xerostomia and sticky saliva in head and neck cancer patients treated with curative radiotherapy?
Radiother Oncol
2005
;
77
:
164
71
.
5.
Jellema
AP
,
Slotman
BJ
,
Doornaert
P
,
Leemans
CR
,
Langendijk
JA
. 
Impact of radiation-induced xerostomia on quality of life after primary radiotherapy among patients with head and neck cancer
.
Int J Radiat Oncol Biol Phys
2007
;
69
:
751
60
.
6.
Nutting
CM
,
Morden
JP
,
Harrington
KJ
,
Urbano
TG
,
Bhide
SA
,
Clark
C
, et al
Parotid-sparing intensity modulated versus conventional radiotherapy in head and neck cancer (PARSPORT): a phase 3 multicentre randomised controlled trial
.
Lancet Oncol
2011
;
12
:
127
36
.
7.
Metcalfe
C
,
Kljavin
NM
,
Ybarra
R
,
de Sauvage
FJ
. 
Lgr5+ stem cells are indispensable for radiation-induced intestinal regeneration
.
Cell Stem Cell
2014
;
14
:
149
59
.
8.
van Luijk
P
,
Pringle
S
,
Deasy
JO
,
Moiseenko
VV
,
Faber
H
,
Hovan
A
, et al
Sparing the region of the salivary gland containing stem cells preserves saliva production after radiotherapy for head and neck cancer
.
Sci Transl Med
2015
;
7
:
305ra147
.
9.
Maimets
M
,
Rocchi
C
,
Bron
R
,
Pringle
S
,
Kuipers
J
,
Giepmans
BN
, et al
Long-term in vitro expansion of salivary gland stem cells driven by wnt signals
.
Stem Cell Reports
2016
;
6
:
150
62
.
10.
Bijl
HP
,
van Luijk
P
,
Coppes
RP
,
Schippers
JM
,
Konings
AW
,
van der Kogel
AJ
. 
Unexpected changes of rat cervical spinal cord tolerance caused by inhomogeneous dose distributions
.
Int J Radiat Oncol Biol Phys
2003
;
57
:
274
81
.
11.
Bijl
HP
,
van Luijk
P
,
Coppes
RP
,
Schippers
JM
,
Konings
AW
,
van Der Kogel
AJ
. 
Regional differences in radiosensitivity across the rat cervical spinal cord
.
Int J Radiat Oncol Biol Phys
2005
;
61
:
543
51
.
12.
Bijl
HP
,
van Luijk
P
,
Coppes
RP
,
Schippers
JM
,
Konings
AW
,
van der Kogel
AJ
. 
Influence of adjacent low-dose fields on tolerance to high doses of protons in rat cervical spinal cord
.
Int J Radiat Oncol Biol Phys
2006
;
64
:
1204
10
.
13.
van Luijk
P
,
Faber
H
,
Schippers
JM
,
Brandenburg
S
,
Langendijk
JA
,
Meertens
H
, et al
Bath and shower effects in the rat parotid gland explain increased relative risk of parotid gland dysfunction after intensity-modulated radiotherapy
.
Int J Radiat Oncol Biol Phys
2009
;
74
:
1002
5
.
14.
Franken
NA
,
Rodermond
HM
,
Stap
J
,
Haveman
J
,
van Bree
C
. 
Clonogenic assay of cells in vitro
.
Nat Protoc
2006
;
1
:
2315
9
.
15.
Wang
JZ
,
Huang
Z
,
Lo
SS
,
Yuh
WT
,
Mayr
NA
. 
A generalized linear-quadratic model for radiosurgery, stereotactic body radiation therapy, and high-dose rate brachytherapy
.
Sci Transl Med
2010
;
2
:
39ra48
.
16.
Marples
B
,
Joiner
MC
. 
The response of chinese hamster V79 cells to low radiation doses: evidence of enhanced sensitivity of the whole cell population
.
Radiat Res
1993
;
133
:
41
51
.
17.
Martin
LM
,
Marples
B
,
Lynch
TH
,
Hollywood
D
,
Marignol
L
. 
Exposure to low dose ionising radiation: molecular and clinical consequences
.
Cancer Lett
2014
;
349
:
98
106
.
18.
Wouters
BG
,
Sy
AM
,
Skarsgard
LD
. 
Low-dose hypersensitivity and increased radioresistance in a panel of human tumor cell lines with different radiosensitivity
.
Radiat Res
1996
;
146
:
399
413
.
19.
Slonina
D
,
Biesaga
B
,
Urbanski
K
,
Kojs
Z
. 
Low-dose radiation response of primary keratinocytes and fibroblasts from patients with cervix cancer
.
Radiat Res
2007
;
167
:
251
9
.
20.
Blanpain
C
,
Fuchs
E
. 
Epidermal homeostasis: a balancing act of stem cells in the skin
.
Nat Rev Mol Cell Biol
2009
;
10
:
207
17
.
21.
Rodrigues-Moreira
S
,
Moreno
SG
,
Ghinatti
G
,
Lewandowski
D
,
Hoffschir
F
,
Ferri
F
, et al
Low-dose irradiation promotes persistent oxidative stress and decreases self-renewal in hematopoietic stem cells
.
Cell Rep
2017
;
20
:
3199
211
.
22.
Kretzschmar
K
,
Clevers
H
. 
Organoids: modeling development and the stem cell niche in a dish
.
Dev Cell
2016
;
38
:
590
600
.
23.
Clevers
H
. 
Modeling development and disease with organoids
.
Cell
2016
;
165
:
1586
97
.
24.
Yui
S
,
Nakamura
T
,
Sato
T
,
Nemoto
Y
,
Mizutani
T
,
Zheng
X
, et al
Functional engraftment of colon epithelium expanded in vitro from a single adult Lgr5(+) stem cell
.
Nat Med
2012
;
18
:
618
23
.
25.
Huch
M
,
Dorrell
C
,
Boj
SF
,
van Es
JH
,
Li
VS
,
van de Wetering
M
, et al
In vitro expansion of single Lgr5+ liver stem cells induced by wnt-driven regeneration
.
Nature
2013
;
494
:
247
50
.
26.
Nagle
PW
,
Hosper
NA
,
Ploeg
EM
,
van Goethem
MJ
,
Brandenburg
S
,
Langendijk
JA
, et al
The in vitro response of tissue stem cells to irradiation with different linear energy transfers
.
Int J Radiat Oncol Biol Phys
2016
;
95
:
103
11
.
27.
Lombaert
IM
,
Brunsting
JF
,
Wierenga
PK
,
Faber
H
,
Stokman
MA
,
Kok
T
, et al
Rescue of salivary gland function after stem cell transplantation in irradiated glands
.
PLoS One
2008
;
3
:
e2063
.
28.
Nanduri
LS
,
Baanstra
M
,
Faber
H
,
Rocchi
C
,
Zwart
E
,
de Haan
G
, et al
Purification and ex vivo expansion of fully functional salivary gland stem cells
.
Stem Cell Reports
2014
;
3
:
957
64
.
29.
Barlow
C
,
Hirotsune
S
,
Paylor
R
,
Liyanage
M
,
Eckhaus
M
,
Collins
F
, et al
Atm-deficient mice: a paradigm of ataxia telangiectasia
.
Cell
1996
;
86
:
159
71
.
30.
Liao
MJ
,
Yin
C
,
Barlow
C
,
Wynshaw-Boris
A
,
van Dyke
T
. 
Atm is dispensable for p53 apoptosis and tumor suppression triggered by cell cycle dysfunction
.
Mol Cell Biol
1999
;
19
:
3095
102
.
31.
Pringle
S
,
Nanduri
LS
,
van der Zwaag
M
,
van Os
R
,
Coppes
RP
. 
Isolation of mouse salivary gland stem cells
.
J Vis Exp
2011
;
48
:
2484
.
32.
Pringle
S
,
Maimets
M
,
van der Zwaag
M
,
Stokman
MA
,
van Gosliga
D
,
Zwart
E
, et al
Human salivary gland stem cells functionally restore radiation damaged salivary glands
.
Stem Cells
2016
;
34
:
640
52
.
33.
Brede
HJ
,
Dietze
G
,
Kudo
K
,
Schrewe
UJ
,
Tancu
F
,
Wen
C
. 
Neutron yields from thick be targets bombarded with deuterons or protons
.
Nucl Instrum Methods Phys Res A
1989
;
274
:
332
44
.
34.
Malmer
CJ
. 
ICRU report 63. nuclear data for neutron and proton radiotherapy and for radiation protection
.
Med Phys
2001
;
28
:
861
-.
35.
Schrewe
UJ
,
Brede
HJ
,
Dietze
G
. 
Investigation of tissue-equivalent proportional counters in mixed neutron photon fields also applying time-of-flight techniques
.
Radiat Prot Dosimetry
1988
;
23
:
239
.
36.
International Commission on Radiation Units and Measurements
. 
ICRU report no. 45; clinical neutron dosimetry, part 1: determination of absorbed dose in a patient treated by external beams of fast neutrons
; 
1989
.
37.
Krueger
SA
,
Joiner
MC
,
Weinfeld
M
,
Piasentin
E
,
Marples
B
. 
Role of apoptosis in low-dose hyper-radiosensitivity
.
Radiat Res
2007
;
167
:
260
7
.
38.
Vignard
J
,
Mirey
G
,
Salles
B
. 
Ionizing-radiation induced DNA double-strand breaks: a direct and indirect lighting up
.
Radiother Oncol
2013
;
108
:
362
9
.
39.
Krueger
SA
,
Collis
SJ
,
Joiner
MC
,
Wilson
GD
,
Marples
B
. 
Transition in survival from low-dose hyper-radiosensitivity to increased radioresistance is independent of activation of ATM Ser1981 activity
.
Int J Radiat Oncol Biol Phys
2007
;
69
:
1262
71
.
40.
Rogakou
EP
,
Pilch
DR
,
Orr
AH
,
Ivanova
VS
,
Bonner
WM
. 
DNA double-stranded breaks induce histone H2AX phosphorylation on serine 139
.
J Biol Chem
1998
;
273
:
5858
68
.
41.
Eccles
LJ
,
Lomax
ME
,
O'Neill
P
. 
Hierarchy of lesion processing governs the repair, double-strand break formation and mutability of three-lesion clustered DNA damage
.
Nucleic Acids Res
2010
;
38
:
1123
34
.
42.
Short
SC
,
Woodcock
M
,
Marples
B
,
Joiner
MC
. 
Effects of cell cycle phase on low-dose hyper-radiosensitivity
.
Int J Radiat Biol
2003
;
79
:
99
105
.
43.
Krueger
SA
,
Wilson
GD
,
Piasentin
E
,
Joiner
MC
,
Marples
B
. 
The effects of G2-phase enrichment and checkpoint abrogation on low-dose hyper-radiosensitivity
.
Int J Radiat Oncol Biol Phys
2010
;
77
:
1509
17
.
44.
Mansour
WY
,
Rhein
T
,
Dahm-Daphi
J
. 
The alternative end-joining pathway for repair of DNA double-strand breaks requires PARP1 but is not dependent upon microhomologies
.
Nucleic Acids Res
2010
;
38
:
6065
77
.
45.
Chalmers
A
,
Johnston
P
,
Woodcock
M
,
Joiner
M
,
Marples
B
. 
PARP-1, PARP-2, and the cellular response to low doses of ionizing radiation
.
Int J Radiat Oncol Biol Phys
2004
;
58
:
410
9
.
46.
Lobrich
M
,
Jeggo
PA
. 
The impact of a negligent G2/M checkpoint on genomic instability and cancer induction
.
Nat Rev Cancer
2007
;
7
:
861
9
.
47.
Sato
T
,
Stange
DE
,
Ferrante
M
,
Vries
RG
,
Van Es
JH
,
Van den Brink
S
, et al
Long-term expansion of epithelial organoids from human colon, adenoma, adenocarcinoma, and barrett's epithelium
.
Gastroenterology
2011
;
141
:
1762
72
.
48.
Shao
L
,
Luo
Y
,
Zhou
D
. 
Hematopoietic stem cell injury induced by ionizing radiation
.
Antioxid Redox Signal
2014
;
20
:
1447
62
.
49.
Eke
I
,
Cordes
N
. 
Radiobiology goes 3D: How ECM and cell morphology impact on cell survival after irradiation
.
Radiother Oncol
2011
;
99
:
271
8
.
50.
Honore
HB
,
Bentzen
SM
. 
A modelling study of the potential influence of low dose hypersensitivity on radiation treatment planning
.
Radiother Oncol
2006
;
79
:
115
21
.
51.
Vissink
A
,
van Luijk
P
,
Langendijk
JA
,
Coppes
RP
. 
Current ideas to reduce or salvage radiation damage to salivary glands
.
Oral Dis
2015
;
21
:
e1
10
.
52.
Li
L
,
Clevers
H
. 
Coexistence of quiescent and active adult stem cells in mammals
.
Science
2010
;
327
:
542
5
.