CoChat
Sign Up
Cited by 343T-cell and B-cell ImmunologySingle-cell and spatial transcriptomicsGenetics, Aging, and Longevity in Model Organisms

The DOI matches its Crossref record; 51 of the 51 references it cites were found in a registry or an index.

Single-cell transcriptomics reveals expansion of cytotoxic CD4 T cells in supercentenarians

Open access
Hashimoto, Kosuke, Kouno, Tsukasa, Ikawa, Tomokatsu, Hayatsu, Norihito, Miyajima, Yurina, Yabukami, Haruka, Terooatea, Tommy, Sasaki, Takashi, Suzuki, Takahiro, Valentine, Matthew, Pascarella, Giovanni, Okazaki, Yasushi, Suzuki, Harukazu, Shin, Jay W., Minoda, Aki, Taniuchi, Ichiro, Okano, Hideyuki, Arai, Yasumichi, Hirose, Nobuyoshi, Carninci, Piero·2019·Proceedings of the National Academy of Sciences
DOI 10.1073/pnas.1907883116PMC PMC6883788

Published version of a preprint

Open this paper in CoChat — free account

Abstract

Supercentenarians, people who have reached 110 y of age, are a great model of healthy aging. Their characteristics of delayed onset of age-related diseases and compression of morbidity imply that their immune system remains functional. Here we performed single-cell transcriptome analysis of 61,202 peripheral blood mononuclear cells (PBMCs), derived from 7 supercentenarians and 5 younger controls. We identified a marked increase of cytotoxic CD4 T cells (CD4 cytotoxic T lymphocytes [CTLs]) as a signature of supercentenarians. Furthermore, single-cell T cell receptor sequencing of 2 supercentenarians revealed that CD4 CTLs had accumulated through massive clonal expansion, with the most frequent clonotypes accounting for 15 to 35% of the entire CD4 T cell population. The CD4 CTLs exhibited substantial heterogeneity in their degree of cytotoxicity as well as a nearly identical transcriptome to that of CD8 CTLs. This indicates that CD4 CTLs utilize the transcriptional program of the CD8 lineage while retaining CD4 expression. Indeed, CD4 CTLs extracted from supercentenarians produced IFN-γ and TNF-α upon ex vivo stimulation. Our study reveals that supercentenarians have unique characteristics in their circulating lymphocytes, which may represent an essential adaptation to achieve exceptional longevity by sustaining immune responses to infections and diseases.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1073/pnas.1907883116

Cites

51 of 51 references verified0 citation errors · 0 not checkable

Statistics

Figures

Fig. 1. Single-cell transcriptome profiling of PBMCs of supercentenarians and controls. (A) Schematic representation of single-cell transcriptome experiments, from blood sample collection to visualization. (B) The number of recovered cells that passed quality control and the median number of genes per cell for each of the donors (7 supercentenarians, SC1−SC7; and 5 controls, CT1–CT5). (C) Two-dimensional tSNE visualization of PBMCs for supercentenarians (Left) and controls (Right). Different colors represent 10 clusters (cell types) defined by the k-means clustering algorithm. (D) Expression of marker genes for 6 major cell types; cell positions are from the tSNE plot in C.
Fig. 1. Single-cell transcriptome profiling of PBMCs of supercentenarians and controls. (A) Schematic representation of single-cell transcriptome experiments, from blood sample collection to visualization. (B) The number of recovered cells that passed quality control and the median number of genes per cell for each of the donors (7 supercentenarians, SC1−SC7; and 5 controls, CT1–CT5). (C) Two-dimensional tSNE visualization of PBMCs for supercentenarians (Left) and controls (Right). Different colors represent 10 clusters (cell types) defined by the k-means clustering algorithm. (D) Expression of marker genes for 6 major cell types; cell positions are from the tSNE plot in C.CC BY · PubMed Central
Fig. 2. Significant reduction of B cells in supercentenarians. (A) Boxplots of the percentage of each cell type (defined by single-cell RNA-Seq) in PBMCs of 7 supercentenarians (SC1–SC7) and 5 controls (CT1–CT5)—the boxes extend from the 25th to 75th percentile and encompass the median (horizontal line). BC, B cell; TC, T cell; NK, natural killer cell; M14, CD14+ monocyte. *P < 0.05 (Wilcoxon rank sum test); no asterisk means not significant. (B) Representative FACS plots showing CD19+ B cells; the plots for other donors are shown in SI Appendix, Fig. S2B. (C) Boxplots of the percentage of each cell type (defined by FACS) in PBMCs of 4 supercentenarians SC1–SC4 and 3 controls CT1–CT3. No asterisk means not significant (Wilcoxon rank sum test).
Fig. 2. Significant reduction of B cells in supercentenarians. (A) Boxplots of the percentage of each cell type (defined by single-cell RNA-Seq) in PBMCs of 7 supercentenarians (SC1–SC7) and 5 controls (CT1–CT5)—the boxes extend from the 25th to 75th percentile and encompass the median (horizontal line). BC, B cell; TC, T cell; NK, natural killer cell; M14, CD14+ monocyte. *P < 0.05 (Wilcoxon rank sum test); no asterisk means not significant. (B) Representative FACS plots showing CD19+ B cells; the plots for other donors are shown in SI Appendix, Fig. S2B. (C) Boxplots of the percentage of each cell type (defined by FACS) in PBMCs of 4 supercentenarians SC1–SC4 and 3 controls CT1–CT3. No asterisk means not significant (Wilcoxon rank sum test).CC BY · PubMed Central
Fig. 3. Expansion of cytotoxic T cells in supercentenarians. (A) Boxplots of percentages of TC1 and TC2 T cells (defined by k-means clustering of single cell RNA-Seq data) in PBMCs of 7 supercentenarians (SC1–SC7) and 5 controls (CT1–CT5). *P < 0.05 (Wilcoxon rank sum test). (B) Two-dimensional tSNE visualization of T cells using the Seurat R package. Different colors represent 2 clusters (Seurat_TC1 and Seurat_TC2), similar to the original TC1 and TC2 clusters. Right (Top and Bottom) show supercentenarians and controls, respectively. (C) Top 20 genes significantly highly expressed in Seurat_TC2 (Left) and Seurat_TC1 (Right). Major cytotoxic effector genes and lymph node homing markers are shown in red. (D) Expression of cytotoxic genes in supercentenarians (Top) and controls (Bottom); cell positions are from the tSNE plot in B. (E) Number of detected genes out of 4 cytotoxic genes (GZMH, GZMB, GZMA, and PRF1) per cell. (F) Percentage of cytotoxic T cells (cells clustered in TC2) among the total T cells. *P < 0.05 (Wilcoxon rank sum test).
Fig. 3. Expansion of cytotoxic T cells in supercentenarians. (A) Boxplots of percentages of TC1 and TC2 T cells (defined by k-means clustering of single cell RNA-Seq data) in PBMCs of 7 supercentenarians (SC1–SC7) and 5 controls (CT1–CT5). *P < 0.05 (Wilcoxon rank sum test). (B) Two-dimensional tSNE visualization of T cells using the Seurat R package. Different colors represent 2 clusters (Seurat_TC1 and Seurat_TC2), similar to the original TC1 and TC2 clusters. Right (Top and Bottom) show supercentenarians and controls, respectively. (C) Top 20 genes significantly highly expressed in Seurat_TC2 (Left) and Seurat_TC1 (Right). Major cytotoxic effector genes and lymph node homing markers are shown in red. (D) Expression of cytotoxic genes in supercentenarians (Top) and controls (Bottom); cell positions are from the tSNE plot in B. (E) Number of detected genes out of 4 cytotoxic genes (GZMH, GZMB, GZMA, and PRF1) per cell. (F) Percentage of cytotoxic T cells (cells clustered in TC2) among the total T cells. *P < 0.05 (Wilcoxon rank sum test).CC BY · PubMed Central
Fig. 4. Expansion of cytotoxic CD4 T cells in supercentenarians. (A) Classification of cytotoxic T cells into 3 subtypes—CD4 CTLs, CD8 CTLs, and γδ T cells—was based on the expression of CD4, CD8, and TRDC (see also SI Appendix, Fig. S4A) in T cells of 7 supercentenarians (SC1–SC7) and 5 controls (CT1–CT5); cell positions are from the tSNE plot in Fig. 3B. (B) Percentages of CD4 CTLs and γδ T cells among the total T cells. *P < 0.05 (Wilcoxon rank sum test); NS, not significant. (C) Percentages of CD4+ T cells and CD8+ T cells in total T cells. NS, not significant (Wilcoxon rank sum test). (D) FACS profiles of 6 supercentenarians (SC1, SC5–7, and SC9) and 1 semisupercentenarian (SC8). Cells gated on CD3+ were profiled using CD4 (x axis) and GZMB or IgG1 κ as an isotype control (y axis). Cells in Top Right corners are CD4 CTLs. (E) Percentages of CD4+ GZMB+ cells among the total T cells of the 6 supercentenarians and 1 semisupercentenarian listed in D and 5 controls (CT4, CT5, and CT6–CT8). *P < 0.05 (Wilcoxon rank sum test). (F) Correlation between percentages of CD4 CTLs determined by RNA-Seq and FACS measurements. Each dot represents 1 donor, shown in green for supercentenarians (SC1, SC5–SC7) and red for controls (CT4, CT5).
Fig. 4. Expansion of cytotoxic CD4 T cells in supercentenarians. (A) Classification of cytotoxic T cells into 3 subtypes—CD4 CTLs, CD8 CTLs, and γδ T cells—was based on the expression of CD4, CD8, and TRDC (see also SI Appendix, Fig. S4A) in T cells of 7 supercentenarians (SC1–SC7) and 5 controls (CT1–CT5); cell positions are from the tSNE plot in Fig. 3B. (B) Percentages of CD4 CTLs and γδ T cells among the total T cells. *P < 0.05 (Wilcoxon rank sum test); NS, not significant. (C) Percentages of CD4+ T cells and CD8+ T cells in total T cells. NS, not significant (Wilcoxon rank sum test). (D) FACS profiles of 6 supercentenarians (SC1, SC5–7, and SC9) and 1 semisupercentenarian (SC8). Cells gated on CD3+ were profiled using CD4 (x axis) and GZMB or IgG1 κ as an isotype control (y axis). Cells in Top Right corners are CD4 CTLs. (E) Percentages of CD4+ GZMB+ cells among the total T cells of the 6 supercentenarians and 1 semisupercentenarian listed in D and 5 controls (CT4, CT5, and CT6–CT8). *P < 0.05 (Wilcoxon rank sum test). (F) Correlation between percentages of CD4 CTLs determined by RNA-Seq and FACS measurements. Each dot represents 1 donor, shown in green for supercentenarians (SC1, SC5–SC7) and red for controls (CT4, CT5).CC BY · PubMed Central
Fig. 5. The differentiation state of T cells for 7 supercentenarians (SC1–SC7) and 5 controls (CT1–CT5). (A) Pseudotime trajectory of T cells estimated using Monocle 2. A continuous value from 0 to 12 was assigned to each cell as a pseudotime. The Bottom shows the general scheme of T cell differentiation. TN, naïve; TCM, central memory; TEM, effector memory; and TEMRA, effector memory reexpressing CD45RA. (B) Expression transition of differentiation-associated genes along the pseudotime. (C) Percentages of T cells along the pseudotime for supercentenarians (SC) and controls (CT). (D) Percentages of CD4 and CT8 CTLs among the total T cells along the pseudotime. (E) Correlation of gene expression between CD4 and CD8 CTLs. (F) Expression transition of selected genes shown separately for CD4 and CD8 CTLs.
Fig. 5. The differentiation state of T cells for 7 supercentenarians (SC1–SC7) and 5 controls (CT1–CT5). (A) Pseudotime trajectory of T cells estimated using Monocle 2. A continuous value from 0 to 12 was assigned to each cell as a pseudotime. The Bottom shows the general scheme of T cell differentiation. TN, naïve; TCM, central memory; TEM, effector memory; and TEMRA, effector memory reexpressing CD45RA. (B) Expression transition of differentiation-associated genes along the pseudotime. (C) Percentages of T cells along the pseudotime for supercentenarians (SC) and controls (CT). (D) Percentages of CD4 and CT8 CTLs among the total T cells along the pseudotime. (E) Correlation of gene expression between CD4 and CD8 CTLs. (F) Expression transition of selected genes shown separately for CD4 and CD8 CTLs.CC BY · PubMed Central
Fig. 6. Single-cell transcriptome and TCR profiles of CD4+ T cells for 2 supercentenarians (SC1, Top; and SC2, Bottom). (A) Schematic representation of experiments for the single-cell transcriptome and TCR analysis. (B) Two-dimensional tSNE visualization of 3 cell types (TC, T cell; NK, natural killer cell; EC, erythroid cell), and CD3D expression (Right). (C) Expression of marker genes for CD4 and CD8 T cells; cell positions are from the tSNE plot in B. (D) T cells recovered in both transcriptome and TCR libraries. Recovered cells were clustered into helper T cells and CD4 CTLs, shown in the tSNE plot. (E) Expression of a marker gene for cytotoxic T cells; cell positions are from the tSNE plot in D. (F) Diversity of TCRs in helper T cells and CD4 CTLs. (G) Cumulative occupancy of the top 50 most abundant clonotypes. (H) Occupancy of the top 10 most abundant clonotypes. (I) Pseudotime and cytotoxicity of clonally expanded CD4 T cells. Cytotoxicity values indicate the mean expression of 5 cytotoxic genes: NKG7, GZMA, GZMB, GZMH, and PRF1. (J) FACS profiles of 1 supercentenarian (SC2) and 1 young control (CT9) upon PMA/ionomycin stimulation. Cells gated on live CD3+ CD4+ CD8− were profiled using GZMB or IgG1 κ (x axis) and IFN-γ, TNF-α, or IgG1 κ (y axis).
Fig. 6. Single-cell transcriptome and TCR profiles of CD4+ T cells for 2 supercentenarians (SC1, Top; and SC2, Bottom). (A) Schematic representation of experiments for the single-cell transcriptome and TCR analysis. (B) Two-dimensional tSNE visualization of 3 cell types (TC, T cell; NK, natural killer cell; EC, erythroid cell), and CD3D expression (Right). (C) Expression of marker genes for CD4 and CD8 T cells; cell positions are from the tSNE plot in B. (D) T cells recovered in both transcriptome and TCR libraries. Recovered cells were clustered into helper T cells and CD4 CTLs, shown in the tSNE plot. (E) Expression of a marker gene for cytotoxic T cells; cell positions are from the tSNE plot in D. (F) Diversity of TCRs in helper T cells and CD4 CTLs. (G) Cumulative occupancy of the top 50 most abundant clonotypes. (H) Occupancy of the top 10 most abundant clonotypes. (I) Pseudotime and cytotoxicity of clonally expanded CD4 T cells. Cytotoxicity values indicate the mean expression of 5 cytotoxic genes: NKG7, GZMA, GZMB, GZMH, and PRF1. (J) FACS profiles of 1 supercentenarian (SC2) and 1 young control (CT9) upon PMA/ionomycin stimulation. Cells gated on live CD3+ CD4+ CD8− were profiled using GZMB or IgG1 κ (x axis) and IFN-γ, TNF-α, or IgG1 κ (y axis).CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Show all 51 references

Cited by

Sorted by citations

See all 343

Related

Metadata

Authors
Hashimoto, K., Kouno, T., Ikawa, T., Hayatsu, N., Miyajima, Y., Yabukami, H., Terooatea, T., Sasaki, T., Suzuki, T., Valentine, M., Pascarella, G., Okazaki, Y., Suzuki, H., Shin, J. W., Minoda, A., Taniuchi, I., Okano, H., Arai, Y., Hirose, N., & Carninci, P.
Affiliations
  • RIKEN Center for Integrative Medical Sciences (JP)
  • Keio University (JP)
Corresponding author
  • Hirose, Nobuyoshi
  • Carninci, Piero
ORCID
Venue
Proceedings of the National Academy of Sciences
Publisher
National Academy of Sciences
Published
2019-11-12
Volume
116
Issue
48
Pages
24242-24251
ISSN
  • 0027-8424 (ISSN-L)
  • 1091-6490
DOI
10.1073/pnas.1907883116
PMC
PMC6883788
OpenAlex
W2983248500
Cited by
343
References
51
Citations by year
2026: 45 · 2025: 49 · 2024: 56 · 2023: 59 · 2022: 51 · 2021: 52 · 2020: 28 · 2019: 3
FWCI
13.36
Type
journal article
Language
en
Open access
Licence
Locations
  • Proceedings of the National Academy of Sciences (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
Funders
  • Keio University
  • Ministry of Education, Culture, Sports, Science and Technology
  • Japan Society for the Promotion of Science (17H04214, 19KK0222, 17H04090, 19H05747, 19H04820, 18K19569)
  • RIKEN
Topics
  • Life Sciences › Immunology and Microbiology › Immunology › T-cell and B-cell Immunology
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Aging › Genetics, Aging, and Longevity in Model Organisms
Keywords
supercentenarians, healthy aging, cytotoxic CD4+ T cells, clonal expansion, T cell receptor repertoire, peripheral blood mononuclear cells, immunosenescence, longevity, interferon-gamma, TNF-α
MeSH terms
Adult, Aged, Aged, 80 and over, B-Lymphocytes, Cell Differentiation, Cells, Cultured, Humans, Interferon-gamma, Leukocytes, Mononuclear, Middle Aged, Tumor Necrosis Factor-alpha, CD4-Positive T-Lymphocytes, Case-Control Studies, Gene Expression Profiling, Single-Cell Analysis, Clonal Evolution

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review