- Research article
- Open Access
Activation of the viral sensor oligoadenylate synthetase 2 (Oas2) prevents pregnancy-driven mammary cancer metastases
Breast Cancer Research volume 24, Article number: 31 (2022)
The interferon response can influence the primary and metastatic activity of breast cancers and can interact with checkpoint immunotherapy to modulate its effects. Using N-ethyl-N-nitrosourea mutagenesis, we found a mouse with an activating mutation in oligoadenylate synthetase 2 (Oas2), a sensor of viral double stranded RNA, that resulted in an interferon response and prevented lactation in otherwise healthy mice.
To determine if sole activation of Oas2 could alter the course of mammary cancer, we combined the Oas2 mutation with the MMTV-PyMT oncogene model of breast cancer and examined disease progression and the effects of checkpoint immunotherapy using Kaplan–Meier survival analysis with immunohistochemistry and flow cytometry.
Oas2 mutation prevented pregnancy from increasing metastases to lung. Checkpoint immunotherapy with antibodies against programmed death-ligand 1 was more effective when the Oas2 mutation was present.
These data establish OAS2 as a therapeutic target for agents designed to reduce metastases and increase the effectiveness of checkpoint immunotherapy.
The interferon response can influence the primary and metastatic activity of breast cancers , and can interact with checkpoint immunotherapy to modulate its effects . Although systemic interferon delivery for cancer produced responses in tumors, it was limited by severe adverse events and has been superseded by newer therapies . To circumvent some of these challenges, intratumor injection of agents that induce an interferon response are currently being trialed in a large number of cancer settings and among these agents are derivatives of polyinosinic–polycytidylic acid (polyIC) . In mice, polyIC greatly improved the ability of anti-programmed death-ligand 1 (PD-L1) therapy to clear or reduce melanomas and cancers of the lung and colon, and this effect was dependent on type 1 interferons .
PolyIC and its degradation products mimic viral RNA, key pathogen-associated molecular-patterns displayed during viral replication. PolyIC binds to and activates toll-like receptor 3 (TLR3), a pathogen-associated molecular pattern receptor located on the cell membrane and within endosomes. This generates a number of immune cell effects, such as an interferon response and innate immune cell activation, the latter effect responsible for the use of polyIC as a vaccine adjuvant .
There is, however, another family of pathogen-associated molecular pattern receptors that recognize double-stranded RNA. The OAS family of proteins are potent inducers of the interferon response . They are pathogen-associated molecular pattern receptors that bind double-stranded RNA in the cytoplasm. These proteins respond by synthesizing 2′–5′-linked RNA from ATP without a template . This in turn binds to and activates latent ribonuclease (RNaseL), causing it to degrade viral RNA and host ribosomal RNA, so disrupting the viral life cycle. Activation of OAS-RNaseL also promotes apoptosis and provokes a robust interferon response, both via poorly understood mechanisms that involve additional detection of RNaseL-degraded RNA by the RIG1-like pathogen recognition receptors and possibly direct interaction of OAS and RNaseL with the mitochondria . The OAS family is routinely ignored by most studies using polyIC, which focus on the action of TLR3. The potential of the OAS family as therapeutic targets in oncology has been ignored because the effects of polyIC are attributed to TLR3, and no effective way to separate the actions of OAS and TLR3 was available.
There is some evidence that RNaseL may be involved in the etiology and progression of cancer. Attention has focused on RNaseL, the down-stream effector of the OAS proteins, activation of which predicted the outcomes in virally induced cancers. An association with RNaseL allelic variation and the risk of human papilloma virus-driven cancers of the cervix and head was demonstrated. Unexpectedly, this study also found an association with cancers that do not show a strong viral etiology such as breast . Activation of RNaseL may occur naturally in cancers, hypothesized to be driven by the recognition of double-stranded RNA species by OAS proteins produced from viral sequences integrated in the genome. These sequences are usually silenced by methylation but are transcribed de novo as a result of the general demethylation of the genome that occurs in cancer cells. This was suggested by the finding that cells deficient in RNaseL were highly resistant to the cytotoxicity of 5-azacytidine, a drug that demethylates the genome. Cells were made resistant to 5-azacytidine via inhibition of RNaseL .
Other cancer studies found that allelic variation in RNaseL was observed to be associated with increased risk of prostate cancer and the increased risk of higher tumor grade, together with an increased level of inflammation markers [12, 13]. Allelic variation of RNaseL may modify the risk of breast cancer in patients carrying high risk mutations . The inherited RNaseL alleles responsible for this are proposed to reduce RNaseL activity, and this proposition is supported by observations in lung cancer cells, which increased the levels of the RNaseL inhibitor RL1 to suppress RNaseL-driven apoptosis, with mitochondrial-mediated apoptosis restored by interferon gamma .
Our laboratories conducted an N-ethyl-N-nitrosourea mutagenesis project in mice to find new genes driving mammary gland development and function, with the intent of also identifying new genes involved in mammary cancer initiation and metastasis. We discovered a dominant mutation in Oas2, (referred to as MT, wild type Oas2 referred to as WT) causing a non-conservative amino acid change that resulted in failed lactation despite normal mammary development during pregnancy. These mice were otherwise completely normal . We observed a robust interferon response associated with apoptosis in the mammary glands of our MT mice. When MT was expressed in human breast cancer cells we found that the interferon response depended on the presence of RNaseL and interferon regulatory factor 7 (IRF7), the key effector molecules at the proximal and distal ends of the OAS-RNaseL signaling pathway . This was the first observation of a regulatory pathway linking activation of viral recognition to the control of lactation. The Oas2 mutant mouse provides a defined and single molecular event that induces an interferon response, with pathological activity only in the mammary gland at the onset of lactation.
The MMTV-PyMT model is the most widely used mouse model of breast cancer. PyMT protein  acts as a scaffold promoting persistent signaling via binding of SRC family kinases, PI3K, SHC, 14-3-3, PLC-gamma and TAZ. It causes early and multifocal onset of estrogen receptor positive hyperplasia that rapidly progresses through a series of well-defined lesions to estrogen receptor negative invasive carcinoma of the luminal molecular subtype . The model is weakly sensitive  or refractory  to PD-L1 or PD-1 as monotherapy. Metastasis to lung is frequent and strongly dependent on innate immune cells previously characterized, as macrophages [21, 22] or myeloid-derived suppressor cells (MDSC)  of monocytic (G-MDSC) or neutrophil (PMN-MDSC) origin, and now thought of as pathologically activated monocytes and neutrophils .
To determine if activation of OAS2 by a future breast cancer therapeutic could provide a new therapeutic strategy, and to distinguish OAS activation from that of TLR3, we combined the Oas2 mutation with expression of the polyoma middle T oncogene (PyMT) driven by the mouse mammary tumor virus promoter (MMTV)  and asked if this altered the initiation or progression of mammary cancer, and the interaction with checkpoint immunotherapy.
The Oas2 colony , on an inbred FVBN background, was crossed with the polyoma middle T antigen mouse (MMTV-PyMT) also on an FVBN background, a kind gift from Dr. William Muller . The colony was then maintained by breeding compound heterozygous males (MT Oas2/WT Oas2/MMTV-PyMT/non-transgenic) with WT Oas2/WT Oas2 and non-transgenic FVBN females. All animals were housed with food and water ad libitum with a 12-h day/night cycle at 22 °C and 80% relative humidity.
One or two 7–9-week-old females were housed with single males overnight and monitored daily for a vaginal plug, and pregnancy was confirmed by weight gain. Females were single housed for parturition and pups culled at birth. Mammary gland tumors were measured with calipers throughout pregnancy and the postpartum period. Tumor volume (mm3) was calculated as (minimum measurement2) * (maximum measurement)/2. Mice were euthanized at specified endpoints, via anesthesia (5% isoflurane, 1 L/min O2) followed by cervical dislocation. At autopsy 4th and combined 2nd/3rd glands or tumor weights were recorded.
PDL-1 antibody treatment in parous mice
Parous mice received 10 µl per gram body weight of a solution of PDL-1 (cat BP0101) or IgG control (cat BP0090) at 1 mg/ml via I.P. injection for six treatments, delivered twice weekly.
Mammary gland whole mounts and lung collection
Mouse mammary glands, tumors and lungs were harvested at specified timepoints, and fixed in 10% buffered formalin for 4 h at room temperature. Tissues were changed into 70% ethanol prior to embedding in paraffin. Mammary glands were defatted in 3–4 changes of acetone before being dehydrated and stained in carmine alum for whole mount photography. Glands were then passed through a series of graded alcohols and embedded in paraffin for sectioning and immunohistochemistry.
Tissue sections were either stained with hematoxylin and eosin for routine histochemistry or stained with antibodies to the following antigens using immunohistochemistry protocols as detailed in Additional file 1: Table S1. For analysis of lung metastasis, two or more sections at 100 µm intervals were cut for staining.
Whole mammary gland or lung sections on glass slides were scanned in Aperio ScanScope XT at Katharina Gaus Light Microscopy Facility, and Garvan Microscope Facility, Sydney, Australia, and analyzed by QuPath . A set of multiple reference images was used to train the DAB pixel detection algorithm. Quantification of phospho-Stat positive nuclei used the positive cell detection QuPath algorithm. Regions that contain positive 3,3′-diaminobenzidine (DAB), hematoxylin or no-stain were classified by batch processing and reported as total tissue area (µm2), DAB-stained area (µm2), non-DAB-stained area (µm2) and non-stained area (µm2). DAB stain area percentage without lumen non-stained area was calculated as;
Kaplan–Meier survival analysis, graph drawing and statistical tests were performed in GraphPad Prism 8.0.2. Strong evidence against the observation being made by chance was accepted at the 90% level.
Immune cell profiling by flow cytometry
Mouse mammary glands (2/3rd & 4th) were crudely macerated and cells disassociated in 1% hyaluronidase, 37 °C shaking, 20 min. Digestion was neutralized with 2% fetal calf serum (FCS)/phosphate-buffered saline (PBS) prior to centrifugation at 1200 rpm, 4 °C. Cell pellets were resuspended (25 µg/ml DNAse1 in 2% FCS/PBS) and strained through 40 µm mesh and the supernatant concentrated by centrifugation. Cell pellets were resuspended 0.8% NH4Cl, DNase1 (25 µg/ml) in PBS and incubated in 37 °C water bath, 1 min. Red blood cells (RBC)/DNA lysis was neutralized in 7 × volume 2% FCS/PBS prior to centrifugation and resuspension in 2% FCS/PBS for immune cell staining using fluorophore conjugated antibodies (refer to antibody panels in Additional file 2: Table S2). Mouse spleen, thymus and lymph nodes (from 4th mammary glands) were collected separately into RPMI media (+ 10% FCS) and strained through 70 µm mesh directly prior to RBC/DNA lysis, cell resuspension and immune cell staining steps (as detailed above). Flow cytometry was performed using BD LSR II SORP and data exported to the FlowJo software (Tree Star Inc.) for data analysis. Routine T [27, 28] and innate immune cell protocols  were followed using antibodies listed in Additional file 2: Table S2.
Effects of pregnancy and mutation of Oas2 in the MMTV-PyMT model
MMTV-PyMT mice were made pregnant at 7–9 weeks of age and analyzed between days 1–3 postpartum. Tumors arose in all glands and were generally first detected by palpation in the combined 2nd and 3rd glands. MMTV-PyMT mice showed normal mammary ductal elongation and branching (Fig. 1A), with alveolar bud formation, but failed lobuloalveoli development (FLAD) during pregnancy. At the end of pregnancy, lobuloalveoli would normally completely fill the fat pad and very few adipocytes would be seen, but in MMTV-PyMT glands this does not occur. As a result, there were no observable differences in mammary gland development between parous MT and WT mice. The failure of MMTV-PyMT dams to feed their pups has been noted previously . This is due to a developmental defect and not physical obstruction of the ductal network as often assumed. The mammary ductal tree showed the same characteristic early lesions identified by Lin et al. , hyperplasia, adenoma/MIN and early carcinoma, with late carcinoma dominant at animals that reached the ethical endpoint. These lesions were arranged in a gradient from youngest to oldest epithelium, hyperplasia to carcinoma, left to right in Fig. 1A. IHC showed that PyMT levels were lowest in hyperplasia and highest in adenoma/MIN (Fig. 1B). We observed no histopathological differences between MT and WT mice. Two novel features were observed in parous mice due to the accumulation of milk within the gland (Fig. 1C). In some mice, engorgement with milk of adenoma/MIN produced highly distended ducts, seen on the left-hand side of the combined 2nd and 3rd mammary gland shown in Fig. 1B, which stained strongly for milk-associated protein (MAP). A second feature we have called engorged carcinoma was characterized as misshapen and fluid-filled lumina, formed by a much-thickened fibrous basement membrane containing fibroblasts and immune cells. Epithelial cells accumulated in multiple layers within the lumina. MAP staining was present but weak, suggesting production of the major milk osmoticon lactose, but relative suppression of milk protein synthesis. Areas of early carcinoma and late carcinoma formed within fields of engorged carcinoma, distinguished by the loss of lumina, poorly defined pseudo-alveoli and milk staining limited to developing necrotic areas. Lung metastases showed a uniform pseudo-alveolar pattern typical of the early carcinoma stage. These lesions expanded and developed necrotic centers, but the pseudo-alveolar pattern never varied (Fig. 1D). The metastasis stained strongly for MAP which accumulated within the pseudo-alveoli and in the interstitial spaces (Fig. 1E). Single PyMT positive cells were never seen, the smallest lesions presented as pseudo-alveoli formed from as few as 4 cells, suggesting that metastasis may occur as groups of mammary epithelial cells.
Adenoma/MINs and carcinomas were not engorged with milk at day 15 postpartum or at later times, demonstrating that these features regress dramatically in volume immediately postpartum by reabsorption of the engorging secretions. This is probably responsible for the regressions in tumor volume that we observed by palpation during the postpartum period (Additional file 3: Fig. S1). Tumor regressions and periods of no tumor expansion occurred in 42% and 42% of mutant MT mice, respectively, compared to 21% and 26% of WT mice, respectively. This is consistent with the complete failure of lactation in MT mice  and thus faster loss of these engorged features.
We previously reported that the Oas2 mutation caused disruption of Signal Transducer And Activator Of Transcription (Stat) phosphorylation  and we determined if this also occurred when combined with the MMTV-PyMT oncogene in the immediate postpartum period, when OAS2 levels are highest having been induced through pregnancy  (Fig. 2A). Phosphorylated Stat1 (pStat1) was seen as small foci of DAB + nuclei in all stages of tumor progression and quantification using QuPath  suggested that some MT animals expressed higher pStat1, as reported in non-PyMT animals . Phosphorylated Stat3 (pStat3) was expressed at lower frequency in carcinoma epithelial cells compared to cells from the other stages and more frequently in MT glands. Phosphorylated Stat5 (pStat5) was more frequently expressed in hyperplasia and showed a reduction in the frequency of expression as the disease progressed to carcinoma, correlating with milk protein levels. pStat5 was less frequently expressed in MT glands, as previously reported for non-PyMT animals . Thus, the effects of MT on Stat activation continue in PyMT animals.
We previously reported increased mammary epithelial cell apoptosis in response to MT in non-PyMT mice following parturition . We used staining for cleaved caspase 3 (CC3) to look for this effect in mice also carrying the PyMT oncogene (Fig. 2B). We found no effect when whole glands were analyzed without regard to lesion type prompting an analysis within lesion types. The highest rates of apoptosis were seen in engorged MIN compared to the other lesion types. There were no significant differences in apoptosis between MT and WT animals in any of the lesion types. This suggests that transformation with PyMT abrogates the apoptotic effects of OAS2 activation.
Effects of MT OAS2 on tumor progression and metastasis
Cohorts of 20 WT or MT PyMT Oas2 mice were kept nulliparous or were made pregnant (80 mice total) from 7 to 9 weeks of age and all pups were removed at birth. Mice were monitored by mammary palpation at twice-weekly intervals. Tumor onset and growth was measured in the major and minor axes using calipers, and mice were euthanized at the estimated ethical endpoint of 10% tumor burden prior to autopsy. Animal survival to various endpoints was analyzed using Kaplan–Meier survival plots with log-rank p values and hazard ratios calculated by the Cox proportional hazards model as implemented in GraphPad-Prism (Fig. 3, and Additional file 4: Fig. S2 for p values and hazard ratios (HR) of two-way comparisons).
Survival to the age of initial tumor detection (Fig. 3A) showed that tumors in WT parous animals were detected at a younger age than WT nulliparous animals, 66 compared to 78 days (p = 0.0038). MT parous animals also had tumors detected at a younger age than MT nulliparous animals, 65.6 days compared to 77.5 days (p = 0.0019). There was no difference in the age of tumor detection between WT and MT animals, either parous or nulliparous. The period from tumor detection to ethical endpoint (Fig. 3B) was not different among WT nulliparous, WT parous or MT parous animals. MT nulliparous animals progressed more quickly than those in the other cohorts, 22.5 days compared to 31 days (p = 0.002) for WT parous, 35 days (p < 0.0001) for MT parous and 25.5 days (p = 0.03) for WT nulliparous. The age that animals reached their ethical endpoint (Fig. 3C) was not different among WT parous, MT parous and MT nulliparous animals. WT nulliparous reached their ethical endpoint later by 5 days (p = 0.0345), 3 days (p = 0.034) and 4 days (p = 0.18), respectively. Overall, parity caused earlier tumor detection but slower tumor progression, opposing effects with the result that parous animals, either MT or WT, reached the ethical endpoint a few days earlier than nulliparous WT animals. The Oas2 mutation had no effect on the age of tumor detection but caused more rapid progression in nulliparous animals, so that MT nulliparous animals reached the ethical endpoint earlier than nulliparous WT, and at the same age as parous animals. These effects showed HR between 1.8 and 2.9 indicating a two- to threefold increase in risk.
We examined metastatic burden in two sections of both lungs per mouse from the cohort of mice above, using PyMT IHC to detect metastases and allow measurement of their stained area using QuPath . In the cohorts, there was no difference in tumor burden at the ethical endpoint, measured at autopsy by weighing tumors and body weights (BW) from animals euthanized at an estimated 10% tumor burden by tumor dimensions (Fig. 4A, each symbol is a mouse).
Results for the quantification of metastases are plotted on a log scale due to the skewed and thus non-Gaussian distribution of this data, p < 0.001 K2 33–70 using the D’Agostino&Pearson normality test (Fig. 4B). Sections with no lung metastases are indicated as “0” at the bottom of the figure as they cannot be plotted on a log axis. Error bars show standard error of the mean, but median values (dashed lines) and the Mann–Whitney nonparametric test of significance are used for comparisons.
Pregnancy caused an increase in the proportion of WT mice with metastases, up from 9 of 19 nulliparous WT mice to 15 of 19 parous WT mice (p = 0.09 Fisher’s exact test). In contrast, pregnancy made no difference to the number of MT Oas2 mice with a metastasis, 10 of 20 nulliparous MT mice compared to 8 of 18 parous MT mice (p = 0.8 Fisher’s exact test). There were half the number of MT mice with metastases than WT mice (p = 0.04 Fisher’s exact test).
The median number of metastases per lung section (Fig. 4B top panel) rose from 0 in nulliparous WT or MT mice to 6.5 in parous WT mice (p < 0.0001 Mann Whitney), but to just 1 in parous MT mice (p = 0.0016). Parous MT mice had sevenfold fewer metastases than WT mice (p = 0.0124). Image analysis was used to quantify the area of lung metastases as a proportion of lung tissue area (Fig. 4B lower panel). The area of lung alveolar space was excluded from the analysis as individual sections showed variation in the degree of lung inflation at fixation. Parity caused the median area of metastasis in WT lung tissue to increase from 0 to 1.3%, but to just 0.18% in MT animals, a sevenfold decrease in metastases area (p = 0.0338). We conclude that activation of OAS2 prevented pregnancy from driving metastases to the lung.
Effects of the OAS2 mutation on the immune system
Lung metastasis in the MMTV-PyMT model is highly dependent on innate immune cells [18, 21,22,23] of monocyte and neutrophil origin which are recruited to tumors where they adopt a pathogenic inflammatory and immunosuppressive state capable of modulating T-cell activity, often called MDSC of granulocyte or polymorphonucleated neutrophil origin (G- or PMN-MDSC), or monocyte origin (M-MDSC). The development of cell surface markers enabling clear distinction of cells that have assumed the pathogenic state remains the key challenge of this field, and the definition of the pathogenic state is currently limited to functional tests .
We undertook a flowcytometric survey of monocytes and granulocytes in tumors (Fig. 5, gating strategy Additional file 6: Fig. S4). Monocyte and neutrophil levels showed individual variation; however, we found a rise in the number of monocytes and neutrophils in more than half of the MT parous tumors compared to WT parous tumors and mean levels increased for both cell types (p = 0.022 and p = 0.0296, respectively). Non-oncogene carrying mice showed no significant changes in monocyte levels with MT or parity. Neutrophil levels increased significantly in parous animals compared to nulliparous animals (p = 0.0571 and p = 0.0043) but MT was without effect. The presence of a tumor appeared to reduce the ability of parity to increase neutrophil levels in some mice. As monocyte and neutrophil levels increase in the tumors of parous MT mice, but not in all MT mice, these results suggest that any action of MT on monocytes and neutrophils to reduce metastasis must be to prevent the transition of monocytes and neutrophils to the pathogenic state.
To determine if the changes in monocytes and neutrophils altered parameters measured in CD4 and CD8 T cells, we surveyed tumors, the draining inguinal lymph node, spleen and thymus. We found great variation among individual mice in most of the measured T-cell parameters, showing that we would require an impractically large cohort size to provide sufficient power to this analysis to detect differences (Additional file 5: Fig. S3).
Effects of PD-L1 on tumor progression and metastasis
To determine if PD-L1 treatment was more effective in parous mice and whether it interacted with MT, we added PD-L1 treatment to the paradigm used above. Mice received either six IgG or PD-L1 treatments (individual details Additional file 6: Fig. S4) when tumors first became palpable and generally coinciding with pregnancy and parity. Mice were monitored with twice weekly measurements of major and minor tumor axis to the ethical end point. Data was analyzed by Kaplan–Meier survival analysis with log-rank p values and Cox proportional hazards estimation of hazard ratios. All tumor growth profiles are shown in in Additional file 6: Fig. S4 and all two-way survival comparisons are shown in Additional file 7: Fig. S5.
PD-L1 or MT had no effect on the age at which tumors were first detected (Fig. 6A, two-way comparisons and p values Additional file 8: Fig. S6). PD-L1 treatment slowed tumor progression in MT mice from detection to endpoint, by 10 days compared to IgG treated mice (Fig. 6B), from a median time of 46 to 56 days (p = 0.0036 HR = 3.8). The effect of PD-L1 in MT mice was also apparent at the ethical endpoint (Fig. 6C), where PD-L1 treated MT mice had a median age of 62.5 days compared to MT IgG treated of 67 days (p = 0.065 HR 1.9). PD-L1 had no effect in WT mice. Tumors in MT IgG treated mice progressed more quickly to the ethical endpoint than in WT mice, and PD-L1 treatment normalized this so that MT PD-L1 treated mice became indistinguishable from WT mice. Since PD-L1 was ineffective in WT mice on tumor progression, we conclude that MT sensitized the mice to PD-L1 treatment, but by promoting more rapid tumor progression which was then ameliorated by PD-L1.
We examined the effect of PD-L1 treatment on metastasis (Fig. 7). Tumor burden was estimated from tumor dimensions and mice were killed at an estimated 10% BW. Actual tumor burden was measured at autopsy. PD-L1 produced slightly lower burden in both WT (p = 0.0813) and MT mice (p = 0.0939). In WT animals, PD-L1 treatment did not significantly reduce the number of metastases, but in MT animals, a significant sevenfold reduction occurred from a median of 29.5 to 4 metastases per section (p = 0.001). The median area of lung metastases as a proportion of lung tissue area was reduced ninefold in WT animals from 6.3 to 0.7 (p = 0.002) and fivefold in MT animals from 3 to 0.56 (p = 0.0002). We conclude that MT increased the ability of PD-L1 treatment to reduce metastases number, a product of the seeding events to lung, but that no effect was seen in tumor area, a product of seeding and additional metastasis growth. Although the interaction of OAS MT and PD-L1 in this model is of small magnitude, it was detected in a model that is at best partially sensitive to PD-L1, hardly the ideal model to examine this question, but never-the-less provides a finding to be followed in other cancer types.
To determine if PD-L1 was having a greater but transient effect during the treatment period, we examined a second cohort of six mice per group immediately after treatment using 2 sections per lung (Fig. 8). Mammary gland weights showed no difference among groups (Fig. 8A). The only group in which any lung sections contained no metastases were MT Oas2 and PD-L1 treated mice, where half the sections were negative (p = 0.0373 Fischer’s exact test compared to WT PD-L1 treated). The median number of metastases in WT mice increased 2.8-fold from 16.5 to 45.5, but in MT mice it fell 9.2-fold from 23 to 2.5 (Fig. 8B upper panel). There was a large difference in median metastasis number per section between WT (45.5) and MT (2.5) animals treated with PD-L1 (p = 0.128).
Metastases as a percentage of tissue area fell 3.5-fold in WT mice from 3.36% to 0.95% (p = 0.0433) and 16-fold in MT mice from 4% to 0.25% (p = 0.0001) (Fig. 8B lower panel). PD-L1 was more effective in MT Oas2 mice. Interestingly, the number of metastases seen two weeks after parity was similar to those detected much later at the ethical endpoint, showing that pregnancy and parity increase metastases within an acute period, rather than over the full length of the postpartum period. These data show that PD-L1 treatment caused a reduction of metastases number and area in MT mice.
In summary, MT prevented the increase in the number and the tissue area of lung metastases caused by parity and MT enhanced the response to PD-L1 therapy.
The publications on allelic variation in the Oas-RNaseL pathway suggest that inherited alleles that are thought to cause pathway activation produce better control of virally induced tumors, but also of non-viral tumors by providing an enhanced response to de novo dsRNA production from integrated viruses. A separate group of Oas-RNaseL pathway alleles is thought to cause pathway suppression and to be involved in the genesis of prostate cancer. Our results suggest, however, that the same Oas2 alleles may have context-dependent effects. In the nulliparous state, activation of Oas2 caused faster primary tumor progression, but in the parous state this effect was lost, and instead suppression of metastases and sensitization to checkpoint immunotherapy was observed. It is likely that distinct effects of Oas-RNaseL activation within the cells of the epithelium and immune system add complexity to the role of this pathway in cancer.
A fully completed pregnancy remains an independent predictor of breast cancer risk , and pregnancy-associated breast cancer is associated with a 40% increase in mortality due to increased metastasis , with even higher rates in cases occurring during pregnancy . The period of increased risk caused by pregnancy may last three or more decades for women who have their first child after their mid 30’s, but less than a decade for women who have their first child before their early 20’s [32, 33]. A degree of life-long protection then follows. This age-dependent period of increased risk defines a large proportion of breast cancer cases as associated with pregnancy; however, current clinical practice uses arbitrary postpartum time periods to define pregnancy associated breast cancer, and does not distinguish pregnancy-associated breast cancer based on the age-dependent period of increased risk. Regardless, this type of breast cancer has no special therapeutic guidelines.
The number of first pregnancies of women over 30 years of age continues to increase especially in developed and developing countries where increasing age at first pregnancy is correlated with increasing female education [34, 35]. Breast cancer associated with pregnancy  is one of the main drivers of the increased prevalence of breast cancer . This effect is compounded by an ongoing decrease in breast feeding, as every 12 months decreases breast cancer risk by 4.3%. In developed countries, decreases in childbirth and lactation have doubled the cumulative risk of breast cancer, with two-thirds of this effect attributable to reduced lactation .
The results of this study have direct relevance to the changes in reproduction practices in developed and developing nations. Activation of OAS2 in the MMTV-PyMT model of breast cancer prevented pregnancy from increasing the number of lung metastases, a finding that establishes the OAS family as therapeutic targets for the development of agents designed to prevent pregnancy from driving breast cancer metastasis. A small interaction with PD-L1 therapy was also detected, albeit in a model at best only partially sensitive to PD-L1. The data provide impetus for further study to determine if the combination of OAS activation and checkpoint inhibition is of therapeutic relevance for the treatment of breast cancer during pregnancy, or other cancers where a more robust effect of PD-L1 occurs. A cautionary finding was that in nulliparous mice activation of OAS2 resulted in more rapid primary tumor progression, suggesting that in nulliparous women the use of interferon-inducing agents should be carefully examined, especially for local recurrences. As the mechanisms responsible for the effects of pregnancy and lactation on postpartum breast cancer risk are poorly understood these findings may also have significance for our understanding of pregnancy associated breast cancer. Importantly, the MT mice show no pathology other than failed lactation, a rare observation in the interferon-response field which is beset with serious adverse events, suggesting that such a therapy would have a favorable toxicity profile.
Our findings may have relevance to breast cancer more generally. We have recently described a pro-metastasis pathology we called involution mimicry  that occurs in response to forced expression of the ETS transcription factor ELF5, the master regulator of lobuloalveolar development during pregnancy, which causes a limited lactation and ongoing involution response in an ELF5 mouse model , but is also seen in breast cancer patients. The ELF5 mouse model shows a similar a engorgement pathology to that we saw following parity in the MT mice and an increase in lung metastasis that is abrogated by neutrophil depletion . We found profound changes in the state of the epithelium, cancer-associated fibroblasts (CAFs), endothelial cells and immune cells . Among the CAFs we found populations previously associated with mammary gland involution following lactation, and an analysis of ligand-receptor interactions among these cell types established the existence of a complex regulatory network linking the induction of ELF5 in the epithelium to profound transcriptional changes in the CAFs and endothelial cells, and to the recruitment and M2 polarization of MDSC . Two linked human molecular signatures were discovered in breast cancer patients, a lactation signature and an involution signature. The lactation signature was enriched in the epithelial cell populations and the involution signature was present in both epithelial and stromal populations. The involution signature was prognostic in the METABRIC cohort . We surmise that elevated ELF5, possibly due to loss of promoter demethylation , drives a limited lactation and involution in these breast cancers.
Our results raise the possibility that the use of an OAS family activator, or less specific inducers of the interferon response, could reduce the pro-metastases pathology that occurs in a significant number of breast cancers during pregnancy, the immediate postpartum period and in cases expressing the involution-mimicry phenotype, and that efficacy could be increased by combination with checkpoint immunotherapy. While the magnitude of some of the effects detailed here are modest, that they occur indicates the existence of mechanistic links that may be exploited by new therapeutic strategies.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Cleaved caspase 3
E-74 like factor 5
E-26 transformation specific
Fetal calf serum
Failed lobuloalveolar development
Immuno gamma globulin
Interferon regulatory factor 7
Myeloid-derived suppressor cells
Mammary intraepithelial neoplasia
- MMTV :
Mouse mammary tumor virus promoter
- Oas2 :
Oligoadenylate synthetase 2 gene
Oligoadenylate synthetase 2 protein
Programmed death receptor 1
Programmed death-ligand 1
Phospho lipase C
Polyoma middle T protein
- PyMT :
Polyoma middle T oncogene
Red blood cells
Retinoic acid induced gene 1
Roswell Park Medical Institute
Src homology 2 domain-containing
Rous sarcoma virus family kinases
Parker BS, Rautela J, Hertzog PJ. Antitumour actions of interferons: implications for cancer therapy. Nat Rev Cancer. 2016;16:131–44.
Benci JL, et al. Opposing functions of interferon coordinate adaptive and innate immune responses to cancer immune checkpoint blockade. Cell. 2019;178:933–48.
Borden EC. Interferons alpha and beta in cancer: therapeutic opportunities from new insights. Nat Rev Drug Discov. 2019;18:219–34.
Kyi C, et al. Therapeutic immune modulation against solid cancers with intratumoral poly-ICLC: a pilot trial. Clin Cancer Res. 2018;24:4937–48.
Nagato T, Lee YR, Harabuchi Y, Celis E. Combinatorial immunotherapy of polyinosinic-polycytidylic acid and blockade of programmed death-ligand 1 induce effective CD8 T-cell responses against established tumors. Clin Cancer Res. 2014;20:1223–34.
Zheng X, Li S, Yang H. Roles of toll-like receptor 3 in human tumors. Front Immunol. 2021;12:667454.
Silverman RH. Viral encounters with 2′,5′-oligoadenylate synthetase and RNase L during the interferon antiviral response. J Virol. 2007;81:12720–9.
Donovan J, Dufner M, Korennykh A. Structural basis for cytosolic double-stranded RNA surveillance by human oligoadenylate synthetase 1. Proc Natl Acad Sci USA. 2013;110:1652–7.
Chakrabarti A, Jha BK, Silverman RH. New insights into the role of RNase L in innate immunity. J Interferon Cytokine Res. 2011;31:49–57.
Madsen BE, et al. Germline mutation in RNASEL predicts increased risk of head and neck, uterine cervix and breast cancer. PLoS ONE. 2008;3:e2492.
Banerjee S, et al. OAS-RNase L innate immune pathway mediates the cytotoxicity of a DNA-demethylating drug. Proc Natl Acad Sci USA. 2019;116:5071–6.
Fesinmeyer MD, Kwon EM, Fu R, Ostrander EA, Stanford JL. Genetic variation in RNASEL and risk for prostate cancer in a population-based case-control study. Prostate. 2011;71:1538–47.
Meyer MS, et al. Genetic variation in RNASEL associated with prostate cancer risk and progression. Carcinogenesis. 2010;31:1597–603.
Nguyen-Dumont T, et al. Is RNASEL:p.Glu265* a modifier of early-onset breast cancer risk for carriers of high-risk mutations? BMC Cancer. 2018;18:165.
Yin H, Jiang Z, Wang S, Zhang P. IFN-gamma restores the impaired function of RNase L and induces mitochondria-mediated apoptosis in lung cancer. Cell Death Dis. 2019;10:642.
Oakes SR, et al. A mutation in the viral sensor 2′–5′-oligoadenylate synthetase 2 causes failure of lactation. PLoS Genet. 2017;13:e1007072.
Fluck MM, Schaffhausen BS. Lessons in signaling and tumorigenesis from polyomavirus middle T antigen. Microbiol Mol Biol Rev. 2009;73:542–63.
Lin EY, et al. Progression to malignancy in the polyoma middle T oncoprotein mouse breast cancer model provides a reliable model for human diseases. Am J Pathol. 2003;163:2113–26.
Allen E, et al. Combined antiangiogenic and anti-PD-L1 therapy stimulates tumor immunity through HEV formation. Sci Transl Med. 2017;9:eaak9679.
Messenheimer DJ, et al. Timing of PD-1 blockade is critical to effective combination immunotherapy with anti-OX40. Clin Cancer Res. 2017;23:6165–77.
Lin EY, et al. Vascular endothelial growth factor restores delayed tumor progression in tumors depleted of macrophages. Mol Oncol. 2007;1:288–302.
Lin EY, Nguyen AV, Russell RG, Pollard JW. Colony-stimulating factor 1 promotes progression of mammary tumors to malignancy. J Exp Med. 2001;193:727–40.
Gallego-Ortega D, et al. ELF5 drives lung metastasis in luminal breast cancer through recruitment of Gr1+ CD11b+ myeloid-derived suppressor cells. PLoS Biol. 2015;13:e1002330.
Veglia F, Sanseviero E, Gabrilovich DI. Myeloid-derived suppressor cells in the era of increasing myeloid cell diversity. Nat Rev Immunol. 2021;21:485–98.
Guy CT, Cardiff RD, Muller WJ. Induction of mammary tumors by expression of polyomavirus middle T oncogene: a transgenic mouse model for metastatic disease. Mol Cell Biol. 1992;12:954–61.
Bankhead P, et al. QuPath: open source software for digital pathology image analysis. Sci Rep. 2017;7:16878.
Bradley LM, Atkins GG, Swain SL. Long-term CD4+ memory T cells from the spleen lack MEL-14, the lymph node homing receptor. J Immunol. 1992;148:324–31.
MacDonald HR, Budd RC, Cerottini JC. Pgp-1 (Ly 24) as a marker of murine memory T lymphocytes. Curr Top Microbiol Immunol. 1990;159:97–109.
Rose S, Misharin A, Perlman H. A novel Ly6C/Ly6G-based strategy to analyze the mouse splenic myeloid compartment. Cytometry A. 2012;81:343–50.
Daling JR, Malone KE, Doody DR, Anderson BO, Porter PL. The relation of reproductive factors to mortality from breast cancer. Cancer Epidemiol Biomark Prev. 2002;11:235–41.
Whiteman MK, et al. Reproductive history and mortality after breast cancer diagnosis. Obstet Gynecol. 2004;104:146–54.
Albrektsen G, Heuch I, Hansen S, Kvale G. Breast cancer risk by age at birth, time since birth and time intervals between births: exploring interaction effects. Br J Cancer. 2005;92:167–75.
Colditz GA, Rosner B. Cumulative risk of breast cancer to age 70 years according to risk factor status: data from the Nurses’ Health Study. Am J Epidemiol. 2000;152:950–64.
Mathews T, Hamilton B. Mean age of mothers is on the rise: United States, 2000–2014. NCHS Data Brief No 232. 2016.
Studies AIoF. Births in Australia—age of new mothers. 2018.
Lee GE, Mayer EL, Partridge A. Prognosis of pregnancy-associated breast cancer. Breast Cancer Res Treat. 2017;163:417–21.
DeSantis CE, et al. International variation in female breast cancer incidence and mortality rates. Cancer Epidemiol Biomark Prev. 2015;24:1495–506.
Collaborative Group on Hormonal Factors in Breast C. Breast cancer and breastfeeding: collaborative reanalysis of individual data from 47 epidemiological studies in 30 countries, including 50302 women with breast cancer and 96973 women without the disease. Lancet. 2002;360:187–95.
Valdes-Mora F, et al. Single-cell RNAseq uncovers involution mimicry as an aberrant development pathway during breast cancer metastasis. Cell Rep. 2021;35:108945.
Lee HJ, et al. Lineage specific methylation of the Elf5 promoter in mammary epithelial cells. Stem Cells. 2011;29:1611–9.
The authors acknowledge the excellent work of the staff of Australian BioResources and the Biological Testing Facility of the Garvan Institute.
This work was supported by grants from the Congress Directed Medical Research Program (BC995364 and DAMD17-01-1-0241), Cure Cancer Australia Foundation, NHMRC Australia (projects 1047149, Fellowships 545805, 481310, 1043400), the Australian Research Council Discovery Project (DP110102288), Fellowships (ECF-13-08 and ECF-16-022) from the National Breast Cancer Foundation, Banque Nationale de Paris-Paribas Australia and New Zealand, Mostyn Family Foundation, Cue Clothing Co., Estee Lauder Australia, Sutton’s Motors Sydney and the RT Hall Trust.
Ethics approval and consent to participate
All mice were housed in specific pathogen-free conditions at the Australian BioResources and the Biological testing Facility of the Garvan Institute, with all animal experiments carried out according to guidelines contained within the NSW (Australia) Animal Research Act 1985, the NSW (Australia) Animal Research Regulation 2010 and the Australian code of practice for the care and use of animals for scientific purposes, (8th Edition 2013, National Health and Medical Research Council (Australia)) and approved by the Garvan/St Vincent’s Animal Ethics and Experimentation Committees (approvals 14/27, 16/26, 17/23, 18/36).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Antibodies, concentration and antigen retrieval conditions for immunohistochemistry. All reagents were from Leica BOND or DAKO for automated or manual IHC (as specified). Visualization of antigen–antibody complexes was performed using the DAB + liquid substrate chromogen system (K3467).
Antibodies used for flow cytometry. Panel name, antibody name by antigen recognized and conjugated fluorochrome, the antibody supplier and catalogue number and final diluted concentration for each is shown.
Effect of MT Oas2 on primary tumor expansion in parous mice. Examples of tumor growth curves of individual mice. The left- and right-hand side 2nd/3rd and 4th mammary glands were palpated twice weekly and tumor growth was estimated by measurement of the major and minor axis of each gland using calipers. Y axis is placed at parity indicated on the x axis as day 0. Numbers show the total number of glands with tumors detected, the number that showed a period of no growth, and the number that showed regression at or following parity. W, wild type; M, MT; W/W, homozygous wild type; M/M, homozygous MT.
Effect of MT Oas2 on primary tumor initiation and expansion in nulliparous and parous mice. Two-way Kaplan–Meier survival analysis for the indicated periods and endpoints, and for the indicated genotypes, WT wildtype Oas2, MT mutant Oas2, and parity status of mouse cohorts (20 per group). P values and hazard ratios (HR) calculated by the log-rank test using GraphPad Prism. n = 19 MT parous, 21 WT parous, 20 MT nulliparous and 19 WT nulliparous.)
Effects of MT Oas2 on T cells in parous mice. A screening of various T-cell parameters was undertaken in mammary tumors, lymph nodes, spleen and thymus. Error bars are standard error of the mean and p values calculated by Student’s t test (all non-significant). Axes represent % of total cells passing the previous gate. Genotypes, WT/WT homozygous wildtype Oas2, MT/MT homozygous mutant Oas2.
Flow cytometry gating strategy for monocytes and neutrophils known as myeloid-derived suppressor cells. Series of gates used to quantify monocytes and neutrophils using flow cytometry.
Effect of PD-L1 and MT Oas2 on primary tumor initiation and expansion in parous mice. Examples of tumor growth curves of individual mice. The 2nd/3rd and 4th mammary glands were palpated twice weekly and tumor growth was estimated by measurement of the major and minor axis of individual tumors using calipers. Y axis is placed at the day of tumor detection on the x axis as day 0. Black diamonds show a treatment with IP IgG or PD-L1 antibodies and red diamonds indicate parity. WT = wild type, MT = MT Oas2 MT/MT = homozygous MT Oas2.
Effect of PD-L1 and MT Oas2 on primary tumor initiation and expansion in parous mice. Two-way Kaplan–Meier survival analysis for the indicated periods and endpoints, and for the indicated genotypes and treatments (PD-L1 or IgG IP), of mouse cohorts. P values and hazard ratios (HR) calculated by the log-rank test using GraphPad Prism. n = 8 WT PD-L1, 8 WT IgG, 14 MT PD-L1, 15 MT IgG.
About this article
Cite this article
Ho, WH.J., Law, A.M.K., Masle-Farquhar, E. et al. Activation of the viral sensor oligoadenylate synthetase 2 (Oas2) prevents pregnancy-driven mammary cancer metastases. Breast Cancer Res 24, 31 (2022). https://doi.org/10.1186/s13058-022-01525-z