METHODS
The genetic variants of CYP2A6 and CYP2A13 genes were investigated by polymerase chain reaction-
restriction fragment length polymorphism, allele-specific PCR and sequencing methods in 492
histopathologically confirmed ESCC cases and an equal number of matched controls. Gene-gene and
gene-environment interactions were calculated by logistic regression analysis.
RESULTS
Inverse association was found between the variant genotypes of CYP2A6 (OR=0.6; 95% CI, 0.4-0.9) and
CYP2A13 (OR=0.5; 95% CI, 0.31-0.8) with ESCC risk. Individually, the inverse association of variant
genotypes of the three studied CYP2A6 genes was retained when harboured by a participant in combination
with CYP2A13 variant genotype (OR=0.3; 95%CI, 0.1-0.8). Participants who were smokers,
consumed alkaline beverage, had used biomass fuel for cooking, lived in adobe houses and had a positive
family history of cancer showed a strong ESCC risk when harbouring homozygous wild genotypes
of CYP2A6 and CYP2A13. Among the different gene environmental interactions, only CYP2A6b
(OR=1.5; 95%CI, 1.1-2.0; Pinteraction=0.018), and CYP2A13 (OR=1.4; 95%CI, 1.1-2.0; Pinteraction=0.021)
genotypes showed statistically significant interactions with smoking.
CONCLUSION
Normal genotypes of CYP2A6 and CYP2A13 considerably increase ESCC risk in subjects who also had
exposure to environmental risk factors.
Keywords: Cytochrome P450 gene variants; gene-environment interaction; Kashmi; oesophagal squamous cell carcinoma
Polymorphisms in xenobiotic metabolising enzymes
(XME) are among the most important genetic
differences in carcinogenesis. The Phase I and Phase
II XMEs biotransform the xenobiotic to make it more
water-soluble, and the Phase III transporter then aids
in the elimination of the changed intermediate from
the body via urine. The intermediates produced during
biotransformation are more reactive and, if not eliminated,
can bind biomolecules, including DNA. Interindividual
genetic variants in XMEs cause differences
in XME expression and activity, as well as the ability to
eliminate reactive intermediates from the body, thereby
modulating the risk of cancer caused by xenobiotics
in food or the environment.[
Kashmir, with a high incidence of ESCC in both
genders,[
Data Collection
Detailed information on age, sex, place of residence,
ethnicity, religion, education, dietary data, including
intake of fresh fruits and vegetables and other potential
confounding factors of interest was collected using
a questionnaire specifically designed for the study
population. Detailed information on the life-long history
of use, with starting and stopping ages and daily
amount of use, was obtained for several tobacco products.
Any change in the type of tobacco products and
amount of use was also recorded. Ever use of traditional
hookah, nass, cigarette and gutkha, was defined
as the use of the respective product (s) daily or at least
weekly for a period of 6 months or more. Information
on family history of any cancer (FHC) was obtained
from all the participants. To assess the socioeconomic
status (SES) of the subjects, information on
potential parameters of SES was obtained including
education level (highest level attained), monthly income
(INR), house type, cooking fuel, and ownership
of several household appliances. Similarly, the information
regarding oral hygiene, house type, second
hand smoking was also acquired from all the subjects.
Genetic analysis
For identification of the CYP2A6 genotypes, PCRRFLP
analysis was performed as described previously.
[
Allele Specific-Polymerase Chain Reaction (ASPCR)
was employed for the CYP2A6b gene, and whole
gene deletion genotyping was based on a 2-step PCR
method. The first PCR reaction produced a 1,961bp
fragment of the CYP2A6b from all individuals with or
without the deleted CYP2A6b gene. The second PCR,
which specifically detected the deleted CYP2A6b gene,
used the product resulting from the first PCR amplification
as a template. 1.5% ethidium bromide-stained
agarose gel was run to check the amplified products.
The presence of the CYP2A6b-specific 1,181bp product
amplified with the first primer pair indicated the CYP2A6b
wild genotype (*1/*1). The presence of the product resulting from amplification with the second primer
pair indicated the deleted CYP2A6b genotype (*4/*4).
The presence of the product in both reactions indicated
the heterozygote genotype (*1/*4). It is important to
mention that 10% of the samples of cases and controls
were randomly tested twice for experimental validation;
however, the results were similar for all duplicate sets.
PCR amplification of CYP2A6c genotypes, including
homozygous for the wild-type (*1A/*1A), heterozygous
type (*1A/*4C) and deletion-type (*4C/*4C)
were determined by PCR-RFLP assay as previously
described.[
The details of PCR conditions, primers, restriction
enzyme, and length of expected fragments on digestion,
mutant alleles and change in nucleotide position
of the above genes are given in Table
It is pertinent to mention that three different SNPs
were simultaneously studied in the case of the CYP2A6
gene based on their substrate specificity with nitrosamines
and PAHs to which the study population is frequently
exposed through various exogenous exposures.
Similarly, in the case of CYP2A13, the three genotypes
wild homozygous (C/C), heterozygous (C/T) and
homozygous mutant (T/T) were determined by PCRRFLP
assay as previously described.[
PCR-RFLP results were validated by sequencing
10% of the randomly picked samples. For sequencing,
unpurified PCR products were directly sent to SciGenome
Private Limited, Cochin Kerala-India. The resulting
sequence chromatograms were then compared
with the original gene sequences for the expected results.
Sequence scanner software (Finch TV Geospiza
1.4.0) was used for comparing sequences for the possible
sequence variations due to gene polymorphisms.
Statistical Analysis
Five millilitres of venous blood was collected from
each subject in sterilised plastic vials containing EDTA
(0.5M; pH=8.0) and stored at -80ºC before DNA extraction.
Genomic DNA was extracted from the collected
blood samples using the phenol-chloroform
method.[
Categorical variables were set for presenting and
calculating numbers and percentages for different
genotypes of CYP2A6 and CYP2A13. Tests for Hardy-
Weinberg Equilibrium (HWE) were conducted by
comparing observed and expected genotype distributions
by the ?2 goodness of fit. Statistical significance
for the departure of a genotype frequency from its expected
frequency under the HWE model was set at
p?0.05. Conditional logistic regression models were
used to calculate odds ratios (ORs) and corresponding
95% confidence intervals (95%CIs) to assess the association of the genotypes with ESCC risk and to
assess the possible gene-gene and gene-environment
interaction (GEI). Subjects were stratified into various
groups based on smoking habit, FHC, type of fuel
used for cooking, house type and various possible
genotypic combinations. Adjustment was made with
known ESCC risk factors like age, sex, residence, education
level, SES, fruit and vegetable consumption,
oral hygiene, animal contact, salted tea consumption
and smoking in different forms. For genetic analysis,
wild homozygous, heterozygous, mutant homozygous
genotypes, as well as a variant group (a combination
of heterozygous and mutant homozygous with at
least one and/or both defective alleles) were analysed
separately. However, for gene-gene or GEI analysis,
genotypes were restricted to homozygous wild and a
variant genotype only. All statistical analysis was done
using STATA software, version 12 (STATA Corp.,
College Station, TX, USA). Two-sided P values <0.05
were considered statistically significant.
Allelic Frequencies and Genotype Analysis
The minor allele frequency differences observed
among cases and controls were statistically significant
(p=<0.05) except for CYP2A6b (p=0.230), and the
genotype frequencies were in agreement with Hardy-
Weinberg Equilibrium (Appendix 1). PCR-RFLP results
of CYP2A6a, CYP2A6c and CYP2A13 genes are
presented in Appendices 2, 3 and 4, respectively.
As shown in Table
On analysing the modulating effect of various known risk factors of ESCC in the study population, the association of the various genotypes with ESCC risk changed reasonably in participants with wild-type genotypes as compared to variant genotype referents (heterozygous and homozygous mutant genotypes clubbed together as variants).
Tobacco Smoking and Nass Use
The smoker participants had higher risk of ESCC when
carry the wild type genotypes of CYP2A6a (OR=2.7;
95% CI, 1.3-5.3); CYP2A6b (OR=2.9; 95% CI, 1.3-6.9);
CYP2A6c (OR=3.0; 95% CI, 1.6-5.6) and CYP2A13
(OR=2.2; 95% CI, 1.2-4.0). Similarly, the use of nass
increased risk in participants when carrying wild-type
genotypes of CYP2A6 or CYP2A13 (Table
Salt Tea Consumption
Unlike the other genotypes, moderately higher ESCC
risk was found in CYP2A6C and CYP2A13 wild genotype
harbouring subjects with salt tea consumption.
Family History
The presence of a family history of cancer was strongly
associated with ESCC risk. The OR and 95% CI for
wild genotypes of CYP2A6a, CYP2A6b, CYP2A6c and
CYP2A13 were (OR=8.8; 95% CI, 4.5-17.9), (OR=5.2;
95% CI, 2.7-9.7), (OR=7.6; 95%CI, 4.0-14.6) and
(OR=7. 95%CI, 4.1-14.1), respectively (Table
Biomass Fuel Use and House Type
High risk of ESCC was found in participants who carried
the wild homozygous genotype of CYP2A6a and used biomass
as cooking fuel (OR=7.8; 95% CI, 2.1-29.0) and lived
in adobe houses (OR=4.7; 95%CI, 2.1-10.7). Similarly,
higher risk was found in participants who lived in adobe
house and used biomass fuel for cooking in other CYP2A6
and CYP2A13 analyzed wild genotypes (Table
Gender Wise Risk
On analysing the gender wise risk, males showed increased
risk while harboring wild genotypes of either
CYP2A6b (OR=2.1; 95% CI, 1.0-4.8); CYP2A6c
(OR=2.5; 95% CI, 1.1-5.4) and CYP2A13 (OR=3.5;
95%CI, 1.6-7.7) (Table
Gene Environment Interaction
Among the various gene environment combinations,
only CYP2A6b (Pinteraction=0.018) and CYP2A13 (Pinteraction
=0.021) genes showed synergistically significant
associations with smoking (Table
In case of the CYP2A6 gene, 23 variants are currently
known that reduce its enzymatic activity, and 5
variants were shown to completely abolish enzymatic
function.[
In this study, we observed an inverse association
of some variant genotypes of CYP2A6 and CYP2A13.
However, the ESCC risk increased significantly with
normal genotypes in the presence of different ESCC
risk factors in the study population.
The reduced risk of ESCC due to CYP2A6 gene
variants in our study is in agreement with the previous
studies on oesophagal and lung malignancies.
[
Similarly, frequencies of the CYP2A6*4 allele
in three regions of China, in Han (N=120), Uighur
(N=100), Bouyei (N=100) and Tibetan (N=100)
(p<0.0001) were 7.9%, 15%, 0% and 2%, respectively.
[
Besides a huge body of literature supporting our
findings, still a few studies that have not shown any significant
relationships between CYP2A6 genotypes and
lung cancer in both never and ever smokers.[
In case of CYP2A13 gene, our results were in line
with the previous reports which have revealed a 2-3
fold reduction in the metabolic activation of tobacco
specific nitrosamine - 4-(methylnitrosamino)-1-(3-
pyridyl)-1-butanone (NNK) in subjects with variant
genotype of CYP2A13 than subjects carrying a
wild genotype, [
Stratification analysis showed that the reduced
risk of lung adenocarcinoma related to the variant
CYP2A13 genotype was limited to smokers, especially
light smokers (OR=0.23; 95% CI, 0.08-0.68) but
not non-smokers or heavy smokers. The two novel
polymorphisms T478C and T494C in the CYP2A13
gene were associated with a significantly reduced risk
of head and neck cancer (OR=0.37; 95%CI, 0.19-
0.71). A CYP2A13 haplotype carrying variant alleles
of T478C/T494C was associated with the reduced risk
of (OR 0.42; 95%CI, 0.22-0.78).[
The retention of inverse effect due to the combination
of CYP2A6 and CYP2A13 variants in our study
is consistent with the earlier studies, but other than
oesophagal malignancy.[
One of the interesting observations of the current
study is that the ESCC risk is more common among
men than women. The biological mechanism for such
association is still not known, but one of the plausible
explanations could be that males have a higher
prevalence of active tobacco smoking as compared to
females. The combination of tobacco smoking exposure
in males with susceptible genotypes puts them
at higher risk and hence male predominance towards
ESCC risk in our population. Interestingly, a reduction
in the consumption of smoking among the carriers
of variant alleles of the CYP2A6 gene could also
be the probable reason for male dominance with CYP2A6
wild genotype carrying subjects, and hence an increased risk than female participants.[
A limited number of studies are available regarding
the modifying effect of the CYP2A6 and CYP2A13
genes in subjects with a positive history of cancer
among relatives. The plausible reason for increased risk
in our study could be either similar exposure to ESCC
risk factors within the family and/or similar genetic
setup among the relatives.[
The synergistic association of the CYP2A13 wild
genotype in the presence of smoking could reflect the
biological feature of CYP2A13-257Cys, which exhibits
a decreased catalytic efficiency toward N-nitrosamines
as compared to CYP2A13?257Arg.[
Limitations
Though, no study has reported the exposure of these
genotypes with ESCC under such environmental combinations
and confounding of the results with the probable
ESCC risk factors, selection or recall bias could be
one of the weak points of this study, although the same
hospital setting and limited number of interviewers
lessen this type of bias to some extent.
Ethics Committee Approval: The study was approved by the Sher-i-Kashmir Institute of Medical Sciences Ethics Committee (no: SIMS 1 31/IEC?SKIMS/2013, date: 17/01/2013).
Informed Consent: Informed consent was obtained from all participants.
Conflict of Interest Statement: The authors have no conflicts of interest to declare.
Funding: The authors declared that this study received no financial support.
Use of AI for Writing Assistance: No AI technologies utilized.
Author Contributions: Concept - N.A.D., G.A.B.; Design - N.A.D., G.A.B; Supervision - N.A.D.; Funding - N.A.D.; Materials - T.R.M., F.A.G.; Data collection and/or processing - G.A.B., G.R., N.A.S.; Data analysis and/or interpretation - N.A.B.; Literature search - N.A.B.; Writing - N.A.B.; Critical review - G.A.B., N.A.S., G.R., F.A.G., T.R.M., M.S., N.A.D.
Peer-review: Externally peer-reviewed.