Investigation of the cell composition and gene expression in the delayed-type hypersensitivity tuberculin skin test
Military Medical Research volume 10, Article number: 14 (2023)
The tuberculin skin test (TST) reagents have continuously improved, with the ESAT6-CFP10 (EC) test having recently been introduced, but are seldom based on the direction of the delayed-type hypersensitivity (DTH) mechanism. Previous studies only partially showed the infiltration and activation of immune cells and the production of cytokines of the skin induration [1, 2], and lack the detailed measurements of cell proportions and gene expression in the DTH response. Therefore, in this study, we revealed the comprehensive characteristics of DTH by single-cell RNA sequencing (scRNA-seq) in the guinea pig tuberculosis (TB) model [Experimental Animal Welfare Ethics Committee, Beijing Tuberculosis and Thoracic Tumor Research Institute (2021–064)].
A denser cell infiltrate was stained in the dermis after exposure to EC or purified protein derivative (PPD) than after exposure to phosphate buffered saline (PBS), with the densest infiltrate observed after exposure to PPD. The dermal thickness after PPD exposure [(1768.3 ± 130.9) μm] was significantly higher than that after PBS [(1442.94 ± 77.73) μm] or EC [(1477.52 ± 95.93) μm] exposure (P < 0.01, Additional file 1: Fig. S1). Moreover, the area of redness and swelling induced by EC in skin induration was broader than that induced by PPD, consistent with previous studies . These results verified the difference in the DTH responses induced by EC and PPD.
To better understand the DTH induced by EC and PPD, scRNA-seq was performed with skin induration (Additional file 1: Methods). After stringent quality control (Additional file 1: Fig. S2), the transcriptomes of 33,074 cells were obtained and analyzed (Fig. 1a, Additional file 1: Table S1). Then, through SingleR annotation, 10 different major cell types were identified (Fig. 1b, Additional file 1: Fig. S3), among which the controversial cell type of hepatocytes was identified due to the limitation of ImmGen databases being based on mice rather than guinea pigs. Based on established marker genes (Fig. 1c), with a definite correlation between different clusters (Additional file 1: Fig. S4), the cells in the merged dataset were identified as fibroblasts, endothelial cells (Endos), monocytes, T cells, macrophages, natural killer (NK) cells, B cells, granulocytes and dendritic cells (DCs), and the composition of the main cell types, such as T cells, in the EC induration samples differed from that in the PPD induration samples (Fig. 1d, Additional file 1: Fig. S5). Although the reactivity of Endos, NK cells, and B cells in guinea pigs differed slightly from that of the corresponding human cells in the induration samples, most were consistent (Fig. 1b), indicating that the guinea pig is a suitable TST model, with skin architecture, physiology, and immunopathology similar to those in human beings .
Moreover, we identified 1615 up-regulated and 1148 down-regulated differentially expressed genes (DEGs) in the EC sample, and 1669 up-regulated and 1184 down-regulated DEGs in the PPD sample (Fig. 1e, Additional file 1: Fig. S6, Additional file 2). Most of the DEGs were common to the EC and PPD samples (Fig. 1f), and 110 up-regulated DEGs were common to at least four identified cell types (Additional file 3). These DEGs were related mainly to the adaptive immune and immune response pathway (Additional file 1: Figs. S7, S8 and S9, Additional file 4), indicating the complex immune interactions during the DTH response between different cell types.
Then, T cells from all samples were clustered unsupervised based on the specific markers as previous reports  and 7 cell subpopulations were identified (Additional file 1: Fig. S10): Th1 cells, Tregs, effector-memory CD8+ T (CD8_EM) cells, naive CD4+ T-like (CD4_NaiveLike) cells, naive CD8+ T-like (CD8_NaiveLike) cells, CD8_EarlyActivation cells and “unknown cells” (unable to be characterized based on specific marker genes) (Fig. 1g, h). The ratio of CD4+ to CD8+ cells was 2 to 1 in the EC sample, but 5 to 1 in the PPD sample (Fig. 1i). The functions of the DEGs were also analyzed (Additional file 1: Fig. S10).
The Endos were significantly increased in the EC and PPD groups (Fig. 1d), and could be categorized into 15 clusters (Fig. 1j–l, Additional file 1: Fig. S11). The functions of the DEGs are analyzed (Additional file 1: Figs. S12 and S13). Moreover, various cytokines were up-regulated in Endos, suggesting that Endos may play more important roles in DTH than previously believed (Additional file 1: Fig. S6).
In conclusion, the heterogeneity of cells and their transcriptomes involved in DTH were first clarified by scRNA-seq in the current study, which determined the immune microenvironment during the DTH response (Additional file 1: Fig. S14). Deciphering the reactivity of different types of cells in DTH could robustly guide the improvement of TST reagents. The most suitable components of the ideal TST reagent are achievable, for example, EC instead of the complex PPD used in TB diagnosis.
Availability of data and materials
The datasets and materials used in this study are available from the corresponding author on reasonable request.
Differentially expressed genes
- NK cells:
Natural killer cells
Phosphate buffered saline
Purified protein derivative
Single-cell RNA sequencing
Tuberculin skin test
Tomlinson GS, Cashmore TJ, Elkington PT, Yates J, Lehloenya RJ, Tsang J, et al. Transcriptional profiling of innate and adaptive human immune responses to mycobacteria in the tuberculin skin test. Eur J Immunol. 2011;41(11):3253–60.
Czarnobilska E, Obtułowicz K, Wsołek K. Type IV of hypersensitivity and its subtypes. Przegl Lek. 2007;64(7–8):506–8.
Abubakar I, Jackson C, Rangaka MX. C-Tb: a latent tuberculosis skin test for the 21st century? Lancet Respir Med. 2017;5(4):236–7.
Orme IM, Ordway DJ. Mouse and guinea pig models of tuberculosis. Microbiol Spectr. 2016. https://doi.org/10.1128/microbiolspec.TBTB2-0002-2015.
Andreatta M, Corria-Osorio J, Müller S, Cubas R, Coukos G, Carmona SJ. Interpretation of T cell states from single-cell transcriptomics data using reference atlases. Nat Commun. 2021;12(1):2965.
The authors thank all the coordinators and technicians for their hard field work.
This work was supported by the National Natural Science Foundation (81871691), and the Beijing Municipal Natural Science Foundation (KZ202110025034).
Ethics approval and consent to participate
This study was approved by the Experimental Animal Welfare Ethics Committee, Beijing Tuberculosis and Thoracic Tumor Research Institute (2021-064).
Consent for publication
The authors declare they have no competing interests.
. Methods. Fig. S1 Characteristics of skin induration. Fig. S2 Quality control of all cells acquired by scRNA-seq. Fig. S3 Cellular diversity and heterogeneity of skin induration. Fig. S4 Correlation analysis on cell clusters and distribution of cell population. Fig. S5 Cellular diversity of skin induration. Fig. S6 DEGs among different cell types caused by EC and PPD. Fig. S7 Functions of up-regulated DEGs shared by EC and PPD. Fig. S8 Functions of down-regulated DEGs shared by EC and PPD. Fig. S9 GO and KEGG enrichment analysis of DEGs in EC and PPD samples. Fig. S10 Expression patterns and transcriptional regulatory networks of T cells. Fig. S11 DEGs in Endos in DTH. Fig. S12 GO and KEGG terms of up-regulated DEGs in Endos. Fig. S13 GO and KEGG terms of down-regulated DEGs in Endos. Fig. S14 Scheme of DTH stimulated by different antigens. Table S1 Results of quality control.
. Up-regulated and down-regulated DEGs in total cells in EC or PPD
. Up-regulated DEGs shared by EC and PPD expressed in at least four cell types
. GO and KEGG terms enriched with up-regulated and down-regulated DEGs in EC or PPD
About this article
Cite this article
Duan, HJ., Chu, HQ., Cao, TM. et al. Investigation of the cell composition and gene expression in the delayed-type hypersensitivity tuberculin skin test. Military Med Res 10, 14 (2023). https://doi.org/10.1186/s40779-023-00450-2