1. Introduction
Despite the development of new therapeutic and diagnostic strategies and even improved patient outcomes, stomach cancer is still a challenge, especially in advanced stages of cancer. In 2018, more than 1 million patients were diagnosed with stomach cancer worldwide, which is 5.7% of all cancer patients. Stomach cancer is aggressive, and 40% of patients have metastatic disease at diagnosis [
1]. Gastric cancer occurs in two main subtypes: intestinal type and diffuse type, which differ molecularly, epidemiologically, and clinically. The intestinal type is characterized by an organized cell structure resembling intestinal epithelium. It occurs mainly in older people and develops gradually, passing through stages of inflammation, atrophic mucositis, and intestinal metaplasia [
2,
3]. The diffuse type, on the other hand, occurs more often in younger people and has a more aggressive course. It is associated with mutations in the CDH1 gene, which leads to a loss of adhesion between cells. Cancer cells of this type are more diffuse and spread throughout the stomach wall, leading to their thickening. For this reason, this type of cancer is more difficult to detect in the early stages and has a worse prognosis [
2]. Therefore, research is constantly being conducted to detect this cancer earlier or to develop more effective treatment regimens. For this purpose, numerous studies are being conducted on the influence of the immune system on the development and treatment of stomach cancer. One of the elements of the immune system taking part in the fight against cancer is Natural killer (NK) cells, which, being an element of the innate response, are also able to influence the adaptive response, thus constituting a kind of bridge between these two types of immunity [
4]. Due to the influence of these cells not only on other cells through the secretion of cytokines, they can also actively destroy infected or degenerated tumor cells. At the same time, it is indicated that in the case of tumors, the number of NK cells infiltrating into the tumor sites is usually reduced, and the function and activation of NK cells are severely inhibited, which promotes the progression and metastasis of tumors [
5]. Reduced NK cell cytolytic activity and cytokine secretion capacity are significantly associated with poor prognosis in various cancers. In turn, increased NK cell activation and high NK cell lysis capacity correlate positively with more prolonged overall survival and a better prognosis for patients [
5].
Therefore, our study aims to evaluate the expression of Toll-like receptors (TLR-2, TLR-3, TLR-4, and TLR-9) on NK and NKT-like cells in patients with gastric cancer (GC), compare these results with healthy volunteers (HV), and investigate variants depending on the cancer subtype. Furthermore, we aimed to evaluate TLR gene expression using RT-qPCR to investigate their potential as diagnostic and prognostic biomarkers, with particular attention to the expression levels of these TLRs in patients with intestinal and diffuse gastric cancer subtypes, analyze their correlation with the stages of cancer progression, and compare the results in different patient groups and with the control group.
2. Materials and Methods
2.1. Research Material and Criteria for Patient Inclusion and Exclusion from the Study
The study included 86 patients with histologically confirmed diagnosis of gastric cancer (GC), regardless of the stage of the disease (I–IV), aged 18 years and older, with a performance status of 0–2 according to the Eastern Cooperative Oncology Group (ECOG) scale. All patients gave written, informed consent to participate in the study after familiarizing themselves with the purpose, scope, and potential risks of the study. Additionally, only peripheral blood in EDTA (10 mL) was collected from them for analysis. Exclusion criteria included the presence of active malignancy other than gastric cancer, except for fully cured basal cell carcinoma, squamous cell carcinoma of the skin, or cervical cancer in situ. Patients who had received chemotherapy, radiotherapy, immunotherapy, or other forms of anticancer treatment within the last 4 weeks before entering the study were also excluded, as were those with active autoimmune disease or chronic inflammatory disease that could affect the results of the study. Patients with severe heart, liver, or kidney failure (e.g., creatinine clearance < 30 mL/min, decompensated cirrhosis, NYHA class III-IV heart failure) were not included in the study, as were those with active, uncontrolled infection, including HIV, HBV, or HCV. Patients who were pregnant, breastfeeding, or unwilling to use effective contraception during the study were also excluded. Additionally, patients who did not provide informed consent to participate in the study were not included.
The control group included 30 healthy volunteers who met the inclusion criteria and did not meet the exclusion criteria. Inclusion criteria included individuals aged 18 years and older, in good health, with no symptoms of any chronic diseases or inflammatory conditions. Volunteers had to have a performance status of 0 according to the Eastern Cooperative Oncology Group (ECOG) scale. All participants gave written, informed consent to participate in the study after familiarizing themselves with its purpose, scope, and potential risks. Only peripheral blood on EDTA (10 mL) was collected from healthy volunteers for analysis. Exclusion criteria included the presence of any chronic disease, including autoimmune, neoplastic, cardiovascular, liver, and kidney diseases, and active infections, including HIV, HBV, or HCV. Volunteers taking any immunosuppressive, antineoplastic, or anti-inflammatory drugs chronically were also excluded from the study. Individuals who were pregnant, breastfeeding, or not using effective contraception were also excluded. Additionally, participants who did not provide informed consent to participate in the study were not included. Detailed information on patient characteristics is presented in
Figure 1 and
Supplementary Materials Table S1.
2.2. Immunophenotyping
To evaluate the expression of TLR2-9 on NK and NKT-cells, the peripheral blood samples from the participants were incubated with a set of monoclonal anti-human antibodies, including CD45 FITC, CD3 BV510, CD4 BV650, CD8 BV605, CD19 PE-Cy7, CD16 BB700, CD56 BUV661, TLR2 PE, TLR4 BV421, TLR3 PE, and TLR9 BV421 (BD, Franklin Lakes, NJ, USA). BD Horizon™ Brilliant Stain Buffer (BD, Franklin Lakes, NJ, USA) was used to enhance antibody stability, thus improving the signal quality for flow cytometry analysis. Erythrocytes were removed using a lysis buffer (BD, Franklin Lakes, NJ, USA), and cells were washed with BD Pharmingen™ Stain Buffer (BSA) (BD, Franklin Lakes, NJ, USA). For intracellular TLR assessment, the BD Cytofix/Cytoperm™ Fixation/Permeabilization Kit (BD, Franklin Lakes, NJ, USA) was applied. The assays were conducted on a CytoFLEX LX flow cytometer (Beckman Coulter, Indianapolis, IN, USA), and daily quality control was ensured with CytoFLEX Daily QC Fluorosphere reagents to minimize instrument-related variability. The acquisition of events was set to collect 15,000 events from the gate expressing CD45, which was the basis for further analyses (
Figure S1). Data were analyzed using Kaluza Analysis software version 2.1 (Beckman Coulter, Indianapolis, IN, USA).
2.3. PBMC Isolation
Peripheral blood was collected from healthy donors/patients into tubes with anticoagulant (EDTA). Blood was diluted 1:1 with physiological saline buffer (PBS) and gently applied onto a Gradisol L density gradient (Aqua-med). Samples were centrifuged at room temperature for 20 min at 400× g without inhibition. The peripheral blood mononuclear cell (PBMC) layer was collected from the interface and washed twice with PBS supplemented with 2% bovine serum (FBS) and 2 mM EDTA, centrifugation for 10 min at 300× g. PBMC was counted and suspended in MACS buffer (PBS, 2% FBS, 2 mM EDTA).
2.4. Separation of NK and NKT-Like Cells
According to the manufacturer’s instructions, NK cells were isolated from PBMC (107 cells/mL) using a commercially available NK Cell Isolation Kit, human (Miltenyi Biotec). The procedure used adverse selection, in which NK cells remained unbound to antibodies. PBMC were incubated with a mixture of antibodies binding markers specific for T cells (CD3), B cells (CD19), monocytes (CD14), granulocytes (CD15), and other nonspecific cells, which were then removed using magnetic beads conjugated with anti-mouse antibodies. The cell suspension was passed through a MACS column in a magnetic field. Purified NK cells were collected in the eluate by washing the column with MACS buffer.
NKT-like cells were isolated from PBMC (107 cells/mL) by a two-step method using a commercially available CD3+CD56+ NKT Cell Isolation Kit, human (Miltenyi Biotec). In the first step, NK cells and monocytes were magnetically labeled using a cocktail of biotin-conjugated antibodies and Anti-Biotin MicroBeads (Miltenyi Biotec). Labeled cells were removed from the sample by passing them through a MACS® column placed in a magnetic field, which allowed the separation of unbound cells as a fraction enriched in NKT cells.
In the second step, CD3⁺CD56⁺ cells were isolated from the obtained NKT cell-enriched fraction by positive selection using CD56 MicroBeads (Miltenyi Biotec) according to the manufacturer’s protocol. Labeled CD3⁺CD56⁺ cells were retained in the MACS® column, and unbound cells were washed away. Then, the column was removed from the magnetic field, and the isolated NKT-like cells were eluted, obtaining a fraction of high purity and viability, ready for further analysis. After separation, cell purity and viability were assessed by flow cytometry (using antibodies against CD3, CD56, and CD16 (BioLegend). NK cells were identified as CD3−CD56bright or CD56dim, and NKT-like cells as CD3+CD56+.
2.5. RNA Isolation Using A&A Biotechnology Kit
Total RNA isolation was performed using an RNA isolation kit (A&A Biotechnology, Gdańsk, Poland) according to the manufacturer’s protocol. Cells (or tissues) were lysed by adding an appropriate lysis buffer containing guanidine, which denatures proteins, including RNases, preventing RNA degradation. The mixture was mixed vigorously and incubated for several minutes at room temperature to lyse the cells completely.
After lysis, the suspension was transferred to an RNA-binding column, placed in a collection tube, and centrifuged for 1 min at 12,000× g, which allowed nucleic acids to bind to the column membrane. The column was then washed multiple times with special wash buffers, including a buffer containing ethanol, which allowed the removal of contaminants such as proteins, lipids, and DNA residues. After each wash, the column was centrifuged to remove excess fluid effectively. In the last step, the RNA was eluted from the column membrane using an elution buffer warmed to room temperature or RNase-treated water (DEPC-treated water). The column was placed in a clean elution tube and centrifuged for 1 min at 12,000× g. The resulting RNA was collected, and its concentration and quality were assessed by spectrophotometry (e.g., NanoDrop) and agarose gel electrophoresis or by microcapillary analyzer (e.g., Agilent Bioanalyzer) to ensure RNA integrity.
The isolated RNA was stored at −80 °C until further use in experiments such as RT-qPCR, RNA sequencing (RNA-seq), or gene expression analysis.
2.6. cDNA Synthesis
Isolated RNA was subjected to cDNA synthesis using the iScript™ cDNA Synthesis Kit (Bio-Rad, Hercules, CA, USA) according to the manufacturer’s instructions. The reaction was performed in a volume of 20 µL containing 1 µg RNA, 4 µL 5X iScript Reaction Mix, 1 µL iScript Reverse Transcriptase, and RNase-free water. The reaction was performed in a thermocycler (e.g., Simply Ampli, Applied Biotechnology) according to the following program: RNA denaturation and primerization for 5 min at 25 °C, cDNA synthesis for 20 min at 46 °C, and then reverse transcriptase inactivation for 1 min at 95 °C. The resulting cDNA was stored at −20 °C until further analysis.
2.7. RT-qPCR (Real-Time Quantitative PCR)
TLR2, TLR3, TLR4, and TLR9 gene expression analysis was performed by RT-qPCR using SYBR Green dye and SsoAdvanced Universal SYBR® Green Supermix® reagents (Bio-Rad). RT-qPCR reactions were performed in a volume of 20 µL containing: 10 µL SYBR Green Supermix, 1 µM PrimePCR Assay for TLR2 (Unique Assay ID: qHsaCED0036567; Assay Design: exonic; Chromosome Location: 4:154624891-154625028; Amplicon Length: 108; Amplicon Context Sequence: GGATTGTTAGAATTAGAGTTTGATGACTGTACCCTTAATGGAGTTGGTAATTTTAG AGCATCTGATAATGACAGAGTTATAGATCCAGGTAAAGTGGAAACGTTAACAATC CGGAGGCTGCATATTCCAAGGTTTTAC), TLR3 (Unique Assay ID: qHsaCED0046212; Assay Design: exonic; Chromosome Location: 4:187003545-187003664; Amplicon Length: 90; Amplicon Context Sequence: CAGCCTTACAGAGAAGCTATGTTTGGAATTAGCAAACACAAGCATTCGGAATCTG TCTCTGAGTAACAGCCAGCTGTCCACCACCAGCAATACAACTTTCTTGGGACTAA AGTGGACAAA), TLR4 (Unique Assay ID: qHsaCED0037607; Assay Design: exonic; Chromosome Location: 9:120474961-120475080; Amplicon Length: 90; Amplicon Context Sequence: CAAGATTCAAAGTATTTATTGCACAGACTTGCGGGTTCTACATCAAATGCCCCTA CTCAATCTCTCTTTAGACCTGTCCCTGAACCCTATGAACTTTATCCAACCAGGTG CATTTAAAGA), TLR9 (Unique Assay ID: qHsaCED0003672; Assay Design: exonic; Chromosome Location: 3:52259820-52259947; Amplicon Length: 98; Amplicon Context Sequence: CCTACATCCCATGAGGGCCTCACACCTGTCCTCTACCAAGCCCAGGGAGGAGCT AAGGCCCAGAGCTCAGGCAGAGAGCAGGGAGAGATGGGAATTCTGGATAGCAC CAGTAGCGGGTACACCTTGCT) and the GAPDH reference gene (Unique Assay ID: qHsaCED0038674; Assay Design: exonic; Chromosome Location: 12:6647267-6647413; Amplicon Length: 117; Amplicon Context Sequence: GTATGACAACGAATTTGGCTACAGCAACAGGGTGGTGGACCTCATGGCCCACAT GGCCTCCAAGGAGTAAGACCCCTGGACCACCAGCCCCAGCAAGAGCACAAGAG GAAGAGAGAGACCCTCACTGCTGGGGAGTCCCTGCCACAC), 2 µL cDNA (approx. 50 ng) and nuclease-free water to make up the volume. The reaction was performed on a CFX96 Real-Time PCR Detection System (Bio-Rad) with the following thermal profile: polymerase activation for 3 min at 95 °C, followed by 40 cycles consisting of denaturation at 95 °C for 10 s, primer annealing at 60 °C for 30 s, and extension at 72 °C for 30 s. The specificity of the amplification was verified by melting curve analysis in the temperature range from 65 °C to 95 °C, increased by 0.5 °C every 5 s. The results of mRNA expression analysis were normalized to the reference gene (GAPDH) using the 2−ΔΔCt method, and the data were presented as relative expression changes compared to the control sample.
2.8. Statistical Analysis
Data were statistically analyzed using Statistica (version 13.5.0.17, TIBCO Software Inc., Palo Alto, CA, USA) and GraphPad Prism (version 5.01 for Windows, GraphPad Software, San Diego, CA, USA). Due to the nonparametric nature of the data, the normality of their distribution was assessed using the Shapiro–Wilk test. In the absence of a normal distribution, appropriate nonparametric tests were used. Comparisons between the two groups were performed using the Mann–Whitney U test, which is a nonparametric test used to compare the medians of two independent groups. The test results were presented as the median and interquartile range (IQR). p values < 0.05 were considered statistically significant. Comparisons between multiple groups were performed using the Kruskal–Wallis test, which is used to assess differences between more than two independent groups. In the case of significant results of the Kruskal–Wallis test, additional post-hoc analysis with Bonferroni correction was performed to identify which groups were significantly different from each other. Bonferroni correction reduces type I error during multiple comparisons. Correlation analysis was performed using Spearman’s rank correlation coefficient, which is used to assess the strength and direction of the relationship between two variables with a nonparametric distribution. Correlation results were presented as Spearman’s correlation coefficient (rho) with the corresponding p value. p values < 0.05 were considered statistically significant. Receiver Operating Characteristic (ROC) curve analysis was performed to assess the ability of selected variables to differentiate between individual study groups. Predictor values for each variable were compared between the two groups, and ROC curves were plotted based on the calculated sensitivity and specificity at different cut-off values. Area Under Curve (AUC) was calculated as a measure of the predictive ability of the variable. AUC values were interpreted as follows: AUC = 0.5 indicated no discrimination, AUC > 0.7 was considered acceptable predictive ability, AUC > 0.8 as good, and AUC > 0.9 as very good discrimination between groups. The AUC value of each ROC curve was compared with the value of 0.5 using a significance test to assess whether the variable was significantly different from random prediction. p values < 0.05 were considered statistically significant.
4. Discussion
Studies conducted by our team have shown significant differences in the expression of TLR receptors between groups of patients with GC and in comparison to the HV group, which may have important diagnostic and prognostic implications. The statistical analysis performed showed a significantly higher expression of the tested TLR receptors (TLR-2, TLR-3, TLR-4, and TLR-9) in the group of patients with GC compared to HV, which suggests that these receptors play an essential role in the immune response to GC.
The signaling pathways involved in NK cell activation by TLRs differ depending on their location—TLR2 and TLR4 are usually present on the cell membrane surface, whereas TLR3 and TLR9 are located intracellularly and are activated upon entry of ligands into the cell. TLR2 recognizes lipoproteins and peptidoglycans derived from Gram-positive bacteria, as well as microbial and fungal components. TLR2 can form heterodimers with TLR1 or TLR6, which broadens the range of ligands it recognizes. Upon activation by a ligand, TLR2 recruits the adaptor MyD88 (myeloid differentiation primary response 88), which leads to the activation of the NF-κB pathway and, consequently, to the transcription of proinflammatory genes and the production of cytokines [
6,
7,
8]. TLR4 is specific for lipopolysaccharides (LPS) present on the surface of Gram-negative bacteria. Its activation initiates a signaling cascade through both MyD88 and TRIF (Toll/IL-1 receptor domain-containing adaptor-inducing interferon-β), leading to the activation of two major signaling pathways: NF-κB and IRF3 (interferon regulatory factor 3). TLR4, in addition to NF-κB, stimulates the production of type I interferons, which enhances antiviral and antitumor responses [
9,
10,
11]. TLR3 is located primarily in endosomes and recognizes double-stranded RNA (dsRNA), which is typical for many viruses. After binding dsRNA, TLR3 mainly activates the TRIF pathway, which leads to the production of proinflammatory cytokines, especially type I interferons. In NK cells, TLR3 activation increases their cytotoxicity and IFN-γ production, supporting antiviral and antitumor responses [
12]. TLR9 is also intracellular and localized in endosomes, where it recognizes unmethylated CpG sequences found in the DNA of bacteria and some viruses. After activation, TLR9 uses the MyD88 adaptors, which leads to the activation of NF-κB and IRF7, resulting in the production of proinflammatory cytokines and interferons. In NK cells, TLR9 activation promotes their cytotoxic activity and cytokine production. Differences in TLR localization affect the way they detect pathogens and activate NK cells [
13,
14]. TLR2 and TLR4, as membrane receptors, enable a rapid response to external threats, while TLR3 and TLR9, as intracellular receptors, detect pathogens after they enter the cell, which is particularly important for antiviral responses and responses against intracellular pathogens.
One of the study’s key findings is the significantly higher expression of TLR in CD3-CD56+ and CD3+CD56+ cells in patients with GC compared to HV. The mean expression of TLR-2 in NK cells of patients with GC was 13.37±10.41, while in the HV group, it was only 1.24 ± 0.84 (
p = 0.000). A similar pattern was observed for other receptors studied, where TLR-3, TLR-4, and TLR-9 were also significantly higher in the GC group. For example, the expression of TLR-9 in CD3-CD56+ cells was, on average, 12.42 ± 9.53 in the GC group and 1.26 ± 0.87 in the HV group (
p = 0.000), which indicates apparent differences in the immune response between patients and healthy individuals. These results are consistent with the literature suggesting that TLRs may be vital in recognizing cancer cells and initiating the inflammatory response [
15,
16,
17,
18,
19]. At the same time, the observed higher expression of TLR receptors in GC patients suggests that these receptors may be associated with disease progression. Studies show that TLR-2, TLR-3, TLR-4, and TLR-9 are involved in the immune response to tumors by stimulating CD3-CD56+ and CD3+CD56+cells, which may lead to increased cytotoxicity and production [
20,
21,
22]. In the case of cancers such as GC, TLR expression may play a dual role: on the one hand, it may promote an anti-tumor response, and on the other hand, it may be associated with a chronic inflammatory state that supports tumor development [
20,
21,
22]. Moreover, differences in TLR expression on CD3-CD56+ and CD3+CD56+cells may be related to adaptive changes in the immune response as the disease progresses. Data suggest that higher levels of TLR expression in mature NK cells are observed in more advanced stages of cancer, which may indicate an adaptive immune response of the organism to disease progression. In advanced stages of cancer, NK cells may increase TLR expression as an adaptive mechanism, enabling a more effective response to danger signals from the tumor microenvironment. It has been suggested that increased TLR expression may support NK cell activation and their ability to produce cytokines and destroy tumor cells, which is particularly important in the face of an increased number of tumor cells and the presence of pathogens or DNA fragments that stimulate TLRs. High TLR expression in mature NK cells, especially in those with the CD56dim phenotype, may indicate an increased readiness to respond to danger signals, which is crucial in antitumor and antiviral responses. NK cells with the CD56 dim phenotype are known for their high cytotoxicity—they are particularly effective in directly destroying cancer or virus-infected cells by releasing perforins and granzymes, which lead to apoptosis of target cells. As GC progresses, increased activation of TLRs may reflect an enhanced immune response, which in turn may impact tumor progression and patient prognosis. In studies of other cancers, such as skin cancer, TLR agonists have been successfully used to stimulate the immune response, suggesting that similar approaches may apply to the treatment of GC [
14,
18,
23,
24,
25]. The differences in TLR expression between CD3-CD56+ and CD3+CD56+cells observed in our study may also provide additional information on the immunological mechanisms involved in gastric cancer. TLR expression was higher in patients with GC in both cell subpopulations, but in CD3-CD56+ cells, more significant differences were observed compared to healthy volunteers. This may indicate a more significant role of CD3-CD56+ cells in the antitumor response in gastric cancer. Moreover, increased TLR expression in CD3+CD56+ cells may suggest that these cells also play an essential role in the immune response to tumors. However, their role may be more related to the modulation of the inflammatory response than direct cytotoxicity [
19,
26,
27].
It is also worth mentioning that the highest expression of TLRs was observed in patients in the advanced stages of the disease. This may indicate a more intense immune response activation as the disease progresses. In particular, the expression values of TLR-2, TLR-3, TLR-4, and TLR-9 on CD3-CD56+ cells increased proportionally with the tumor stage, suggesting that these receptors may be functional as diagnostic and prognostic biomarkers in assessing the progression of GC. Such a phenomenon may reflect the increasing demand for an immune response in more advanced stages of the disease [
28,
29]. Our studies suggest TLR receptors may become attractive therapeutic targets for treating GC. There are reports in the literature on the use of TLR agonists in treating other types of cancers, which indicates the possibility of using similar strategies in GC [
14,
26,
27,
28,
30]. TLR9 agonists are currently being investigated in anticancer therapy, both as monotherapy and in combination with other treatments such as chemotherapy and immunotherapy. Moreover, increased expression of TLRs in CD3-CD56+ and CD3+CD56+ cells in GC patients suggests that stimulation of these receptors could contribute to an enhanced immune response to the tumor [
14,
26,
27,
28,
30].
Limitations of the Study
Although the results of the study presented in this publication are statistically significant, they only show a small fragment of the immunopathogenesis of GC. It should be noted that the study was conducted on a relatively small number of patients with GC and HV, which may limit the possibility of generalizing the results to a broader population. The study did not include a detailed analysis of potential confounding factors, such as lifestyle, diet, or exposure to environmental factors that may affect TLR expression in CD3-CD56+ and CD3+CD56+. Therefore, this aspect should also be analyzed in future studies. Moreover, the study participants differed not only in terms of the type of cancer (diffuse vs. intestinal), but also in the advancement of the disease. The insufficient number of these groups may limit the possibility of precise assessment of the relationship between TLR expression and specific clinical features of the patients. Additionally, the presented analyses were cross-sectional and did not include long-term monitoring of patients with GC, which makes it impossible to assess changes in TLR expression over time and their potential impact on disease progression, response to treatment, or prognosis. The obtained results of the study show us an insight into the deregulation of the immune system by CD3-CD56+ and CD3+CD56+ cells at the time of GC diagnosis and may be a good starting point for further detailed analyses. All these limitations indicate the need for further, more detailed studies with larger cohorts, taking into account the control of confounding factors and functional analysis of TLR receptors to better understand their role in the pathogenesis and potential treatment of GC.