Background: Chronic inflammation is an important pathological factor in colorectal tumorigenesis. Lymphotoxin-alpha (LT-α), a pleiotropic pro-inflammatory cytokine has been shown to possess both cancer promoting and cancer inhibiting activities. Several studies have analyzed the association of intronic LT-α+252A/G single-nucleotide polymorphism (SNP) in human LT-α gene with various cancers including colorectal cancer (CRC), but the outcome have been mixed and inconclusive. Materials and Methods: The present case-control study analyzed the association of LT-α+252A/G SNP with CRC risk in the ethnic Kashmiri population. The genotype frequencies of LT-α+252A/G intronic SNP were compared between 142 CRC patients and 184 individually matched healthy controls by polymerase chain reaction-restriction fragment length polymorphism method. The association between the LT-α+252A/G SNP and CRC risk was investigated through conditional logistic regression models adjusted for multiple possible confounding (third) variables. Further, the association between CRC risk and various clinico-pathological parameters, demographic variables, and environmental factors within the case group subjects and the SNP under study was also analyzed. Results: The overall association between the LT-α+252A/G SNP and the modulation of CRC risk was found to be significant (P = 0.013). Further, we found a significant effect modification of the association between the LT-α+252A/G SNP genotypes and CRC risk by gender (P = 0.046). We also found that the LT-α+252A/G SNP within the case group was significantly associated with gender (P = 0.0014) and lymph node status (P < 0.0001). Conclusion: This study has demonstrated that LT-α+252A/G SNP is significantly associated with risk of CRC in the ethnic Kashmiri population although the nature of this association could not be deciphered further in a statistically significant manner. However, this study needs to be replicated with larger sample size and if possible in other ethnically defined populations that exhibit comparable incidence of CRC to substantiate and elaborate our findings in a more comprehensive manner.
Introduction
Colorectal cancer (CRC) also known as colon cancer or large bowel cancer includes the neoplasia of the colon, rectum, and appendix and is one of the most common cancers and leading cause of cancer-related mortality and morbidity worldwide. It ranks third among the cancers in men and second in women globally, is the third-most frequent cause of cancer-related mortality with nearly 1.2 million new cases diagnosed each year and approximately 600,000 deaths reported annually.[1],[2] The Valley of Kashmir, part of Jammu and Kashmir State, located in Northern India has a unique ethnically defined population with highly distinctive food habits compared with the rest of the country and the world. In Kashmir Valley, CRC represents the third-most common gastrointestinal tract cancer after esophageal cancer and gastric cancer[3] and is the fourth-common cancer among the males and the third among the females.[3],[4]
The pathogenesis of CRC is strongly associated with innate immune processes and intestinal inflammation. The chronic inflammation plays a vital role in the promotion of colorectal tumorigenic processes. The chronic inflammation can also stimulate and complement the noninflammatory pro-tumorigenic processes through supply of various mediators including cytokines, which promote tumor growth, survival, progression, angiogenesis, invasion, and metastasis.[5],[6] Cytokines, both pro-inflammatory and anti-inflammatory, constitute the vital component in the mediation and regulation of inflammatory immune responses which play a key role in the pathogenesis of colorectal tumorigenesis.[7]
Human lymphotoxin-alpha (LT-α) is a pleiotropic pro-inflammatory cytokine, which plays an essential role in the regulation of various immune responses including inflammation, immuno-stimulation, antiviral responses, and cytotoxic activities on infected cells and tumor cells.[8],[9] LT-α is a vital mediator of lymphoid neogenesis and plays a key role in the development, organization, and maintenance of secondary and tertiary lymphoid organs including Peyer's patches, lymph nodes, the lymphoid follicles in tonsils, adenoids, and the spleen,[10],[11] which are crucial for the initiation and maintenance of adaptive immune responses.
The role of LT-α in colorectal and other tumors is controversial, as it has been shown to possess both cancer promoting and cancer inhibiting activities. LT-α is involved in anti-tumor surveillance through its cytotoxic effect on tumor cells. LT-α is involved in the activation of natural killer (NK) cells, which represent an important component of nonspecific host defense mechanism involved in tumor rejection and inhibition of metastasis. LT-α promotes the differentiation, maturation and recruitment of NK cells to tumor lesions.[12] However, various studies point towards the involvement of LT-α in the development of various cancers including CRC.[13]
The LT-α gene is located within the class III region of the major histocompatibility complex on the short arm of chromosome 6 at position 21.3. The gene encoding LT-α is highly polymorphic and has been reported to contain several different polymorphisms mostly in the form of single nucleotide polymorphisms (SNPs) which are found within the gene promoter besides the intronic and exonic regions. Of all the polymorphisms, the LT-α SNPs designated as rs1041981, rs2239704, rs2229094, rs746868, and rs909253 have been reported to be functionally relevant with regard to the regulation of gene expression and therefore may influence LT-α expression and contribute to the inter-individual variability in the expression and tissue/serum levels of this cytokine.[14],[15] This may further result in differences in relevant biological functions including the modulation of inflammatory response and subsequently the differences in susceptibility among different individuals to several diseases including various cancers.[16],[17],[18],[19],[20]
The LT-α+252A/G SNP (rs909253) also known as IVS1 +90 A/G SNP representing an adenine (A) to guanine (G) substitution at +252 nucleotide (nt) position is located within the first intron of LT-α gene and has been reported to be significantly associated with differential LT-α expression and activity. The AG and more potently GG genotypes and G allele of LT-α+252A/G SNP have been reported to be associated with increased expression, tissue/serum levels and higher activity of LT-α in comparison to AA genotype and A allele of this SNP.[21],[22],[23],[24],[25] The more common +252A allele is sometimes referred to as LT10.5 or LT-α (10.5 kb) whereas the less common +252G allele is referred to as LT5.5 or LT-α (5.5 kb).[26] It has been reported that the less common or variant genotype LT-α+252GG and less common or variant allele LT-α+252G resulted in a 1.5-fold increase in LT-α gene expression in comparison to the more common LT-α+252AA genotype and more common LT-α+252A allele and same has been reflected through the considerably higher tissue and serum levels of LT-α observed in subjects with GG genotype in comparison to the AA genotype.[23],[27],[28],[29] The LT-α+252A/G SNP also regulates the expression and activity of tumor necrosis factor-alpha (TNF-α) and LT-α+252G allele has also been reported to be associated with increased expression, enhanced tissue/serum levels and higher activity of TNF-α in comparison to LT-α+252A allele.[28],[29] Although the exact functional significance of LT-α+252A/G SNP is not known yet fully, it has been proposed that the differential LT-α expression and activity stems from the transcriptional regulatory activity of this SNP. The A to G substitution at +252 bp position has been reported to influence the binding of an unknown nuclear factor (NF) such that the protein binds with higher affinity to the +252G allele, resulting in 1.5-fold higher expression of LT-α protein in comparison to the expression associated with +252A allele.[23] The haplotype specific increased expression, tissue/serum levels, and higher activity of LT-α associated with LT-α+252A/G SNP may result in the constitutive activation of NF-κB and enhanced NF-κB signaling may modulate the inflammatory response and further stimulate or promote the activation of pathogenic pathways involved in several diseases including various types of cancers.
Several studies have analyzed the possible association of LT-α+252A/G SNP with risk of various cancer types including gastric cancer,[17],[30] breast cancer,[20] lung cancer,[31] non-Hodgkin's lymphoma,[32],[33] oral cancer,[34] endometrial cancer,[35] myeloma,[36] leukemia,[37] bladder cancer,[38] and cervical cancer.[39] However, we found only few studies that have analyzed the possible association of LT-α+252A/G SNP with CRC risk.[18],[40]
In this study, we systematically analyzed the possible association between LT-α+252A/G SNP and susceptibility to or risk of CRC in the Kashmiri population through a case-control study design. We also evaluated the possible effect modification of CRC risk by age, gender, and smoking status. Further, we investigated the possible association of LT-α+252A/G SNP with various clinico-pathological parameters, demographic variables and environmental factors including smoking habit and studied their role in modulating the risk of CRC in the population under study.
Materials and Methods
Study subjects.
The present study included two subject groups: Case and control. The case group included 142 individuals recruited consecutively irrespective of their age and gender with primary CRC who underwent surgical resection for primary CRC tumors at the Department of General Surgery, Sher-I-Kashmir Institute of Medical Sciences (SKIMS), Srinagar, Kashmir. The diagnosis of CRC was confirmed histopathologically. The tumor stage and the tumor grade were classified according to the 8th ed.ition of TNM classification of Union International Control of Cancer. Only those cases who had not received any neoadjuvant chemo or radiotherapy were chosen for this study. All the cases were >18 years old and had no prior history of any malignancy. Blood and tissue samples were obtained from these CRC patients. The control group included 184 healthy individuals with no history or prior diagnosis of any malignant disorder or any other serious disease, which were recruited during the same time period and from the same geographic area and from whom blood was collected and used as control for the present study. The control group included both general population-based subjects and hospital based subjects. The control group subjects were matched to the case group subjects individually for age (±5 years), sex, place of residence (rural/urban), smoking habit, and ethnicity to minimize the confounding effect of these various relevant factors. Both the case and the controls chosen for this study were ethnic Kashmiris.
Data collection
The data relevant to the study concerning all the CRC patients including various clinico-pathological parameters, demographic variables and the environmental factors was obtained and evaluated from the patient medical records, pathology reports and also from the personal interviews with the patients and/or their guardians (for those who were illiterate or unable to communicate). The interviews were conducted in local language for easy and direct communication, which also helped to gather maximum possible relevant information. The data collected included tumor location, Dukes Stage, lymph node status, age, sex, place of residence, ethnicity, smoking habit, and the family history of cancer among several other potential confounding parameters. The relevant data were also obtained for each of the recruited controls mostly through personal interviews and included parameters such as age, sex, place of residence, ethnicity, and smoking habit. The data collection was carried out by research professionals only, to fulfill the requisite quality standards during the course of this study. All the patients and/or their guardians were informed about the study and their willingness to participate in this study was documented using a predesigned questionnaire and same procedure was followed for the controls. All the procedures concerning the study participants including sample procurement and the data collection were carried out in accordance with the ethical standards laid down by the Institutional Ethics Committee, SKIMS and the World Health Organization and the Code of Ethics of the World Medical Association (Declaration of Helsinki, 1964 and its seventh amendment, 2013) for experiments in humans.[41]
Sample preparation and DNA extraction
The tumor tissue samples collected after surgical resection were immediately snap frozen in liquid nitrogen and then stored at −80°C until further use for DNA extraction and other experimental purposes. Peripheral blood sample, 3–5 ml from each case and control group individual was collected by venopuncture into ethylene diamine tetra acetic acid coated blood vacutainer collection tubes (purple capped tubes; ADS Hitech Polymers, India) and stored at −80°C until further use. The genomic DNA was extracted from both the tumor tissue and blood specimens using DNeasy™ Blood and Tissue Kit (catalog no. 69504; Qiagen, Germany) and Quick-gDNA™ MiniPrep kit (catalog no. D3024; Zymo Research, USA) according to the manufacturers' instructions. The extracted DNA was stored at −20°C until further use. The qualitative and the quantitative assessments of the extracted genomic DNA samples were carried out by absorbance measurements at 260 nm and 280 nm using ultraviolet-visible Spectrophotometric analysis and also by agarose gel electrophoresis. The DNA extracted from blood samples of case and control group subjects was used for this study.
Single nucleotide polymorphism analysis or genotyping
The LT-α+252A/G SNP was genotyped using the polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) assay.
Lymphotoxin-alpha +252A/G SNP polymerase chain reaction
The PCR for the amplification of LT-α gene region containing the LT-α+252A/G SNP was carried out in a total volume of 25 μL containing 100 ng – 1 &#s956;g of genomic DNA, 0.7–1 U Taq DNA polymerase with lX Standard Taq reaction buffer (New England Biolabs, UK), 1.8 mM MgCl2; 0.28 mM deoxynucleotide triphosphate mix (New England Biolabs, UK); 0.56 μM forward and revere oligonucleotide primers (Integrated DNA Technologies, India) and nuclease-protease free water (Qiagen, Germany) added up to a final volume of 25 μL. Alternatively and randomly Phusion DNA Polymerase with Phusion HF Buffer (New England Biolabs, Inc., UK) were used instead of Taq DNA polymerase with Standard Taq reaction buffer to check for any Taq polymerase induced amplification errors.
The PCR conditions used for the amplification of LT-α gene region containing the LT-α+252A/G SNP were as follows; initial denaturation at 95°C for 6 min followed by 35 cycles of denaturation at 95°C for 45 s; annealing at 64°C for 60 s and extension at 72°C for 45 s followed by a single final extension step at 72°C for 10 min. The oligonucleotide primers used for the amplification of the specific gene region containing the LT-α+252A/G SNP were 5'-CCGTGCTTCGTGCTTTGGACTA-3' (Forward) and 5'-AGAGCTGGTGGGGACATGTCTG-3' (Reverse). The desired PCR product obtained for LT-α+252A/G SNP was 741 bp in size [[Figure 1] for PCR gel picture].{Figure 1}
Genotyping
The LT-α+252A/G SNP was genotyped using the restriction enzyme NcoI. The restriction enzyme was procured from Thermo Fisher Scientific, USA. The digestion was carried out according to the manufacturers' instructions in a 30 μL reaction volume containing 10 μL of PCR product and 10 U of appropriate restriction enzyme and incubated at 37°C overnight. Wild genotype (AA) is not cleaved by NcoI enzyme thereby yielding a single 741 bp fragment whereas variant genotype (GG) yields two fragments 545 bp and 196 bp in size. Heterozygous genotype (AG) yields three fragments 741 bp, 545 bp, and 196 bp in size. The digestion products of LT-α+252A/G SNP were separated on a 3% agarose gel stained with ethidium bromide (HiMedia, India) to a final concentration of 0.5 μg/ml [[Figure 2] for RFLP digestion gel picture].{Figure 2}
Quality control
The quality control included the assessment of genotyping errors including the false estimates of a particular allele or genotype frequency and the evaluation of the reproducibility of the genotyping done. For these assessments, approximately 10% of the patient and control samples selected randomly were re-genotyped. In addition, in each PCR-RFLP setup, previously amplified and genotyped samples representing different genotypic scenarios were included as a positive control. The genotyping reproducibility of the samples was very high for the LT-α SNP with a weighted kappa coefficient of 0.99, which meant a high concordance rate of 99%.[42],[43]
Statistical analyses
The frequencies of genotypes and alleles for the SNP under study were obtained through direct counting. The numbers and percentages were calculated and presented for each of the categorical variables along with means, standard deviations, median, and inter-quartile range for continuous variables. Conditional logistic regression analysis was carried out to calculate unadjusted and adjusted odds ratios (ORs) and corresponding 95% confidence intervals (CIs) to assess the possible association of the relevant SNP genotypes with CRC risk and to assess the possible gene-environment interactions if applicable. In order to eliminate the possible confounding (third) variables, the conditional logistic regression models were adjusted for the known risk factors such as gender, age and smoking habit and with the place of residence. The possible effect measure modification of the association between various genotypes relevant to the SNP under study and CRC risk by various CRC risk factors including age, gender, and smoking status was also included in the conditional logistic regression models and analyzed. The correlation between the genotypes and the clinico-pathological parameters, demographic variables, and environmental factors including smoking habit within the case group was analyzed using Fisher's exact test. The fitness of the genotype distributions to Hardy-Weinberg equilibrium (HWE) for the allele and the genotype frequencies in the population under study was tested using the Chi-square test. A two sided probability value of or <5% (P ≤ 0.05) was considered statistically significant for all types of analyses. All statistical analyses were performed using IBM, Armonk, New York, United States.
The effective sample size and the statistical power were computed using the “Genetic Power Calculator” developed by Purcell et al. (http://zzz.bwh.harvard.edu/gpc/). We obtained a healthy power score of about 77% for the SNP under study in the present case-control study involving 142 case subjects and 184 controls.
Results
The frequencies of various clinico-pathological parameters, demographic variables, and environmental factors in CRC case subjects and relevant parameters in controls representing the general characteristics of study participants are given in [Table 1]. The frequencies of the genotypes of LT-α+252A/G SNP for both the case and the control groups are listed in [Table 2]. The more common AA genotype of LT-α+252A/G SNP was almost equal in frequency among the case group (59.86% [85/142]) and the control group (60.87% [112/184]). The frequency of the heterozygous genotype (AG) in the case group (40.14% [57/142]) was higher than that of the control group (33.7% [62/184]). The variant genotype (GG) was altogether absent in the case group but was present in the control group (5.43% [10/184]). Further, the frequency of the more common LT-α+252A allele was found to be 79.93% (227/284) among the case group subjects and 77.72% (286/368) among the control group subjects. The frequency of the less common LT-α+252G allele was found to be 20.07% (57/284) among the case group subjects and 22.28% (82/368) among the control group subjects. The frequency of the combined variant genotype (AG + GG) in the case group (40.14% [57/142]) was not much different from that of the control group (39.13% [72/184]) [Table 2]. The overall association between the LT-α+252A/G SNP and the modulation of CRC risk was found to be significant (P = 0.013) [Table 2]. However, the nature of this association could not be deciphered further in a statistically significant manner. Further, the genotype frequencies for LT-α+252A/G SNP were found to be in agreement with HWE among the control group subjects (χ2 = 0.135; P = 0.713) but not among the case group subjects (χ2 = 8.953; P = 0.002).{Table 1}{Table 2}
The possible effect measure modification or effect modification of the association between LT-α+252A/G genotypes and CRC risk by various CRC risk factors including age, gender, and smoking status are summarized in [Table 3]. On analyzing the effect modification of LT-α+252A/G genotypes by age, gender, and smoking status, it was found that the effect of the combined variant genotype (AG + GG) on CRC risk was significantly influenced by gender (P = 0.046). A decreased CRC risk was observed in females (OR, 0.42 [95% CI, 0.19–0.92]; P = 0.030).{Table 3}
The numbers and the frequencies of the subsets of various characteristics of the case group subjects under study, i.e., age, gender, dwelling, smoking status, tumor location, tumor grade, and lymph node status for LT-α+252A/G SNP are listed in [Table 4]. We analyzed the correlation of the LT-α+252A/G promoter SNP with the subsets of these various characteristics of the case group participants. The LT-α+252A/G SNP was significantly associated with gender (P = 0.0014). The male participants who carried the heterozygous genotype (AG) were at an increased risk of developing CRC in comparison to females (OR, 3.07 [95% CI, 1.52–6.19]; P = 0.0017). The LT-α+252A/G SNP was also significantly associated with smoking status (P = 0.0141). The participants who carried the heterozygous genotype (AG) and were smokers were at an increased risk of developing CRC (OR, 2.35 [95% CI, 1.18–4.66]; P = 0.0163). The LT-α+252A/G SNP also showed an overall strongly significant association with lymph node status (P < 0.0001). Further, the participants who carried the heterozygous genotype (AG) of LT-α+252A/G SNP were at an increased risk of lymph node infiltration (OR, 4.54 [95% CI, 2.21–9.31]); P < 0.0001. Some statistical parameters mentioned here are not shown in [Table 4].{Table 4}
Discussion
In the present study, we evaluated the role of functional LT-α+252A/G SNP located within the first intron of LT-α gene as a potential modulator of risk of CRC in ethnic Kashmiri population in a case-control study design with 142 case subects and 184 controls.
In the present study, we evaluated the distribution of the LT-α+252A/G SNP genotypes in CRC patients and controls, and found that this LT-α promoter SNP showed an overall significant association with the modulation of the CRC risk in our population. However, the nature of this association could not be deciphered further in a statistically significant manner. The factors that led to this inconclusive result could be both multiple and complex. The more expressive variant genotype LT-α+252GG was totally absent from the case group subjects which possibly led to the haploinsufficiency effect whereby the observed overall association of CRC risk with LT-α+252G allele was not strong enough to express itself clearly. In other words, the association was there but due to lack of LT-α+252GG genotype, the effect of LT-α+252G allele was less pronounced which possibly led to this inconclusively. At the protein level, this possibly means that the increased LT-α expression associated with the heterozygous genotype (AG) in case group subjects will be high compared to the low expressive wild genotype (AA) and this increased expression resulting from a single LT-α+252G allele though was high enough to exhibit that somehow the differential LT-α expression could modulate the CRC risk but not enough to depict this modulation clearly enough. Further studies involving a large number of CRC patients and healthy controls may possibly help to decipher this association in a statistically significant manner and may explain its influence on the CRC risk in a more conclusive manner and substantiate our hypothesis. Our results though inconclusive regarding the nature of risk modulation are still in agreement with one study that reported the association of the LT-α+252A/G SNP with modulation of CRC risk.[40]
In the present study, we also evaluated the possible effect modification of the association between LT-α+252A/G SNP genotypes and CRC risk by age, gender, and smoking status. We found a significant effect modification of association between the combined variant genotype (AG + GG) of LT-α+252A/G SNP and CRC risk by gender. The female gender decreased the CRC risk. The females carrying the variant allele, LT-α+252A/G in heterozygous form (AG) or variant homozygous form (GG) were at a decreased risk of developing CRC in comparison to males. A plausible explanation of this finding is that the regulation of immune response including the innate and adaptive immune responses and cytokine expression and activity in humans has been naturally programmed to exhibit a gender-based dimorphism.[44],[45],[46],[47] It is now well established that the immune response is regulated in a gender specific manner in which sex hormones, the androgens and the estrogens play an essential role but affect the immune system in opposite ways.[48],[49] The androgens have been reported to favor the T Helper 1 (Th1)-type response,[50] whereas as estrogens promote T Helper 2 (Th2)-type immune response.[48] It is also known now through various studies that the cytokines, interferon gamma, a T Helper 1 (Th1)-type cytokine, in males and interleukin-6 (IL-6), a T Helper 2 (Th2)-type cytokine, in females are the main players and regulators of immune system in a gender specific manner.[51],[52],[53] In other words, the T Helper 1 (Th1)-type immune response is predominant in males whereas T Helper 2 (Th2)-type immune response is predominant in females. LT-α is a T Helper 1 (Th1)-type cytokine and it is possible that LT-α may possibly be more active and may possibly have a more profound immunomodulatory role in males in comparison to females. This enhanced immunomodulatory role may in part be achieved through increased expression of LT-α in males in comparison to females. Further, the increased expression and circulating levels of LT-α have been associated with the development of various cancers including CRC.[13],[54],[55],[56] Therefore, with respect to our findings, it is reasonable to argue that the affect of LT-α+252G allele resulting in the increased expression, serum levels and activity of LT-α will be more profound in males in comparison to the females considering the predominance of T Helper 1 (Th1)-type response in males. This may explain our finding that the females carrying LT-α+252G allele were at a decreased risk of developing CRC in comparison to males. However, this hypothesis needs to be mechanistically validated to obtain conclusive evidence. Further studies involving comprehensive mechanistic evaluation of the possible gender specific differences in the expression and circulating levels of LT-α protein are warranted to arrive at a conclusive explanation, which may also possibly substantiate our hypothesis. Further, it is important to emphasize here that the immune response may be regulated differently in a gender-specific manner but the outcome, that is, the preservation of immune homeostasis is similar in both males and females.
In this study, we also evaluated the association of LT-α+252A/G SNP with the numbers and the frequencies of the subsets of various characteristics of the case group subjects under study, i.e., age, gender, dwelling, smoking status, tumor location, tumor grade, and lymph node status.
When stratifying by gender, we found that the male participants who carried the heterozygous genotype (AG) of LT-α+252A/G SNP were at an increased risk of developing CRC in comparison to females. This finding is similar to the one observed in case of effect modification analysis and as such has already been explained.
When stratifying by smoking status, we found that the subjects, who carried the heterozygous genotype (AG) of LT-α+252A/G SNP and were smokers, had an increased risk of developing CRC in comparison to non-smokers. Smoking is an established risk factor of CRC.[51],[52],[53] Smoking has been reported to induce inflammatory response[60] by stimulating and increasing the production of various pro-inflammatory cytokines such as TNF-α, LT-α IL-1, IL-6, IL-8 and GM-CSF and by inhibiting and/or decreasing the production of anti-inflammatory cytokines such as IL-10.[60],[61],[62],[63] Smoking has also been reported to decrease NK cell function and impair regulatory T cell (Treg) activity.[63] In a broader sense, the smokers in comparison to nonsmokers are exposed to an increased burden of the carcinogenic agents from cigarette/tobacco smoking which increases the risk of developing CRC in comparison to nonsmokers. Further, the heterozygous genotype (AG) is associated with moderate increase in LT-α production in comparison to low producer genotype (AA). Also increased expression and circulating levels of LT-α have been associated with the development of various cancers including CRC.[13],[54],[55],[56] Thus, in response to our finding, it is reasonable to argue that for the same genotype LT-α+252AG and the level of other environmental exposures, smokers compared to nonsmokers were at an increased risk of developing CRC.
When stratifying by lymph node status, we found that the subjects who carried the heterozygous genotype (AG) of LT-α+252A/G SNP were at an increased risk of lymph node infiltration. As discussed already, the heterozygous genotype (AG) is associated with a moderate increase in LT-α production in comparison to low producer genotype (AA).Further, the increased expression and circulating levels of LT-α have been associated with the development of various cancers including CRC.[13],[54],[55],[56] and more severe outcome in some cancers,[20],[64] which manifests itself through an increased tumor invasion and metastasis and is partly reflected through increased lymph node infiltration. Therefore, with reference to our finding that the carriers of the heterozygous genotype (AG) of LT-α+252A/G SNP were at an increased risk of lymph node infiltration, it is reasonable to argue that the increased LT-α expression, tissue/serum levels and eventually the increased LT-α activity associated with the LT-α+252AG genotype may be responsible for the increased colorectal tumor invasion and metastasis which encompasses increased lymph node infiltration. However, further studies involving comprehensive mechanistic evaluation of the expression and the circulating levels of LT-α protein in various phases of colorectal tumorigenesis and their possible role in tumor progression are warranted to prove our hypothesis.
In this study, we found that the genotype frequencies of LT-α+252A/G SNP among the case group subjects deviated significantly from the HWE. The main reasons responsible for the deviations from the HWE include inbreeding, consanguinity, population stratification, small population size, migration, mutations, and genotyping errors[65],[66],[67] and some of these are possibly responsible for the deviations from HWE observed in our study. The Kashmiri population consisting mostly Muslims represents an almost pure ethnic population[68] among whom consanguineous marriages are quite common and often traditional.[69] Further, due to overall genetic isolation from the rest of the world, the study population also exhibits a considerable degree of stratification and the population is also relatively small in size.[68] The genotyping errors resulting from technological and human reasons or both also represent an important reason behind the deviation from HWE in association studies but the genotyping errors are unlikely to be responsible for the deviations from HWE observed in the present study because we consistently obtained a high degree of the genotyping reproducibility for the samples studied. We obtained a weighted kappa coefficient of 0.99, which corresponded to a high concordance rate of 99% for the SNP studied.
To the best of our knowledge, the present study is the first of its kind with regard to examining the association of LT-α+252A/G SNP in the Kashmiri population. The use of clinically diagnosed and histopathologically confirmed CRC samples, involvement of population-based controls in addition to the hospital-based controls which were recruited from the same geographical area during the same time period and were matched to the case group subjects individually for age, sex, place of residence (rural/urban), smoking habit and ethnicity to minimize the confounding effect of these relevant factors and the adjustment of results for multiple potential confounding (third) variables can be described as the major strengths of this study. The major limitation of this study is the modest sample size to detect comprehensively the gene-gene and gene-environment interactions which usually require much larger sample size. Therefore, more of the similar studies but with a larger sample size and if possible incorporating other ethnic populations are needed to substantiate our findings or elaborate our findings in a more comprehensive manner regarding the association of the SNP under study with risk of CRC. However, these limitations are unlikely to affect the outcome of this study.
Conclusion
We have demonstrated through this study that the LT-α+252A/G intronic SNP is significantly associated with risk of CRC in the ethnic Kashmiri population. However, the nature of this association could not be deciphered further in a statistically significant manner. Further studies involving a large number of CRC patients and healthy controls may possibly help to decipher this association in a statistically significant manner and may explain its influence on the CRC risk in a more conclusive manner. We have also demonstrated that there is a significant effect modification of the association between LT-α+252A/G SNP genotypes and CRC risk by gender. Further, we have also demonstrated that there is a significant association between LT-α+252A/G SNP and some characteristics of the case group subjects including gender, smoking status and lymph node status.
Acknowledgments
MZB acknowledges the Department of Biotechnology, University of Kashmir for providing all the research facilities and University Grants Commission (UGC), GOI for UGC-BSR fellowship. MZB would like to thank his colleagues Taseen Gul, Insha Zahoor, Zeenat Hamid, Henah Mehraj Balkhi and Ashaq Hussain Mir from Department of Biotechnology, University of Kashmir, and Idrees Ayoub Shah from the Department of Biochemistry, University of Kashmir, Hazratbal, Srinagar.
The authors wish to thank each and every CRC patient and control subject who took part in this study and cooperated during the interview and sample collection. We also thank the head and technical staff of the operating theatre in the Department of General Surgery, SKIMS, Soura, Kashmir who helped us with tissue procurement, and the anonymous pathologists at the Department of Pathology, SKIMS, Soura, Kashmir for the histopathological assessment of the tumor tissues.
Financial support and sponsorship
MZB acknowledges University Grants Commission (UGC), GOI for UGC-BSR fellowship.
Conflicts of interest
There are no conflicts of interest.
MZB acknowledges University Grants Commission (UGC), GOI for UGC-BSR fellowship.
Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. Global cancer statistics. CA Cancer J Clin 2011;61:69-90.
Torre LA, Siegel RL, Ward EM, Jemal A. Global cancer incidence and mortality rates and trends – An update. Cancer Epidemiol Biomarkers Prev 2016;25:16-27.
Sameer AS. Colorectal cancer: Molecular mutations and polymorphisms. Front Oncol 2013;3:114.
Wani MA, Jan FA, Khan NA, Pandita KK, Khurshid R, Khan SH. Cancer trends in Kashmir; common types, site incidence and demographic profiles: National Cancer Registry 2000-2012. Indian J Cancer 2014;51:133-7.
Qian BZ, Pollard JW. Macrophage diversity enhances tumor progression and metastasis. Cell 2010;141:39-51.
Wogan GN, Dedon PC, Tannenbaum SR, Fox JG. Infection, inflammation and colon carcinogenesis. Oncotarget 2012;3:737-8.
Dranoff G. Cytokines in cancer pathogenesis and cancer therapy. Nat Rev Cancer 2004;4:11-22.
Aggarwal BB. Signalling pathways of the TNF superfamily: A double-edged sword. Nat Rev Immunol 2003;3:745-56.
Locksley RM, Killeen N, Lenardo MJ. The TNF and TNF receptor superfamilies: Integrating mammalian biology. Cell 2001;104:487-501.
Drayton DL, Liao S, Mounzer RH, Ruddle NH. Lymphoid organ development: From ontogeny to neogenesis. Nat Immunol 2006;7:344-53.
Ware CF. Network communications: Lymphotoxins, LIGHT, and TNF. Annu Rev Immunol 2005;23:787-819.
Smyth MJ, Johnstone RW, Cretney E, Haynes NM, Sedgwick JD, Korner H, et al. Multiple deficiencies underlie NK cell inactivity in lymphotoxin-alpha gene-targeted mice. J Immunol 1999;163:1350-3.
Popivanova BK, Kitamura K, Wu Y, Kondo T, Kagaya T, Kaneko S, et al. Blocking TNF-alpha in mice reduces colorectal carcinogenesis associated with chronic colitis. J Clin Invest 2008;118:560-70.
Boraska V, Rayner NW, Groves CJ, Frayling TM, Diakite M, Rockett KA, et al. Large-scale association analysis of TNF/LTA gene region polymorphisms in type 2 diabetes. BMC Med Genet 2010;11:69.
Tan JH, Temple SE, Kee C, Waterer GW, Tan CR, Gut I, et al. Characterisation of TNF block haplotypes affecting the production of TNF and LTA. Tissue Antigens 2011;77:100-6.
Chae YS, Kim JG, Sohn SK, Moon JH, Kim SN, Lee SJ, et al. Lymphotoxin alfa and receptor-interacting protein kinase 1 gene polymorphisms may correlate with prognosis in patients with diffuse large B cell lymphoma treated with R-CHOP. Cancer Chemother Pharmacol 2010;65:571-7.
Lu R, Dou X, Gao X, Zhang J, Ni J, Guo L. A functional polymorphism of lymphotoxin-alpha (LTA) gene rs909253 is associated with gastric cancer risk in an Asian population. Cancer Epidemiol 2012;36:e380-6.
Sainz J, Rudolph A, Hoffmeister M, Frank B, Brenner H, Chang-Claude J, et al. Effect of type 2 diabetes predisposing genetic variants on colorectal cancer risk. J Clin Endocrinol Metab 2012;97:E845-51.
Zhang Y, Wang MY, He J, Wang JC, Yang YJ, Jin L, et al. Tumor necrosis factor-alpha induced protein 8 polymorphism and risk of non-Hodgkin's lymphoma in a Chinese population: A case-control study. PloS One 2012;7:e37846.
Zhou P, Huang W, Chu X, Du LF, Li JP, Zhang C. The lymphotoxin-α 252A>G polymorphism and breast cancer: A meta-analysis. Asian Pac J Cancer Prev 2012;13:1949-52.
Clarke R, Xu P, Bennett D, Lewington S, Zondervan K, Parish S, et al. Lymphotoxin-alpha gene and risk of myocardial infarction in 6,928 cases and 2,712 controls in the ISIS case-control study. PLoS Genet 2006;2:e107.
dos Santos M, Stur E, Maia LL, Agostini LP, Peterle GT, Mendes SO, et al. Genetic variability of inflammatory genes in the Brazilian population. Genet Test Mol Biomarkers 2013;17:844-8.
Hansson GK, Robertson AK, Söderberg-Nauclér C. Inflammation and atherosclerosis. Annu Rev Pathol 2006;1:297-329.
Temple SE, Cheong KY, Almeida CM, Price P, Waterer GW. Polymorphisms in lymphotoxin alpha and CD14 genes influence TNFalpha production induced by Gram-positive and Gram-negative bacteria. Genes Immun 2003;4:283-8.
Wang Q. Molecular genetics of coronary artery disease. Curr Opin Cardiol 2005;20:182-8.
Messer G, Spengler U, Jung MC, Honold G, Blömer K, Pape GR, et al. Polymorphic structure of the tumor necrosis factor (TNF) locus: An NcoI polymorphism in the first intron of the human TNF-beta gene correlates with a variant amino acid in position 26 and a reduced level of TNF-beta production. J Exp Med 1991;173:209-19.
Temple SE, Almeida CM, Cheong KY, Wunderink RG, Waterer GW. A diplotype in the lymphotoxin alpha gene is associated with differential expression of LTA mRNA induced by Gram-positive and Gram-negative bacteria. Int J Immunogenet 2007;34:157-60.
Pociot F, Briant L, Jongeneel CV, Mölvig J, Worsaae H, Abbal M, et al. Association of tumor necrosis factor (TNF) and class II major histocompatibility complex alleles with the secretion of TNF-alpha and TNF-beta by human mononuclear cells: A possible link to insulin-dependent diabetes mellitus. Eur J Immunol 1993;23:224-31.
Stüber F, Petersen M, Bokelmann F, Schade U. A genomic polymorphism within the tumor necrosis factor locus influences plasma tumor necrosis factor-alpha concentrations and outcome of patients with severe sepsis. Crit Care Med 1996;24:381-4.
García-González MA, Nicolás-Pérez D, Lanas A, Bujanda L, Carrera P, Benito R, et al. Prognostic role of host cyclooxygenase and cytokine genotypes in a Caucasian cohort of patients with gastric adenocarcinoma. PLoS One 2012;7:e46179.
Rausch SM, Gonzalez BD, Clark MM, Patten C, Felten S, Liu H, et al. SNPs in PTGS2 and LTA predict pain and quality of life in long term lung cancer survivors. Lung Cancer 2012;77:217-23.
Aissani B, Ogwaro KM, Shrestha S, Tang J, Breen EC, Wong HL, et al. The major histocompatibility complex conserved extended haplotype 8.1 in AIDS-related non-Hodgkin lymphoma. J Acquir Immune Defic Syndr 2009;52:170-9.
Skibola CF, Bracci PM, Nieters A, Brooks-Wilson A, de Sanjosé S, Hughes AM, et al. Tumor necrosis factor (TNF) and lymphotoxin-alpha (LTA) polymorphisms and risk of non-Hodgkin lymphoma in the InterLymph Consortium. Am J Epidemiol 2010;171:267-76.
Yapijakis C, Serefoglou Z, Vylliotis A, Nkenke E, Derka S, Vassiliou S, et al. Association of polymorphisms in Tumor Necrosis Factor Alpha and Beta genes with increased risk for oral cancer. Anticancer Res 2009;29:2379-86.
Niwa Y, Ito H, Matsuo K, Hirose K, Ito N, Mizuno M, et al. Lymphotoxin-alpha polymorphisms and the risk of endometrial cancer in Japanese subjects. Gynecol Oncol 2007;104:586-90.
Davies FE, Rollinson SJ, Rawstron AC, Roman E, Richards S, Drayson M, et al. High-producer haplotypes of tumor necrosis factor alpha and lymphotoxin alpha are associated with an increased risk of myeloma and have an improved progression-free survival after treatment. J Clin Oncol 2000;18:2843-51.
Demeter J, Porzsolt F, Rämisch S, Schmidt D, Schmid M, Messer G. Polymorphism of the tumour necrosis factor-alpha and lymphotoxin-alpha genes in chronic lymphocytic leukaemia. Br J Haematol 1997;97:107-12.
Nonomura N, Tokizane T, Nakayama M, Inoue H, Nishimura K, Muramatsu M, et al. Possible correlation between polymorphism in the tumor necrosis factor-beta gene and the clinicopathological features of bladder cancer in Japanese patients. Int J Urol 2006;13:971-6.
Nieves-Ramirez ME, Partida-Rodriguez O, Alegre-Crespo PE, Tapia-Lugo Mdel C, Perez-Rodriguez ME. Characterization of single-nucleotide polymorphisms in the tumor necrosis factor α promoter region and in lymphotoxin α in squamous intraepithelial lesions, precursors of cervical cancer. Transl Oncol 2011;4:336-44.
de Jong MM, Nolte IM, te Meerman GJ, van der Graaf WT, de Vries EG, Sijmons RH, et al. Low-penetrance genes and their involvement in colorectal cancer susceptibility. Cancer Epidemiol Biomarkers Prev 2002;11:1332-52.
World Medical Association. World Medical Association Declaration of Helsinki: Ethical principles for medical research involving human subjects. JAMA 2013;310:2191-4.
Bonin A, Bellemain E, Bronken Eidesen P, Pompanon F, Brochmann C, Taberlet P. How to track and assess genotyping errors in population genetics studies. Mol Ecol 2004;13:3261-73.
Wang DG, Fan JB, Siao CJ, Berno A, Young P, Sapolsky R, et al. Large-scale identification, mapping, and genotyping of single-nucleotide polymorphisms in the human genome. Science 1998;280:1077-82.
Aulock SV, Deininger S, Draing C, Gueinzius K, Dehus O, Hermann C. Gender difference in cytokine secretion on immune stimulation with LPS and LTA. J Interferon Cytokine Res 2006;26:887-92.
Cannon JG, St Pierre BA. Gender differences in host defense mechanisms. J Psychiatr Res 1997;31:99-113.
Moxley G, Posthuma D, Carlson P, Estrada E, Han J, Benson LL, et al. Sexual dimorphism in innate immunity. Arthritis Rheum 2002;46:250-8.
Schuurs AH, Verheul HA. Effects of gender and sex steroids on the immune response. J Steroid Biochem 1990;35:157-72.
McCarthy M. The “gender gap” in autoimmune disease. Lancet 2000;356:1088.
Whitacre CC. Sex differences in autoimmune disease. Nat Immunol 2001;2:777-80.
Beagley KW, Gockel CM. Regulation of innate and adaptive immunity by the female sex hormones oestradiol and progesterone. FEMS Immunol Med Microbiol 2003;38:13-22.
Berghella AM, Contasta I, Del Beato T, Pellegrini P. The discovery of how gender influences age immunological mechanisms in health and disease, and the identification of ageing gender-specific biomarkers, could lead to specifically tailored treatment and ultimately improve therapeutic success rates. Immun Ageing 2012;9:24.
Contasta I, Totaro R, Pellegrini P, Del Beato T, Carolei A, Berghella AM. A gender-related action of IFNbeta-therapy was found in multiple sclerosis. J Transl Med 2012;10:223.
Pellegrini P, Contasta I, Del Beato T, Ciccone F, Berghella AM. Gender-specific cytokine pathways, targets, and biomarkers for the switch from health to adenoma and colorectal cancer. Clin Dev Immunol 2011;2011:819724.
Haybaeck J, Zeller N, Wolf MJ, Weber A, Wagner U, Kurrer MO, et al. A lymphotoxin-driven pathway to hepatocellular carcinoma. Cancer Cell 2009;16:295-308.
Hehlgans T, Stoelcker B, Stopfer P, Müller P, Cernaianu G, Guba M, et al. Lymphotoxin-beta receptor immune interaction promotes tumor growth by inducing angiogenesis. Cancer Res 2002;62:4034-40.
Or YY, Chung GT, To KF, Chow C, Choy KW, Tong CY, et al. Identification of a novel 12p13.3 amplicon in nasopharyngeal carcinoma. J Pathol 2010;220:97-107.
Botteri E, Iodice S, Bagnardi V, Raimondi S, Lowenfels AB, Maisonneuve P. Smoking and colorectal cancer: A meta-analysis. JAMA 2008;300:2765-78.
Hannan LM, Jacobs EJ, Thun MJ. The association between cigarette smoking and risk of colorectal cancer in a large prospective cohort from the United States. Cancer Epidemiol Biomarkers Prev 2009;18:3362-7.
Otani T, Iwasaki M, Yamamoto S, Sobue T, Hanaoka T, Inoue M, et al. Alcohol consumption, smoking, and subsequent risk of colorectal cancer in middle-aged and elderly Japanese men and women: Japan Public Health Center-based prospective study. Cancer Epidemiol Biomarkers Prev 2003;12:1492-500.
Yanbaeva DG, Dentener MA, Creutzberg EC, Wesseling G, Wouters EF. Systemic effects of smoking. Chest 2007;131:1557-66.
Arnson Y, Shoenfeld Y, Amital H. Effects of tobacco smoke on immunity, inflammation and autoimmunity. J Autoimmun 2010;34:J258-65.
Sopori M. Effects of cigarette smoke on the immune system. Nat Rev Immunol 2002;2:372-7.
Vassallo R, Tamada K, Lau JS, Kroening PR, Chen L. Cigarette smoke extract suppresses human dendritic cell function leading to preferential induction of Th-2 priming. J Immunol 2005;175:2684-91.
Seidemann K, Zimmermann M, Book M, Meyer U, Burkhardt B, Welte K, et al. Tumor necrosis factor and lymphotoxin alfa genetic polymorphisms and outcome in pediatric patients with non-Hodgkin's lymphoma: Results from Berlin-Frankfurt-Münster Trial NHL-BFM 95. J Clin Oncol 2005;23:8414-21.
Hosking L, Lumsden S, Lewis K, Yeo A, McCarthy L, Bansal A, et al. Detection of genotyping errors by Hardy-Weinberg equilibrium testing. Eur J Hum Genet 2004;12:395-9.
Leal SM. Detection of genotyping errors and pseudo-SNPs via deviations from Hardy-Weinberg equilibrium. Genet Epidemiol 2005;29:204-14.
Wigginton JE, Cutler DJ, Abecasis GR. A note on exact tests of Hardy-Weinberg equilibrium. Am J Hum Genet 2005;76:887-93.
Ayub SG, Ayub T, Khan SN, Rasool S, Mahboob-ul-Hussain , Wani KA, et al. Epidemiological distribution and incidence of different cancers in Kashmir valley 2002-2006. Asian Pac J Cancer Prev 2011;12:1867-72.
Fareed M, Afzal M. Estimating the inbreeding depression on cognitive behavior: A population based study of child cohort. PLoS One 2014;9:e109585.
|