(a) Age-dependent changes in thymic architecture, as shown by representative H and E staining of thymic sections. Scale bars represent 150 µm. Medullary islands stain as light purple while cortical …
(a) Thymocyte cellularity fluctuations over age. (b) Ratio of %TEC:%thymocytes over age. (c) TEC cellularity over age for mTEC and cTEC FAC-sorted populations. Error bars represent mean +/- standard …
(a) FACS gating strategy to separate different T-cell subtypes. (b) Maturation trajectory and negative selection checkpoints for TCR-positive thymocytes. (c) Frequency of different thymocyte …
(a) A schematic of T cell receptor rearrangements used to design simulations. (b) Proportions of productive TCRs (y-axis) simulated at different sample sizes (x-axis). (c–f) Proportions of TCR alpha …
(a) A schematic showing the experimental design and FACS phenotypes of sorted cells for single-cell RNA-sequencing. Right panel shows cell composition fluctuations as a relative fraction of all …
(a) FACS gating strategy for TEC isolation. Identical gating strategies were used for all time points. (b) Filtering strategy to identify high-quality TEC libraries. (c–f) Fractions of libraries …
This figure depicts the assignment of single TEC from the Bornstein, the Park or our Ageing dataset to our nine TEC subtypes (left side) or the Bornstein/Park mTECI-IV nomenclature.
(a) Observed percentages (%) of TEC based on pre-scoring into FAC sorted subpopulations. (b) Estimated contributions of each FACS sort type to each single cell subtype through age. Each coloured dot …
Each subplot is a SPRING-layout of the shared nearest-neighbour graph of single TEC, derived from scRNA-seq transcriptional profiles as in Figure 2. Graph nodes are coloured according to the …
The X-axis shows the fraction of marker genes that overlap the specified pathway, the size of the dot represents the number of marker genes in the enriched pathway, and the colour of the dot …
The X-axis shows the fraction of marker genes that overlap the specified pathway, the size of the dot represents the number of marker genes in the enriched pathway, and the colour of the dot …
The heatmap shows the fraction of times that the libraries are co-clustered based on a variety of transformations, clustering methods and the number of features (Materials and methods). The heatmap …
(a) Enrichment of MSigDB biological pathways with age in mature cTEC, intertypical TEC and mature mTEC, annotated as in (Figure 2b). Bars denote normalised enrichment score (NES) for significant …
Each panel shows a heatmap and a MA-plot depicting the expression of the differentially expressed genes within the (a) mature cTEC, (b) intertypical TEC, or (c) mature mTEC subtypes. The heatmaps …
(a) TRA expression enriched in Mature mTEC and Post-Aire mTEC. Shown is the percentage of expressed genes that are classed as TRAs (see Materials and methods), for each TEC subtype across mouse …
(a) A schematic representing the transgenic Dox-inducible ZsGreen (ZsG) lineage tracing of β−5t-expressing mTEC precursors (top), and lineage tracing experiment in 1-week-old thymi (bottom). The …
(a) Representative FACS plots depicting the labelling efficiency of cTEC and mTEC of 1-week-old mice, 48 hr after treatment with 0.3 mg of doxycycline per mouse. (b) FACS gating strategy employed …
(a) The top genes that discriminate between TEC subtypes trained using single TEC from 1-week-old mice. Shown is the mean decrease in accuracy when dropping each gene in the random forest (left) and …
(a) The mean percentage of cTEC (left) and mTEC (right) labelled either 2 days (blue points) or 28 days (orange points) after doxycycline treatment of 3xtgβ5t mice. Error bars are the mean +/- the …
(a) cTEC/mTEC distribution as traditionally defined by flow cytometry using the surface markers Ly51 and reactivity to UEA1. (b) Proportion of ZsGreen-labelled total TEC detected at each time-point, …
(a) FAC-Sorting gating strategy for the isolation of ZsGreen +/- TEC subpopulations. (b) Multiplet detection using multiplexed hashtag oligos (HTO). Coloured bars denote the number of cells in each …
(a) Boxplots showing the distributions of the estimated deconvolution size factors (left) and number of detected genes for each single-cell cluster (right). (b) Uniform manifold approximation and …
(a) A mapping of single-cell clusters onto TEC subtype phenotypes. (b) A confusion matrix showing the proportions of single-cells in each TEC cluster (columns) and their classification into Ageing …
Boxplots showing the distribution of marker gene expression (y-axis) for TEC subtypes across TEC clusters (x-axis). Boxes are coloured by the inferred TEC subtype to which they belong.
Each panel is coloured according to the TEC subtype annotation and corresponds to Figure 6c. Arrows point to the positions of the relevant clusters in the UMAP.
(a) UMAP as in (Figure 6c) coloured by diffusion pseudotime (DPT) distance in the medullary lineage (see Figure 7—figure supplement 1 for details). (b) Changes in ZsGreen+ labelled TEC …
(a–b) Diffusion maps of single mTEC and Intertypical TEC, showing two medulla branches on DC1 vs DC2 (a) and DC2 vs DC3 (b). (c) A UMAP subset to mTEC and Intertypical TEC coloured by DPT on both …
(a–b) Diffusion maps of single cTEC and Intertypical TEC showing diffusion components (DC) (a) 1 vs. 2 and (b) 2 vs. 3. (c) A rug plot showing the positions of TEC sub-clusters (y-axis) along the …
Boxplot showing the proportion of TEC in each subtype cluster which express key genes linked to thymic involution and TEC identity: (a) Psmb11 (β5-t), (b) Ly6a (Sca1), (c) Bmp4, (d) Inhba (Activin …
Expression in single TEC of (a) Psmb11 , (b) Ccl21a , (c) both (Psmb11+/Ccl21a+; blue), and (d) expression levels of Psmb11 vs. Ccl21a, with TEC subtype indicated by colour. TEC with no detected …
Tables listing relevant experiment related information.
(Table 1) Numbers of TEC isolated in the single-cell experiment. (Table 2) Marker genes for each subtype of TEC. (Table 3) Tests conducted for marker proteins in TEC. (Table 4) Single-cell defined TEC subtypes and known concordant phenotypes. (Table 5) Details of antibodies used in flow cytometry staining panels to identify TEC and thymocytes undergoing negative selection. (Table 6) ADT primers for Droplet sequencing. (Table 7) HTO primers for Droplet sequencing. (Table 8) Tissue-specific gene classification via FANTOM5 Cage-Seq:Tissue samples were grouped into 27 broad groups based on the annotation data.