![]() |
|
|
| Abstract |
|---|
|
|
|---|
Genetic polymorphisms that affect the metabolic activity or expression of
biotransformation enzymes may be important contributors to interindividual
differences in susceptibility to environmental diseases (e.g.,
Ariyoshi et al., 2002
;
Thier et al., 2002
;
Xu et al., 2002
). For CYP2A13,
variations in its expression level or metabolic capacity may significantly
alter the extent of metabolic activation of NNK and other xenobiotic
substrates in the respiratory tract, which may lead to altered susceptibility
to tobacco-related tumorigenesis. In a previous study, we used polymerase
chain reaction-single-strand conformation polymorphism (PCR-SSCP) analysis to
investigate the genetic polymorphisms of the CYP2A13 gene
(Zhang et al., 2002
). A total
of seven variant alleles was detected, but only four represented single
nucleotide polymorphisms (SNPs) (i.e., variants having allele frequencies
greater than 1%), and only one SNP was detected in the coding region, in exon
5, leading to an Arg257Cys amino acid change. The Arg257Cys variant has been
characterized following heterologous expression; it was 37 to 56% less active
than the 257Arg protein toward all substrates tested in a reconstituted
system, and it displayed a >2-fold decrease in catalytic efficiency toward
NNK (Zhang et al., 2002
).
Because PCR-SSCP may not detect all mutations, it is possible that additional variant alleles occur in the CYP2A13 gene, which remain to be identified by the use of more sensitive techniques. In the present study, we have developed a long-PCR protocol to amplify the CYP2A13 gene for sequence analysis. A total of 32 samples from Chinese donors was sequenced, which led to the identification of 13 previously unidentified SNPs in the CYP2A13 gene. These include three SNPs in the coding region, one of which leads to a nonfunctional allele. The frequencies of the coding region mutations were further examined using random samples of white, black, Hispanic, and Asian newborns from New York. In addition, haplotype analysis was performed to identify linkages between the various mutations.
| Materials and Methods |
|---|
|
|
|---|
Amplification of CYP2A13 Gene by Long-Distance PCR. Primers shown in Table 1 were designed according to the CYP2A13 gene sequence from the completed human genome data base (GenBank accession no. AC008962 [GenBank] ). Primers 5'-F1 and 3'-R2 were used to amplify the full-length CYP2A13 gene. PCR amplification was carried out in a PerkinElmer thermal cycler 9600 instrument (Applied Biosystems, Foster City, CA). The reaction mixtures, in a total volume of 50 µl, contained about 50 to 200 ng of genomic DNA or 50 to 200 pg of PCR product as a template, 1.8 mM MgSO4, 60 mM Tris-SO4 (pH 9.1), 18 mM (NH4)2SO4, 200 µM each dNTP, 200 nM each primer, and 2 µl of ELONGase Enzyme Mix (Invitrogen, Carlsbad, CA). After an initial denaturation at 94°C for 2 min, 35 cycles of amplification, each consisting of a denaturation at 94°C for 45 s, annealing at 64.5°C for 45 s, and an extension at 68°C for 10 min, were carried out, followed by a final extension at 68°C for 10 min. An H2O blank (no template) control was routinely used for detecting potential contamination of reagents.
|
DNA Sequencing. All long-PCR products were gel-purified using a QIAquick Gel Extraction kit (QIAGEN) and then subjected to direct sequencing using an automated DNA sequencer from Applied Biosystems (model 3100) at the Molecular Genetics Core of the Wadsworth Center. Each sample was initially sequenced with a total of 10 different primers, as follows: 5'-F1, 5'-R1, Exon1R, Exon2F2, Exon3F2, Exon5F, Exon6R, Exon7R, Exon8R, and Exon9R2 (see Table 1 for sequences). When PCR product from a single amplification was insufficient for all sequencing reactions, combined PCR products from two or more amplifications using the same genomic DNA sample or the first PCR product as a template were used. Mutations identified were confirmed by a second, independent PCR and by sequencing experiments using primers shown in Table 1, for detection of potential PCR and sequencing errors.
Determination of Allele Frequencies. The frequencies of exon 1, 2,
and 3 variant alleles were determined using anonymous human genomic DNA
samples isolated from newborn blood spots from the New York State Newborn
Screening Program. Use of these specimens was approved by the Institutional
Review Board at the New York State Department of Health. The methods of DNA
preparation and of PCR amplification of DNA fragments covering exon 1, 2, and
3 have been described previously (Zhang et
al., 2002
). The lower and upper limits of the 95% confidence
interval for the allele frequency were calculated, with a correction for
continuity, using a program
(http://faculty.vassar.edu/lowry/prop1.html)
that is based on methods described by Newcombe
(1998
).
Haplotype Analysis. DNA samples containing heterozygous mutations at both exon 1 and exon 5 SNP sites were used as a template to amplify the full-length CYP2A13 gene using primers 5'-F1 and 3'-R2, as described above. The 9456-bp PCR product was used as a template in a nested PCR using primers 5'-F2 and 3'-R1, which amplify an 8723-bp fragment. The nested PCR mixtures contained, in a total volume of 50 µl, about 50 to 200 pg of template, 1.5 mM MgSO4, 60 mM Tris-SO4 (pH 9.1), 18 mM (NH4)2SO4, 200 µM each dNTP, 200 nM each primer, and 1 µl of ELONGase Enzyme Mix. After an initial denaturation at 94°C for 2 min, 35 cycles of amplification, each consisting of a denaturation at 94°C for 20 s, annealing at 66°C for 1 min, and an extension at 68°C for 8 min, were carried out, followed by a final extension at 72°C for 30 min. The nested PCR products were purified by agarose gel electrophoresis, with use of crystal violet for staining DNA (Invitrogen), and cloned into the pCR-XL-TOPO plasmid vector (Invitrogen). Sequence analysis of individual clones, which represent single alleles, was carried out using primers 5'-R1, Exon1R, Exon5F, and Exon9R2.
| Results |
|---|
|
|
|---|
|
A total of 14 SNPs were detected, as shown in Table 2; all variant alleles were detected as heterozygotes. Of the 14 SNPs, four are located in the coding region, whereas the other 10 are located in noncoding regions: two in the 5'-flanking region, two in the 3'-untranslated region, and six in the introns. Five of the 14 SNPs had apparent allele frequencies of greater than 10%, whereas the other nine had apparent allele frequencies of 1.8 to 4.8% in this small sample.
|
Among the four SNPs detected in the coding region, one (Arg257Cys, caused
by a 3375C > T missense mutation, with an apparent allele frequency of
10.9%) was the same as had been previously identified
(Zhang et al., 2002
), whereas
the other three (one each in exons 1, 2, and 3) had not been previously
reported. The SNP in exon 1 was a 74G > A missense mutation, leading to a
predicted amino acid change from Arg to Gln at position 25. This SNP was
detected in seven subjects, with an apparent allele frequency of 10.9%. The
SNP detected in exon 2 was a 578C > T missense mutation, resulting in an
amino acid change from an Arg to a stop codon at position 101. This SNP was
detected in two subjects, with an apparent allele frequency of 3.2%. The SNP
in exon 3, which was detected in only one subject (1.8%), was a missense
mutation (1706C > G), which is expected to cause a conserved Asp158Glu
substitution.
None of the 10 SNPs found in the noncoding regions had been detected in the
previous study (Zhang et al.,
2002
). The two SNPs in the 5'-flanking region, -729C > T
(at 729 bp before ATG) and -411G > A, had apparent allele frequencies of
4.7 and 15.6%, respectively, whereas the two SNPs in the 3'-untranslated
region, 7520C > G and 7571G > C, had apparent allele frequencies of 4.7
and 10.9%, respectively. The other six SNPs were detected in introns 1, 2, 3,
5, 6, and 8 (Table 2), with
apparent allele frequencies between 1.8 and 15.6%.
A rare mutation in exon 3, 1662G > C (Gly144Arg), previously detected in
a Hispanic newborn (Zhang et al.,
2002
), and a silent mutation in exon 2, 523C > T, previously
detected in an Asian newborn (Zhang et
al., 2002
), were not detected in these Chinese samples. The
occurrence of the intron-5 and intron-7 mutations reported earlier
(Zhang et al., 2002
) was not
examined in the present study.
The frequencies of the Arg25Gln, Arg101Stop, and Asp158Glu alleles were further examined using random samples of white, black, Hispanic, and Asian newborns in the New York State Newborn Screening Program. For Arg25Gln, nine heterozygotes and two homozygotes for the 74G > A mutation were detected among the 102 samples analyzed, with apparent allele frequencies of 2.1% (95% confidence interval, 0.112.5%) in white (n = 24), 11.5% (4.824.1%) in black (n = 26), 1.9% (0.111.6%) in Hispanic (n = 26), and 9.6% (3.621.8%) in Asian (n = 26) newborns. The frequency in the Asian newborns (9.6%) was very similar to that found in the Chinese samples (10.9%). On the other hand, the Arg101Stop mutation was not detected in 136 newborn samples examined (23 white, 21 black, 19 Hispanic, and 73 Asian), suggesting that this mutation may be unique for the Chinese patient population. For Asp158Glu, one heterozygote (Hispanic) was detected among 86 newborns examined (17 white, 21 black, 25 Hispanic, and 23 Asian), indicating that this allele, although infrequent, is not exclusively found in Chinese individuals.
The seven Chinese samples heterozygous for the Arg25Gln allele were also
heterozygous for -411G > A, 3375C > T (exon 5, Arg257Cys), 7233T > G
(intron 8), and 7571G > C (3'-untranslated region), suggesting that
these mutations may be parts of a common haplotype. The potential linkage
between the Arg25Gln and the Arg257Cys alleles was also apparent when the
distribution of the Arg25Gln allele in the 102 newborn DNA samples described
above was compared with that of the Arg257Cys allele
(Table 3), which had been
determined previously in the same samples by PCR-SSCP
(Zhang et al., 2002
). The data
in Table 3 also indicate the
occurrence of a haplotype (found in one individual) that consists of the
25Gln, but not the 257Cys allele. The linkage between the Arg25Gln and the
Arg257Cys alleles was subsequently confirmed by sequencing subcloned long-PCR
products that correspond to single chromosomes. Results from four Chinese DNA
samples containing heterozygous mutations at both exon 1 and exon 5 SNP sites
consistently showed that the 74G > A and 3375C > T alleles were on the
same chromosome, as were the -411G > A, 7233T > G, and 7571G > C
alleles.
|
| Discussion |
|---|
|
|
|---|
The Arg257Cys mutation is known to cause a general reduction in CYP2A13
activity for all substrates tested, including NNK,
2'-methoxyacetophenone, coumarin, hexamethylphosphoramide,
N,N-dimethylaniline, and N-nitrosomethylphenylamine
(Zhang et al., 2002
), but it
is not yet known whether the Arg25Gln mutation has any functional
consequences. This mutation is unlikely to affect membrane retention or
protein synthesis, since both Arg and Gln appear in orthologous CYP2A proteins
(e.g., Arg in CYP2A13 and Gln in CYP2A6). Moreover, in CYP2C1, mutation of
Lys21 [which appears to correspond to Arg25 of CYP2A13 according to a sequence
alignment between CYP2A and CYP2C (Gotoh,
1992
)] to Asn did not affect membrane retention of a green
fluorescent protein reporter
(Szczesna-Skorupa and Kemper,
2000
). On the other hand, the effects of specific mutations in
this linker region on catalytic activity have not been examined. Although this
region is not part of the surface for interaction with NADPH-cytochrome P450
reductase or for substrate binding
(Williams et al., 2000
), and
although the Arg25Gln change is unlikely to affect gross structure (since
Gln25 occurs in functional CYP2As), this mutation does involve the loss of a
positive charge. It remains to be determined whether such a change could have
subtle effects on protein orientation and substrate access through the lipid
bilayer. It will also be important to determine the impact of the Arg25Gln and
Arg257Cys double mutations on CYP2A13 function since these mutations are
linked in most cases. However, the Asp158Glu mutation in exon 3 is a
conservative change and is less likely to have a significant impact on
function.
The Arg101Stop mutation represents a null allele, since the premature stop
codon will most likely lead to the synthesis of a truncated protein containing
only the amino-terminal 100 residues, which could be unstable and would not be
functional as a cytochrome P450 monooxygenase. The clear-cut functional
consequence of this mutation will make it particularly interesting for future
studies that correlate CYP2A13 genotype with respiratory tract
diseases. Individuals with the null allele may be at a reduced risk of
chemical toxicity from compounds metabolically activated by CYP2A13 in the
respiratory tract. Nonfunctional alleles have also been found in the
CYP2A6 gene (e.g., Nunoya et al.,
1998
; Oscarson et al.,
1999
), which has been shown to play a major role in the metabolic
disposition of nicotine and may thus influence the levels of exposure to
tobacco-related carcinogens, such as NNK, in active smokers (e.g.,
Oscarson, 2001
;
Xu et al., 2002
;
Yoshida et al., 2002
). In
fact, male smokers homozygous for the CYP2A6*1 allele seem
to have an elevated risk for tobacco-induced lung cancers
(Ariyoshi et al., 2002
). It is
likely that individuals defective in both CYP2A6 and CYP2A13
genes may be further protected against tobacco-related toxicity in the
lung.
Functional consequences of the noncoding region SNPs are more difficult to
predict. Of the two SNPs detected in the 5'-flanking region, neither is
located in a binding site for any known transcription factor, as determined by
a search of the TRANSFAC data base
(Heinemeyer et al., 1998
) using
the TFSEARCH program (Y. Akiyama: TFSEARCH: Searching Transcription Factor
Binding Sites,
http://www.crbc.jp/papia.html).
None of the intron mutations occur at a known splicing site
(Senapathy et al., 1990
) or
would generate a new splice site, as determined by an analysis using the Omiga
program, version 2.0 (Accelrys, Cambridge, UK). The two SNPs in the
3'-untranslated region are not expected to alter RNA folding, according
to an analysis using mfold, version 3.1
(http://www.bioinfo.rpi.edu/applications/mfold)
(Mathews et al., 1999
;
Zuker et al., 1999
). Further
studies on the impact of these mutations on the expression of the
CYP2A13 gene are warranted, since significant interindividual
differences in the level of CYP2A proteins have been found in microsomes from
fetal nasal mucosa (Gu et al.,
2000
).
Xiuling Zhang
Ying Chen
Yiqin Liu
Xiang Ren
Qing-Yu Zhang
Michele Caggana
Xinxin Ding
Wadsworth Center, New York State Department of Health, School of Public Health, State University of New York at Albany, New York (X.Z., Q-Y.Z., M.C., X.D.); and Qilu Hospital of Shandong University, Jinan, Shandong, China (Y.C., Y.L., X.R.)
| Acknowledgments |
|---|
| Footnotes |
|---|
1 Abbreviations used are: P450, cytochrome P450; NNK,
4-(methylnitrosamino)-1-(3-pyridyl)-1-butanone; PCR-SSCP, polymerase chain
reaction-single-strand conformation polymorphism; SNP, single nucleotide
polymorphism; bp, base pair(s). ![]()
Address correspondence to: Dr. Xinxin Ding, Wadsworth Center, New York State Department of Health, Empire State Plaza, Box 509, Albany, New York 12201-0509. E-mail: xding{at}wadsworth.org
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. D'Agostino, X. Zhang, H. Wu, G. Ling, S. Wang, Q.-Y. Zhang, F. Liu, and X. Ding Characterization of CYP2A13*2, a Variant Cytochrome P450 Allele Previously Found to Be Associated with Decreased Incidences of Lung Adenocarcinoma in Smokers Drug Metab. Dispos., November 1, 2008; 36(11): 2316 - 2323. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Zhang, J. D'Agostino, H. Wu, Q.-Y. Zhang, L. von Weymarn, S. E. Murphy, and X. Ding CYP2A13: Variable Expression and Role in Human Lung Microsomal Metabolic Activation of the Tobacco-Specific Carcinogen 4-(Methylnitrosamino)-1-(3-pyridyl)-1-butanone J. Pharmacol. Exp. Ther., November 1, 2007; 323(2): 570 - 578. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Misra, J.-Y. Hong, and S. Kim Multiplex Genotyping of Cytochrome P450 Single-Nucleotide Polymorphisms by Use of MALDI-TOF Mass Spectrometry Clin. Chem., May 1, 2007; 53(5): 933 - 939. [Abstract] [Full Text] [PDF] |
||||
![]() |
S.-L. Wang, X.-Y. He, J. Shen, J.-S. Wang, and J.-Y. Hong The Missense Genetic Polymorphisms of Human CYP2A13: Functional Significance in Carcinogen Activation and Identification of A Null Allelic Variant Toxicol. Sci., November 1, 2006; 94(1): 38 - 45. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. Hukkanen, P. Jacob III, and N. L. Benowitz Metabolism and Disposition Kinetics of Nicotine Pharmacol. Rev., March 1, 2005; 57(1): 79 - 115. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Zhang, M. Caggana, T. L. Cutler, and X. Ding Development of a Real-Time Polymerase Chain Reaction-Based Method for the Measurement of Relative Allelic Expression and Identification of CYP2A13 Alleles with Decreased Expression in Human Lung J. Pharmacol. Exp. Ther., October 1, 2004; 311(1): 373 - 381. [Abstract] [Full Text] [PDF] |
||||
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||