Mutant mice lacking alternatively spliced p53 isoforms unveil Ackr4 as a male-specific prognostic factor in Myc-driven B-cell lymphomas

  1. Anne Fajac
  2. Iva Simeonova
  3. Julia Leemput
  4. Marc Gabriel
  5. Aurélie Morin
  6. Vincent Lejour
  7. Annaïg Hamon
  8. Jeanne Rakotopare
  9. Wilhelm Vaysse-Zinkhöfer
  10. Eliana Eldawra
  11. Marina Pinskaya
  12. Antonin Morillon
  13. Jean-Christophe Bourdon
  14. Boris Bardot  Is a corresponding author
  15. Franck Toledo  Is a corresponding author
  1. Genetics of Tumor Suppression, Institut Curie, France
  2. CNRS UMR3244, France
  3. Sorbonne University, France
  4. PSL Research University, France
  5. Non Coding RNA, Epigenetic and Genome Fluidity, Institut Curie, France
  6. School of Medicine, Ninewells Hospital, University of Dundee, United Kingdom

Abstract

The Trp53 gene encodes several isoforms of elusive biological significance. Here, we show that mice lacking the Trp53 alternatively spliced (AS) exon, thereby expressing the canonical p53 protein but not isoforms with the AS C-terminus, have unexpectedly lost a male-specific protection against Myc-induced B-cell lymphomas. Lymphomagenesis was delayed in Trp53+/+Eμ-Myc males compared to Trp53ΔAS/ΔAS Eμ-Myc males, but also compared to Trp53+/+Eμ-Myc and Trp53ΔAS/ΔAS Eμ-Myc females. Pre-tumoral splenic cells from Trp53+/+Eμ-Myc males exhibited a higher expression of Ackr4, encoding an atypical chemokine receptor with tumor suppressive effects. We identified Ackr4 as a p53 target gene whose p53-mediated transactivation is inhibited by estrogens, and as a male-specific factor of good prognosis relevant for murine Eμ-Myc-induced and human Burkitt lymphomas. Furthermore, the knockout of ACKR4 increased the chemokine-guided migration of Burkitt lymphoma cells. These data demonstrate the functional relevance of alternatively spliced p53 isoforms and reveal sex disparities in Myc-driven lymphomagenesis.

eLife assessment

This important study using engineered mouse models provides a first and compelling demonstration of a pathogenic phenotype associated with lack of expression of p53AS isoforms, isoforms of the p53 protein with a different C-terminus than canonical p53. The role of these isoforms has been elusive so far and this first demonstration represents a substantial advance in our understanding of the complex role(s) of p53 isoforms. The revised article adequately addresses previous concerns.

https://doi.org/10.7554/eLife.92774.3.sa0

eLife digest

Human cells divide many times during a lifetime, a process that requires careful regulation to avoid uncontrolled cell division, which can lead to various disorders, including cancer. For example, TP53, which encodes multiple proteins, is the most commonly mutated gene in cancers.

TP53 carries the instructions to make a tumor suppressor protein, known as p53, which can stop cancers from forming and spreading. In humans and mice, TP53 (and the mouse analogue Trp53) can also be read by the cell to make several slightly different versions of the p53 protein, known as isoforms. The p53 isoforms are much less studied and their role in an organism is still unclear.

To address this, Fajac et al. used genome editing to make mouse strains that were still able to express p53, but were only able to create a specific subset of p53 isoforms. In these mice, part of the Trp53 gene had been mutated to remove the cell’s ability to make isoforms with an alternative C-terminal end.

Fajac et al. then allowed these mice to breed with mice that were model organisms for a cancer called B-cell lymphoma. This revealed that male offspring that lacked alternative p53 isoforms were more susceptible to B-cell lymphoma and that they had decreased levels of the protein ACKR4, a receptor for signaling proteins that regulate cellular movement. Human datasets showed that having higher levels of ACKR4 could be linked to a better disease prognosis in male patients with Burkitt lymphoma, a rare but aggressive form of B-cell lymphoma. The same effect was not observed in females, suggesting that measuring ACKR4 gene expression in male patients with Burkitt lymphoma might be useful to identify the patients at higher risk.

The study from Fajac et al. provides a new perspective on p53 – one of the most studied proteins. It highlights specific p53 isoforms and the ACKR4 protein as a potential way to identify male patients at higher risk from a type of B-cell lymphoma.

Introduction

TP53, the human gene for tumor suppressor p53, encodes several isoforms owing to distinct promoters, alternative splicing, and multiple translation initiation sites (Bourdon et al., 2005; Courtois et al., 2002; Flaman et al., 1996; Yin et al., 2002). p53 alternative isoforms can be abnormally expressed in cancer cells and some may regulate the canonical p53 protein (Anbarasan and Bourdon, 2019; Bourdon et al., 2005; Mondal et al., 2013; Senturk et al., 2014). However, aberrant RNA splicing is a common feature of cancer cells (Graubert et al., 2012; Martin et al., 2013; Pajares et al., 2007; Sette and Paronetto, 2022) and to which extent alternative splicing generates functionally relevant proteins is controversial (Abascal et al., 2015; Bardot and Toledo, 2015; Blencowe, 2017; Tress et al., 2017a; Tress et al., 2017b; Ule and Blencowe, 2019; Weatheritt et al., 2016). Thus, the biological importance of many p53 isoforms remains elusive.

Like its human TP53 homolog, the murine Trp53 gene encodes multiple isoforms differing in their N- or C-termini (Arai et al., 1986; Marcel et al., 2011). Mouse models to evaluate the role of p53 isoforms differing in their N-terminus revealed that Δ40-p53 overexpression leads to accelerated ageing (Maier et al., 2004; Steffens Reinhardt et al., 2020). However, the potential role of p53 isoforms with an alternative C-terminus was not analyzed in vivo. p53 isoforms with distinct C-termini result from the splicing of two mutually exclusive final exons: exon 11, encoding the canonical ‘α’ C-terminal domain, and the alternatively spliced (AS) exon, encoding another C-terminus (Arai et al., 1986). In adult mice, isoforms with the canonical C-terminus are predominant in all tissues (Figure 1—figure supplement 1A). Two models (Trp53Δ31 and Trp53ΔCTD), designed to study the consequences of a loss of the canonical p53 C-terminus, exhibited signs of increased p53 activity, leading to a rapidly lethal anemia (Hamard et al., 2013; Simeonova et al., 2013). To determine the role of p53-AS isoforms in vivo, we created Trp53ΔAS, a mouse model with a specific deletion of the AS exon (Figure 1—figure supplement 1B). In mouse embryonic fibroblasts (MEFs), the Trp53ΔAS allele prevented the expression of isoforms with the AS C-terminus, whereas it did not affect RNA levels for p53 isoforms with the canonical C-terminus (Figure 1—figure supplement 1C). We previously used this model to show that p53-AS isoforms had no role in the anemia affecting Trp53Δ31/Δ31 mice (Simeonova et al., 2013). However, a detailed phenotyping of Trp53ΔAS/ΔAS mice remained to be performed. The detailed phenotyping, presented here, yielded surprising information on lymphomagenesis.

Results

Stress responses in WT and Trp53ΔAS/ΔAS cells

We analyzed cellular stress responses in thymocytes, known to undergo a p53-dependent apoptosis upon irradiation (Lowe et al., 1993), and in primary fibroblasts, known to undergo a p53-dependent cell cycle arrest in response to various stresses; for example, DNA damage caused by irradiation or doxorubicin (Kastan et al., 1992), and the Nutlin-mediated inhibition of Mdm2, a negative regulator of p53 (Vassilev et al., 2004). We first compared thymocytes from irradiated wild-type (WT) and Trp53ΔAS/ΔAS mice. In WT thymocytes, isoforms with the AS C-terminus were five times less abundant than isoforms with the α C-terminus at the RNA level (Figure 1A), and in western blots the p53-AS protein appeared as a faint band running just ahead of, and often hard to separate from, the band specific for p53-α, the canonical full-length p53 (Figure 1B). In Trp53ΔAS/ΔAS thymocytes, mRNA levels for α isoforms were slightly decreased, if at all (Figure 1A), whereas p53-α protein levels appeared markedly decreased (Figure 1B), raising the possibility that p53-AS isoforms might contribute to p53-α abundance. Nevertheless, the transactivation of classical p53 target genes (Figure 1C) and apoptotic response (Figure 1D, Figure 1—figure supplement 1D) were not significantly altered by the loss of AS isoforms. Likewise, no significant difference was observed between WT and Trp53ΔAS/ΔAS fibroblasts in assays for cell cycle control (Figure 1E, Figure 1—figure supplement 1E), expression of well-known p53 target genes (Figure 1F, Figure 1—figure supplement 1F and G), proliferation under hyperoxic conditions (Figure 1G), or the growth of tumor xenografts (Figure 1H).

Figure 1 with 1 supplement see all
The loss of p53-AS isoforms does not alter cellular stress responses or survival to spontaneous tumors.

(A) mRNAs for p53-α and p53-AS isoforms from thymocytes of irradiated mice were quantified by RT-qPCR, and p53-α levels in Trp53+/+ mice were assigned a value of 1. Means ± SEM (n = 3). (B) Protein extracts from thymocytes of Trp53+/+ (WT) or Trp53ΔAS/ΔAS (ΔAS) irradiated mice were immunoblotted with p53 or actin antibodies. After normalization to actin, full-length (FL) p53-α levels in WT thymocytes were assigned a value of 1. (C) mRNA levels of p53 target genes in thymocytes, before or after γ-irradiation. Means ± SEM (n = 3). (D) Thymocyte apoptotic response to γ-irradiation. Means ± SEM (n = 6). (E) Cell cycle control in mouse embryonic fibroblasts (MEFs) after γ-irradiation. Asynchronous MEFs were exposed to 0–10 Gy γ-irradiation, and after 24 hr, cells were labeled with BrdU for 1 hr and analyzed by FACS. Means ± SEM from >3 independent experiments with at least two independent MEF clones per genotype. (F) mRNA levels of p53 target genes in MEFs untreated or treated with 0.5 μg/ml of clastogenic doxorubicin (Doxo) or 10 μM of the Mdm2 antagonist Nutlin. Means ± SEM from >3 experiments with ≥2 independent MEF clones. (G) MEFs proliferation under hyperoxic conditions. Cells were grown according to a 3T3 protocol. Each point is the mean from four independent MEF clones, the value for each clone resulting from triplicates. (H) Growth of tumor xenografts. E1A+Ras-expressing MEFs were injected into the flanks of nude mice and tumor volumes were determined after 1–25 days. Means ± SD (n = 4 per timepoint and genotype). (I) Tumor-free survival of Trp53+/+ and Trp53ΔAS/ΔAS mice (n = cohort size). (J) Incidence of the indicated tumor types, determined at death after macroscopic examination and histological analysis of Trp53+/+ (WT) and Trp53ΔAS/ΔAS (ΔAS) mice. In (A, D, E), ns = non-significant in Student’s t-test.

Lymphomagenesis in WT and Trp53ΔAS/ΔAS mice

We compared spontaneous tumor onset in WT and Trp53ΔAS/ΔAS littermates for over 2 years and observed no significant difference in tumor-free survival (Figure 1I). Because lymphoma is a common neoplasm in C57Bl/6J WT mice (Brayton et al., 2012) and our mouse cohorts resulted from >10 generations of backcrosses with C57Bl/6J mice, we searched for evidence of lymphoma in the lymph nodes and spleen, by macroscopic examination at autopsy and histological analyses. B-cell lymphomas were observed in about 30% of mice of either genotype (Figure 1J). In WT mice, a higher incidence of B-cell lymphomas was observed in females, in agreement with previous observations (Brayton et al., 2012). By contrast, no obvious sex-specific bias was observed for B-cell lymphomas in Trp53ΔAS/ΔAS mice (Figure 1J), raising the possibility that the loss of p53-AS isoforms affected B-cell lymphomagenesis. However, the numbers of lymphoma-bearing mice were too small to be conclusive.

We next used Eμ-Myc transgenic mice, prone to highly penetrant B-cell lymphomas (Adams et al., 1985). Trp53+/+ Eμ-Myc and Trp53ΔAS/ΔAS Eμ-Myc mice developed B-cell lymphomas (Figure 2—figure supplement 1A) with similar survival curves when sexes were not considered (Figure 2—figure supplement 1B). Importantly, however, death was accelerated, and tumor lymph nodes were larger, in Trp53ΔAS/ΔAS Eμ-Myc males compared to their Trp53+/+ Eμ-Myc male counterparts, whereas no difference in lymphomagenesis was noticeable between Trp53ΔAS/ΔAS Eμ-Myc and Trp53+/+ Eμ-Myc female mice (Figure 2A and B). Our data (Figure 2A and B, Figure 2—figure supplement 1C), together with the fact that B-cell lymphomas occur with a higher incidence in WT C57Bl/6J female mice (Brayton et al., 2012), suggested that Trp53+/+ Eμ-Myc male mice are more refractory to B-cell lymphomas, and that p53-AS isoforms might confer this male-specific protection against lymphomagenesis.

Figure 2 with 1 supplement see all
Male-specific acceleration of Myc-induced B-cell lymphomagenesis in mice lacking p53-AS isoforms.

(A) Tumor-free survival of Trp53+/+ Eμ-Myc and Trp53ΔAS/ΔAS Eμ-Myc mice, classified according to sex (n = cohort size). (B) Tumor volumes upon dissection of Trp53+/+ Eμ-Myc and Trp53ΔAS/ΔAS Eμ-Myc mice, classified according to sex. Means ± SEM from 100 lymph nodes from Trp53+/+ Eμ-Myc males, 96 from Trp53+/+ Eμ-Myc females, 148 from Trp53ΔAS/ΔAS Eμ-Myc males, and 124 from Trp53ΔAS/ΔAS Eμ-Myc females. (C) Myc mRNA and protein levels in lymph node tumors. (D) Levels of p53-α and p53-AS transcripts and p53 protein levels in lymph node tumors. (E) Transcript levels of the indicated p53 target genes. Means ± SEM (n = 6 per genotype). (F, G) Apoptosis (F) and cell proliferation (G) in tumor lymph nodes from Eμ-Myc males were determined by immunohistochemistry with antibodies against cleaved caspase-3 and ki67, respectively. Positive cells were counted and normalized to the analyzed areas. Means ± SEM (n = 6 mice per assay and genotype). In (F, G) scale bars = 50 μm. Statistical analyses with Mantel–Cox (A) and Student’s t (B–G) tests. ***p<0.001, *p<0.05, ns: nonsignificant.

Cause for accelerated lymphomagenesis in Trp53ΔAS/ΔAS Eμ-Myc males

We next aimed to determine the mechanisms underlying the accelerated lymphomagenesis in Trp53ΔAS/ΔAS Eμ-Myc males. Inactivating p53 mutations were not more frequent in tumors from Trp53ΔAS/ΔAS Eμ-Myc males than in those from Trp53+/+ Eμ-Myc males, ruling out additional mutations at the Trp53 locus as potential causes for accelerated lymphomagenesis in Trp53ΔAS/ΔAS Eμ-Myc males (Figure 2—figure supplement 1D and E). We next analyzed a subset of tumors with no detectable Trp53 mutation in males of both genotypes and found that Myc was expressed at similar RNA and protein levels in all tumors (Figure 2C). No difference in p53-α mRNA levels was observed in tumors from both genotypes, although a decrease at the protein level was detected in most tumors from Trp53ΔAS/ΔAS Eμ-Myc males (Figure 2D). Nevertheless, similar transcript levels for classical p53 target genes were observed in tumor cells of both genotypes (Figure 2E). To test whether a higher tumor volume in Trp53ΔAS/ΔAS Eμ-myc males might result from lower apoptosis and/or higher cell proliferation, we next analyzed tumors by immunohistochemistry with antibodies against cleaved caspase-3 or ki-67, respectively. Similar apoptotic and proliferation indexes were observed for both genotypes (Figure 2F and G). In sum, classical assays for p53 activity in tumors failed to account for differences between the two male genotypes.

The speed of lymphomagenesis in Eμ-Myc mice correlates with the extent of B-cell expansion in the first stages of B-cell differentiation (Langdon et al., 1986) and p53 was proposed to control the pool of early B cells (Slatter et al., 2010). Therefore, we determined the levels of the early pre-B/ immature B cells in 6-week-old mice, before any sign of tumor. We analyzed the spleen, a preferential site of B-cell expansion (Langdon et al., 1986) with a relatively high AS/α isoform ratio (Figure 1—figure supplement 1A). Flow cytometry with a combination of markers was used to discriminate the pre-B, immature, transitional, and mature B subpopulations. As expected (Langdon et al., 1986), we observed high numbers of pre-B and immature B cells in Eμ-Myc mice. In males, pre-B and immature B cells were more abundant in Trp53ΔAS/ΔAS Eμ-Myc animals, while no difference was observed for transitional and mature B cells (Figure 3A, Figure 3—figure supplement 1A). By contrast, in the spleen of Trp53+/+ and Trp53ΔAS/ΔAS 6-week-old male mice without the Eμ-Myc transgene, most B cells were mature B cells (Figure 3—figure supplement 1B). In Eμ-Myc females, the numbers of pre-B and immature B cells were similar between genotypes, as were the numbers of transitional and mature B cells (Figure 3A). Interestingly, Trp53+/+ Eμ-myc males, which develop lymphomas less rapidly, exhibited the lowest number of immature B cells (Figure 3A), suggesting a direct correlation between the level of immature B cell expansion and the speed of lymphomagenesis. Together, these data suggested that p53-AS isoforms may not be required to control the pool of early B cells under normal conditions, but that in an Eμ-Myc context they would limit the expansion of pre-tumor early B cells, specifically in males.

Figure 3 with 1 supplement see all
The loss of p53-AS isoforms affects Ackr4 expression in Eμ-Myc male mice.

(A) B-cell subpopulations in spleens of 6-week-old Trp53+/+ Eμ-Myc and Trp53ΔAS/ΔAS Eμ-Myc mice. Means ± SEM (n = 6 per genotype). (B, C) RNAseq analysis of spleens from Trp53+/+ Eμ-Myc (n = 3) and Trp53ΔAS/ΔAS Eμ-Myc (n = 4) 4–6-week-old male mice. Volcano plot (B), with differentially expressed genes (DEGs) in red. Unsupervised clustering heatmap plot (C), with DEGs ranked according to mean fold changes, and protein-coding genes in bold. (D) RT-qPCR analysis of candidate DEGs from spleens of Trp53+/+ (p53WT) Eμ-Myc males and Trp53ΔAS/ΔAS (p53ΔAS) Eμ-Myc males. Means ± SEM (n = 3–4 per genotype). (E) RT-qPCR analysis of indicated DEGs from spleens of 4–6-week-old Trp53+/+ Eμ-Myc males, Trp53ΔAS/ΔAS Eμ-Myc males, Trp53+/+ Eμ-Myc females, and Trp53ΔAS/ΔAS Eμ-Myc females. Means ± SEM (n = 3–4 per sex and genotype). (F) Ackr4 is transactivated by p53 in response to stress. Ackr4 mRNAs in untreated or doxorubicin-treated WT and Trp53-/- mouse embryonic fibroblasts (MEFs). Data from 2 to 3 MEFs per genotype (Younger et al., 2015). (G) A putative p53 response element in Ackr4 intron 1. Top: map of the Ackr4 gene. (boxes: exons [brown box: translated region]; black line: intron 1); middle: p53 ChIP in doxorubicin-treated MEFs according to ChIP-Atlas (SRX270554) (Oki et al., 2018); bottom: p53 Response Element (p53RE) consensus sequence (R = G or A, W = A or T, Y = C or T), the putative p53RE and its mutated counterpart. (H) Luciferase assays of the candidate p53RE. A 1.5 kb fragment containing the WT or mutant p53RE was cloned upstream a luciferase reporter, then transfected into Trp53-/- MEFs together with an expression plasmid for full length p53 (FL), p53-AS or the DNA-binding mutant p53R270H (RH). Means ± SEM (n = 4–6). (I) In MEFs, p53 activation leads to an increased Ackr4 expression attenuated by estradiol. Ackr4 and Cdkn1a mRNAs were quantified by RT-qPCR from Trp53+/+ and Trp53ΔAS/ΔAS MEFs, untreated or treated with 10 μM Nutlin and/or 5 μg/ml 17-β estradiol (E2). Means ± SEM from four independent experiments. (J) Evidence for Myc binding at the Mt2 promoter in B cells. ChIP-Atlas reports Myc binding to Mt2 promoter sequences in primary B cells from the lymph nodes of Eμ-Myc mice (SRX353785, SRX353783) and in Eμ-Myc-induced lymphoma cells (SRX522383). Chr: chromosome. (K) Gene set enrichment analysis (GSEA). GSEA, performed in Trp53+/+ Eμ-Myc (WT_Myc) and Trp53ΔAS/ΔAS Eμ-Myc (p53DAS_Myc) male splenic cells, indicated an enrichment of hallmark Myc targets in Trp53ΔAS/ΔAS Eμ-Myc cells. In (A, D, E, F, H, I) **p<0.001, **p<0.01, *p<0.05, °p≤0.057, ns: nonsignificant in Student’s t or Mann–Whitney tests.

Transcriptomes from Trp53+/+ Eμ-Myc and Trp53ΔAS/ΔAS Eμ-Myc male spleens

We next performed bulk RNA-seq and differential expression analyses comparing the spleens from 4- to 6-week-old Trp53ΔAS/ΔAS Eμ-Myc males to spleens from age-matched Trp53+/+ Eμ-Myc males. This revealed a limited number of significantly deregulated genes (Figure 3B), including 13 protein-coding genes and 11 pseudogenes (Figure 3C). Out of the 13 protein-coding genes, we focused on the 10 genes not encoding an immunoglobulin and analyzed the same samples by RT-qPCR (Figure 3D). For 6 of the 10 genes, expression levels were too low to be quantified (Tcstv3), or differences in expression were not statistically significant (Masp2, Akr1c19, Cd300lh, Il5ra, Slc26a1). Of note, Il5ra belonged to this group, although it is regulated by p53 (Zhu et al., 2022), which illustrates the difficulty to analyze subtle effects in our experiments. Taking this into account, we considered as potentially interesting the four remaining genes, exhibiting differences in mRNA levels with statistical significance (p<0.05) or borderline statistical significance (p=0.057): Ackr4, less expressed in Trp53ΔAS/ΔAS Eμ-Myc males, and Fam132a, Mt2, and Prss50, with an increased expression in Trp53ΔAS/ΔAS Eμ-Myc males (Figure 3D). Importantly, survival curves indicated that the Myc-induced lethality was delayed in Trp53+/+ Eμ-Myc males compared to Trp53ΔAS/ΔAS Eμ-Myc males, Trp53+/+ Eμ-Myc females, and Trp53ΔAS/ΔAS Eμ-Myc females (Figure 2—figure supplement 1C). Thus, we quantified transcripts for Ackr4, Fam132a, Mt2, and Prss50 in the spleen of 4–6-week-old Trp53+/+ Eμ-Myc and Trp53ΔAS/ΔAS Eμ-Myc females, then compared mRNA levels in Trp53+/+ Eμ-Myc males versus the three other groups. Significantly higher expression of Ackr4 and lower expression of Mt2 were found in Trp53+/+ Eμ-Myc males (Figure 3E).

Ackr4 (also known as Ccrl1) encodes the atypical chemokine receptor 4, a decoy receptor promoting the degradation of chemokines that modulate cancer cell proliferation and metastasis (Chow and Luster, 2014; Marcuzzi et al., 2018; Müller et al., 2001). Our data suggested that Ackr4 might be a gene transactivated by p53-α and/or p53-AS isoforms. Consistent with this, by extracting data from a transcriptome-wide study in MEFs (Younger et al., 2015), we found evidence for a p53-dependent transactivation of Ackr4 in response to doxorubicin (Figure 3F). Furthermore, ChIP-Atlas, the database of chromatin immunoprecipitation experiments (Oki et al., 2018), indicated p53 binding to sequences within the intron 1 of Ackr4 in doxorubicin-treated MEFs, and we identified a candidate p53 responsive element in this intron (Figure 3G). We next used luciferase assays to show that this p53 responsive element can be bound and regulated by both p53-α and p53-AS (Figure 3G and H, Figure 3—figure supplement 1C). Together, these data show that Ackr4 is indeed a p53 target gene, although RNAseq data indicated that it is expressed at much lower levels than classical p53 targets like Cdkn1a or Mdm2 in the splenic cells of Eμ-Myc mice (Supplementary file 1). Furthermore, Ackr4 was shown to be regulated by Foxl2 and estrogen signaling in ovarian cells (Georges et al., 2014) and 17-β estradiol was recently found to regulate ACKR4 expression in meniscal cells from both sexes, albeit differentially (Knewtson et al., 2020). Accordingly, we observed, in both WT and Trp53ΔAS/ΔAS MEFs, that p53 activation with the Mdm2 antagonist Nutlin led to the transactivation of Ackr4, but that a concomitant treatment with 17-β estradiol markedly decreased, or completely abrogated, Ackr4 transactivation (Figure 3I). By contrast, Cdkn1a was efficiently transactivated under both conditions in mutant and WT cells (Figure 3I). These data indicate that Ackr4 is a p53 target gene whose p53-mediated transactivation can be inhibited by estrogens.

The Mt2 gene, encoding the potentially oncogenic metallothionein-2 (Si and Lang, 2018), was less expressed in Trp53+/+ Eμ-Myc male pre-tumoral splenic cells, which raised the possibility of its direct or indirect repression by p53, potentially through the binding of p53 or the DREAM complex at its promoter (Engeland, 2018; Peuget and Selivanova, 2021). However, ChIP-Atlas reported no binding of these proteins at the Mt2 promoter. Alternatively, evidence that Myc may impact on Mt2 expression was obtained previously (Qin et al., 2021), and ChIP-Atlas reported Myc binding at the Mt2 promoter in primary B cells from lymph nodes of Eμ-Myc mice as well as Eμ-Myc-induced lymphoma cells (Figure 3J). This may suggest that the lower expression of Mt2 in pre-tumoral splenic cells from Trp53+/+ Eμ-Myc males (Figure 3E) might result from a subtle difference in Myc signaling. Consistent with this, the comparison of transcriptomes of splenic cells from Trp53ΔAS/ΔAS Eμ-Myc males and Trp53+/+ Eμ-Myc males, when analyzed by gene set enrichment analysis (Subramanian et al., 2005), revealed an enrichment of hallmark Myc target genes in Trp53ΔAS/ΔAS Eμ-Myc male splenic cells (Figure 3K, Supplementary file 2).

Relevance of ACKR4 expression in Burkitt lymphomas

Murine and human alternatively spliced p53 isoforms exhibit structural differences (Marcel et al., 2011) and the ChIP-Atlas database (Oki et al., 2018) does not report p53 binding to the human ACKR4 intron 1. Nevertheless, we found that p53 activation in human cells also led to an increased ACKR4 expression abrogated by 17-β estradiol, whereas 17-β estradiol had no significant effect on the p53-mediated transactivation of CDKN1A (Figure 4A). This led us to investigate the potential relevance of ACKR4 expression in human B-cell lymphomas. We analyzed public databases of B-cell lymphomas patients with clinical and gene expression information. We first analyzed #GSE4475 (Hummel et al., 2006), a dataset of mature aggressive B-cell lymphomas previously used to define Burkitt lymphoma-like specific transcriptomes, comprising 159 patients (91 men, 68 women) with clinical follow-up. Overall, ACKR4 gene expression was not significantly different between male and female patients (Figure 4B, left). However, average mRNA levels appeared higher in males when we considered the 30% patients of each sex with the highest ACKR4 expression (Figure 4B, right). Strikingly, when we compared the survival of 30% patients with the highest ACKR4 mRNA levels to the survival of the 30% patients with the lowest ACKR4 mRNA levels, high ACKR4 expression correlated with a better prognosis in men, but not in women (Figure 4C). By contrast, for MT2A, the human homolog of Mt2, differences in mRNA levels did not correlate with significant differences in survival for either sex (Figure 4—figure supplement 1A).

Figure 4 with 4 supplements see all
ACKR4 is a male-specific prognostic factor in Burkitt lymphoma.

(A) In human cells, p53 activation leads to an increased ACKR4 expression abrogated by estradiol. ACKR4 and CDKN1A mRNAs were quantified by RT-qPCR from p53-proficient (MRC5) and p53-deficient (MRC5-SV40) human fibroblasts, untreated or treated with Nutlin and/or estradiol (E2). Means ± SEM from four independent experiments. (B, C) Analysis of lymphoma dataset #GSE4475. ACKR4 gene expression was plotted for all lymphoma patients with clinical follow-up (91 men [M], 68 women [W]), classified according to sex (B, left). Gene expression (B, right) or survival curves (C) were plotted for the 30% patients (27 men, 20 women) with the highest or lowest ACKR4 expression, classified according to sex. (D, E) Analysis of Burkitt lymphoma-specific dataset #phs00235. ACKR4 gene expression was plotted for all patients with a Burkitt lymphoma diagnosed at age 0–17 (48 males, 29 females), classified according to sex (D, left). Gene expression (D, right) or survival curves (E) were plotted for the 30% patients (15 men, 9 women) with the highest or lowest ACKR4 expression, classified according to sex. (F) The knockout of ACKR4 in Burkitt lymphoma Raji cells increases their CCL21-guided migration. Chemotaxis was assayed by using Boyden chambers with bare polycarbonate membranes as previously described (Calpe et al., 2011). Equal number of cells were deposited on the membrane of a transwell insert, then migration was determined by counting cells in the lower compartment, after 15 hr of culture with or without CCL21 added to the lower chamber. Statistical analyses by Student’s t or Mann–Whitney tests (A, B, D, F) and Mantel–Cox (C, E) test. ***p<0.001, **p<0.01, *p<0.01, °p=0.054, ns: nonsignificant.

We also analyzed dataset #GSE181063 (Lacy et al., 2020), comprising mostly diffuse large B-cell lymphomas (DLBCL; 613 men, 536 women) and a few Burkitt lymphomas (65 men, 18 women). We found no difference in survival curves of DLBCL patients with low versus high ACKR4 levels, neither in men nor in women. However, there was again an increased survival for Burkitt lymphoma male patients with high ACKR4 expression, but not for women (Figure 4—figure supplement 1B). Next, we analyzed #phs000235 (Morin et al., 2011), a Burkitt lymphoma-specific dataset (65 men, 37 women) comprising mostly patients diagnosed at 0–17 years of age, hence providing cohorts homogeneous for both tumor type and age of onset. Again, ACKR4 was expressed at higher levels in a subset of male patients (Figure 4D) and high ACKR4 expression correlated with a better prognosis only in males (Figure 4E). Finally, we analyzed dataset #GSE136337 (Danziger et al., 2020), comprising data from the malignant plasma cells of patients with multiple myeloma (260 men, 166 women), and found that ACKR4 is not a prognostic factor for this cancer type (Figure 4—figure supplement 1C). Altogether, these analyses led us to conclude that, as in Eμ-Myc mice, ACKR4 is a male-specific positive prognostic factor in Burkitt lymphoma, the archetype of MYC-driven B-cell lymphomas.

We next considered the possibility that ACKR4 might be more than a biomarker, if it acts as a suppressor of MYC-driven B-cell lymphomagenesis. In support of this hypothesis, ACKR4 scavenges the chemokine CCL21, a ligand of the chemokine receptor CCR7 (Bastow et al., 2021; Ulvmar et al., 2014), and the ACKR4-mediated sequestration of CCL21 may impair the CCR7 signaling cascade, which might lead to decreased MYC activity (Shi et al., 2015). Consistent with this, a Ccr7 deficiency was shown to delay the lymphomagenesis induced by Eμ-Myc in mice (Rehm et al., 2011). In addition, Ccr7 is required for lymphoma cell lodging to secondary lymphoid organs (Rehm et al., 2011) and ACKR4 expression was inversely correlated with the metastasis capacity of different types of cancer cells (Shi et al., 2015; Zhu et al., 2014). These data suggested that ACKR4 might regulate the behavior of Burkitt lymphoma cells. To test this hypothesis, we designed a CRISPR-Cas9 approach to perform the knockout of ACKR4 in Raji cells. Raji is a Burkitt lymphoma cell line isolated from a 11-year-old male patient (Pulvertaft, 1964). Although p53 is mutated in Raji cells (Duthu et al., 1992), these cells overexpress MYC (Hamlyn and Rabbitts, 1983) and express both ACKR4 and CCR7 (Ferreira et al., 2014), suggesting that they might be suitable to evaluate the impact of ACKR4 on Burkitt lymphoma cell behavior, particularly their migratory capacities. Raji cells were transfected with a vector expressing Cas9, a puromycin resistance gene and either of two guide RNAs targeting ACKR4 (or no guide RNA for control), then puromycin-resistant cells were selected and recovered as cellular pools or diluted to isolate cellular clones (Figure 4—figure supplement 2). Both guide RNAs targeted sequences in ACKR4 mapping upstream an encoded DRY motif essential for signal transduction (Watts et al., 2013), so that Cas9-induced DNA breaks would generate knockout alleles. With this strategy, we obtained puromycin-resistant clones from Raji cells, two of which were verified to be ACKR4 KO clones by DNA sequencing (Figure 4—figure supplement 3). We compared the proliferation and migration capacities of Raji control cells (transfected with the Cas9 expression vector without guide RNAs) and the two independent ACKR4 KO Raji clones identified, cultured for 15 hr in medium supplemented or not with the CCL21 chemokine. Under these conditions, Raji cells of all genotypes appeared to proliferate similarly (Figure 4—figure supplement 4). Strikingly however, the KO of ACKR4 led to a fourfold increase in CCL21-mediated cell migration, consistent with the hypothesis that ACKR4 may hinder MYC-driven B-cell lymphomagenesis (Figure 4F).

Discussion

Here, we analyzed a mouse model with a specific deletion of the AS exon of the Trp53 gene. Despite a subtle phenotype, this model revealed that a male-specific protective effect against Eμ-Myc-induced B-cell lymphomas is lost in the absence of p53-AS isoforms. Trp53ΔAS/ΔAS males also appeared more prone to develop spontaneous lymphomas, suggesting that the sex-specific protective effect conferred by p53-AS isoforms might not be restricted to the Eμ-Myc model.

Our transcriptomic data from splenic cells of Trp53+/+ Eμ-Myc and Trp53ΔAS/ΔAS Eμ-Myc males disclosed very few differentially expressed genes and highlighted Ackr4 as a male-specific positive prognostic factor in Eμ-Myc-induced lymphomas. Mechanistically, we identified Ackr4, expressed at low levels in splenic cells, as a p53 target gene that may be transactivated by p53-α and/or p53-AS according to luciferase assays. That Ackr4 might be regulated by both types of p53 isoforms was expected because their DNA binding domains are identical. In fact, if one considers that p53-α isoforms appear more abundant than p53-AS isoforms in wild-type cells, and that the loss of p53-AS isoforms correlated with a decrease in p53-α levels in the thymocytes and tumor lymph nodes of mutant mice, then it seems likely that the reduced transactivation of Ackr4 in the splenic cells of Trp53ΔAS/ΔAS Eμ-Myc males could mainly result from decreased p53-α levels, rather than the loss of p53-AS isoforms per se. In addition, we observed that 17-β estradiol can inhibit the p53-mediated transactivation of Ackr4. Together, our data suggest that Ackr4 may be regulated by p53-α, p53-AS, and estrogens, likely accounting for sex-specific and p53-status-dependent differences in gene expression.

Our analyses reveal that in both mice and humans Ackr4/ACKR4 is a male-specific prognostic factor in Burkitt-like lymphomas. Furthermore, several lines of evidence suggest that Ackr4 might act as a tumor suppressor of Myc-driven B-cell lymphomas. As mentioned before, the ACKR4-mediated sequestration of CCL21 may impair the CCR7 signaling cascade, which might lead to decreased MYC activity (Shi et al., 2015). In Trp53+/+ Eμ-Myc male splenic cells, the observed lower expression of Mt2, known to be regulated by Myc (Qin et al., 2021), and of many genes that are hallmark Myc targets (as revealed by GSEA), appear consistent with this hypothesis. In addition, Ccr7 is required for lymphoma cell lodging to secondary lymphoid organs (Rehm et al., 2011) and ACKR4 expression was inversely correlated with the metastasis capacity of different types of cancer cells (Shi et al., 2015; Zhu et al., 2014). Consistent with this, we found that the KO of ACKR4 in Raji Burkitt lymphoma cells led to a dramatic increase in CCL21-guided cell migration. Finally, Ackr4 regulates B cell differentiation (Kara et al., 2018), which raises the possibility that an alteration of the p53-Ackr4 pathway in Trp53ΔAS/ΔAS Eμ-Myc male splenic cells might contribute to increase the pools of pre-B and immature B cells that may be prone to lymphomagenesis. In sum, a decrease in Ackr4 expression might promote B-cell lymphomagenesis through several non-exclusive mechanisms. Importantly, ACKR4 was previously found to inhibit the growth and metastasis of breast, cervical, colorectal, hepatocellular, and nasopharyngeal cancer cells (Feng et al., 2009; Hou et al., 2013; Ju et al., 2019; Shi et al., 2015; Zhu et al., 2014), although no report mentioned any sex-specific bias for cancers occurring in both sexes. Our data provide evidence that sex-specific differences in Ackr4 expression may have prognostic value. This suggests that measuring ACKR4 gene expression in male patients with Burkitt lymphoma could be useful to identify the patients at higher risk, for whom specific therapeutic regimens might be required.

Interestingly, our data suggested that Mt2 might be a male-specific negative prognostic factor in murine Eμ-Myc-induced lymphomas, but MT2A expression levels had no prognostic value in human lymphomas. A possible explanation for this discrepancy is suggested by the fact that Mt2 is regulated by Myc. A translocation leading to MYC overexpression drives oncogenesis in all Burkitt lymphomas, but half of them exhibit additional missense MYC mutations enhancing its tumorigenicity (Chakraborty et al., 2015). The transcriptional program of a WT and two lymphoma-associated Myc mutants were recently compared, and we noticed that one of the mutants led to an alteration in Mt2 expression (Mahani et al., 2021), which would abrogate any potential prognostic value.

Finally, a polymorphism in the MDM2 gene promoter provided evidence that sex-specific hormones may affect p53 signaling and tumorigenesis (Bond and Levine, 2007). More recently, a higher frequency of TP53 mutations in men, together with an increased vulnerability to alterations of X-linked genes encoding p53 regulators, was proposed to explain a higher cancer incidence and death in male patients (Haupt et al., 2019). Here, on the contrary, male mice and a subset of male patients were more efficiently protected against Burkitt-like lymphomas, which adds another layer of complexity to sex-specific differences in tumorigenesis. The p53 pathway thus underlies cancer sex-disparities through multiple mechanisms, which may notably include variations in p53 isoforms or Ackr4 expression.

Materials and methods

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Gene (Mus musculus)Trp53GenBankENSMUSG00000059552
Gene (M. musculus)Ackr4GenBankENSMUSG00000079355
Gene (M. musculus)Mt2GenBankENSMUSG00000031762
Gene (Homo sapiens)ACKR4GenBankENSMUSG00000129048
Strain, strain background (M. musculus, both sexes)Trp53ΔAS, C57Bl/6JSimeonova et al., 2013
Strain, strain background (M. musculus, both sexes)Eμ−Myc, C57Bl/6JJackson LabsB6.Cg-Tg(IghMyc)22Bri/J
Strain, strain background (M. musculus, females)CD-1 Nude
CD-1
Charles River LabsCrl:CD1-Foxn1nu
Strain, strain background (M. musculus, both sexes)C57Bl/6JCharles River Labs
Cell line (M. musculus, both sexes)WT, Trp53+/ΔAS, Trp53ΔAS/ΔAS, Trp53+/-, Trp53-/- fibroblastsThis paperPrimary fibroblasts prepared from E13.5 days embryos‘Cells and cell culture reagents’
Cell line (H. sapiens, male)MRC5Sigma-AldrichMRC5 PD19 (#05072101)
Cell line (H. sapiens, male)MRC5-SV40Sigma-AldrichMRC5-SV2 (#84100401)
Cell line (H. sapiens, female)HEK293TATCCCRL-3216
Cell line (H. sapiens, male)RajiATCCCCL-86
Cell line (H. sapiens, male)Raji ACKR4 KOThis paperACKR4 KO 4.14 & 5.2 Raji derivativesFigure 4—figure supplement 3
Transfected construct (Adenoviral E1A)pWZL-E1A12SAddgenepWZL hygro 12S E1A (#18748)
Transfected construct (human Ras)pBabe-Hrasv12AddgenepBabe-puro Ras v12 (# 1768)
Antibodyp53 (rabbit polyclonal)NovocastraLeica NCL-p53-CM5p1/2000
AntibodyMyc (mouse monoclonal)Santa Cruz9E-10 sc401/1000
Antibodyp21 (mouse monoclonal)Santa CruzF-5 sc62461/200
AntibodyActin (mouse monoclonal)Santa CruzActin-HRP sc477781/5000
AntibodyCD45R/B220 APC (rat, monoclonal)BD BiosciencesAnti-mouse CD45R/B220 APC (#561880)1/200
AntibodyIgD (rat, monoclonal)BD BiosciencesAnti-mouse IgD BV 605 (#563003)1/100
AntibodyCD43 (rat, monoclonal)BD BiosciencesAnti-mouse CD43 FITC (#561856)1/200
AntibodyIgM (rat, monoclonal)BD BiosciencesAnti-mouse IgM PE (#562033)1/50
Recombinant DNA reagentpSpCas9(BB)–2A-PuroAddgenePX459 (#48139)
Sequence-based reagentTrp53α-FThis paperqPCR primerAAAGGATGCCCATGCTACAGA; Figure 1—figure supplement 1
Sequence-based reagentTrp53α-RThis paperqPCR primerTCTTGGTCTTCAGGTAGCTGGAG; Figure 1—figure supplement 1
Sequence-based reagentTrp53AS-FThis paperqPCR primerAAAGGATGCCCATGCTACAGA; Figure 1—figure supplement 1
Sequence-based reagentTrp53AS-RThis paperqPCR primerTGAAGTGATGGGAGCTAGCAGTT; Figure 1—figure supplement 1
Sequence-based reagentAckr4-FThis paperqPCR primerGCACCTCTCCCAGCTTAAACA; Figure 3
Sequence-based reagentAckr4-RThis paperqPCR primerAATAGTATTCCGCTGACTGGTTCAG; Figure 3
Sequence-based reagentACKR4-FThis paperqPCR primerACTGCTCCTCTCTGCCGACTAC; Figure 4
Sequence-based reagentACKR4-RThis paperqPCR primerGCCATTCATTTCATTTTCCTCAT; Figure 4
Sequence-based reagentACKR4-g4This paperGuide for CRISPR #4TGGTAGTGGCAATTTATGCC; Figure 4—figure supplement 3
Sequence-based reagentACKR4-g5This paperGuide for CRISPR #5GGGCTGTTAATGCAGTTCAT; Figure 4—figure supplement 3
Peptide, recombinant proteinCCL21Preprotech#300-35A
Peptide, recombinant proteinSuperscript IVInvitrogenTF #18090010
Commercial assay or kitNucleospin RNA IIMacherey-NagelFS #NZ74095520
Commercial assay or kitPower SYBR GreenApplied Biosystems# 4367659
Commercial assay or kitSupersignal West FemtoThermo Fisher# 34096
Commercial assay or kitAnnexinV-FITC apoptosis staining/ detection kitAbcam# Ab14085
Commercial assay or kitTruseq stranded Total RNAIllumina#20020596
Commercial assay or kitNucleofector Amaxa kit VLonza# VCA-1003
Chemical compound, drugDoxorubicinSigma-Aldrich# D1515
Chemical compound, drugEtoposideSigma-Aldrich# E1383
Chemical compound, drugNutlin 3aSigma-Aldrich# SML-0580
Chemical compound, drug17β-estradiolSigma-Aldrich# E2758
Software, algorithmFlowJoBeckton-DickinsonRRID:SCR_008520v 10.10
Software, algorithmfeatureCountsLiao et al., 2014
Software, algorithmDESeq2 R packageLove et al., 2014
Software, algorithmGSEA softwareSubramanian et al., 2005
Software, algorithmPWMScanAmbrosini et al., 2018
Software, algorithmCRISPORHaeussler et al., 2016
Software, algorithmPrismGraphPadRRID:SCR_002798v 5.0

Mice

Design and construction of the Trp53ΔAS mouse model were previously described (Simeonova et al., 2013). A minimum of 10 backcrosses with C57Bl/6J mice of both sexes (Charles River Laboratories) were performed before establishing the cohorts of Trp53+/+ and Trp53ΔAS/ΔAS littermate mice used in this study. Mouse genotyping with multiple primer sets confirmed >99% C57Bl/6J genetic background after 10 backcrosses (primer sequences available upon request). Cohorts of Trp53+/+ Eμ-Myc and Trp53ΔAS/ΔAS Eμ-Myc mice were established with identical parental origin of the Eμ-Myc transgene. For all experiments, mice housing and treatment were conducted according to the Institutional Animal Care and Use Committee of the Institut Curie (approved project #03769.02).

Cells and cell culture reagents

Request a detailed protocol

MEFs were isolated from 13.5 days post-coitum embryos and cultured in a 5% CO2 and 3% O2 incubator, in Dulbecco's Modified Eagle Medium (DMEM) GlutaMAX (Gibco), with 15% Fetal Bovine Serum (FBS) (PAN Biotech), 100 μM 2-mercaptoethanol (Millipore), 0.1 mM non-essential amino acids and penicillin/streptomycin (Gibco) for less than five passages, except for 3T3 experiments, performed in a 5% CO2 incubator for nine passages. Cells were treated for 24 hr with 0.5 μg/ml doxorubicin (Sigma-Aldrich), 15 μM etoposide (Sigma-Aldrich), 10 μM Nutlin 3a (Vassilev et al., 2004) (Sigma-Aldrich), and/or 5 μg/ml 17β-estradiol (Sigma-Aldrich). At least three independent experiments with at least two independent littermate MEF clones of each genotype and each sex were performed to measure DNA damage responses. For estradiol assays, four independent experiments with three independent MEF male clones of each genotype were performed. Human lung fibroblast MRC5 and its SV40-transformed derivatives were cultured in a 5% CO2 and 3% O2-regulated incubator in Minimum Essential Medium (MEM) medium without Phenol Red (Gibco), completed with 10% FBS, 2 mM l-glutamine (Gibco), 1 mM pyruvate, 0.1 mM non-essential amino acids, and penicillin/streptomycin, and treated for 24 hr with 10 μM Nutlin 3a and/or 5 μg/ml 17β-estradiol (Merck). Four independent experiments were performed. Burkitt lymphoma cells (Raji or Raji ACKR4 KO derivatives) were cultured in a 5% CO2 incubator in RPMI GlutaMAX (Gibco) supplemented with 10% FBS, 2 mM l-glutamine, 1 mM pyruvate, 0.1 mM non-essential amino acids, 0.45% glucose (Sigma-Aldrich), and penicillin/streptomycin. In proliferation and migration assays, Raji control cells and two independent ACKR4 KO derivatives were treated or not for 15 hr with 1 μg/ml CCL21 (PreproTech) then counted in triplicates. For all cell lines, cell culture supernatants were found negative for mycoplasma contamination.

Quantitative RT-PCR

Request a detailed protocol

Total RNAs were extracted using nucleospin RNA II (Macherey-Nagel), reverse-transcribed using superscript IV (Invitrogen), and real-time quantitative PCRs were performed on an ABI PRISM 7500 using Power SYBR Green (Applied Biosystems) as previously described (Simeonova et al., 2012). For quantification of p53 isoforms in healthy tissues, a forward primer in exon 10 and a reverse primer encompassing the boundary between exons 10 and 11 were used for p53-α amplification, whereas the same forward primer and a reverse primer located in exon AS were used for p53-AS amplification (see Supplementary file 3 for primer sequences). To determine the AS/α mRNA ratios, expression levels were compared with a standard curve generated by serial dilutions of a plasmid containing both p53-AS and p53-α cDNAs.

Western blots

Request a detailed protocol

Thymocytes were lysed in RIPA buffer (50 mM Tris–HCl pH 8, 150  mM NaCl, 5  mM EDTA, 0.5% deoxycholic acid, 0.1% SDS, 1% NP-40) with a cocktail of protease inhibitors (Roche) and 1  mM PMSF (Sigma). Whole-cell extracts were sonicated three times for 10 s and centrifuged at 13,000  r.p.m. for 30  min to remove cell debris. MEFs or B-cell lymphomas were lysed in Giordano’s buffer (50  mM Tris–HCl pH 7.4, 250  mM NaCl, 5  mM EDTA, 0.1% Triton X-100) with a cocktail of protease inhibitors (Roche) and 1 mM PMSF (Sigma). Protein lysate concentration was determined by bicinchoninic acid (BCA) assay (Thermo Scientific) and 30  μg of each lysate was fractionated by SDS–PAGE on a 4–12% polyacrylamide gel and transferred onto polyvinylidene difluoride (PVDF) membrane (Amersham). Membranes were incubated with antibodies against p53 (CM5, Novocastra), myc (9E-10, Santa Cruz), p21 (F-5, Santa Cruz), and actin (actin-HRP sc47778, Santa Cruz) and revealed with SuperSignal West femto detection reagent (Thermo Scientific).

Apoptosis assays

Request a detailed protocol

Six-week-old Trp53+/+ and Trp53ΔAS/ΔAS male mice were whole-body irradiated with 5 Gy of γ-irradiation. Mice were sacrificed 4 hr later and thymocytes were recovered, stained with AnnexinV-FITC Apoptosis detection kit (Abcam), then analyzed by flow cytometry using FlowJo.

Cell-cycle assays

Request a detailed protocol

Log-phase MEFs were irradiated at room temperature with a CS γ-irradiator at doses of 3 or 10 Gy, incubated for 24 hr, then pulse-labeled for 1 hr with 10 μM BrdU, fixed in 70% ethanol, double-stained with FITC anti BrdU and propidium iodide, and sorted by flow cytometry using a BD FACSort. Data were analyzed using FlowJo.

Oncogene-induced tumor xenografts

Request a detailed protocol

MEFs with the indicated genotypes were sequentially infected with pWZL-E1A12S and pBABE-Hrasv12 viruses as previously described (Toledo et al., 2006). In total, 5 × 106 E1A- and Ras- (E1A+Ras) expressing MEFs of each genotype were injected subcutaneously into the flanks of 7-week-old female athymic nude mice (at least four mice per genotype) and tumor volumes were determined 1, 8, 15, 21, and 25 days after injection. Importantly, populations of (E1A+Ras)-expressing cells were used to minimize potential differences in expression levels that could result from independent viral insertion sites.

Cell sorting of B-cell subpopulations

Request a detailed protocol

Splenic cells were recovered from 6-week-old asymptomatic mice and incubated with DAPI and the following antibodies: APC rat anti-mouse CD45R/B220, FITC rat anti-mouse CD43, PE rat anti-mouse IgM, and BV605 rat anti-mouse IgD (BD Pharmingen). First, the B220+ CD43 cells were selected by flow cytometry from DAPI-negative living cells, yielding subsequently four different B subpopulations based on IgM and IgD labeling: IgM-/IgD- preB lymphocytes, IgM low/IgD- immature B lymphocytes, IgM high/IgD- transitional B lymphocytes, and IgM+/IgD+ mature B lymphocytes.

RNA-seq analysis

Request a detailed protocol

Total RNA was extracted from the spleen of 4–6-week--old asymptomatic mice using nucleospin RNA II (Macherey-Nagel). The quality of RNA was checked with Bioanalyzer Agilent 2100 and RNAs with a RIN (RNA integrity number) > 6 were retained for further analysis. RNA was depleted from ribosomal RNA, then converted into cDNA libraries using a TruSeq Stranded Total Library preparation kit (Illumina). Paired-end sequencing was performed on an Illumina MiSeq platform. Reads were mapped to the mouse genome version GRCm38 and counted on gene annotation gencode.vM18 with featureCounts (Liao et al., 2014). Differentially expressed genes of C57Bl/6J genetic background with an adjusted p-value<0.05 were identified using the DESeq2 R package (Love et al., 2014). Gene set enrichment analysis was performed using the GSEA software with canonical pathway gene sets from the Mouse Molecular Signature Database (Subramanian et al., 2005).

Luciferase assays

Request a detailed protocol

The candidate p53 responsive element (p53 RE) in the Ackr4 promoter was identified using the JASPAR database of binding profiles (Fornes et al., 2020) with the position weight matrix scanner PWMscan (Ambrosini et al., 2018). A 1.5 kb fragment from Ackr4 intron 1, containing a WT or mutant p53 RE at its center, was cloned upstream an SV40 minimal promoter and a luciferase reporter gene in the backbone of a PGL3 plasmid (Promega). We used lipofectamine 2000 to transfect Trp53-/- MEFs with 2 μg of either luciferase expression vector, 2 μg of an expression vector for p53WT, p53AS or the DNA-binding mutant p53R270H, and 30 ng of a renilla luciferase expression plasmid (pGL4.73, Promega) for normalization. Transfected cells were incubated for 24 hr, then trypsinized, resuspended in 75 μl culture medium with 7.5% FBS, and transferred into a well of an optical 96-well plate (Nunc). The dual-glo luciferase assay system (Promega) was used according to the manufacturer’s protocol to lyse the cells and read firefly and renilla luciferase signals. Results were normalized, then the average luciferase activity in cells transfected with the WT p53RE luciferase reporter and the p53R270H expression plasmid was assigned a value of 1.

Generation of ACKR4 knockout Burkitt lymphoma cells

Request a detailed protocol

Six guide RNAs (gRNAs) were designed to target the human ACKR4 gene using the web tool CRISPOR (Haeussler et al., 2016). For each gRNA, two reverse complementary oligonucleotides were designed (with added backbone sequences, including a 5′ G nucleotide for gRNA sequences without a 5′ G to improve transcription efficiency from a U6 promoter), annealed, and cloned in the Cas9 expression vector pSpCas9(BB)–2A-Puro (PX459) (Addgene). Preliminary tests for efficiency to induce cleavage in the targeted DNA regions were performed in HEK293T cells, which led to select the gRNAs #4 (5′-TGGTAGTGGCAATTTATGCC-3′) and #5 (5′-GGGCTGTTAATGCAGTTCAT-3′) for further experiments. Raji Burkitt lymphoma cells were transfected with a PX459 vector expressing gRNA #4 or #5 (or no gRNA for control) using Nucleofector Amaxa kit V (Lonza), and selection was carried out with 1 μg/ml puromycin for 3 weeks to obtain stably transfected cells. Clones were obtained by seeding cells at low density (1 cell/5 wells) in 96-well plates, expanded, then subdivided into two parts, half for storage in liquid nitrogen and half for ACKR4 genotyping by Sanger DNA sequencing (see Supplementary file 3 for primer sequences). For each cell clone, the ACKR4 target DNA regions were amplified by PCR, PCR products were cloned in a PGL3 plasmid, and eight transformed bacteria colonies were sequenced (Figure 4—figure supplement 2). This led to identify two homozygous ACKR4 KO clones, one (clone 4.14) obtained using gRNA #4 and one (clone 5.2) using gRNA #5 (Figure 4—figure supplement 3).

Cell migration assay

Request a detailed protocol

Cell migration assays toward the chemokine CCL21 were performed with Boyden chambers essentially as described (Calpe et al., 2011) by measuring transwell migration across bare polycarbonate membranes with a pore size of 5 μm (Corning). A total of 100 μl of culture medium (RPMI GlutaMAX, 10% FBS, 2 mM l-glutamine, 1 mM pyruvate, 0.1 mM non-essential amino acids, 0.45% glucose, and penicillin/streptomycin) containing 5 × 105 cells was added to a 6.5 mm diameter transwell insert, and 600 μl of culture medium with or without 1 μg/ml CCL21 were added to the lower compartment. After 15 hr at 37°C in 5% CO2, the number of migrated cells in the lower chamber was determined using a Coulter Counter.

Statistical analyses

Request a detailed protocol

Student’s unpaired t-tests were used in most figures to analyze the differences between WT vs. ΔAS values. Log-rank (Mantel–Cox) tests were used to analyze Kaplan–Meier tumor-free survival curves. Analyses were performed using GraphPad Prism 5, and values of p<0.05 were considered significant.

Data availability

RNA sequencing data have been deposited in the Gene Expression Omnibus (GEO) under the accession code GSE209708. All other data are available within the article and its supplementary information.

The following data sets were generated
    1. Fajac A
    2. Bardot B
    3. Toledo F
    4. Gabriel M
    (2024) NCBI Gene Expression Omnibus
    ID GSE209708. Splenocyte mRNA profiles of 4-6 weeks-old p53+/+ Em-Myc and p53DAS/DAS Em-Myc mice.
The following previously published data sets were used
    1. Care MA
    2. Barrens SL
    (2021) NCBI Gene Expression Omnibus
    ID GSE181063. Whole genome expression profiling based on paraffin embedded tissue of a large DLBCL cohort.

References

    1. Duthu A
    2. Debuire B
    3. Romano J
    4. Ehrhart JC
    5. Fiscella M
    6. May E
    7. Appella E
    8. May P
    (1992)
    p53 mutations in Raji cells: Characterization and localization relative to other Burkitt’s lymphomas
    Oncogene 7:2161–2167.
    1. Flaman JM
    2. Waridel F
    3. Estreicher A
    4. Vannier A
    5. Limacher JM
    6. Gilbert D
    7. Iggo R
    8. Frebourg T
    (1996)
    The human tumour suppressor gene p53 is alternatively spliced in normal cells
    Oncogene 12:813–818.

Article and author information

Author details

  1. Anne Fajac

    1. Genetics of Tumor Suppression, Institut Curie, Paris, France
    2. CNRS UMR3244, Paris, France
    3. Sorbonne University, Paris, France
    4. PSL Research University, Paris, France
    Contribution
    Conceptualization, Formal analysis, Investigation, Visualization, Writing – original draft, Writing – review and editing
    Competing interests
    No competing interests declared
  2. Iva Simeonova

    1. Genetics of Tumor Suppression, Institut Curie, Paris, France
    2. CNRS UMR3244, Paris, France
    3. Sorbonne University, Paris, France
    4. PSL Research University, Paris, France
    Contribution
    Investigation, Initiated the project
    Competing interests
    No competing interests declared
  3. Julia Leemput

    1. Genetics of Tumor Suppression, Institut Curie, Paris, France
    2. CNRS UMR3244, Paris, France
    3. Sorbonne University, Paris, France
    4. PSL Research University, Paris, France
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  4. Marc Gabriel

    1. CNRS UMR3244, Paris, France
    2. Sorbonne University, Paris, France
    3. PSL Research University, Paris, France
    4. Non Coding RNA, Epigenetic and Genome Fluidity, Institut Curie, Paris, France
    Contribution
    Formal analysis
    Competing interests
    No competing interests declared
  5. Aurélie Morin

    1. Genetics of Tumor Suppression, Institut Curie, Paris, France
    2. CNRS UMR3244, Paris, France
    3. Sorbonne University, Paris, France
    4. PSL Research University, Paris, France
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  6. Vincent Lejour

    1. Genetics of Tumor Suppression, Institut Curie, Paris, France
    2. CNRS UMR3244, Paris, France
    3. Sorbonne University, Paris, France
    4. PSL Research University, Paris, France
    Contribution
    Investigation
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-2797-1507
  7. Annaïg Hamon

    1. Genetics of Tumor Suppression, Institut Curie, Paris, France
    2. CNRS UMR3244, Paris, France
    3. Sorbonne University, Paris, France
    4. PSL Research University, Paris, France
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  8. Jeanne Rakotopare

    1. Genetics of Tumor Suppression, Institut Curie, Paris, France
    2. CNRS UMR3244, Paris, France
    3. Sorbonne University, Paris, France
    4. PSL Research University, Paris, France
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  9. Wilhelm Vaysse-Zinkhöfer

    1. Genetics of Tumor Suppression, Institut Curie, Paris, France
    2. CNRS UMR3244, Paris, France
    3. Sorbonne University, Paris, France
    4. PSL Research University, Paris, France
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  10. Eliana Eldawra

    1. Genetics of Tumor Suppression, Institut Curie, Paris, France
    2. CNRS UMR3244, Paris, France
    3. Sorbonne University, Paris, France
    4. PSL Research University, Paris, France
    Contribution
    Investigation
    Competing interests
    No competing interests declared
  11. Marina Pinskaya

    1. CNRS UMR3244, Paris, France
    2. Sorbonne University, Paris, France
    3. PSL Research University, Paris, France
    4. Non Coding RNA, Epigenetic and Genome Fluidity, Institut Curie, Paris, France
    Contribution
    Formal analysis
    Competing interests
    No competing interests declared
  12. Antonin Morillon

    1. CNRS UMR3244, Paris, France
    2. Sorbonne University, Paris, France
    3. PSL Research University, Paris, France
    4. Non Coding RNA, Epigenetic and Genome Fluidity, Institut Curie, Paris, France
    Contribution
    Formal analysis
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0575-5264
  13. Jean-Christophe Bourdon

    School of Medicine, Ninewells Hospital, University of Dundee, Dundee, United Kingdom
    Contribution
    Resources, Initiated the project
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4623-9386
  14. Boris Bardot

    1. Genetics of Tumor Suppression, Institut Curie, Paris, France
    2. CNRS UMR3244, Paris, France
    3. Sorbonne University, Paris, France
    4. PSL Research University, Paris, France
    Contribution
    Conceptualization, Formal analysis, Supervision, Investigation, Visualization, Writing – original draft, Writing – review and editing
    For correspondence
    boris.bardot@curie.fr
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-4976-9593
  15. Franck Toledo

    1. Genetics of Tumor Suppression, Institut Curie, Paris, France
    2. CNRS UMR3244, Paris, France
    3. Sorbonne University, Paris, France
    4. PSL Research University, Paris, France
    Contribution
    Conceptualization, Formal analysis, Supervision, Funding acquisition, Investigation, Visualization, Writing – original draft, Project administration, Writing – review and editing, Initiated the project
    For correspondence
    Franck.Toledo@curie.fr
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3798-4106

Funding

Fondation de France (Comité tumeurs)

  • Franck Toledo

Ligue Contre le Cancer (Labellisation)

  • Franck Toledo

Ligue Contre le Cancer (Comité Ile de France)

  • Franck Toledo

Fondation ARC pour la Recherche sur le Cancer

  • Franck Toledo

European Research Council (875532-Prostator-ERC-2019-PoC)

  • Antonin Morillon

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

This project was supported by grants from the Comité Tumeurs of Fondation de France (to FT), the Comité Ile-de-France and Comité national (Labellisation 2014–18) of the Ligue Nationale Contre le Cancer (to FT), and the Fondation ARC pour la recherche sur le Cancer (to FT). IS, JR, and EE were PhD fellows of the Ministère de la Recherche; JL and AMorin were postdoctoral fellows of Cancéropôle Ile-de-France and Institut National du Cancer. MG was paid by European the Research Council 875532-Prostator-ERC-2019-PoC attributed to AMoril; JCB was supported by Cancer Research-UK. We thank K Fernandes for his participation in making the ΔAS mutation, and members of the Institut Curie platforms: I Grandjean, H Gautier, C Daviaud, M Garcia, M Verlhac, A Fosse, and P Bureau (animal facility); S Baulande and S Lameiras (NGS); M Huerre, A Nicolas, and R Leclere (histopathology); Z Maciorowski, A Viguier, and S Grondin (flow cytometry).

Ethics

For all experiments, mice housing and treatment were conducted according to Institutional Animal Care and Use Committee of the Institut Curie (approved project #03769.02).

Version history

  1. Sent for peer review:
  2. Preprint posted:
  3. Reviewed Preprint version 1:
  4. Reviewed Preprint version 2:
  5. Version of Record published:
  6. Version of Record updated:

Cite all versions

You can cite all versions using the DOI https://doi.org/10.7554/eLife.92774. This DOI represents all versions, and will always resolve to the latest one.

Copyright

© 2024, Fajac et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 750
    views
  • 102
    downloads
  • 0
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Anne Fajac
  2. Iva Simeonova
  3. Julia Leemput
  4. Marc Gabriel
  5. Aurélie Morin
  6. Vincent Lejour
  7. Annaïg Hamon
  8. Jeanne Rakotopare
  9. Wilhelm Vaysse-Zinkhöfer
  10. Eliana Eldawra
  11. Marina Pinskaya
  12. Antonin Morillon
  13. Jean-Christophe Bourdon
  14. Boris Bardot
  15. Franck Toledo
(2024)
Mutant mice lacking alternatively spliced p53 isoforms unveil Ackr4 as a male-specific prognostic factor in Myc-driven B-cell lymphomas
eLife 13:RP92774.
https://doi.org/10.7554/eLife.92774.3

Share this article

https://doi.org/10.7554/eLife.92774

Further reading

    1. Cancer Biology
    2. Genetics and Genomics
    Joakim W Karlsson, Vasu R Sah ... Jonas A Nilsson
    Research Article

    Uveal melanoma (UM) is a rare melanoma originating in the eye’s uvea, with 50% of patients experiencing metastasis predominantly in the liver. In contrast to cutaneous melanoma, there is only a limited effectiveness of combined immune checkpoint therapies, and half of patients with uveal melanoma metastases succumb to disease within 2 years. This study aimed to provide a path toward enhancing immunotherapy efficacy by identifying and functionally validating tumor-reactive T cells in liver metastases of patients with UM. We employed single-cell RNA-seq of biopsies and tumor-infiltrating lymphocytes (TILs) to identify potential tumor-reactive T cells. Patient-derived xenograft (PDX) models of UM metastases were created from patients, and tumor sphere cultures were generated from these models for co-culture with autologous or MART1-specific HLA-matched allogenic TILs. Activated T cells were subjected to TCR-seq, and the TCRs were matched to those found in single-cell sequencing data from biopsies, expanded TILs, and in livers or spleens of PDX models injected with TILs. Our findings revealed that tumor-reactive T cells resided not only among activated and exhausted subsets of T cells, but also in a subset of cytotoxic effector cells. In conclusion, combining single-cell sequencing and functional analysis provides valuable insights into which T cells in UM may be useful for cell therapy amplification and marker selection.

    1. Cancer Biology
    2. Neuroscience
    Jeffrey Barr, Austin Walz ... Paola D Vermeer
    Research Article

    Cancer patients often experience changes in mental health, prompting an exploration into whether nerves infiltrating tumors contribute to these alterations by impacting brain functions. Using a mouse model for head and neck cancer and neuronal tracing, we show that tumor-infiltrating nerves connect to distinct brain areas. The activation of this neuronal circuitry altered behaviors (decreased nest-building, increased latency to eat a cookie, and reduced wheel running). Tumor-infiltrating nociceptor neurons exhibited heightened calcium activity and brain regions receiving these neural projections showed elevated Fos as well as increased calcium responses compared to non-tumor-bearing counterparts. The genetic elimination of nociceptor neurons decreased brain Fos expression and mitigated the behavioral alterations induced by the presence of the tumor. While analgesic treatment restored nesting and cookie test behaviors, it did not fully restore voluntary wheel running indicating that pain is not the exclusive driver of such behavioral shifts. Unraveling the interaction between the tumor, infiltrating nerves, and the brain is pivotal to developing targeted interventions to alleviate the mental health burdens associated with cancer.