Circulating immunological transcriptomic profile identifies DDX3Y and USP9Y on the Y chromosome as promising biomarkers for predicting response to PD-1/PD-L1 blockade

封锁 转录组 染色体 计算生物学 PD-L1 生物 免疫学 医学 免疫系统 遗传学 基因 免疫疗法 受体 基因表达
作者
Liting You,Zhaodan Xin,Feifei Na,Min Chen,Wen Yang,Jin Li,Jiajia Song,Ling Bai,Jianzhao Zhai,Xiaohan Zhou,Binwu Ying,Juan Zhou
出处
期刊:Chinese Medical Journal [Lippincott Williams & Wilkins]
标识
DOI:10.1097/cm9.0000000000003403
摘要

To the Editor: Immune checkpoint inhibitors (ICIs) have shown remarkable clinical responses; however, their efficacy remains limited to a small subset of patients. Peripheral blood mononuclear cells (PBMCs) offer an effective, accessible, and minimally invasive approach to assess tumor immune status and identify ICI responders. In this study, we aimed to elucidate the role of PBMC gene expression in ICI treatment response and prognosis. This study received approval from the Biomedical Ethics Committee of West China Hospital, Sichuan University (No. 2019 [1045]) and the requirement to obtain informed consent was waived. We enrolled patients with histologically and clinically confirmed solid tumors who were scheduled to undergo anti-programmed death 1 (PD-1)/programmed death ligand 1 (PD-L1) immunotherapy at West China Hospital of Sichuan University between May 1, 2020, and May 31, 2021 for the discovery and validation cohort [Supplementary Figure 1, https://links.lww.com/CM9/C240]. Of these, 55 were included in the discovery cohort and 98 in the validation cohort. According to Response Evaluation Criteria in Solid Tumours version 1.1 (RECIST version 1.1), patients with no disease progression (PD) or tumor-induced deaths within 6 months of anti-PD-1/PD-L1 treatment, including those with complete response (CR), partial response (PR), or stable disease (SD), were grouped as ICI responders. Patients experiencing PD within 6 months were classified as ICI non-responders. We identified 27 ICI responders and 28 non-responders in the discovery cohort and 47 responders and 51 non-responders in the validation cohort. No significant differences in sex, age, cancer type, clinical stage, treatment line, or treatment regimen were observed between responders and non-responders in either the discovery or validation cohort (P >0.05, Supplementary Table 1, https://links.lww.com/CM9/C240). Transcriptome sequencing of PBMCs in the discovery cohort (N = 55) identified a total of 206 differentially expressed genes (DEGs) (P <0.05, fold change [FC] >1.3 or <0.76), with 116 upregulated and 90 downregulated genes [Figure 1A, B]. Principal component analysis (PCA) based on these DEGs accurately distinguished ICI responders and non-responders as separate groups [Figure 1C]. Gene Ontology analyses revealed DEGs enrichment in the positive regulation of protein kinase B signaling, epithelial cell development, and external encapsulating structure organization [Supplementary Figure 2, https://links.lww.com/CM9/C240]. Kyoto Encyclopedia of Genes and Genomes pathway analysis revealed significant DEGs enrichment in the nitrogen metabolism, calcium, and Rap1 signaling pathways [Supplementary Figure 2, https://links.lww.com/CM9/C240]. Through reverse transcription-polymerase chain reaction (RT-PCR), we confirmed significant differences in DDX3Y, USP9Y, UTY, KDM5D, and RPS4Y1 expression, between responders and non-responders (P <0.05, Figure 1D, Supplementary Figure 3, https://links.lww.com/CM9/C240).Figure 1: Circulating immunological transcriptomic profile identifies Y chromosome genes as promising biomarkers for efficacy of PD-1/PD-L1 blockade. (A–C) DEGs determined by transcriptome sequencing between ICI-responsive and ICI-non-responsive patients in the discovery cohort. (A) Difference ranking chart of DEGs. (B) Volcano plot of DEGs. The dashed line indicates the screening threshold (P = 0.05, FC = 1.3). (C) PCA plot based on DEGs between ICI-responsive and ICI-non-responsive patients. (D–I) DDX3Y and USP9Y predict the therapeutic response and survival outcomes after PD-1/PD-L1 blockade. (D, G) Expression levels of DDX3Y and USP9Y between responsive and non-responsive patients in the discovery (D) and validation cohorts (G). Differences in survival outcomes based on different DDX3Y (E, H) and USP9Y (F, I) expression levels in the discovery and validation cohorts. Data are presented as mean ± standard deviation. Statistical significance in (D, G) was determined by a two-sided unpaired t-test. Survival data in (E, F, H, I) were analyzed by log-rank test. * P <0.05. CI: Confidence interval; DEGs: Differentially expressed genes; FC: Fold change; HR: Hazard ratio; ICI: Immune checkpoint inhibitor; PCA: Principal component analysis; PD-1: Programmed death 1; PD-L1: Programmed death ligand 1; PFS: Progression-free survival.Further, in an independent external validation cohort (N = 98), we demonstrated that DDX3Y and USP9Y could predict the therapeutic efficacy and survival outcomes of patients receiving PD-1/PD-L1 blockade. ICI responders exhibited higher expression levels of DDX3Y and USP9Y than non-responders [Figure 1D, G]. Both DDX3Y and USP9Y were associated with progression-free survival (PFS) in patients receiving PD-1/PD-L1 blockade. For DDX3Y, the median PFS was shorter in the DDX3Y-low group than in the DDX3Y-high group (DDX3Y low vs. high: 77.0 days vs. 207.0 days, log-rank P = 0.0304, HR [95% CI]: 2.137 [1.109–6.305] in the discovery cohort; DDX3Y low vs. high: 164.5 days vs. 591.0 days, log-rank P = 0.0151, HR [95% CI]: 2.055 [1.148–3.205] in the validation cohort) [Figure 1E, H]. For USP9Y, the median PFS was shorter in the USP9Y-low group than in the USP9Y-high group (USP9Y low vs. high: 94.0 days vs. 625.0 days, log-rank P = 0.0299, HR [95% CI]: 2.084 [1.084–4.109] in the discovery cohort; USP9Y low vs. high: 160.0 days vs. 375.0 days, log-rank P = 0.0083, HR [95% CI]: 2.045 [1.202–3.218] in the validation cohort) [Figure 1F, I]. Then, we conducted subgroup analyses to assess the potential of DDX3Y and USP9Y in predicting responses within specific subgroups. We identified DDX3Y as a predictor of ICI efficacy in various subgroups (log-rank P <0.05), with a shorter median PFS in the DDX3Y-low group than in the DDX3Y-high group in both the discovery and validation cohorts; however, this difference was significant, except for the NSCLC subgroup [Supplementary Figure 4, https://links.lww.com/CM9/C240]. For USP9Y, patients in the USP9Y-high group demonstrated improved survival outcomes after PD-1/PD-L1 blockade treatment compared to those in the USP9Y-low group across both the discovery and validation cohorts; however, this finding was significant, except in the NSCLC subgroup of the validation cohort [Supplementary Figure 4, https://links.lww.com/CM9/C240]. Overall, our findings strengthen the evidence supporting DDX3Y and USP9Y as predictors of PD-1/PD-L1 blockade efficacy. As DDX3Y and USP9Y are located on the Y chromosome, it is necessary to eliminate any sex-related bias. First, we confirmed no significant difference in sex-related distribution between responders and non-responders [Supplementary Table 1, https://links.lww.com/CM9/C240]. Subsequently, we conducted a subgroup analysis, specifically for male patients [Supplementary Figure 5, https://links.lww.com/CM9/C240], which revealed a significantly higher expression of DDX3Y in ICI-responder males than in non-responder males (P = 0.003, Supplementary Figure 5A, https://links.lww.com/CM9/C240), and individuals with high DDX3Y expression exhibited longer PFS than those with low DDX3Y expression (log-rank P = 0.0188, HR [95% CI]: 1.630 [1.100–2.683], Supplementary Figure 5C, https://links.lww.com/CM9/C240). Similarly, the expression levels of USP9Y were higher in responder males (P = 0.0009, Supplementary Figure 5B, https://links.lww.com/CM9/C240) than in non-responder males, and individuals with low USP9Y expression had shorter PFS than those with high USP9Y expression (log rank P = 0.0059, HR [95% CI]: 1.812 [1.188–2.702], Supplementary Figure 5D, https://links.lww.com/CM9/C240). These findings further support the potential roles of DDX3Y and USP9Y as reliable predictors of the PD-1/PD-L1 blockade. The RNA-sequencing results from the discovery cohort revealed a significant number of DEGs located on the Y chromosome, including DDX3Y, USP9Y, UTY, KDM5D, and RPS4Y1. These genes exhibited higher expression levels in ICI responders than in non-responders (as depicted in Figure 1A, B). Additionally, low expression of UTY, KDM5D, and RPS4Y1 was closely associated with a shorter PFS following PD-1/PD-L1 blockade (KDM5D, HR [95% CI]: 2.677 [1.572–10.250], log-rank P = 0.0044; RPS4Y1, HR [95% CI]: 2.495 [1.430–7.706], log-rank P = 0.0062; UTY, HR [95% CI]: 3.257 [1.619–6.238], log-rank P = 0.0010; Supplementary Figures 3D–F, https://links.lww.com/CM9/C240) in the discovery cohort. In the independent external validation cohort, patients with higher KDM5D and RPS4Y1 expression showed longer PFS (KDM5D, HR [95% CI]: 1.928 [1.222–3.289], log-rank P = 0.0067; RPS4Y1, HR [95% CI]: 2.061 [1.148–3.158], log-rank P = 0.0135), although we did not obtain results for the UTY gene because of insufficient samples [Supplementary Figure 3G, H, https://links.lww.com/CM9/C240]. Additionally, our study revealed that patients with a higher average expression of the aforementioned genes on the Y chromosome in PBMCs exhibited a higher percentage of whole-blood lymphocytes and a lower percentage of whole-blood neutrophils in both the discovery [Supplementary Figure 6A, https://links.lww.com/CM9/C240] and validation cohort [Supplementary Figure 6B, https://links.lww.com/CM9/C240]. Although the differences in the percentage of both lymphocytes and neutrophils based on different DDX3Y and RPS4Y1 expression levels were not statistically significant in the validation cohort [[Supplementary Figure 6B, https://links.lww.com/CM9/C240], their potential roles in immune cells should not be disregarded. Currently, the cause-and-effect relationship is not fully understood; however, our study proposes that the loss of chromosome Y (LOY) or extreme downregulation of chromosome Y (EDY) in immune cells may lead to a decrease in the proportion of lymphocytes, thereby resulting in primary resistance to PD-1/PD-L1 blockade immunotherapy. In summary, our findings suggest that LOY or EDY in PBMCs may serve as potential unfavorable predictors for PD-1/PD-L1 blockade response. This mechanism may function by inhibiting lymphocyte proliferation and activation, thereby weakening the immune response to ICI immunotherapy. In summary, our data showed that the expression of a cluster of Y-linked genes, including DDX3Y, USP9Y, KDM5D, UTY, and RPS4Y1, was significantly lower in non-responders to PD-1/PD-L1 blockade than in responders. The higher expression of these genes located on the Y chromosome was closely associated with better survival outcomes. Sex-related differences in immune responses are well acknowledged.[1] It was also observed that ICIs are more effective in males than females.[2,3] We propose that LOY and EDY may serve as potential biomarkers and play a role in the immune response to PD-1/PD-L1 blockade. LOY and EDY are common in patients with solid tumors, occurring at varying frequencies.[4] A recent study reported that LOY in cancer cells can drive immune escape and tumor growth.[5] LOY and EDY also occur in immune cells, and their interactions with tumor cells are crucial for determining tumor progression, patient prognosis, and the efficacy of immunotherapy. Changes in Y chromosome in immune cells play an important role in tumor immune responses and ICI efficacy. Our findings suggest that gene expression on Y chromosome may regulate lymphocyte proliferation; however, further research is needed to understand these mechanisms. In conclusion, DDX3Y and USP9Y emerge as useful biomarkers for predicting PD-1/PD-L1 blockade response. This study identified Y chromosome immunosurveillance as a key factor contributing to sex differences that affect the efficacy of PD-1/PD-L1 blockade. Funding This work was supported by the National Natural Science Foundation of China (No. NSFC82372331), the Project of Science and Technology Department of Sichuan Province (No. 2024NSFSC1551 and No. 2023NSFSC0716). Conflicts of interest None.
最长约 10秒,即可获得该文献文件

科研通智能强力驱动
Strongly Powered by AbleSci AI
更新
PDF的下载单位、IP信息已删除 (2025-6-4)

科研通是完全免费的文献互助平台,具备全网最快的应助速度,最高的求助完成率。 对每一个文献求助,科研通都将尽心尽力,给求助人一个满意的交代。
实时播报
罗一完成签到,获得积分10
2秒前
4秒前
丘比特应助wu采纳,获得10
7秒前
俏皮芷蕊发布了新的文献求助30
7秒前
称心的菲鹰完成签到,获得积分10
8秒前
碧蓝问安发布了新的文献求助10
9秒前
9秒前
打打应助ZZZ采纳,获得10
11秒前
15秒前
呆萌板凳发布了新的文献求助10
15秒前
hp关闭了hp文献求助
16秒前
17秒前
都选C完成签到,获得积分10
18秒前
壮观以松完成签到,获得积分10
18秒前
Liufgui应助郭小宝采纳,获得20
18秒前
heli完成签到,获得积分10
20秒前
如意枫叶发布了新的文献求助10
21秒前
都选C发布了新的文献求助10
22秒前
英俊的铭应助淡烟流水采纳,获得10
23秒前
23秒前
Miracle完成签到,获得积分10
25秒前
29秒前
wu发布了新的文献求助10
29秒前
忧心的听双完成签到,获得积分10
29秒前
Timon完成签到,获得积分10
30秒前
深情安青应助Miracle采纳,获得10
31秒前
李健应助猪猪hero采纳,获得10
31秒前
kingwill应助Harlotte采纳,获得60
31秒前
张雯思发布了新的文献求助10
32秒前
西瓜刀发布了新的文献求助10
32秒前
寒冷的发箍完成签到,获得积分10
34秒前
34秒前
34秒前
35秒前
wu完成签到,获得积分10
35秒前
orixero应助张大英采纳,获得10
35秒前
37秒前
heli发布了新的文献求助10
37秒前
38秒前
猪猪hero发布了新的文献求助10
38秒前
高分求助中
A new approach to the extrapolation of accelerated life test data 1000
ACSM’s Guidelines for Exercise Testing and Prescription, 12th edition 500
‘Unruly’ Children: Historical Fieldnotes and Learning Morality in a Taiwan Village (New Departures in Anthropology) 400
Indomethacinのヒトにおける経皮吸収 400
Phylogenetic study of the order Polydesmida (Myriapoda: Diplopoda) 370
基于可调谐半导体激光吸收光谱技术泄漏气体检测系统的研究 350
Robot-supported joining of reinforcement textiles with one-sided sewing heads 320
热门求助领域 (近24小时)
化学 材料科学 医学 生物 工程类 有机化学 生物化学 物理 内科学 纳米技术 计算机科学 化学工程 复合材料 遗传学 基因 物理化学 催化作用 冶金 细胞生物学 免疫学
热门帖子
关注 科研通微信公众号,转发送积分 3989242
求助须知:如何正确求助?哪些是违规求助? 3531393
关于积分的说明 11253753
捐赠科研通 3270010
什么是DOI,文献DOI怎么找? 1804868
邀请新用户注册赠送积分活动 882084
科研通“疑难数据库(出版商)”最低求助积分说明 809136