The salivary glands are often damaged during head and neck cancer radiotherapy. This results in chronic dry mouth, which adversely affects quality of life and for which there is no long-term cure. Mouse models of salivary gland injury are routinely used in regenerative research. However, there is no clear consensus on the radiation regime required to cause injury. Here, we analysed three regimes of γ-irradiation of the submandibular salivary gland. Transcriptional analysis, immunofluorescence and flow cytometry was used to profile DNA damage, gland architecture and immune cell changes 3 days after single doses of 10 or 15 Gy or three doses of 5 Gy. Irrespective of the regime, radiation induced comparable levels of DNA damage, cell cycle arrest, loss of glandular architecture, increased pro-inflammatory cytokines and a reduction in tissue-resident macrophages, relative to those observed in non-irradiated submandibular glands. Given these data, coupled with the fact that repeated anaesthetic can negatively affect animal welfare and interfere with saliva secretion, we conclude that a single dose of 10 Gy irradiation is the most refined method of inducing acute salivary gland injury in a mouse model.

Therapeutic radiation remains a predominant treatment for targeting various cancers. Although recent advances in the delivery of radiation, such as intensity-modulated radiation therapy, aim to minimise off-target side effects, healthy tissues that lie within the therapeutic field often also receive high doses of radiation, leading to cellular damage and organ dysfunction (Bentzen, 2006). The salivary glands are often irreversibly damaged following radiotherapy for head and neck cancer, and the resulting salivary gland dysfunction results in xerostomia, or chronic dry mouth (Lechelt et al., 2018). These adverse side effects of radiation therapy substantially decrease the quality of life in survivors and there is currently no long-term regenerative therapy for this debilitating condition (Bentzen, 2006). Therefore, to better the quality of life in patients, the effects of radiation injury on salivary glands and the mechanisms by which the tissue repairs itself or fails to repair itself must be better understood.

Mice provide excellent preclinical models owing to their anatomical similarities to humans (Emmerson and Knox, 2018), coupled with their genetic tractability and the availability of transgenic strains. Research into salivary gland radiation injury routinely uses mice to model the tissue injury that patients undergoing radiotherapy experience. However, there is no clear consensus on the dose that is required to mirror the injury evident in patients. Patients with head and neck cancer generally receive between 50 and 70 Gy radiation (Bourhis et al., 2012; Thomson et al., 2020) targeted to the tumour, which is split over a 6 week period (termed fractionated radiotherapy). Although radiotherapy is often not purposefully targeted at the salivary glands, as the proportion of head and neck cancer that arises from salivary gland tissue is very low (<5%) (Lawal et al., 2015; Mezi et al., 2020), the salivary glands often inadvertently fall into the irradiation field and can experience doses of anywhere between 8 and 55 Gy (Hawkins et al., 2018; van Timmeren et al., 2021). This results in cell disorganisation, collagen deposition and adipocyte accumulation (Nam et al., 2016), as well as a loss of saliva production (Eisbruch et al., 2001; Henson et al., 1999) and/or a change in saliva consistency (Winter et al., 2021; Funegård et al., 1994; Makkonen et al., 1986). Acute hyposalivation occurs as early as 1 week after radiotherapy (Franzén et al., 1992), whereas longer-term changes are associated with acinar cell loss, ductal metaplasia and fibrosis (Sullivan et al., 2005) and lead to xerostomia in the majority of patients (Wijers et al., 2002). However, as human biopsy samples can only be accessed when a donor is already undergoing a surgical intervention, there is a dearth of data on the acute effects of radiation injury on human salivary gland tissue. As such, animal models are crucial to understand these early effects at the cellular level.

To date, various doses and regimes of radiation have been used in mouse models to mirror the human clinical situation, including single doses ranging from 1 Gy (Limesand et al., 2006; Humphries et al., 2006; Avila et al., 2009) to 15 Gy (Peng et al., 2020; Lombaert et al., 2008; Bralic et al., 2005; Takeda et al., 2003; Cotrim et al., 2007, 2005; Urek et al., 2005; Limesand et al., 2010), and fractionated doses ranging from 4 Gy (Limesand et al., 2010) to 30 Gy (Cotrim et al., 2007). Table 1 outlines the various radiation regimes used in the field and includes the type of radiation used in each study. However, it is widely known that different mouse strains or genetic backgrounds exhibit differences in their sensitivity to radiation in other organs (Jackson et al., 2010; Down and Yanch, 2010; Down et al., 1986; Dileto and Travis, 1996), and, as such, we have also included details of the genetic background used in all published studies. Regardless of the dose and type of irradiation injury, there are some commonalities reported across studies, including DNA damage, changes in acinar cells, fibrosis and loss of saliva production.

In this study, we compared three regimes of γ-irradiation, namely, single doses of 10 Gy and 15 Gy, and three doses of 5 Gy delivered over a period of 5 days (Fig. 1C), in order to ascertain and compare the levels of cellular damage elicited, focusing on the submandibular salivary gland (SMG), the best studied of the three major salivary glands in mice. We examined SMGs at 3 days after irradiation, as we have previously shown that, at this timepoint, we see significant cellular damage, disruption of the epithelial architecture and transcriptional changes associated with DNA damage and stress response to injury (Emmerson et al., 2018; McKendrick et al., 2023). Using transcriptional analysis, immunofluorescence staining and flow cytometry, we characterised epithelial and immune cell changes following the three different regimes of radiation injury, and we show that although there are subtle differences among different approaches, all regimes elicit comparable damage relative to that in unirradiated SMGs.

We used a well-established model of targeted radiation to compare three regimes of γ-irradiation (Fig. 1A). This model exposes the three pairs of salivary glands, while protecting the nose and the body with lead (Fig. 1B). We delivered a single dose of either 10 or 15 Gy γ-irradiation, or three doses of 5 Gy irradiation over a period of 5 days (hereafter 3×5 Gy), followed by analysis 3 days after the final dose (Fig. 1C). We and others have previously demonstrated that radiation to either the SMG, sublingual gland or parotid gland causes transcriptional changes in the hours and days immediately after irradiation (4 and 8 h and 1-7 days post irradiation) that are indicative of DNA and cellular damage, including elevation in the expression of genes involved in apoptosis and cell cycle arrest [Bax and Cdkn1a (encoding p21)] (Avila et al., 2009, Emmerson et al., 2018, McKendrick et al., 2023) and a reduction in the expression of genes involved in cell cycle (Mki67 and Ccnd1) and salivary gland function [Aqp5 and Mist1 (also known as Bhlha15)] (Emmerson et al., 2018; McKendrick et al., 2023). In this study, we first analysed gene expression changes at 3 days post irradiation in unirradiated SMGs (0 Gy) and SMGs exposed to 10 Gy, 15 Gy or 3×5 Gy. We found that all three radiation regimes elicited a similar and significant increase in the expression of Bax, a gene involved in the apoptotic and cell death pathway (Westphal et al., 2014), and a reduction in the expression of cyclin D1 (Ccnd1) and the water channel-encoding gene aquaporin 5 (Aqp5) (Fig. 1D). However, the expression of Cdkn1a was unchanged between unirradiated and irradiated SMGs (Fig. 1D).

We subsequently performed immunofluorescence analysis for markers of tissue architecture, DNA damage and cell cycle. We first analysed the epithelial architecture by staining for the epithelial marker E-cadherin (ECAD, encoded by Cdh1) and the acinar water channel aquaporin 5 (AQP5). AQP5 is necessary for water transfer and plays a major role in the salivary secretion process (Delporte and Steinfeld, 2006); thus, a loss of expression or aberrant localisation leads to salivary gland dysfunction. We found that, compared to the unirradiated control, the irradiated samples exhibited disruption to the tissue architecture, as previously reported (McKendrick et al., 2023), regardless of the radiation regime used (Fig. 2A,A′). We quantified AQP5 and ECAD fluorescence intensity (Fig. 2C,D) and found AQP5 levels to be in line with gene expression data (Fig. 1C). AQP5 fluorescence intensity was significantly reduced in all three radiation regimes (Fig. 2C), in line with the observed acinar cell disruption (Fig. 2A), whereas ECAD intensity was unchanged (Fig. 2A,D). This suggests that acinar cells are more sensitive than ductal epithelial cells to radiation at this early timepoint. Furthermore, we analysed acinar cell atrophy by measuring average acinar cell diameter as previously described (Lombaert et al., 2020). We found that acinar cell diameter was significantly reduced following irradiation, but it was not different among irradiation regimes (Fig. 2A,E). We also performed immunofluorescence analysis for CD31 (PECAM1), which is found on endothelial cells and is used as a marker of vasculature (Lertkiatmongkol et al., 2016). In line with previously published studies (Cotrim et al., 2007; McKendrick et al., 2023), we found that the fluorescence intensity of CD31 was modestly but significantly reduced following irradiation (Fig. 2B,F), suggesting that irradiation directly impacts vessel arrangement and architecture.

We next examined cellular damage in response to radiation. We first performed immunofluorescence analysis for 53BP1 (encoded by Trp53bp1) to analyse the extent of DNA damage, as it forms visible foci in nuclei following exposure to ionising radiation (Gupta et al., 2014; Shibata and Jeggo, 2020). Indeed, we have previously shown that 53BP1+ foci are evident in the nuclei of ductal cells of the murine SMG at 1 and 3 days post irradiation (McKendrick et al., 2023). Here, we found evidence of 53BP1+ foci in the nuclei of cells in all three radiation regimes (Fig. 3A,A′), and the numbers of foci were significantly elevated compared to those in non-irradiated controls, in which there were no foci evident (Fig. 3D), consistent with previous findings for 10 and 15 Gy irradiation (Marmary et al., 2016; McKendrick et al., 2023). Although we found that the expression of Cdkn1a (encoding p21) was unchanged following radiation (Fig. 1C), we found that the numbers of p21+ cells in the irradiated tissues were significantly higher than those in non-irradiated controls (Fig. 3B,E). In contrast, we found that the numbers of proliferating cells, marked by Ki67 (encoded by Mki67), were significantly reduced following radiation, irrespective of the regime used (Fig. 3C,F). Overall, these findings demonstrate cellular damage caused by irradiation, leading to cell cycle changes.

Radiation injury leads to inflammation and immune cell influx in numerous tissues (Cheng et al., 2021; Li et al., 2020; Nolan et al., 2022; Meziani et al., 2018; Groves et al., 2023). Conversely, tissue-resident immune cell populations and their niche can be altered following radiation injury (Groves et al., 2023; Hudson et al., 2023; Soysa et al., 2019; Venkateswaran et al., 2019). Although these changes can be beneficial if at the site of the tumour (Shinde-Jadhav et al., 2021), they can lead to negative effects on healthy surrounding tissues that are inadvertently irradiated during therapy. We found that radiation induced expression of the genes encoding the cytokines IL-6 (Il6) and TGFβ (Tgfb1) in whole SMGs (Fig. 4A). Sustained expression of IL-6 has been associated with salivary gland senescence and long-term hypofunction (Marmary et al., 2016), so the early upregulation seen in this study is not surprising and sustained expression could be causative in the longer-term dysfunction described in this previously published study. Furthermore, TGFβ plays a pivotal role in fibrosis after radiation injury (Dai et al., 2019) and accelerates the DNA damage response in epithelial cells in vitro via SMAD signalling (Lee et al., 2016). Interestingly, although significant upregulation of Il6 and Tgfb1 was observed in response to a single dose of 15 Gy radiation, this was not the case in the fractionated group, suggesting that low-dose-sustained injury does not have such a significant effect on cytokine expression, at least at this early timepoint after injury (Fig. 4A).

Finally, we have recently demonstrated the crucial role that macrophages play in salivary gland regeneration after irradiation injury (McKendrick et al., 2023). In this study, we examined the number of CD45 (encoded by Ptprc)-positive cells (i.e. total leukocytes) and macrophages at 3 days after radiation. We have previously shown that macrophages are the most abundant leukocyte in the salivary gland and that they express high levels of F4/80 (encoded by Adgre1) and medium to high levels of CD11b (encoded by Itgam) (McKendrick et al., 2023). Thus, we used these two markers to analyse macrophages in the current study, using both immunofluorescence staining coupled with confocal microscopy and flow cytometry. Although the number of total CD45+ leukocytes, measured by immunofluorescence analysis, was unchanged between unirradiated and radiated SMGs (Fig. 4B,D), the number of F4/80+ macrophages was significantly reduced after radiation in all regimes examined (Fig. 4C,E). This is in line with our previous findings that total SMG macrophages are significantly reduced at 3 days after 10 Gy targeted γ-irradiation, but that this loss is recovered by 7 days (McKendrick et al., 2023). This is likely due to the observed loss of a proliferative subpopulation of macrophages, which is eventually recovered (McKendrick et al., 2023). In order to explore this further, we confirmed this result by flow cytometry (Fig. 4F-I). The gating strategy for flow cytometry is shown in Fig. 4F. Here, we found that the number of CD45+ cells (as a percentage of total live cells) was significantly reduced with 3×5 Gy irradiation, but not with 10 or 15 Gy irradiation (Fig. 4G). However, we found that the total number of CD11b+ F4/80+ macrophages was significantly reduced following irradiation in all regimes (Fig. 4H,I), in line with results from immunofluorescence analyses.

There is currently no established radiation regime to model radiotherapy-induced salivary gland injury in mouse models, and approaches range from low, single doses to large, fractionated doses. Here, we compared three different radiation regimes in order to directly compare the early injury that is induced, and show that although there are some subtle differences in different regimes, all approaches elicit similar levels of epithelial disruption, DNA damage and changes in cell cycle, and alterations to tissue-resident macrophages.

It is well established that off-target radiation causes a decline in saliva production (Eisbruch et al., 2001; Henson et al., 1999) changes in saliva consistency (Winter et al., 2021; Funegård et al., 1994; Makkonen et al., 1986), acinar cell loss, ductal metaplasia and fibrosis (Sullivan et al., 2005) in patients with head and neck cancer undergoing radiotherapy. However, it is challenging to explore the very early effects of radiation injury in human biopsies as material is scarce and the confounding factors of human studies make real data difficult to interpret. Thus, animal models provide a valuable resource to understand the cellular effects.

In this study, we observed transcriptional changes indicative of DNA damage in all three radiation regimes. We found elevated Bax expression at 3 days post irradiation, in line with previous studies that describe elevation of Bax expression at 3, 24 and 48 h post irradiation (Uchida et al., 2022), and elevation of Bax and Puma (also known as Bbc3) expression at 4 and 8 h post irradiation (Avila et al., 2009). However, these data are at odds with previous histological analysis of irradiated rat and mouse SMGs at 48 h post irradiation, which did not reveal detectable apoptosis (Marmary et al., 2016; Choi et al., 2009). This suggests that although the tissue exhibits DNA damage, either it does not undergo cell death via apoptosis, or the window of apoptosis is missed in these studies. Although apoptosis has been well described in the developing embryonic salivary gland (Teshima et al., 2016), reports on protein levels in the adult salivary gland after irradiation are scarce, with some exceptions (Ren et al., 2022; Gilman et al., 2019). However, we find that targeted radiation to the murine SMG leads to DNA damage foci, in line with previous reports in mice (Marmary et al., 2016; Ren et al., 2022; Uchida et al., 2022; Chang et al., 2019; Varghese et al., 2018). Given that radiation also causes DNA damage (marked by γH2AX and/or 53BP1 foci) in salivary gland cells in vitro (Uchida et al., 2022; Chang et al., 2019; Affandi et al., 2021) and salivary gland tissues ex vivo (Song et al., 2021), we speculate that this is a direct effect on the cells. Indeed, lymphocytes exhibit 53BP1 and γH2AX foci with increasing doses of irradiation (Sorokina et al., 2013). It is interesting that DNA damage is readily observed but apoptosis is not evident; this may mean that the DNA damage is not in fact leading to cell death and apoptosis. This phenomenon has been demonstrated in mammary gland and hair follicles, where stem cells exhibit DNA damage but are resistant to apoptosis, and instead exhibit enhanced DNA repair activity and continue to proliferate (Chang et al., 2015; Insinga et al., 2013; Sotiropoulou et al., 2010).

It has been proposed that the mechanism of action for the observed acinar cell loss following radiation is due to sublethal DNA damage, which manifests over time and finally becomes lethal at a delayed phase, resulting in acinar cell loss upon proliferation and replacement (Nagler, 2002). Given that we saw evidence of DNA damage in all three radiation regimes, coupled with the evident loss of AQP5 intensity, as previously described (Uchida et al., 2022; Choi et al., 2009), and reduction in acinar cell diameter, we propose that all three models tested here are capable of recapitulating the observed acinar cell loss in human salivary glands (Nam et al., 2016). Furthermore, we saw a reduction in the expression of Aqp5 after 3 days, similar to what has been previously reported at 3 h post irradiation (Uchida et al., 2022). However, the previous study found that expression recovered by 24 h, whereas we observed that expression remained reduced at 3 days post irradiation across all three radiation regimes. Furthermore, we found that the expression of E-cadherin was not changed following irradiation, as has been previously described (Uchida et al., 2022). Although expression of Cdkn1a, the gene encoding p21, was unchanged at 3 days post irradiation, we observed significantly more p21+ cells in the SMG following irradiation, irrespective of the radiation regime analysed.

To date, most studies in which inflammatory cell infiltrate or resident inflammatory cells have been reported in post-radiation salivary glands have relied on Hematoxylin- and Eosin-stained images (Nam et al., 2016; Muhvic-Urek et al., 2006; Kim et al., 2019). However, in order to gain an accurate representation of inflammation, immunofluorescence and/or flow cytometric analysis is necessary. We found that although total leukocyte numbers were largely unchanged following irradiation injury, there was a reduction in the number of tissue-resident (F4/80+ CD11b+) macrophages, in line with our previous findings (McKendrick et al., 2023). Furthermore, the expression of pro-inflammatory cytokines was upregulated, in line with previous reports (Uchida et al., 2022). However, we cannot rule out that this may be an indirect effect: for example, Cxcl2, encoding a pro-inflammatory cytokine produced by macrophages, is upregulated in both oral tissues and SMGs after irradiation (Uchida et al., 2022; Shen et al., 2018). Moreover, bystander effects have been reported in non-irradiated contralateral salivary glands (Uchida et al., 2022), suggesting that irradiation injury also has systemic, long-range effects on salivary gland. This could explain why patients lose significant salivary gland function even when radiation is targeted on one side. We have previously shown that exposure to a single dose of 10 Gy γ-irradiation led to a substantial reduction in a proliferating subset of tissue-resident macrophages after 3 days (McKendrick et al., 2023). In all three radiation regimes tested in this study, we found a significant reduction in total macrophage numbers, which implies that exposure to either a single dose or multiple fractionated doses has a profound effect on tissue-resident macrophages.

In summary, we have demonstrated that, although there are subtle differences, the three radiation regimes analysed here are comparable in terms of DNA damage, cellular changes and macrophage abundance, and all three regimes mirror the early effects of radiotherapy injury. Furthermore, in order to induce injury in mouse salivary glands, most researchers will anaesthetise mice to immobilise them and directly irradiate the neck area only. However, a recent study has reported that repeated administration of ketamine affects salivary gland secretion and interferes with the reported loss following irradiation injury (Uchida et al., 2022). Thus, a model that involves minimum anaesthesia but recapitulates the injury evident in humans is both the most refined method and likely the most scientifically robust. However, we recognise that there will be differences in the source of ionising radiation and, as such, studies using X-rays may differ from those using γ-irradiation. That said, recent evidence suggests that γ-irradiation and X-ray irradiation elicit comparable effects in a model of bone marrow depletion for reconstitution (Wittenborn et al., 2021). Comparative research is required to ascertain whether this is also the case in salivary gland radiation injury. However, overall and based on the data presented, we propose that the lowest single dose tested here (10 Gy) is scientifically sound and the most appropriate model for future studies into acute cellular salivary gland damage.

Limitations of the study

Although this study is informative, it has some caveats that should be recognised. Firstly, we only analysed SMGs from female mice; however, it is becomingly increasingly evident that sexual dimorphism exists in murine salivary gland structure and transcriptome (Mukaibo et al., 2019; Pinkstaff, 1998; Chai et al., 1993; Brown et al., 2020). Indeed, we have recently shown that there are subtle differences in immune cell populations between male and female mouse SMGs (McKendrick et al., 2023). Thus, male mice may exhibit subtle differences to female mice in response to radiation. This is not uncommon: male and female mice respond differently to radiation-induced arthrofibrosis (Rodman et al., 2022), lung toxicity (Bilal et al., 2014) and memory impairment (Perez et al., 2018); and male and female patients differ in sensitivity to irradiation-induced heart toxicity (Snyder et al., 2012) and laryngeal/pharyngeal side effects (Page et al., 2018). However, numerous studies, including our own, have used both male and female mice in salivary gland regeneration and have seen equivalent results (May et al., 2018; Emmerson et al., 2018; McKendrick et al., 2023); thus, we predict that male mice would also elicit comparable responses to the different regimes of irradiation tested here.

Secondly, we only assessed one timepoint following radiation in this study. However, we have previously shown that samples at 3 days post irradiation exhibit hallmarks of injury, including the expression of genes associated with apoptosis, loss of functional proteins, initiation of a DNA damage response and changes in macrophage proliferation (Emmerson et al., 2018; McKendrick et al., 2023); thus, this is a timepoint of significant interest. However, we cannot rule out the possibility that early and late effects are not differentially regulated across the three regimes tested. For example, following thoracic radiotherapy, fibrosis is not typically observed until 4-6 months post irradiation (Grundmann et al., 2009; Liu et al., 2018; De la Cal et al., 2012), and the extent of lung fibrosis increases as the dose of radiation increases (Rosen et al., 2001). Importantly, TGFβ signalling is a crucial mediator of the fibrotic response (reviewed by Meng et al. 2016; Frangogiannis 2020). However, it is becoming evident that expression of TGFβ in tissues including the lung and intestine varies depending on the dose and duration of radiation (Meyer et al., 2017b). Although the role of TGFβ in radiation-induced salivary gland fibrosis is yet to be elucidated, TGFβ expression is elevated in human salivary glands at 3 to 156 months after radiotherapy, and in mouse salivary glands at 2 weeks after radiation injury (Anscher et al., 1995; Hakim et al., 2011; Spiegelberg et al., 2014). However, a recent study has demonstrated elevated expression of TGFβ following ductal ligation (Woods et al., 2015), another salivary gland injury model that aims to recapitulate the damage caused by a ductal blockage (i.e. chronic obstructive sialadenitis or salivary gland stones) and exhibits significant fibrosis. Elevated TGFβ is also evident in clinical samples of human chronic obstructive sialadenitis (Teymoortash et al., 2003; Kizu et al., 1996). Given these data, it seems likely that elevated salivary gland TGFβ expression is a hallmark and possible cause of tissue fibrosis, but it is not unique to radiation injury. Cellular senescence is another key component of radiation-induced salivary gland dysfunction, which is mediated, at least in part, by IL-6 signalling (Marmary et al., 2016), but does not become apparent until 8-12 weeks post radiation. Although we report that the acute injury phase of irradiation-induced salivary gland damage is comparable between single doses of 10 Gy and 15 Gy and a fractionated dose of 3×5 Gy and, as such, 10 Gy is the most refined model, we cannot rule out the possibility that the models may differ when exploring different stages of the repair process and, in particular, the roles of TGFβ and IL-6, and the development of fibrosis and accumulation of cellular senescence.

Finally, as radiation therapy becomes more targeted in human patients, it is important that radiation-induced changes are recapitulated in animal models. Thus, it is important to recognise and respond to changes in preclinical testing, as reviewed by Verhaegen et al. (2023). However, to date, there has been no evidence that the benefits of proton therapy outweigh the costs for patients with head and neck cancer, and in the UK and USA, the majority of patients are still treated with conventional photon therapy (https://www.england.nhs.uk/publication/proton-beam-therapy-for-head-and-neck-cancer-in-adults/; Blanchard et al., 2018; Leeman et al., 2017; Sio et al., 2016). Thus, in this study, we have explored and validated the acute cellular responses of an animal model that most closely represents the injury that humans experience, highlighting the most refined experimental radiation regime.

Mouse studies

All procedures were approved by the UK Home Office and performed under Project Licence PB5FC9BD (Emmerson). For the purposes of this study, only female mice were used.

Radiation-induced salivary gland injury

Twenty-week-old C57BL/6 mice were anaesthetised using 1 mg/kg medetomidine hydrochloride (Dormitor) and 75 mg/kg ketamine (Ketavet) prepared in 0.9% saline (Thermo Fisher Scientific). The necks of anaesthetised mice were irradiated using a single 137Cs source in a Shepherd Mark-I-68A irradiator (JL Shepherd & Associates), as previously described (Emmerson et al., 2018; McKendrick et al., 2023). Non-irradiated control mice were not anaesthetised. After a 20 min period of anaesthesia, mice were given 1 mg/kg of the reversal agent Antisedan (Vetoquinol UK) and were allowed to recover in cages with additional heat support before returning to normal housing. Mice were provided with soft diet and DietGel (ClearH2O) ad libitum. Mice were euthanised at 3 days post irradiation and SMGs collected for analysis.

Tissue digestion

Individual SMGs (approximately 80 mg) were mechanically minced in 2 ml of RPMI-1640 (Gibco) containing 5% foetal calf serum (FCS; Sigma-Aldrich), 25 µl collagenase-II (23 mg/ml) (Sigma-Aldrich), 25 µl hyaluronidase (40 mg/ml) (Sigma-Aldrich) and 125 µl of 0.1 M CaCl2 using a GentleMACS machine (Miltenyi Biotec) using program A.01 in a GentleMACS C-tube. Tissues were incubated in a shaking incubator at 100 rpm at 37°C for 60 min and subsequently centrifuged at 400 g for 5 min at 4°C. The supernatant was discarded and the pellet was resuspended in 2 ml of RPMI-1640 containing 5% FCS. The solution was filtered through a 70 µm nylon mesh (Thermo Fisher Scientific) and centrifuged at 400 g for 5 min at 4°C. The cell pellet was resuspended in 1 ml of 1× red blood cell lysis buffer (Abcam) for 5 min on ice before centrifugation at 400 g for 5 min at 4°C. The cell pellet was resuspended in 1 ml of pre-warmed (37°C) trypsin containing 0.25% EDTA and incubated at 37°C for 5 min before trituration. This step was repeated three times until single cells were evident. The solution was centrifuged at 400 g for 5 min at 4°C. Then, the cell pellet was resuspended in 1 ml of FACS buffer [Hank's Balanced Salt Solution (Lonza) containing 1% bovine serum albumin (BSA; Sigma-Aldrich) and 2 mM EDTA (Sigma-Aldrich)] and passed through a 20 µm filter-capped 5 ml fluorescence-activated cell sorting (FACS) tube (Falcon) using a 5 ml syringe with a 25G needle.

Flow cytometry and FACS

Equal numbers of cells were resuspended in 100 µl of FACS buffer. Cells were stained with 1:1000 anti-CD16/32 (2.4G2; Biolegend) for 15 min at 4°C to reduce non-specific antibody binding to receptors for IgG. Cells were subsequently stained with conjugated antibodies (Table S1) for 60 min at 4°C in the dark. Samples were washed with FACS buffer and centrifuged at 400 g for 5 min at 4°C before resuspension in 300 µl of FACS buffer. Single stain controls were prepared using OneComp Beads (Thermo Fisher Scientific). Fluorescence minus one (FMO) controls were prepared using cells. SYTOX Green (Thermo Fisher Scientific) was used as a dead cell marker. Samples were analysed using an LSR II flow cytometer (BD Biosciences). Data were analysed using FlowJo software (v9). The dump gate (FITC) comprised SiglecF+ cells, Ly6G+ cells and dead (SYTOX Green+) cells.

Tissue processing for histology

SMGs were fixed for 6-8 h in 4% paraformaldehyde (Thermo Fisher Scientific) at room temperature with constant mixing, followed by three washes with phosphate-buffered saline (PBS; Merck). After fixation, SMGs were processed for the generation of frozen sections by incubating in increasing concentrations of sucrose (10% and 30%) before embedding in OCT (Leica). 12 µm sections were cut using a cryostat (Leica) and stored at −20°C.

Immunofluorescence analysis

Immunofluorescence analysis of tissue sections has been previously described (Emmerson et al., 2018). In brief, tissues were fixed with 4% paraformaldehyde if not previously fixed, and permeabilised with ice-cold acetone/methanol (1:1) for 1 min. Tissues were blocked for 2 h at room temperature with 5% BSA, 5% donkey serum (Merck) and 0.01% Tween 20 in PBS. Salivary glands were incubated with primary antibodies overnight at 4°C. Antibodies are listed in Table S2. Antibodies were detected using donkey Cy2-, Cy3- or Cy5-conjugated secondary Fab fragment antibodies (Jackson Laboratories, 711-546-152, 712-166-150, 705-166-147 and 712-546-153; 1:300), nuclei were stained using Hoechst 33342 (1:1000, Sigma-Aldrich), and slides were mounted using Prolong Gold anti-fade mounting medium. Fluorescence was analysed using a Leica SP8 confocal microscope and ImageJ software (National Institutes of Health).

Histological analysis and cell counts

For immunofluorescence analysis, cells positively stained for markers were counted using ImageJ. Three random fields of view per sample were acquired on a Leica SP8 microscope at 25× or 40× magnification (Nyquist). Images were run through an ImageJ cell counting macro either as single images or, if the file was a z-stack, using the middle image of the stack. Using the macro, images were split into individual channels and the appropriate channel was extracted. Positive cells, such as macrophages, were thresholded and counted according to their size using the ‘Analyse particles’ command. As a quality control, an output file was saved where the counted macrophages were highlighted in green and could be manually checked and confirmed. This macro has previously been published in McKendrick et al. (2023). Fluorescence intensity was measured in ImageJ using the ‘Threshold’ function. Acinar cell size was measured as previously described (Lombaert et al., 2020) using the Image J ‘Measure’ function.

Quantitative PCR analysis

RNA was isolated from a piece (∼3 mm3) of whole SMG tissue by sonicating in lysis buffer from the RNAqueous Micro Kit (Life Technologies) and following the manufacturer's instructions. Total RNA was DNase treated (Life Technologies) prior to cDNA synthesis (First Strand Synthesis Kit, Thermo Fisher Scientific). Quantitative PCR (qPCR) with SYBR Green was performed using 5 ng cDNA and primers [designed using Primer3 (https://www.primer3plus.com/index.html) and Beacon Designer software (https://www.premierbiosoft.com/molecular_beacons/) or described on PrimerBank (http://pga.mgh.harvard.edu/primerbank/)]. Primer sequences are listed in Table S3. Melt curves and primer efficiency were determined as previously described (Hoffman et al., 2002). Gene expression was normalised to that of the housekeeping gene Gapdh and to a selected unirradiated control mouse sample. Reactions were run in triplicate.

Statistical tests

Normal distribution was assessed using the D'Agostino–Pearson omnibus test. Data were analysed for statistical significance using one-way ANOVA (multiple groups) with post hoc testing performed using Tukey’s Q test (GraphPad Prism). For multiple testing, we used a false discovery rate of 0.05. All graphs show the mean±standard deviation (s.d.) as indicated in the figure legends. Statistical tests used for each experiment are also indicated in all figure legends.

The authors would like to thank Dr Fiona Rossi and the Centre for Regenerative Medicine (CRM) Flow Cytometry Facility; Dr Matthieu Vermeren and the CRM Imaging Facility; and the University of Edinburgh Little France animal facility.

Author contributions

Conceptualization: E.E.; Methodology: A.L.J., S.S.E., J.G.M., E.E.; Software: S.S.E.; Validation: A.L.J., S.S.E., J.G.M., E.E.; Formal analysis: A.L.J., S.S.E., J.G.M., L.M.H., E.E.; Investigation: A.L.J., S.S.E., J.G.M., E.M., E.E.; Resources: E.E.; Writing - original draft: E.E.; Writing - review & editing: A.L.J., S.S.E., L.M.H., E.E.; Visualization: A.L.J., S.S.E., E.E.; Supervision: E.E.; Project administration: E.E.; Funding acquisition: E.E.

Funding

S.S.E. and E.M. are funded by the Wellcome Trust grant 108906/Z/15/Z. J.G.M., L.M.H. and E.E. are funded by UK Research and Innovation (UKRI)/Medical Research Council (MRC) grants (MR/W004763/1, MR/X018733/1 and MR/S005544/1, respectively). J.G.M. is also funded by the University of Edinburgh Chancellor's Fellowship studentship, and E.E. is funded by a Chancellor's Fellowship from the University of Edinburgh. Open Access funding provided by University of Edinburgh. Deposited in PMC for immediate release.

Data availability

All relevant data can be found within the article and its supplementary information.

Affandi
,
T.
,
Ohm
,
A. M.
,
Gaillard
,
D.
,
Haas
,
A.
and
Reyland
,
M. E.
(
2021
).
Tyrosine kinase inhibitors protect the salivary gland from radiation damage by increasing DNA double-strand break repair
.
J. Biol. Chem.
296
,
100401
.
Anscher
,
M. S.
,
Kong
,
F.-M.
,
Murase
,
T.
and
Jirtle
,
R. L.
(
1995
).
Short communication: normal tissue injury after cancer therapy is a local response exacerbated by an endocrine effect of TGF beta
.
Br. J. Radiol.
68
,
331
-
333
.
Avila
,
J. L.
,
Grundmann
,
O.
,
Burd
,
R.
and
Limesand
,
K. H.
(
2009
).
Radiation-induced salivary gland dysfunction results from p53-dependent apoptosis
.
Int. J. Radiat. Oncol. Biol. Phys.
73
,
523
-
529
.
Bentzen
,
S. M.
(
2006
).
Preventing or reducing late side effects of radiation therapy: radiobiology meets molecular pathology
.
Nat. Rev. Cancer
6
,
702
-
713
.
Bilal
,
B. U.
,
Kocak
,
Z.
,
Cosar
,
R.
,
Altaner
,
S.
and
Kocak
,
M.
(
2014
).
Histopathological investigation of differences between female and male rats on radiation-induced acute lung toxicity
.
J. Cancer Res. Ther.
2
,
84
-
88
.
Blanchard
,
P.
,
Gunn
,
G. B.
,
Lin
,
A.
,
Foote
,
R. L.
,
Lee
,
N. Y.
and
Frank
,
S. J.
(
2018
).
Proton therapy for head and neck cancers
.
Semin. Radiat. Oncol.
28
,
53
-
63
.
Bourhis
,
J.
,
Sire
,
C.
,
Graff
,
P.
,
Grégoire
,
V.
,
Maingon
,
P.
,
Calais
,
G.
,
Gery
,
B.
,
Martin
,
L.
,
Alfonsi
,
M.
,
Desprez
,
P.
et al. 
(
2012
).
Concomitant chemoradiotherapy versus acceleration of radiotherapy with or without concomitant chemotherapy in locally advanced head and neck carcinoma (GORTEC 99-02): an open-label phase 3 randomised trial
.
Lancet Oncol.
13
,
145
-
153
.
Bralic
,
M.
,
Muhvic-Urek
,
M.
,
Stemberga
,
V.
,
Golemac
,
M.
,
Jurkovic
,
S.
,
Borcic
,
J.
,
Braut
,
A.
and
Tomac
,
J.
(
2005
).
Cell death and cell proliferation in mouse submandibular gland during early post-irradiation phase
.
Acta Med. Okayama
59
,
153
-
159
.
Brown
,
C. T.
,
Nam
,
K.
,
Zhang
,
Y.
,
Qiu
,
Y.
,
Dean
,
S. M.
,
Dos Santos
,
H. T.
,
Lei
,
P.
,
Andreadis
,
S. T.
and
Baker
,
O. J.
(
2020
).
Sex-dependent regeneration patterns in mouse submandibular glands
.
J. Histochem. Cytochem.
68
,
305
-
318
.
Chai
,
Y.
,
Klauser
,
D. K.
,
Denny
,
P. A.
and
Denny
,
P. C.
(
1993
).
Proliferative and structural differences between male and female mouse submandibular glands
.
Anat. Rec.
235
,
303
-
311
.
Chang
,
C. H.
,
Zhang
,
M.
,
Rajapakshe
,
K.
,
Coarfa
,
C.
,
Edwards
,
D.
,
Huang
,
S.
and
Rosen
,
J. M.
(
2015
).
Mammary stem cells and tumor-initiating cells are more resistant to apoptosis and exhibit increased DNA Repair activity in response to DNA damage
.
Stem Cell Reports
5
,
378
-
391
.
Chang
,
S.
,
Hu
,
L.
,
Xu
,
Y.
,
Li
,
X.
,
Ma
,
L.
,
Feng
,
X.
,
Wang
,
J.
,
Zhang
,
C.
and
Wang
,
S.
(
2019
).
Inorganic nitrate alleviates total body irradiation-induced systemic damage by decreasing reactive oxygen species levels
.
Int. J. Radiat. Oncol. Biol. Phys.
103
,
945
-
957
.
Cheng
,
J. N.
,
Luo
,
W.
,
Sun
,
C.
,
Jin
,
Z.
,
Zeng
,
X.
,
Alexander
,
P. B.
,
Gong
,
Z.
,
Xia
,
X.
,
Ding
,
X.
,
Xu
,
S.
et al. 
(
2021
).
Radiation-induced eosinophils improve cytotoxic T lymphocyte recruitment and response to immunotherapy
.
Sci. Adv.
7
,
eabc7609
.
Choi
,
J. H.
,
Wu
,
H.-G.
,
Jung
,
K. C.
,
Lee
,
S. H.
and
Kwon
,
E. K.
(
2009
).
Apoptosis and expression of AQP5 and TGF-beta in the irradiated rat submandibular gland
.
Cancer Res. Treat
41
,
145
-
154
.
Cotrim
,
A. P.
,
Sowers
,
A. L.
,
Lodde
,
B. M.
,
Vitolo
,
J. M.
,
Kingman
,
A.
,
Russo
,
A.
,
Mitchell
,
J. B.
and
Baum
,
B. J.
(
2005
).
Kinetics of tempol for prevention of xerostomia following head and neck irradiation in a mouse model
.
Clin. Cancer Res.
11
,
7564
-
7568
.
Cotrim
,
A. P.
,
Sowers
,
A.
,
Mitchell
,
J. B.
and
Baum
,
B. J.
(
2007
).
Prevention of irradiation-induced salivary hypofunction by microvessel protection in mouse salivary glands
.
Mol. Ther.
15
,
2101
-
2106
.
Dai
,
Y. H.
,
Li
,
X.-
,
Dong
,
D.-
,
Gu
,
H.-
,
Kong
,
C.-
and
Xu
,
Z.
(
2019
).
P27 Promotes TGF-β-mediated pulmonary fibrosis via interacting with MTORC2
.
Can. Respir. J.
2019
,
7157861
.
De la Cal
,
C.
,
Fernández-Solari
,
J.
,
Mohn
,
C. E.
,
Prestifilippo
,
J. P.
,
Pugnaloni
,
A.
,
Medina
,
V. A.
and
Elverdin
,
J. C.
(
2012
).
Radiation produces irreversible chronic dysfunction in the submandibular glands of the rat
.
Open Dent J.
6
,
8
-
13
.
Delporte
,
C.
and
Steinfeld
,
S.
(
2006
).
Distribution and roles of aquaporins in salivary glands
.
Biochim. Biophys. Acta
1758
,
1061
-
1070
.
Dileto
,
C. L.
and
Travis
,
E. L.
(
1996
).
Fibroblast radiosensitivity in vitro and lung fibrosis in vivo: comparison between a fibrosis-prone and fibrosis-resistant mouse strain
.
Radiat. Res.
146
,
61
-
67
.
Down
,
J. D.
and
Yanch
,
J. C.
(
2010
).
Identifying the high radiosensitivity of the lungs of C57L mice in a model of total-body irradiation and bone marrow transplantation
.
Radiat. Res.
174
,
258
-
263
.
Down
,
J. D.
,
Nicholas
,
D.
and
Steel
,
G. G.
(
1986
).
Lung damage after hemithoracic irradiation: dependence on mouse strain
.
Radiother. Oncol.
6
,
43
-
50
.
Eisbruch
,
A.
,
Kim
,
H. M.
,
Terrell
,
J. E.
,
Marsh
,
L. H.
,
Dawson
,
L. A.
and
Ship
,
J. A.
(
2001
).
Xerostomia and its predictors following parotid-sparing irradiation of head-and-neck cancer
.
Int. J. Radiat. Oncol. Biol. Phys.
50
,
695
-
704
.
Emmerson
,
E.
and
Knox
,
S. M.
(
2018
).
Salivary gland stem cells: A review of development, regeneration and cancer
.
Genesis
56
,
e23211
.
Emmerson
,
E.
,
May
,
A. J.
,
Nathan
,
S.
,
Cruz-Pacheco
,
N.
,
Lizama
,
C. O.
,
Maliskova
,
L.
,
Zovein
,
A. C.
,
Shen
,
Y.
,
Muench
,
M. O.
and
Knox
,
S. M.
(
2017
).
SOX2 regulates acinar cell development in the salivary gland
.
Elife
6
,
e26620
.
Emmerson
,
E.
,
May
,
A. J.
,
Berthoin
,
L.
,
Cruz-Pacheco
,
N.
,
Nathan
,
S.
,
Mattingly
,
A. J.
,
Chang
,
J. L.
,
Ryan
,
W. R.
,
Tward
,
A. D.
and
Knox
,
S. M.
(
2018
).
Salivary glands regenerate after radiation injury through SOX2-mediated secretory cell replacement
.
EMBO Mol. Med.
10
,
e8051
.
Frangogiannis
,
N.
(
2020
).
Transforming growth factor-β in tissue fibrosis
.
J. Exp. Med.
217
,
e20190103
.
Franzén
,
L.
,
Funegård
,
U.
,
Ericson
,
T.
and
Henriksson
,
R.
(
1992
).
Parotid gland function during and following radiotherapy of malignancies in the head and neck. A consecutive study of salivary flow and patient discomfort
.
Eur. J. Cancer
28
,
457
-
462
.
Funegård
,
U.
,
Franzén
,
L.
,
Ericson
,
T.
and
Henriksson
,
R.
(
1994
).
Parotid saliva composition during and after irradiation of head and neck cancer
.
Eur. J. Cancer B Oral. Oncol.
30b
,
230
-
233
.
Gilman
,
K. E.
,
Camden
,
J. M.
,
Klein
,
R. R.
,
Zhang
,
Q.
,
Weisman
,
G. A.
and
Limesand
,
K. H.
(
2019
).
P2X7 receptor deletion suppresses γ-radiation-induced hyposalivation
.
Am. J. Physiol. Regul. Integr. Comp. Physiol.
316
,
R687
-
r696
.
Groves
,
A. M.
,
Misra
,
R.
,
Clair
,
G.
,
Hernady
,
E.
,
Olson
,
H.
,
Orton
,
D.
,
Finkelstein
,
J.
,
Marples
,
B.
and
Johnston
,
C. J.
(
2023
).
Influence of the irradiated pulmonary microenvironment on macrophage and T cell dynamics
.
Radiother. Oncol.
183
,
109543
.
Grundmann
,
O.
,
Mitchell
,
G. C.
and
Limesand
,
K. H.
(
2009
).
Sensitivity of salivary glands to radiation: from animal models to therapies
.
J. Dent. Res.
88
,
894
-
903
.
Grundmann
,
O.
,
Fillinger
,
J. L.
,
Victory
,
K. R.
,
Burd
,
R.
and
Limesand
,
K. H.
(
2010
).
Restoration of radiation therapy-induced salivary gland dysfunction in mice by post therapy IGF-1 administration
.
BMC Cancer
10
,
417
.
Gupta
,
A.
,
Hunt
,
C. R.
,
Chakraborty
,
S.
,
Pandita
,
R. K.
,
Yordy
,
J.
,
Ramnarain
,
D. B.
,
Horikoshi
,
N.
and
Pandita
,
T. K.
(
2014
).
Role of 53BP1 in the regulation of DNA double-strand break repair pathway choice
.
Radiat. Res.
181
,
1
-
8
.
Hakim
,
S. G.
,
Ribbat
,
J.
,
Berndt
,
A.
,
Richter
,
P.
,
Kosmehl
,
H.
,
Benedek
,
G. A.
,
Jacobsen
,
H. C.
,
Trenkle
,
T.
,
Sieg
,
P.
and
Rades
,
D.
(
2011
).
Expression of Wnt-1, TGF-β and related cell-cell adhesion components following radiotherapy in salivary glands of patients with manifested radiogenic xerostomia
.
Radiother. Oncol.
101
,
93
-
99
.
Hawkins
,
P. G.
,
Lee
,
J. Y.
,
Mao
,
Y.
,
Li
,
P.
,
Green
,
M.
,
Worden
,
F. P.
,
Swiecicki
,
P. L.
,
Mierzwa
,
M. L.
,
Spector
,
M. E.
,
Schipper
,
M. J.
et al. 
(
2018
).
Sparing all salivary glands with IMRT for head and neck cancer: Longitudinal study of patient-reported xerostomia and head-and-neck quality of life
.
Radiother. Oncol.
126
,
68
-
74
.
Henson
,
B. S.
,
Eisbruch
,
A.
,
D'Hondt
,
E.
and
Ship
,
J. A.
(
1999
).
Two-year longitudinal study of parotid salivary flow rates in head and neck cancer patients receiving unilateral neck parotid-sparing radiotherapy treatment
.
Oral Oncol.
35
,
234
-
241
.
Hoffman
,
M. P.
,
Kidder
,
B. L.
,
Steinberg
,
Z. L.
,
Lakhani
,
S.
,
Ho
,
S.
,
Kleinman
,
H. K.
and
Larsen
,
M.
(
2002
).
Gene expression profiles of mouse submandibular gland development: FGFR1 regulates branching morphogenesis in vitro through BMP- and FGF-dependent mechanisms
.
Development
129
,
5767
-
5778
.
Hudson
,
W. H.
,
Olson
,
J. J.
and
Sudmeier
,
L. J.
(
2023
).
Immune microenvironment remodeling after radiation of a progressing brain metastasis
.
Cell Rep. Med.
4
,
101054
.
Humphries
,
M. J.
,
Limesand
,
K. H.
,
Schneider
,
J. C.
,
Nakayama
,
K. I.
,
Anderson
,
S. M.
and
Reyland
,
M. E.
(
2006
).
Suppression of apoptosis in the protein kinase Cdelta null mouse in vivo
.
J. Biol. Chem.
281
,
9728
-
9737
.
Insinga
,
A.
,
Cicalese
,
A.
,
Faretta
,
M.
,
Gallo
,
B.
,
Albano
,
L.
,
Ronzoni
,
S.
,
Furia
,
L.
,
Viale
,
A.
and
Pelicci
,
P. G.
(
2013
).
DNA damage in stem cells activates p21, inhibits p53, and induces symmetric self-renewing divisions
.
Proc. Natl. Acad. Sci. USA
110
,
3931
-
3936
.
Jackson
,
I. L.
,
Vujaskovic
,
Z.
and
Down
,
J. D.
(
2010
).
Revisiting strain-related differences in radiation sensitivity of the mouse lung: recognizing and avoiding the confounding effects of pleural effusions
.
Radiat. Res.
173
,
10
-
20
.
Kim
,
J. W.
,
Kim
,
J. M.
,
Choi
,
M. E.
,
Kim
,
S.-K.
,
Kim
,
Y.-M.
and
Choi
,
J.-S.
(
2019
).
Adipose-derived mesenchymal stem cells regenerate radioiodine-induced salivary gland damage in a murine model
.
Sci. Rep.
9
,
15752
.
Kizu
,
Y.
,
Sakurai
,
H.
,
Katagiri
,
S.
,
Shinozaki
,
N.
,
Ono
,
M.
,
Tsubota
,
K.
and
Saito
,
J.
(
1996
).
Immunohistological analysis of tumour growth factor beta 1 expression in normal and inflamed salivary glands
.
J. Clin. Pathol.
49
,
728
-
732
.
Lawal
,
A. O.
,
Adisa
,
A. O.
,
Kolude
,
B.
and
Adeyemi
,
B.
(
2015
).
Malignant salivary gland tumours of the head and neck region: a single institutions review
.
Pan. Afr. Med. J.
20
,
121
.
Lechelt
,
L. A.
,
Rieger
,
J. M.
,
Cowan
,
K.
,
Debenham
,
B. J.
,
Krewski
,
B.
,
Nayar
,
S.
,
Regunathan
,
A.
,
Seikaly
,
H.
,
Singh
,
A. E.
and
Laupacis
,
A.
(
2018
).
Top 10 research priorities in head and neck cancer: Results of an Alberta priority setting partnership of patients, caregivers, family members, and clinicians
.
Head Neck
40
,
544
-
554
.
Lee
,
J.
,
Kim
,
M.-R.
,
Kim
,
H.-J.
,
An
,
Y. S.
and
Yi
,
J. Y.
(
2016
).
TGF-β1 accelerates the DNA damage response in epithelial cells via Smad signaling
.
Biochem. Biophys. Res. Commun.
476
,
420
-
425
.
Leeman
,
J. E.
,
Romesser
,
P. B.
,
Zhou
,
Y.
,
McBride
,
S.
,
Riaz
,
N.
,
Sherman
,
E.
,
Cohen
,
M. A.
,
Cahlon
,
O.
and
Lee
,
N.
(
2017
).
Proton therapy for head and neck cancer: expanding the therapeutic window
.
Lancet Oncol.
18
,
e254
-
e265
.
Lertkiatmongkol
,
P.
,
Liao
,
D.
,
Mei
,
H.
,
Hu
,
Y.
and
Newman
,
P. J.
(
2016
).
Endothelial functions of platelet/endothelial cell adhesion molecule-1 (CD31)
.
Curr. Opin Hematol.
23
,
253
-
259
.
Li
,
H.
,
Chen
,
X.
,
Zeng
,
W.
,
Zhou
,
W.
,
Zhou
,
Q.
,
Wang
,
Z.
,
Jiang
,
W.
,
Xie
,
B.
and
Sun
,
L.-Q.
(
2020
).
Radiation-enhanced expression of CCL22 in nasopharyngeal carcinoma is associated with CCR4(+) CD8 T cell recruitment
.
Int. J. Radiat. Oncol. Biol. Phys.
108
,
126
-
139
.
Li
,
J.
,
Sudiwala
,
S.
,
Berthoin
,
L.
,
Mohabbat
,
S.
,
Gaylord
,
E. A.
,
Sinada
,
H.
,
Cruz Pacheco
,
N.
,
Chang
,
J. C.
,
Jeon
,
O.
,
Lombaert
,
I. M. A.
et al. 
(
2022
).
Long-term functional regeneration of radiation-damaged salivary glands through delivery of a neurogenic hydrogel
.
Sci. Adv.
8
,
eadc8753
.
Limesand
,
K. H.
,
Schwertfeger
,
K. L.
and
Anderson
,
S. M.
(
2006
).
MDM2 is required for suppression of apoptosis by activated Akt1 in salivary acinar cells
.
Mol. Cell. Biol.
26
,
8840
-
8856
.
Limesand
,
K. H.
,
Avila
,
J. L.
,
Victory
,
K.
,
Chang
,
H.-H.
,
Shin
,
Y. J.
,
Grundmann
,
O.
and
Klein
,
R. R.
(
2010
).
Insulin-like growth factor-1 preserves salivary gland function after fractionated radiation
.
Int. J. Radiat. Oncol. Biol. Phys.
78
,
579
-
586
.
Liu
,
X.
,
Ong
,
H. L.
and
Ambudkar
,
I.
(
2018
).
TRP Channel involvement in salivary glands-some good, some bad
.
Cells
7
,
74
.
Lombaert
,
I. M.
,
Brunsting
,
J. F.
,
Wierenga
,
P. K.
,
Faber
,
H.
,
Stokman
,
M. A.
,
Kok
,
T.
,
Visser
,
W. H.
,
Kampinga
,
H. H.
,
de Haan
,
G.
and
Coppes
,
R. P.
(
2008
).
Rescue of salivary gland function after stem cell transplantation in irradiated glands
.
PLoS ONE
3
,
e2063
.
Lombaert
,
I. M. A.
,
Patel
,
V. N.
,
Jones
,
C. E.
,
Villier
,
D. C.
,
Canada
,
A. E.
,
Moore
,
M. R.
,
Berenstein
,
E.
,
Zheng
,
C.
,
Goldsmith
,
C. M.
,
Chorini
,
J. A.
et al. 
(
2020
).
CERE-120 Prevents irradiation-induced hypofunction and restores immune homeostasis in porcine salivary glands
.
Mol. Ther. Methods Clin. Dev.
18
,
839
-
855
.
Makkonen
,
T. A.
,
Tenovuo
,
J.
,
Vilja
,
P.
and
Heimdahl
,
A.
(
1986
).
Changes in the protein composition of whole saliva during radiotherapy in patients with oral or pharyngeal cancer
.
Oral. Surg. Oral. Med. Oral. Pathol.
62
,
270
-
275
.
Marmary
,
Y.
,
Adar
,
R.
,
Gaska
,
S.
,
Wygoda
,
A.
,
Maly
,
A.
,
Cohen
,
J.
,
Eliashar
,
R.
,
Mizrachi
,
L.
,
Orfaig-Geva
,
C.
,
Baum
,
B. J.
et al. 
(
2016
).
Radiation-induced loss of salivary gland function is driven by cellular senescence and prevented by IL6 modulation
.
Cancer Res.
76
,
1170
-
1180
.
Martin
,
K. L.
,
Hill
,
G. A.
,
Klein
,
R. R.
,
Arnett
,
D. G.
,
Burd
,
R.
and
Limesand
,
K. H.
(
2012
).
Prevention of radiation-induced salivary gland dysfunction utilizing a CDK inhibitor in a mouse model
.
PLoS One
7
,
e51363
.
May
,
A. J.
,
Cruz-Pacheco
,
N.
,
Emmerson
,
E.
,
Gaylord
,
E. A.
,
Seidel
,
K.
,
Nathan
,
S.
,
Muench
,
M. O.
,
Klein
,
O.
and
Knox
,
S. M.
(
2018
).
Diverse progenitor cells preserve salivary gland ductal architecture after radiation-induced damage
.
Development
145
,
dev166363
.
McKendrick
,
J. G.
,
Jones
,
G.-R.
,
Elder
,
S. S.
,
Watson
,
E.
,
T'Jonck
,
W.
,
Mercer
,
E.
,
Magalhaes
,
M. S.
,
Rocchi
,
C.
,
Hegarty
,
L. M.
,
Johnson
,
A. L.
et al. 
(
2023
).
CSF1R-dependent macrophages in the salivary gland are essential for epithelial regeneration after radiation-induced injury
.
Sci. Immunol.
8
,
eadd4374
.
Meng
,
X. M.
,
Nikolic-Paterson
,
D. J.
and
Lan
,
H. Y.
(
2016
).
TGF-β: the master regulator of fibrosis
.
Nat. Rev. Nephrol.
12
,
325
-
338
.
Meyer
,
S.
,
Chibly
,
A. M.
,
Burd
,
R.
and
Limesand
,
K. H.
(
2017a
).
Insulin-like growth factor-1-mediated DNA repair in irradiated salivary glands is sirtuin-1 dependent
.
J. Dent. Res.
96
,
225
-
232
.
Meyer
,
J. E.
,
Finnberg
,
N. K.
,
Chen
,
L.
,
Cvetkovic
,
D.
,
Wang
,
B.
,
Zhou
,
L.
,
Dong
,
Y.
,
Hallman
,
M. A.
,
Ma
,
C.-M. C.
and
El-Deiry
,
W. S.
(
2017b
).
Tissue TGF-β expression following conventional radiotherapy and pulsed low-dose-rate radiation
.
Cell Cycle
16
,
1171
-
1174
.
Mezi
,
S.
,
Pomati
,
G.
,
Botticelli
,
A.
,
De Felice
,
F.
,
Musio
,
D.
,
della Monaca
,
M.
,
Amirhassankhani
,
S.
,
Vullo
,
F.
,
Cerbelli
,
B.
,
Carletti
,
R.
et al. 
(
2020
).
Primary squamous cell carcinoma of major salivary gland: “Sapienza Head and Neck Unit” clinical recommendations
.
Rare Tumors
12
,
2036361320973526
.
Meziani
,
L.
,
Mondini
,
M.
,
Petit
,
B.
,
Boissonnas
,
A.
,
Thomas de Montpreville
,
V.
,
Mercier
,
O.
,
Vozenin
,
M.-C.
and
Deutsch
,
E.
(
2018
).
CSF1R inhibition prevents radiation pulmonary fibrosis by depletion of interstitial macrophages
.
Eur. Respir. J.
51
,
1702120
.
Muhvic-Urek
,
M.
,
Bralic
,
M.
,
Curic
,
C.
,
Pezelj-Ribaric
,
S.
,
Borcic
,
J.
and
Tomac
,
J.
(
2006
).
Imbalance between apoptosis and proliferation causes late radiation damage of salivary gland in mouse
.
Physiol. Res.
55
,
89
-
95
.
Mukaibo
,
T.
,
Gao
,
X.
,
Yang
,
N. Y.
,
Oei
,
M. S.
,
Nakamoto
,
T.
and
Melvin
,
J. E.
(
2019
).
Sexual dimorphisms in the transcriptomes of murine salivary glands
.
FEBS Open Bio.
9
,
947
-
958
.
Nagler
,
R. M.
(
2002
).
The enigmatic mechanism of irradiation-induced damage to the major salivary glands
.
Oral Dis.
8
,
141
-
146
.
Nam
,
K.
,
Maruyama
,
C. L.
,
Trump
,
B. G.
,
Buchmann
,
L.
,
Hunt
,
J. P.
,
Monroe
,
M. M.
and
Baker
,
O. J.
(
2016
).
Post-Irradiated Human Submandibular Glands Display High Collagen Deposition, Disorganized Cell Junctions, and an Increased Number of Adipocytes
.
J. Histochem. Cytochem.
64
,
343
-
352
.
Nolan
,
E.
,
Bridgeman
,
V. L.
,
Ombrato
,
L.
,
Karoutas
,
A.
,
Rabas
,
N.
,
Sewnath
,
C. A. N.
,
Vasquez
,
M.
,
Rodrigues
,
F. S.
,
Horswell
,
S.
,
Faull
,
P.
et al. 
(
2022
).
Radiation exposure elicits a neutrophil-driven response in healthy lung tissue that enhances metastatic colonization
.
Nat. Cancer
3
,
173
-
187
.
Page
,
B.
,
Han
,
P.
,
Peng
,
L. C.
,
Cheng
,
Z.
,
Harkness
,
J.
,
Shen
,
C.
,
Choflet
,
A.
,
Cecil
,
E.
,
Hui
,
X.
,
Schmitt
,
N. C.
et al. 
(
2018
).
Gender differences in radiation therapy effects in male and female patients with head and neck cancer
.
Int. J. Radiat. Oncol. Biol. Phys.
100
,
1400
-
1401
.
Peng
,
X.
,
Wu
,
Y.
,
Brouwer
,
U.
,
van Vliet
,
T.
,
Wang
,
B.
,
Demaria
,
M.
,
Barazzuol
,
L.
and
Coppes
,
R. P.
(
2020
).
Cellular senescence contributes to radiation-induced hyposalivation by affecting the stem/progenitor cell niche
.
Cell Death Dis.
11
,
854
.
Perez
,
E. C.
,
Rodgers
,
S. P.
,
Inoue
,
T.
,
Pedersen
,
S. E.
,
Leasure
,
J. L.
and
Gaber
,
M. W.
(
2018
).
Olfactory memory impairment differs by sex in a rodent model of pediatric radiotherapy
.
Front. Behav. Neurosci.
12
,
158
.
Pinkstaff
,
C. A.
(
1998
).
Salivary gland sexual dimorphism: a brief review
.
Eur. J. Morphol.
36
,
31
-
34
.
Ren
,
J.
,
Huang
,
R.
,
Li
,
Y.
,
Chen
,
R.
,
Tian
,
H.
and
Liu
,
C.
(
2022
).
Radioprotective effects and mechanism of HL-003 on radiation-induced salivary gland damage in mice
.
Sci. Rep.
12
,
8419
.
Rodman
,
S. N.
,
Kluz
,
P. N.
,
Hines
,
M. R.
,
Oberley-Deegan
,
R. E.
and
Coleman
,
M. C.
(
2022
).
Sex-based differences in the severity of radiation-induced arthrofibrosis
.
J. Orthop. Res.
40
,
2586
-
2596
.
Rosen
,
I. I.
,
Fischer
,
T. A.
,
Antolak
,
J. A.
,
Starkschall
,
G.
,
Travis
,
E. L.
,
Tucker
,
S. L.
,
Hogstrom
,
K. R.
,
Cox
,
J. D.
and
Komaki
,
R.
(
2001
).
Correlation between lung fibrosis and radiation therapy dose after concurrent radiation therapy and chemotherapy for limited small cell lung cancer
.
Radiology
221
,
614
-
622
.
Shen
,
Z.
,
Wang
,
J.
,
Huang
,
Q.
,
Shi
,
Y.
,
Wei
,
Z.
,
Zhang
,
X.
,
Qiu
,
Y.
,
Zhang
,
M.
,
Wang
,
Y.
,
Qin
,
W.
et al. 
(
2018
).
Genetic modification to induce CXCR2 overexpression in mesenchymal stem cells enhances treatment benefits in radiation-induced oral mucositis
.
Cell Death Dis.
9
,
229
.
Shibata
,
A.
and
Jeggo
,
P. A.
(
2020
).
Roles for 53BP1 in the repair of radiation-induced DNA double strand breaks
.
DNA Repair (Amst)
93
,
102915
.
Shinde-Jadhav
,
S.
,
Mansure
,
J. J.
,
Rayes
,
R. F.
,
Marcq
,
G.
,
Ayoub
,
M.
,
Skowronski
,
R.
,
Kool
,
R.
,
Bourdeau
,
F.
,
Brimo
,
F.
,
Spicer
,
J.
et al. 
(
2021
).
Role of neutrophil extracellular traps in radiation resistance of invasive bladder cancer
.
Nat. Commun.
12
,
2776
.
Sio
,
T. T.
,
Lin
,
H.-K.
,
Shi
,
Q.
,
Gunn
,
G. B.
,
Cleeland
,
C. S.
,
Lee
,
J. J.
,
Hernandez
,
M.
,
Blanchard
,
P.
,
Thaker
,
N. G.
,
Phan
,
J.
et al. 
(
2016
).
Intensity modulated proton therapy versus intensity modulated photon radiation therapy for oropharyngeal cancer: first comparative results of patient-reported outcomes
.
Int. J. Radiat. Oncol. Biol. Phys.
95
,
1107
-
1114
.
Snyder
,
M.
,
Burmeister
,
J.
,
Joiner
,
M.
,
Meyer
,
J.
,
Tait
,
L.
,
Cohen
,
S.
,
McSpadden
,
E.
and
Konski
,
A.
(
2012
).
Potential gender differences in a normal tissue complication probability model for heart toxicity during radiation therapy for esophageal Cancer
.
Int. J. Radiat. Oncol. Biol. Phys.
84
,
S757
-
S758
.
Song
,
Y.
,
Uchida
,
H.
,
Sharipol
,
A.
,
Piraino
,
L.
,
Mereness
,
J. A.
,
Ingalls
,
M. H.
,
Rebhahn
,
J.
,
Newlands
,
S. D.
,
DeLouise
,
L. A.
,
Ovitt
,
C. E.
et al. 
(
2021
).
Development of a functional salivary gland tissue chip with potential for high-content drug screening
.
Commun Biol
4
,
361
.
Sorokina
,
S.
,
Markova
,
E.
,
Gursky
,
J.
,
Dobrovodsky
,
J.
and
Belyaev
,
I.
(
2013
).
Relative biological efficiency of protons at low and therapeutic doses in induction of 53BP1/γH2AX foci in lymphocytes from umbilical cord blood
.
Int. J. Radiat. Biol.
89
,
716
-
723
.
Sotiropoulou
,
P. A.
,
Candi
,
A.
,
Mascré
,
G.
,
De Clercq
,
S.
,
Youssef
,
K. K.
,
Lapouge
,
G.
,
Dahl
,
E.
,
Semeraro
,
C.
,
Denecker
,
G.
,
Marine
,
J.-C.
et al. 
(
2010
).
Bcl-2 and accelerated DNA repair mediates resistance of hair follicle bulge stem cells to DNA-damage-induced cell death
.
Nat. Cell Biol.
12
,
572
-
582
.
Soysa
,
R.
,
Lampert
,
S.
,
Yuen
,
S.
,
Douglass
,
A. N.
,
Li
,
W.
,
Pfeffer
,
K.
and
Crispe
,
I. N.
(
2019
).
Fetal origin confers radioresistance on liver macrophages via p21(cip1/WAF1)
.
J. Hepatol.
71
,
553
-
562
.
Spiegelberg
,
L.
,
Swagemakers
,
S. M. A.
,
van Ijcken
,
W. F. J.
,
Oole
,
E.
,
Wolvius
,
E. B.
,
Essers
,
J.
and
Braks
,
J. A. M.
(
2014
).
Gene expression analysis reveals inhibition of radiation-induced TGFβ-signaling by hyperbaric oxygen therapy in mouse salivary glands
.
Mol. Med.
20
,
257
-
269
.
Sullivan
,
C. A.
,
Haddad
,
R. I.
,
Tishler
,
R. B.
,
Mahadevan
,
A.
and
Krane
,
J. F.
(
2005
).
Chemoradiation-induced cell loss in human submandibular glands
.
Laryngoscope
115
,
958
-
964
.
Takeda
,
I.
,
Kizu
,
Y.
,
Yoshitaka
,
O.
,
Saito
,
I.
and
Yamane
,
G.-
(
2003
).
Possible role of nitric oxide in radiation-induced salivary gland dysfunction
.
Radiat. Res.
159
,
465
-
470
.
Teshima
,
T. H.
,
Wells
,
K. L.
,
Lourenço
,
S. V.
and
Tucker
,
A. S.
(
2016
).
Apoptosis in early salivary gland duct morphogenesis and lumen formation
.
J. Dent. Res.
95
,
277
-
283
.
Teymoortash
,
A.
,
Tiemann
,
M.
,
Schrader
,
C.
,
Hartmann
,
O.
and
Werner
,
J. A.
(
2003
).
Transforming growth factor beta in chronic obstructive sialadenitis of human submandibular gland
.
Arch. Oral Biol.
48
,
111
-
116
.
Thomson
,
D. J.
,
Palma
,
D.
,
Guckenberger
,
M.
,
Balermpas
,
P.
,
Beitler
,
J. J.
,
Blanchard
,
P.
,
Brizel
,
D.
,
Budach
,
W.
,
Caudell
,
J.
,
Corry
,
J.
et al. 
(
2020
).
Practice recommendations for risk-adapted head and neck cancer radiotherapy during the COVID-19 pandemic: An ASTRO-ESTRO consensus statement
.
Radiother. Oncol.
151
,
314
-
321
.
Uchida
,
H.
,
Ingalls
,
M. H.
,
Maruyama
,
E. O.
,
Johnston
,
C. J.
,
Hernady
,
E.
,
Faustoferri
,
R. C.
and
Ovitt
,
C. E.
(
2022
).
Short-term and bystander effects of radiation on murine submandibular glands
.
Dis Model. Mech.
15
,
dmm049570
.
Urek
,
M. M.
,
Bralic
,
M.
,
Tomac
,
J.
,
Borcic
,
J.
,
Uhac
,
I.
,
Glazar
,
I.
,
Antonic
,
R.
and
Ferreri
,
S.
(
2005
).
Early and late effects of X-irradiation on submandibular gland: a morphological study in mice
.
Arch. Med. Res.
36
,
339
-
343
.
van Timmeren
,
J. E.
,
Chamberlain
,
M.
,
Bogowicz
,
M.
,
Ehrbar
,
S.
,
Dal Bello
,
R.
,
Garcia Schüler
,
H.
,
Krayenbuehl
,
J.
,
Wilke
,
L.
,
Andratschke
,
N.
,
Guckenberger
,
M.
et al. 
(
2021
).
MR-Guided adaptive radiotherapy for head and neck cancer: prospective evaluation of migration and anatomical changes of the major salivary glands
.
Cancers (Basel)
13
,
5404
.
Varghese
,
J. J.
,
Schmale
,
I. L.
,
Mickelsen
,
D.
,
Hansen
,
M. E.
,
Newlands
,
S. D.
,
Benoit
,
D. S. W.
,
Korshunov
,
V. A.
and
Ovitt
,
C. E.
(
2018
).
Localized Delivery of Amifostine Enhances Salivary Gland Radioprotection
.
J. Dent. Res.
97
,
1252
-
1259
.
Venkateswaran
,
K.
,
Shrivastava
,
A.
,
Agrawala
,
P. K.
,
Prasad
,
A. K.
,
Devi
,
S. C.
,
Manda
,
K.
,
Parmar
,
V. S.
and
Dwarakanath
,
B. S.
(
2019
).
Mitigation of radiation-induced gastro-intestinal injury by the polyphenolic acetate 7, 8–diacetoxy-4-methylthiocoumarin in mice
.
Sci. Rep.
9
,
14134
.
Verhaegen
,
F.
,
Butterworth
,
K. T.
,
Chalmers
,
A. J.
,
Coppes
,
R. P.
,
de Ruysscher
,
D.
,
Dobiasch
,
S.
,
Fenwick
,
J. D.
,
Granton
,
P. V.
,
Heijmans
,
S. H. J.
,
Hill
,
M. A.
et al. 
(
2023
).
Roadmap for precision preclinical x-ray radiation studies
.
Phys. Med. Biol.
68
,
06RM01
.
Westphal
,
D.
,
Kluck
,
R. M.
and
Dewson
,
G.
(
2014
).
Building blocks of the apoptotic pore: how Bax and Bak are activated and oligomerize during apoptosis
.
Cell Death Differ.
21
,
196
-
205
.
Wijers
,
O. B.
,
Levendag
,
P. C.
,
Braaksma
,
M. M. J.
,
Boonzaaijer
,
M.
,
Visch
,
L. L.
and
Schmitz
,
P. I. M.
(
2002
).
Patients with head and neck cancer cured by radiation therapy: a survey of the dry mouth syndrome in long-term survivors
.
Head Neck
24
,
737
-
747
.
Winter
,
C.
,
Keimel
,
R.
,
Gugatschka
,
M.
,
Kolb
,
D.
,
Leitinger
,
G.
and
Roblegg
,
E.
(
2021
).
Investigation of changes in saliva in radiotherapy-induced head neck cancer patients
.
Int. J. Environ. Res. Public Health
18
,
1629
.
Wittenborn
,
T. R.
,
Fahlquist Hagert
,
C.
,
Ferapontov
,
A.
,
Fonager
,
S.
,
Jensen
,
L.
,
Winther
,
G.
and
Degn
,
S. E.
(
2021
).
Comparison of gamma and x-ray irradiation for myeloablation and establishment of normal and autoimmune syngeneic bone marrow chimeras
.
PLoS ONE
16
,
e0247501
.
Woods
,
L. T.
,
Camden
,
J. M.
,
El-Sayed
,
F. G.
,
Khalafalla
,
M. G.
,
Petris
,
M. J.
,
Erb
,
L.
and
Weisman
,
G. A.
(
2015
).
Increased expression of TGF-β signaling components in a mouse model of fibrosis induced by submandibular gland duct ligation
.
PLoS One
10
,
e0123641
.

Competing interests

The authors declare no competing or financial interests.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution and reproduction in any medium provided that the original work is properly attributed.

Supplementary information