The pipeline begins with one of the most detailed single-cell RNA sequencing atlases of triple-negative breast cancer available, dataset GSE176078, which profiles more than 100,000 cells from 26 treatment-naive patients. After rigorous quality control that removed cells with excessive mitochondrial content or too few detected genes, the researchers normalized expression data, identified 2,000 highly variable genes, and applied principal component analysis, k-nearest neighbor graph construction, and Leiden clustering to organize the cellular landscape. The resulting atlas resolved nine major cell compartments, including T cells, myeloid cells, B cells, cancer epithelial cells, and cancer-associated fibroblasts. Validation against the curated reference annotations produced adjusted rand index values between 0.288 and 0.311 and normalized mutual information values between 0.616 and 0.671, confirming that the unsupervised clustering substantially recovered known biology despite the inherent difficulty of matching 35 data-driven clusters to nine curated labels.

With the atlas established, the analysis zeroed in on the T cell compartment, extracting roughly 15,000 to 30,000 T cells and re-processing them independently to resolve functional sub-states. Two gene-set module scores quantified exhaustion and cytotoxicity across the population. The centerpiece of the study, however, is a deceptively simple quantity: the per-cell ratio of PDCD1, the gene encoding the inhibitory receptor PD-1, to CD2, the adhesion receptor that stabilizes the immune synapse between T cells and their targets. When PD-1 expression dominates and CD2 is low, inhibitory signaling suppresses activation so thoroughly that even dense CD58 ligand on an opposing surface cannot overcome it. When CD2 is moderately expressed but PD-1 only mildly elevated, calibrated CD2/CD58 ligand density can restore synapse stability and partial T cell activation. Prior work had shown that CD2 expression on tumor-infiltrating CD8 T cells correlates negatively with exhaustion markers and that CD2/CD58 co-stimulation is less sensitive to PD-1 inhibition than canonical CD28 signaling.

Aggregating these measurements per patient revealed three distinct immune phenotypes across the cohort. Roughly eight patients formed a high-exhaustion group in which more than 65 percent of CD8 T cells were classified as exhausted and PD-1 expression exceeded CD2 by three-and-a-half to five-fold. Their cells showed the lowest cytotoxic scores, consistent with a terminally exhausted state in which effector function is epigenetically silenced, making CD2-directed surface engineering alone an insufficient strategy. A second group of about ten patients retained measurable cytotoxic capacity alongside moderate exhaustion; their principal barrier appeared to be inadequate synapse stabilization rather than irreversible transcriptional shutdown, making them the most tractable candidates for calibrated CD2/CD58 ligand density optimization. A third group of roughly eight patients showed sparse CD8 infiltration altogether, representing immune-cold tumors where T cell recruitment must precede any activation engineering. Across all patients, mean PDCD1 expression spanned a seven-fold range and the resulting PDCD1/CD2 ratios ranged from approximately 0.3 to 4.1, capturing clinically meaningful variability invisible to bulk gene expression summaries.

To test whether the ratio carries clinical weight beyond the single-cell cohort, the team turned to TCGA-BRCA bulk RNA sequencing data of roughly 1,100 samples. In unadjusted Cox regression, higher bulk PDCD1/CD2 ratios were associated with lower mortality, with a hazard ratio of 0.47 and a confidence interval of 0.28 to 0.79, reaching statistical significance below 0.005. That direction appears paradoxical at first, since a high ratio reflects exhaustion at the single-cell level. The authors are careful to explain the discrepancy: in bulk tissue, PDCD1 transcripts derive from the entire immune infiltrate rather than from exhausted CD8 cells alone, so a high bulk ratio primarily signals strong immune infiltration, itself a favorable prognostic factor in breast cancer. This cell-type composition effect inverts the ratio’s meaning across measurement platforms, and the study explicitly frames the survival association as exploratory and hypothesis-generating rather than as validation of a clinical biomarker.

Cross-modal comparison of the 24 patients with matched bulk and single-cell data reinforced this interpretation with appropriate caution. Spearman correlations between single-cell-derived and bulk-derived ratios were weakly negative at approximately minus 0.30, consistent with the infiltration confound compressing the denominator, while exhaustion scores showed a weakly positive correlation of about 0.28, directionally consistent with the hypothesis that higher single-cell exhaustion burden corresponds to higher bulk exhaustion signal. Neither correlation reached conventional statistical significance given the small matched sample, but the directional agreement across platforms suggests the single-cell signals are not artifacts of normalization. Kaplan-Meier analysis stratified by the bulk ratio median showed longer median overall survival in the high-ratio group, with median survival differences of roughly 12 to 18 months and curve separation emerging after about two years of follow-up.

Beyond stratification, the pipeline delivers what previous tools have not: engineering output. A targeted ligand-receptor proxy screen across five immune axis pairs placed PD-1/PD-L1 highest in both T cell-tumor and T cell-myeloid pairings, reflecting the co-elevation of PD-1 on exhausted T cells and PD-L1 across tumor and myeloid compartments. CD2/CD58 ranked moderately at the tumor interface but low with myeloid cells, indicating that CD2-directed ligand optimization is specific to the direct tumor-T cell boundary. LAG-3/HLA-DRA showed the inverse profile, strongest with myeloid and B cell ligand sources, while CD28 co-stimulation was broadly downregulated throughout the tumor microenvironment, consistent with exhaustion-associated loss of CD28 expression. The authors emphasize these rankings are comparative heuristics from product-of-means calculations, not statistically validated communication events, and require experimental confirmation through co-culture or blocking assays.

The final translation step is the DesignPriorityScore, a rule-based metric combining normalized exhaustion burden with 50 percent weight, PDCD1/CD2 axis imbalance with 30 percent weight, and CD8 infiltration with 20 percent weight. Each of the 26 patients receives a ranked score and a corresponding recommendation: recruitment-first strategies for immune-desert tumors, combined PD-1 blockade with CD2 reinforcement for the highest-exhaustion group, CD2/CD58 axis optimization for the moderate group, and CD28 co-stimulation otherwise. Bootstrap resampling with 200 iterations, threshold sensitivity testing across exhaustion quantiles from 0.60 to 0.90, and weight sensitivity analysis confirmed ranking stability, with mean Spearman correlations above 0.85 and top-quartile patients retaining their position more than 90 percent of the time. Adjacent patients were separated by an average score gap of 0.03 normalized units, sufficient to avoid ambiguous recommendations.

The authors are candid about limitations. The pipeline was built and tested on a single cohort of 26 patients, exhaustion states were assigned by quantile-based gene-set scoring rather than experimental annotation, and the survival associations remain unadjusted for tumor purity, stage, and molecular subtype. The scoring weights are biologically motivated rather than data-trained, and no functional experiments have yet linked the recommendations to actual T cell activation outcomes. Future directions include integrating the spatial transcriptomics data that accompany the atlas to map where exhausted cells reside relative to tumor cells, training the score on labeled immunotherapy response cohorts, and connecting patient-level ratios to biophysical models of ligand density and clustering on synthetic cell membranes. As a prototype, the study demonstrates that single-cell immune phenotyping can be pushed beyond description into prescriptive design, offering synthetic immunologists a patient-specific starting point for building the engineered immune cells of tomorrow.

Subject of Research: A computational framework translating single-cell PDCD1/CD2 immune axis measurements into patient-specific synthetic immune cell engineering priorities for triple-negative breast cancer.

Article Title: Computational stratification and engineering framework of the PDCD1/CD2 immune axis in triple-negative breast cancer

Article References: Chowdhury, K. (2026). Computational stratification and engineering framework of the PDCD1/CD2 immune axis in triple-negative breast cancer. Clinical Cancer Bulletin, 5(1), Article 11. https://doi.org/10.1007/s44272-026-00063-5

Image Credits: AI Generated

DOI: 10.1007/s44272-026-00063-5

Keywords: triple-negative breast cancer, single-cell RNA sequencing, T cell exhaustion, PDCD1/CD2 axis, immune checkpoint blockade, synthetic immunology, tumor microenvironment, immune stratification, CD2/CD58 co-stimulation, DesignPriorityScore, TCGA-BRCA, ligand-receptor interactions