Single-cell differential splicing analysis reveals high heterogeneity of liver tumor-infiltrating T cells

Shang Liu,Biaofeng Zhou,Liang Wu,Yan Sun,Jie Chen,Shiping Liu

DOI: https://doi.org/10.1038/s41598-021-84693-w

IF: 4.6

2021-03-05

Scientific Reports

Abstract:Abstract Recent advances in single-cell RNA sequencing (scRNA-seq) have improved our understanding of the association between tumor-infiltrating lymphocyte (TILs) heterogeneity and cancer initiation and progression. However, studies investigating alternative splicing (AS) as an important regulatory factor of heterogeneity remain limited. Here, we developed a new computational tool, DESJ-detection, which accurately detects differentially expressed splicing junctions (DESJs) between cell groups at the single-cell level. We analyzed 5063 T cells of hepatocellular carcinoma (HCC) and identified 1176 DESJs across 11 T cell subtypes. Interestingly, DESJs were enriched in UTRs, and have putative effects on heterogeneity. Cell subtypes with a similar function closely clustered together at the AS level. Meanwhile, we identified a novel cell state, pre-activation with the isoform markers ARHGAP15-205. In summary, we present a comprehensive investigation of alternative splicing differences, which provided novel insights into T cell heterogeneity and can be applied to other full-length scRNA-seq datasets.

multidisciplinary sciences

What problem does this paper attempt to address?

The main problem this paper attempts to address is the study of alternative splicing heterogeneity in tumor-infiltrating T cells (TILs) in hepatocellular carcinoma (HCC) using single-cell RNA sequencing (scRNA-seq) technology. Specifically, the paper develops a new computational tool—DESJ-detection, which is used to accurately detect differentially expressed splicing junctions (DESJs) at the single-cell level. By analyzing data from 5063 T cells, researchers identified 1176 differentially spliced junctions and found that these junctions are enriched in untranslated regions (UTRs), potentially impacting T cell heterogeneity. Additionally, the paper identifies a new cell state—pre-activation state, with ARHGAP15-205 as a marker. In summary, this study provides new insights into T cell heterogeneity and can be applied to other full-length scRNA-seq datasets.