New genetic biomarkers predicting azathioprine blood concentrations in combination therapy with 5-aminosalicylic acid.

BACKGROUND AND AIMS: Azathioprine (AZA) is widely used for the treatment of inflammatory bowel disease (IBD) patients. AZA is catabolized by thiopurine S-methyltransferase (TPMT), which exhibits genetic polymorphisms. It has also been reported that 5-aminosalicylic acid (5-ASA) inhibits TPMT activity, and that increased 6-thioguanine nucleotide (6-TGN, a metabolite of AZA) blood concentrations result in an increased number of ADRs. In this study, single nucleotide polymorphisms (SNPs) related to differential gene expression affecting AZA drug metabolism in combination therapy with 5-ASA were examined.
METHODS: To identify genetic biomarkers for the prediction of 6-TGN blood concentration, ExpressGenotyping analysis was used. ExpressGenotyping analysis is able to detect critical pharmacogenetic SNPs by analyzing drug-induced expression allelic imbalance (EAI) of premature RNA in HapMap lymphocytes. We collected blood samples on 38 patients with inflammatory bowel disease treated with AZA and corroboration of the obtained SNPs was attempted in clinical samples.
RESULTS: A large number of SNPs with AZA/5-ASA-induced EAI within the investigated HapMap lymphocytes was identified by ExpressGenotyping analysis. The respective SNPs were analyzed in IBD patients' blood samples. Among these SNPs, several that have not yet been described to be induced by AZA/5-ASA were found. SNPs within SLC38A9 showed a particular correlation with patients' 6-TGN blood concentrations.
CONCLUSIONS: Based on these results, ExpressGenotyping analysis and genotyping of patients appears to be a useful way to identify inter-individual differences in drug responses and ADRs to AZA/5-ASA. This study provides helpful information on genetic biomarkers for optimized AZA/5-ASA treatment of IBD patients.

Introduction

The immune-modulating thiopurines, 6-mercaptopurine (6-MP) and its pro-drug azathioprine (AZA), are widely used in inflammatory bowel disease (IBD) treatment [1]–[3]. AZA is a pro-drug of 6-MP carrying a methyl imidazole group attached to a sulfur carbon of 6-MP and metabolized to 6-MP by cleavage require the action of gluthatione S transferase. For the maintenance of remission in IBD, AZA and 6-MP are the first-line immunomodulators [4], [5]. Both AZA and 6-MP must be extensively metabolized before the pharmacologically active metabolites, 6-thioguaninenucleotides (6-TGN), are generated. The clinical efficacy of 6-MP has been reported to be correlated with erythrocyte levels of 6-TGN [6], [7]. The usefulness of monitoring 6-TGN blood concentrations in patients receiving AZA and 6-MP has also been suggested [8]–[10]. Although these thiopurine derivatives are considered to be relatively safe for maintenance therapy, several studies have reported discontinuation of thiopurine derivatives in up to 50% of patients during long-term therapy, mainly due to the development of adverse drug reactions (ADRs) [4], [11], [12]. AZA is catabolized by thiopurine S-methyltransferase (TPMT), which exhibits single nucleotide polymorphisms (SNPs). These SNPs have been reported to contribute to bone marrow suppression as ADRs [9], [13], [14]. Patients with mutant alleles are reported to exhibit lower enzymatic activity of TPMT and to be more likely to show bone marrow suppression. In the Japanese population, TPMT A719G is the most common allele [15]–[18] and is associated with low TPMT enzyme activity [19], [20]. However, in the Japanese population, the frequency of TPMT mutations is reported to be lower (2%) than that reported from Western countries [15]–[17], [21], [22]. It has also been reported that 5-aminosalicylic acid (5-ASA) inhibits TPMT activity, and that an increased 6-TGN blood concentration results in an increased number of ADRs [3]. However, the mechanism of ADRs induced by AZA has not yet been completely understood. Therefore, to identify new genetic biomarkers to predict 6-TGN blood concentrations before the administration of thiopurine is very valuable and important.

In this study, SNPs related to differential gene expression affecting AZA drug metabolism in AZA/5-ASA combination therapy were investigated. This genetic information is important to predict the 6-TGN blood concentrations of IBD patients treated with AZA/5-ASA. Recently, ExpressGenotyping analysis has been developed to explore the SNPs affecting gene expression changes of each allele (these gene expression changes are called Expression Allelic Imbalance (EAI)) after drug exposure (Figure 1). It is clear from a previous report that the concept of allelic variation such as this EAI is also important [23]–[25]. Such SNPs may be associated with inter-individual differences in drug response and ADRs. Therefore, ExpressGenotyping analysis can be used as a novel technology to detect critical genetic biomarkers that predict 6-TGN blood concentrations of IBD patients.

Figure 1

ExpressGenotyping analysis.

ExpressGenotyping analysis can detect critical pharmacogenetic SNPs by analyzing drug-induced expression allelic imbalance (EAI) of premature RNA in HapMap lymphocytes. This new analysis technique has the potential to accelerate the realization of personalized medicine, which is closely associated with inter-individual differences in drug response and ADRs.

Materials and Methods

Patients

Patients with established IBD receiving AZA were enrolled between April 1ST 2009 and August 31ST 2009. All patients were Japanese persons attending the gastroenterology outpatient clinic at the Hospital of the Kyoto Prefectural University of Medicine. The characteristics of the enrolled patients are shown in Table 1. Of the 38 patients, 35 were also treated with 5-ASA. Clinical information related to drug administration for these patients is shown in Table 2. The protocol of this study was approved by the Ethics Committee of the Kyoto Prefectural University of Medicine and written informed consent was obtained from all patients before the original procedures were conducted. This study was performed in accordance with the ethical principles associated with the Declaration of Helsinki and registered in the University Hospital Medical Network Clinical Trials Registry (UMIN 000010197).

Table 1

Characteristics of enrolled patients.

Total (n)	38
Sex (n; male/female)	29/9
Age (years; median±range)	44.0±10.6
Disease (n; UC/CD/BD)	21/14/3
Location of UC patients (n = 21)
Pancolitis	19
Distal colitis	2
Location of CD patients (n = 14)
Colon	4
Colon and ileum	6
Ileum	4
Duration of IBD (years; median±range)	12.8±8.97

UC: ulcerative colitis, CD: Crohn's disease, BD: Behcet's disease.

Table 2

Clinical information of drug administration for patients.

Sample number	Catalog ID	Dose of azathioprine (mg/day)	Dose of 5-ASA (mg/day)	Concentration of 6-TGN (pmol/8×108 RBCs)	Risk ratios (6-TGN/AZA)
1	5842535	75	3600	195	2.60
2	5779006	50	2250	474	9.48
3	770916	25	1500	423	16.92
4	4872959	25	3000	130	5.20
5	5459355	50	3000	287	5.74
6	6232762	50	2250	243	4.86
7	2639084	50	2250	164	3.28
8	0156158	50	1500	144	2.88
9	5571916	50	3000	50	1.00
10	3281892	25	2250	464	18.56
11	2852065	100	3000	353	3.53
12	162980	50	2250	572	11.44
13	0777935	100	0	253	2.53
14	5445356	100	2250	178	1.78
15	6007418	50	3600	152	3.04
16	5294675	25	2250	144	5.76
17	4549910	100	3000	87	0.87
18	5999831	25	2400	163	6.52
19	2002116	25	2400	54	2.16
20	6158273	50	4000	365	7.30
21	6176417	50	3600	293	5.86
22	4255178	50	3600	557	11.14
23	6230023	50	3600	99	1.98
24	0185310	75	2000	595	7.93
25	0437466	50	3000	488	9.76
26	1830237	50	2250	359	7.18
27	1996388	50	3000	146	2.92
28	2240525	50	1500	301	6.02
29	2366881	50	3000	287	5.74
30	4782476	50	2000	193	3.86
31	4916624	50	3000	201	4.02
32	5478242	100	2250	380	3.80
33	5480939	25	3000	292	11.68
34	5615377	100	4000	113	1.13
35	5796199	100	0	50	0.50
36	5892802	75	3000	366	4.88
37	5906785	50	2000	561	11.22
38	6176174	50	0	288	5.76

Clinical information related to drug administration for patients.

Since the 5-ASA dose was less than 2000 mg/day in sample numbers 3, 8, 13, 28, 35, and 38, they were not used for data analysis in this study.

Measurement of 6-TGN

The blood 6-TGN concentration was measured by Mitsubishi Chemical Medience Corporation using HPLC (Inertsil ODS-3 (3.0×100 mm, 4 µm). This is the modified method from the previous report [23].

Exposure of HapMap lymphocytes to drugs

AZA and 5-aminosalicylic acid (5-ASA) were obtained from MP Biomedicals, Inc. (Tokyo, Japan). They were dissolved in RPMI medium 1640 (Invitrogen, Van Allen Way Carlsbad, CA, USA) including 15% fetal bovine serum (Invitrogen) as vehicle for in vitro experiments. Suspensions of 30 different HapMap lymphocytes from Japanese subjects (Coriell, Camden, NJ, USA), at a concentration of 2.5×106 cells in 0.5 mL of RPMI medium 1640 (Invitrogen) including 15% fetal bovine serum, were incubated with 4.5 mL of AZA/5-ASA (100 µM AZA and 10 mM 5-ASA). The suspensions were sampled after 24 hours and used for the following assay. As controls, the respective HapMap lymphocytes were treated in the same way, but only with the vehicle RPMI medium 1640 including 15% fetal bovine serum (no drugs).

mRNA expression analysis

Total RNA was extracted from the AZA/5-ASA-treated and control cells using TRIzol Reagent (Invitrogen). Total RNA concentration was determined using a NanoDrop ND-1000 instrument (NanoDrop Technologies, Wilmington, DE, USA), and RNA quality was checked with an Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). Starting with 250 ng of total RNA, GeneChip U133 Plus 2.0 Array analysis (Affymetrix, Santa Clara, CA, USA) was conducted following the manufacturer's protocol, using the GeneChip 3′IVT Express Kit (Affymetrix) and the GeneChip Hybridization Wash and Stain Kit (Affymetrix). The raw expression signal obtained through a GeneChip Instrument System (Affymetrix) was normalized by scaling of the target signal to 100. For further data processing, probes were selected showing a P call in the detection step of the U133 Plus 2.0 Array among these selected probes; any probe set showing a ≥2-fold increase or decrease (AZA/5-ASA- treated vs. vehicle) of gene expression for a given sample was considered a significant change.

SNP 6.0 array analysis

DNA was extracted from the 30 different untreated HapMap lymphocytes using a QIAamp DNA mini Kit (QIAGEN), and the quality and concentration of the DNA samples were determined with a NanoDrop ND-1000 (NanoDrop Technologies). On the other hand, DNA of IBD patients was extracted with FALCO Biosystems (Kyoto, Japan). Starting with 500 ng of DNA, the 6.0 SNP Mapping Array assay was performed using the GeneChip 6.0 SNP Mapping Assay Kit (Affymetrix) following the manufacturer's protocol. Data for SNP typing were determined using the Genotyping Analysis Software (GTYPE) Version 4.1 of the GeneChip Instrument System (Affymetrix). Based on the 6.0 SNP Mapping Array genotyping information, an AA, AB, or BB call was obtained. In this analysis, a Birdseed v1 or v2 algorithm was included in the calculation [26].

MassARRAY analysis

DNA of IBD patients was extracted with FALCO Biosystems. Genotyping of SNPs (rs1142345, rs1800584, rs56161402, rs6921269, rs2842951, rs2842934, rs1800460, rs2518463, rs72552739, rs1800462, and rs3898137) in the TPMT gene was performed using the Sequenom MassARRAY technology platform (Sequenom, San Diego, CA, USA), with the iPLEX Gold Reagent Kit (Sequenom), SpectroCHIP ARRAYS, Clean Resin Kit (Sequenom), and HotStarTaq DNA Polymerase (QIAGEN, Hilden, Germany) based on the single-base extension reaction. Locus-specific PCR primers and allele-specific detection primers were designed using the MassARRAY Assay Design 3.1 software (Sequenom). Allele detection was performed using a MassARRAY Analyzer Compact MALDI-TOF mass spectrometer (Sequenom), a MassARRAY Nanodispenser (Sequenom), and the MassARRAY TYPER 3.4 software (Sequenom).

Correlation analysis of SNPs and 6-TGN risk ratios in IBD patients

We difined the risk ratio as the proportion of blood 6-TGN concentration and the dose of AZA intake (6-TGN concentration/AZA dose). The correlations between the AA, AB, and BB genotypes and 6-TGN risk ratios of clinical data in IBD patients were analyzed. For the correlation analysis, the genotypes of each sample from 32 individuals were plotted on the X-axis as AA = 0, AB = 1, and BB = 2, and 6-TGN risk ratios for the corresponding genes were plotted on the Y-axis. These for which all AA, AB, and BB genotypes were observed in at least 2 samples were selected for subsequent analysis, and r2 was calculated for the combination of each SNP and 6-TGN risk ratio. With this correlation analysis, SNPs with r2≥0.4 (p-values<0.01) with 6-TGN risk ratios were identified, and these were considered genetic biomarkers.

ExpressGenotyping analysis

Premature RNA expression analysis was performed using the ExpressGenotyping reactor and ExpressGenotyping analyzer (Japanese Patent Number: 4111985, Method for identifying a gene with allelic variation in gene expression (FILING DATA: December 15, 2005)). This new analysis technique is performed using the GIM algorithm [27]–[29]. The ExpressGenotyping reactor consists of a group of reaction steps where cDNA is synthesized preferentially from the premature RNA in the total RNA. The subsequent assay using the synthesized cDNA is essentially similar to the Affymetrix 6.0 SNP Mapping Array's assay, but instead of ordinary DNA typing, quantitative allele-specific expression information is obtained. The ExpressGenotyping analyzer combines the information from the DNA and cDNA typing and detects variations of premature RNA expression from each allele at particular SNPs. Consequently, comprehensive premature RNA expression assays using the ExpressGenotyping reactor and ExpressGenotyping analyzer can detect genes with differential allelic expressions.

Detection of drug-induced EAI by ExpressGenotyping analysis

To identify genetic biomarkers involved in the AZA/5-ASA-response, ExpressGenotyping analysis was performed as follows: 1) The ExpressGenotyping data outputs of drug-treated and vehicle-treated cells were obtained as a text file; 2) The respective vehicle and drugs files were combined with probe ID; 3) Probes with a heterozygous (AB) call within a gene region (3UTR, 5UTR, CDS, exon, and intron) were selected; 4) The expression ratios of premature RNA (vehicle sample (A allele/B allele or B allele/A allele) “vehicle EAI (max)”, drug sample (A allele/B allele or B allele/A allele) “drug EAI (max)” were calculated; 5) SNPs with an EAI ratio [calculated as drug EAI (max)/vehicle EAI (max)] of ≥2 were considered to show an EAI specific for drug exposure; and 6) EAI data of 30 test samples that were cell lines were collected. Next, AB heterozygous genotype analysis based on statistics was performed to identify responsive SNPs from the cells exhibiting the AB heterozygous genotype. SNPs with 3 or more heterozygous calls among the 30 test samples were used for this analysis. In this method, SNPs showing p-values<0.01 in either analysis of covariance of expression values or the t-test of the difference of inter-allelic expression ratios in cells under stimulation with test compounds and also exhibiting inter-allelic expression ratios of ≥2-fold were added to the list of test compound-responsive SNPs.

Comparative analysis between clinical and ExpressGenotyping data

The genetic biomarkers analyzed in IBD patients' blood samples were correlated with drug-induced EAI obtained by ExpressGenotyping analysis. In this analysis, a coincidence between genetic biomarkers of clinical data (r2 = 0.4 or more, p-value = 0.01 or less) of IBD patients and drug-induced EAI values of ≥2 was detected based on the statistics of the AB heterozygous genotype analysis.

Statistics

Risk ratio was assessed for each genotype group of patients. The data were compared using correlation analysis and t test. Prizm version 5.0 d (GraphPad Software Inc, La Jolla, CA) and SciPy (statistical library of Python) (http://www.scipy.org/) were used for all statistical analysis. All statistical significance was defined as P<.05.

Results

6-TGN blood concentration and 5-ASA dose

The relationship between 6-TGN blood concentration and the 5-ASA dose was analyzed. There was no positive correlation between 6-TGN blood concentration and low-dose 5-ASA (Figure 2a) or high-dose 5-ASA (Figure 2b) administration.

Figure 2

Frequency distribution of 6-TGN levels with the treatment of AZA and 5-ASA.

(A) 6-TGN levels of patients taking less than 3000 mg/day of 5-ASA. (B) 6-TGN levels of patients taking more than 3000 mg/day of 5-ASA. Relationships between 6-TGN levels and MCV (C) and between 6-TGN levels and WBC (D).

As shown in Figure 2c, mean corpuscular volume (MCV) was highly correlated with 6-TGN blood concentration. The white blood cell (WBC) level showed a significantly negative correlation with 6-TGN levels (Figure 2d).

Gene expression changes of mRNA were observed in mRNA expression analysis of HapMap lymphocytes from Japanese subjects. The scatter plots (Figure S1) showed significant changes in mRNA expression for the 100 µM AZA and 10 mM 5-ASA-treated vs. vehicle-treated HapMap lymphocytes for 24 hours. As the probe numbers of gene expression changed, by observing an average of 1257 up (change ≥2×)-regulated-probes (796 genes) and 815 down (change ≤0.5×)-regulated-probes (621 genes), 100 µM AZA and 10 mM 5-ASA treatment caused a marked effect on mRNA expression. GOTERM analysis of the Database by Annotation Visualization and Integrated Discovery (DAVID) showed the up-related terms (apoptosis, programmed cell death, and death etc.) and the down-related-terms (M phase, cell cycle phase, and cell cycle etc.) (http://david.abcc.ncifcrf.gov/). Based on the above results, exposure with 100 µM AZA and 10 mM 5-ASA for 24 hours was selected as the experimental condition for the subsequent ExpressGenotyping analysis of 30 different HapMap lymphocytes (Figure S2).

SNPs analysis in IBD patients

Correlation analysis of SNPs and 6-TGN risk ratios in IBD patients with SNP 6.0 Array analysis detected 1286 SNPs (516 genes) with r≥0.4 (p-values<0.01) (Figure S3). The GOTERM analysis of the database by DAVID identified cell adhesion, biological adhesion, and cell projection morphogenesis, etc. The top most significant SNPs were genotypes for rs3846502, rs3761769, rs16884434, and rs6897117 in SLC38A9 (Figure 3). SLC38A9 has been reported as a transporter gene. On the other hand, no correlations with SNPs in TPMT gene were obtained. Therefore, in order to perform a more detailed analysis of the TPMT gene, by MassARRAY analysis, the genotypes (AA, AB, or BB) of various TPMT SNPs were identified. TPMT is an enzyme that metabolizes AZA to 6-TGN, and some SNPs cause differences in metabolic activity. Four SNPs of TPMT have been analyzed (rs1800462, TPMT*2 [30]–[32]; rs1800460, TPMT*3B [33]; rs1142345 “A719G”, TPMT*3C [34]; rs1800584, TPMT*4 [35]) that affect 6-TGN blood concentration and, therefore, hematological toxicity. No difference in the distribution of SNPs was found, and genotype effects of the risk allele on hematological toxicity were not seen in all clinical samples measured (Figure S4). Other taq SNPs (rs56161402, rs6921269, rs2842951, rs2842934, rs2518463, rs72552739 and rs3898137) in TPMT showed no differences in these genotypes, as with the SNPs noted above. These SNPs of rs2842951, rs2518463, and rs3898137 observed in different distributions did not show significant correlations (rs2842951 (r = 0.08014, p-values = 0.66281), rs2518463 (r = 0.10855, p-values = 0.55427), and rs3898137 (r = 0.04124, p-values = 0.82268)). From the above results, differences in the 6-TGN risk ratios were found to not be effected by TPMT SNPs of IBD patients.

Figure 3

SLC38A9 SNPs of IBD patients.

SLC38A9 SNPs were determined by correlation analysis with genotypes and 6-TGN risk ratios in IBD patients.

Detection of drug-induced EAI

For the detection of AZA/5-ASA-induced EAIs, 8373 SNPs (3338 genes) with an EAI ratio ≥2 in ExpressGenotyping analysis (Figuer S5) were selected, and the numbers of AZA/5-ASA-induced EAIs for each HapMap lymphocyte are shown in Figure 4a. The dispersion diagram of the AZA/5-ASA-induced EAIs in the AB heterozygous genotype analysis for each chromosome is shown in Figure 4b. The most significant genes showing ≥3 fold-change of EAI were the SNPs of chromosome 11: GALNTL4 (rs2129523), chromosome 14: NPAS3 (rs12433161), chromosome 20: PHACTR3 (rs6015561) and chromosome 20: MACROD2 (rs8116068). In addition, the AB heterozygous genotype analysis of ExpressGenotyping data revealed reliable AZA/5-ASA-induced EAIs (64 SNPs (55 genes). The results are shown in (Figure S6). GOTERM analysis of the database with DAVID identified the regulation of apoptosis, regulation of programmed cell death, and regulation of cell death, etc. in the AB heterozygous genotype analysis. These EAI-related SNPs on ExpressGenotyping analysis may be important as genetic biomarkers affecting gene expression changes of each allele in predicting 6-TGN blood concentrations in patients taking AZA/5-ASA medication.

Figure 4

Number of drug-induced EAIs for 30 different HapMap lymphocytes by ExpressGenotyping analysis.

(A) For the drug-induced EAIs, SNPs with an EAI ratio ≥2 were selected. Dispersion diagram of drug-induced EAIs for each chromosome. (B) The drug-induced EAI of fold-change in the dispersion diagram detected with AB heterozygous genotype analysis by ExpressGenotyping. Details of the AB heterozygous genotype analysis are described below. SNPs with the AB heterozygous genotype in >3 of 30 HapMap lymphocytes were used to search for SNPs responsible for producing the drug-induced intracellular allelic imbalance. These analyses were performed using the SNPs for which the AB heterozygous genotype was observed in at least 3 samples.

Comparative data analysis by clinical data and ExpressGenotyping analysis

An attempt was made to identufy genetic biomarkers by combining 6-TGN risk ratios of clinical data and AZA/5-ASA-induced EAI of ExpressGenotyping analysis. In this comparative data analysis, SNPs in SLC38A9 showed a correlation with the 6-TGN risk ratios and AZA/5-ASA-induced EAI. The SNP rs6897117 of SLC38A9 was significant in the clinical data, and ExpressGenotyping analysis showed an EAI ratio ≥2 in GM18975 (Figure 5). The EAI ratio of SLC38A9 showed that the gene expression level of allele B was reduced, and it may increase the 6-TGN risk ratios. Therefore, these SNPs identified by ExpressGenotyping analysis might serve as genetic biomarkers for prediction of 6-TGN blood concentrations in patients taking AZA/5-ASA medication.

Figure 5

Results of drug-induced EAI for SNP (rs6897117) of SLC38A9.

The SNP of SLC38A9 shows a consistent relationship between 6-TGN risk ratios and drug-induced EAI.

Discussion

In the present study, new SNPs on SLC38A9 as genetic biomarkers to predict blood 6-TGN concentrations in patients taking AZA were identified. As the usage of AZA for IBD patients is now increasing, and the efficacy of AZA to maintain remission has been reported to be superb, it is very important that complications do not occur with AZA. In the analysis of TPMT SNPs of IBD patients in this study, the difference in the 6-TGN risk ratios was found to not have the risk allele related to hematological toxicity. Compared to the previous genetic biomarker, the new biomarker, SNPs of SLC38A9, was more sensitive.

In this study, 6-TGN concentration was not correlated with the dose of AZA treatment. Most of the studies have also reported no relationship between AZA dose and 6-TGN concentration [36]. This means that weight-based dosing of AZA does not guarantee therapeutic metabolite concentrations, and the efficacy of AZA metabolism might determine 6-TGN concentrations. As the number of WBCs and 6-TGN concentration showed a reverse correlation in the present study, prediction of blood 6-TGN concentration before AZA treatment will contribute to prevent leukopenia induced by 6-TGN.

Solute carriers (SLCs) are the largest group of transporters in the human genome [37], [38]. Many of the SLCs have important functions in cellular processes, including roles as transporters for amino acids and neurotransmitter cycling, while others are passive transporters or exchangers [38]–[40]. The SLCs have been categorized into at least 46 different families with varied biochemical properties [41], and two more families have recently been identified [42]. In mammals, the SLC proteins form four major phylogenetic groups, α-, β-, γ-, and δ-groups, with proteins in each group having a common evolutionary origin. The β-group includes the SLC32, SLC36, and SCL38 families, and it is the second largest phylogenetic cluster of amino acid transporters.

The SLC38 family consists of 11 members. The orphans SLC38A7-11 were only recently identified, and their substrate of transport has thus far not been characterized [43], while the other members, SLC38A1-6, are functionally classified as sodium-coupled neutral amino acid transporters, also known as SNATs. SNAT1 (SLC38A1, [44], [45]), SNAT2 (SLC38A2, [46]–[49]), and SNAT4 (SLC38A4, [50], [51]) have been further classified into system A, while SNAT3 (SLC38A3) and SNAT5 (SLC38A5) belong to system N transporters [52]–[54]. The system N/A classification depends on the functional properties and patterns of substrate recognition. SNAT6 (SLC38A6) [55] and the other more recently discovered members of the family have not been classified according to the N/A system. All SNAT genes, except SNAT4 [50], are expressed in the brain, with SNAT1 [44], [56] and SNAT2 [57] expressed both in astrocytes and neurons, while SNAT3 [52] and SNAT5 [58] are located in glial cells. No exact cellular localization of the protein or substrate profile has been reported for the orphan members of the SLC38 family. L-glutamine is a favored substrate for most of the members throughout the SLC38 family, and it is suggested that the transporters are part of the glutamine-glutamate cycle in the CNS [52], [54]. In this study, the genotype of SLC38A9 was found to be related to the blood 6-TGN concentration. The direct effect of SLC38A9 SNPs on 6-TGN blood concentration is still unclear, and further study will be necessary to reveal the pathway from SLC38A9 SNPs to the metabolism of AZA.

ExpressGenotyping analysis developed by HaploPharma Inc. is an innovative technology to detect critical pharmacogenetic SNPs and to provide genetic biomarkers to analyze them. ExpressGenotyping analysis enables the discovery of genetic biomarkers with a much smaller number of samples than with conventional approaches, such as genome-wide association study (GWAS), because these conventional approaches require complex procedures that usually need samples of hundreds or thousands of cases and a long time period, as well as being quite costly. In comparison to these approaches, ExpressGenotyping analysis can analyze the expression of each allele in a sample separately and, thus, reveal differences in allelic gene expression within the sample. Therefore, ExpressGenotyping analysis can be used to perform practical genetic biomarker discovery in samples with smaller numbers, at low cost. The basic functional principle of ExpressGenotyping is as follows. 1) Cell lines (HapMap lymphocytes) are treated with the drug of interest. 2) RNA and DNA are prepared from drug-treated and untreated cell lines. 3) Combined SNP typing/expression profiling on Affymetrix GeneChip microarrays is performed. 4) Analysis of expression profiles (drug-treated vs. vehicle-treated) and drug-induced EAI is then done. 5) Candidate genetic biomarkers for drug response prediction are identified. 6) Candidate genetic markers are validated with sample material from clinical studies. These technologies combined allow the genome-wide detection of variations in gene expression of individual alleles within one sample and eliminate the requirement for extensive inter-sample statistical analyses to identify a correlation between gene expression and genotype. In this report, the genetic biomarkers were identified by combining clinical data of IBD patients with AZA/5-ASA-induced EAI obtained by ExpressGenotyping analysis. Many SNPs in different genes were identified as candidate genetic biomarkers showing the relationship biologically (GOTERM analysis of DAVID). Of the genetic biomarkers related to the clinical 6-TGN blood concentration, SNPs in SLC38A9 were the most significant. This SLC38A9 gene might have an important role as a transporter for various drugs. Therefore, this SNP identified by ExpressGenotyping analysis might serve as a genetic biomarker for prediction of the 6-TGN blood concentration in patients on AZA/5-ASA.

These results show that ExpressGenotyping analysis is a useful way to identify individual differences in the incidence of drug response and ADRs, and it is expected to be a new method to accelerate the achievement of personalized medicine. In the future, the EAI EG method can be applied to the analysis of ex ante evaluation of immunosuppressive agents in the treatment of IBD, such as AZA.

Scatter plot analysis. In order to determine the experimental conditions for ExpressGenotyping analysis, scatter plot analysis was performed with the U133 Plus 2.0 Array. Data analysis was used to average signal values of *GM18940, GM18942, and GM18943 (*Sample ID in HapMap lymphocytes from Japanese individuals).

(TIFF)

Click here for additional data file.

30 different HapMap lymphocytes. The catalog ID indicates the database number for HapMap lymphocytes from Japanese subjects.

(DOCX)

Click here for additional data file.

Correlation analysis of SNPs and 6-TGN risk ratios in IBD patients. Criteria: r≥0.4, p-values<0.01.

(XLSX)

Click here for additional data file.

SNP genotypes of TPMT. TPMT 11 SNP genotypes were analyzed for 32 IBD patients.

(DOCX)

Click here for additional data file.

AZA/5-ASA-induced EAI detected with ExpressGenotyping analysis. Criteria: AZA/5-ASA-induced EAI≥2.

(XLSX)

Click here for additional data file.

AZA/5-ASA-induced EAI by AB heterozygous genotype analysis of ExpressGenotyping data. This list shows the AZA/5-ASA-induced EAI in more than one (AB genotype in >3 of 30) HapMap lymphocyte.

(DOCX)

Click here for additional data file.