HLA-G 3’UTR Polymorphisms Are Linked to Susceptibility and Survival in Spanish Gastric Adenocarcinoma Patients

HLA-G is a non-classical class I HLA molecule that induces tolerance by acting on receptors of both innate and adaptive immune cells. When overexpressed in tumors, limits surveillance by the immune system. The HLA-G gene shows several polymorphisms involved in mRNA and protein levels. We decided to study the implication of two polymorphisms (rs371194629; 14bp INS/DEL and rs1063320; +3142 C/G) in paired tissue samples (tumoral and non-tumoral) from 107 Spanish patients with gastric adenocarcinoma and 58 healthy control individuals, to assess the possible association of the HLA-G gene with gastric adenocarcinoma susceptibility, disease progression and survival. The presence of somatic mutations involving these polymorphisms was also analyzed. The frequency of the 14bp DEL allele was increased in patients (70.0%) compared to controls (57.0%, p=0.025). In addition, the haplotype formed by the combination of the 14bp DEL/+3142 C variants is also increased in patients (54.1% vs 44.4%, p=0.034, OR=1.74 CI95% 1.05-2.89). Kaplan-Meier analysis revealed that 14bp DEL/DEL patients showed lower 5-year life-expectancy than INS/DEL or INS/INS (p=0.041). Adjusting for TNM staging (Cox regression analysis) disclosed a significant difference in death risk (p=0.03) with an expected hazard 2.6 times higher. Finally, no somatic mutations were found when comparing these polymorphisms in tumoral vs non-tumoral tissues, which indicates that this is a preexisting condition in patients and not a de novo, tumor-restricted, event. In conclusion, the variants predominant in patients were those increasing HLA-G mRNA stability and HLA-G expression, clearly involving this molecule in gastric adenocarcinoma susceptibility, disease progression and survival and making it a potential target for immunotherapeutic approaches.

Introduction

Gastric epithelial adenocarcinomas are the most common form of stomach tumors, representing 90% of the cases. Adenocarcinomas are usually located in the cardia (31.0%), antrum (26.0%) and body (14.0%), and their development is associated to the infection by Helicobacter pylori, as 84.0% of the patients have been infected with this bacterium (1). In 2020, this tumor was the 6th most frequent cancer in the world (11.1 new cases per 100.000 persons per year), and the 5th with highest mortality rate (7.7 patients per 100.000 persons per year), which made it one of the most aggressive tumors. The 5-year survival rate for this type of cancer was 29.5%, a situation that has been maintained for the last 30 years (2, 3).

Tumors with high frequency of somatic mutations, such as stomach malignancies (4), express tumor neo-antigens on their surface and are, thus, potential targets of the immune response. However, even in an immunocompetent organism, neoplastic cells develop mechanisms, such as the expression of immunomodulatory molecules as HLA-G, to evade the action of the immune system. Published works describe that HLA-G possesses a key role in the cancer immunoediting mechanism, attenuating the elimination of tumor cells (5–8).

HLA-G is a non-classical class I HLA molecule composed of a heavy chain bound to β2 microglobulin. The HLA-G gene, located on chromosome 6, exhibits 7 introns and 8 exons that codifies for the heavy chain. Exon 1 codes for the signal peptide, exons 2, 3 and 4 the extracellular domains α1, α2 and α3, respectively, and exons 5 and 6 the transmembrane region and the cytoplasmic domain, respectively (9). Exon 7 is transcribed in the pre-mRNA molecule, but not present in the mature-mRNA, whereas exon 8 is not translated., but in the latter lays the 3’UTR region, involved in the transcriptional regulation of the gene (10, 11).

HLA-G shows tolerogenic functions by inhibiting immunocompetent cells, through the interaction of HLA-G (whether membrane or soluble isoforms) with cognate receptors. These receptors are the immunoglobulin-like transcription receptor type 2 (ILT2, CD85j), present on NK, B, T cells and antigen-presenting cells (APCs) (12, 13); the ILT4 (CD85d) receptor (that along with ILT2, displays ITIM motifs), unique to APCs (12, 14); the killer immunoglobulin-like receptor (KIR) 2DL4, present on NK cells (15); and, finally, CD8 (16). Because of this tolerogenic role, HLA-G is involved in a wide variety of processes, such as maternal-fetal tolerance (17), organ transplantation (17), viral infections (18, 19), autoimmunity (18, 20), and cancer progression (5).

The 3’UTR region of the HLA-G gene shows numerous variations that can have an impact on the mRNA and, therefore, on the protein levels (10). Of these variants, rs371194629 and rs1063320 have been studied in several types of cancer and pathologies (21–23).

The rs371194629 polymorphism (14bp INS/DEL) is caused by the deletion (DEL) from the ancestral variant (INS) (10, 24) of a 14bp segment (5′-ATTTGTTCATGCCT-3’) located at the +2960 position of the 3’UTR region. This 14bp segment has been associated with both, the splicing and the stability of the mRNA (10, 25, 26), as it contains an AUUUG domain putatively exerting an AU-pentamer-like effect, decreasing mRNA stability (27). Therefore, the DEL allele provides a higher stability of the mRNA (25), associated with a high expression of HLA-G (26).

The rs1063320 polymorphism (+3142C/G), consists of the transversion of a cytosine (C, the ancestral variant) to guanine (G) at position +3142 of the 3’UTR region, modulating the affinity of 148a, 148b and 152 miRNAs (known to either favor the direct mRNA degradation or block mRNA to protein translation) for this region (28–30). Should a C be found at position +3142, miRNAs affinity will decrease, increasing the mRNA availability and the production of HLA-G (10, 11).

Other authors have described the expression of HLA-G in gastric tumors (31), therefore, studying the polymorphisms related to the expression levels (such as 14bp INS/DEL and +3142 C/G) of this molecule could serve to identify new genetic markers involved in the risk and evolution of this pathology.

Therefore, we decided to study the influence of the aforementioned polymorphisms and the combined haplotypes thereof in paired tissue samples (tumoral and non-tumoral) from a group of patients with gastric adenocarcinoma, and compare the frequencies achieved with that of a control population. This approach will allow us to determine whether the HLA-G variants may be adequate gastric adenocarcinoma risk markers and whether somatic mutations take place in tumoral, but not healthy, gastric tissue.

Patients, Materials, and Methods

Samples used in this study were obtained from the Servicio de Cirugía General y Digestiva, Hospital Príncipe de Asturias (Alcalá de Henares - Madrid) and sent to the Departamento de Inmunología (Facultad de Medicina. Universidad Complutense de Madrid - Madrid), where they were processed.

Patients

One hundred and seven Spanish patients diagnosed with gastric adenocarcinoma were included in this study. Patients were classified according to the TNM staging criteria (stages I through IV) (32) (Table 1).

TABLE 1

Table 1 Demographic and clinical characteristics of patients.

The inclusion criteria for this study were patients (women or men) over 18 years old with gastric adenocarcinoma, stratified according to UICC/AJCC criteria (7th edition 2009) (32), in any TNM stage and Karnofsky index >70% or performance status ≤ 2. All patients signed an informed consent before their inclusion in the study, and all samples were anonymized upon arrival to the laboratory.

The exclusion criteria included resecability of the primary tumor, coexistence with other neoplastic diseases, pregnancy and severe alterations of hepatic, cardiovascular o renal function.

Comorbidities included in this study were ULCUS (21%), hypertension (33%), cardiopathy (8%) pulmonary disease (21%), liver disease (2%), renal disease (2%), cerebrovascular disease (4%), diabetes (17%), smoker (19%) and alcoholism (4%).

Tissue Samples

Tumoral (T) and distal, non-tumoral (NT), gastric tissue samples were obtained from each patient upon surgery. A total of 214 tissue samples (T+NT) were available for the study.

Controls

A total of 58 sex-and age-matched Spanish healthy donors, from the same geographic region as the patients, were included as controls. DNA was obtained from blood or saliva samples, as described in the following section.

Genomic DNA Extraction

DNA isolation, both from blood and tissue, was carried out by using the Illustra Nucleon BACC (GE Healthcare) kit, following the manufacturer’s instructions. In the case of tissue specimens, fragments of 25 mg were mechanically disrupted, subjected to proteinase K treatment, and followed by the DNA precipitation protocol included in the kit. DNA from saliva control samples were extracted using the Oragene DNA 500 kit (DNA Genotek) and purified with PrepIT-L2P (DNA Genotek).

The concentration and quality of DNA extracted per sample was determined by spectrophotometric methods in a NanoDropOne (ThermoScientific).

Analysis of the 14bp Polymorphism

The region of exon 8 containing the 14bp polymorphism was amplified by PCR, using primers and conditions previously published (33, 34) (Table 2A), further confirmed using the NCBI Blast tool. Amplified products (224bp, INS variant or 210bp, DEL variant) were resolved by electrophoretic analysis in 3% agarose gels for 80 min at 90V (Figure 1A).

TABLE 2

Table 2 Primers and conditions used (A) for the 14bp INS/DEL polymorphism and (B) for the +3142 C/G polymorphism (33–35).

FIGURE 1

Figure 1 (A) Electrophoretic and (B) sequencing analysis of the 14bp INS/DEL polymorphism. For each patient, the tumoral (T) DNA sample was analyzed paired to its non-tumoral (NT) counterpart. In this example, Patient 1 is INS homozygous, Patient 2 is DEL homozygous and Patient 3 is heterozygous (INS/DEL). In all the instances the T and NT samples from each patient have the same polymorphic variant, indicating that there are no somatic mutations. INS = 224bp, DEL = 210bp. Red box= 14bp insertion. Arrow = 14bp deletion.

Analysis of the +3142 C/G Polymorphism

PCR primers and conditions used have been already published (34, 35) (Table 2B). The PCR product was subjected to BseSI digestion (PCR-RFLP) and resolved in a 2% agarose gel for 50 min at 90V. Effective digestion will disclose a “G” at position +3142, yielding two bands of 316bp and 90bp, whereas an undigested amplicon will indicate a “C” at this position and a single 406bp band (Figure 2A).

FIGURE 2

Figure 2 (A) Electrophoretic and (B) sequencing analysis of the +3142 C/G polymorphism. For each patient, the tumoral (T) DNA sample was analyzed paired to its non-tumoral (NT) counterpart. In this example, Patient 1 is G homozygous, Patient 2 is heterozygous (C/G) and Patient 3 is C homozygous. In all the instances the T and NT samples from each patient have the same polymorphic variant, indicating that there are no somatic mutations. C = 406bp, G = 316 and 90bp. Blue tint = +3142 SNP. Red box= BseSI restriction target. When a G is placed in the +3142 position the restriction occurs.

Sequencing

DNA sequencing (Sanger method) was carried out in 272 samples (107 patient paired samples plus 58 controls), to confirm the PCR and PCR-RFLP results (Figures 1B, 2B).

Statistics

The data of the polymorphisms studied by PCR or PCR-RFLP were analyzed with the software SNPStats. This software allows to assess Hardy-Weinberg Equilibrium (exact test), chi-square test, OR estimation analysis of association between polymorphisms and disease applying logistic regression models, that consider the dominant, recessive, codominant and log-additive models of inheritance. SNPStats also allows the analysis of linkage disequilibrium, using the D statistic and a correlation coefficient, and the analysis of haplotypes (EM algorithm) (36, 37).

Kaplan-Meier method was used to estimate the 5-year survival function of patients with gastric cancer and different genetic factors (GraphPad Prism 8.0 software). Multivariate Cox regression models were used to simultaneously assess the effect of genetic factors and other factors such as comorbidities, clinical features and demographic characteristics on 5-year survival of patients (software R). For all the Cox regression fits, the individual and global Schoenfeld test indicated that no covariate in the model nor the model as a whole violate the Proportional Hazard assumption, meaning that the hazard ratio stays constant over time (38).

P-values below 0.05 were considered statistically significant. As two polymorphisms were considered, the Holm-Bonferroni (HB) sequential correction method for multiple testing was applied to the statistical analyses when required. The HB method compares the k-ranked p-value to the nominal significance level (0.05) divided by (n-k+1), where in this case n = 2 (the number of polymorphisms) and k = 1 and 2.

Results

The clinical characteristics of patients are shown in Table 1.

Genetic Analysis

Hardy-Weinberg equilibrium was confirmed for both polymorphisms and groups of individuals included in the study (data not shown).

14bp INS/DEL Polymorphism

Comparing the frequency of the 14bp INS/DEL polymorphism in patients and our control group yielded statistically significant differences.

The 14bp DEL variant was the most frequent in the group of patients (70.0%) compared to the control group (57.0%, p=0.025 see Table 3A). Genotype distribution (INS/INS, INS/DEL, DEL/DEL) showed no statistically significant differences between both groups (P=0.089, Table 3B), although DEL/DEL individuals were more abundant (50.9%) in patients than in controls (34.5%).

TABLE 3

Table 3 Polymorphism 14bp INS/DEL.

The best fit model for this polymorphism, based on the Akaike Information Criterion (AIC) and Bayesian Information Criterion (BIC) values, was the log-additive model, where having each copy of DEL modifies the risk of developing gastric cancer in an additive form. DEL-allele bearers showed a higher risk of developing gastric cancer (p=0.034 OR 1.65, CI95% 1.04-2.61, see Table 3C).

As already mentioned, the DEL variant provides greater stability to the HLA-G mRNA, and, hence, yields higher levels of the protein, favoring tumor progression.

+3142 C/G Polymorphism

Although no significant difference was observed in the distribution of the alleles of this polymorphism (Table 4A), when considering the inheritance models of the different genotypes, significant differences were found. Based on the AIC and BIC values, the recessive inheritance model is the one that best fits the data obtained, with an increased frequency of C/C individuals in the group of patients (30.8%) compared to controls (14.3%, p=0.017, OR=2.68, CI95% 1.14-6.28, see Table 4C).

TABLE 4

Table 4 Polymorphism +3142 C/G.

Like the DEL variant of the previous polymorphism, the +3142 C variant would yield higher levels of the final protein, possibly favoring tumor progression.

Based on the Holm-Bonferroni sequential correction, significant associations are still detected in both +3142 (through the recessive model with the first-ranked p-value = 0.017<0.025) and 14bp (through the log-additive model with the second-ranked p-value=0.034<0.05) polymorphisms.

Haplotype Formed by the 14bp INS/DEL and +3142 C/G Polymorphisms

These two polymorphisms lie in close vicinity (182bp) and are in linkage disequilibrium, according to the values obtained with the SNPStats software (D’=0.98, r=0.73, Table 5), forming haplotypes.

TABLE 5

Table 5 Haplotype frequencies estimation and distribution.

The frequency of the 14bp-DEL/+3142-C (DEL/C) haplotype (Table 5) is higher in patients (54.1%) than in control individuals (44.4%, p=0.034, OR=1.74 CI95% 1.05-2.89). These results are consistent with the allele frequencies obtained for each polymorphism, as the variants producing higher expression of HLA-G are more frequent in patients with gastric adenocarcinoma. Therefore, this supports the results obtained at the individual level for each polymorphism and denotes DEL/C as a possible risk haplotype.

Somatic Mutations

The distribution of the polymorphisms studied was compared in paired (T and NT) gastric tissue samples in all patients. In no instance was a difference found between the paired tissue samples analyzed, irrespective of the polymorphisms (14bp INS/DEL or +3142 C/G) considered. To further confirm these data, blood samples (EDTA) from patients were drawn, and the results obtained matched those of tissues (results not shown).

Survival Curves

Five-year survival rate of patients with gastric cancer was calculated. Only patients enrolled in the study for at least 5 years (52 out of the 107 patients) were included. Patients were divided into two groups, according to their 14bp genotype: DEL/DEL and INS/DEL+ INS/INS, the former expressing more HLA-G than the latter, according to previous published work (39). Kaplan-Meier analysis (Figure 3) revealed significantly lower survival rate in patients bearing the 14bp DEL/DEL genotype (N=23) as compared to patients with either INS/DEL or INS/INS genotypes (N=29) (28% vs 54%, p=0.041). After adjustment for stage, multivariate Cox regression analysis (Figure 4) indicated a statistically significant difference in death risk (Wald test p=0.03, HR=2.6 CI95% 1.10-6.00), with an expected hazard 2.6 times higher among DEL/DEL patients as compared to patients with either INS/DEL or INS/INS genotypes. Besides, patient inclusion in stage III (Wald test p=0.018, HR 7.1 CI95% 1.39-35.56) or stage IV (Wald test p< 0.001, HR 14.5 CI95% 3.05-68.80) group was a clear statistically significant risk factor. Compared with stage I patients, the risk of death was 7.95 times higher in stage III patients and 14.48 times higher in stage IV patients.

FIGURE 3

Figure 3 5-year disease specific survival curve (Kaplan-Meier) obtained for the 14bp polymorphism (patients were grouped in two categories: DEL/DEL, N=23 vs INS/DEL+INS/INS, N=29).

FIGURE 4

Figure 4 Multivariate Cox regression carried out to assess simultaneously the effect of other factors (sex, age, cancer stage, cancer location and cancer type) on 5-year survival of patients. Our analyses indicate that polymorphism 14bp is significantly associated with 5-year survival when adjusting by cancer stage. The “*”symbol means p < 0.05, and the “***” symbol means p < 0.005.

The +3142C/G polymorphism variants yielded no significant effect on the survival of these patients (data not shown).

Discussion

The association between HLA-G, immunoediting and cancer has been of great interest in recent years (40–42). In fact, HLA-G-mediated signaling pathway is now considered as a new therapeutic immune checkpoint, in addition to other well-established ones (41).

We thus decided to analyze the involvement of HLA-G in the development of gastric adenocarcinoma. To this end, we examined in a group of 107 patients the influence of two polymorphisms (14bp INS/DEL and +3142 C/G) that affect the stability and availability of the HLA-G mRNA and, therefore, the level of the corresponding protein. A higher frequency of the variants increasing HLA-G mRNA stability (14bp DEL and +3142 C) is expected in patients, as they will favor the development and progression of gastric adenocarcinoma. Our results clearly show that HLA-G polymorphisms are linked to gastric cancer susceptibility.

Genetic Analysis

14bp INS/DEL Polymorphism

A significant difference was found in the frequency of the DEL allele in patients compared to controls (Table 3A). Since already published work (9–11) showed that this variant increases HLA-G expression, we suggest that patients bearing the DEL allele will extensively display HLA-G, allowing tumor progression. These differences were also confirmed using a larger cohort of controls of Italian origin (N=245) obtained in a bibliographic search (43) (data not shown): the DEL allele was present in 57.0% of individuals, significantly different from our diseased group (p=0.001). We used an Italian population as they have a similar HLA genetic background (44). Likewise, the 14bp DEL variant was linked to the development of breast, esophageal and colorectal cancer (39, 45–47). We may conclude that the 14bp DEL variant is a potential new gastric carcinoma risk marker.

As far as inheritance models are concerned, and according to the AIC, BIC and OR values (Table 3C), the most appropriate model to describe the distribution of these genotypes is the log-additive, whereby only the risk variant (DEL), has an effect in disease susceptibility, which, furthermore, is additive. Thus, the presence of a single copy of DEL (as in INS/DEL individuals) implies 1.65 times the risk of suffering gastric cancer, and the presence of two copies (DEL/DEL individuals; 50.9% in patients vs 34.5% in controls) will further increase this risk.

Previous published research articles revealed that the DEL allele exerts an effect on the stability and levels of HLA-G mRNA (25, 48) and on the levels of soluble HLA-G protein (26, 49), and reinforce the notion that the log-additive model is the one that best fits our study, since the high frequency of DEL/DEL patients (50.9%) implies a strong genetic predisposition of this genotype in the development of gastric cancer, tightly related to a high expression of HLA-G: the more DEL is present in the transcript the more the HLA-G mRNA is stabilized and translated.

Moreover, similar findings have been reported by other authors. Eskandari-Nasab E et al. (47) reported a higher frequency of the 14bp DEL allele and of the DEL/DEL genotype in breast cancer patients compared to the control group. Likewise, Jiang Y et al. (45) reported in a meta-analysis that the HLA-G 14bp DEL allele and the DEL/DEL genotype were associated with increased cancer risk.

+3142 C/G Polymorphism

As for the polymorphism +3142 C/G, no statistically significant differences were detected between patients and control subjects although there is a slight increase in the frequency of the C allele in patients compared to controls (Table 4A), and an increase in the frequency of the C/C genotype (Table 4B); again, this result fits in the proposed hypothesis, since bearing the C allele favors a higher expression of HLA-G (9–11, 28, 30).

In this case, the best inheritance model, based on AIC and BIC values, is the recessive one (Table 4C). According to this model, C/C individuals present 2.68 times increased risk of developing gastric cancer. Previous works described that +3142C leads to higher expression of HLA-G (28, 30), due to a lower affinity of different miRNA (148a, 148b and 152) for the HLA-G mRNA. In a similar way to 14bp DEL, which it is described to increase mRNA stability, +3142C/C genotype is overrepresented in patients (30.8% vs 14.3%, p=0.017), possibly involving this polymorphism with the development of gastric cancer.

Again, results lending support to our data have been published elsewhere. Jiang Y et al. (50) described in a meta-analysis that the HLA-G +3142 C>G mutation significantly decreased cancer risk, both in the allelic and recessive comparison models.

Haplotype Formed by the 14bp INS/DEL and +3142 C/G Polymorphisms

Disease susceptibility has long been linked to extended haplotypes of the HLA system (51, 52). The combination of the alleles of these two polymorphisms renders different haplotypes. Assuming the role the variants may have in HLA-G levels and cancer risk, the combination of the DEL and C alleles would pose the highest cancer risk. In fact, and according to the calculations carried out by the SNPStats software, the INS/G combination (considered as the reference value by the software) would be underrepresented in patients, whereas the DEL/C haplotype is significantly more frequent in patients (54.1%) than in healthy controls (44.4%) (Table 5). This indicates that the former haplotype is a protective factor, while the latter is a risk factor for this disease. The increase in the frequency of DEL/C haplotype matches the results obtained in the analysis of the individual polymorphisms, suggesting that both 14bp DEL and +3142C variants (that presumably lead to a higher HLA-G expression)are associated with susceptibility to gastric adenocarcinoma. This is, to the best of our knowledge, the first time that an association between the HLA-G 3’UTR region and the development of gastric cancer is disclosed in our population.

Somatic Mutations

Random mutations take place in cancer, conferring cells a proliferative and invasive capacity and allowing them to escape immune surveillance. In our case, an increase in variants favoring HLA-G expression and thus, abrogating immune response (i.e.: 14bp DEL and +3142C) could be expected in tumoral (T) but not in non-tumoral (NT) distal cells. However, after analysis of paired (T+NT) tissue samples from the 107 patients studied, no somatic mutations were found. The polymorphism present in a T sample analyzed matched that of the paired NT sample in every single patient tested. We can then confidently conclude that the polymorphisms here studied, and their influence on gastric cancer is a pre-existing condition in these patients.

Survival

Further to mediating disease risk, we measured whether these variants were involved in life expectancy. None of the comorbidities studied were linked to patient survival, whereas TNM staging and 14bp polymorphism revealed a clean association. Disease-specific survival rate (Figures 3 and 4) is significantly diminished, and hazard ratio increased (by up to 14-fold) in patients expressing higher levels of HLA-G according to published works (25, 26, 48, 49) (in our case, bearers of the DEL/DEL genotype). Therefore, tumor cells evade the immune system and proliferate unchecked, leading to disease dissemination and death. The 14bp DEL variant has been already associated with worse survival in a cohort of patients with colorectal cancer (39), a tumor with similar histological features to gastric adenocarcinoma, which supports the relevance of this polymorphism in the progression of cancer.

This finding suggests the possibility of considering this molecule as a potential target for therapeutic approaches. Downregulating HLA-G expression with miRNAs (as has been done in other clinical settings) (53–55) or blocking (with monoclonal antibodies) its interaction with cognate receptors (in a way similar to PD-1/PD-L1 current immunotherapies) (56–58), will make tumors visible to immunocompetent cells, eliciting an active immune response.

Although the functional effect of these polymorphisms and HLA-G expression could explain their linkage to disease susceptibility and progression, we cannot, nevertheless, exclude the possibility that these polymorphisms be in LD with other genes (i.e., other class I HLA genes) that could truly mediate the development and prognostic of this disease.

A limitation for this study is that our research is focused on genetic polymorphisms and, although we have confirmed HLA-G expression in tissue and sHLA-G in plasma in part of our cohort of patients (Supplementary Figure 1), we could not evaluate the association of HLA-G expression with the risk of developing gastric cancer or the survival of patients.

Further studies, focusing on HLA-G expression on the patients, are required to precisely assess the role of HLA-G in gastric cancer.

Besides these limitations, a more extensive cohort of patients and including other 3’UTR polymorphism (such as +3003 T/C or +3184 A/G) in a larger study would increase the reliability of the associations herein proposed. To ease comparisons between groups, future studies with larger cohorts will include ancestry informative markers (AIMs).

We conclude that the polymorphisms 14bp INS/DEL and +3142 C/G of the HLA-G gene mediate gastric cancer risk and survival, and suggest the possibility of establishing new therapeutic approaches aiming at counterbalancing the negative role of this protein in tumors.

Data Availability Statement

The nucleotide sequences have been deposited to Genbank - accession MZ130952-MZ130955.

Ethics Statement

The studies involving human participants were reviewed and approved by Comité ético de investigación clínica, Hospital Clínico San Carlos, Madrid, Spain. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

CV-Y: manuscript writing, investigation, and analysis. IJ: design, analysis, and supervision. MM-A: investigation support and validation. EM-L: analysis. AL-N: investigation support. FS-T: manuscript revision. AG-C, AL-G, IL, and RG: patient follow-up, sample and data collection. AA-V: critical review, project administration, and funding acquisition. JM-V: manuscript writing and revision, supervision, project administration, and funding acquisition. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by grants from Instituto de Salud Carlos III (PI18/00626 and PI18/00721), with funds from the European Union (Fondo Europeo de Desarrollo Regional FEDER). IJ is a grant recipient of a Universidad Complutense de Madrid—Real Colegio Complutense Harvard grant, (Ayudas para contratos predoctorales de personal investigador en formación CT18/16).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

We thank Prof. E. D. Carosella (Saint Louis Hospital, Paris, France) for scientific support and Darío Martínez Martínez for help with data analysis.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2021.698438/full#supplementary-material

Supplementary Figure 1 | Membrane bound and soluble HLA-G expression in patients with gastric cancer. (A) Representative stainings of HLA-G in tumoral (Tumor 1-4) and distal (non-tumoral) tissue of the patients included in the study. (B) Representative measure of sHLA-G expression by ELISA in the patients included in the study (1-37) and supernatants of two HLA-G negative (-) and positive (+) cell lines.

References

1. Pucułek M, Machlowska J, Wierzbicki R, Baj J, Maciejewski R, Sitarz R. Helicobacter Pylori Associated Factors in the Development of Gastric Cancer With Special Reference to the Early-Onset Subtype. Oncotarget (2018) 9:31146–62. doi: 10.18632/oncotarget.25757

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2021) 7(3):209–49. doi: 10.3322/caac.21660

CrossRef Full Text | Google Scholar

3. WHO. Cancer Today: Data Visualization Tools for Exploring the Global Cancer Burden in 2020. Available at: https://gco.iarc.fr/today/home.

Google Scholar

4. Schumacher TN, Schreiber RD. Neoantigens in Cancer Immunotherapy. Science (2015) 348:69–74. doi: 10.1126/science.aaa4971

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Urosevic M, Dummer R. Human Leukocyte Antigen–G and Cancer Immunoediting. Cancer Res (2008) 68:627–30. doi: 10.1158/0008-5472.CAN-07-2704

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Loumagne L, Baudhuin J, Favier B, Montespan F, Carosella ED, Rouas-Freiss N. In Vivo Evidence That Secretion of HLA-G by Immunogenic Tumor Cells Allows Their Evasion From Immunosurveillance. Int J Cancer (2014) 135:2107–17. doi: 10.1002/ijc.28845

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Oucherif O, Naimi D. Function of HLA-G in Cancer Immunoediting and its Clinical Benefits. African J Cancer (2015) 7:132–9. doi: 10.1007/S12558-015-0361-6

CrossRef Full Text | Google Scholar

8. Rodríguez JA. HLA-Mediated Tumor Escape Mechanisms That may Impair Immunotherapy Clinical Outcomes via T-Cell Activation. Oncol Lett (2017) 14:4415–27. doi: 10.3892/ol.2017.6784

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Alegre E, Rizzo R, Bortolotti D, Fernandez-Landázuri S, Fainardi E, González A. Some Basic Aspects of HLA-G Biology. J Immunol Res (2014) 2014:657625–10. doi: 10.1155/2014/657625

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Donadi E, Castelli E, Arnaiz-Villena A, Roger M, Rey D, Moreau P. Implications of the Polymorphism of HLA-G on its Function, Regulation, Evolution and Disease Association. Cell Mol Life Sci (2011) 68:369–95. doi: 10.1007/s00018-010-0580-7

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Arnaiz-Villena A, Juarez I, Suarez-Trujillo F, López-Nares A, Vaquero C, Palacio-Gruber J, et al. HLA-G: Function, Polymorphisms and Pathology. Int J Immunogenet (2020) 48(2):172–92. doi: 10.1111/iji.12513

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Shiroishi M, Tsumoto K, Amano K, Shirakihara Y, Colonna M, Braud VM, et al. Human Inhibitory Receptors Ig-Like Transcript 2 (ILT2) and ILT4 Compete With CD8 for MHC Class I Binding and Bind Preferentially to HLA-G. Proc Natl Acad Sci USA (2003) 100:8856–61. doi: 10.1073/pnas.1431057100

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Wan R, Wang Z, Li H, Peng X, Liu G, Ou J, et al. Human Leukocyte Antigen-G Inhibits the Anti-Tumor Effect of Natural Killer Cells via Immunoglobulin-Like Transcript 2 in Gastric Cancer. Cell Physiol Biochem (2017) 44:1828–41. doi: 10.1159/000485819

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Gao A, Sun Y, Peng G. ILT4 Functions as a Potential Checkpoint Molecule for Tumor Immunotherapy. Biochim Biophys Acta (BBA) - Rev Cancer (2018) 1869:278–85. doi: 10.1016/j.bbcan.2018.04.001

CrossRef Full Text | Google Scholar

15. Faure M, Long EO. KIR2DL4 (CD158d), an NK Cell-Activating Receptor With Inhibitory Potential. J Immunol (2002) 168:6208–14. doi: 10.4049/jimmunol.168.12.6208

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Contini P, Ghio M, Poggi A, Filaci G, Indiveri F, Ferrone S, et al. Soluble HLA-A,-B,-C and -G Molecules Induce Apoptosis in T and NK CD8+ Cells and Inhibit Cytotoxic T Cell Activity Through CD8 Ligation. Eur J Immunol (2003) 33:125–34. doi: 10.1002/immu.200390015

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Rebmann V, da Silva Nardi F, Wagner B, Horn PA. HLA-G as a Tolerogenic Molecule in Transplantation and Pregnancy. J Immunol Res (2014) 2014. doi: 10.1155/2014/297073

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Rizzo R, Bortolotti D, Bolzani S, Fainardi E. HLA-G Molecules in Autoimmune Diseases and Infections. Front Immunol (2014) 5:592. doi: 10.3389/fimmu.2014.00592

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Amiot L, Vu N, Samson M. Immunomodulatory Properties of HLA-G in Infectious Diseases. J Immunol Res (2014) 2014. doi: 10.1155/2014/298569

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Lafon M, Prehaud C, Megret F, Lafage M, Mouillot G, Roa M, et al. Modulation of HLA-G Expression in Human Neural Cells After Neurotropic Viral Infections. J Virol (2005) 79:15226–37. doi: 10.1128/JVI.79.24.15226-15237.2005

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Jiang Y, Chen S, Jia S, Zhu Z, Gao X, Dong D, et al. Association of HLA-G 3' UTR 14-Bp Insertion/Deletion Polymorphism With Hepatocellular Carcinoma Susceptibility in a Chinese Population. DNA Cell Biol (2011) 30:1027–32. doi: 10.1089/dna.2011.1238

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Dias FC, Castelli EC, Collares CVA, Moreau P, Donadi EA. The Role of HLA-G Molecule and HLA-G Gene Polymorphisms in Tumors, Viral Hepatitis, and Parasitic Diseases. Front Immunol (2015) 6:9. doi: 10.3389/fimmu.2015.00009

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Gautam S, Kumar U, Kumar M, Kanga U, Dada R. Association of HLA-G 3'utr Polymorphisms With Soluble HLA-G Levels and Disease Activity in Patients With Rheumatoid Arthritis: A Case-Control Study. Immunol Invest (2020) 49:88–105. doi: 10.1080/08820139.2019.1657146

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Castelli EC, Ramalho J, Porto IOP, Lima THA, Felício LP, Sabbagh A, et al. Insights Into HLA-G Genetics Provided by Worldwide Haplotype Diversity. Front Immunol (2014) 5:476. doi: 10.3389/fimmu.2014.00476

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Hviid TVF, Hylenius S, Rørbye C, Nielsen LG. HLA-G Allelic Variants are Associated With Differences in the HLA-G mRNA Isoform Profile and HLA-G mRNA Levels. Immunogenetics (2003) 55:63–79. doi: 10.1007/s00251-003-0547-z

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Chen X-, Yan W-, Lin A, Xu H-, Zhang J-, Wang X-. The 14 Bp Deletion Polymorphisms in HLA-G Gene Play an Important Role in the Expression of Soluble HLA-G in Plasma. Tissue Antigens (2008) 72:335–41. doi: 10.1111/j.1399-0039.2008.01107.x

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Rousseau P, Le Discorde M, Mouillot G, Marcou C, Carosella ED, Moreau P. The 14 Bp Deletion-Insertion Polymorphism in the 3' UT Region of the HLA-G Gene Influences HLA-G mRNA Stability. Hum Immunol (2003) 64:1005–10. doi: 10.1016/j.humimm.2003.08.347

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Tan Z, Randall G, Fan J, Camoretti-Mercado B, Brockman-Schneider R, Pan L, et al. Allele-Specific Targeting of microRNAs to HLA-G and Risk of Asthma. Am J Hum Genet (2007) 81:829–34. doi: 10.1086/521200

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Zhu X, Han T, Wang X, Li Y, Yang H, Luo Y, et al. Overexpression of miR-152 Leads to Reduced Expression of Human Leukocyte Antigen-G and Increased Natural Killer Cell Mediated Cytolysis in JEG-3 Cells. Am J Obstetrics Gynecol (2010) 202:592.e1–7. doi: 10.1016/j.ajog.2010.03.002

CrossRef Full Text | Google Scholar

30. Porto IOP, Mendes-Junior CT, Felício LP, Georg RC, Moreau P, Donadi EA, et al. MicroRNAs Targeting the Immunomodulatory HLA-G Gene: A New Survey Searching for microRNAs With Potential to Regulate HLA-G. Mol Immunol (2015) 65:230–41. doi: 10.1016/j.molimm.2015.01.030

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Murdaca G, Calamaro P, Lantieri F, Pigozzi S, Mastracci L, Grillo F, et al. HLA-G Expression in Gastric Carcinoma: Clinicopathological Correlations and Prognostic Impact. Virchows Arch (2018) 473:425–33. doi: 10.1007/s00428-018-2379-0

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Santiago JMR, Sasako M, Osorio J. TNM-7th Edition 2009 (UICC/AJCC) and Japanese Classification 2010 in Gastric Cancer. Towards Simplicity and Standardisation in the Management of Gastric Cancer. Cir Esp (2011) 89:275–81. doi: 10.1016/j.ciresp.2010.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Hviid TV, Hylenius S, Hoegh AM, Kruse C, Christiansen OB. HLA-G Polymorphisms in Couples With Recurrent Spontaneous Abortions. Tissue Antigens (2002) 60:122–32. doi: 10.1034/j.1399-0039.2002.600202.x

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Ben Fredj N, Sakly K, Sakly N, Bortolotti D, Aissi M, Frih-Ayed M, et al. The Association Between Functional HLA-G 14bp Insertion/Deletion and +3142 C>G Polymorphisms and Susceptibility to Multiple Sclerosis. Immunol Lett (2016) 180:24–30. doi: 10.1016/j.imlet.2016.10.006

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Cordero EAA, Veit TD, Silva M. A. L. DA, Jacques SMC, Silla LMDR, Chies JAB. HLA-G Polymorphism Influences the Susceptibility to HCV Infection in Sickle Cell Disease Patients. Tissue Antigens (2009) 74:308–13. doi: 10.1111/j.1399-0039.2009.01331.x

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Iniesta R, Guinó E, Moreno V. Análisis Estadístico De Polimorfismos Genéticos En Estudios Epidemiológicos. Gaceta Sanitaria (2005) 19:333–41. doi: 10.1157/13078029

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Solé X, Guinó E, Valls J, Iniesta R, Moreno V. SNPStats: A Web Tool for the Analysis of Association Studies. Bioinformatics (2006) 22:1928–9. doi: 10.1093/bioinformatics/btl268

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Schoenfeld D. Chi-Squared Goodness-Of-Fit Tests for the Proportional Hazards Regression Model. Biometrika (1980) 67:145–53. doi: 10.2307/2335327

CrossRef Full Text | Google Scholar

39. Garziera M, Bidoli E, Cecchin E, Mini E, Nobili S, Lonardi S, et al. HLA-G 3’utr Polymorphisms Impact the Prognosis of Stage II-III CRC Patients in Fluoropyrimidine-Based Treatment. PloS One (2015) 10:e0144000. doi: 10.1371/journal.pone.0144000

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Rouas-Freiss N, LeMaoult J, Verine J, Tronik-Le Roux D, Culine S, Hennequin C, et al. Intratumor Heterogeneity of Immune Checkpoints in Primary Renal Cell Cancer: Focus on HLA-G/Ilt2/Ilt4. Oncoimmunology (2017) 6:e1342023. doi: 10.1080/2162402X.2017.1342023

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Lin A, Yan W. Heterogeneity of HLA-G Expression in Cancers: Facing the Challenges. Front Immunol (2018) 9:2164. doi: 10.3389/fimmu.2018.02164

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Lázaro-Sánchez AD, Salces-Ortiz P, Velásquez LI, Orozco-Beltrán D, Díaz-Fernández N, Juárez-Marroquí A. HLA-G as a New Tumor Biomarker: Detection of Soluble Isoforms of HLA-G in the Serum and Saliva of Patients With Colorectal Cancer. Clin Transl Oncol (2020) 22:1166–71. doi: 10.1007/s12094-019-02244-2

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Dardano A, Rizzo R, Polini A, Stignani M, Tognini S, Pasqualetti G, et al. Soluble Human Leukocyte Antigen-G and Its Insertion/Deletion Polymorphism in Papillary Thyroid Carcinoma: Novel Potential Biomarkers of Disease? J Clin Endocrinol Metab (2012) 97:4080–6. doi: 10.1210/jc.2012-2231

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Arnaiz-Villena A, Dimitroski K, Pacho A, Moscoso J, Gómez-Casado E, Silvera-Redondo C, et al. HLA Genes in Macedonians and the Sub-Saharan Origin of the Greeks. Tissue Antigens (2001) 57:118–27. doi: 10.1034/j.1399-0039.2001.057002118.x

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Jiang Y, Lu J, Wu Y, Zhao X, Li L. Genetic Variation in the HLA-G 3′UTR 14–Bp Insertion/Deletion and the Associated Cancer Risk: Evidence From 25 Case–Control Studies. Biosci Rep (2019) 39. doi: 10.1042/BSR20181991

CrossRef Full Text | Google Scholar

46. Haghi M, Ranjbar M, Karari K, Samadi-Miandoab S, Eftekhari A, Hosseinpour-Feizi MA. Certain Haplotypes of the 3'-UTR Region of the HLA-G Gene are Linked to Breast Cancer. Br J BioMed Sci (2021) 78:87–91. doi: 10.1080/09674845.2020.1856495

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Eskandari-Nasab E, Hashemi M, Hasani S, Omrani M, Taheri M, Mashhadi M. Association Between HLA-G 3'UTR 14-Bp Ins/Del Polymorphism and Susceptibility to Breast Cancer. Cancer Biomark (2013) 13:253–9. doi: 10.3233/CBM-130364

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Verloes A, Spits C, Vercammen M, Geens M, LeMaoult J, Sermon K, et al. The Role of Methylation, DNA Polymorphisms and microRNAs on HLA-G Expression in Human Embryonic Stem Cells. Stem Cell Res (2017) 19:118–27. doi: 10.1016/j.scr.2017.01.005

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Durmanova V, Bandzuchova E, Zilinská Z, Tirpakova J, Kuba D, Buc M, et al. Impact of HLA-G 14 Bp Polymorphism and Soluble HLA-G Level on Kidney Graft Outcome. Open Life Sci (2016) 11:372–9. doi: 10.1515/biol-2016-0049

CrossRef Full Text | Google Scholar

50. Jiang Y, Li W, Lu J, Zhao X, Li L. HLA-G +3142 C>G Polymorphism and Cancer Risk: Evidence From a Meta-Analysis and Trial Sequential Analysis. Medicine (2019) 98:e16067. doi: 10.1097/MD.0000000000016067

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Tait BD, Harrison LC, Drummond BP, Stewart V, Varney MD, Honeyman MC. HLA Antigens and Age at Diagnosis of Insulin-Dependent Diabetes Mellitus. Hum Immunol (1995) 42:116–22. doi: 10.1016/0198-8859(94)00075-2

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Bitti PP, Murgia BS, Ticca A, Ferrai R, Musu L, Piras ML, et al. Association Between the Ancestral Haplotype HLA A30B18DR3 and Multiple Sclerosis in Central Sardinia. Genet Epidemiol (2001) 20:271–83. doi: 10.1002/1098-2272(200102)20:2<271::AID-GEPI9>3.0.CO;2-L

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Kaminski V, Ellwanger JH, Chies JAB. Down-Regulation of HLA-G Gene Expression as an Immunogenetic Contraceptive Therapy. Med Hypotheses (2017) 102:146–9. doi: 10.1016/j.mehy.2017.03.030

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Hanna J, Hossain GS, Kocerha J. The Potential for microRNA Therapeutics and Clinical Research. Front Genet (2019) 10:478. doi: 10.3389/fgene.2019.00478

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Fu J, Mao J, Wang C. The microRNA-152/Human Leukocyte Antigen-G Axis Affects Proliferation and Immune Escape of Non-Small Cell Lung Cancer Cells. J Int Med Res (2020) 48:300060520970758. doi: 10.1177/0300060520970758

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Hendriks D, Choi G, de Bruyn M, Wiersma VR, Bremer E. Chapter Seven - Antibody-Based Cancer Therapy: Successful Agents and Novel Approaches. In: Galluzzi L, editor. International Review of Cell and Molecular Biology. Cambridge, MA, USA: Academic Press (2017). p. 289–383. doi: 10.1016/bs.ircmb.2016.10.002

CrossRef Full Text | Google Scholar

57. Attia JVD, Dessens CE, van de Water R, Houvast RD, Kuppen PJK, Krijgsman D. The Molecular and Functional Characteristics of HLA-G and the Interaction With Its Receptors: Where to Intervene for Cancer Immunotherapy? Int J Mol Sci (2020) 21. doi: 10.3390/ijms21228678

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Carosella ED, Gregori S, Tronik-Le Roux D. HLA-G/LILRBs: A Cancer Immunotherapy Challenge. Trends Cancer (2021) 7(5):389–92. doi: 10.1016/j.trecan.2021.01.004

PubMed Abstract | CrossRef Full Text | Google Scholar