Abstract
Despite the essential role of thymic epithelial cells (TEC) in T cell development, the signals regulating TEC differentiation and homeostasis remain incompletely understood. In this study, we show a key in vivo role for the vitamin A metabolite, retinoic acid (RA), in TEC homeostasis. In the absence of RA signaling in TEC, cortical TEC (cTEC) and CD80loMHC class IIlo medullary TEC displayed subset-specific alterations in gene expression, which in cTEC included genes involved in epithelial proliferation, development, and differentiation. Mice whose TEC were unable to respond to RA showed increased cTEC proliferation, an accumulation of stem cell Ag-1hi cTEC, and, in early life, a decrease in medullary TEC numbers. These alterations resulted in reduced thymic cellularity in early life, a reduction in CD4 single-positive and CD8 single-positive numbers in both young and adult mice, and enhanced peripheral CD8+ T cell survival upon TCR stimulation. Collectively, our results identify RA as a regulator of TEC homeostasis that is essential for TEC function and normal thymopoiesis.
Introduction
The generation of functional T cells requires interactions between developing thymocytes and cortical and medullary thymic epithelial cells (TEC) (1). Cortical TEC (cTEC) promote the commitment, expansion, and differentiation of CD4−CD8− double-negative (DN) thymocytes to CD4+CD8+ double-positive (DP) thymocytes (2) and positive selection of DP thymocytes with TCR specificities able to recognize peptide–MHC complexes (3). In contrast, medullary TEC (mTEC), through their ability to express a broad range of tissue-restricted Ags (3), ensure the elimination of autoreactive TCR specificities within the single-positive (SP) thymocyte pool, are required for the generation of Foxp3+ regulatory T cells, and promote SP thymocyte maturation (4).
Mature cTEC and mTEC originate from a common bipotent progenitor present in both embryonic (5–7) and postnatal (8) thymus via lineage-committed cTEC and mTEC precursor intermediates (9, 10). The environmental signals regulating TEC development and differentiation are incompletely understood but include receptor activator for NF-κB ligand, CD40L, and lymphotoxin derived from developing thymocytes (11) and fibroblast growth factor-7 and -10 (12, 13), and insulin-like growth factor-1 and -2 (14, 15) from thymic mesenchymal cells. We recently found that gp38+ thymic mesenchymal cells produce the vitamin A metabolite, retinoic acid (RA), and that addition of a pan–RA receptor (RAR) antagonist to fetal thymic organ cultures enhances TEC proliferation, which was associated with increased expression of cTEC-lineage genes (16). Although these observations suggest that RA signaling may regulate TEC homeostasis, the in vivo role of RA signaling specifically in TEC in thymic development, homeostasis, and thymopoiesis remains unknown.
RA signals through RAR/retinoid X receptor heterodimers that function as ligand-induced transcription factors binding to RA-responsive elements in the regulatory regions of target genes to regulate gene transcription (17). To investigate the in vivo role of RA signaling in TEC, we generated mice whose TEC express a dominant-negative form of RARα (dnRARα) and are thus unable to respond to RA signals. Our results provide in vivo evidence that RA signaling in TEC regulates postnatal TEC development and expansion with downstream impacts on thymopoiesis.
Materials and Methods
Mice
Foxn1Cre.dnRARlsl/lsl, dnRARlsl/lsllo) and control cells (VCThi) at a 1:1 ratio. Mice were sacrificed 4 h postinjection and CTL activity in spleen was calculated as the ratio of target/control cells. All animal procedures were approved by the Lund/Malmö Animal Ethics Committee.
Cell preparation
For TEC analysis, thymus tissue was digested as previously described (16) and TEC were enriched by depleting CD45+ cells with CD45 microbeads (Miltenyi Biotec). For thymocyte and splenocyte analysis, organs were mashed through a 70-μm cell strainer and RBCs lysed using ACK solution.
Flow cytometry analysis and cell sorting
Cells were stained as previously described (16). All analyses were pregated on live singlet cells. Data acquisition was performed on an LSR II (BD Biosciences) and analyzed using FlowJo software (Tree Star). Sorting was performed on a FACSAria II (BD Biosciences). Postsort purity was >95%. For intracellular staining, cells were fixed and permeabilized using a Foxp3 staining buffer set (eBioscience) or 4% paraformaldehyde and 0.1% saponin. CCL21 staining was performed as previously described (18). TCR Vβ usage was assessed using the anti-mouse TCR Vβ screening panel (BD Biosciences). For analysis of BrdU incorporation, mice were injected i.p. with BrdU (1 mg; BD Biosciences), sacrificed after 16 h, and TEC suspensions were stained with FITC-conjugated BrdU (FITC BrdU flow kit; BD Biosciences).
Analysis of thymic sections
Frozen thymic sections were stained with H&E and analyzed as described (19
Quantitative PCR
Indicated samples were analyzed for the presence of the genomic stop cassette using primers neomycin_fw (5′-TGCCGAATATCATGGTGGAAAAT-3′) and neomycin_rev (5′-TCGTCAAGAAGGCGATAGAAGGCG-3′) and compared with the presence of the dnRAR construct. Expression of dnRAR was assessed using primers dnRAR_fw (5′-GACTGCAGAAGTGCTTTGAA-3′) and dnRAR_rev (5′-GCACTCGGGCTTGGGCA-3′) using Kapa SYBR Fast quantitative PCR master mix (Kapa Biosystems). In all cases, mRNA levels were normalized to β-actin.
T cell stimulation
Naive splenic T cells were incubated (5 × 104 cells per well) on anti-CD3 (145-2C11; Bio X Cell) precoated U-bottom 96-well plates with or without soluble anti-CD28 (37.51, 1 μg/ml; Bio X Cell) at 37°C and 5% CO2. Thymidine (1 μCi [3H]thymidine; PerkinElmer) was added after 48 h, and 18 h later cells were harvested with a cell harvester (Tomtec). [3H]thymidine was assessed on a 1450 MicroBeta TriLux (PerkinElmer). For assessing proliferation, cells were VCT labeled prior to incubation and harvested after 66 h. For Ca2+ flux, cells (4 × 106) were loaded with Fluo-4 (7 μM; Life Technologies) and Fura Red (35 μM; Molecular Probes) and incubated with anti-CD3 (10 mg/ml, 145-2C11; Bio X Cell) for 10 min at 37°C. Crosslinking Ab (1–20 μg, goat anti–Armenian hamster IgG; Jackson ImmunoResearch Laboratories) was added and Ca2+ flux monitored on a FACSAria II for 250 s.
RNA sequencing
RNA was prepared using an RNeasy Micro Kit (Qiagen) and amplified using an Ovation RNA-Seq System V2 (NuGEN) kit. cDNA libraries, generated with an Ovation Ultralow System V2 kit (NuGEN), were fragmented using a Bioruptor Pico (Diagenode) and sheared cDNA end repaired to generate blunt ends, ligated to Illumina adaptors with indexing tags, and purified using AMPure XP beads. Libraries were sequenced on a HiSeq 2500 System (Illumina) using 150-bp paired-end reads (Center for Biological Sequence Analysis, Technical University of Denmark, Copenhagen, Denmark). Sample reads were aligned to the mouse Ensembl genome assembly using HISAT-2 (version 2.0.1), and expression counts were quantified using htseq-count. For normalization, gene expression data of each sample were combined in R and genes with <25 counts in two samples within a single condition were excluded. Principal component analysis (PCA) was performed using the R function prcomp and PCA plots were generated using ggplot2. The R package Limma was used for differentially expressed (DE) analysis. Genes with a q value < 0.05 (Benjamini–Hochberg) were considered significant. The complete sets of DE genes (DEG) were tested for gene ontology (GO) biological process enrichment using gProfileR (version 0.6.1) (20). Heat maps were created using gplots function heatmap.2. RNA sequencing (RNA-seq) data were deposited in the National Center for Biotechnology Information Sequence Read Archive (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA454173) under accession number PRJNA454173.
Statistical analysis
All statistical analyses apart from RNA-seq experiments were performed using GraphPad Prism 7.0. Statistical analysis was performed using unpaired two-tailed Student t test. A p value < 0.05 was considered significant (*p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001).
Results
Adult mice whose TEC cannot respond to RA display enhanced numbers of cTEC
To investigate the role of RA signaling in TEC in vivo, dnRARlsl/lsl mice (21) that carry a loxP-flanked stop cassette preceding the coding sequence for a dnRARα (RARα403), previously used to assess the role of RA signaling in hematopoietic cells (22, 23), were crossed with Foxn1Cre mice (24) to generate Foxn1Cre.dnRARlsl/lsl mice. Deletion of the stop cassette in TEC from Foxn1Cre.dnRARlsl/lsl mice (hereafter termed Cre+) but not dnRARlsl/lsl littermates (hereafter termed Cre−) was confirmed at the genomic level by quantitative PCR (Fig. 1A), which resulted in expression of dnRAR by these cells (Fig. 1B).
Adult Foxn1Cre.dnRARlsl/lsl mice have increased numbers of cTEC. (A) Presence of the genomic STOP cassette and (B) relative expression of dnRAR mRNA in sorted TEC and CD45+ thymocytes from Cre− (open bar) and Cre+ (filled bar) mice. Fold levels represent the mean (SD) of (A) four and (B) five biological replicates. (C) Representative flow cytometry plots and (D) total number of cTEC (UEA-1−Ly51+), mTEClo (UEA-1+Ly51−CD80loMHC-IIlo), and mTEChi (UEA-1+Ly51−CD80hiMHC-IIhi) subsets (pregating: live single CD45−EpCAM+) from female (♀) and male (♂) Cre− and Cre+ mice. Data are from three experiments. **p < 0.01, ***p < 0.001. EpCAM, epithelial cell adhesion molecule; UEA-1, Ulex europeaus agglutinin-1.
Adult (8 wk) Cre+ and Cre− littermates showed similar thymic architecture, including cortical and medullary size (Supplemental Fig. 1A, 1B). However, quantitative flow cytometry analysis demonstrated a significant increase in cTEC numbers in both female and male Cre+ mice (Fig. 1C, 1D), whereas CD80lo MHC class II (MHC-II)lo mTEC (mTEClo) and CD80hiMHC-IIhi mTEC (mTEChi) as well as proportions of AIRE or CCL21 expressing mTEC were unchanged (Fig. 1D, Supplemental Fig. 1C). Thus, RA signaling in TEC appears to be required for cTEC homeostasis in adult mice.
RA regulates numerous distinct biological pathways in cTEC and mTEClo subsets
To further understand the role of RA signaling in TEC biology, cTEC, mTEClo, and mTEChi were sorted in triplicate from 8-wk-old Cre+ and Cre− littermates, and their gene expression profiles were assessed by RNA-seq. PCA confirmed reproducibility between the datasets and clearly separated cTEC, mTEClo, and mTEChi subsets (Fig. 2A). Only four genes were DE in mTEChi between Cre+ and Cre− mice (data not shown, q value < 0.05, Benjamini–Hochberg FDR), indicating that RA signaling has little impact on these cells. In contrast, 382 genes in cTEC and 408 genes in mTEClo were identified as DE between Cre+ and Cre− mice (Fig. 2B, 2C). Only 81 DEG were shared by cTEC and mTEClo, indicating that most direct or indirect RA-controlled genes are regulated in a TEC subset–specific manner (Fig. 2C).
Altered gene expression in TEC from Foxn1Cre.dnRARlsl/lsl mice. Comparative gene expression analysis of TEC subsets from 8-wk-old Cre− and Cre+ mice. cTEC, mTEClo, and mTEChi subsets were sorted in three experiments from five to six pooled female mice per experiment. (A) PCA of all TEC subsets and (B) volcano plots showing DEG in Cre+ versus Cre− mice from RNA-seq analysis of cTEC and mTEClo subsets. DE (FDR < 0.05, colored dots) and non-DE (gray dots) are shown. (C) Venn diagram of DEG in indicated TEC subsets from Cre+ mice. ↑, upregulated in the absence of RA signaling; ↓, downregulated in the absence of RA signaling; ↕, downregulated in one TEC subset but upregulated in the other. (D and E) Lists of most significant GO terms in cTEC (D) and mTEClo (E) from Cre+ mice.
GO enrichment analysis of the DEG for cTEC and mTEClo demonstrated significantly enriched GO terms for epithelial cell proliferation and epithelial cell differentiation in cTEC (Fig. 2D) whereas those for mTEClo included processes of secretion, cell adhesion, and ion transport (Fig. 2E).
Impaired RA signaling in TEC results in increased cTEC proliferation
Given the enhanced cTEC numbers in adult Cre+ mice (Fig. 1D) and the RA-regulated GO term related to epithelial cell proliferation in cTEC (Figs. 2D, 3A, 3B), we assessed TEC proliferation in adult Cre+ and Cre− littermates. Sixteen hours after BrdU injection cTEC, but not mTEC, from Cre+ mice showed increased BrdU incorporation compared with cTEC from Cre− mice (Fig. 3C, Supplemental Fig. 1D). Thus, enhanced cTEC proliferation may contribute, at least in part, to the increased numbers of cTEC in adult Cre+ mice.
Increased proliferation of cTEC in Foxn1Cre.dnRARlsl/lsl mice. (A) Cell proliferation GO term network and (B) heat maps of DEG in the epithelial cell proliferation GO term. ↑, upregulated in the absence of RA signaling; ↓, downregulated in the absence of RA signaling. Circle size is relative to number of DEG in GO term. (C) Percentage of BrdU+ cells in indicated TEC subsets 16 h postinjection. Data are from three experiments. *p < 0.05.
Impaired RA signaling in TEC results in increased numbers of stem cell Ag-1hi cTEC and reduced numbers of mTEC during postnatal TEC expansion
The identification of RA-regulated GO terms related to epithelial development and differentiation in Cre+ cTEC (Figs. 2D, 4A, 4B) prompted us to investigate alterations in TEC in early life, during the period of TEC expansion (25). The Ly51+ cTEC compartment of adult mice contains bipotent TEC progenitors (9, 26), and stem cell Ag-1 (Sca-1) has been implicated as a marker that is enriched in cells with bipotent progenitor potential (9, 26). We thus assessed surface expression of Sca-1 within the cTEC compartment over time. Cre+ mice had a significant increase in Sca-1hi cTEC numbers from birth into adulthood compared with Cre− littermates, whereas Sca-1lo cTEC numbers were comparable (Fig. 4C, 4D). Furthermore, Cre+ mice had reduced numbers of mTEC during the first 3 wk of life (Fig. 4E), a time period where mTEC derive primarily from bipotent progenitors displaying a cTEC-specific phenotype (27–30). These results demonstrate that RA signaling modulates cTEC/mTEC ratios in early life, potentially by regulating the ability of bipotent cTEC-like progenitors to give rise to mTEC.
Distinct developmental kinetics of TEC and increase of Sca-1hi cTEC in Foxn1Cre.dnRARlsl/lsl mice. (A) Developmental process GO term network and (B) heat maps of DEG in the epithelium development GO term. ↑, upregulated in the absence of RA signaling; ↓, downregulated in the absence of RA signaling. Circle size is relative to number of DEG in GO term. (C) Representative flow cytometry plots of Sca-1hi and Sca-1lo cTEC and (D and E) total number of indicated TEC subsets over time. Mean (SD) of three experiments with three to six mice per time point and experiment. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
Increased TEC regeneration after thymic injury in the absence of RA signaling
To investigate the role of RA signaling in TEC during thymic regeneration, 8-wk-old Cre+ mice and Cre− littermates were injected i.p with a single dose of dexamethasone (31). Such treatment results in rapid thymus involution, followed by an extended recovery phase (32). In line with the above studies, dexamethasone treatment resulted in a dramatic reduction in thymocyte and TEC numbers 2 d postinjection (Fig. 5). Thymocyte as well as mTEClo and mTEChi numbers were significantly lower in Cre+ compared with Cre− mice (Fig. 5A, 5C), whereas cTEC loss remained the same between both groups (Fig. 5B). In contrast, during the recovery phase 12 d postinjection, cTEC and mTEClo numbers were significantly higher in Cre+ mice compared with Cre− littermates, despite similar levels of thymocyte restoration (Fig. 5). Thus, RA signaling in TEC appears to influence TEC recovery during injury-induced thymic regeneration.
Enhanced TEC regeneration in Foxn1Cre.dnRARlsl/lsl mice upon glucocorticoid-induced thymic injury. (A–C) Impact of glucocorticoid-induced thymic injury on (A) thymic cellularity, (B) Sca-1hi or Sca-1lo cTEC, and (C) mTEClo and mTEChi numbers. Data are from two experiments. **p < 0.01, ***p < 0.001, ****p < 0.0001.
RA signaling in TEC impacts conventional T cell development
Given the importance of TEC in T cell development, we next assessed the impact of RA signaling in TEC on thymopoiesis. Cre+ mice had reduced thymic cellularity from 1 d to 5 wk after birth compared with Cre− littermates (Fig. 6A), whereas adult Cre+ mice displayed normal thymocyte cellularity (Fig. 6A). Young (12-d-old) Cre+ mice had reduced numbers of DN, DP, and SP thymocytes (Fig. 6B). Among DN subsets, DN3 and DN4, but not ETP or DN2 numbers, were reduced (Supplemental Fig. 2A), whereas localization of CD25+ thymocytes, which in the outer cortex and subcapsular zone almost exclusively represent DN2 and DN3 thymocytes (33), to the subcapsular zone appeared normal (Fig. 6C). Total DP numbers were similar in 5- and 8-wk-old Cre+ and Cre− mice (Fig. 6B); however, numbers of CD69intTCRβint DP (DPint) and CD69hiTCRβhi DP (DPhi), which represent minor DP subpopulations undergoing and immediately subsequent to positive selection (34, 35), were reduced in Cre+ mice at 5 (Fig. 6D, 6E) but not 8 wk (data not shown). CD4SP and CD8SP numbers were reduced in Cre+ mice at all analyzed time points (Fig. 6B), including both immature CD24hiCD62Llo and mature CD24loCD62Lhi SP subsets (Supplemental Fig. 2C, 2D). Expression of TEC genes previously implicated in regulating thymocyte development and survival did not differ between Cre+ and Cre− mice (Supplemental Fig. 2E, 2F).
Foxn1Cre.dnRARlsl/lsl mice display reduced thymic cellularity in early life and maintain decreased SP numbers into adulthood. (A) Thymic cellularity and (B) total number of DN (CD4−CD8−), DP (CD4+CD8+), CD4SP (CD4+CD8−), and CD8SP (CD4−CD8+) thymocytes at indicated time points. Data are from (A) three to five or (B) three experiments. (C) Confocal images of thymus sections of 12-d-old mice stained for CD25 (green) and CD8 (purple). Images are representative of four mice per strain. Scale bars, 100 μm. (D) Representative flow cytometry plots and (E) total number of DP subsets before (DPlo), during (DPint) and after selection (DPhi) in 5-wk-old mice. Data are from three experiments. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
Finally, the proportions of Foxp3+ cells among CD4SP thymocytes and numbers of γδ T cells, NKT cells, and intraepithelial lymphocyte precursors (36) were similar between Cre+ and Cre− mice (Supplemental Fig. 2G–L). Taken together, these results suggest that in the absence of RA signaling in TEC, the generation of conventional CD4SP and CD8SP thymocytes is moderately impaired.
Normal positive selection in mice that lack RA signaling in TEC
The reduction in postselection DP and SP subsets in Cre+ mice prompted us to determine whether the absence of RA signaling in TEC impacted positive selection. Cre+ and Cre− mice were backcrossed onto a TCR transgenic OT-I and OT-II background to assess selection in the setting of a fixed monoclonal TCR repertoire, allowing for detection of subtle changes in thymocyte selection (37, 38). Similar frequencies of OT-I CD8SP and OT-II CD4SP were found in Cre+ and Cre− mice (Fig. 7A, Supplemental Fig. 3A). Additionally, CD5 expression on DP thymocytes, as a measure of TCR signaling strength (39), was similar in Cre+ and Cre− mice (Fig. 7B). Finally, the overall TCR repertoire of CD4SP and CD8SP thymocytes, as assessed by flow cytometry analysis of Vβ gene usage, was similar in Cre+ and Cre− mice, although we did observe significant alterations in the gene usage of Vβ2 and Vβ5.1/5.2 in CD8SP thymocytes (Supplemental Fig. 3B). Collectively, these findings suggest that RA signaling in TEC does not impact positive selection in a major way.
Normal positive selection and thymic egress in Foxn1Cre.dnRARlsl/lsl mice. (A) Percentage of CD8SP or CD4SP cells of live thymocytes in Cre− (○) and Cre+ (●) OT-I and OT-II mice. (B) CD5 expression of DP thymocytes. Data are representative of two experiments with four mice per group and experiment. (C) Total number of splenic naive CD4+ and CD8+ cells. (D) Frequencies of GFP+ RTE within splenic naive (CD62LhiCD44lo) CD4+ and CD8+ T cells in mice reconstituted with Rag2pGFP bone marrow and (E and F) mean fluorescence intensity (MFI) of (E) CD24 and (F) Qa2 expression on naive splenic CD4+ and CD8+ T cells. Data are from three (A and C) or two (D–F) experiments.
RA signaling in TEC is not required for thymic egress of mature thymocytes
Given the reduction in SP thymocytes, we next assessed the peripheral T cell compartment of Cre+ mice. The total number of naive (CD62LhiCD44lo) CD4+ and CD8+ splenic T cells did not differ between 8-wk-old Cre+ and Cre− littermates (Fig. 7C). Given that recent thymic emigrants (RTEs) represent only a minor subset of the naive T cell pool in adult mice (40), we next assessed thymic egress in irradiated Cre+ and Cre− littermates reconstituted with bone marrow cells from Rag2p-GFP mice (41). GFP signal in these mice peaks in DP thymocytes and remains detectable 1−2 wk after cessation of Rag expression, thus allowing for the identification of GFP+ RTEs in the periphery (41). The frequencies of GFP+ RTEs among naive splenic CD4+ and CD8+ T cells and GFP expression levels on RTEs did not differ between Cre+ and Cre− BM chimeras (Fig. 7D, Supplemental Fig. 3C, 3D). Additionally, surface expression of CD24 and Qa2, markers used to discriminate RTEs from the mature naive T cell pool (41), was similar in unmanipulated Cre+ and Cre− mice (Fig. 7E, 7F). Taken together, these results suggest that the absence of RA signaling in TEC does not impact SP thymocyte residence time in the medulla or thymic egress.
Impaired RA signaling in TEC alters the response of naive CD8+ T cells to TCR stimulation
Finally, to assess whether the absence of RA signaling in TEC impacted peripheral naive T cell functionality, sorted splenic naive CD4+ and CD8+ T cells from Cre+ and Cre− mice were stimulated with anti-CD3 or anti-CD3/anti-CD28 and their proliferation was assessed by thymidine incorporation (Fig. 8A, 8B). CD8+ T cells from Cre+ mice incorporated thymidine to a greater extent than did their Cre− counterparts over a wide range of anti-CD3 concentrations (Fig. 8B). Consistent with this, more CD8+ T cells were recovered from anti-CD3 and anti-CD3/anti-CD28 cultures when cells originated from Cre+ mice (Fig. 8C). In contrast, induction of CD69 (Fig. 8D, Supplemental Fig. 4A), calcium flux in response to anti-CD3 stimulation across a range of crosslinking Ab concentrations (Fig. 8E, data not shown), and cellular division assessed by VCT dilution (Supplemental Fig. 4B, 4C) were similar in CD8+ T cells from Cre+ and Cre− mice. Finally, OVA and anti-CD40 immunized Cre+ and Cre− mice displayed a similar ability to kill OVA-pulsed target cells in vivo (Fig. 8F, 8G), indicating that CD8+ T cells from Cre+ mice are capable of differentiating into functional CTL. Collectively, these results suggest that although naive CD8+ T cells show enhanced survival following anti-CD3 stimulation in vitro, their ability to respond to anti-CD3 and to differentiate into functional CTL appears normal.
Increased numbers of naive CD8+ T cells from Foxn1Cre.dnRARlsl/lsl mice following anti-CD3 stimulation in vitro. (A and B) [3H]thymidine incorporation of anti-CD3– or anti-CD3/anti-CD28–stimulated naive (A) CD4+ or (B) CD8+ T cells. Mean (SD) of three experiments with four mice per group and anti-CD3 condition in each experiment. Results were normalized between experiments by setting the mean cpm for Cre− at 5 μg/ml anti-CD3 to 100% for each experiment (cpm counts for the three experiments: CD4+, 89,942, 58,209, 233,127; CD8+, 89,513, 38,643, 43,635). (C) Total number of CD8+ T cells and (D) frequencies of CD8+CD69+CD62Llo T cells after indicated stimulation. Mean (SD) of four mice per group and anti-CD3 condition in each experiment. (E) Ca2+ response of naive splenic CD8+ T cells to anti-CD3 stimulation. Data are representative of two experiments with four mice per group. (F) Representative flow cytometry plots mice and (G) percentage of killed target cells 4 h postinjection into OVA and anti-CD40–immunized mice. Data are from three experiments. *p < 0.05, **p < 0.01, ***p < 0.001, ****p < 0.0001.
Discussion
In this study, we assessed the impact of RA signaling in TEC for TEC and thymocyte development in vivo. Our results confirm and extend our previous in vitro observations in fetal thymic organ cultures (16) that inhibition of RA signaling enhances cTEC expansion. Additionally, we show that the absence of RA signaling alters the expression of numerous, yet largely nonoverlapping, genes in cTEC and mTEClo, with DEG in cTEC associated with epithelial cell proliferation and development. In the absence of RA signaling in TEC, immature Sca-1hi cTEC numbers increased whereas mTEC numbers decreased during early postnatal thymic expansion, resulting in reduced thymopoiesis. Inhibition of RA signaling in TEC also led to reduced SP numbers, and naive CD8+ T cells generated under these conditions displayed increased survival following anti-CD3 stimulation in vitro. Collectively, our results suggest that RA signaling in TEC in vivo is required for normal TEC development and function.
Lack of RA signaling in TEC had little impact on the gene expression profile of mTEChi, indicating that mature mTEChi in vivo are not exposed to or are unresponsive to RA. In contrast, significant changes were observed in the gene expression profiles of cTEC and mTEClo. DEG between cTEC and mTEClo were largely nonoverlapping, consistent with previous studies in embryonic fibroblasts and stem cells (42), suggesting that gene regulation by RA is highly cell type specific. Indeed, the biological processes controlled by RA in cTEC but not mTEClo included cell development and differentiation, GO terms previously identified to be overrepresented in cTEC when compared with mTEC (43). Which proportion of these DEG are direct RA targets or regulated indirectly by RA through the regulation of, for example, transcription factors, nuclear receptors, or histone-modifying enzymes (44) is difficult to assess using current technologies given limitations in TEC numbers.
GO analysis of DEG in cTEC in the absence of RA signaling predicted alterations in epithelial cell proliferation and included upregulation of FOXN1 and insulin-like growth factor-1, known drivers of TEC proliferation (11), indicating a potential role for RA signaling in limiting cTEC cycling. Consistent with this, as well as our prior in vitro studies (16), we observed that cTEC but not mTEC proliferation was increased in the absence of RA signaling. This increase in proliferation likely contributes, at least in part, to the increased numbers of cTEC observed in Cre+ mice, although we cannot exclude additional effects.
Consistent with the GO analysis predicting alterations in epithelial differentiation, absence of RA signaling in TEC resulted in increased numbers of Sca-1hi cTEC, which in early postnatal life was associated with decreased numbers of mTEC. mTEC derive from bipotent progenitors that bear a cTEC phenotype during embryonic and early postnatal thymic expansion (27, 28, 45), and Sca-1hi cTEC have previously been described to contain bipotent progenitors in the postnatal thymus (9, 26). Thus, although the identity of bipotent progenitors within the postnatal cTEC compartment remains unknown, our results indicate that RA signaling in such progenitors may serve to promote mTEC development. Although the underlying direct or indirect RA-regulated genes, or sets of genes, responsible for these alterations remain unclear, these results are in line with the role of RA in regulating the differentiation of intestinal epithelial transit amplifying progenitors into secretory and enteroendocrine cells in vivo (46) and various epithelial cell lines in vitro (47, 48). Future experiments assessing the gene expression profile and turnover of Sca-1hi and Sca-1lo cTEC subsets in Cre+ and Cre− mice may help identify potential additional roles of RA in cTEC development from bipotent precursors.
In contrast to early postnatal life, mTEC in adult mice are maintained through mTEC-restricted progenitors (10, 30, 49). Cre+ mice showed a normalization of mTEC numbers with age when mice transition to an adult mTEC developmental program, indicating that RA may promote the differentiation of mTEC from bipotent, but not mTEC-restricted, precursors. The full recovery of mTEC numbers in adult Cre+ mice further suggests that RA signaling may confer a proliferative or survival advantage for mTEC-committed cells. In support of this, we found enhanced restoration of mTEC numbers in the absence of RA signaling following thymic injury, a process driven by mTEC-committed precursors (49).
In the absence of RA signaling in TEC we observed age-dependent alterations in thymopoiesis, with a significant reduction in DN, DP, and SP thymocytes at 12 d of age, which became restricted to SP thymocytes with age. At 12 d, ETP and DN2 numbers appeared normal in Cre+ mice whereas DN3 numbers were reduced, indicating that RA signaling in TEC is not required for thymic seeding, but potentially provides an optimal environment for DN3 survival and development. RA deficiency impairs expression of the Notch ligand Delta-like 1 in the embryonic spinal cord (50). Although expression of the Notch ligand Delta-like 4, a key factor for survival and differentiation of DN3 thymocytes (51), was similar in 8-wk-old Cre+ and Cre− mice, further analysis of the transcriptional profile of cTEC from 12-d-old Cre+ and Cre- mice might shed light on why DN3 thymocyte numbers are reduced at this age. Finally, although reduced DN3 numbers in 12-d-old Cre+ mice may partially underlie the reduction in subsequent thymocyte subsets, mature thymocyte numbers were more significantly reduced at this age, indicating potential additional roles of RA signaling in TEC in early thymopoiesis.
Total DP thymocyte numbers recovered with age in Cre+ mice, potentially as a result of niche filling; however, postselection DP subsets were reduced at 5 wk of age whereas SP thymocyte numbers remained reduced in adults. Although we cannot rule out subtle changes in the selecting peptide repertoire, our results indicate that positive selection is largely unaffected in Cre+ mice. However, alterations in the gene usage of several Vβ gene families might suggest minor effects on the resulting TCR repertoire, the significance of which remains to be determined. The reduced numbers of SP thymocytes in Cre+ mice appear not to be due to alterations in thymic egress or SP maturation. Moreover, expression of TEC genes previously implicated in thymopoiesis did not differ between Cre+ and Cre− mice. We thus speculate that mTEC of Cre+ mice, despite normal numbers, are qualitatively less capable of supporting SP thymocyte survival or, alternatively, induce stronger negative selection. Such effects are likely mediated by as yet unidentified RA-regulated genes or sets of genes impacting TEC–thymocyte interactions or through alterations in niche availability.
A final finding from our studies was that absence of RA signaling in TEC impacted naive CD8+ T cells. Despite normal numbers in vivo, naive splenic CD8+ T cells from Cre+ mice incorporated more thymidine following anti-CD3 stimulation in vitro. Anti-CD3–induced calcium flux, CD69 expression and proliferation, and the in vivo CTL activity of CD8+ T cells from Cre+ and Cre− mice were similar, indicating normal responsiveness to TCR stimulation. In contrast, the number of CD8+ T cells recovered from anti-CD3 cultures was higher when cells originated from Cre+ mice. A possible explanation for these findings is that T cell development within an RA signaling–deficient TEC environment alters the threshold by which CD8+ T cells undergo activation-induced cell death, perhaps through inhibition of Nur77 or its family members, known to induce TCR-mediated apoptosis during negative selection but also in mature T cells (52). Analysis of signaling pathways downstream of the TCR, as well as assessment of primary and recall CD8+ T cell responses in Cre+ and Cre− mice, should help elucidate the extent to which RA signaling deficiency in TEC impacts peripheral T cell function.
In summary, in this study, we identify RA signaling in TEC as a key regulator of TEC homeostasis and function. Future studies focusing on the identification of direct and indirect RA target genes in TEC subsets should help provide a better understanding of the molecular mechanisms underlying these findings.
Disclosures
The authors have no financial conflicts of interest.
Acknowledgments
We thank Dr. Julien Vandamme (Technical University of Denmark) and Dr. Kirstine G. Belling (University of Copenhagen, Copenhagen, Denmark) for technical support and Ann-Charlotte Selberg (Lund University) for animal care and genotyping.
Footnotes
This work was supported by Swedish Natural Science Research Council Grant 621-2013-4048, Danish Medical Research Council Grant DFF-1331-00024, and by funds from the IngaBritt and Arne Lundbergs Foundation. G.A. is supported by Medical Research Council Programme Grant MR/N000919/1 and Cancer Research UK Cancer Immunology Grant CRUK A20753. A.J.W. is supported by Wellcome Trust Seed Award 204375/Z/16/Z.
The sequences presented in this article have been submitted to the National Center for Biotechnology Information Sequence Read Archive (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA454173) under accession number PRJNA454173.
The online version of this article contains supplemental material.
Abbreviations used in this article:
- cTEC
- cortical TEC
- DE
- differentially expressed
- DEG
- DE gene
- DN
- double-negative
- dnRARα
- dominant-negative form of RARα
- DP
- double-positive
- GO
- gene ontology
- MHC-II
- MHC class II
- mTEC
- medullary TEC
- PCA
- principal component analysis
- RA
- retinoic acid
- RAR
- RA receptor
- RNA-seq
- RNA sequencing
- RTE
- recent thymic emigrant
- Sca-1
- stem cell Ag-1
- SP
- single-positive
- TEC
- thymic epithelial cell
- VCT
- violet cell tracer.
- Received March 20, 2018.
- Accepted May 7, 2018.
- Copyright © 2018 by The American Association of Immunologists, Inc.



















