Fig 1: scATAC-seq, scRNA-seq and scTCR-seq analysis of the changes in T cell subsets in GO.a UMAP projections of T cell subsets of scATAC-seq dataset. b Differences in the proportions of Th1, Th17, CD8 Te and CD4 CTL among HC (n = 10) and GO (n = 21). Data are represented as the median [IQR]. *P < 0.05 by Mann-Whitney U test. c Proportion of each T subsets cell type in scRNA-seq after processing with the Milo algorithm. d UpSet plot showing the integrated comparative analysis of upregulated DARs in effector T cells of the GO group compared to HC. e UpSet plot showing the integrated comparative analysis of upregulated DEGs in effector T cells of the GO group compared to HC. f Enriched gene oncology terms in upregulated DARs of T cell subsets in the GO group. g Enriched gene oncology terms in upregulated DEGs of T cell subsets in the GO group. h Bar chart illustrating the proportion of unique TCR clones in HC and GO groups. Data are represented as median [IQR]. *P < 0.05 by Mann-Whitney U test. i Bar chart showing expression abundance of GO-specific clones across T cell subsets in the GO group. j UMAP plot mapping GO-specific clones in T cells. Contour plots (k) and bar chart (l) indicating CD8 + GNLY + KLRC1- cell proportions from HC and GO donors. m Immunofluorescence images showing the infiltration of CD8 + GNLY+ cells in the orbits of HC and GO donors (left). Bar chart illustrating the difference in the proportion of CD8 + GNLY+ cells infiltration in HC and GO groups (Right). Data are represented as the median [IQR]. *P < 0.05 by Mann-Whitney U test. scale bar: 20 μm.