Abstract
Background/Aim: The matrix metalloproteinases (MMP) play an important role in the physiological and pathological remodeling of tissues including carcinogenesis. The study's aim was to assess the relations between MMP-2(-735C/T), MMP-2(-1306C/T), MMP-9(-1562C/T), and MMP-3(-11715A/6A) polymorphisms, and clinical/morphological manifestation of laryngeal squamous cell carcinoma (LSCC) and benign vocal fold lesions (BVFL). Patients and Methods: Two hundred and seventeen patients with LSCC and BVFL and 458 controls were included in this study. The genotyping was performed using the real-time polymerase chain reaction method. Results: The MMP-2(-1306C/T) C/T genotype was significantly rarer among the patients with moderate-poorly differentiated LSCC compared to the control group, however the MMP-3(-11715A/6A) 6A/6A genotype was significantly more frequent compared to controls. Smoking and 6A/6A genotype of MMP-3(-11715A/6A) polymorphism were associated with increased odds of LSCC risk. No associations between MMP genotypes and BVFL were found. Conclusion: Smoking and MMP-3 (-11715A/6A) 6A/6A genotype may cause a higher risk for developing LSCC.
- Laryngeal squamous cell carcinoma
- matrix metalloproteinases
- gene polymorphism
- benign vocal fold lesions
Head and neck squamous cell carcinoma (HNSCC), which includes tumors of the oral cavity, pharynx, hypopharynx, and larynx is one of the most common cancers worldwide (1). Laryngeal squamous cell carcinoma (LSCC) constitutes approximately 2-3% of all malignant tumors representing the most common neoplasm of the head and neck region. According to the morbidity rate, LSCC ranks second among the neoplasms of the respiratory system after lung cancer (2). Despite the current advances in diagnostics, surgery, radiotherapy and chemotherapy, there was no significant increase in the 5-year survival rate in the last few decades (3). Following LSCC treatment, the quality of life of patients in advanced stages of the disease often remains significantly impaired, manifesting with dysphagia and loss of voice and speech (4).
The laryngeal carcinoma does not occur in unaltered tissues. It is formed under the influence of harmful habits when hyperplastic (dysplasia, atypia) changes emerge in squamous cell epithelium. Many factors, such as smoking, alcohol use, HPV infection, and genetic predisposition are associated with an increased risk of LSCC (5). Although tobacco and alcohol addiction play a critical role in LSCC carcinogenesis, only a small proportion of smokers and drinkers are finally diagnosed with LSCC. This leads to presuming that genetic susceptibility to LSCC may vary in the general population (6).
The important role of matrix metalloproteinases (MMPs) in HNSCC development has been established in the last decades (7). MMPs belong to a family of zinc-containing enzymes which are responsible for degrading all extracellular matrix components, which is a very important event during the cancer cell invasion process and metastasis of most neoplasms (8). MMPs represent the endopeptidases that contain an active site Zn2+ and are divided into subfamilies based on relationships and structure of the catalytic domain. MMPs encompass a family of currently 25 related, yet distinct vertebrate gene products of which 23 are found in humans (9). Under normal conditions, MMPs are involved in both wound repair, tissue regeneration, and reproduction (10). MMPs may also be implicated in the carcinogenesis process, as previous studies have shown that MMPs may play a role in several steps of cancer development, including cancer cell growth, differentiation, migration, apoptosis, invasion, and metastasis (11). Typically, MMPs may be divided into the following five groups: i) collagenases (MMP-1,-8,-13), ii) gelatinases (MMP-2,-9), iii) stromelysins (MMP-3,-10), iv) membrane MMPs (MMP-14,-15,-16,-17), and v) others (MMP-19,-21,-23A,-23B,-27,-28) (9).
MMP-2, MMP-9, and MMP-3 represent the key members of the MMPs family and have been found to play a crucial role in physiological and pathological events associated with the metabolism and remodelling of the connective tissue (12, 13). Overexpression of these MMPs has been found to be associated with the development of different cancer types, including HNSCC (7), thereby indicating that these MMPs may also be implicated in the development of LSCC. The role of MMP-2 and MMP-9 in LSCC tumorigenesis has been investigated in numerous studies (14-16), however there is only one study analysing laryngeal carcinoma association with MMP-3 (17). To the best of our knowledge, there is no literature data analysing MMP-2 (-735 C/T), MMP-2 (-1306 C/T), MMP-9 (-1562 C/T), and MMP-3 (-1171 5A/6A) gene polymorphism associations with LSCC and non-malignant laryngeal lesions.
The aim of the study was to assess the relations between MMP-2 (-735 C/T), MMP-2 (-1306 C/T), MMP-9 (-1562 C/T), and MMP-3 (-1171 5A/6A) gene polymorphisms and the LSCC and non-malignant laryngeal lesions as well as their clinical and morphologic manifestations.
Patients and Methods
Patients. The research protocol of the present study was approved by the Kaunas Regional Biomedical Research Ethics Committee (Protocol No. BE-2-34). The total study group consisted of 217 LSCC and benign vocal fold lesions (BVFL) patients who underwent treatment at the Department of Otolaryngology of Lithuanian University of Health Sciences (LUHS) and the Oncological Hospital of the LUHS.
The BVFL group consisted of 105 patients and included the following clinical entities: i) vocal fold polyps (VFP) (n=39, age: 43, range=21-69), ii) laryngeal keratosis (LK) (n=36, age:54, range=34-78), and recurrent respiratory papillomatosis (RRP) subgroup (n=30, age:38, range=21-83). The clinical diagnosis of BVFL was based on typical clinical signs revealed during video laryngostroboscopy (VLS) and direct microlaryngoscopy. All patients underwent endolaryngeal microsurgery; therefore, the final diagnosis was confirmed using the results of histological examination of the removed mass lesions from the larynx.
The LSCC group consisted of 105 males and 7 females who received no treatment prior to being diagnosed with LSCC (n=112). The age of LSCC group patients ranged from 36 to 79 years (median 62 years). Clinical diagnosis of LSCC was based on the findings of VLS and direct microlaryngoscopy. The final diagnosis was confirmed following histopathological investigation. The patients of the LSCC group underwent standard additional diagnostic tests, including: chest x-ray, abdomen ultrasonography, and computed tomography (CT) of the neck. The subjects with previously diagnosed malignant tumors in other body locations were not involved into the study. The most important clinical and pathological features of BVFL and LSCC patients are presented in Tables I and II.
The reference group consisted of random research samples from Kaunas city (Lithuania) residents who were selected during the HAPPIE (The international Health, Alcohol and Psychosocial Factors in Eastern Europe project) study (18). This reference group was used in the present study to compare the genotypic distribution of MMP gene polymorphisms. The samples were collected in the Population Study Laboratory at the Institute of Cardiology of LUHS by carrying out a selection based on the lists of residents of Kaunas city. Totally, 1128 residents of Kaunas within the age group from 45 to 73 years (median: 60 years) participated in the HAPPIE study. The reference group for this study was formed by selecting 428 men and 30 women within the age group from 45 to 73 years (median: 62 years) from the HAPPIE study. The reference group was created by taking into consideration the distribution of age and gender in the LSCC group. Therefore, the medians of the patient age of the reference group and the LSCC group did not differ significantly (p>0.05).
DNA extraction. The DNA extraction and genotyping were carried out at the Laboratory of Molecular Cardiology at the Institute of Cardiology of LUHS. Briefly, 10 ml of blood were withdrawn from a peripheral vein for each study participant, and DNA was extracted using the Genomic DNA Purification Kit (Fermentas, Vilnius, Lithuania) according to the recommendations of the manufacturer or the silica gel column method utilizing the genomic DNA extraction kit SorpoClean™ Genomic DNA Extraction Module (SORPO Diagnostics, Vilnius, Lithuania) according to the manufacturer's protocol (Metabion, Planegg, Germany). The genotyping of MMP-3 (−1171 5A/6A), MMP-2 (−735 C/T), MMP-2 (−1306 C/T), and MMP-9 (−1562 C/T) was carried out using the real-time polymerase chain reaction (PCR) method. For the genotyping of MMP-2 (−735 C/T) (rs2285053) and MMP-2 (−1306 C/T) (rs243865) Applied Biosystems (Thermo Fisher Scientific, San Jose, CA, USA) kits were used. Genotyping was performed using the HT 7900 real-time PCR quantification system (Applied Biosystems). The real-time PCR reagents (2X Maxima™ Probe/ROX qPCR Master mix buffer, fluorescent dye-labelled markers, sterile ddH2O) were taken out of −20°C and were thawed at room temperature. The thawed reagents were centrifuged (10,000 rpm) and were stored in an ice tub for setting up the RT PCR reactions. The PCR mixture (9 μl) was poured into each well of a 96-well microtiter plate, and then 1 μl of matrix DNA of samples (~10 ng) and 1 μl of a negative control (−K) were added. An optic film was pasted on the microtiter with the 96 wells, and the microtiter was centrifuged for 15 s at 10,000 rpm. The remaining procedure was continued according to the manual provided by the manufacturer.
A detailed description of the genotyping methods used in the present research have been presented in previous studies (19, 20).
Statistical analysis. For our statistical analysis, the program PASW Statistics 20.0 (SPSS Inc., Chicago, IL, USA) was used. The Hardy – Weinberg analysis was performed to compare the observed and expected MMP's genotype frequencies using χ2 test in all groups. The distribution of the MMP's single nucleotide polymorphism (SNP) between the patients and the reference groups were compared using either χ2 test or Fisher's exact test. Relation of smoking risk and genotype in the LSCC patients' group was analysed by odds ratios (OR) and 95% confidence intervals (95%CI). The significance level of 0.05 was chosen for testing statistical hypotheses.
Results
The distributions of MMP-3 (-1171 5A/6A), MMP-2 (−735 C/T), MMP-2 (-1306 C/T), and MMP-9 (−1562 C/T) genotypes promoter polymorphisms in BVFL, LSCC and reference groups were in agreement with the Hardy – Weinberg equilibrium (χ2<3.84, p>0.05).
The distribution of MMP-3 (-1171 5A/6A), MMP-2 (−735 C/T), MMP-2 (-1306 C/T), and MMP-9 (−1562 C/T) gene polymorphisms among the patients from BVFL group did not differ significantly (p>0.05) from the reference group. The MMP-2 (-1306 C/T) C/T genotype was significantly more rare among the LSCC patients compared to the reference group (p=0.04). However, in the LSCC patients group only a tendency of more frequent MMP-3 (-1171 5A/6A) 6A/6A genotype polymorphism was revealed compared to the reference group (p=0.059). For the patients who were diagnosed with advanced stages of LSCC (III-IV), their distribution of genotype MMP-3 (-1171 5A/6A) showed a significantly more frequent 6A/6A genotype (p=0.004) and a more frequent presence of the 6A allele (p=0.041) compared to the genotype and allele distribution of the reference group. A similar outcome was revealed after determining a significantly more frequent 6A/6A genotype in the moderate and poorly (G2-3) differentiated LSCC group compared to the reference group (p=0.017). The MMP-2 (-1306 C/T) C/T genotype was significantly more rare among the patients with moderate and poorly (G2-3) differentiated LSCC compared to the reference group (p=0.007). Compared to genotypes of the reference group, MMP-2 (-735 C/T) and MMP-9 (-1562 C/T) genotypes of the LSCC group did not differ. The results of genotypic distributions among LSCC and reference groups are presented in Table III.
Having applied the logistic regression model, the possibility to develop LSCC appeared to significantly depend on the habit of smoking and the MMP-3 (-1171 5A/6A) genotype 6A/6A. The possibility of the LSCC development for smoking individuals was 7.69 fold higher compared to non-smokers (OR=7.69, 95%CI=4.63-1.77, p=0.0001), and for those with MMP-3 6A/6A genotype it increases 1.78 fold compared to the 5A/5A or 5A/6A genotypes (OR=1.78 95%CI=1.08-2.93, p=0.022).
Discussion
Expression of MMPs has already been investigated in numerous studies and in a large variety of different tumors, including different locations of squamous cell carcinoma. Nevertheless, in LSCC, the correlations between the expression of MMPs, morphological parameters of the tumor and clinical signs still remain controversial; moreover, prognostic roles of different MMPs in LSCC have not been established as yet (21-25).
Activation of MMPs expression can be caused by various gene polymorphisms in promoter regions, when their binding to transcription factors or other regulating elements is disrupted. MMP polymorphisms can be induced by nucleotide changes in promoter regions by insertions, substitutions or microsatellite instability (12). In around 90% of cases, a SNP is determined in the locus where one of the basic changes appears in the DNA strain. Polymorphisms can also be detected in MMP gene promoter regions but most of them are not biologically active. Only a small part of SNPs that can change gene transcription intensity is biologically active and can have an impact on the predisposition of a disease (26).
It has been previously established that MMP-2 -1306 C/T transition in promoter regions consequently changes their activity. Experiments provided by Price’ team have shown that the C allele increases MMP-2 expression more significantly compared to the T allele (27). Similarly, C/T transition located in the MMP-2 -735 promoter region can lead to the same effect (28). Yu et al., have studied Chinese patients with esophageal squamous cell carcinoma carrying the -1306 C/T and -735 C/T polymorphisms in the MMP-2 promoter region. The authors describe the direct important associations between the haplotype C (-1306)-C (-735) and an increased risk of esophageal cancer (29). Penget et al., have confirmed this association in a meta-analysis study (30). In this study, the MMP-2 -735 CC genotype carriers had a higher risk of esophageal and lung cancer compared to those with TT and CT genotypes; the MMP-2 -1306 CC genotype carriers were in higher risk of esophageal, lung, gastric, and head-neck cancer comparing to TT and CT carriers (30). In 2015, Kim et al. conducted a large meta-analysis of sixteen case-control studies with 11792 individuals of MMP-2 (-735 C/T) SNP. The results of this study demonstrated significant associations of the C/C genotype with increased cancer risk in the Caucasian population, however, overall MMP-2 (-735 C/T) polymorphism has not been found to be associated with cancer risk in any calculated model (31). In the present study, analysis of MMP-2 (-735 C/T) genotypes distribution did not show any statistically significant relations with LSCC risk, however the MMP-2 (-1306) C/T genotype was significant rarer in LSCC patients compared to healthy controls.
The MMP-9 (-1562 C/T) polymorphism was described for the first time by Zhang et al. in 1999. The authors found that the C to T allele substitution in the MMP-9 promotor region may increase the promotor activity by 1.5-fold (32). In later study, increased T allele frequency in 152 oral squamous cell carcinomas (OSCC) was detected by Vairaktaris et al. (33). Tu et al. from Taiwan, have investigated the possible associations between the MMP-9 (-1562 C/T) polymorphism, OSCC and oral submucosal fibrosis (OSF) in areca users. The authors concluded that the T allele influenced the OSCC risk only in younger areca chewers (34). Similarly, in 2011, Chaudhary et al. described their results from a genetic investigation in OSCC and OSF patients in an Indian population. They evaluated MMP-2 (-1306 C/T) and MMP-9 (-1562 C/T) polymorphisms in 1260 individuals; 422 patients with HNSCC, 412 patients with OSF, and 426 controls, and identified a strong association between the T allele of both gene polymorphisms and HNSCC risk in comparison to controls, however, this was not the case with OSF (35). The interesting results were presented by Liu et al. in meta-analysis study. Significant associations between the MMP-9 (-1306 C/T) polymorphism and the risk of different cancer types metastases were found only with the C/T genotype. The results from the Asian population showed that the C/T and T/T genotypes had a higher risk of cancer metastases in comparison to C/C. In contrast, no such association was found in the European population (36). On the other hand, in study done by Nasr et al. in 2007, no significant association between the MMP-9 (-1306 C/T) polymorphism and nasopharyngeal carcinoma was found (37). This is in agreement with the results of the present study, where the distribution of the MMP-9 (-1306 C/T) genotypes in patients with LSCC did not differ significantly compared to the reference group.
The literature concerning the effect of both MMP-3 expression and its polymorphisms on HNSCC carcinogenesis is scarce and inconsistent. Until now only a few studies have investigated the MMP-3 (-1171 5A/6A) polymorphisms concerning the different head and neck region tumors. A common variant in the promoter region of the human MMP-3 gene with 1 allele having a run of 5 adenines (5A) and the other having 6 adenines (6A) have been revealed as manifesting an impact on gene expression (38). MMP-3 gene is located in chromosome 11 (11q22.2-11q22.3 region). Insertion of one adenine (A) in the -1171 base-pair position of MMP-3 promoter causes a formation of 6 adenines (6A) instead of 5 (5A). It has been shown that the 6A allele has a higher binding affinity for the ZBP-89 transcription factor, which decreases both promoter transcription activity and MMP-3 gene expression (38). Chaudhary et al., have evaluated the MMP-3 (-1171 5A/6A) polymorphism in 362 participants: i) 101 cases of OSF, ii) 135 cases HNSCC and iii) 126 healthy controls. This study found that the 5A allele might play an important role in the susceptibility to HNSCC, as individuals with the 5A/5A genotype had nearly a two-fold risk of HNSCC (OR=1.94) when compared to controls (39). Tu et al., have analysed OSF and OSCC, which are prevalent among the Asian areca users, and found that the 5A/5A genotype was associated with high risk of OSF but not of OSCC (40).
Furthermore, in studies by Hashimoto et al. and Nishizawa et al., no significant association between the MMP-3 (-1171 5A/6A) polymorphism and a risk of HNSCC was found (41-42). On the other hand, Zhang et al., have performed a meta-analysis, that included thirteen studies analysing associations between MMP-2, 3 and 9 promoter polymorphisms and ethnic, tumor site, staging and other clinical-morphological features. For MMP-3 (-1171 5A/6A) and MMP-9 (-1562 C/T) polymorphisms, no association was found in total meta-analysis cohort; however, in subgroup analyses based on ethnicity and tumor site, significant associations were detected between the MMP-3 (-1171 5A/6A) polymorphism and HNSCC risk in a European population and pharyngeal/laryngeal cancer. Authors have drawn attention to urge for completion of future studies with larger samples that should further evaluate the role of the MMP-3 (-1171 5A/6A) polymorphism on HNSCC risk in different populations (43).
To the best of our knowledge, the associations between MMP-3 (-1171 5A/6A) polymorphism and single LSCC have not been investigated before. Therefore, in the present study, the MMP-3 (-1171 5A/6A) gene polymorphisms and their associations with clinical and morphological manifestations in patients with LSCC and non-malignant laryngeal lesions were investigated for the first time. Our study revealed that the possibility to acquire LSCC for smokers is 7.69 times higher compared to non-smokers and for those who have MMP-3 6A/6A genotype this chance increases by 1.78 times compared to subjects who have 5A/5A or 5A/6A genotypes. Also, in patients with advanced LSCC the distribution of genotype MMP-3 (-1171 5A/6A) showed a significantly more frequent presence of the 6A/6A genotype and a more frequent allele 6A compared to the genotype and allele distribution of the reference group. A similar outcome was revealed after determining a significantly more frequent 6A/6A and a rarer 5A/5A genotype in the moderate-poorly (G2-3) differentiated LSCC patients' subgroup compared to the reference group.
It is important to note that some experimental studies have indicated an opposite understanding of the significance of MMP-3 during the carcinogenesis process. In a series of studies in mice by McCawley et al., the authors have found that MMP-3 may play a protective role in the development of SCC (44, 45), as treatment in MMP-3-null mice using chemical carcinogens enhanced tumor growth as compared to the control mice. MMP-3-null mice demonstrated a faster initial tumor growth associated with increased cell proliferation and they had more undifferentiated or early metastatic tumors with more frequent distant lung metastases. From these results, the authors proposed that MMP-3 is expressed as a protective response and plays an important role in host defence during SCC tumorigenesis. The results of these studies could partly explain the more frequent presence of MMP-3 (-1171) 6A/6A genotype and allele 6A in advanced and moderate/poorly differentiated LSCC patients' group in our series. It may be presumed that a more frequent allele 6A may decrease the promoter transcription activity and MMP-3 gene expression, thus removing the possible protective role of MMP-3 during LSCC development.
In conclusion, the results of the present study demonstrate that combination of smoking and 6A/6A genotype of MMP-3 (-1171 5A/6A) polymorphism are associated with significantly increased odds of LSCC risk, thus confirming a presumably important role of MMP-3 (-1171) 6A/6A genotype in LSCC carcinogenesis. These results seem promising, as they may improve our understanding of the molecular events leading to LSCC formation and development. Further investigation in this field on a larger data sample is advocated.
Acknowledgements
The Authors thank Professor Abdonas Tamosiunas for the reference group (HAPPIE) data.
Footnotes
Authors' Contributions
All Authors (VL, VL, RL, PV, and VU) were involved in the conception and design of the present investigation. VL, RL, and PV performed laboratory tests. VL, RL, PV and VU analyzed the data. VL, RL and PV drafted the article. VU and VL revised the last version of manuscript. The final manuscript was read and approved by all authors of the study.
This article is freely accessible online.
Conflicts of Interest
The Authors have no competing interests to declare.
- Received September 22, 2019.
- Revision received October 10, 2019.
- Accepted October 18, 2019.
- Copyright© 2020, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved