{'fraction_non_alphanumeric': 0.045705928058869236, 'fraction_numerical': 0.05343022990081814, 'mean_word_length': 5.253215204344098, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 1, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Differences in UGT1A1, UGT1A7, and UGT1A9 Polymorphisms between Uzbek and Japanese Populations

Background and Objectives Uridine-diphosphate glucuronosyltransferase 1A (UGT1A) is a key enzyme involved in irinotecan metabolism, and polymorphisms in the UGT1A gene are associated with irinotecan-induced toxicity. The aim of this study was to elucidate the allele frequencies of UGT1A polymorphisms in healthy Uzbek volunteers, and to compare them with those of the Japanese population. Method A total of 97 healthy volunteers from Uzbekistan were enrolled and blood samples were collected from each participant. Genotyping analysis was performed by fragment size analysis for UGT1A1*28, direct sequencing for UGT1A7*3 and UGT1A9*22, and TaqMan assays for UGT1A1*93, UGT1A1*6, UGT1A1*27, UGT1A1*60, and UGT1A7*12. The frequencies of polymorphisms were compared with the Japanese population by using the data previously reported from our study group. Results When the Uzbek and Japanese populations were compared, heterozygotes or homozygotes for UGT1A1*28, UGT1A1*60, and UGT1A1*93 were significantly more frequent in the Uzbek population (P < 0.01). The rate of UGT1A7*12 was not significantly different between the two populations, whereas UGT1A1*6 and UGT1A9*22 were significantly less frequent in the Uzbek population (P < 0.05). UGT1A7*1 were less prevalent in the Uzbek population than in the Japanese population (P < 0.01). Conclusion The Uzbek population has different frequencies of polymorphisms in UGT1A genes compared with the Japanese population. A comprehensive study of the influence of UGT1A1 polymorphisms on the risk of irinotecan-induced toxicity is necessary for optimal use of irinotecan treatment.


Introduction
Irinotecan in combination with fluorouracil and leucovorin prolongs the progression-free survival and overall survival of patients with metastatic colorectal cancer compared with fluorouracil and leucovorin alone. However, grade 3/4 diarrhea and neutropenia can occur during the course of the treatment, jeopardizing approximately 20-25 and 25-55 % of patients, respectively [1,2]. The toxicity of irinotecan is influenced by several non-genetic factors [3,4]. A number of recent studies have found a strong relationship between the adverse effects of irinotecan-based treatment and polymorphisms of the uridine-diphosphate glucuronosyltransferase 1A (UGT1A) gene, which encodes a crucial enzyme in irinotecan metabolism [5][6][7][8][9][10][11][12][13].
Administered irinotecan is activated and transformed through hydrolysis into 7-ethyl-10-hydroxycamptothecin (SN-38), which is a potent inhibitor of topoisomerase I and is the metabolite responsible for adverse reactions to irinotecan-based treatment [14]. SN-38 is then converted into an inactive form, 10-O-glucuronyl-SN-38 (SN-38G), by various UGT1A isoforms including UGT1A1, 1A6, 1A7, 1A8, 1A9, and 1A10 [15][16][17]. Therefore, irinotecaninduced toxicities could be due to a reduced glucuronidation rate of SN-38 caused by genetic variants of UGT1A. One of the representative polymorphisms that strongly influences irinotecan clearance is UGT1A1*28, which has (TA) 7 in the promoter region instead of (TA) 6 as in UGT1A1*1 [18]. Individuals with UGT1A1*28 demonstrate reduced irinotecan clearance [10,12] due to reduced transcriptional activity of UGT1A1 [18]. The UGT1A1*28 polymorphism may increase the clinical risk of severe irinotecan-induced toxicity in colorectal cancer patients [5,6]. In contrast to the effect of UGT1A1*28 on transcriptional efficacy, UGT1A1*6 has G replaced by A at position ?211 relative to the UGT1A1 transcription start site, which results in decreased irinotecan metabolism [19].
Several polymorphisms, especially UGT1A1, 1A7, and 1A9, are associated with an alteration in irinotecan metabolism and are receiving increasing attention in clinical settings [8,20], but only UGT1A1*28 and UGT1A1*6 have been extensively studied. We recently examined irinotecan metabolism-related UGT1A polymorphisms and reported their frequency in three geographically different regions in Japan [21]. We consider that studying the interethnic diversity of these polymorphisms will lead to an understanding of the inconsistent drug responses observed between different populations, and will ultimately help establish the best strategies for application of new drugs in these populations. In this context, this study was designed to determine the allele frequency of UGT1A1*27, UGT1A1*60, UGT1A1*93, UGT1A7*2 (N129K), UGT1A7*3 (co-occurrence of N129K and W208R), UGT1A7*12, and UGT1A9* 22 as well as UGT1A1*28 and UGT1A1*6 in healthy volunteers of the Republic of Uzbekistan, a country located between East Asia and Europe (Table 1; Fig. 1).

Subjects
This study was approved by the ethics committee of Yamaguchi University, and written informed consent was obtained from all participants. We recruited 97 volunteers from Uzbekistan. A 7 ml sample of peripheral blood was collected from each individual and stored in ethylenediaminetetraacetic acid (EDTA). The 150 healthy Japanese volunteers were from three districts in Japan, and we have investigated and reported on their polymorphisms previously [21].

Genotyping
Genomic DNA from peripheral blood anti-coagulated with EDTA was extracted by using the conventional NaI method. UGT1A1*28 was identified using fragment size analysis, and UGT1A7*3 and UGT1A9*22 were detected using direct sequencing. A TaqMan assay was performed for the identification of UGT1A1*93, UGT1A1*6, UGT1A1*27, UGT1A1*60, and UGT1A7*12. Primers and probes used in the study are listed in Table 1.
For fragment size analysis, polymerase chain reaction (PCR) reactions using 800 ng of template DNA were performed according to the manufacturer's instructions (Ex Taq; Takara, Tokyo, Japan). The PCR products of TA6 (6 TA repeats; 94 bp) and TA7 (7 TA repeats; 96 bp) were mixed with Hi-Di formamide and the internal size standard (GeneScan LIZ-500, Applied Biosystems, CA, USA) at a ratio of 1:10 (v/v). After electrophoresis in the ABI Prism 3100 Genetic Analyzer (Applied Biosystems), fragment sizes were compared with the internal size standard and determined by the local Southern algorithm and by GeneMapper software version 3.5 (Applied Biosystems).
For direct sequencing, PCR was conducted using the Gene Amp PCR System PC808 (ASTEC, Tokyo, Japan) with Ex Taq polymerase. Purification of PCR products was performed with ExoSAP-IT (Amersham Bioscience, Tokyo, Japan) for 20 min at 37°C followed by 20 min at 80°C. Sequencing reactions were conducted with a BigDye Terminator Cycle Sequencing Kit (Applied Biosystems, Tokyo, Japan), and, after the reaction products were purified with ethanol, analysis was performed using an ABI 3100-Avant Genetic Analyzer (Applied Biosystems). Table 1, the transposition of T to G at position 387 (or N129K), C to A at position 391, and G to A at position 392 are termed UGT1A7*2. These three polymorphisms are genetically linked, and only N129K is described in this manuscript. The transposition of T to C at position 622 is termed W208R, and co-occurrence of N129K and W208R is termed UGT1A7*3, whereas the occurrence of W208 alone is termed UGT1A7*4. Thus, in the case of heterozygous N129K and heterozygous W208R in the same individual, the possible genotypes are UGT1A7*1/*3 and *2/*4. However, because UGT1A7*4 is known to be quite rare in non-Uzbek populations [9,17], the statistical analysis was performed on the assumption that all the genotypes were UGT1A7*1/*3.

As described in
For TaqMan assays for UGT1A1*93, forward and reverse PCR primers and TaqMan probes were customsynthesized by Applied Biosystems. Primers and probes for UGT1A1*6, UGT1A1*27, UGT1A1*60, and UGT1A7*12 were obtained from Applied Biosystems (TaqMan SNP Genotyping Assays). Reaction mixtures were loaded into 384-well plates and analyzed in an ABI Prism 7900HT Sequence Detection System (Applied Biosystems).

Statistical Analysis
Each nucleotide polymorphism was evaluated for the Hardy-Weinberg equilibrium, and by linkage disequilibrium (LD) analysis using Haploview 4.2 software (Massachusetts Institute of Technology, Cambridge, MA, USA; http://www. broadinstitute.org/scientific-community/science/programs/ medical-and-population-genetics/haploview/haploview). The Lewontin's coefficient D 0 and correlation coefficient, r 2 , were calculated as measures of LD. The proportions of subjects with homozygous wild-type alleles or heterozygous or homozygous variant alleles were calculated with 95 % Agresti-Coull confidence intervals (95 % CI).
Nucleotide positions are relative to the transcription start site Transposition of T to G at position 387 (or N129K) is termed UGT1A7*2. Transposition of T to C at position 622 is termed W208R, and cooccurrence of N129K and W208R is termed UGT1A7*3 TA6 6 TA repeats, TA7 7 TA repeats, F-FAM forward primer-labeled reporter 1 probe, F forward primer, R reverse primer, FAM reporter 1 probe, VIC reporter 2 probe, WT wild type Fisher's exact test with a two-sided significance level of 0.05 was used for comparing the areas. For a two-sided 95 % CI for a binomial proportion whose true value is varied from 0.5 to 0.1, a sample size of 50 yields a halfwidth of, at most, 14 % in any situations of the true value. These analyses were performed using SAS v9.3 (SAS Institute Inc., Cary, NC, USA).

Baseline Characteristics, Hardy-Weinberg Equilibrium, and Linkage Disequilibrium
The average age of the participants was 41.1 (range 19-78) years in Uzbekistan, and 39.9 (range 18-67) years in Japan. The 97 volunteers from Uzbekistan comprised 37 female and 60 males, while 125 of the 150 Japanese volunteers were female. As described previously [21], the Japanese volunteers were mainly recruited from nurses, which is an occupation traditionally chosen by women in Japan.
All of the UGT1A polymorphisms were in Hardy-Weinberg equilibrium (P [ 0.05). Perfect linkages were observed between W208 and UGT1A7*12 in both the Uzbek and the Japanese populations (data not shown). UGT1A1*28 and UGT1A1*93 showed very strong linkage in the Uzbek population (r 2 = 0.95, D 0 = 1.00) and they matched completely in the Japanese population (Fig. 2). Close linkage was also observed between N129K and UGT1A9*22, both in the Uzbek (r 2 = 0.94, D 0 = 1.00), and in the Japanese (r 2 = 0.92, D 0 = 1.00) populations. Table 2 lists the frequencies of UGT1A1 polymorphisms, UGT1A1*6, *28, *60, *27, and *93. The two polymorphisms that have been extensively studied for their role in irinotecan metabolism, UGT1A1*28 and UGT1A1*6, showed different incidences in the Uzbek and Japanese populations. UGT1A1*28 showed a significantly higher incidence and UGT1A1*6 showed a significantly lower incidence in the Uzbek and the Japanese populations (P \ 0.01, P = 0.012, respectively). Neither UGT1A1*36 (5TA repeat in the promoter region) nor UGT1A1*37 (8TA repeat in the promoter region) were identified among the Uzbek volunteers. Among the other UGT1A1 polymorphisms examined, UGT1A*27 was very rare in both Uzbek and Japanese populations, occurring at a rate of 0.005 and 0.003, respectively, and UGT1A1*60 and UGT1A1*93 were significantly more prevalent in the Uzbek population than in the Japanese population (both P \ 0.01).

Frequencies of UGT1A7 and UGT1A9
Polymorphisms Table 3 lists the frequencies of UGT1A7*3 (N129K and W208R), UGT1A7*12, and UGT1A9*22. The Uzbek population had a significantly higher incidence of homozygous N129K than did the Japanese population (P \ 0.01). However, the incidences of W208R and UGT1A7*12 were not significantly different between the two populations; furthermore, W208R and UGT1A7*12 showed evidence of genetic linkage, which is reflected in their identical values in Table 3. No Uzbek individual carrying heterozygous UGT1A7*1/*4 or *3/*4, or homozygous UGT1A7*4 was detected. UGT1A7*1 is less prevalent in Uzbekistan than in Japan (P \ 0.01). UGT1A9*22 was significantly less prevalent in the Uzbek population than in the Japanese population (P \ 0.01).

Discussion
UGT1A is an essential enzyme for the elimination of numerous endogenous and exogenous compounds through glucuronidation [22]. Here, we analyzed the allele frequencies of several UGT1A polymorphisms in 97 healthy Uzbek volunteers and compared them with those of the Japanese population. The polymorphism UGT1A1*28 ([TA] 7 in the UGT1A1 promoter region) has been extensively investigated [23,24], particularly in terms of inter-individual variability of irinotecan-induced gastrointestinal and hematological toxicity [5][6][7][8][9][10][11][12]. Because inactivation of the active form of irinotecan (SN-38) protects against severe irinotecan-related toxicity, decreased transcriptional activity of UGT1A1 due to UGT1A1*28 is considered to be a risk factor for adverse reactions to irinotecan treatment. Here, we showed that the incidence of UGT1A1*28 homozygotes in the Uzbek population was 9.3 %, which was significantly higher than that in the Japanese population that we studied previously [21]. In published reports, the incidence of UGT1A1*28 homozygotes is high in Europe (5-14.8 %), Africa (5.9-17.9 %), and India (12.8-19.3 %); and less frequent in East Asia and Japan (0-5.9 %) [5][6][7][8][25][26][27][28][29][30][31]; the allele frequency of UGT1A1*28 shows a similar trend to the incidence of homozygotes (Table 4). Horsfall et al. [23] and Premawardhena et al. [24] reported further diversity in the numbers of TA repeats (i.e. [TA] 5 termed UGT1A1*36, and [TA] 8 termed UGT1A1*37) among individuals from Africa and those with varying degrees of African ancestry in North and Central America. UGT1A1*36 and UGT1A1*37 are also present in Caucasian and Indian populations at a very low rate [24,30]; however, to our knowledge, they have not been detected in Japan. Neither UGT1A1*36 nor UGT1A1*37 were identified among the 97 Uzbek volunteers.
The polymorphism UGT1A1*6 has also been extensively studied and has been shown to be associated with decreased irinotecan metabolism [19]. The current study showed the rates of heterozygotes and homozygotes for UGT1A1*6 in Uzbekistan were 13.4 and 2.1 %, respectively, which was significantly lower than the equivalent rates in the Japanese population we studied previously [21]. In a previous study, only two heterozygotes among 150 Caucasians and no individuals among 150 African-Americans were found to carry UGT1A1*6 [25]. The incidence of UGT1A1*6 polymorphisms in the Uzbek population was more frequent that that previously reported in Caucasian and African populations ( Table 4).
The variant 686C[A (UGT1A1*27) is a very rare polymorphism that is associated with a significant reduction in UGT1A activity in vitro [19]. No UGT1A1*27 allele was detected in 150 Caucasian or 150 African-American healthy volunteers [25], or in more than 300 Indian healthy volunteers [30,31] (Table 4). In contrast, studies from Japan identified heterozygous UGT1A1*27 in 2.6 % of Japanese colorectal cancer patients [5], in 0.7 % of Japanese or African-American healthy volunteers [21,25], and in 2.3 % of a mixed population of Japanese arrhythmic patients and cancer patients [26]. The current study demonstrated that the incidence of UGT1A1*27 heterozygotes was also low in the Uzbek population (1 of 97 individuals). UGT1A1*60 and UGT1A*93 are the polymorphisms that result in reduced UGT1A1 enzyme efficacy and may be related to the hematological toxicity observed after irinotecan treatment [8,32]. The allele frequency of UGT1A1*60 in the Uzbek population (0.5) was higher than that in various reports of the Japanese population (0.23-0.26), including our previous study. The frequency was similar to that reported for Caucasian (0.44-0.62) and Indian (0.43) populations, and lower than that reported for the African population (0.85) ( Table 4). Heterozygous or homozygous UGT1A1*93 was observed significantly more frequently in the Uzbek population than in the Japanese population that we studied previously [21]. In a small study of the Chinese population [27], UGT1A1*93 was detected at a frequency of 0.12, which is similar to that reported in various studies of the Japanese population, including our previous study. In contrast, UGT1A1*93 is relatively common in Europe, Africa, India, and Uzbekistan.
UGT1A9*22 is characterized by a one-base thymidine insertion in the promoter region of UGT1A9, which increases the UGT1A9 levels and causes a higher clearance rate of SN-38G than that observed for UGT1A9*1 [33]. Our previous study [21] revealed that the allele frequency of UGT1A9*22 in healthy Japanese volunteers was 0.66, which is comparable to that in Japanese colorectal cancer patients (0.54-0.66) [9,33,34]. A study of colorectal cancer patients who were mainly of Caucasian origin Heterozygous and homozygous UGT1A1*28, *60, and *93 were observed significantly more frequently in the Uzbek population than in the Japanese population. In contrast, heterozygous and homozygous UGT1A1*6 were significantly less prevalent in the Uzbek population than in the Japanese population. The occurrence of UTG1A1*27 was rare in both populations reported a UGT1A9*22 allele frequency of 0.4 [8,35], which is similar to that observed here for the healthy Uzbek volunteers.
The relationship between UGT1A polymorphisms and irinotecan-induced toxicities remains controversial. For instance, several early studies reported a strong association between UGT1A1*28 with severe (grade 3 or 4) toxicity in Japanese [5,6,11] and Caucasian colorectal patients [10]; however, another study found that UGT1A1*28 had little or no effect on the rate of adverse effects in such patients For each genotype, the values represent the number of subjects (percentage of subjects, 95 % CI). The occurrence of N129 alone is termed UGT1A7*2, and the occurrence of W208 alone is termed UGT1A7*4, whereas co-occurrence of N129K and W208R is termed UGT1A7*3. In the case of heterozygous N129K and heterozygous W208R, the possible genotypes are UGT1A7*1/*3 and *2/*4; however, because UGT1A7*4 is known to be quite rare in non-Uzbek populations, the statistical analysis was performed on the assumption that all the genotypes were UGT1A7*1/*3 [7,8,27]. Inconsistent patient responses to irinotecan imply that factors such as age, organ function, concomitant therapy [3,4], and dose of irinotecan [7] might modulate the therapeutic and adverse effects of this drug. However, the variable responses might also be explained by the synergic effect of UGT1A1*28 and other UGT1A polymorphisms, such as UGT1A1*6 and UGT1A7*3 [12,13,20,28]. Our finding that there is a significantly different distribution of UGT1A polymorphisms in Uzbekistan and Japan strengthens this hypothesis, and provides further evidence that extensive studies of the relationships between UGT1A haplotypes and the risk of irinotecan toxicity would be worthwhile. One limitation of this study was that the number of the Uzbek healthy volunteers examined was small. Thus, the rate of detection of rare polymorphisms, such as UGT1A1*36, UGT1A1*37, and UGT1A1*27, might not be representative of the whole Uzbek population. However, for other polymorphisms, the numbers of volunteers was sufficient to show very strong statistical differences between the Uzbek population and the Japanese population. Another limitation of the study was that a statistical comparison of the frequency of UGT1A polymorphisms in the Uzbek population and other populations that were studied using different genotyping methods was not possible. Our ongoing study of various ethnic groups from Asia to Europe will allow such direct comparisons and may clarify the movement of people throughout history [21].
LD analysis in the present study revealed strong linkage between some UGT1A genotypes. Recently, strong linkage between N129K and UGT1A9*22, as we have observed in Uzbek and Japanese populations, was also reported in a Caucasian population [8]. In contrast, the linkage between UGT1A1*28 and UGT1A1*93 was stronger in both Uzbek and Japanese populations (r 2 = 0.94 and 0.92, respectively) than in a Caucasian population (r 2 = 0.6782) [8]. This suggests that the difference in haplotype frequencies between different populations reflects both genotype (polymorphism) frequencies and the strength of LD. It is also plausible that the different strength of linkage is one of the reasons for the inconsistency observed between the adverse effects of irinotecan and the presence of UGT1A polymorphisms in treated patients. Thus, further investigation is required of UGT1A LD and haplotypes.
Our results revealed that UGT1A1*28 is more frequent in the Uzbek population than in the Japanese population, and is present at a similar rate to that reported for Caucasian and African populations. In contrast, UGT1A1*6 is less prevalent in the Uzbek population than in the Japanese population, but has a higher incidence than previous reports in Caucasians. Although UGT1A1*27 has not been reported in African and Caucasian populations, it was identified in a single Uzbek volunteer. UGT1A1*60 and *93 were detected more frequently in the Uzbek population than in the Japanese population, and the rate was similar to that reported for Caucasians. Finally, UGT1A9*22 was identified in the Uzbek population at a lower rate than that in the Japanese population. These findings confirm the inter-ethnic diversity of polymorphisms in UGT1A and lead us to speculate that the Uzbekistan population has genetic characteristic between those of East Asia and European countries, consistent with its geographical location (Table 4). We anticipate that further comprehensive studies revealing the relationship between UGT1A1 polymorphisms and the risk of irinotecan-induced toxicity will facilitate individual treatment option decisions.