Transient adenovirus-Cre infection causes long-lasting remodeling of the mammary gland immune landscape.

Introduction
Development of cancers, such as breast cancer, from their corresponding cellular origins occurs within a dynamic immunological ecosystem in which immune cells can both restrain and promote tumorigenesis. A “fight” between anti-tumorigenic and pro-tumorigenic immune and inflammatory mechanisms can determine the course of tumor development and its phenotype1,2. In breast cancer, during its early stages of development, innate and adaptive immune surveillance mechanisms, such as cytotoxic T cells, NK cells, dendritic cells, and macrophages, can recognize and eliminate transformed mammary epithelial cells (MECs), limiting malignant progression3–5. However, as oncogenic signaling, genomic instability, and tissue remodeling intensify, the evolving tumor microenvironment can reshape these immune populations, such as during the transition from the precancer to cancer stages6. Myeloid-derived immune cells, particularly tumor-associated macrophages and neutrophils, may adopt pro-tumorigenic phenotypes that support immunosuppression, angiogenesis, and invasion, while regulatory T cells and exhausted T-cell states dampen effective antitumor immunity7–9. Defining how these immune populations are recruited, activated or suppressed, and reprogrammed during the early phases of mammary tumorigenesis is essential for designing strategies that bolster endogenous immune defenses. By understanding the immunological events that tip the balance from elimination to escape, immune-based preventive approaches may be developed to intercept breast cancer at its roots, before clinically detectable disease emerges.
As access to early stages of human cancer development, including that of breast cancer, is very difficult if not impossible, in order to achieve the goal of immune-based cancer interception, animal models that can recapitulate this early phase of tumorigenesis are essential. In breast cancer, several elegant approaches have been developed to model mammary tumor initiation from a defined subpopulation of MECs and to follow their progression in mice10. One of the prevalent approaches is the Cre/loxP-based strategy that enables precise temporal and spatial control over oncogenic events within the mammary epithelium11. By placing Cre recombinase under the control of mammary-specific promoters, such as MMTV, Wap, or Krt812,13, and activating it through ligand-dependent systems such as CreER (Cre-estrogen receptor fusion, inducible by tamoxifen14,15), researchers can trigger defined genetic alterations, including oncogene activation or tumor suppressor loss, at selected MEC subpopulations or developmental stages. This approach faithfully recapitulates the stepwise evolution of breast cancer, permitting the study of early cellular transformation, clonal expansion, and cancer-immune cell interaction within an intact immune microenvironment16–19. To complement the CreER/tamoxifen-based inducible approach, my group developed a novel approach based on intraductal injection of Cre-expressing adenovirus to mouse mammary glands20–22. Adenovirus is a DNA virus that does not integrate into the host genome and the adenoviral vector we use only leads to transient expression of the Cre recombinase; thus, it serves a similar purpose as an inducible Cre system. The cell type-specificity (for inducing Cre-mediated recombination) is achieved by using different MEC subpopulation-specific promoters (e.g., Krt8, Wap) to drive adenovirus-Cre expression.
Breast cancer mouse models based on both CreER/tamoxifen and adenovirus-Cre-based approaches offer excellent in vivo systems to characterize immune cell phenotypes and their dynamic changes at different stages of mammary tumorigenesis. In these induced mice, mammary gland immune cells can be modulated by signals from mutated MECs. However, it is unclear whether the injected tamoxifen or adeno-Cre virus (to induce Cre-mediated recombination and initiation of mammary tumorigenesis) can affect the immune microenvironment. This is possible as administration of tamoxifen in vivo can potentially affect estrogen signaling in the mammary gland23,24, while estrogen/ER (estrogen receptor) signaling plays a key role in regulating mammary gland immune cell populations25. Similarly, injection of adenovirus to mammary glands should trigger virus-induced immune reactions, at least transiently20. However, it is unclear whether any of these induction approach-related modulations of mammary gland immune cells is transient or more long-lasting and whether such change would affect immune cell phenotypes and dynamics, and ultimately, mammary tumor development. A thorough understanding of these should have important implications for the proper use of these modeling approaches to study immune cells during mammary tumorigenesis.

Results

Adenoviral infection increases leukocyte infiltration in the mammary gland
In this study, we used a breast cancer mouse model that we previously developed based on induced loss of p53 in luminal MECs17. The experimental mice (PY) carry Trp53 conditional knockout alleles (Trp53fl/fl)26 along with a conditional Cre-reporter [Rosa26-LSL-YFP (R26Y)]27, while R26Y reporter-only mice served as wild-type (WT) controls (Supplementary figure 1a). To induce mammary tumor initiation specifically in luminal MECs, we utilized an adenovirus expressing Cre under the control of the luminal keratin 8 (Krt8, or K8) promoter (Ad-K8-Cre)20. Intraductal injection of Ad-K8-Cre into the abdominal (#4) mammary glands of PY female mice induced simultaneous YFP reporter expression and Trp53 deletion in luminal MECs, resulting in the development of mammary tumors (predominantly of the Claudin-low subtype) with 100% penetrance17. For comparison, we also used a CreER/tamoxifen-based approach to induce p53-loss in the same luminal MEC population. We generated KPY experimental mice and KY controls by crossing the K8-CreER allele, a transgenic construct expressing CreERT2 under the same K8 promoter12, with PY and R26Y mice, respectively (Supplementary figure 1b). Tamoxifen (TAM) induction in KPY female mice led to development of mammary tumors similar to those observed in Ad-K8-Cre-induced PY mice, also with full penetrance17.
To investigate how p53-loss in luminal MECs alters the mammary gland immune microenvironment during the premalignant stage, we induced p53 deletion by administering Ad-K8-Cre or TAM to PY (and R26Y) or KPY (and KY) female mice at ~8 weeks of age, respectively. Based on our previous findings, mammary tumors begin to emerge in induced mice around 5–6 months after induction17. Accordingly, in this study we defined 3 and 5 months post-induction as the mid- and late-premalignant stages, respectively (Fig. 1a). To determine whether intraductally delivered adenovirus itself elicits additional immunogenicity capable of altering the mammary immune microenvironment, we evaluated leukocyte infiltration one week after injection. FACS analysis revealed that adenovirus-injected R26Y mice exhibited a significantly higher proportion of CD45+ leucocytes in the injected abdominal mammary glands (55.5%±7.1%) compared with non-injected thoracic glands from the same mice (35.76%±4.74%) (Fig. 1b–c). In contrast, the K8-CreER-based model relies on systemic TAM administration via intraperitoneal injection, which could in principle exert global effects on all mammary glands through its modulation of estrogen signaling. To assess this possibility, we compared leukocyte infiltration in mammary glands from TAM-treated KY mice with that of age-matched, untreated KY controls. The data indicate that TAM treatment does not significantly alter leukocyte infiltration (28.36%±3.16% vs. 30.98%±3.32%) (Fig. 1d–e).

Transient adenoviral infection in the mammary gland increases T cell infiltration
To determine whether transient adenoviral infection in PY and R26Y mice has any long-term effects on the immune microenvironment of the injected mammary glands, we examined leucocyte infiltration at mid- and late-premalignant stages using FACS analysis. At the mid-stage, both adenovirus-infected PY and R26Y mice displayed higher levels of leukocyte infiltration compared with their counterparts induced via TAM (i.e., KPY and KY) (Fig. 2a–b). Induced p53 loss in luminal MECs modestly increased immune activity in the mammary glands of KPY mice relative to KY controls; however, no such difference was detected between adenovirus-induced PY and R26Y mice (Fig. 2a–b). By the late stage, leukocyte infiltration in induced KPY mice rose to levels comparable to those of adenovirus-induced PY and R26Y mice, whereas KY controls remained substantially lower. In contrast, leukocyte infiltration in adenovirus-induced PY and R26Y mice remained similar to each other throughout both stages (Fig. 2a–b). These findings indicate that transient adenoviral infection causes a sustained increase in mammary gland leukocyte infiltration. Furthermore, loss of p53 in luminal MECs likely triggers an anti-tumor immune response that enhances mammary gland immune activity. This effect is readily detectable in TAM-induced KPY mice but may be obscured in adenovirus-induced PY mice due to the elevated baseline immune infiltration caused by prior adenoviral exposure.
To exclude the possibility that the increased immune activation observed in the KPY model following TAM induction was simply due to a higher recombination efficiency compared with adenoviral delivery, we evaluated the initial induction and subsequent clonal expansion of premalignant MECs in both models. We first quantified the proportion of premalignant cells (YFP+) within the lineage-negative compartment (Lin−, i.e., CD45−TER119−CD31−) two weeks after induction. The frequency of YFP+ premalignant cells was comparable between experimental and control groups across both induction strategies (Supplementary figure 2a-b). Moreover, more than 90% of YFP+ cells localized to the luminal MEC gate (Lin−CD24highCD29low), indicating that both induction approaches achieved similar levels of p53 deletion specifically within luminal MECs (Supplementary figure 2c). We next assessed the clonal expansion capacity of p53-null MECs during premalignancy by FACS. In the adenovirus-induced model, p53-null MECs (YFP+) represented 0.29%±0.13% of the Lin− compartment at two weeks post-induction, increasing to 0.85%±0.3% at the mid-stage and 2.89%±0.92% at the late-stage of premalignancy. In contrast, YFP+ MECs in induced R26Y mice (which retain WT p53) remained at a steady percentage (~0.6% YFP+) throughout the premalignant period (Supplementary figure 2d). A similar pattern was observed in the TAM-induced model (Supplementary figure 2e). Together, these results demonstrate that both induction strategies generate comparable initial recombination efficiencies in luminal MECs, and that p53-loss confers a consistent clonal expansion advantage during premalignancy.
To evaluate how transient adenoviral infection reshapes the immune landscape of the mammary gland, we quantified major immune cell populations, including T cells (CD45+CD3+), NK cells (CD45+CD3−NKp46+), macrophages (CD45+CD11b+F4/80+), dendritic cells (CD45+CD11b+CD11c+), B cells (CD45+CD19+), and neutrophils (CD45+CD11b+Ly-6G+), by FACS analysis (Supplementary figure 3). Across all time points, adenovirus-infected mammary glands in both PY and R26Y mice exhibited a marked increase in T cells (Fig. 2c), indicating that the elevated baseline immune infiltration observed following adenoviral injection is driven primarily by T cell enrichment.

Adenoviral infection induces CD8+ tissue-resident T cells in the mammary gland
We next examined the CD3+ T cell compartment in greater detail. Regardless of the p53 mutagenesis status, adenoviral infection markedly increased CD3+ T cell infiltration, both in frequency and absolute cell number, in virally infected mammary glands at both mid- and late-premalignant stages, compared with TAM-induced mice (Fig. 3a). This elevation in T cell abundance was driven primarily by an increase in CD8+ T cells, rather than CD4+ T cells (Fig. 3b–c). Given that transient adenoviral infection produced a sustained elevation in CD8+ T cells, we investigated whether this reflected an expansion of tissue-resident memory T cells. To do so, we assessed CD103 and CD73 expression, canonical markers of tissue-resident memory T cells28,29, on both CD8+ and CD4+ T cell populations. In mice analyzed 3 months post-induction, ~90% of infiltrating CD8+ T cells in adenovirus-infected mammary glands (both PY and R26Y) expressed CD103, compared with only ~30% in TAM-induced KPY and KY mice (Fig. 3d). CD73 expression was relatively high on CD8+ T cells across all groups (Fig. 3d). In contrast, neither CD103 nor CD73 showed substantial differences on CD4+ T cells across induction methods (Fig. 3e). Collectively, these findings indicate that transient adenoviral infection elicits a robust CD8+ T cell response, promoting their infiltration and long-term residency within the mammary gland microenvironment.

CD8+ tissue-resident T cell phenotype masks p53 loss-driven CD8+ T cell activation
We next characterized CD8+ and CD4+ T cell subsets and functional states during premalignancy by assessing expression of IFN-γ, as well as the activation markers PD-1 and CD69, using FACS analysis. In the TAM-induced model (KPY), induced p53-loss in luminal MECs elicited a pronounced increase in IFN-γ, PD-1, and CD69 expression in CD8+ T cells at the late-premalignant stage. At the mid-stage, CD8+ T cells in KPY mice showed increased CD69 expression but no significant elevation in IFN-γ or PD-1 (Fig. 4a–c). In contrast, in the adenovirus-induced model, induced p53-loss in luminal MECs from PY mice did not significantly alter IFN-γ or PD-1 expression in CD8+ T cells when comparing to induced R26Y controls (Fig. 4a–c). Instead, CD8+ T cells in both adenovirus-induced PY and R26Y mice displayed consistently elevated activation states at both mid- and late-stages relative to their TAM-induced counterparts (Fig. 4a–c). This sustained activation profile appears to result from the large pool of tissue-resident CD8+ T cells resulted from the prior adenoviral infection. Supporting this, the majority of PD-1+CD8+ and CD69+CD8+ cells in adenovirus-induced PY and R26Y mice co-expressed CD103, whereas most PD-1+CD8+ and CD69+CD8+ cells in TAM-induced KPY and KY mice lacked CD103 expression (Supplementary figure 4a). In contrast, the CD103 expression in the PD-1+CD4+ and CD69+CD4+ cells from both models was low and not significantly different (Supplementary figure 4b).
In the CD4+ T cell compartment, adenovirus-induced PY and R26Y mice exhibited a transient increase in IFN-γ-producing, activated CD4+ T cells at the mid-premalignant stage, which declined to levels comparable to those of TAM-induced mice by the late stage (Fig. 4d). In contrast, mammary glands from TAM-induced KPY mice displayed a selective increase in IFN-γ-producing, activated CD4+ T cells during the late-premalignant stage compared with KY controls, a feature not observed in the adenovirus-induced model (Fig. 4d). A similar late-stage increase in PD1+CD4+ and CD69+CD4+ T cells was also uniquely observed in induced KPY mice (Fig. 4e–f). Taken together, these findings indicate that transient adenoviral infection recruits CD8+ T cells into the injected mammary gland, where they differentiate into tissue-resident T cells. These virus-induced CD8+ T cells display a persistently activated phenotype that likely obscures the CD8+ T cell response elicited by the expanding p53-null mutant MEC population. In parallel, viral infection appears to drive otherwise quiescent CD4+ T cells into a transient, IFN-γ-producing, virus-specific state, which may limit their engagement in immune responses against p53-null MECs.

Adenoviral infection-induced CD8+ T cell influx condenses the myeloid compartment
Lastly, we examined the myeloid compartment in the mammary gland. In the TAM-induced model, induced p53-loss in luminal MECs resulted in increased macrophage infiltration (both in percentage and absolute number) at mid- and late-premalignant stages (in KPY vs. KY, Fig. 5a–b). In contrast, adenovirus-induced PY and R26Y mice exhibited a pronounced reduction in the macrophage frequency, whereas absolute macrophage numbers remained more comparable to those observed in TAM-induced mice (Fig. 5a–b). These macrophage changes appeared to be driven primarily by adenoviral infection rather than p53 status. Nevertheless, related to p53-loss, a modest late-stage increase in macrophage abundance was also detected in adenovirus-induced PY mice (Fig. 5a–b). Despite these quantitative differences, macrophage polarization marker expression (e.g., CD206, CD86, and I-A/I-E) remained largely similar between adenovirus-induced and TAM-induced models (Fig. 5c–e). However, we observed a trend toward reduced CD206+, MHC-II+ (i.e., I-A/I-E+), or CD86+ macrophage subpopulations in mice with induced loss of p53 (PY or KPY) compared with their respective controls (R26Y or KY) (Fig. 5c–e).

Discussion
In this work, by characterizing and comparing mammary gland immune cells during premalignant stages in breast cancer mouse models with the same driver and tumor outcomes but different induction approaches, we observed that transient adenoviral infection not only led to an expected initial immune response, but also reshaped the mammary gland immune microenvironment long-term. The main viral infection-induced long-lasting change was the profound increase in CD8+ tissue-resident T cells, thereby confounding the detection of p53 loss-induced CD8+ T cell responses. Moreover, transient adenoviral exposure appeared to skew CD4+ T cells toward IFN-γ-producing anti-virus-specific states at the mid-stage. Whether such CD4+ T cell state affects their response to p53-deficient MECs during premalignancy is unclear. Although compared to T cells, the myeloid compartment was less affected by the adenoviral exposure, the expansion of T cell population in the adenovirus-induced model may have compressed the myeloid compartment. In the TAM-induced KPY (vs. KY) mice, induced loss of p53 in luminal MECs led to significant increases in both the percentage and absolute number of macrophages at both mid- and late-stages of premalignancy. However, such changes were less profound in adenovirus-induced PY (vs. R26Y) mice. Nevertheless, both modeling approaches revealed a similar trend of reductions in CD206+, MHC-II+ (i.e., I-A/I-E+), or CD86+ macrophage subpopulations, correlating with induced p53-loss. Overall, this study highlights the complexity of immune cell changes when using different cancer-induction approaches and offers practical insights for selecting an appropriate induction system to achieve more accurate interpretations of immune regulation during mammary tumorigenesis and beyond.
Despite adenovirus-induced changes in the mammary gland immune landscape, adenovirus-induced PY and TAM-induced KPY mice (under similar genetic background) developed mammary tumors with a similar latency17. This does not necessarily suggest that virus-induced changes in mammary gland immune cells have no influence on the course of tumor development. Other factors, such as the time required for p53-deficient MECs to acquire secondary driver mutations, may play a more dominant role in determining the tumor latency. Alternatively, as a tightly regulated host defense system, virus-induced upregulation of CD8+ tissue-resident T cells may be off-set by changes in other immune cell subsets (e.g., macrophages), making the overall immune reactions toward p53-null mutant MECs comparable to those without the virus-induced changes.
Viral infections have been found in both normal and neoplastic human breast epithelial cells. Although infections by viruses such as Epstein-Barr virus (EBV), human cytomegalovirus (HCMV), certain human papillomaviruses (HPV), and mouse mammary tumor viruses (MMTV) are more frequently found in human breast cancers30–33, viral infections of normal human breast epithelial cells have also been reported32,34. For instance, it was reported that HCMV infection could be detected in glandular epithelium in 63% of normal adult breast cases32. In another study, it was shown that EBV infection occurs in breast epithelial cells but not in breast cancer cells34. Interestingly, similar to adenovirus, HCMV is a DNA virus that does not typically integrate its DNA into the host genome and persists in the host cell episomally35. It is possible that viral infections of human breast epithelial cells, even only transient, could also lead to long-last changes in the immune cell landscape in the affected breast tissues, which may have influences on the risk of breast cancer development.
Adenovirus-Cre infection is a commonly used approach to induce cancer initiation in genetically engineered mice, such as those for lung cancer36,37, ovarian/gynecologic cancers38,39, sarcoma40,41, and bladder cancer42,43. It is anticipated that transient adenoviral infection in their corresponding tissues of origin may also affect their immune landscapes, potentially leading to long-lasting immune cell changes locally. Thus, caution should be taken when interpreting immune-related data from cancer mouse models when this induction approach is used.

Methods

Animals
All animal experiments were conducted in accordance with the approved animal protocol (2020N000122) and overseen by the Institutional Animal Care and Use Committee (IACUC) of Brigham and Women’s Hospital (BWH). The reporting in the manuscript follows the recommendations in the ARRIVE guidelines. Trp53fl/fl mice (B6.129P2- Trp53tm1Brn/J, strain# 008462), Rosa26-stop-YFP (R26Y) mice (B6.129X1-Gt(ROSA)26Sortm1(EYFP)Cos/J, strain# 006148), and K8-CreER transgenic mice (Tg(Krt8-cre/ERT2)17Blpn/J, strain# 017947) were purchased from The Jackson Laboratory (JAX). Homozygous Trp53fl/fl;R26Yhomo mice (PY) were obtained by breeding Trp53fl/fl mice with R26Y mice. K8-CreER;Trp53fl/fl;R26Yhomo mice (KPY) were produced by further crossing PY mice with K8-CreER transgenic mice. K8-CreER;R26Yhomo control mice (KY) were generated by breeding the K8-CreER allele into R26Yhomo mice. All mouse lines have been backcrossed into the FVB/NJ background for at least six generations.

Mouse modeling and lineage tracing
Female mice at 8 weeks of age (~15–20g of body weight) were used for lineage tracing. Cre/loxP recombination in luminal mammary epithelial cells (MECs) was induced by two approaches: adenovirus-Cre infection or tamoxifen (TAM) injection. Induced mice were analyzed at 3 or 5 months later (i.e., mid-premalignant stage: ~3 months post-induction, ~20–25g of body weight; late-premalignant stage: ~5 months post-induction, ~25–35g of body weight).
The Keratin 8 (K8) promoter-driven Cre-expressing adenoviral vector was generated in-house and deposited at the University of Iowa Viral Vector Core Facility. The Ad-K8-Cre adenovirus was subsequently obtained from the same facility (WC-Li-535). The concentrated adenovirus was diluted in sterile advanced DMEM/F12 culture medium supplemented with 0.1% Bromophenol blue and 0.01 M CaCl2. 5 μL of the diluted adenovirus prep (109 pfu/mL) was injected into each abdominal mammary gland, as previously described22.
For tamoxifen-induced recombination, tamoxifen (T5648, Sigma-Aldrich, Chicago, USA) was dissolved in corn oil by rotating at 37 °C overnight to a final concentration of 50 mg/mL. Mice received intraperitoneal (i.p.) injections of 5 mg tamoxifen (100 μL) every other day for two doses.

Mammary gland single cell preparation
Mice were euthanized by carbon dioxide (CO2) overdose in Euthanex multi-cage chamber units, which were set to introduce 100% CO2 at a fill rate of 30% displacement of the chamber volume per minute with CO2, added to the existing air in the chamber. Both abdominal (#4) mammary glands were harvested from each individual mouse and then pooled as a single sample. Tissues were minced into fine fragments and digested in digestion media (5 mL advanced DMEM/F12 culture medium containing 2% Fetal Bovine Serum (FBS), 1% HEPES, 16.5 μg Collagenase I, and 75 μg DNase I) at 37 °C for 1 hour with rotation. The resulting single-cell suspension was treated with Red Blood Cell (RBC) lysis buffer for 5 minutes at room temperature (RT) to deplete RBCs. After washing out the lysis buffer, the remaining cells were subjected to flow cytometry analysis.

Flow cytometry (FACS) analysis
To assess cytokine expression via intracellular staining, an aliquot of single-cell suspension from the mammary gland digestion was incubated with a cell activation cocktail in RMPI 1640 supplemented with 10% FBS and 1% Penicillin-Streptomycin for 5 hours in the presence of Brefeldin A. The remaining cells were immediately subjected to antibody staining for the surface markers.
Briefly, the cells were first incubated with fixable viability dye (1:500 dilution) for 20 minutes at RT in the dark to label the dead cells. After centrifugation and removal of the dye solution, the cells were incubated with anti-mouse CD16/32 antibody (1:200; Cat# 156603, clone S17011E; BioLegend, San Diego, USA) at 4 °C for 5 minutes to block non-specific binding. Subsequently, cells were stained with a master mix of fluorophore-conjugated antibodies against surface markers (see Table 1 for details) at 4 °C for 30 minutes, with dilutions as specified by the manufacturers’ protocols. After washing out the antibody dilution, the cells were resuspended and immediately analyzed by flow cytometry.
For intracellular staining, pre-activated cells were harvested and processed for surface marker staining first, as described above. Then, the cells were fixed and permeabilized using the True-Nuclear Transcription Factor Buffer Set (Cat# 424401, BioLegend), and subsequently labeled with the corresponding FACS antibodies.
FACS data were acquired on a BD LSR II or BD Symphony A5 analyzer (BD Biosciences, San Jose, USA) and analyzed using Flowjo software (10.10.0). Gating strategies are provided in Supplementary Fig. 3.

Statistics
The results were presented as mean ± S.E.M. Student’s t tests were used for the comparisons between two groups. Comparisons of multiple groups were performed by two-way ANOVA, and p-values were adjusted using Tukey’s method. No randomization or blinding was used in the in vivo studies. Differences are considered to be significant for *p<0.05, **p<0.01, and ***p<0.005.

Supplementary Material
Supplementary Fig. 1 Schematic diagrams of the genetic approaches used in different mouse models. (a) Ad-K8-Cre intraductal injection model; (b) K8-CreER/TAM model.
Supplementary Fig. 2 Initial induction rates and clonal expansion trajectories based on different induction schemes. Representative FACS plots (a) and statistical analysis (b) of the proportion of YFP+ mutant MECs within the lineage negative (Lin−, i.e., CD31−TER−119−CD45−) compartment. (c) FACS plots represent the CD24/CD29 expression pattern of Lin− cells (right) or Lin−YFP+ mutant cells (left), which defines the luminal MECs (CD24+CD29low), basal MECs (CD24+CD29high), and mesenchymal-like cells (CD24−CD29+). (d, e) The proportion (%) of mutant MECs (YFP+) within the Lin− compartment from each model.
Supplementary Fig. 3 FACS gating strategy for major immune cell populations in the mammary glands.
Supplementary Fig. 4 Tissue-resident CD8+
T cells exhibited an activation phenotype in mammary glands of the virus-induced model. Expression of CD69, PD-1, and CD103 in CD8+ (a) and CD4+ (b) T cells in mammary glands of mice from both models with the indicated genotypes.