Integrative analysis of scRNA-seq and scATAC-seq revealed transit-amplifying thymic epithelial cells expressing autoimmune regulator

  1. Takahisa Miyao
  2. Maki Miyauchi
  3. S Thomas Kelly
  4. Tommy W Terooatea
  5. Tatsuya Ishikawa
  6. Eugene Oh
  7. Sotaro Hirai
  8. Kenta Horie
  9. Yuki Takakura
  10. Houko Ohki
  11. Mio Hayama
  12. Yuya Maruyama
  13. Takao Seki
  14. Hiroto Ishii
  15. Haruka Yabukami
  16. Masaki Yoshida
  17. Azusa Inoue
  18. Asako Sakaue-Sawano
  19. Atsushi Miyawaki
  20. Masafumi Muratani
  21. Aki Minoda
  22. Nobuko Akiyama  Is a corresponding author
  23. Taishin Akiyama  Is a corresponding author
  1. Laboratory for Immune Homeostasis, RIKEN Center for Integrative Medical Sciences, Japan
  2. Immunobiology, Graduate School of Medical Life Science, Yokohama City University, Japan
  3. Laboratory for Cellular Epigenomics, RIKEN Center for Integrative Medical Sciences, Japan
  4. YCI Laboratory for Immunological Transcriptomics, RIKEN Center for Integrative Medical Sciences, Japan
  5. YCI Laboratory for Metabolic Epigenetics, RIKEN Center for Integrative Medical Sciences, Japan
  6. Laboratory for Cell Function Dynamics, RIKEN Center for Brain Science, Japan
  7. Transborder Medical Research Center, and Department of Genome Biology, Faculty of Medicine, University of Tsukuba, Japan
8 figures, 1 table and 9 additional files

Figures

Figure 1 with 3 supplements
Droplet-based single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) analysis of thymic epithelial cells (TECs) in 4-week mice.

(A) Uniform manifold approximation and production (UMAP) plot of scATAC-seq data from TECs (EpCAM+ CD45 TER119) from 4-week mice. Cell clusters are separated by colors and numbers in the plot. The two datasets were integrated using the Seurat package. The graph on the right shows percentages of each cluster in the total number of cells detected (15,255 cells). (B) Chromatin accessibility of typical marker genes of TECs. Accessibility in each gene region is represented in red. (C) Violin plot depicting chromatin accessibility in Aire and Cd80 gene regions in each cluster. (D) Pseudo-bulk accessibility tracks of the Aire gene region in each cluster (upper panels) and frequency of sequenced fragments within the Aire gene region of individual cells in cluster 0, 2, and 4 (lower panels).

Figure 1—figure supplement 1
Violin plot of chromatin accessibility in thymic epithelial cell (TEC) marker gene regions in each cluster.
Figure 1—figure supplement 2
Comparison of uniform manifold approximation and production (UMAP) plots of single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) between two experiments.

(A) Uniform manifold approximation and production (UMAP) plots of single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) data of murine thymic epithelial cells (TECs). Two scATAC-seq data (Experiment #1 and Experiment #2) were separately shown after the integration. The graph on the right shows percentages of each cluster. (B) Chromatin accessibility of typical marker genes of TECs. Two scATAC-seq data (Experiment #1 and Experiment #2) were separately shown after the integration. Accessibility in each gene region is represented in red.

Figure 1—figure supplement 3
Comparison of single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) data between two experiments.

(A) Violin plot depicting chromatin accessibility in Aire gene regions in each cluster. Two single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) data (Experiment #1 and Experiment #2) were separately shown after the integration. (B) Pseudo-bulk accessibility tracks of the Aire gene region in each cluster (upper panels) and frequency of sequenced fragments within the Aire gene region of individual cells in cluster 0, 2, and 4 (lower panels). Two scATAC-seq data (Experiment #1 and Experiment #2) were separately shown after the integration.

Figure 2 with 6 supplements
Droplet-based single-cell RNA sequencing (scRNA-seq) analysis of thymic epithelial cells (TECs) in 4-week mice.

(A) Uniform manifold approximation and production (UMAP) plot of scRNA-seq data from TECs (EpCAM+ CD45 TER119) from 4-week mice. Cell clusters (R0 to R17) are indicated by colors and numbers in the plot. The graph on the right shows the percentages of each cluster in the total number of cells detected (11,792 cells). (B) Violin plots depicting expression levels of typical TEC marker genes in each cluster.

Figure 2—figure supplement 1
Heatmap of the top five genes selectively expressed in each subcluster.

Yellow color indicates high expression.

Figure 2—figure supplement 2
Comparison of single-cell RNA sequencing data between two experiments.

(A) Uniform manifold approximation and production (UMAP) plots of thymic epithelial cell (TEC) single-cell RNA sequencing (scRNA-seq) data. Two scRNA-seq data (Experiment #1 and Experiment #2) were separately shown after the integration using Seurat package. The lower left graph shows percentages of each cluster. The lower right graph shows cumulative plot for percentages of clusters. (B) UMAP plots of TEC scRNA-seq data. Two scRNA-seq data (Experiment #1 and Experiment #2) were separately shown after merging two data without using integration functions. The lower left graph shows percentages of each cluster. The lower right graph shows cumulative plot for percentages of clusters.

Figure 2—figure supplement 3
Violin plots of marker genes and cell cycle-related genes.

(A) Violin plots for expression level of typical thymic epithelial cell (TEC) marker genes in single-cell RNA sequencing (scRNA-seq) analysis of TECs. (B) Violin plots for expression level of cell cycle-related genes in scRNA-seq analysis of TECs.

Figure 2—figure supplement 4
Comparison of single-cell RNA sequencing data between this study and previous studies.

(A) Integration of single-cell RNA sequencing (scRNA-seq) data derived from a previously reported well-based study (Bornstein et al., 2018) and scRNA-seq data derived from the present droplet-based study. scRNA-seq data from the two studies were integrated. (B) Uniform manifold approximation and production (UMAP) projections of the two scRNA-seq datasets are shown. Data from a previous study (Dhalla et al., 2020) were re-analyzed and integrated with data from the present study. Expression of typical marker genes in each data. (C) UMAP projections of our scRNA-seq data (upper left) and re-analyzed scRNA-seq data from a previous study (Wells et al., 2020) (upper right) are shown. In lower graphs, cell types in our single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) data were assigned by transferring the cell type annotation of Well’s scRNA-seq data.

Figure 2—figure supplement 5
Dot plots of expression of marker genes for our single-cell RNA sequencing (scRNA-seq) clusters in clusters of other scRNA-seq datasets (Bornstein et al., 2018; Dhalla et al., 2020; Wells et al., 2020).
Figure 2—figure supplement 6
Dot plots of expression of marker genes for Baran-Gales’s single-cell RNA sequencing (scRNA-seq) clusters (Baran-Gale et al., 2020) in clusters of our scRNA-seq datasets.
Figure 3 with 4 supplements
Integrative analysis of single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) data and single-cell RNA sequencing (scRNA-seq) data of thymic epithelial cells (TECs).

(A) Gene expression was predicted from scATAC-seq data using Signac. Individual cells in the cluster from scATAC data (clusters 0 to 11) were assigned and transferred to the uniform manifold approximation and production (UMAP) plot of scRNA-seq cluster (R0 to R17). (B) Correlation between clusters derived from scATAC-seq and scRNA-seq datasets of TECs. Cell types were annotated in scATAC dataset of TECs by transferring clusters from an scRNA-seq dataset. (C) Ratios of cells assigned to each scRNA-seq cluster in each scATAC cluster.

Figure 3—figure supplement 1
Integrative analysis of single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) data and single-cell RNA sequencing (scRNA-seq) data of thymic epithelial cells (TECs).

Gene expression was predicted from individual cells in scATAC-seq data. Individual cells in the scATAC data were assigned to an scRNA-seq cluster (R0 to R17).

Figure 3—figure supplement 2
Comparison of Integrative analysis of single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) data and single-cell RNA sequencing (scRNA-seq) data of thymic epithelial cells (TECs) between two experiments.

After cell types were annotated in single-cell assays for transposase-accessible chromatin (scATAC) dataset of thymic epithelial cells (TECs) by transferring clusters from an single-cell RNA sequencing (scRNA-seq) dataset, two scATAC-seq data were split.

Ratio of cells assigned to each scRNA-seq cluster in each scATAC cluster was shown in the lower graph.

Figure 3—figure supplement 3
Pseudo-bulk accessibility tracks and frequency of sequenced fragments.

Typical differentially accessible regions between clusters 0 and 4 are depicted.

Figure 3—figure supplement 4
Trajectory analysis of single-cell RNA sequencing data.

(A) RNA velocity analysis of single-cell RNA sequencing (scRNA-seq) data. (B) Monocle 3 trajectory analysis of single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) data. The trajectory was manually started from cluster 3 or 4.

Figure 4 with 1 supplement
Subcluster analysis of the thymic epithelial cell (TEC) subset expressing a high level of cell cycle-related genes.

(A) Uniform manifold approximation and production (UMAP) plot of single-cell RNA sequencing (scRNA-seq) data of each subcluster (R1A to R1G) in R1. Cell subclusters (R1A to R1G) are separated by colors and numbers in the plot. The graph on the right shows the percentages of each cluster in the parent R1 cluster. (B) Heatmap of the top five genes selectively expressed in each subcluster. Yellow color indicates high expression. (C) Expression levels of Aire and Ccl21a in the subcluster are exhibited as violin plots. (D) Expression levels of Aire, Ccl21a, and Cd80 in the subcluster are shown in dot plots. (E) Expression levels of marker genes for G2/M phase (upper, Mki67 and Hmgb2) and S phase (lower, Tyms and Slbp) in the subcluster are exhibited as violin plots. (F) Percentage of cells predicted as each cell cycle (G1, G2M, and S phases) in the subclusters.

Figure 4—figure supplement 1
Integrative analysis of cluster 4 from single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) data and subclusters in R1 from single-cell RNA sequencing (scRNA-seq) data.

Gene expression was predicted from individual cells in cluster 4 of scATAC-seq data. Individual cells in the cluster 4 were assigned to the subclusters (R1A to R1G) in R1.

Figure 5 with 3 supplements
A highly proliferative subset of Aire+ CD80hi medullary thymic epithelial cells (mTECs).

(A) Schematic depiction of cell cycles and Fucci fluorescence. (B) Flow cytometric analysis of TECs from Fucci2a mice crossed with Aire-GFP-reporter mice. The gating strategy is shown. Intensities of GFP to monitor Aire expression in each subset (Venus+ CD80hi mTEC, Venus CD80hi mTEC, and CD80lo mTEClo) are shown in the right panels. Left, Airegfp/+:: Fucci2a; right, control::Fucci2a. Typical figures of three independent experiments are exhibited. (C) Immunostaining of a sorted Venus+ CD80hi mTEC subset via anti-Aire antibody and DAPI (nucleus staining). Typical panels of three independent experiments are exhibited. Scale bars, 10 μm. (D) Immunostaining of thymic sections from Fucci2a mice with anti-Aire and anti-keratin-5 (Krt5) antibodies. Typical panels of three independent experiments are exhibited. Scale bars, 10 μm. (E) Scatter plots of RNA sequencing data from Venus+ CD80hi mTEC and Venus CD80hi mTEC subsets. The left panel shows a plot of all detected genes and the right panel shows tissue-specific antigen (TSA) genes detected. N = 3. (F) Atypical RNA sequencing tracks of Aire, typical Aire-dependent TSA genes (Ins and Sst), Fezf2, and Top2a (a marker of G2/M phase). (G) Scatter plots and volcano plots of RNA sequencing data from Venus+ CD80hi mTEC and Venus CD80hi mTEC subsets. Upper panels show Aire-dependent TSAs, lower panels show Aire-independent TSAs. Red dots in volcano plots indicate genes for which expression differed significantly (twofold change and FDR p < 0.05) in Venus+ and Venus CD80hi mTEC subsets. Numbers of differentially expressed genes are shown in the panels. N = 3. Y axis is log10 of FDR p-value.

Figure 5—figure supplement 1
Flow cytometric analysis of medullary thymic epithelial cells (mTECs) from Fucci mouse.

Ly51UEA-1+ TECs were separated by expression of Venus+ and CD80 (left). CD80 expression was detected in mCherryhi mTECs (red) and Venus+ mTECs (green) in right panels.

Figure 5—figure supplement 2
Immunostaining of thymic sections from Fucci2a mice with anti-GFP (for Venus staining, green) and anti-keratin-5 (Krt5, blue) antibodies.

Typical panels of three independent experiments are exhibited. Scale bars, 100 μm.

Figure 5—figure supplement 3
Gene set enrichment analysis of differentially expressed genes between Venus and Venus + cells.

Gene sets of Aire-dependent tissue-specific antigens (TSAs) (Aire-dep TSAs), Aire-independent TSAs (Aire-indep TSAs), and all TSAs were tested. NES indicates normalized enrichment score.

Fate mapping study with in vivo BrdU pulse-labeling of Fucci thymic epithelial cells (TECs).

(A) Schematic procedure of in vivo BrdU pulse labeling of Fucci mice, and analysis of BrdU staining in mCherryhiCD80hi and mCherryloCD80hi medullary TECs (mTECs) by flow cytometiric analysis. BrdU staining was done after sorting each cell fraction. (B) Typical flow cytometric profile of BrdU staining in mCherryloCD80hi mTECs (upper panels) and mCherryhiCD80hi mTECs (lower panels) at days 0.5, 3.5, and 6.5 after the BrdU injection. Data for the ratio of BrdU+ cells in each mTEC fraction are summarized in right-hand figures. N = 7 for 0.5 day after the BrdU injection, N = 3 for 3.5 and 6.5 days after the injection. Two-tailed Student’s t-tests. **p < 0.01 and *p < 0.05. NS, not significant (p > 0.05). p = 1.5 × 10–3 for the upper figure and p = 0.033 for the lower figure. Original data were shown in Figure 6—source data 1. (C) Cell number of BrdU+mCherryloCD80hi mTECs and BrdU+mCherryhiCD80hi mTECs at days 0.5, 3.5, and 6.5 after the BrdU injection. Two-tailed Student’s t-tests. **p < 0.01. NS, not significant (p > 0.05). p = 4.3 × 10–3 for the left figure and p = 5.1 × 10–3 for the right figure. Original data were shown in Figure 6—source data 1. (D) Mean fluorescence intensity (MFI) of BrdU staining in mCherryloCD80hi at day 0.5 and mCherryhiCD80hi at days 3.5 and 6.5. MFIs of other time points were difficult to evaluate because of very low cell numbers. Two-tailed Student’s t-tests. *p = 0.015 and **p = 6.5 × 10–3. NS, not significant (p > 0.05). Original data were shown in Figure 6—source data 1.

Figure 7 with 3 supplements
Fate mapping study of proliferating Aire+ medullary thymic epithelial cells (mTECs) in in vitro reaggregated thymic organ culture (RTOC).

(A) RTOC experiment to test the differentiation capacity of proliferating Aire+ mTECs. Proliferating Aire+ mTECs (mCherrylo) and E15.5 embryonic thymic cells were reaggregated and subsequently cultured for 5 days. Reaggregated thymic organ (RTO) was analyzed by flow cytometry. Representative flow cytometric profiles of RTOC are shown. N = 5. The ratio of mCherryhi cells in TECs is summarized in the right-hand figure. *p < 0.05. p = 0.027 between CD80hi and CD80lo in mCherrylo and p = 0.024 between CD80hi mCherrylo and CD80hi RTOC control. (B) Volcano plots of RNA-seq data from mCherrylo CD80hi mTECs (mCherrylo), mCherryhi CD80hi mTECs (mCherryhi), and mCherryhi CD80hi mTECs in RTOC (mCherryhi in RTOC). Red dots in volcano plots indicate genes for which expression differed significantly between the two subsets. Numbers of differentially expressed genes are shown in the panels. N = 3. Y axis is log10 of FDR p-value. (C) Expression levels of Aire and Mki67 in mCherrylo, mCherryhi, and mCherryhi in RTOC. (D) Scatter plot of normalized expression values of TA-TEC marker candidates in mCherrylo and mCherryhi in RTOC. TA-TEC marker candidate genes were selected from bulk RNA-seq data and scRNA-seq data in Figure 7—figure supplement 2 (E) Integration of well-based single-cell random displacement amplification sequencing (scRamDA-seq) data (mCherrylo, mCherryhi, and mCherryhi in RTOC) with the droplet-based scRNA-seq data in Figure 2. (F) Frequency of each cell cluster in scRamDA-seq data of mCherrylo, mCherryhi, and mCherryhi -RTOC. (G) Volcano plot of tissue-specific antigen (TSA) expression in each cell cluster in scRamDa-seq data of mCherryhi-RTOC as compared to mCherrylo. Red dots indicate significantly changed TSA genes.

Figure 7—figure supplement 1
Flow cytometric analysis and gene set enrichment analysis of reaggregation thymic organ culture experiments.

(A) Flow cytometric analysis for purity check for sorted mCherrylo medullary thymic epithelial cells (mTECs). (B) Ratio of Venus+ cells in sorted mCherrylo and RTOC. (C) Reaggregated thymic organ culture (RTOC) experiment to test the survival of contaminated mCherryhi mTECs. Sorted mCherryhi mTECs (200, 500, or 1000 cells) and E15.5 embryonic thymic cells were reaggregated and subsequently cultured for 5 days. Reaggregated thymic organ was analyzed by flow cytometer. (D) Gene set enrichment analysis (GSEA) of differentially expressed genes between mCherryhi and mCherrylo for pro-apoptotic gene set (HALLMARK_APOPTOSIS, M5902) (Liberzon et al., 2015). (E) Flow cytometric analysis of RTOC using allogenic fetal thymus (Balb/cA).

Figure 7—figure supplement 2
Analysis of differentially expressed tissue specific antigen (TSA) genes between mCherryhi thymic epithelial cells (TECs) in reaggregated thymic organ culture (RTOC) and mCherrylo TECs.

(A) Gene set enrichment analyses of differentially expressed genes between mCherryhi thymic epithelial cells (TECs) in reaggregated thymic organ culture (RTOC) and mCherrylo TECs (left panels), mCherryhi TECs and mCherrylo TECs (middle panels), and mCherryhi TECs and mCherryhi TECs in RTOC are exhibited. Gene sets of Aire-dependent tissue-specific antigens (TSAs) (Aire-dep TSAs), and Aire-independent TSAs (Aire-indep TSAs) were tested. NES indicated normalized enrichment score. FDR q indicated the false discovery rate q value. (B) Venn diagram of gene lists expressed in proliferating TECs at higher level in the three different RNA-seq datasets (mCherrylo vs. mCherryhi in Figure 6, Venus + vs. Venus in Figure 4 and cluster R1 in Figure 2). TA-TEC gene candidates were selected from the Venn diagram. TA-TEC marker gene candidates were selected by omitting cell cycle-related genes (GO:0007049 and Tirosh et al., 2016) form the TA-TEC gene candidates. The list of genes is summarized in Supplementary file 4. (C) Volcano plot for TA-TEC marker candidate expression in mCherrylo and mCherryhi in RTOC.

Figure 7—figure supplement 3
Integration of droplet-based single-cell RNA sequencing and well-based single-cell random displacement amplification sequencing data.

(A) Uniform manifold approximation and production (UMAP) plot of droplet-based single-cell RNA sequencing (scRNA-seq) and well-based single-cell random displacement amplification sequencing (scRamDA-seq) data after their integration. Cell clusters are separated by colors and numbers in the plot. The graph on the right shows the percentages of each cluster in the total number of cells detected. Each cluster was assigned based on gene expression profile and corresponded with clusters in Figure 2. (B) Violin plots depicting expression level of typical thymic epithelial cell (TEC) marker genes in each cluster.

Figure 8 with 2 supplements
Proliferating Aire+ CD80hi medullary thymic epithelial cells (mTECs) persist in older mice.

(A) Flow cytometry analysis of CD80hi mTEC subsets from Fucci2a mice aged 4, 8, and 19 weeks. Representative data are shown. Percentages of Venus+ cells in CD80hi mTEC subsets are summarized in the graph in the right panel. N = 4 each for Airegfp/+:: Fucci2a (circle) and control::Fucci2a (closed triangles). (B) Schematic depiction of the proposed process of Aire+ mTEC development in the adult thymus. Transit-amplifying TSAlo Aire+ TECs give rise to mature mTECs. Precursor cells to the transit-amplifying TECs have not been determined yet. Cluster numbers in Figure 1 are shown together with the model of mTEC subsets I to IV.

Figure 8—figure supplement 1
Integration of single-cell RNA sequencing data of 4-week-old and fetal thymic epithelial cell.

(A) Single-cell RNA sequencing (scRNA-seq) data in this study (4-week-old mice) were integrated with scRNA-seq data reported by others (Kernfeld et al., 2018) by Seurat package (Find integration anchors function). (B) Expression levels of Mki67 (upper panels) and Aire (lower panels) in the subcluster are shown in violin plots (left) and dot plots (right). (C) scRNA-seq data in this study (4-week-old mice) were integrated with scRNA-seq data reported by others (Kernfeld et al., 2018) by Harmony algorithm (Korsunsky et al., 2019).

Figure 8—figure supplement 2
Single-cell RNA sequencing (scRNA-seq) data of feta thymus (Kernfeld et al., 2018) were integrated with other adult scRNA-seq data (Bornstein et al., 2018; Dhalla et al., 2020; Wells et al., 2020).

Tables

Appendix 1—key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Genetic reagent (Mus musculus)B6.Cg-Aire < tm2Mmat>/RbrcRIKEN BioResource Research CenterBRC No:RBRC03515
Genetic reagent (Mus. musculus)B6;129-Gt(ROSA)26Sor < tm1(Fucci2aR)Jkn>RIKEN BioResource Research CenterBRC No:RBRC06511
Genetic reagent (Mus. musculus)B6(Cg)-Foxn1tm3(cre)Nrm/JJackson LaboratoryIMSR Cat#JAX:018448, RRID:IMSR_JAX:018448
Genetic reagent (Mus. musculus)CAG-Cre transgenic miceProvided by Jun-ichi Miyazaki
AntibodyPurified anti-mouse CD16/32(Rat monoclonal)BioLegendCat#101302, RRID:AB_312801FACS(1:200)
AntibodyAPC/Cyanine7 anti-mouse CD45(Rat monoclonal)BioLegendCat#103116, RRID:AB_312981FACS (1:200)
AntibodyPE Rat anti-mouse CD45(Rat monoclonal)eBioscienceCat#12-0451-82, RRID:AB_465668FACS (1:200)
AntibodyAPC/Cyanine7 anti-mouse TER-119/Erythroid Cells(Rat monoclonal)BioLegendCat#116223, RRID:AB_2137788FACS (1:200)
AntibodyPE anti-mouse TER-119/Erythroid Cells(Rat monoclonal)eBioscienceCat#12-5921-82, RRID:AB_466042FACS (1:200)
AntibodyBrilliant Violet 510 anti-mouse CD326 (Ep-CAM)(Rat monoclonal)BioLegendCat#118231, RRID:AB_2632774FACS (1:400)
AntibodyFITC anti-mouse CD326 Ep-CAM (Rat monoclonal)BioLegendCat#118208, RRID:AB_1134107FACS (1:400)
AntibodyAlexa Fluor 647 anti-mouse Ly-51 (Rat monoclonal)BioLegendCat#108312, RRID:AB_2099613FACS (1:400)
Chemical compound, drugBiotinylated Ulex Europaeus Agglutinin I (UEA I)Vector LaboratoriesCat#B-1065–2FACS (1:800)
Chemical compound, drugStreptavidin PE/Cyanine7 ConjugateeBioscienceCat#25-4317-82FACS (1:800)
Chemical compound, drugStreptavidin APC/Cyanine7 ConjugateBD PharmingenCat#554063 RRID:AB_10054651FACS (1:400)
AntibodyPE anti-mouse CD80 (Armenian hamster monoclonal)eBioscienceCat#12-0801-81, RRID:AB_465751FACS (1:300)
AntibodyPacific Blue anti-mouse CD80 Antibody (Armenian hamster monoclonal)BioLegendCat#104724, RRID:AB_2075999FACS (1:300)
AntibodyAlexa Fluor 647 anti-mouse Aire (Rat monoclonal)eBioscienceCat#51-5934-80IHC (1:100)
AntibodyPurified Rabbit anti-Keratin 5 (rabbit polyclonal)BioLegendCat#905504, RRID:AB_2616956IHC (1:400)
AntibodyAlexa Fluor 647 anti-Rabbit IgG (H + L)(Donkey polyclonal)InvitrogenCat#A-31573, RRID:AB_2536183IHC (1:1000)
Chemical compound, drugLiberase TMRoche DiagnosticsCat#5401127001
Chemical compound, drug7-Aminoactinomycin DCalbiochemCat#129935-1MGCN
Chemical compound, drugSYTOX Blue Nucleic Acid StainInvitrogenCat#S11348
Software, algorithmFlowJo version 10BDFlowJo, RRID:SCR_008520
Software, algorithmCell Ranger v3.0.010× GenomicsCell Ranger, RRID:SCR_017344
Software, algorithmSEURAT version 4.1.0https://github.com/satijalab/seurat/blob/master/vignettes/install.RmdSEURAT, RRID:SCR_007322
Software, algorithmVelocyto version 0.6https://github.com/velocyto-team/velocyto.RVelocyto, RRID:SCR_018167
Software, algorithmpagoda2 version 1.0.9https://github.com/kharchenkolab/pagoda2pagoda2, RRID:SCR_017094
Software, algorithmCell Ranger ATAC version1.1.010× GenomicsCell Ranger ATAC, RRID:SCR_021160
Software, algorithmSignac version 1.5.0https://github.com/timoast/signac/blob/master/vignettes/install.RmdSignac, RRID:SCR_021158
Software, algorithmMonocle3, version 0.2.3https://cole-trapnell-lab.github.io/monocle3/docs/installation/Monocle3, RRID:SCR_018685
Software, algorithmCLC Genomics Workbench Version 7.5.1QIAGENCLC Genomics Workbench, RRID:SCR_011853
Commercial assay or kitChromium Single Cell 3’ Library & Gel Bead Kit v210× GenomicsCat#PN-120237
Commercial assay or kitChromium Single Cell A Chip Kit10× GenomicsCat#PN-120236
Commercial assay or kitChromium i7 Multiplex Kit10× GenomicsCat#PN-120262
Commercial assay or kitChromium Next GEM Single Cell ATAC Library & Gel Bead Kit10× GenomicsCat#PN-1000176
Commercial assay or kitChromium Next GEM Chip H Single Cell Kit10× GenomicsCat#PN-1000162
Commercial assay or kitSingle Index Kit N, Set A10× GenomicsCat#PN-1000212
Commercial assay or kitNEBNext rRNA Depletion KitNew England BiolabsCat#E6310
Commercial assay or kitNEBNext Ultra Directional RNA Library Prep Kit for IlluminaNew England BiolabsCat#E7420
Commercial assay or kitKAPALibraryQuantificationKits Illumina/UniversalNippon GeneticsCat#KK4824
Chemical compound, drugKAPAHiFi DNA PolymeraseNippon GeneticsCat#KK2102
Commercial assay or kitAgilent High Sensitivity DNA KitAgilent TechnologiesCat#5067–4626
Commercial assay or kitMultina DNA-12000SHIMADZUCat#S292-36600-91

Additional files

Supplementary file 1

List of genes specifically expressed in each cluster.

https://cdn.elifesciences.org/articles/73998/elife-73998-supp1-v2.xlsx
Supplementary file 2

Percentage of single-cell RNA sequencing (scRNA-seq) clusters in single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) clusters after the integration.

https://cdn.elifesciences.org/articles/73998/elife-73998-supp2-v2.xlsx
Supplementary file 3

Differentially chromatin-accessible regions between cluster 0 and cluster 4.

https://cdn.elifesciences.org/articles/73998/elife-73998-supp3-v2.xlsx
Supplementary file 4

Percentage of single-cell RNA sequencing (scRNA-seq) subclusters of R1 in single-cell assays for transposase-accessible chromatin sequencing (scATAC-seq) clusters after the integration.

https://cdn.elifesciences.org/articles/73998/elife-73998-supp4-v2.xlsx
Supplementary file 5

Gene ontology (GO) analysis of genes differentially expressed in Venus + cells.

https://cdn.elifesciences.org/articles/73998/elife-73998-supp5-v2.xlsx
Supplementary file 6

List of all, Aire-dependent, Aire-independent tissue-specific antigen genes.

https://cdn.elifesciences.org/articles/73998/elife-73998-supp6-v2.xlsx
Supplementary file 7

Summary for assignment of individual single cells in single-cell random displacement amplification sequencing (scRamDa-seq) of mCherryhi, mCherrylo, and mCherryhi-RTOC.

https://cdn.elifesciences.org/articles/73998/elife-73998-supp7-v2.xlsx
Supplementary file 8

Summary for assignment of individual single cells in single-cell random displacement amplification sequencing (scRamDa-seq) of mCherryhi, mCherrylo, and mCherryhi-RTOC.

https://cdn.elifesciences.org/articles/73998/elife-73998-supp8-v2.xlsx
Transparent reporting form
https://cdn.elifesciences.org/articles/73998/elife-73998-transrepform1-v2.docx

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. Takahisa Miyao
  2. Maki Miyauchi
  3. S Thomas Kelly
  4. Tommy W Terooatea
  5. Tatsuya Ishikawa
  6. Eugene Oh
  7. Sotaro Hirai
  8. Kenta Horie
  9. Yuki Takakura
  10. Houko Ohki
  11. Mio Hayama
  12. Yuya Maruyama
  13. Takao Seki
  14. Hiroto Ishii
  15. Haruka Yabukami
  16. Masaki Yoshida
  17. Azusa Inoue
  18. Asako Sakaue-Sawano
  19. Atsushi Miyawaki
  20. Masafumi Muratani
  21. Aki Minoda
  22. Nobuko Akiyama
  23. Taishin Akiyama
(2022)
Integrative analysis of scRNA-seq and scATAC-seq revealed transit-amplifying thymic epithelial cells expressing autoimmune regulator
eLife 11:e73998.
https://doi.org/10.7554/eLife.73998