Skip to main content

MassARRAY analysis of twelve cancer related SNPs in esophageal squamous cell carcinoma in J&K, India



MassARRAY (Agena Bioscience™) combines competitive PCR with MALDI-TOF mass spectrometry (MS) analysis that gives highly accurate, sensitive, and high-throughput methods for the quantitative analysis of variation of gene expression in multiple samples. SNPs (Single Nucleotide Polymorphisms) have a very high potential of discovering disease-gene relationships. SNP-genotyping through MassARRAY is not only a cost-effective genotyping method but also provides a platform to validate variants observed through a high-throughput Next-generation sequencing (NGS).


In the present study, we have incorporated the use of matrix-assisted laser desorption/ionization-time of flight, mass spectrometry (MALDI-TOF) as a tool for differentiating genotypes based on the mass of variant. We have performed multiplex PCR and genotyped 12 SNPs in 758 samples (166 cases and 592 controls). The 12 studied SNPs were chosen with a rationale for their association with multiple cancers in literature.


This is the first study to explore these SNPs with esophageal cancer within the J&K population. Out of 12 SNPs, two SNPs rs12190287 of TCF21 and rs10046 of CYP19A1 were significantly associated with esophageal cancer with Odds Ratio (OR) 1.412 (1.09–1.8 at 95% CI, p = 0.008) and 1.54 (1.21–2.072 at 95% CI, p = 0.0007) within the population of Jammu and Kashmir.


We explored 12 SNPs that were found to be associated with multiple cancers in literature with esophageal cancer within the population of J&K. This is the first study to find the relation of these SNPs with ESCC within the studied population. This study explores the relation of genetic and environmental factors with the ESCC susceptibility.

Peer Review reports


Esophageal cancer (EC) is the most common type of cancer worldwide but is least studied with poor survival and highly aggressive nature. Various risk factors have been associated with EC which include lifestyle, dietary habits, low socio-economic status, poor oral hygiene, and genetics [1]. According to a survey by GLOBOCAN in 2018 about 572,034 of new cases and 508,585 deaths were reported worldwide about EC. About 70% of cases have been observed in men, and there is a 2–3 fold difference in incidence and mortality rates between the sexes worldwide [2]. In India, the Northeastern states like Assam, Meghalaya, Mizoram, and Nagaland tops the chart about EC, for both men and women. Kashmir valley is another state with the highest incidence of EC [3]. Genetics is one of the major risk factors associated with EC which includes loss and gain of chromosomes, gene amplification, and microsatellite instability [4]. 90% of esophageal cancers are ESCC (esophageal squamous cell carcinoma) and about 5% are EAC (esophageal adenocarcinoma) throughout the world. The remaining 5% include rare malignancies [4]. We attempted to investigate the role of cancer-related genetic variants in ESCC within the population of J&K. J&K is an ignored state owing to its peculiar geographical background and political circumstances. Although the incidence rate of EC is very high, yet the data related to the genetics of ESCC in the studied population is meager. Replication studies are important to enhance the credibility of a study. Genes that have been explored about ESCC within the population of J&K are summarized in the (Supplementary Table 1).



A total of 758 individuals (166 cancerous and 592 healthy controls) were recruited for the study after approval from the Institutional Ethical Review Board (IERB) - Shri Mata Vaishno Devi University (SMVDU) under notification number (SMVDU/IERB/16/41). All the patients recruited were not on chemotherapy, radiotherapy, and do not have any other form of cancer. All the control samples recruited were without any family history of ESCC. A 2 ml of a venous blood sample was collected. A well-written consent was taken from both cases and controls (Supplementary file 1). The clinical parameters of both cases and controls are provided in Supplementary Table 2.

DNA isolation

The genomic DNA was isolated by using the phenol-chloroform method [5] and the Qiagen DNA isolation Kit (Catalogue No. 51206).

The quality of the genomic DNA was checked by agarose gel electrophoresis (Bio-Rad Gel Doc™ EZ imager) and quantification was performed using Eppendrof’s BioSpectrometer™ at wavelength 260 nm and 280 nm. The amount of DNA was calculated using the following formula, DNA μg/ml = OD at 260 nm x 50X Dilution factor. The ratio at 260/280 was taken as criteria to check the purity of DNA.


Agena MassARRAY platform was used for SNP genotyping. It is a robust tool with high accuracy and is cost-effective as it involves multiplex PCR [6]. Genotyping was performed in the Central MassARRAY Analyzer facility at Shri Mata Vaishno Devi University. Agena Design Suite v2.0 was used to design-forward, reverse, and single base extension primer (customized) (Supplementary Table 3). Multiplex PCR was used to amplify the targeted region. 1 μl volume of genomic DNA (10 ng/μl concentration) was loaded in 384 well PCR plates which were dried at 85 °C for 10 min. After drying, the reaction mixture was prepared. The reaction mixture includes dNTPs, primer pool (forward and reverse primers pool), reaction buffers, and DNA polymerase. After the first PCR is done, the reaction mixture was treated with shrimp alkaline phosphatase (SAP), the multiplex PCR reaction is then performed with modified ddNTPs and primers (pooled single extension primers). PCR conditions were used using the Gabriel 2009 protocol [6].

After the final PCR, the product was subjected to resin treatment and then is transferred to 384 well Spectro-CHIP using MassARRAY Nanodispenser. The product is then subjected to MT analyzer and the data obtained is analyzed by pre-installed Typer Analyzer v4.0. The genotyping results were replicated in 10% of random samples and the concordance rate was 98.6%. In the reaction of 384 well plates, one positive and one negative control was added for quality check.

Genotyping quality control

SNPs with a call rate above 90% were only included for further the statistical analysis. All the 12 SNPs were following HWE (Hardy Weinberg Equilibrium) in both cases and controls.

Statistical analysis

The statistical analysis was done using Plink v1.09 [7] with a maximum of 1000 permutations. All the 12 SNPs were following HWE (Hardy Weinberg Equilibrium). The significance level was calculated based on 3 × 2 chi-square tests for genotype frequencies between ESCC cases and healthy controls. Further logistic regression analysis was also done using SPSS V.23. The analysis was done based on Odds Ratio (OR), confidence interval (CI) and p-value as the level of significance from confounding factors like age, gender, and BMI.

GTEx portal was used to determine NES (normalized effect size) value {Low, 2017 #12} and the gene expression of associated SNPs and variant effect prediction (Supplementary Figures 1 and 2). NES value for variants rs10046 is negative and rs12190287 was positive. Both are significantly associated with the esophageal tissue. A positive beta means that the tested allele is associated with an increased expression of that gene; a negative beta means that the tested allele is associated with a reduction in gene expression.


This is a case-control study that includes 166 ESCC cases and 592 healthy controls belonging to the Northern region of India. The mean age and standard deviation (SD) for cases were 60.4 ± 12.6 and for controls it was 58.4 ± 18.4. A total of 86 males and 80 females (cases) and 192 males and 400 females (controls) were recruited in the study. The BMI which is one of the main risk factors for malignancies was also recorded. The BMI of patients was 21.1 ± 5 and in the case of controls it was 27.6 ± 5.1. In about 68 cases there was metastasis while in 98 cases were non-metastatic. While recording family history, it was observed that about 22 cases had a family history of ESCC, 46 cases did not have a family history and 98 cases were not aware of the history of ESCC in their families. In the present study, two SNPs were showing an increased risk of ESCC within the population. These SNPs are rs12190287 of TCF21 and rs10046 of CYP19A1. In our study, the genetic variant rs12190287 has been evaluated concerning ESCC and it was observed that the variant under study was associated with the higher risk of ESCC within the population of Jammu and Kashmir with OR 1.412 (1.09–1.8, at 95% CI, p = 0.008). The genetic variant in CYP19A1 (Cytochrome P450 family 19 sub-familyA1) rs10046 (C > T) has been associated with a higher risk of ESCC with OR 1.584 (1.21–2.072, at 95% CI, p = 0.007) as shown in Table 1.

Table 1 Association analyses of the variants with Esophageal Cancer

In the present study, all the genetic variants were tested individually in a standard way to find its association with the disease. Both significant and non-significant SNPs are equally important for SNP analysis in the case-control study design. The SNP is said to be significantly associated or not associated with the risk of disease based on the p-value. If the p-value is less than or equal to a specific threshold (0.05), the SNP is said to be significantly associated with the higher risk of the disease provided OR (odd’s ratio) is above 1 at 95% CI (confidence interval) and if the p-value is greater than 0.05 then the variant is considered to be not associated [8]. The details of both significant and non-significant SNPs have been given in Table 1. All the SNPs summarized in the table below were following HWE. In any genetic association study, HWE is an essential tool to find genotyping errors [9].


Esophageal cancer is a very belligerent type of carcinoma with a very low survival rate. It’s because of its truculent nature, it is least studied worldwide. There are numerous risk factors associated with esophageal carcinoma which include genetic factors and environmental factors [10]. In Asia, ESCC is the most prevalent type of cancer while as in Caucasians there is a higher incidence of EAC. It is because of the different heritage backgrounds that these two populations show different susceptibilities towards esophageal carcinoma [11]. Epidemiological studies have shown that esophageal carcinoma is highly prevalent in the state of Assam located in the eastern region of India [12]. Jammu and Kashmir are one of the high incidence areas of “Central Asian esophageal cancer belt”. Central Asian esophageal cancer belt includes northern China, republics of Kazakhstan, Uzbekistan, and Turkmenistan. Kashmir Valley also falls under this high incidence region at the southernmost end [4]. The State of Jammu and Kashmir is geographically divided into three distinct regions Jammu, Kashmir, and Ladakh. They are culturally and ethnically distinct. These regions have historically/geographically married within their communities thus preserving the genetic pool. Though various esophageal carcinoma studies from the Kashmir region have explored demography and genetics but it remains uninvestigated from the Jammu region. There is a major difference between the dietary habits of people from Jammu as compared to people from the Kashmir or the Ladakh region [13]. A detailed questionnaire was designed for the patients in the present study which consisted of information like family history, age, gender, BMI, etc.

In the present study, genetic elucidation among cases and controls was explored. Genetics is an important risk factor associated with esophageal carcinoma. In the SNP analysis we have replicated the newly identified variants (rs12190287 of TCF21, rs10046 of CYP19A1, rs2735940 of TERT, rs751402 of ERCC5, rs2699887 of SLC14A2, rs3792152 of REV1, rs10069690 of TERT, rs2981582 of FGFR2, rs1695 of GSTP1, rs251796 of TERF2, rs2229080 of DCC, rs1801010 of BCL2. Two SNPs that were showing an increased risk of ESCC within the population of our study are rs12190287 of TCF21 and rs10046 of CYP19A1.

TCF21 (Transcription factor 21) helps in epithelial differentiation and thus plays a specific role in the differentiation of one or more subsets of epicardial cell types {Copyright© 1996–2018#97}. It is a candidate tumor suppressor located at chromosome 6 and it has been associated with lung, head, and neck cancers. It was found that TCF21 in gastric carcinoma cells was highly methylated than in the normal adjacent cells. rs12190287 (C > G) of TCF21 has been studied with Osteosarcoma risk in the Chinese population and was showing a higher risk of Osteosarcoma in the studied population [14].

CYP19A1 (Cytochrome P450 family 19 sub-familyA1) gene plays an important role in the biosynthesis of estrogen and it has been associated with the progression of breast cancer in the Chinese population [15]. Several association studies have been done with CYP19A1 and it was observed that there are enormous inter-population differences. The genetic variant rs10046 (C > T) has been significantly associated with Iranian women, Inuit women, and Xinjiang Uygur women [16].

Both the genetic variants rs12190287 of TCF21 and rs10046 of CYP19A1 have been studied in multiple cancers which include breast, lung, neck, and gastric cancer. Our study associated the above genetic variants with ESCC within the population of Jammu and Kashmir. The functional aspects of the variants understudy are yet to be evaluated. A recent study has observed that genetic variants that have been previously associated with one cancer are associated with other cancers too. Cross cancer analysis has been done for the identification of variants to estimate the genetic correlation between them by using the data from GWAS. Substantial evidence has been found which has identified pleiotropy among loci with strong association. A new high throughput study must be used to evaluate the functional implication of the variants [17].

The genes (Table 1) that were not associated with the ESCC in the studied population include TERT(Telomerase Reverse Transcriptase) which is highly active in cancer cells but shows low or inactivity in normal somatic cells. rs2735940 was significantly associated with a higher risk of ESCC in the Chinese population. This genetic variant has previously been associated with breast cancer in the European population [18, 19]. rs751402 of ERCC5 (Excision Repair Cross Complementing Group 5) plays an important role in DNA damage and repair and its deficiency can lead to genomic instability and carcinogenesis [20]. There are two genetic variants rs751402 (C > T) and rs2298881 (A/C/T) rs751402 that have been associated with ESCC in the Chinese population and rs2298881 has been associated with other cancers like lung, gastric and laryngeal cancers across the population [21,22,23]. rs2699887 of PIK3CA(Phosphatidylinositol-4-5-Biphosphate 3-Kinase Catalytic Sub- Unit Alpha) has been associated with ESCC in the USA (Texan population). rs3792152 of REV1(REV1, DNA Directed Polymerase) recruits DNA polymerases involved in translesion synthesis (TLS) of damaged DNA. Genetic variant rs3792152 of REV1 is associated with breast cancer risk in Thai women [24] and rs2981582 ofFGFR2(Fibroblast Growth Factor Receptor 2) plays a vital role in the cell proliferation, migration, and apoptosis [25]. Genetic variant rs2981582 of FGFR2 has been associated with breast cancer in different women population (Dutch, Arabic, and West Siberia) [26,27,28]. rs1695 of GSTP1(Glutathione S-transferase Pi 1) belongs to a family of enzymes that play an important role in detoxification by catalyzing the conjugation of many hydrophobic and electrophilic compounds. The genetic variant rs1695 of GSTP1 has shown an increased risk of ESCC and EAC in the population of Kashmir Valley [29]. In the present study, though the genetic variant rs1695 of GSTP1 has shown no significant association with ESCC. Similar results were obtained after stratification by race (Caucasian/Asian). Our results were in agreement with the above-mentioned meta-analysis of 52 studies [30]. rs2229080 of DCC(Deleted in Colorectal Carcinoma Netrin1 Receptor) has been studied in the Chinese population concerning breast cancer and it was found that this variant has been associated with reduced breast cancer risk [31]. The genetic variant rs2229080 (C > G) has been studied in ESCC and gastric cancer patients in Kashmir Valley also and it has been observed that rs2229080 of DCC has shown no association with the risk of ESCC. Various studies have shown the role of DCC in colorectal carcinoma, gastric, breast, esophageal and prostate cancer [32]. rs1801018 of BCL2(B-cell lymphoma 2) contributes to programmed cell death and apoptosis and has been studied in ESCC in the Chinese population but there was no significant association observed between the variant and the risk of the disease [33]. All of the above-mentioned SNPs have shown non-significant association with ESCC within the population of J&K in the present study. Although, these variants need to be investigated in the large cohort of population for a conclusive statement.


The Present study explores the link between genetics, environmental factors, and ESCC. This is the first study to investigate these genetic variants within the population of Jammu and Kashmir about esophageal squamous cell carcinoma. The present study identified important regions of genetic variation associated with risk for the development of the disease. Understanding of these biomarkers will help in elucidating the biological pathways and possible new strategies for identification and prevention of the malignancies. Though it has to be replicated in the large cohort for a conclusive statement but it is an important study that can establish the clinical relevance of novel biomarkers.

Availability of data and materials

Data and material are available. The datasets generated or analyzed during the current study are not publicly available but are available with the corresponding author and can be provided on reasonable request.



Body Mass Index


Confidence Interval


Esophageal Cancer


Esophageal Squamous Cell Carcinoma


Hardy Weinberg equilibrium


Jammu and Kashmir


Next generation sequencing


Normalized effect Size


Odds Ratio


Polymerase chain reaction


Single nucleotide polymorphism


Shrimp alkaline phosphatase


  1. 1.

    Zhang Y. Epidemiology of esophageal cancer. World J Gastroenterol: WJG. 2013;19(34):5598.

    Article  Google Scholar 

  2. 2.

    Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.

    Article  Google Scholar 

  3. 3.

    Samarasam I. Esophageal cancer in India: current status and future perspectives. Int J Adv Med Health Res. 2017;4(1):5.

    Article  Google Scholar 

  4. 4.

    Mir MM, Dar NA. Esophageal cancer in Kashmir (India): an enigma for researchers. Int J Health Sci. 2009;3(1):71.

    Google Scholar 

  5. 5.

    Ghatak S, Muthukumaran RB, Nachimuthu SK. A simple method of genomic DNA extraction from human samples for PCR-RFLP analysis. J Biomol Tech. 2013;24(4):224.

    PubMed  PubMed Central  Google Scholar 

  6. 6.

    Gabriel S, Ziaugra L, Tabbaa D. SNP genotyping using the Sequenom MassARRAY iPLEX platform. Curr Protoc Hum Genet. 2009;60(1):2–12.

    Google Scholar 

  7. 7.

    Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81(3):559–75.

    CAS  Article  Google Scholar 

  8. 8.

    O'dushlaine C, Kenny E, Heron EA, Segurado R, Gill M, Morris DW, et al. The SNP ratio test: pathway analysis of genome-wide association datasets. Bioinformatics. 2009;25(20):2762–3.

    CAS  Article  Google Scholar 

  9. 9.

    Yong Zou G, Donner A. The merits of testing hardy-Weinberg equilibrium in the analysis of unmatched case-control data: a cautionary note. Ann Hum Genet. 2006;70(6):923–33.

    Article  Google Scholar 

  10. 10.

    Yang CS. Research on esophageal cancer in China: a review. Cancer Res. 1980;40(8 Part 1):2633–44.

    CAS  PubMed  Google Scholar 

  11. 11.

    Zhang H-Z, Jin G-F, Shen H-B. Epidemiologic differences in esophageal cancer between Asian and Western populations. Chin J Cancer. 2012;31(6):281.

    Article  Google Scholar 

  12. 12.

    Kumar Phukan R, Kanta Chetia C, Shahadat Ali M, Mahanta J. Role of dietary habits in the development of esophageal cancer in Assam, the north-eastern region of India. Nutr Cancer. 2001;39(2):204–9.

    Article  Google Scholar 

  13. 13.

    Sehgal S, Kaul S, Gupta B, Dhar M. Risk factors and survival analysis of the esophageal cancer in the population of Jammu, India. Indian J Cancer. 2012;49(2):245.

    CAS  Article  Google Scholar 

  14. 14.

    Jiang Z, Zhang W, Chen Z, Shao J, Chen L, Wang Z. Transcription factor 21 (TCF21) rs12190287 polymorphism is associated with osteosarcoma risk and outcomes in east Chinese population. Med Sci Monit. 2017;23:3185.

    Article  Google Scholar 

  15. 15.

    Long J-R, Kataoka N, Shu X-O, Wen W, Gao Y-T, Cai Q, et al. Genetic polymorphisms of the CYP19A1 gene and breast cancer survival. Cancer Epidemiol Prevent Biomark. 2006;15(11):2115–22.

    CAS  Article  Google Scholar 

  16. 16.

    Warsy AS, Almukaynizi FB, AlDaihan S, Alam S, Daghastani M. Genetic Polymorphisms in Aromatase (CYP19) Gene and Cancer, Genetic Polymorphisms, Narasimha Reddy Parine, IntechOpen; 2017. Available from:,

  17. 17.

    Wu Y-H, Graff RE, Passarelli MN, Hoffman JD, Ziv E, Hoffmann TJ, et al. Identification of pleiotropic cancer susceptibility variants from genome-wide association studies reveals functional characteristics. Cancer Epidemiol Prevent Biomark. 2018;27(1):75–85.

    CAS  Article  Google Scholar 

  18. 18.

    Zhou L, Fu G, Wei J, Shi J, Pan W, Ren Y, et al. The identification of two regulatory ESCC susceptibility genetic variants in the TERT-CLPTM1L loci. Oncotarget. 2016;7(5):5495.

    Article  Google Scholar 

  19. 19.

    Bernert M. Investigating telomere dynamics in oesophageal squamous carcinoma cells using standard and gold nanoparticle-based assays; 2017.

    Google Scholar 

  20. 20.

    Zhu M-L, Shi T-Y, Hu H-C, He J, Wang M, Jin L, et al. Polymorphisms in the ERCC5 gene and risk of esophageal squamous cell carcinoma (ESCC) in eastern Chinese populations. PLoS One. 2012;7(7):e41500.

    CAS  Article  Google Scholar 

  21. 21.

    Gao H, Ge R, Liu H, Wang Y, Yan S. Effect of ERCC1 polymorphism on the response to chemotherapy and clinical outcome of non-small cell lung cancer. Genet Mol Res. 2014;13(4):8997–9004.

    CAS  Article  Google Scholar 

  22. 22.

    Lu Z-m, Luo T-h, M-m N, G-e F, Ma L-y, Xue X-c, et al. Influence of ERCC1 and ERCC4 polymorphisms on response to prognosis in gastric cancer treated with FOLFOX-based chemotherapy. Tumor Biol. 2014;35(4):2941–8.

    CAS  Article  Google Scholar 

  23. 23.

    Couto P, Miranda D, Vieira R, Vilhena A, De Marco L, Bastos-Rodrigues L. Association between CLOCK, PER3 and CCRN4L with non-small cell lung cancer in Brazilian patients. Mol Med Rep. 2014;10(1):435–40.

    CAS  Article  Google Scholar 

  24. 24.

    Sangrajrang S, Sato Y, Sakamoto H, Ohnami S, Khuhaprema T, Yoshida T. REV1 genetic polymorphisms and breast cancer risk in Thai women. Genes Environ. 2011;33(4):167–72.

    CAS  Article  Google Scholar 

  25. 25.

    Stelzer G, Rosen R, Plaschkes I, Zimmerman S, Twik M, Fishilevich S, Iny Stein T, Nudel R, Lieder I, Mazor Y, Kaplan S, Dahary D, Warshawsky D, Guan - Golan Y, Kohn A, Rappaport N, Safran M, and Lancet D. The GeneCards Suite: From Gene Data Mining to Disease Genome Sequence Analysis, Current Protocols in Bioinformatics, GeneCards – the human gene database, 2016;54:1.30.1-1.30.33.[PDF].

  26. 26.

    Shan J, Mahfoudh W, Dsouza SP, Hassen E, Bouaouina N, Abdelhak S, et al. Genome-wide association studies (GWAS) breast cancer susceptibility loci in Arabs: susceptibility and prognostic implications in Tunisians. Breast Cancer Res Treat. 2012;135(3):715–24.

    CAS  Article  Google Scholar 

  27. 27.

    Boyarskikh UA, Zarubina NA, Biltueva JA, Sinkina TV, Voronina EN, Lazarev AF, et al. Association of FGFR2 gene polymorphisms with the risk of breast cancer in population of West Siberia. Eur J Hum Genet. 2009;17(12):1688.

    CAS  Article  Google Scholar 

  28. 28.

    Huijts PE, Vreeswijk MP, Kroeze-Jansema KH, Jacobi CE, Seynaeve C, Krol-Warmerdam EM, et al. Clinical correlates of low-risk variants in FGFR2, TNRC9, MAP3K1, LSP1 and 8q24 in a Dutch cohort of incident breast cancer cases. Breast Cancer Res. 2007;9(6):R78.

    Article  Google Scholar 

  29. 29.

    Malik MA, Upadhyay R, Mittal RD, Zargar SA, Mittal B. Association of xenobiotic metabolizing enzymes genetic polymorphisms with esophageal cancer in Kashmir Valley and influence of environmental factors. Nutr Cancer. 2010;62(6):734–42.

    CAS  Article  Google Scholar 

  30. 30.

    Wang S, Zhang J, Jun F, Bai Z. Glutathione S-transferase pi 1 variant and squamous cell carcinoma susceptibility: a meta-analysis of 52 case-control studies. BMC Med Genet. 2019;20(1):22.

    Article  Google Scholar 

  31. 31.

    Liu X, Wang X, Fu SW, Wang M, Kang H, Guan H, et al. Genetic association of deleted in colorectal carcinoma variants with breast cancer risk: a case-control study. Oncotarget. 2016;7(22):32765.

    Article  Google Scholar 

  32. 32.

    Malik MA, Gupta A, Zargar SA, Mittal B. Role of genetic variants of deleted in colorectal carcinoma (DCC) polymorphisms and esophageal and gastric cancers risk in Kashmir Valley and meta-analysis. Tumor Biol. 2013;34(5):3049–57.

    CAS  Article  Google Scholar 

  33. 33.

    Li Q, Yin J, Wang X, Wang L-M, Shi Y-J, Zheng L, et al. B-cell lymphoma 2 rs17757541 C> G polymorphism was associated with an increased risk of gastric cardiac adenocarcinoma in a Chinese population. Asian Pac J Cancer Prev. 2013;14(7):4301–6.

    Article  Google Scholar 

Download references


RS would like to acknowledge Dr. Swarkar Sharma, for his guidance and support. He conceptualised the study and guided throughout the experimental work.


No Funding.

Author information




RS designed the study, performed experiments, extrapolated results and has written the manuscript. RS acknowledges SV, AB and GR in performing the experimental work, BS and DB in sample collection. RS would like to acknowledge VS, IS, HS and RK in editing the manuscript. RS acknowledges NAD for providing samples. All authors have read and approved the manuscript.

Corresponding authors

Correspondence to Ruchi Shah or Nazir Ahmed Dar.

Ethics declarations

Ethics approval and consent to participate

A total of 758 individuals (166 cancerous and 592 healthy control were recruited for the study. The study was approved by the Institutional Ethics Review board (IERB) of Shri Mata Vaishno Devi University (SMVDU) vide IERB Serial No: SMVDU/IERB/16/43. The consent form was designed in three languages (English, Hindi, Urdu) and an informed and written consent was obtained from participants (cases and controls) before the sample collection.

Consent for publication

All data is coded, no personal information, images have been used in the study.

Competing interests

No competing interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Supplementary Table 1. Details of the genes selected for the study. Supplementary Table 2. Details and clinical features of the cases and controls of J&K Population. Supplementary Table 3. SNPs and primers associated.

Additional file 2: Supplementary Figure 1. Gene expressions of CYP19A1 (fig. 1a) and TCF21 (fig. 1b) in different tissues. Violin plots showing expression of CYP19A1 and TCF21 in different tissues which includes esophageal tissue. On the basis of TPM value (Transcript per million) CYP19A1 is less expressed in esophageal tissue whereas TCF21 is well expressed in Esophageal tissue and NES value for variants rs10046 of CYP19A1 is negative and rs12190287of TCF21 is positive. Supplementary Figure 2. Gene location and variant effect prediction. TSS is transcription start site and TES is transcription end site predicted by using The Genotype-Tissue Expression (GTEx) locus. It aims to characterize variation in gene expression levels across individuals and diverse tissues of the human body. Both CYP19A1 (fig. 2a) and TCF21 (fig. 2b) has a strong eQTL (expression quantitative trait loci) signal in esophageal tissue.

Additional file 3: Supplementary file 1. Datasheet and consent form for sample collection.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Shah, R., Sharma, V., Bhat, A. et al. MassARRAY analysis of twelve cancer related SNPs in esophageal squamous cell carcinoma in J&K, India. BMC Cancer 20, 497 (2020).

Download citation


  • Single nucleotide polymorphism
  • Esophageal cancer
  • Jammu and Kashmir