$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
The research involving human subjects was reviewed and approved by the Ethics Committee of the First Affiliated Hospital of Soochow University (Approval No. 2019026). The ethical review board of the First Affiliated Hospital of Soochow University has approved or granted an exemption for the ethical review of this data usage. Additionally, informed consent was obtained from all participants.
Bioinformatics data analysis
After a thorough Pan-cancer expression analysis examination of the TCGA and Tumor IMmune Estimation Resource (TIMER) databases, the behavior of LLT1 was investigated across various cancer types, with a particular focus on gastric cancer. The TISIDB database was utilized to investigate the correlation between LLT1 expression levels and different molecular subtypes of gastric cancer. The correlation between patients' clinical characteristics and LLT1 expression was further analyzed with TCGA data through R programming. GO and KEGG analyses were employed to investigate the functional role of LLT1 in gastric carcinogenesis. The ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues using Expression data) algorithm was additionally employed to quantify immune and stromal cell infiltration levels in tumor tissues according to LLT1 expression patterns. Further analysis utilizing the CIBERSORT algorithm was conducted to evaluate correlations between LLT1 expression levels and specific immune cell infiltration patterns within the tumor microenvironment. Comprehensive bioinformatic analysis was performed using TCGA datasets and R statistical programming (v4.4.1) to interrogate potential correlations between LLT1 expression profiles and immunoregulatory biomarkers. Finally, the TCIA (The Cancer Immunome Atlas) database was utilized to investigate whether LLT1 expression could predict response to PD-1/PD-L1 immune checkpoint blockade therapy.
GAC specimens and tissue microarray (TMA) construction
Two sets of tissue microarray data were generated as follows:
The first set retrospectively analyzed specimens from 268 patients with GAC who underwent surgery in 2011. Pathologists selected paraffin blocks, identified tumor core regions, and constructed them into tissue microarrays. Every individual sample site was allocated a distinct identifier to meticulously monitor associated clinical parameters such as sex, age, neoplasm dimensions, degree of neoplastic differentiation, phase of cancer progression, invasion depth, lymphatic involvement, and the presence of distant spread. The duration of overall survival (OS) was defined as the interval from the initiation of therapy until the occurrence of mortality or the cutoff date of December 2017. It should be highlighted that comprehensive OS data was only available for a subset of 198 participants.
The second set included 115 treatment-naïve patients with advanced GC, who underwent surgical/endoscopic resection at the First Affiliated Hospital of Soochow University between 2016 and 2023. Pathologists evaluated the paraffin blocks and constructed tissue microarrays with close follow-up of disease progression. Notably, all specimens were collected before any preoperative treatment, and postoperative treatment primarily consisted of PD-1 monoclonal antibodies combined with chemotherapy.
All the surgical resection specimens were collected from the antrum of the stomach. All samples were fixed in neutral buffered formalin and embedded in paraffin. Based on the evaluation of H&E-stained sections, two senior pathologists independently identified the most representative tissue regions from the paraffin-embedded specimens. The tissue microarrays (TMAs) were constructed using an automated tissue arrayer, with tissue cores precisely positioned in the TMA blocks according to predefined modules. Finally, the completed TMA blocks were sectioned at 4 µm thickness by an experienced histotechnologist.
Immunohistochemistry (IHC) and multiplex immunofluorescence (mIHC)
Immunohistochemical and multiplex immunofluorescence assays were conducted on a constructed tissue microarray. IHC staining was carried out in accordance with a previously established methodology26. The first set of microarray paraffin sections was transferred onto glass slides and baked at 72°C for 1 h. After the sections were deparaffinized with xylene and rehydrated in successively increasing ethanol dilutions, endogenous peroxidase activity was blocked by the application of peroxidase blocking solution for 30 min. Antigen retrieval was performed by boiling the slides at 110 °C in Citrate Antigen Retrieval Solution (pH 6.0) for 3 min. After blocking with 5% BSA blocking buffer, the sections were incubated with mouse anti-human CLEC2D monoclonal, mouse anti-human CD3 mono clonal (ready to use), mouse anti-human CD4 monoclonal (ready to use), rabbit anti-human CD8 monoclonal (ready to use), mouse anti-human CD11c monoclonal (1:50) mouse anti-human CD19 monoclonal (ready to use), mouse anti-human CD31 monoclonal (ready to use), mouse anti-human CD56 monoclonal (ready to use), mouse anti-human CD68 monoclonal (ready to use), rabbit anti-human Foxp3 monoclonal (1:50), mouse anti-human IL-17 monoclonal (1:50), mouse anti-human IFN-γ monoclonal (1:50), mouse anti-human Ki67 monoclonal (1:50), mouse anti-human PD-L1 monoclonal (ready to use) at 4 °C overnight27,28. The slides were then washed and incubated with horseradish peroxidase-labeled secondary antibodies (anti-mouse and rabbit) at room temperature for 30 min. After incubation, the samples were washed, and the EnVision detection system was used for immunodetection. Under the microscope, when brown areas begin to appear on the chip, careful observation is required. The chromogenic reaction should be stopped once the brown areas no longer increase. Slides were counterstained with Mayer's hematoxylin and mounted with neutral glue. In the negative control, mouse or rabbit IgG was used in place of the primary antibody. The positive control was human tonsil tissue. The second set of microarray sections is used for multiplex immunohistochemistry (mIHC). Reusing the aforementioned antigen retrieval and blocking conditions, the antibody panel is spectrally optimized as follows: PD-L1 (480 nm), CD161 (520 nm), LLT1 (570 nm), PD-1 (620 nm), CD8 (690 nm), and PAN-CK (780 nm). Primary antibodies were incubated under the same conditions as in the original protocol29 (except for PANCK and CD161, which were newly introduced), at the recommended 1:100 dilution, followed by tyramide signal amplification (TSA) fluorescence staining. After each round of fluorescence staining, antibody stripping is performed, followed by repeated antigen retrieval, endogenous peroxidase blocking, and non-specific target blocking before incubating the next primary antibody. When performing multiplex immunohistochemistry, three points are crucial: pre-experiments must be conducted for each fluorescent dye, strongly expressed antibodies should be paired with weak fluorescent signals, while weakly expressed antibodies should be paired with strong fluorescent signals, and after staining each fluorescent marker, the staining results must be observed under a microscope under light-protected conditions. After staining with the 6th fluorescent dye, DAPI counterstaining is applied to label cell nuclei, and the slides are mounted. Finally, all tissue microarrays were scanned using a digital pathology scanner and analyzed with an AI-powered digital pathology image analysis software.
IHC/mIHC staining evaluation
The staining outcomes were interpreted by two expert pathologists in a double-blind fashion. They gauged the depth of staining using a high-power objective (200x magnification) or matching fluorescent setting (200x), assigning scores ranging from 0 (no staining) to 3 (intense staining). The fractions of LLT1-positive cells within gastric carcinoma samples were divided into four stratifications: 0 (<5%), 1 (5-25%), 2 (26-50%), and 3 (>50%). The Immunoreactivity Score (IRS) was derived from the product of the staining depth and proportion of LLT1-positive cells. Specimens scoring zero were regarded as negative, whereas scores of 1-4 indicated weak expression, 5-7 moderate, and 8-9 strong expression; the same criteria were applied to evaluate the expression of other markers such as Ki67 and CD3. Based on the evaluation and grouping methodology, the samples were labeled as having strong (IRS8-9), moderate (IRS5-7), weak (IRS1-4), or no expression (IRS0). Patients with the lowest scores (IRS≤4) were assigned to the low-expression group, and those with higher scores (IRS≥5) were assigned to the high-expression group.
However, for the PD-L1 marker, a different approach was adopted using the CPS. This score was determined by the ratio of PD-L1-positive cells (including tumor, immune, and stromal cells) to the total number of viable tumor cells scaled up by a factor of 100. Individuals with a CPS exceeding five were placed in the high-expression cohort, whereas those with a CPS of five or below were grouped into the low-expression cohort.
Immunotherapy response grouping
In line with the Response Evaluation Criteria in Solid Tumors version 1.1 (RECIST 1.1), the evaluation of target lesions in patients falls into four distinct categories: (a) Complete Response (CR): all target lesions have vanished entirely, and the short axis of all pathological lymph nodes (both target and non-target) must measure less than 10 mm; (b) Partial Response (PR): the cumulative diameter of all measurable target lesions has reduced by at least 30% from the baseline measurements; (c) Progressive Disease (PD): this category is applicable when the cumulative diameter of target lesions has increased by at least 20% compared to the smallest recorded value during the study (using the baseline as a reference if it represents the smallest value), with the absolute increase being no less than 5 mm; additionally, the emergence of new lesions is also categorized as Progressive Disease (PD); (d) Stable Disease (SD): the cumulative diameter of the target lesions neither achieves the reduction criteria for PR nor the increase criteria for PD throughout the study. Two stratification techniques were implemented. The first method involved segregating patients based on their response status, where PD (n=53) and SD (n=46) patients comprised the non-responder group (n=99), whereas PR (n=15) and CR (n=1) patients formed the responder group (n=16). The second approach entailed categorizing patients according to disease progression, with CR, PR, and SD patients classified as the immunotherapy-effective cohort (n=62) and PD patients as the immunotherapy non-effective cohort (n=53).
Statistical analysis
Empirical evaluations were performed using the R computing platform (version 4.4.1) and the SPSS Statistical software package (version 27.0). Disparities among study cohorts were investigated by the application of the chi-square contingency table analysis. Spearman's non-parametric correlation technique was employed to probe the strength and direction of the associations between variables. Survival time variations were analyzed through the construction of Kaplan-Meier estimators complemented by log-rank testing for significance. Subsequently, Cox proportional hazards regression models were employed to identify factors potentially influencing patient survival outcomes. Statistical significance was predefined as a two-sided p-value threshold of <0.05. To ensure the robustness of the findings, comprehensive sensitivity analyses were conducted to evaluate potential confounding variables, with multivariate models adjusted accordingly. Furthermore, explicit interpretations of hazard ratios were provided to facilitate clinical understanding, emphasizing the practical implications of identified risk factors on patient prognosis.