Evaluating NAT2PRED for inferring the individual acetylation status from unphased genotype data
© Sabbagh et al; licensee BioMed Central Ltd. 2009
Received: 11 August 2009
Accepted: 31 December 2009
Published: 31 December 2009
Genetically determined differences in N-acetylation capacity have proved to be important determinants of both the effectiveness of therapeutic response and the development of adverse drug reactions and toxicity during drug treatment. NAT2PRED is a web-server that allows a fast determination of NAT2 acetylation phenotype from genotype data without taking the extra step of reconstructing haplotypes for each individual (publicly available at http://nat2pred.rit.albany.edu). However, the classification accuracy of NAT2PRED needs to be assessed before its application can be advocated at a large scale.
The ability of NAT2PRED to classify individuals according to their acetylation status (slow, intermediate and rapid acetylators) was evaluated in a worldwide dataset composed of 56 population samples (8,489 individuals) from four continental regions.
NAT2PRED correctly identified slow acetylators with a sensitivity above 99% for all populations outside sub-Saharan Africa. Nevertheless, NAT2PRED showed a poor ability to distinguish between intermediate and rapid acetylators, with a classification error rate reaching up to 10% in the non-African samples.
NAT2PRED is an excellent tool to infer the individual acetylation status from NAT2 genotype data when the main interest is to distinguish slow acetylators from the others. This should facilitate the determination of the individual acetylation status in routine clinical practice and lead to better monitoring of risks associated with cancer and adverse drug reactions.
The human acetylation polymorphism is one of the most intensively studied pharmacogenetic traits that underlie interindividual and interethnic differences in response to xenobiotics. Genetically determined differences in N-acetylation capacity have proved to be important determinants of both the effectiveness of therapeutic response and the development of adverse drug reactions and toxicity during drug treatment [1, 2]. Some of the drugs excreted by acetylation are crucial in the treatment of diseases representing a worldwide concern, such as tuberculosis, AIDS-related complex diseases, and hypertension. Moreover, numerous association studies have linked the acetylation phenotype to susceptibility to a variety of complex human diseases, the most consistent findings being those regarding urinary bladder cancer, asthma and other allergic disorders [3–6].
Single nucleotide polymorphisms (SNPs) in the coding region of the NAT2 gene determine the acetylation phenotype. Thanks to the well-established genotype-phenotype correlation at this locus, the individual acetylation status can be reliably predicted from the haplotype combination at NAT2, according to the acknowledged classification of NAT2 haplotypes into either low-activity or fully functional alleles (see the consensus NAT2 gene nomenclature website: http://www.louisville.edu/medschool/pharmacology/NAT.html). Slow, intermediate and rapid acetylators are defined as carriers of zero, one or two functional haplotypes, respectively. Many studies did not distinguish between rapid and intermediate acetylators, categorizing both types of subjects as rapid acetylators.
Determining the linkage phase of SNP alleles along the coding sequence of NAT2 is crucial to unequivocally assign an individual's multi-site NAT2 genotype to a particular combination of two multilocus haplotypes and to correctly infer his acetylation status. Ambiguous NAT2 genotyping data that may lead to patient misclassification appear indeed to be common in human populations . Special analytical techniques have been designed to unambiguously determine the allocation of SNP alleles to either DNA strand. However, molecular haplotyping methods are labour-intensive and expensive and are not suitable for routine clinical applications. A cheap and straightforward alternative for haplotype reconstruction is the use of computational algorithms. Several studies have demonstrated the high effectiveness of haplotype reconstruction algorithms in the particular case of the NAT2 gene by comparing the computationally inferred haplotypes to the real ones resolved through the use of molecular haplotyping techniques [7–9]. But such in silico approaches may not be convenient to clinicians not familiar with computational haplotype analysis. Moreover, algorithmic techniques are statistical and require the analysis of a population rather than a single or a few individuals.
To circumvent these limitations, Kuznetsov et al.  have developed a web-server NAT2PRED that allows a fast determination of NAT2 acetylation phenotypes (slow, intermediate, and rapid) from the unphased genotype data at six polymorphic positions in NAT2. These six SNPs are the most commonly reported ones in surveys of NAT2 sequence variation in human populations . NAT2PRED alleviates the need of reconstructing haplotypes by implementing a supervised pattern recognition method that was trained on the NAT2 genotyping data of 1,377 subjects of known acetylation status. This tool is publicly available http://nat2pred.rit.albany.edu and has a simple intuitive user interface. It showed a nearly perfect classification accuracy of 99.9% in a sample mostly composed of Caucasians . However, it remains unclear to what extent this tool can be applied to individuals from any ethnicity since it was developed on a dataset where 94% of subjects were Caucasian. The accuracy of NAT2PRED needs to be assessed before its application can be advocated at a large scale.
The objective of the present study was to empirically evaluate the performance of NAT2PRED in a wide collection of population samples worldwide. To this end, we performed an extensive survey of the literature to identify those samples that were adequately genotyped for all the common SNPs in NAT2. In total the collected data consisted of 8,489 individuals from 56 human populations representing major geographic regions: sub-Saharan Africa (12 samples), Europe and North Africa (23), Central/South Asia (5), East Asia (13), and America (3). In each sample, NAT2 haplotypes were reconstructed using either molecular or computational methods. Therefore, these data provided an opportunity to compare, for each subject, the acetylation status inferred by NAT2PRED from the unphased NAT2 genotype data with the "true" one predicted from the pair of haplotypes reconstructed through molecular or computational haplotyping. NAT2PRED correctly identified slow acetylators with a sensitivity above 99% for all populations outside sub-Saharan Africa.
We selected from published reports up to November 2009 all the population samples that were genotyped for the seven most common SNPs in NAT2 (191G>A (rs1801279), 282C>T (rs1041983, rs59855457), 341T>C (rs1801280, rs56935242), 481C>T (rs1799929, rs60310310), 590G>A (rs1799930, rs60190029), 803A>G (rs1208, rs56599719, rs58999469), and 857G>A (rs1799931, rs58803786)) and for which genotype frequency data were available in the paper. We also included those non-African samples that were genotyped for all SNPs except SNP 191G>A as this is rare in non-African populations. The final dataset used for the present study describe 56 population samples from throughout the world, representing 8,489 individuals. A summary description of the data set is provided in Additional File 1.
In each sample, NAT2 haplotypes were either directly resolved using molecular-haplotyping techniques (through allele-specific PCR and restriction mapping) or were computationally inferred from the unphased multilocus genotypes using statistical algorithms (based either on a parsimony, maximum-likelihood, or Bayesian approach). For some samples, a combination of the two approaches was used. The acetylation phenotype of each individual was then predicted by assuming that the homozygous or compound heterozygous genotype for two haplotypes of the series NAT2*4, NAT2*11, NAT2*12 or NAT2*13 results in the rapid acetylator status, the occurrence of one of these haplotypes in combination with a haplotype of the series NAT2*5, NAT2*6, NAT2*7 or NAT2*14 results in the intermediate acetylator status and the occurrence of two haplotypes of the series NAT2*5, NAT2*6, NAT2*7 or NAT2*14 results in the slow acetylator phenotype.
Besides, NAT2PRED was applied to the unphased genotype data of each sample (at the six polymorphic positions 282, 341, 481, 590, 803, and 857) to provide an inferred acetylation status that was compared to the "true" one predicted from the diplotype configuration at NAT2. The web-server was used in the standard way, by specifying the genotype at each polymorphic position using radio buttons and clicking the 'Submit' button. We retained for each subject the final prediction (i.e., the phenotype with the highest probability). The performance of NAT2PRED was measured by computing in each sample the classification error rate, which is the proportion of individuals whose acetylation status was incorrectly inferred by the web-server. The classification error rate of NAT2PRED was evaluated by considering either two (slow vs. others) or three (slow vs. intermediate vs. rapid) phenotypic classes for the acetylation status.
Results and discussion
Classification error rates of NAT2PRED in each geographic regiona.
Classification error rate (%)
Two acetylation phenotypes b
Three acetylation phenotypes c
Sub-Saharan Africa (n = 617)
Europe and North Africa (n = 4,391)
Central/South Asia (n = 757)
East Asia (n = 2,346)
America (n = 378)
It is important to note that for most of the samples investigated, NAT2 haplotypes were not actually determined by a molecular resolution of mutation linkage patterns but were rather inferred through computational methods (see Additional File 1). As real NAT2 haplotypes may be quite different from those calculated by computational inference , it is of high relevance to evaluate the performance of NAT2PRED in only those samples where the real NAT2 haplotypes were determined. There are only six population samples where molecular techniques were used to reconstruct haplotypes and applied to the whole set of individuals. In the single African sample (Tswana from South Africa), NAT2PRED showed a poor classification accuracy with a classification error rate of 10% and 16% when considering two (slow vs. others) or three (slow vs. intermediate vs. rapid) phenotypic classes for the acetylation status, respectively. In the five remaining non-African samples, the ability of NAT2PRED to distinguish slow acetylators from the others was very high with an average classification error rate of only 0.5%, ranging from 0% in UK Caucasians and Nicaraguans to 1.2% in Dravidian Indians. This means that in 99.5% of cases, NAT2PRED correctly identified slow acetylators in the same way the usual procedure did (that is, by reconstructing the real haplotypes through molecular techniques and by inferring the individual acetylation status from the haplotype combination at NAT2). On the other hand, a poorer performance was observed when intermediate acetylators were to be distinguished from rapid acetylators in these five samples (average error rate of 2.6%, ranging from 0% to 10.2%).
Our results show that NAT2PRED is an excellent tool to infer the individual acetylation status from NAT2 genotype data when the main interest is to distinguish slow acetylators from the others. Its classification error rate is indeed too high to ensure a reliable discrimination between intermediate and rapid acetylators. However, since the slow acetylator phenotype is the one which exhibits today the strongest association with susceptibility to bladder cancer and adverse drug reactions [2, 3], the web-server NAT2PRED should be of great use in clinical laboratories. This should facilitate the determination of the individual acetylation status in routine clinical practice and lead to better monitoring of risks associated with cancer and adverse drug reactions. On the other hand, NAT2PRED must be used with caution for patients whose population of origin contains a sizeable frequency of SNP 191G>A. This concerns populations of African origin but also several non-African populations, notably South-European and South-American populations, where this SNP has been identified, probably as a consequence of genetic admixture with Africans . In particular, a frequency above 1% has been reported for the 191A allele in Spanish (1.4%) , in Coyaima from Colombia (3.5%) , in Brazilians from Rio de Janeiro (4.7%) , and in Brazilians from Goiàs (2.4%) .
- Meisel P: Arylamine N-acetyltransferases and drug response. Pharmacogenomics. 2002, 3 (3): 349-366. 10.1517/14622418.104.22.1689.View ArticlePubMedGoogle Scholar
- Ladero JM: Influence of polymorphic N-acetyltransferases on non-malignant spontaneous disorders and on response to drugs. Curr Drug Metab. 2008, 9 (6): 532-537. 10.2174/138920008784892038.View ArticlePubMedGoogle Scholar
- Agúndez JA: Polymorphisms of human N-acetyltransferases and cancer risk. Curr Drug Metab. 2008, 9 (6): 520-531. 10.2174/138920008784892083.View ArticlePubMedGoogle Scholar
- García-Closas M, Malats N, Silverman D, Dosemeci M, Kogevinas M, Hein DW, Tardón A, Serra C, Carrato A, García-Closas R, Lloreta J, Castaño-Vinyals G, Yeager M, Welch R, Chanock S, Chatterjee N, Wacholder S, Samanic C, Torà M, Fernández F, Real FX, Rothman N: NAT2 slow acetylation, GSTM1 null genotype, and risk of bladder cancer: results from the Spanish Bladder Cancer Study and meta-analyses. Lancet. 2005, 366 (9486): 649-659. 10.1016/S0140-6736(05)67137-1.View ArticlePubMedPubMed CentralGoogle Scholar
- Hein DW: N-acetyltransferase 2 genetic polymorphism: effects of carcinogen and haplotype on urinary bladder cancer risk. Oncogene. 2006, 25 (11): 1649-1658. 10.1038/sj.onc.1209374.View ArticlePubMedPubMed CentralGoogle Scholar
- Batra J, Ghosh B: N-acetyltransferases as markers for asthma and allergic/atopic disorders. Curr Drug Metab. 2008, 9 (6): 546-553. 10.2174/138920008784892074.View ArticlePubMedGoogle Scholar
- Agúndez JA, Golka K, Martínez C, Selinski S, Blaszkewicz M, García-Martín E: Unraveling ambiguous NAT2 genotyping data. Clin Chem. 2008, 54 (8): 1390-1394. 10.1373/clinchem.2008.105569.View ArticlePubMedGoogle Scholar
- Sabbagh A, Darlu P: Inferring haplotypes at the NAT2 locus: the computational approach. BMC Genet. 2005, 6 (1): 30-10.1186/1471-2156-6-30.View ArticlePubMedPubMed CentralGoogle Scholar
- Golka K, Blaszkewicz M, Samimi M, Bolt HM, Selinski S: Reconstruction of N-acetyltransferase 2 haplotypes using PHASE. Arch Toxicol. 2008, 82 (4): 265-270. 10.1007/s00204-007-0247-9.View ArticlePubMedGoogle Scholar
- Kuznetsov IB, McDuffie M, Moslehi R: A web-server for inferring the human N-acetyltransferase-2 (NAT2) enzymatic phenotype from NAT2 genotype. Bioinformatics. 2009, 25 (9): 1185-1186. 10.1093/bioinformatics/btp121.View ArticlePubMedPubMed CentralGoogle Scholar
- Sabbagh A, Langaney A, Darlu P, Gérard N, Krishnamoorthy R, Poloni ES: Worldwide distribution of NAT2 diversity: implications for NAT2 evolutionary history. BMC Genet. 2008, 9: 21-10.1186/1471-2156-9-21.View ArticlePubMedPubMed CentralGoogle Scholar
- Patin E, Harmant C, Kidd KK, Kidd J, Froment A, Mehdi SQ, Sica L, Heyer E, Quintana-Murci L: Sub-Saharan African coding sequence variation and haplotype diversity at the NAT2 gene. Hum Mutat. 2006, 27 (7): 720-10.1002/humu.9438.View ArticlePubMedGoogle Scholar
- García-Martín E: Interethnic and intraethnic variability of NAT2 single nucleotide polymorphisms. Curr Drug Metab. 2008, 9 (6): 487-497. 10.2174/138920008784892155.View ArticlePubMedGoogle Scholar
- Fuselli S, Gilman RH, Chanock SJ, Bonatto SL, De Stefano G, Evans CA, Labuda D, Luiselli D, Salzano FM, Soto G, Vallejo G, Sajantila A, Pettener D, Tarazona-Santos E: Analysis of nucleotide diversity of NAT2 coding region reveals homogeneity across Native American populations and high intra-population diversity. Pharmacogenomics J. 2007, 7 (2): 144-152. 10.1038/sj.tpj.6500407.View ArticlePubMedGoogle Scholar
- Teixeira RL, Miranda AB, Pacheco AG, Lopes MQ, Fonseca-Costa J, Rabahi MF, Melo HM, Kritski AL, Mello FC, Suffys PN, Santos AR: Genetic profile of the arylamine N-acetyltransferase 2 coding gene among individuals from two different regions of Brazil. Mutat Res. 2007, 624 (1-2): 31-40.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2350/10/148/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.