Methods: The rs4644 and rs4652 variants of galectin-3 were genotyped by TaqMan single nucleotide polymorphism assay using genomic deoxyribonucleic acid isolated from the peripheral blood of 65 (54 males, 11 females; mean age: 60.1±11.9 years; range, 34 to 83 years) with Stage IA-IIIA non-small cell lung cancer who underwent primary surgical treatment and 95 healthy individuals (48 males, 47 females; mean age: 53.9±13.5 years; range, 32 to 87 years) between March 2017 and September 2018. Circulating galectin-3 levels in serum samples of the patient and control groups were assessed by enzyme-linked immunosorbent assay. Messenger ribonucleic acid expression of galectin-3 in tumor and surrounding tissues of the patient group was examined by real-time quantitative polymerase chain reaction. Both predictive and prognostic significance of the results were analyzed.
Results: The presence of angiolymphatic invasion was significant in the patients with rs4652 AA genotype (p=0.04). Serum galectin-3 levels were significantly higher in the patients than the controls (p<0.0001). The patients with rs4644 CA/CC (p<0.0001 and p<0.0001) and rs4652 AA/AC (p=0.001 and p<0.0001) genotypes had higher serum galectin-3 levels than their corresponding controls. Serum galectin-3 levels increased in the presence of vascular invasion in patients with both rs4644 AC (p=0.03) and rs4652 AC (p=0.019) genotypes. The receiver operating characteristic curve suggested serum galectin-3 level as a strong predictive marker for the patient group with a cut-off value of 17.089 ng/mL (area under the curve: 0.910±0.04; 95% confidence interval: 0.832-0.988; p<0.001). Univariate analysis revealed the association of lower serum galectin-3 levels with better survival (p=0.048). Multivariate survival analysis showed that only high serum galectin-3 levels tended to be related to survival of the patients (hazard ratio: 5.106; 95% confidence interval: 0.956-27.267; p=0.056).
Conclusion: The presence of galectin-3 gene variants may lead to histopathological differences among patients with non-small cell lung cancer. Serum galectin-3 level may be a valuable diagnostic biomarker and be associated with survival of these patients.
Galectin-3 is an approximately 30 kDa beta-galactoside-binding protein which is encoded by the LGALS3 gene located on chromosome 14q21-q22 in human and is composed of six exons and five introns.[6] As a multifunctional protein, galectin-3 plays a pivotal role in cancer contributing to tumor progression, angiogenesis, and metastasis.[7-9] Furthermore, alterations in expression and plasma levels of galectin-3 have been considered as potential biomarkers with a clinical value, particularly in cancer.[10,11]
There are two most common single nucleotide polymorphism (SNP) sites located in exon 3 of LGALS3, namely rs4644 and rs4652 variants. The variant of rs4644 +191 C>A substitutes histidine to proline at residue 64, whereas the variant of rs4652 +292 A>C changes threonine at residue 98 to proline.[12] Double-sided roles of different variants of galectin-3 have been studied in different types of cancer. In a study, a polymorphism at rs4644 variant resulted in susceptibility to matrix metalloproteinase (MMP) cleavage and the acquisition of resistance to drug-induced apoptosis; thus, it was correlated with the incidence of breast cancer.[13] In another study, a possible protective role of rs4644 variant of galectin-3 was suggested in prostate cancer.[14] In addition, rs4652 variant, but not rs4644 variant, was associated with the tumor grade and prognosis of glioma in Chinese population.[15] Besides, rs4652 variant was thought to contribute to the susceptibility and chemotherapeutic drug resistance in gastric carcinoma.[16] Recently, Fang et al.[17] recorded that allele C of rs4652 variant might be risk factor for cervical cancer.
Despite the fact that SNPs play a key role in the occurrence and development of lung cancer, there has yet no a report showing specifically a correlation of galectin-3 gene variants with the risk and prognosis of NSCLC. In the present study, we, for the first time, aimed to investigate the association of galectin-3 and its genetic variants with the risk development and prognosis for NSCLC.
TaqMan SNP genotyping
The peripheral blood samples from all participants
were collected into the ethylenediaminetetraacetic acid
(EDTA)-containing tubes. Genomic deoxyribonucleic
acid (DNA) was isolated from peripheral circulating
lymphocytes by salting-out process.[19] The distribution
of SNPs rs4644 and rs4652 were performed using the
TaqMan SNP genotyping assays C___7593635_1_ and C___7593636_30, respectively. Briefly, approximately
50 ng of each DNA sample was subjected to a
5'-nuclease allelic discrimination assay using TaqMan
Universal PCR Master Mix (Applied Biosystems Inc.,
CA, USA) and the primer-probe mix consisting of two
SNP region-specific pre-validated primers, as well as
mutation-specific fluorescence-labeled TaqMan minor
groove binder (MGB) probes. In this study, one probe
was labeled with VIC dye to detect the allele A; the
second probe was labeled with fluorescein amidite dye
to detect the allele C. The reaction was amplified by
an initial denaturation step of 2 min at 95°C, followed
by 40 cycles of 15 sec at 95°C and 1 min at 60°C.
Endpoint analysis was performed on an ABI Prism
7900 Sequence Detection System (Applied Biosystems
Inc., CA, USA).
Enzyme-linked immunosorbent assay (ELISA)
Blood serum was obtained from participants
by centrifugation and stored -80°C. Serum levels
of galectin-3 concentrations were measured in
patients and controls by using human galectin-3
platinum ELISA kit (Invitrogen, Bender MedSystems,
Vienna, Austria) (sensitivity: 0.29 ng/mL; intraassay
coefficient of variation [CV]: 7.5% and interassay
CV: 5.4%) according to the manufacturer?s
instructions.
Real-time quantitative polymerase chain
reaction (RT-qPCR)
Total ribonucleic acid (RNA) was isolated from
the tumor and non-malignant surrounding lung
tissue specimens of the same patient by the TRizol
(Invitrogen, CA, USA) method. Validation of RNA
concentrations were performed by NanoDrop2000
(Thermo Fisher Scientific Inc., MA, USA). Firststrand
complementary DNA (cDNA) was synthesized
using high-capacity cDNA reverse transcription kit
(Applied Biosystems Inc., CA, USA). The quantitative
real-time PCR was conducted by Stratagene Mx3005P
(Agilent Technologies Inc., CA, USA) system in the
presence of SensiFAST TM SYBR® N o-ROX g reen
dye (Bioline Reagents Ldt., TN, USA). All PCR
primers used were as follow: LGALS3 forward,
5'GTGCCTCGCATGCTGATAAC3' and reverse,
5'GCAACCTTGAAGTGGTCAGG3'; GAPDH primer,
5'TGCACCACCAACTGCTTAGC3' and reverse,
5'GGCATGGACTGTGGTCATGAG3'. Gene cards
were analyzed using the threshold cycle (CT) relative
quantification method. The CT values were normalized
for endogenous reference [ΔCT=CT (GAPDH)-CT
(LGALS3)] and compared with the control using
the ΔΔCT formula [ΔΔCT=ΔCT (tumor tissue)-ΔCT (surrounding healthy tissue as control)]. Relative
messenger RNA (mRNA) expression (fold change)
was calculated by 2-ΔΔCt method (GAPDH for the
internal control). Each sample was tested in duplicate.
Statistical analysis
Statistical analysis was performed using the PASW
version 17.0 software (SPPS Inc., Chicago, IL, USA)
and GraphPad Prism version 8.0.0 (GraphPad, La
Jolla, CA, USA). The distribution of rs4644 and rs4652
genotypes and the frequency of alleles between the
patients and control groups were compared using the
chi-square, Fisher"s exact, and Student"s t-tests. Serum
galectin-3 levels and galectin-3 mRNA expression
analysis results between NSCLC patients and controls
were analyzed using the Mann-Whitney U test. The
inter-group differences of ELISA and pathological
parameters with different genotypes were compared
using the Kruskal-Wallis test. Multivariate logistic
regression analysis was adjusted by risk factors in
terms of galectin-3, sex, and smoking status. The
receiver operating characteristic (ROC) curve analysis
calculated the area under the curve (AUC) to check the
diagnostic performance of serum galectin-3 level in
NSCLC. The optimal cut-off value on ROC curve was
determined with respect to AUC value by maximizing
both specificity and sensitivity, as well as their
summation.[20] Prognostic factors were also evaluated
in NSCLC patients. Differences in survival were
investigated using the log-rank test in the univariate
analysis, and multivariate analysis was done using the
Cox proportional hazards regression model. A p value
of <0.05 was considered statistically significant.
Table 1: Baseline characteristics of study population
The clinical features of the patients with genotypes of rs4644 and rs4652 variants were presented in Table 3. Regarding disease`s characteristics including tumor stage, lymph node metastasis, perineural invasion, angiolymphatic invasion and vascular invasion, no association was observed among NSCLC patients with genotypes of rs4644 variant. However, there was a markedly suggestive evidence of the presence of angiolymphatic invasion for NSCLC patients who had AA genotype of rs4652 variant (p=0.04).
The mean serum galectin-3 levels were significantly higher in NSCLC patients (26.05±1.77 ng/mL), compared to healthy controls (11.62±1.30 ng/mL) (p<0.0001) (Figure 1a). The patients with CA and CC genotypes of rs4644 variant were associated a higher serum galectin-3 level compared to the controls (for CA genotype: 26.02±1.72 vs. 12.08±2.05 ng/mL, respectively; p<0.0001; for CC genotype: 28.8±3.43 vs. 12.27±1.80 ng/mL; respectively; p<0.0001) (Figure 1b). Moreover, the patients with AC and AA genotypes of rs4652 variant were correlated with higher serum galectin-3 levels, compared to controls (for AC genotype: 26.50±1.86 vs. 12.08±1.76 ng/mL, respectively; p<0.0001; for AA genotype: 28.76±3.43 vs. 1 1.70±2.36 n g/mL, r espectively; p =0.001) (Figure 1c). Moreover, serum galectin-3 levels of NSCLC patients with AC genotype of both rs4644 and rs4652 variants were found to be higher in the presence of vascular invasion than in the absence of it (for rs4644 variant: 30.81±2.04 ng/mL vs. 23.03±1.85 ng/mL, respectively; p=0.03; for rs4652 variant: 31.94±2.20 ng/mL vs. 22.87±1.37 ng/mL, respectively; p=0.019, respectively) (Tables 4 and 5).
Multivariate logistic regression analysis was also tested for a statistical interaction in which the association of serum galectin-3 level, sex, and smoking status with NSCLC occurrence. It was shown that serum galectin-3 level (odds ratio [OR]: 57.431; 95% confidence interval [CI]: 9.085-363.046; p=0.000) and smoking (OR: 21.730; 95% CI: 1.357-348.056; p=0.03) could be predictive risk factors of NSCLC influencing the diagnosis independent from variant distribution of galectin-3 (Table 6).
Table 6: Multivariate logistic regression analysis for the risk of development of NSCLC
To further evaluate the diagnostic value of serum galectin-3 level in NSCLC patients, the ROC curve was plotted (Figure 2). The cut-off value of serum galectin-3 level was calculated as 17.089 ng/mL with the maximal sensitivity (88%) and specificity (84%) (AUC: 0.910±0.04; 95% CI: 0.832-0.988; p<0.001).
Galectin-3 mRNA expression level was determined in the tumor tissue, compared to the corresponding surrounding healthy tissue of the same participant for each NSCLC patient by the RT-qPCR. The mean galectin-3 mRNA expression was higher in surrounding healthy tissue referred to as control tissue (0.09±0.039 fold) (0.09±0.039 fold) than tumor tissue (0.04± 0.012), although this difference was not statistically significant (p=0.4959) (Figure 3).
The five-year overall survival rate of operated NSCLC patients was 68.2% and the median survival time was 62.7 months (95% CI: 53.4-72.1). While analyzing survival using log-rank comparison of survival data with regards to genotypes of galectin-3 gene variants, there was no statistically significantly different survival rates between the patients with genotypes of rs4644 variant (p=0.625) (Figure 4). Similarly, no statistically significant difference in the survival rates was found between patients with regards to rs4652 genotypes (p=0.349) (Figure 5). In terms of galectin-3 serum levels, the five-year survival rate of patients who had higher serum galectin-3 levels (>17.089 ng/mL) was 44.4%, whereas it was 64.3% in the patients with lower serum galectin-3 levels (p=0.048) (Figure 6). The N factor, T factor, angiolymphatic involvement, and serum galectin-3 levels were subjected to the multivariate Cox proportional analysis. It revealed that only high levels of serum galectin-3 tended to be associated with survival of the patients (hazard ratio [HR]=5.106; 95% CI: 0.956-27.267; p=0.056) (Table 7).
Table 7: Multivariate cox regression model for death hazard estimation in NSCLC patients
In terms of SNPs, the distribution of rs4644 and rs4652 variants were not found to be exclusive between the patients and controls. However, rs4652 variant of galectin-3 was reported to be associated with platinumbased chemotherapy response and prognosis in patients with NSCLC rather than rs4644 variant.[22] However, we did not analyze the response to chemotherapy of our patients, since all patients underwent resection surgery. Regarding the possible roles of rs4644 and rs4652 variants in the incidence of several cancer types,[13-17] the sample size of the present study could be enlarged to reach statistical significance for the genotype and allele distribution between the patients and controls.
In our study, the frequency of patients with AA homozygous genotype of rs4652 variant was significantly increased in the presence of angiolymphatic invasion. In the studies including NSCLC patients with Stage I-IIA, having angiolymphatic invasion was found to be a negative prognostic factor for the development of long-term survival and relapse.[23] However, we obtained no association between the genotypes of galectin-3 variants (rs4644 and rs4652) and survival in NSCLC patients.
Serum galectin-3 levels were higher in patients with AC and CC genotypes of rs4644 variant than in healthy controls. Moreover, rs4652 variant showed higher serum galectin-3 levels in patients with AA and AC genotypes, compared to the control group. There has been no study exemplifying the distribution of serum galectin-3 levels with respect to genotypes of galectin-3 variants in NSCLC patients. However, it has been evidenced that polymorphic amino acid changes in both the rs4644 and rs4652 positions along with the distribution of genotypes affect the secretion of serum galectin-3 in the body through taking an active role in various pathological processes which contribute to the cancer and other diseases.[12] Nevertheless, the exact mechanism of how amino acid change alters the protein expression and function is still unclear. It is known that rs4644 in the 64th residue which harbors proline than histidine has a more stable structure which is evolutionarily conserved, any change in this sequence is likely to affect the function of this molecule.[24] Besides, proline amino acid at the 98th residue of rs4652 variant is also known as critical in protein transportation.[22] Studies on hamsters, N-terminal sequence containing residues 89-96 (Tyr-Pro-Ser-Ala-Pro-Gly-Ala-Tyr) in the galectin-3 gene was found to play a critical role in protein secretion.[25] Also, N-terminal region of galectin-3 is evolutionary conserved in mammalians.[22] Therefore, human N-terminal region located on the human rs4652 variant of galectin-3, which is highly homologous to hamsters, may alter the intracellular galectin-3 level in the body.
Considering the pathological features of patients, in the presence of vascular invasion, serum galectin-3 levels were found to be higher in patients with AC genotypes of both rs4644 and rs4652 variants. The main molecular mechanisms proposed for the galectin-3 promoted angiogenesis include the binding of galectin-3 carbohydrate recognition domain (CRD) to αv3 integrins on endothelial cells, and subsequently inducing the signaling pathways related to angiogenic activity such as vascular endothelial growth factor, fibroblast growth factor, activation of focal adhesion kinase phosphorylation.[8] Shi et al.[16] hypothesized that an A>C alteration at rs4652 site located on the CRD of galectin-3 might strengthen the interaction between CRD and fibroblast surface receptors, thereby, leading to pathway the myofibroblast promotion of angiogenesis and the growth of gastric cancer cells,[16] as activated fibroblasts, also known as myofibroblasts, were thought to increase angiogenesis by interacting with cancer cells.[26]
Multivariate regression analyses have confirmed that smoking contributed to the lung cancer risk matching variables such as age and sex.[27] To extend the work of literature, serum galectin-3 level was also identified as a promising predictive risk factor for NSCLC in our study. Furthermore, the cut-off value of serum galectin-3 level (17.089 ng/mL) was proposed for assessing the discriminatory accuracy between NSCLC patients and healthy individuals.
Finally, in the survival analysis, we found that rs4644 and rs4652 genotypes of galectin-3 were not associated with survival. However, we found that high serum galectin-3 level was marginally associated with lower survival among NSCLC patients. Multivariate analysis indicated that high serum galectin-3 level tended to be an independent surrogate marker for worse survival. This finding may help to select patients for possible adjuvant, despite having earlystage tumors. It is also plausible to propose that the patients with higher serum galectin-3 levels could be followed more frequently for possible earlier recurrences. Similarly, Kataoka et al.[28] found that higher expressions of galectin-3 in tumor cells were associated with tumor recurrence. However, they did not perform multivariate analysis to prove that whether high galectin-3 level was independently related with worse survival.
The main limitations to our study include relatively low number of patients and the matching problems for sample from the same study population, resulting from the exclusion of certain patients, particularly for the expression analysis. Therefore, a heterogeneous gene expression profile of galectin-3 was obtained among tumor tissues of NSCLC patients. A larger series with the analysis of higher number of patients could have revealed a statistically significant association between serum galectin-3 level and survival. Thus, the outcomes of this study warrant further confirmation with larger groups.
In conclusion, the present study may provide a basis for the literature investigating the role of galectin-3 variants in lung cancer by shedding light into whether galectin-3 may be one of the important biomarkers to be suggested in disease risk and progression. Serum galectin-3 level may be utilized to better identify patients with early-stage lung cancer who can benefit from adjuvant therapy after surgical resection. Further in vitro studies are needed to elucidate cellular functional changes in cancer signaling pathways related to rs4644 and rs4652 variants of galectin-3.
Declaration of conflicting interests
The authors declared no conflicts of interest with respect to
the authorship and/or publication of this article.
Funding
This study was supported by the grant from Istanbul
University Scientific Research Found (IU BAP; TDK-2018-
27977).
1) Siegel RL, Miller KD, Jemal A. Cancer statistics, 2019. CA
Cancer J Clin 2019;69:7-34.
2) Goldstraw P, Crowley J, Chansky K, Giroux DJ, Groome PA,
Rami-Porta R, et al. The IASLC Lung Cancer Staging Project:
Proposals for the revision of the TNM stage groupings in the
forthcoming (seventh) edition of the TNM classification of
malignant tumours. J Thorac Oncol 2007;2:706-14.
3) Liao Z, Lin SH, Cox JD. Status of particle therapy for lung
cancer. Acta Oncol 2011;50:745-56.
4) Detterbeck FC, Boffa DJ, Tanoue LT. The new lung cancer
staging system. Chest 2009;136:260-71.
5) Chang A. Chemotherapy, chemoresistance and the
changing treatment landscape for NSCLC. Lung Cancer
2011;71:3-10.
6) Kadrofske MM, Openo KP, Wang JL. The human LGALS3
(galectin-3) gene: Determination of the gene structure and
functional characterization of the promoter. Arch Biochem
Biophys 1998;349:7-20.
7) Liu FT, Rabinovich GA. Galectins as modulators of tumour
progression. Nat Rev Cancer 2005;5:29-41.
8) Markowska AI, Liu FT, Panjwani N. Galectin-3 is an
important mediator of VEGF- and bFGF-mediated
angiogenic response. J Exp Med 2010;207:1981-93.
9) Takenaka Y, Fukumori T, Raz A. Galectin-3 and metastasis.
Glycoconj J 2002;19:543-9.
10) Saussez S, Lorfevre F, Lequeux T, Laurent G, Chantrain
G, Vertongen F, et al. The determination of the levels of
circulating galectin-1 and -3 in HNSCC patients could be
used to monitor tumor progression and/or responses to
therapy. Oral Oncol 2008;44:86-93.
11) Chung LY, Tang SJ, Wu YC, Sun GH, Liu HY, Sun
KH. Galectin-3 augments tumor initiating property and
tumorigenicity of lung cancer through interaction with
?-catenin. Oncotarget 2015;6:4936-52.
12) Hu CY, Chang SK, Wu CS, Tsai WI, Hsu PN. Galectin-3
gene (LGALS3) +292C allele is a genetic predisposition
factor for rheumatoid arthritis in Taiwan. Clin Rheumatol
2011;30:1227-33.
13) Balan V, Nangia-Makker P, Schwartz AG, Jung YS, Tait
L, Hogan V, et al. Racial disparity in breast cancer and
functional germ line mutation in galectin-3 (rs4644): A pilot
study. Cancer Res 2008;68:10045-50.
14) Meyer A, Coinac I, Bogdanova N, Dubrowinskaja N,
Turmanov N, Haubold S, et al. Apoptosis gene polymorphisms
and risk of prostate cancer: A hospital-based study of
German patients treated with brachytherapy. Urol Oncol
2013;31:74-81.
15) Chen HJ, Zheng ZC, Yuan BQ, Liu Z, Jing J, Wang SS. The
effect of galectin-3 genetic variants on the susceptibility and
prognosis of gliomas in a Chinese population. Neurosci Lett
2012;518:1-4.
16) Shi Y, Lin X, Chen G, Yan J, Ying M, Zheng X. Galectin-3
rs4652 A>C polymorphism is associated with the risk of
gastric carcinoma and P-glycoprotein expression level. Oncol
Lett 2017;14:8144-9.
17) Fang SQ, Feng YM, Li M. Correlations of galectin-3 gene
polymorphisms with risk and prognosis of cervical cancer in
Chinese populations: A case-control study. Oncol Res Treat
2017;40:533-9.
18) Groome PA, Bolejack V, Crowley JJ, Kennedy C, Krasnik
M, Sobin LH, et al. The IASLC Lung Cancer Staging
Project: validation of the proposals for revision of
the T, N, and M descriptors and consequent stage
groupings in the forthcoming (seventh) edition of the
TNM classification of malignant tumours. J Thorac
Oncol 2007;2:694-705.
19) Miller SA, Dykes DD, Polesky HF. A simple salting out
procedure for extracting DNA from human nucleated cells.
Nucleic Acids Res 1988;16:1215.
20) Akobeng AK. Understanding diagnostic tests 3:
Receiver operating characteristic curves. Acta Paediatr
2007;96:644-7.
21) Kacan M, Turkyilmaz M, Karabulut F, Altun O, Baran Y.
Complexation, thermal and catalytic studies of N-substituted
piperazine, morpholine and thiomorpholine with some
metal ions. Spectrochim Acta A Mol Biomol Spectrosc
2014;118:572-7.
22) Wu F, Hu N, Li Y, Bian B, Xu G, Zheng Y. Galectin-3 genetic
variants are associated with platinum-based chemotherapy
response and prognosis in patients with NSCLC. Cell Oncol
(Dordr) 2012;35:175-80.
23) Sung SY, Kwak YK, Lee SW, Jo IY, Park JK, Kim KS, et
al. Lymphovascular invasion increases the risk of nodal and
distant recurrence in node-negative stage I-IIA non-smallcell
lung cancer. Oncology 2018;95:156-62.
24) Kaur T, Sodhi A, Singh J, Arora S, Kamboj SS, Kaur M.
Evaluation of galectin-3 genetic variants and its serum levels
in rheumatoid arthritis in North India. International Journal
of Human Genetics 2015;15:131-8.
25) Menon RP, Hughes RC. Determinants in the N-terminal
domains of galectin-3 for secretion by a novel pathway
circumventing the endoplasmic reticulum-Golgi complex.
Eur J Biochem 1999;264:569-76.
26) Webber JP, Spary LK, Sanders AJ, Chowdhury R, Jiang
WG, Steadman R, et al. Differentiation of tumour-promoting
stromal myofibroblasts by cancer exosomes. Oncogene
2015;34:290-302.