Regenerative medicine holds promise to cure radiation-induced salivary hypofunction, a chronic side effect in patients with head and neck cancers, therefore reliable preclinical models for salivary regenerative outcome will promote progress towards therapies. In this study, our objective was to develop a cone beam computed tomography-guided precision ionizing radiation-induced preclinical model of chronic hyposalivation using immunodeficient NSGSGM3 mice. Using a Schirmer's test based sialagogue-stimulated saliva flow kinetic measurement method, we demonstrated significant differences in hyposalivation specific to age, sex, precision-radiation dose over a chronic (6 months) timeline. NSG-SMG3 mice tolerated doses from 2.5 Gy up to 7.5 Gy. Interestingly, 5–7.5 Gy had similar effects on stimulated-saliva flow (∼50% reduction in young female at 6 months after precision irradiation over sham-treated controls), however, >5 Gy led to chronic alopecia. Different groups demonstrated characteristic saliva fluctuations early on, but after 5 months all groups nearly stabilized stimulated-saliva flow with low-inter-mouse variation within each group. Further characterization revealed precision-radiation-induced glandular shrinkage, hypocellularization, gland-specific loss of functional acinar and glandular cells in all major salivary glands replicating features of human salivary hypofunction. This model will aid investigation of human cell-based salivary regenerative therapies.

Head and neck cancer accounts for ∼4% of all cancers in the United States (1 ). The current standard of care for management can involve surgery, radiation therapy, chemotherapy, targeted therapy or a combination of these treatment modalities depending on various factors including the stages of cancer, the surgical accessibility of the tumor, and morbidity associated with each modality (2 ). In nearly 70% of the patients undergoing radiation therapy due to location of many of these oral cancers, non-diseased adjacent tissues such as salivary parenchyma is irreparably damaged due to its high sensitivity to radiation. This often results in radiation-induced “objective” hyposalivation or a “subjective” clinical condition referred to as xerostomia or dry mouth (2, 3 ). Complications of long-term radiation-induced hyposalivation include impaired taste, difficulties in speech, mastication and denture retention, increased risk of dental caries and an overall compromised quality-of-life (4). Management of hyposalivation often involves the use of cholinergic drugs such as pilocarpine (5 ) or cevimeline (6 ). The effectiveness of these systemic sialagogues often depends on the presence of functionally intact glandular tissue units (5 ). Recently, putative stem cells, that can restore the salivary gland functions upon orthotopic transplantation in radio-ablated salivary glands of recipient mice, have been reported in the submandibular glands of mice (710 ) and humans (11 ). To compare the quality and quantity of regenerative outcomes in transplant experiments, a standardized in vivo assay needs to be developed using preclinical animal models. However, there are several limitations with currently available preclinical models of hyposalivation. For example, the radiation dose used for many of these preclinical models range widely from 2 Gy to 30 Gy (12 ). Although there appears to be some consensus on 15 Gy for immunocompetent C57BL/6 mice (7–9, 13–15), but the radiation doses are not well characterized for immunodeficient and radiation-sensitive models such as NSG mice that are more relevant to xenotransplantation studies.

Furthermore, several methods of saliva collection from preclinical models after cholinergic stimulation have been reported including direct cannulation of excretory salivary ducts (16 ), vacuum pump suction of saliva from oral cavity (17 ), glass capillaries (18 ), micropipette (1921 ), gravimetry using pre-weighed filter paper strips (22 ) or a swab method (23 ). However, the methods such as direct cannulation are associated with risk of injury precluding repeated sampling of the gland, while gravimetry, glass capillary and micropipette-based methods show inter-operator variability (23 ). Therefore, an instantaneous method to measure kinetics of stimulated reflex saliva reliably and repetitively with little inter-operator variability would be advantageous.

The influence of age on the functionality of salivary glands is often debatable. Some studies have highlighted progressive and minor decline in saliva production with increasing age (2426 ), while others show aging does not have an impact on functionality of salivary glands in healthy individuals (27, 28 ). However, gender and hormone status in women are known to have an impact on salivary flow rates (2931 ). Nevertheless, age and gender specific differences in the context of stimulated salivary flow rates are not well understood in immunodeficient mouse models.

A small population c-KIT receptor expressing cells were identified in human salivary glands as putative stem cells and their transplantation restored salivary gland functions in NSG mice (11 ). Since the NSG mice do not express hc-KITL and therefore may not fully support the engraftment of hc-KIT expressing putative stem cells.

To address these lacunae, we comprehensively characterized a transgenic line of highly immunodeficient NSG-SGM3 mice that produces hc-KITL and identified that 20-week-old female NSG-SGM3 mice serves as a novel and robust model for chronic hyposalivation studies and preclinical assessment of experimental therapies. Further, we validated the Schirmer's test as a simple and reliable method to assess kinetics of saliva flow in mice. Moreover, characterization of glands at necropsy revealed highly complex differences in sex, age, radiation dose and salivary gland subtype-specific effects of salivary-ablative precision irradiation in NSG-SGM3 mice.

Animals

Female and male NOD.Cg-Prkdcscid Il2rgtm1Wjl Tg(CMV-IL3,CSF2, KITLG)1Eav/MloySzJ (NSG-SGM3) mice were housed in barrier facility and fed ad libitum with food pellets and acidified water. A total of 59 mice i.e., 40 young mice (∼20 weeks old, equal sex representation) and 19 old mice (∼1 year old, female only), were used in this pilot study. One male mouse died due to serious injury incurred during cage fighting at approximately 3 months postirradiation and data was censored. All procedures were approved by Mayo Clinic's Institutional Animal Care and Use Committee.

Salivary-Ablative Precision-Irradiation Set-up

Each mouse was anesthetized via isoflurane for the duration of the procedure via a nose cone delivery system on a bite block and placed in a feet-first prone (FFP) position in the irradiation chamber. A commercially available stereotactic stage (Model 900M, Kopf Instruments, Tujunga, CA) was modified to facilitate mounting to the X-RAD small animal irradiator system (SmART) (Precision X-ray Inc., North Bradford, CT) (32 ) (Fig. 1A). Opposing beam stereotactic radiation doses of a total 0 (or sham), 2.5, 5 or 7.5 Gy was delivered at 5 Gy/min dose rate using X-RAD SmART irradiator to the head and neck in a brain-sparing manner. The stereotactic X-RAD SmART irradiator set-up doesn't require additional blocking shield. Each mouse was volumetrically imaged with cone beam computed tomography (CBCT) using aluminum filter (60 kVp, 0.3 mA, 2.0-mm Al filter, 256 projections 0.2 mm 3-voxel size) and target was defined to ensure complete coverage of major salivary glands. Subsequently, the filter was switched to high-energy copper (irradiation) filter and beam is focused on the target using 10-mm-circular collimator for brain-sparing salivary-ablative radiation treatment.

Kinetics of Sialagogue-Stimulated Saliva Flow Measurements

Sham-treated and precision-irradiated animals were anesthetized using 2% isoflurane and injected subcutaneously with 2 mg/kg pilocarpine, a sialagogue drug. The reflex saliva was collected every 5 min for up to 30 min from the floor of the mouth using pre-weighed Schirmer's test strip (Tearflo, HUB Pharmaceuticals LLC, Rancho Cucamonga, CA) or filter paper. The distance of saliva migration under the capillary action (mm) was recorded on Schirmer's test strips. Subsequently, the weight of blotted strips was determined gravimetrically. Saliva was measured at three-week intervals in all groups for up to 6 months after precision irradiation.

Histology and Immunostaining

Whole salivary glands were fixed in 10% buffered formalin phosphate (Fisher Scientific; SF-1004), paraffin-embedded, sectioned (5-µm thick) and stained for Hematoxylin and Eosin (H&E), trichrome or Sirius red. For immunostaining, tissue sections were deparaffinized and antigen retrieval was performed in antigen unmasking solution (Vector Laboratories Inc, Burlingame, CA) following the manufacturer's instructions. Sectioned tissues were then stained with primary and secondary antibodies described in Supplementary Table S2 (https://doi.org/10.1667/RADE-21-00237.1.S1) Nuclei were counter-stained with DAPI Fluoromount-G® (Southern Biotech, Birmingham, AL). All stained slides were imaged with Cytation5 (BioTek Instruments Agilent Technologies).

Statistical Analysis

All analyses were performed using GraphPad Prism 8.0. For parametric comparisons between sham-treated and precision-irradiated groups, an unpaired two-tailed Student's t-test was performed. Multiple group comparisons were performed by one-way or two-way ANOVA. Differences were considered significant at P < 0.05. Data is displayed as mean ± SEM unless otherwise stated.

To generate a preclinical model of chronic hyposalivation, we utilized 3D image-guided stereotactic X-RAD SmART to radioablate the major salivary glands of immunodeficient NSG-SGM3 mice (33 ). Details of the functionality of this irradiator are described elsewhere (3436 ). Dosimetric calculations were performed with an in-house 1D “point-dose calculator” (PDC) tool developed by our medical physics group to enable efficient dose calculations (single prescription reference point on the central beam axis) (37 ). This simple 1D PDC was verified using the vendor supplied 3D Monte Carlo treatment planning system, SmART-ATP version 1.1 (SmART Scientific Solutions B.V., Maastricht, Netherlands), which is based on the open-source Monte Carlo code (EGSnrc/BEAMnrc)(38, 39 ) (Fig. 1AC). The chosen beam arrangement was parallel opposed beams using a 10-mm-circular collimator, with tipped laterals (89° and 271°) aligned to match beam divergence along the brain edge to maximize brain sparing and achieve specific targeting of major salivary glands along the midline of the mouse in the neck region [Fig. 1C (i-iii)]. The individual young (approximately 20-week-old male or female) or old (more than 1-year-old female) mouse was restrained on a bite block using an anesthesia nose cone and imaged with CBCT. Subsequently, CBCT image-guidance defined target and single radiation dose ranging from 0 (sham) to 10 Gy (old mice) or 7.5 Gy (young mice) was delivered at target site, i.e., the three major salivary glands (Fig. 1B). The high accuracy of CBCT guided targeting in each mouse and the beam arrangement enabled precision dose delivery to whole-salivary-gland structures. Both sham-treated and precision-radiation treated mice were followed for 6 months.

After evaluating all strains of mice and irradiated anatomical sites, we found many differences on the used radiation dose, sex, age and/or duration of the study period (Supplementary Table S1; https://doi.org/10.1667/RADE-21-00237.1.S1). For example, a single 15 Gy radiation dose on immunocompetent female C57BL/6 mice reduced ∼50% of saliva levels after 3 months postirradiation (40 ), and showed inflammatory lesions characteristic of human gland pathology (41 ). In contrast, Pringle et al. (11 ) showed that local irradiation with only a single 5 Gy dose was needed to achieve the same 50% reduction in saliva in immunodeficient NSG mice, implying that radiosensitivity between immunocompetent and immunodeficient mice is different. It is also important to note that a 5 Gy whole-body dose of radiation is lethal to NSG mice. While immunocompetent animals are informative and desirable models for studying damage to radio-ablated salivary glands, there is still an urgent need for an immunodeficient model that optimally supports the survival of human cells as part of regenerative therapies. In this study, we characterized the effect of precision radiation in a preclinical immunodeficient NSG-SGM3 model. Humanized NSG-SGM3 mice (33 ) are one of the most immunodeficient mouse strains. They were originally derived by crossing NOD/LtSz-scid IL2RG (NSG) mice (42 ) with NOD/LtSz-scid/scid (NOD/ SCID) mice, are engineered to express three human cytokines (43 ), and are extremely radiosensitive owing to the presence of the Prkdcscid mutation (44 ). Using our CBCT precision dose of radiation, we found that NSG-SGM3 mice (young and old) tolerated doses up to 7.5 Gy (Fig. 1DF). Increased mortality at 10 Gy in older mice precluded further life-style studies using this dose in NSG-SGM3 mice. We also observed that doses above 5 Gy led to permanent alopecia in the mandibular region, whereas 5 Gy itself showed transient alopecia. The dose of 2.5 Gy had no effect on hair follicles (Fig. 1G).

Next, we evaluated differences in saliva flow among the various groups. Saliva flow stimulation is a nerve-mediated reflex action modulated by the central nervous system, and this activity is an important indicator of functional salivary glands (45 ). In hyposalivation models, cholinergic drugs such as pilocarpine or carbachol are used to stimulate saliva secretion (11 , 46, 47 ). Although there are several methods that have been reported to collect stimulated reflex saliva and evaluate changes in flowrate in vivo or ex vivo (Supplementary Table S1; https://doi.org/10.1667/RADE21-00237.1.S1), we found that the field lacked an independent, easy and objective method for instantaneous kinetic measurement of stimulated-saliva flow in model organisms. Therefore, we adopted the Schirmer's test strip method to assess quantitative changes in salivary flowrate in NSG-SGM3 mice. Schirmer's test strips are routinely used in clinical practice to measure eye dryness (48 ) and hyposalivation in patients (49 ). Thus, we tested this method against the commonly used gravimetric method in our mice (Supplementary Table S1; https://doi.org/10.1667/RADE21-00237.1.S1). We injected 2 mg/kg of pilocarpine subcutaneously into the dorsal flanks of the animal and started continuously collecting saliva 5 min after injection from the floor of mouth using pre-weighed filter paper or Schirmer's test strips for 30 min while changing the strips every 5 min. Saliva migration under the capillary action (mm) was recorded on Schirmer's test strips by making a permanent pencil mark, and weight of blotted filter paper strips was determined gravimetrically (Fig. 2A). Measurements of saliva fraction over the chronic timeline demonstrated a high degree of correlation between the Schirmer's test strip method and the conventionally used gravimetric method in both young male (R2 = 0.90, P < 0.0001; Fig. 2B) and female (R2 = 0.88, P < 0.0001; Fig. 2C) NSG-SGM3 mice. Thus, we continued with the Schirmer's test for further analysis. The cumulative percent distribution analysis demonstrated a significant difference in saliva flowrate between young- and old-female mice (Fig. 2D), but not between young male and female mice (Fig. 2E), suggesting age-dependent changes in salivary physiology. Indeed, our kinetic analysis of flowrate in sham-treated mice suggested that reflex-saliva flowrate peaked a few minutes earlier in old females compared to young females, but no such difference was noted in young males vs. females after stimulation (Fig. 2F and G). In all mice, regardless of sex and age, the saliva flowrate declined after peaking approximately 10 min, but the decline was more dramatic in old-female mice (Fig. 2FI). Schirmer's test demonstrated age-specific (Fig. 2F) and sex-specific (Fig. 2G) differences in kinetics of stimulated-saliva flow captured over 30 min, as well as characteristic differences in flowrate at 10- and 25-min after stimulation (Fig. 2H and I). This trend was consistent throughout the chronic timeline in sham-treated mice (Fig. 3A and B). Taken together, our results suggest that Schirmer's test strips offer a rapid and reliable quantitative method to assess real-time changes in saliva flowrate in mice. Additionally, sexual dimorphism, as previously reported in submandibular glands due to the unique presence of granular convoluted tubular (GCT) structures in male mice (50 ), was also noticeable in NSG-SGM3 mice (Supplementary Fig. S1; https://doi.org/10.1667/RADE-21-00237.1.S1).

Then, we compared the kinetics of pilocarpine-stimulated saliva flow over the course of 30 min using Schirmer's test strips in young-male and -female mice, and old-female mice precision irradiated with 0, 2.5, 5 and 7.5 Gy doses. At 6 months postirradiation, we observed significantly reduced salivary flowrates at 10 min in mice precision irradiated with doses of 5 and 7.5 Gy compared to sham-treated controls (Fig. 3A and B). At 5 min after stimulation, saliva flowrate after 2.5 Gy precision irradiation was significantly lower in young males and old females, but not in young females compared to sham-treated mice. However, this difference disappeared at 10 min in males (Fig. 3A and B), suggesting important sex- and age-specific complex salivary physiology linked with different doses of precision radiation.

To define the chronic hyposalivation timeline in NSG-SGM3 model, we determined changes in stimulated-saliva fraction in precision-irradiated mice compared to sham-treated mice over 6 months postirradiation. We observed marked reduction in stimulated-saliva fraction 3 weeks after 2.5 Gy precision dose in young males and old females, but not in young females when compared to sham-irradiated animals (Fig. 3C). In young mice treated with 5 and 7.5 Gy, after about 12 weeks of fluctuations, saliva fraction stabilized at 50% compared to sham-treated animals, and this was sustained until the endpoint of 6 months (Fig. 3C and D). Interestingly, 2.5 Gy precision-irradiated mice showed a decline in stimulated-saliva fraction from 3 (males and old females) to 9 weeks (young females) after precision irradiation, but then showed partial restoration (Fig. 3C and D). In the case of old females, the stimulated-saliva fraction started marginally declining again after approximately 9 weeks (Fig. 3C and D). Of note, the sex-specific early responses of 2.5 Gy doses were not noticeable when the secreted fraction of saliva data from young males and females were combined (Supplementary Figure S2; https://doi.org/10.1667/RADE-21-00237.1.S1).

Radiation-induced salivary gland injury is also associated with glandular shrinkage, loss of acinar cell area, inflammation, fibrosis, microvascular injury and atrophy (51 ). Therefore, at the endpoint of 6 months, we sacrificed precision-irradiated as well as sham-treated mice and harvested major salivary glands. In general, the weight of major salivary glands in old-female mice was significantly greater relative to their body weight than young mice (Fig. 4A). The weight of major salivary glands relative to their body weight reduced significantly in young- and old-female mice with increasing precision-radiation dose (Fig. 4B and C). Interestingly, we did not observe such correlation between precision-radiation dose and salivary gland weight in male mice across the radiation doses (Fig. 4B). Next, we prepared tissue sections from harvested parotid (PG), submandibular (SMG) and sublingual (SLG) glands of sham-treated old- and young-female mice or 7.5 Gy precision irradiated and performed histological and immunohistochemical analysis. H&E staining revealed gross changes in tissue morphology and loss of acinar and ductal cells in all major salivary glands of precision-irradiated mice compared to sham-treated animals [Fig. 4D (i); Supplementary Fig. 3-(i); https://doi.org/10.1667/RADE-21-00237.1.S1]. Similarly, Sirius red [Fig. 4D (ii); Supplementary Fig. S3 (ii)] and trichrome [Fig. 4D (iii); Supplementary Fig. S3 (iii)] staining showed relative thickening of collagen fibers and fibrosis in tissue sections from precision-irradiated mice compared to sham-treated animals. Moreover, to assess the extent of radiation-induced damage to salivary epithelium, we stained whole salivary tissue sections for epithelial cell adhesion molecule (EpCAM) by immunostaining (10 ). We observed reduced expression of EpCAM in salivary epithelium of precision-irradiated mice compared to sham controls, which could be attributed to the loss of functional ducts and acini [Figure 4D (iv); Supplementary Figure S3 (iv)]. To validate this observation, we performed immunostaining for aquaporin-5, a membrane protein that facilitates water movement across the basal/lateral/apical membrane in acinar cells (52 ). Consistently, our results revealed reduction in expression of aquaporin-5 in acinar cells of precision-irradiated mice compared to sham controls [Fig. 4D (v); Supplementary Fig. S3 (v)]. Interestingly, we observed higher expression of β-catenin in ductal cells of precision-irradiated mice relative to controls [Fig. 4D (v); Supplementary Figure s3 (v)], suggesting ongoing tissue remodeling activity (10, 53, 54). In addition, we performed immunostaining of sodium-potassium-chloride cotransporter 1 (NKCC1), which is highly expressed in baso-lateral membrane of acinar cells. We observed reduced expression of NKCC1 in SMG of precision-irradiated mice when compared to controls [Fig. 4D (vi); Supplementary Fig. s3 (vi)]. However, we did not observe any changes in NKCC1 expression in PG and SLG between precision-irradiated and control mice.

Here, we report a novel and robust immunodeficient NSG-SGM3 mouse model of precision-radiation-induced chronic hyposalivation suitable for assessment of safety and efficacy of experimental human salivary cell therapies.

The precision dose of 5 Gy used in NSG-SGM3 mice not only selectively ablated salivary function in the long term, but also served to partially eliminate endogenous salivary cells to make room for cell engraftment. Interestingly, doses higher than 5 Gy did not guarantee a more complete salivary ablation, and only increased other complications, such as chronic alopecia or morbidity. Our kinetic analysis revealed that a 5 Gy dose was optimal for a long-term ∼50% reduction in saliva of NSG-SGM3 mice, consistent with observations reported elsewhere with NSG mice (11). Additionally, the effect of sex, age and radiation doses on reflex saliva kinetics had not been previously reported in NSG mice. Aging is a natural physiological phenomenon that has functional impact on most organs including salivary glands. The aging salivary glands of mice are characterized by acinar cell atrophy (26), apoptosis and loss of salivary gland function (55 ). Consistently, in our kinetic analysis we demonstrate that the reflex-saliva flowrate peaks a few minutes earlier in old females, compared to young females, but drops sharply to levels significantly lower than that of young female mice. Therefore, the older female mice may not be a preferred choice to measure human salivary regeneration and saliva production.

We also compared precision dose specific reflex saliva flowrate by gravimetric and Schrimer's strip methods over a period of 6 months. The two methods highly correlated with each other in both sexes. Each Schrimer's strip costs approximately $0.30 (U.S.). Single kinetic measurement of a flowrate over 30 min in a single animal requires 6 strips. The reading can be noted on-site in real time, or a pencil mark is sufficient to store and analyze the strips later with no concerns of evaporation of saliva.

SMG gland is the largest of three major salivary glands in rodents and is relatively more radiosensitive than PG (50, 56, 57 ), and therefore precision-radiation induced net loss in saliva may be due to higher loss of SMG function compared to PG or SLG. This is supported by more structural abnormalities in the SMG observed in our histological analysis. However, it is to be noted that this contrasts with human salivary glands, where the PG is the largest and most radiosensitive than the SMG and SLG (51 ). In fact, SMG contributes to the majority of produced saliva in rodents and irradiation of salivary glands may cause imbalance in overall salivary composition. However, these gland-specific dynamics did not affect our overall salivary readout. Sexual dimorphisms in SMG in mice are also well documented and differ with humans. There are larger GCT ducts in male SMG compared to female SMG, and such sex specific differences do not exist in human SMG (5860 ). In fact, the SMG of female mouse morphologically resembles to that of human SMG due to their lack of GCT.

We also noted that after 5 and 7.5 Gy precision doses of radiation, reflex saliva levels spontaneously, but transiently, improved before it declined again. Interestingly, the onset, duration and level of this transient host response appeared to vary based on age, sex and radiation dose.

Currently, the models available for human salivary stem cell testing are limited (Supplementary Table S1; https://doi.org/10.1667/RADE-21-00237.1.S1). The strains with incomplete immunodeficiency due to residual NK cell activity after single genetic alteration include Nude (61 ), SCID (62, 63 ), Rag1–/– (64 ), and Rag2–/– mice (65 ), and those with multiple genetic alterations include non-obese diabetic (NOD)-SCID mice (66 ). SCID mice carry a spontaneous Prkdc mutation that affects V(D)J recombination during T- and B-cell maturation and NHEJ DNA repair in all cells rendering derivative mice like NSG and NSG-SGM3 sensitive to radiation (44, 67 ). Whereas Rag2–/– mice derived with a germline deletion of the recombination-activating gene 2 (Rag2), an enzyme specifically required for V(D)J recombination during T- and B-cell maturation enables the cells to tolerate “normal” murine related radiation doses compared to NSG mice (68 ). While highly immunodeficient mice lacking NK cells was made possible with loss-of-function genetic alteration of gamma chain/ Il2rg gene (6971 ) in NSG and Rag2/Il2rg (R2G2) double knockout (68, 72 ), only NSG mice model, to our knowledge, has been reported for human salivary stem cell transplants (11). Future comparative investigations of the extent of host salivary dysfunction after radiation injury and efficiency of human salivary stem cell engraftment in Il2rg–/– mice with SCID or Rag2–/– are needed.

The NSG-SGM3 mice described here has the following features; first, it expresses physiological levels of 3 human cytokines (hIL-3, hGM-CSF and hc-KITL) and has been previously shown to promote better engraftment of hc-KIT receptor expressing hematopoietic stem cells than NSG mice (73, 74 ). Second, NSG-SGM3 mice has all the attributes of NSG mice including radiosensitivity and immunodeficiency intact. However, it is important to underline that the expression of hc-KIT receptor is not exclusive to a sub-population of salivary gland cells but is widely expressed by other tissue cells including hematopoietic cells, germ cells, mast cells, melanoma cells, gastrointestinal cajal cells, breast luminal cells and small neurons of brain (75, 76 ). Therefore, our planned future studies involve comparison of salivary transplants in NSG and NSG-SGM3 mice to delineate the functionality of hc-KIT in salivary stem cell engraftment. Indeed, hc-KIT as salivary stem cell marker remains disputed and the number of c-Kit expressing cells was extremely low and decreases with age in salisphere cultures (11 ). Further research into the phenotype to prospectively isolate salivary stem cells is urgently needed.

The use of isoflurane anesthesia during saliva measurement is a limitation of the current study, as it has been reported to reduce pilocarpine-induced saliva production (77 ). In fact, most other methods of saliva measurement use ketamine anesthesia that is known to increase salivary secretions through sympathetic stimulation (78 ). However, since the isoflurane is used similarly for sham-treated and radiation-treated groups of mice in this study, any influence of anesthesia on saliva secretion is normalized. Although, the method that involves the use of vacuum pump suction for saliva collection does not require anesthesia (17 ), but in the absence of anesthesia animals tend to struggle and this may induce varying level of anxiety that in turn may have an influence on secreted saliva. In the current study, we utilized precision stereotactic doses to radio-ablate salivary glands in NSG-SGM3 young male, female and old female mice to model hyposalivation. The high accuracy of CBCT-guided targeting and the beam arrangement enabled precision dose delivery to whole salivary gland structures in each mouse in a brain sparing manner. This method may overcome issues associated with head and neck radiotherapy which may lead to cognitive impairment (79 ).

In conclusion, we propose 5 Gy precision-irradiated female immunodeficient NSG-SGM3 mice as a suitable preclinical model for chronic hyposalivation studies, and the Schirmer's test as a reliable quantitative assay for salivary regenerative outcome. We expect the new preclinical model to facilitate rapid progress towards development of salivary regenerative therapies.

Supplementary Fig. S1: H&E staining of the submandibular gland of male and female NSG-SGM3 mice showing unique presence of GCT structures indicated by yellow arrowheads. Male GCT is markedly hypertrophic composed of number of large secretory granules compared to females. Scale bars = 200 µm. H&E staining showed absence of GCT structures in male and female human submandibular glands. Scale bars = 20 µm.

Supplementary Fig. S2: Plot showing data from young NSG-SGM3 mice (male and female combined) from FIG. 3C. Data is plotted as mean±SEM.

Supplementary Fig. S3: Necropsy reveals precision-IR-associated changes on salivary tissue, stroma, epithelial cell polarity, and integrity in old female NSG-SGM3 mice. Representative photomicrograph of stained parotid, submandibular and sublingual glands of old female (n=3) NSG-SGM3 mice at 6 month post-IR (0 and 7.5 Gy). Sections were stained with H&E (i), Sirius red (ii), trichrome (iii), and immunostained for EpCAM (iv), Aquaporin 5/beta-catenin (v) and NKCC1 (iv). Scale bars = 20µm.

Supplementary Table S1: Summary of methods reported for modelling radiation induced hyposalivation in animals.

Supplementary Table S2: Details of antibodies used for immunohistochemistry.

This study was supported by Mayo Clinic and in part by grants to JJ and NK from Mayo Clinic Center for Regenerative Medicine. We thank Dr. James E. Melvin (NIH/NIDCR) for helpful discussions on salivary measurement methods, and Drs. Geng Xian Shi and Gang Liu for their technical assistance. All procedures were approved by Mayo Clinic Institutional Animal Care and Use Committee and Institutional Review Board.

1.
Siegel
RL,
Miller
KD,
Jemal
A.
Cancer Statistics, 2017.
CA: a cancer journal for clinicians.
2017
;
67
(
1
):
7
30
.
2.
Vissink
A,
Burlage
FR,
Spijkervet
FK,
Jansma
J,
Coppes
RP.
Prevention and treatment of the consequences of head and neck radiotherapy.
Critical reviews in oral biology and medicine: an official publication of the American Association of Oral Biologists.
2003
;
14
(
3
):
213
25
.
3.
Sun
HB,
Gao
XJ,
Deng
J,
Li
NY,
Lu
HJ.
Progress of oral sequelae during head-neck radiotherapy.
The Chinese journal of dental research: the official journal of the Scientific Section of the Chinese Stomatological Association.
2010
;
13
(
1
):
51
5
.
4.
Villa
A,
Connell
CL,
Abati
S.
Diagnosis and management of xerostomia and hyposalivation.
Ther Clin Risk Manag.
2015
;
11
:
45
51
.
5.
Takakura
AC,
Moreira
TS,
Laitano
SC,
De Luca Junior
LA,
Renzi
A,
Menani
JV.
Central muscarinic receptors signal pilocarpine-induced salivation.
J Dent Res.
2003
;
82
(
12
):
993
7
.
6.
Chambers
MS,
Jones
CU,
Biel
MA,
Weber
RS,
Hodge
KM,
Chen
Y,
et al.
Open-label, long-term safety study of cevimeline in the treatment of postirradiation xerostomia.
Int J Radiat Oncol Biol Phys.
2007
;
69
(
5
):
1369
76
.
7.
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
(
4
):
e2063
.
8.
Mishima
K,
Inoue
H,
Nishiyama
T,
Mabuchi
Y,
Amano
Y,
Ide
F,
et al.
Transplantation of side population cells restores the function of damaged exocrine glands through clusterin.
Stem Cells.
2012
;
30
(
9
):
1925
37
.
9.
Nanduri
LS,
Lombaert
IM,
van der Zwaag
M,
Faber
H,
Brunsting
JF,
van Os
RP,
et al.
Salisphere derived c-Kit+ cell transplantation restores tissue homeostasis in irradiated salivary gland.
Radiother Oncol.
2013
;
108
(
3
):
458
63
.
10.
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
(
1
):
150
62
.
11.
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
(
3
):
640
52
.
12.
Weng
PL,
Aure
MH,
Ovitt
CE.
Concise Review: A critical evaluation of criteria used to define salivary gland stem cells.
Stem Cells.
2019
;
37
(
9
):
1144
50
.
13.
Varghese
JJ,
Schmale
IL,
Mickelsen
D,
Hansen
ME,
Newlands
SD,
Benoit
DSW,
et al.
Localized delivery of amifostine enhances salivary gland radioprotection.
J Dent Res.
2018
;
97
(
11
):
1252
9
.
14.
Zeidan
YH,
Xiao
N,
Cao
H,
Kong
C,
Le
QT,
Sirjani
D.
Botulinum toxin confers radioprotection in murine salivary glands.
Int J Radiat Oncol Biol Phys.
2016
;
94
(
5
):
1190
7
.
15.
Yamamura
Y,
Yamada
H,
Sakurai
T,
Ide
F,
Inoue
H,
Muramatsu
T,
et al.
Treatment of salivary gland hypofunction by transplantation with dental pulp cells.
Arch Oral Biol.
2013
;
58
(
8
):
935
42
.
16.
Marmary
Y,
Fox
PC,
Baum
BJ.
Fluid secretion rates from mouse and rat parotid glands are markedly different following pilocarpine stimulation.
Comp Biochem Physiol A Comp Physiol.
1987
;
88
(
2
):
307
10
.
17.
Lin
AL,
Johnson
DA,
Wu
Y,
Wong
G,
Ebersole
JL,
Yeh
CK.
Measuring short-term gamma-irradiation effects on mouse salivary gland function using a new saliva collection device.
Arch Oral Biol.
2001
;
46
(
11
):
1085
9
.
18.
Scofield
RH,
Asfa
S,
Obeso
D,
Jonsson
R,
Kurien
BT.
Immunization with short peptides from the 60-kDa Ro antigen recapitulates the serological and pathological findings as well as the salivary gland dysfunction of Sjogren's syndrome.
J Immunol.
2005
;
175
(
12
):
8409
14
.
19.
Deshmukh
US,
Nandula
SR,
Thimmalapura
PR,
Scindia
YM,
Bagavant
H.
Activation of innate immune responses through Toll-like receptor 3 causes a rapid loss of salivary gland function.
J Oral Pathol Med.
2009
;
38
(
1
):
42
7
.
20.
Deshmukh
US,
Ohyama
Y,
Bagavant
H,
Guo
X,
Gaskin
F,
Fu
SM.
Inflammatory stimuli accelerate Sjogren's syndrome-like disease in (NZB x NZW)F1 mice.
Arthritis Rheum.
2008
;
58
(
5
):
1318
23
.
21.
Szczerba
BM,
Kaplonek
P,
Wolska
N,
Podsiadlowska
A,
Rybakowska
PD,
Dey
P,
et al.
Interaction between innate immunity and Ro52-induced antibody causes Sjogren's syndrome-like disorder in mice.
Ann Rheum Dis.
2016
;
75
(
3
):
617
22
.
22.
Boyd
A,
Aragon
IV,
Abou Saleh
L,
Southers
D,
Richter
W.
The cAMP-phosphodiesterase 4 (PDE4) controls beta-adrenoceptor- and CFTR-dependent saliva secretion in mice.
Biochem J.
2021
;
478
(
10
):
1891
906
.
23.
Bagavant
H,
Trzeciak
M,
Papinska
J,
Biswas
I,
Dunkleberger
ML,
Sosnowska
A,
et al.
A method for the measurement of salivary gland function in mice.
J Vis Exp.
2018
(
131
).
24.
Turner
MD,
Ship
JA.
Dry mouth and its effects on the oral health of elderly people.
J Am Dent Assoc.
2007
;
138
Suppl:
15S
20S
.
25.
Xu
F,
Laguna
L,
Sarkar
A.
Aging-related changes in quantity and quality of saliva: Where do we stand in our understanding?
J Texture Stud.
2019
;
50
(
1
):
27
35
.
26.
Enoki
N,
Kiyoshima
T,
Sakai
T,
Kobayashi
I,
Takahashi
K,
Terada
Y,
et al.
Age-dependent changes in cell proliferation and cell death in the periodontal tissue and the submandibular gland in mice: a comparison with other tissues and organs.
J Mol Histol.
2007
;
38
(
4
):
321
32
.
27.
Gupta
A,
Epstein
JB,
Sroussi
H.
Hyposalivation in elderly patients.
J Can Dent Assoc.
2006
;
72
(
9
):
841
6
.
28.
Lasisi
TJ,
Shittu
ST,
Oguntokun
MM,
Tiamiyu
NA.
Aging affects morphology but not stimulated secretion of saliva in rats.
Ann Ib Postgrad Med.
2014
;
12
(
2
):
109
14
.
29.
Inoue
H,
Ono
K,
Masuda
W,
Morimoto
Y,
Tanaka
T,
Yokota
M,
et al.
Gender difference in unstimulated whole saliva flow rate and salivary gland sizes.
Arch Oral Biol.
2006
;
51
(
12
):
1055
60
.
30.
Rukmini
JN,
Sachan
R,
Sibi
N,
Meghana
A,
Malar
CI.
Effect of Menopause on Saliva and Dental Health.
J Int Soc Prev Community Dent.
2018
;
8
(
6
):
529
33
.
31.
Minicucci
EM,
Pires
RB,
Vieira
RA,
Miot
HA,
Sposto
MR.
Assessing the impact of menopause on salivary flow and xerostomia.
Aust Dent J.
2013
;
58
(
2
):
230
4
.
32.
Walb
MC,
Carlson
BL,
Sarkaria
JN,
Tryggestad
EJ.
Quantifying the setup uncertainty of a stereotactic murine micro-image guided radiation therapy system.
Br J Radiol.
2019
;
92
(
1095
):
20180487
.
33.
Wunderlich
M,
Chou
FS,
Link
KA,
Mizukawa
B,
Perry
RL,
Carroll
M,
et al.
AML xenograft efficiency is significantly improved in NOD/SCID-IL2RG mice constitutively expressing human SCF, GM-CSF and IL-3.
Leukemia.
2010
;
24
(
10
):
1785
8
.
34.
Verhaegen
F,
Granton
P,
Tryggestad
E.
Small animal radiotherapy research platforms.
Phys Med Biol.
2011
;
56
(
12
):
R55
83
.
35.
Stewart
JM,
Ansell
S,
Lindsay
PE,
Jaffray
DA.
Online virtual isocenter based radiation field targeting for high performance small animal microirradiation.
Phys Med Biol.
2015
;
60
(
23
):
9031
46
.
36.
Clarkson
R,
Lindsay
PE,
Ansell
S,
Wilson
G,
Jelveh
S,
Hill
RP,
et al.
Characterization of image quality and image-guidance performance of a preclinical microirradiator.
Med Phys.
2011
;
38
(
2
):
845
56
.
37.
Smith
EJ,
Tryggestad
EJ,
Carlson
BL,
Walb
MC,
Sarkaria
JN.
Imaging and dosimetry study of inter-fraction setup error in a murine xenograft flank tumor radiation model.
Radiat Res.
2020
;
193
(
2
):
161
70
.
38.
van Hoof
SJ,
Granton
PV,
Verhaegen
F.
Development and validation of a treatment planning system for small animal radiotherapy: SmART-Plan.
Radiother Oncol.
2013
;
109
(
3
):
361
6
.
39.
Rogers
DW,
Faddegon
BA,
Ding
GX,
Ma
CM,
We
J,
Mackie
TR.
BEAM: a Monte Carlo code to simulate radiotherapy treatment units.
Med Phys.
1995
;
22
(
5
):
503
24
.
40.
Peng
X,
Wu
Y,
Brouwer
U,
van Vliet
T,
Wang
B,
Demaria
M,
et al.
Cellular senescence contributes to radiation-induced hyposalivation by affecting the stem/progenitor cell niche.
Cell Death Dis.
2020
;
11
(
10
):
854
.
41.
Urek
MM,
Bralic
M,
Tomac
J,
Borcic
J,
Uhac
I,
Glazar
I,
et al.
Early and late effects of X-irradiation on submandibular gland: a morphological study in mice.
Arch Med Res.
2005
;
36
(
4
):
339
43
.
42.
Shultz
LD,
Lyons
BL,
Burzenski
LM,
Gott
B,
Chen
X,
Chaleff
S,
et al.
Human lymphoid and myeloid cell development in NOD/ LtSz-scid IL2R gamma null mice engrafted with mobilized human hemopoietic stem cells.
J Immunol.
2005
;
174
(
10
):
6477
89
.
43.
Nicolini
FE,
Cashman
JD,
Hogge
DE,
Humphries
RK,
Eaves
CJ.
NOD/SCID mice engineered to express human IL-3, GM-CSF and Steel factor constitutively mobilize engrafted human progenitors and compromise human stem cell regeneration.
Leukemia.
2004
;
18
(
2
):
341
7
.
44.
Biedermann
KA,
Sun
JR,
Giaccia
AJ,
Tosto
LM,
Brown
JM.
scid mutation in mice confers hypersensitivity to ionizing radiation and a deficiency in DNA double-strand break repair.
Proc Natl Acad Sci U S A.
1991
;
88
(
4
):
1394
7
.
45.
Proctor
GB.
The physiology of salivary secretion.
Periodontol 2000.
2016
;
70
(
1
):
11
25
.
46.
Nanduri
LS,
Maimets
M,
Pringle
SA,
van der Zwaag
M,
van Os
RP,
Coppes
RP.
Regeneration of irradiated salivary glands with stem cell marker expressing cells.
Radiother Oncol.
2011
;
99
(
3
):
367
72
.
47.
Gilman
KE,
Camden
JM,
Klein
RR,
Zhang
Q,
Weisman
GA,
Limesand
KH.
P2X7 receptor deletion suppresses gamma-radiation-induced hyposalivation.
Am J Physiol Regul Integr Comp Physiol.
2019
;
316
(
5
):
R687
R96
.
48.
Senchyna
M,
Wax
MB.
Quantitative assessment of tear production: A review of methods and utility in dry eye drug discovery.
J Ocul Biol Dis Infor.
2008
;
1
(
1
):
1
6
.
49.
López-Jornet
P,
Camacho-Alonso
F,
Bermejo-Fenoll
A.
A simple test for salivary gland hypofunction using Oral Schirmer's test.
J Oral Pathol Med.
2006
;
35
(
4
):
244
8
.
50.
Maruyama
CL,
Monroe
MM,
Hunt
JP,
Buchmann
L,
Baker
OJ.
Comparing human and mouse salivary glands: A practice guide for salivary researchers.
Oral Dis.
2019
;
25
(
2
):
403
15
.
51.
Grundmann
O,
Mitchell
GC,
Limesand
KH.
Sensitivity of salivary glands to radiation: from animal models to therapies.
J Dent Res.
2009
;
88
(
10
):
894
903
.
52.
Hosoi
K,
Yao
C,
Hasegawa
T,
Yoshimura
H,
Akamatsu
T.
Dynamics of Salivary Gland AQP5 under Normal and Pathologic Conditions.
Int J Mol Sci.
2020
;
21
(
4
).
53.
Hai
B,
Yang
Z,
Millar
SE,
Choi
YS,
Taketo
MM,
Nagy
A,
et al.
Wnt/beta-catenin signaling regulates postnatal development and regeneration of the salivary gland.
Stem Cells Dev.
2010
;
19
(
11
):
1793
801
.
54.
Hai
B,
Yang
Z,
Shangguan
L,
Zhao
Y,
Boyer
A,
Liu
F.
Concurrent transient activation of Wnt/beta-catenin pathway prevents radiation damage to salivary glands.
Int J Radiat Oncol Biol Phys.
2012
;
83
(
1
):
e109
16
.
55.
Choi
JS,
Park
IS,
Kim
SK,
Lim
JY,
Kim
YM.
Analysis of age-related changes in the functional morphologies of salivary glands in mice.
Arch Oral Biol.
2013
;
58
(
11
):
1635
42
.
56.
Coppes
RP,
Vissink
A,
Konings
AW.
Comparison of radiosensitivity of rat parotid and submandibular glands after different radiation schedules.
Radiother Oncol.
2002
;
63
(
3
):
321
8
.
57.
Kim
JY,
An
CH,
Kim
JY,
Jung
JK.
experimental animal model systems for understanding salivary secretory disorders.
Int J Mol Sci.
2020
;
21
(
22
).
58.
Pinkstaff
CA.
Salivary gland sexual dimorphism: a brief review.
Eur J Morphol.
1998
;
36
Suppl:
31
4
.
59.
Gresik
EW.
The granular convoluted tubule (GCT) cell of rodent submandibular glands.
Microsc Res Tech.
1994
;
27
(
1
):
1
24
.
60.
Gresik
EW,
Hosoi
K,
Kurihara
K,
Maruyama
S,
Ueha
T.
The rodent granular convoluted tubule cell–an update.
Eur J Morphol.
1996
;
34
(
3
):
221
4
.
61.
Flanagan
SP.
'Nude', a new hairless gene with pleiotropic effects in the mouse.
Genet Res.
1966
;
8
(
3
):
295
309
.
62.
Shultz
LD,
Ishikawa
F,
Greiner
DL.
Humanized mice in translational biomedical research.
Nat Rev Immunol.
2007
;
7
(
2
):
118
30
.
63.
Bosma
GC,
Custer
RP,
Bosma
MJ.
A severe combined immunodeficiency mutation in the mouse.
Nature.
1983
;
301
(
5900
):
527
30
.
64.
Mombaerts
P,
Iacomini
J,
Johnson
RS,
Herrup
K,
Tonegawa
S,
Papaioannou
VE.
RAG-1-deficient mice have no mature B and T lymphocytes.
Cell.
1992
;
68
(
5
):
869
77
.
65.
Shinkai
Y,
Rathbun
G,
Lam
KP,
Oltz
EM,
Stewart
V,
Mendelsohn
M,
et al.
RAG-2-deficient mice lack mature lymphocytes owing to inability to initiate V(D)J rearrangement.
Cell.
1992
;
68
(
5
):
855
67
.
66.
Shultz
LD,
Schweitzer
PA,
Christianson
SW,
Gott
B,
Schweitzer
IB,
Tennent
B,
et al.
Multiple defects in innate and adaptive immunologic function in NOD/LtSz-scid mice.
J Immunol.
1995
;
154
(
1
):
180
91
.
67.
Fulop
GM,
Phillips
RA.
The scid mutation in mice causes a general defect in DNA repair.
Nature.
1990
;
347
(
6292
):
479
82
.
68.
Goldman
JP,
Blundell
MP,
Lopes
L,
Kinnon
C,
Di Santo
JP,
Thrasher
AJ.
Enhanced human cell engraftment in mice deficient in RAG2 and the common cytokine receptor gamma chain.
Br J Haematol.
1998
;
103
(
2
):
335
42
.
69.
Cao
X,
Shores
EW,
Hu-Li
J,
Anver
MR,
Kelsall
BL,
Russell
SM,
et al.
Defective lymphoid development in mice lacking expression of the common cytokine receptor gamma chain.
Immunity.
1995
;
2
(
3
):
223
38
.
70.
DiSanto
JP,
Muller
W,
Guy-Grand
D,
Fischer
A,
Rajewsky
K.
Lymphoid development in mice with a targeted deletion of the interleukin 2 receptor gamma chain.
Proc Natl Acad Sci U S A.
1995
;
92
(
2
):
377
81
.
71.
Ohbo
K,
Suda
T,
Hashiyama
M,
Mantani
A,
Ikebe
M,
Miyakawa
K,
et al.
Modulation of hematopoiesis in mice with a truncated mutant of the interleukin-2 receptor gamma chain.
Blood.
1996
;
87
(
3
):
956
67
.
72.
Mazurier
F,
Fontanellas
A,
Salesse
S,
Taine
L,
Landriau
S,
Moreau-Gaudry
F,
et al.
A novel immunodeficient mouse model–RAG2 x common cytokine receptor gamma chain double mutants–requiring exogenous cytokine administration for human hematopoietic stem cell engraftment.
J Interferon Cytokine Res.
1999
;
19
(
5
):
533
41
.
73.
Billerbeck
E,
Barry
WT,
Mu
K,
Dorner
M,
Rice
CM,
Ploss
A.
Development of human CD4+FoxP3+ regulatory T cells in human stem cell factor-, granulocyte-macrophage colony-stimulating factor-, and interleukin-3-expressing NOD-SCID IL2Rgamma(null) humanized mice.
Blood.
2011
;
117
(
11
):
3076
86
.
74.
Wunderlich
M,
Chou
FS,
Sexton
C,
Presicce
P,
Chougnet
CA,
Aliberti
J,
et al.
Improved multilineage human hematopoietic reconstitution and function in NSGS mice.
PLoS One.
2018
;
13
(
12
):
e0209034
.
75.
Lennartsson
J,
Ronnstrand
L.
Stem cell factor receptor/c-Kit: from basic science to clinical implications.
Physiol Rev.
2012
;
92
(
4
):
1619
49
.
76.
Liang
J,
Wu
YL,
Chen
BJ,
Zhang
W,
Tanaka
Y,
Sugiyama
H.
The C-kit receptor-mediated signal transduction and tumor-related diseases.
Int J Biol Sci.
2013
;
9
(
5
):
435
43
.
77.
Knudsen
J,
Nauntofte
B,
Josipovic
M,
Engelholm
SA,
Hyldegaard
O.
Effects of isoflurane anesthesia and pilocarpine on rat parotid saliva flow.
Radiat Res.
2011
;
176
(
1
):
84
8
.
78.
Kohrs
R,
Durieux
ME.
Ketamine: Teaching an old drug new tricks.
Anesth Analg.
1998
;
87
(
5
):
1186
93
.
79.
Zer
A,
Pond
GR,
Razak
ARA,
Tirona
K,
Gan
HK,
Chen
EX,
et al.
Association of neurocognitive deficits with radiotherapy or chemoradiotherapy for patients with head and neck cancer.
JAMA Otolaryngol Head Neck Surg.
2018
;
144
(
1
):
71
9
.

Author notes

Editor's note. The online version of this article (DOI: https://doi.org/10.1667/RADE-21-00237.1) contains supplementary information that is available to all authorized users.

1 These authors contributed equally.

2 Current affiliation: Department of Radiation Oncology and Winship Cancer Institute, Emory University, Atlanta, GA.

3 Senior authors.

Supplementary data