Genetic Variant Rs2301721 in the HOXA7 Gene Is Associated with Leukemia Risk in Jammu Region of India
Total Page:16
File Type:pdf, Size:1020Kb
RESEARCH Genetic variant rs2301721 in the HOXA7 gene is associated with leukemia risk in Jammu region of India a b a a a a Amrita Bhat , Deepak Abrol , Gh. Rasool Bhat , Ruchi Shah , Sonali Verma , Divya Bakshi , Bhanu Sharma a c d a , Supinder Singh , Rakesh Raina , Rakesh Kumar * a Shri Mata Vaishno Devi University, Katra b Department of Radiotherapy, GMC, Kathua, J&K c Department of Medicine, ASCOMS, Sidhra. d Siddhi Vinayak Dental and Health Care Centre, Rupnagar, Jammu, India *Corresponding author e-mail: [email protected] ABSTRACT Background: Leukemia is a heterogeneous disorder, characterized by high proliferation of white blood cells. Various genetic studies have tried to assess the role of SNPs contributing to the development of leukemia. The HOX is implicated in normal as well as abnormal hematopoiesis. The role of the present variant rs2301721 has been previously studied in various population groups; however, the role of this variant in Leukemia in the Jammu region/population is unclear. Aim: In the present study, we investigated the genetic variant rs2301721 in the HOXA7 gene in leukemia patients from the Jammu region of India. Methods: The variant was genotyped by using Sanger sequencing in 180 individuals (90 leukemic cases and 90 healthy controls). The association of SNP with the disease was evaluated by using logistic regression. Results: It was observed that the variant rs2301721 showed significant association with Leukemia risk in Jammu population The allelic OR of variant rs2301721 was 1.54 (1.01 – 2.34) with p-value = 0.004 and under H.W.E it was 0.838 when corrected for age, gender, BMI, smoking and alcohol. Conclusion: The present study concludes that the variant rs2301721 in the HOXA7 gene acts as a risk factor in the development of leukemia in the study population. KEYWORDS: Homeobox proteins (HOXA7); Linkage Disequilibrium (LD); Jammu region; SNP Citation: Bhat et al. Genetic variant rs2301721 in the HOXA7 gene is associated with Leukemia risk in Jammu region of India. Polymorphism 2020; 5: 68-76. POLYMORPHISM 68 RESEARCH INTRODUCTION may regulate gene expression, morphogenesis, and Leukemia consists of a group of heterogeneous differentiation in hematopoietic stem cells (Inoue et malignancies in which immature and dysfunctional al., 2013). The present study aimed at exploring the hematopoietic progenitors proliferate and association of variant rs2301721 of HOXA7 with the accumulate in the bone marrow (Alharbi, risk of Leukemia in Jammu population. Pettengell, Pandha, & Morgan, 2013). In the hematopoietic stem cells, a disruption of the cellular MATERIAL AND METHODS processes, including proliferation, differentiation, Ethics statement and cell growth ultimately causes Leukemia The study was approved by the Institutional Ethics ("Leukemia," 2017). In the United States, leukemia is Review Board (IERB) of Shri Mata Vaishno Devi ranked 8th for being the most commonly diagnosed University (SMVDU) under notification number cancer with approximately 61,789 new cases and (SMVDU/IERB/16/41). All the details have been 23,100 deaths in 2019. According to the Population- recorded in a predesigned proforma and written Based Cancer Registry of India, males are more informed consent was obtained from each affected than females by the ratio 2:1 (Asthana, participant before conducting the study. All 2018). The incidence of leukemia in Jammu has shot experimental protocols were conducted according up rapidly from the previous decade to the guidelines and regulations set by IERB, (Chokkalingam & Buffler, 2008). One of India’s SMVDU. highly diverse and predominately inbred population groups, the population of Jammu and Kashmir (J&K) region, remained unexplored for Sample collection and DNA isolation blood disorders. Leukemia is multifactorial in origin, A total of 180 subjects were recruited for the study, which can be caused by both genetic as well as out of which 90 were cases (leukemic patients) and non-genetic factors (Sinnett, 2000). Various genetic 90 were healthy controls. All the cases were alterations, including genetic polymorphisms histopathologically confirmed. The genomic DNA resulting in gain-of-function mutation in proto- was isolated from the blood samples using a oncogenes as well as inactivating tumor suppressor Qiagen DNA Isolation kit (Catalogue No. 51206). genes can lead to the progression of leukemia Agarose gel electrophoresis was used to analyze (Münger, 2002). Genome-wide association studies the quality of the genomic DNA and quantification (GWAS) have advanced our understanding of was performed using UV spectrophotometer susceptibility to leukemia; however, much of the (Nanodrop). heritable risk remains unidentified. This factor has been increasingly documented as an inherited Exome sequencing genetic factor that can contribute to the progression of leukemia. GWAS data indicate that The exome sequencing for 05 cases was done on much of the heritable risk of leukemia is ascribable Illumina NGS Platform. The alignment against the to common genetic variations which need to be hg GrCh37 genome was carried out using BWA and discovered. Recently, GWAS has found the the variant calling was done using GATK pipeline. association of variant rs2301721 of HOXA7 with The alignment against the hg19 GrCh37 genome depression disease in Russian population (Shadrina, was carried out using BWA, and the variant calling Bondarenko, & Slominsky, 2018) and colon cancer was done using GATK pipeline. The variants were in Korean population (Kan et al., 2007). The genetic identified using GATK tool and annotated using variant rs2301721 of HOXA7 has cytogenetic ANNOVAR. The databases used for study include localization to 7p15.2 with the 1st exonic position and 1000 genome, DBNSFP, dbSNP, GWAS and ClinVar. encodes a DNA-binding transcription factor which POLYMORPHISM 69 RESEARCH After filtration of variants against various filters by the use of Exome Capture Kits targeted region. RESULT The total data was generated 45.63 GB with more Sanger Sequencing than 100X sequencing coverage. HOXA7 variant Polymerase chain reaction (PCR) was done to rs2301721 was frequently observed in all samples. amplify PAX5 using Mastercycler® pro, Eppendorf Detail of gene is given in Table 2 and Figure 1 AG, Hamburg, Germany. Overall the target region HOXA7 was not been reported in leukemia from was amplified using the PCR program with 95ºC for Jammu Kashmir population of North India. As the 1 minute, 95ºC for 1 min, 60ºC for 45 seconds with variant of HOXA7 was shortlisted from WES, it was 36 cycles, 72ºC for 1 min, 72ºC for 2min, 4ºC infinity imperative to further validate in a larger number of (∞) shown in Table 1 and sequencing was samples and the same was done by sanger performed by Aggrigenome Pvt. Limited using sequencing. The main aim was to explore the Sequence Scanner Software (ABI 3730 XL DNA association of variant rs2301721 of HOXA7 with the Analyser) and Chromas v2.6.6 (Technelysium Pty Leukemia risk. Ltd, South Brisbane, Australia). This case-control association study included 180 samples (90 leukemic cases and 90 healthy Statistical Analysis controls). The clinical characteristic distributions of Statistical Analysis of the data was done by using the cases and the controls are given in Table 2. In SPSS software (v.20; Chicago, IL). Chi-square (χ2) Table 1, 60% of Leukemia patients were males, and was performed, and genotype frequencies were 40% were females, while in controls, 70% were tested for total Hardy–Weinberg equilibrium. Binary males, and 30% were females. The mean ages were Logistic Regression was also used to estimate OR at 40.51(±14.67) for leukemia patients and 50.76 95% CI, and the respective level of significance was (±13.30) years for controls. BMI in cases was 21.2 estimated as p-value. Moreover, various (±6.10) and 24.9 (±5.3) in Controls. interactions of genes were analyzed by the string software tool. Table: 1 Showing the details of forwarding and reversed primers of variant rs2301721 of HOXA7. Gene Zygosity SNP ID cDNA Protei Exon Function Type of Alteration Change n No. chang e HOXA7 Heterozygosit rs2301721 g.2862429 p.A18T 01 Mis-sense Non synonymous y 4 G>A variant SNV POLYMORPHISM 70 RESEARCH Table 2: Clinical Characteristics of Cases and Controls. Characteristics Cases Controls p-value Age * 40.7 (±14.5) 49.1(±15.30) <0.01 Gender (in %) F=40 F=30 - M=60 M=70 BMI** 21.2 (±6.10) 24.9 (±5.3) <0.01 Histopathology (Leukemia Types) Acute Lymphoblastic Leukemia (ALL) 0.17 Acute Myeloid Leukemia (AML) 0.23 Chronic Myeloid Leukemia (CML) 0.50 Chronic Lymphoid Leukemia (CLL) 0.10 Smoker 0.63 0.22 - Non-Smoker 0.37 0.78 - Alcoholic 0.53 0.17 - Non-Alcoholic 0.47 0.83 - *age in years and **BMI in kg/m2 GENES TERT HOXA7 8% ATM 17% 8% TP63 9% PAX5 14% XRCC1 11% ASXL1 11% IKZF1 11% ARID5B 11% Figure 1: Showing the percentage frequency of genes targeted in exome sequencing, where HOXA7 had a higher percentage (17%) frequency among all genes. The C>T genetic variant (rs10899750) in the HOXA7 HOXA7 shows homozygous wild (CC), gene results in an amino acid change from Alanine heterozygous (CT) and Homozygous mutant (TT) to tryptophan. An allelic frequency distribution is (Figure 3). shown in Figure 2 and the chromatogram of POLYMORPHISM 71 RESEARCH CASES CONTROLS 40 35 30 25 20 15 Frequency 10 5 0 CC CT TT Allele Figure 2: Showing Allelic frequency distribution of variant rs2301721 of HOXA7 Figure 3: Chromatogram of HOXA7 with variant rs2301721 having missense mutation C >T. Showing A) homozygous wild (CC) B) heterozygous (CT) and C) Homozygous mutant (TT). The allelic OR of variant rs2301721 was 1.54 (1.01 – HOXA7 shows a significant association with the risk 2.34) with p-value = 0.004 and under H.W.E it was of leukemia in the population of Jammu region of 0.838.