Azathioprine with Allopurinol: Lower Deoxythioguanosine in DNA and Transcriptome Changes Indicate Mechanistic Differences to Azathioprine Alone

Article first published online 27 April 2017 Supplemental Digital Content is Available in the Text.

effects to standard-dose azathioprine or 6-MP. It is also used in a few patients with a poor response to azathioprine/6-MP, and a specific metabolite profile characterized by high methylated metabolite (MMP) and low thioguanine nucleotide (TGNs) levels in red-blood cells (RBC). 3,4 This "hypermethylation" profile may result from preferential methylation of 6-MP by thiopurine methyltransferase (TPMT) to MMP (Fig. 1), giving a poor clinical response and increased chance of hepatotoxicity. [5][6][7][8] Methylated thiopurine metabolites make an important contribution to thiopurine toxicity 9 ; patients with a hypermethylator profile or hepatotoxicity respond to LDAA with an increase in RBC TGNs, a reduction in RBC MMP, and disease remission. 4,10,11 In a recent open-label study, Pavlidis et al 12 treated a cohort of thiopurine-naive IBD patients with LDAA, dose-adjusted for TPMT status, and without prior selection based on RBC metabolite levels. Clinical response in that study was seen in 45/57 (79%) patients with Crohn's disease and 34/53 (64%) patients with ulcerative colitis/IBD-unclassified, with remission rates of 61% and 53% of patients, respectively. This is further evidence that the use of LDAA both in hypermethylators 3,4,10,11 and when used as a firstline therapy without metabolic screening 12 is an effective clinical strategy comparable to established guidelines for IBD. 13 The aim of this pilot study was to test the hypothesis that the AZA and LDAA treatment strategies are equivalent in terms of biological markers of thiopurine effect. As markers, we have investigated the utility of DNA-incorporated deoxythioguanosine (dTG DNA ) in comparison with RBC TGN measurement and reduction in neutrophil and white blood-cell counts. In addition, we have used transcriptome analysis to test the hypothesis that similar biological processes are altered by the different thiopurinedosing strategies.

Chemicals and Enzymes
The deoxythioguanosine (dTG) standard was from Carbosynth (Compton, United Kingdom); deoxyadenosine (dA) was from Sigma-Aldrich and deuterated methylmercaptopurine (internal standard) from Toronto Research Chemicals (Ontario, Canada). HPLC-grade formic acid was from Fisher Scientific (Loughborough, United Kingdom). All other reagents were from Sigma (Gillingham, United Kingdom). Calf intestinal alkaline phosphatase and nuclease P1 from Penicillium citrinum were from Sigma.

Study Population and Sample Collection
Clinical samples were provided from adult IBD patients treated with either AZA or LDAA at 8 centers in the UK (see Table 1, Supplemental Digital Content 1, http://links.lww.com/IBD/B503). Patients treated with LDAA were from East Surrey Hospital (ESH) where patients were normally treated with LDAA as a first-line therapy, and without preselection based on RBC metabolites. 12 The numbers of adult patients (ESH patient numbers in parentheses) being treated for IBD comprized 81 (43) with Crohn's Disease, 57 (28) with ulcerative colitis, 4 (1) with indeterminate colitis, and 3 (3) with microscopic colitis. At all centers, patients treated with AZA were treated according to standard clinical guidelines. 13 Of the 145 patients, 84 (50 AZA and 34 LDAA patients) were not prescribed other medication; the remaining 61 patients were receiving a variety of additional drugs, including mesalamine (5-aminosalicylic acid; 34 patients), infliximab (4 patients, 1 of whom was receiving mesalamine as well) and/or steroids (prednisolone; 3 in the LDAA group and 7 in the AZA group); other drugs included statins, antidepressants, and proton pump inhibitors. Data were analyzed only with respect to thiopurine treatment status, using patients treated with LDAA at ESH, and patients treated with AZA across all centers; to avoid bias, LDAA/AZA data from ESH were also analyzed separately. As part of routine clinical practice some participating centers requested RBC thiopurine metabolite measurements (RBC-TGN and RBC-MMP) at the same time as samples were taken for dTG DNA analysis; most (92%) of the samples with these data were from ESH.
Ethics Committee Centre, and each patient gave written informed consent.
Sample Preparation and LC-MS/MS Analysis of Deoxythioguanosine Incorporated in DNA (dTG DNA ) DNA was isolated from whole blood collected in EDTA tubes, using previously published methods, 14 and the ratio of dTG relative to dA in DNA (expressed as the number of moles of dTG/ 1 · 10 6 moles of dA) was analyzed using a recently developed HPLC-mass spectrometry method. 15 Briefly, chromatographic separation of dTG and dA was with a XSelect HSS T3 3.5 mm 4.6 · 100 mm (Waters) column with a VanGuard 3.5 mm 3.9 · 5guard column maintained at 308C. Analytes were eluted with mobile phases of 0.01 M aqueous 0.05% formic acid (A) and 0.05% formic acid in acetonitrile (B). The flow rate was 0.5 mL/minute, and the mobile phase system consisted of a starting condition of 1% buffer B increasing to 3% at 1.1 minutes, 8% at 2.4 minutes, and increasing to a maximum of 30% at 4.1 minutes then decreasing to 5% at 4.5 minutes, maintained until 5.5 minutes then decreasing to 1% for an equilibration period of 2.5 minutes. Standards ranged from 100 mg/mL dA and 2000 ng/mL dTG and were serially diluted in water to 0.39 mg/mL dA and 3.9 ng/mL dTG, respectively. High (40,000 ng/mL dA:400 ng/mL dTG), medium (4000 ng/mL dA:40 ng/mL dTG) and low (400 ng/mL dA:4 ng/mL dTG) matrix controls were prepared in triplicate. Standards, matrix controls, and samples were injected in a volume of 50 mL; for sample injections this is equivalent to 0.2 mg patient DNA per injection.

RNA-seq Methods
Patients for RNA-seq analysis were being treated with the appropriate drug, were in clinical remission at 12 weeks, and for whom a pretreatment blood sample was available; samples for analysis were selected "blind" from the pool of patient samples that met these criteria. At the time of selection, other clinical data for the patients were not available. RNA-seq of total RNA from whole blood, collected into Tempus RNA tubes, was performed by AROS Applied Biotechnology (Aarhus, Denmark). Sequencing libraries were prepared in 2 batches, LDAA patients and AZA patients, each batch consisting of 10 samples representing 1 pretreatment and 1 post-treatment sample (12 weeks after initiation of thiopurine therapy) from each patient. Input RNA was 400 ng from each sample, and libraries were prepared using the Illumina TruSeq Stranded mRNA LT Sample preparation kit after fragmentation using the alternative Illumina fragmentation protocol and RiboZero rRNA depletion. Libraries were quality controlled with respect to concentration and size profile. Size profiles (Agilent Bioanalyzer) from samples tested were 300 to 700 bp; library and library pool concentrations were determined using the KAPA qPCR kit in triplicate reactions with 3 independent 1 · 10 6 -fold dilutions of libraries. RNA-seq libraries were sequenced on a Hi-Seq 2500 using Illumina SBS v4 chemistry using a paired-end 100-bp protocol. Base calls were converted and demultiplexed according to the sequencing primer barcodes, using Illumina CASAVA software, into FASTQ files. Paired-end reads (100 bp) in FASTQ files were aligned to UCSC hg19 using Bowtie 2.2 16 and Tophat 2.0.13 17 running under Linux (Ubuntu). Alignments were counted using Bioconductor R packages "Genomi-cRanges" and "GenomicAlignments."

Statistical Analyses
All statistical analyses were conducted using R 3.1.2 and 3.3.0. 18 We regard these analyses as exploratory and have not applied corrections for experiment-wise error rate in rank correlation and rank sum tests applied to different variable and data subsets. Nonparametric methods (2-sided; Kendall's Rank Correlation, Asymptotic Wilcoxon Mann-Whitney Rank Sum Test, Kruskal-Wallis Test) were used for basic analyses of thiopurine metabolites and clinical correlates. Parametric methods (1-and 2-sample [Welch's] t tests) were used where data were normally distributed (Shapiro-Wilks test). For analysis of dTG DNA at different times after treatment initiation, patients were divided into 2 groups: treatment time for 4 to 167 days and .167 days; the boundary at 167 days (24 wk) represents times to achieve remission (#167 d) and maintenance of remission (.167 d) used in a recent clinical study. 19 Gene-count data from RNA-seq studies were analyzed using a mixed-effect approach in "edgeR" to derive treatment effects (post/pre) across all patients in each group (standard-dose azathioprine, AZA, and low-dose azathioprine with allopurinol, LDAA). 20 Genes .1.5-fold over-or underexpressed in response to drug treatment (Benjamini-Hochberg False Discovery Rate [BH FDR] , 0.05) were mapped separately to Gene Ontology (GO) biological processes using "goseq" in R. For the goseq analysis, a BH FDR of ,0.001 was applied to the LDAA data and ,0.05 to the AZA data. The lists of statistically significant GO terms returned by goseq were summarized graphically using the web resource Revigo, 21 using a similarity cutoff of 0.7, and Cytoscape. 22 Robustness of the goseq analyses was tested by applying a jack-knife procedure, running the goseq analysis on the edgeR output after excluding each patient in turn. Summary lists of genes differentially expressed in common with both treatments were mapped to GO terms using "BiomaRt" in R and the GO terms summarized using Revigo using a similarity score of 0.4.

Incorporation of dTG into DNA of Nucleated Blood Cells
Eighteen samples with data on white blood-cell counts were received from newly diagnosed patients with IBD at initiation of treatment with AZA or LDAA (9 patients per group); these patients were not prescribed additional medication. Twelve weeks after initiation of therapy, there was a significant reduction in neutrophil and white-blood cell counts ( Fig. 2; 1-sample t test, P # 0.003, n ¼ 18); 3 patients in the AZA group were leukopenic (white-blood cell counts , 4 · 10 9 /L). There was no significant difference in reduction of neutrophil or white-blood cell count between the 2 drug regimens (Welch's t test; P $ 0.64), although dTG DNA for patients in the LDAA group was significantly lower (Mann-Whitney test, P ¼ 0.03; Fig. 2).
For the dataset as a whole, dTG DNA was also lower in the patients receiving LDAA compared with those on AZA (Mann-Whitney test, Z ¼ 2.95, P ¼ 0.003; Fig. 3), but there was no significant difference in white-blood cell counts between the 2 treatments (Mann-Whitney, P ¼ 0.65). With respect to time of treatment, dTG DNA levels in patients receiving either AZA or LDAA for more than 24 weeks (.167 d; maintenance of remission) were lower than for patients who had received therapy for up to 24 weeks (4-167 d; achievement of clinical remission) FIGURE 2. Changes in white blood-cell counts (A), neutrophil counts (B), and levels of dTG DNA (C) in newly diagnosed patients at 12 weeks after therapy with AZA or LDAA. The data show a significant reduction in mean white-cell and neutrophil counts (neutrophils in brackets) from 9.3 (6.4) to 6.6 (4.4) as a result of therapy with AZA or LDAA (n ¼ 18; one-sample t-test on difference per patient, P # 0.001), but no significant difference between the AZA or LDAA treatments (t-test; P $ 0.64). However, there were significantly lower dTG DNA levels in the LDAA group (Mann-Whitney Rank Sum Test, Z ¼ 2.21, P ¼ 0.03). Units: neutrophil and white blood-cell counts · 10 9 /L; dTG DNA , moles dTG/10 6 moles dA.  19 Initial assessment of dTG DNA by time periods at 4 to 83, 84 to 167, 168 to 252, or .252 days after initiation of therapy suggested a decrease after 167 days (see Fig. 1, Supplemental Digital Content 1, http://links. lww.com/IBD/B503); therefore, periods 1 and 2, and 3 and 4 were combined for analysis. The distributions of treatment times were similar for AZA and LDAA-treated patients. Units: dTG DNA , moles dTG/10 6 moles deoxyadenosine.
(P # 0.01 for all data and ESH data alone; Fig. 3). There were no patients deficient in TPMT activity (,10 mU/L), and no relationship between dTG DNA and patient TPMT activity, either for the whole dataset or each drug treatment group alone (Kendall, P . 0.5).
Patients on LDAA received azathioprine at approximately one-quarter the dose of those on AZA (median 40.5 mg [range , and 125 mg [range  per day, averaged across a week, for LDAA and AZA, respectively). For patients receiving AZA, there was a significant correlation between dose and dTG DNA (P ¼ 0.012, Table 1), but this was marginal for LDAA patients (P ¼ 0.073). In some patients within both AZA and LDAA treatment groups, dTG DNA was below the level of detection; for the 12 patients with undetectable dTG DNA where RBC TGN measurements were available, RBC TGN levels ranged between 194 and 1016 pmol/8 · 10 8 RBC, indicating that undetectable dTG DNA was not a result of poor compliance. Furthermore, there were no differences in white-blood cell counts between patients with detectable dTG DNA and those without (Mann-Whitney, P ¼ 0.96, AZA and LDAA groups combined).

RBC Metabolites and Relationships with dTG DNA
With respect to the relationships between RBC-TGNs and -MMPs, there was no correlation between RBC-TGN and RBC-MMP for LDAA patients, but the converse was true for AZA. However, there was a significant positive correlation between prescribed azathioprine dose and RBC-MMP levels for LDAA but not for AZA. Conversely, RBC-TGN showed no correlation with thiopurine dose for either AZA or LDAA. Finally, no significant relationship was apparent between dTG DNA and RBC-TGNs or RBC-MMPs for patients on AZA or LDAA (Table 1).
RBC-TGN levels differed between patients in the 2 drug groups with patients on LDAA having higher median RBC-TGN levels (Mann-Whitney test, P ¼ 0.008), although the effect was not significant for the ESH data alone (P ¼ 0.062; Table 2).   Note that for one patient, the "pretreatment" sample for RNA was taken 11 days after initiation of therapy, but the analysis was unaffected by inclusion of this as a time-zero sample. Units: dTG DNA , moles dTG/10 6 moles deoxyadenosine (dA); RBC-TGN, Units: pmoles/8 · 10 8 RBC. LDAA patients were not receiving additional medication; in contrast, of the AZA patients, one was also being treated with cyclosporin, mesalamine and prednisolone, one with mesalamine and one was also being treated with infliximab. CD, Crohn's Disease; nd, not determined; UC, ulcerative colitis.
Although median and mean RBC-MMP levels were lower in the LDAA group, this difference was not statistically significant (P $ 0.26). In contrast to dTG DNA measurements, there was no difference between the early and late time periods of treatment for either RBC metabolite (Mann-Whitney Rank Sum Test, P $ 0.35). . GO profiling from RNA-seq data for thiopurine-treated patients: (A) biological processes significantly associated with genes underexpressed as a consequence of LDAA therapy; (B) sampled (jack-knife: each patient excluded in turn) data, biological processes significantly associated with underexpressed genes in LDAA-treated patients which were common to all subsampled sets; (C) biological processes significantly associated with over-expressed genes in LDAA-treated patients; (D) biological processes significantly associated with overexpressed genes which were common to all subsampled sets of LDAA-treated patients; in A to D, the BH FDR was set to P , 0.001 for GO categories. E, GO profiling from RNA-seq data for AZA-treated patients, showing biological processes significantly associated with underexpressed genes (BH FDR was set to P , 0.05 for GO categories). No biological processes were significantly associated with increases in gene expression in response to AZA, but only the GO cell compartments "cell periphery" and "plasma membrane," were significantly associated with increases in gene expression in response to AZA and these were common to all subsamples.

Transcriptome Changes as Markers of Biological Response to Thiopurine Therapy
For 10 newly diagnosed patients, blood samples for transcriptome analysis by RNA-seq were collected prior to thiopurine therapy, and a second sample collected from each patient after 12 weeks of treatment with either AZA alone or LDAA; for these patients, 5 treated with AZA and 5 with LDAA, no adverse drug reactions were recorded, and all were in clinical remission. The time of sampling was consistent with achievement of maximal levels of dTG DNA and 9 of the 10 patients had detectable dTG DNA ; mean and median dTG DNA levels were similar in the 2 groups. The patient without detectable dTG DNA was receiving LDAA and RBC-TGN levels indicated compliance ( Table 3). All LDAA patients were from 1 site (ESH) and were not on additional medication, whereas, the AZA patients were from 3 separate sites and some were receiving additional drugs (Table 3). For the LDAA patient group overall, 619 genes were repressed and 614 had higher expression in response to therapy (seeTables S2 and S3, Supplemental Digital Content 1, http://links.lww.com/IBD/B503); applying a BH FDR of ,0.001 for Gene Ontology (GO) analysis of gene lists, underexpressed genes were significantly associated with biological processes of immune defense ( Fig. 4; see Table S4, Supplemental Digital Content 1, http://links.lww.com/IBD/B503). Subsampling by reanalysis after exclusion of each patient in turn revealed a core set of biological-process GO categories comprising "defence response" (GO:0006952), "defence response to other organism" (GO:0098542) and "immune response" (GO:0006955). GO biological-process categories significantly associated with overexpressed genes in the LDAA-treated patients were more diverse, with subsampling giving RNA metabolic processes and processes of macromolecular localization and organization as the GO category core sets for biological processes ( Fig. 4; see Table S5, Supplemental Digital Content 1, http://links.lww.com/IBD/B503). Within the top 5 genes repressed by LDAA (lowest P values), SNRK-AS1 is potentially involved in colon-cell proliferation, PLTP modulates inflammation and immune responses and C1orf56 (MENT) is a tumor modifier associated with lymphoma; the top 5 overexpressed genes were associated with activation of host-defense responses and hematopoiesis (see Tables S2 and S3, Supplemental Digital Content 1, http://links.lww.com/IBD/B503).
For the AZA-treated group, 464 genes were repressed and 165 had higher expression in response to therapy (see Tables S6 and S7, Supplemental Digital Content 1, http://links.lww.com/IBD/B503). In these gene lists, there was no significant representation in any GO categories at a BH FDR of ,0.001. At an FDR of ,0.05, underexpressed genes were enriched in diverse GO categories, particularly those concerned with ion transport and homeostasis, proteolysis and protein catabolism, and cell-cycle regulation ( Fig. 4; see Table S8, Supplemental Digital Content 1, http://links.lww.com/IBD/B503). The RNA-seq data for the AZA patient group as a whole were less consistent than the LDAA group; 3 of the 5 jack-knife subsamples for underexpressed genes had no significant GO categories. Overexpressed genes in the complete AZA dataset fell into only 2 GO categories: plasma membrane and cell periphery cell compartments, with no representation in biological-process GO categories; these GO cell-compartment categories were also common to all subsamples. However, overexpressed genes in 3 of the 5 subsamples were involved in regulation of cell activation and immune system biological processes (see Table S9, Supplemental Digital Content 1, http:// links.lww.com/IBD/B503).
The RAC1 GTPase is thought to be a target of thiopurines 23 but was not represented in the lists of genes differentially expressed in response to therapy (see Tables S2, S3, S6, and S7, Supplemental Digital Content 1, http://links.lww.com/IBD/B503). Furthermore, GO terms associated with RAC1 function (http://www.genecards.org/cgibin/carddisp.pl?gene¼RAC1) were not included in the lists of GO terms significantly associated with differentially expressed genes.

DISCUSSION
These data confirm the measurable incorporation of thioguanine nucleotides into the DNA of nucleated blood cells of patients with IBD on thiopurine therapy. A clinical effect of thiopurines is the reduction of neutrophil and white cell counts [25][26][27] which was clearly demonstrated with no difference between AZA and LDAA treatment groups; in contrast, dTG DNA levels were lower in LDAA-treated patients, as shown by the treatment group sampled after 12 weeks of therapy and the full study dataset. Given correlations between thiopurine dose and dTG DNA within treatments, lower dTG DNA levels in the LDAA group may reflect the lower thiopurine dose. These data suggest that the biological effects of AZA or LDAA are independent of dTG DNA levels in nucleated blood cells, in agreement with other studies showing no correlation between dTG DNA and white-blood cell counts. 28 Maximum dTG DNA levels were in patients sampled within 24 weeks of the start of therapy, with lower levels in patients sampled at longer times after therapy initiation; lower levels of dTG DNA on maintenance of remission may result from adaptation to thiopurines with respect to pathways of metabolism and rates of cell renewal in relevant target tissues.
In contrast to a previous study, 28 we found no significant correlation between dTG DNA and RBC-TGNs or RBC-MMPs. Furthermore, RBC-TGN levels were higher in LDAA-treated patients, in contrast to the lower dTG DNA levels in this patient group. Our current understanding is that inhibition of xanthine oxidase by allopurinol reduces the metabolic inactivation of thiopurines, allowing more conversion to TGNs. 29 MMPs, methylation products of TPMT activity, decrease in the presence of allopurinol by an unknown mechanism; TPMT inhibition by an allopurinol metabolite 30 or the 6-MP metabolite thioxanthine, 29,31,32 have been suggested, and depletion of phosphoribosylpyrophosphate or competitive inhibition of membrane transport processes are possible alternative mechanisms. 33 Clearly, in our study the RBC-TGN and RBC-MMP levels varied between the AZA and LDAA treatment groups as expected.
It may be paradoxical that dTG DNA levels in nucleated blood cells of LDAA-treated patients were low despite higher levels of RBC TGNs compared to patients treated with AZA. Nucleoside triphosphates are toxic to cells, 34 and homeostatic mechanisms may counter raised TGTP levels by conversion of nucleoside triphosphates to monophosphates, as we have found with TGTP in vitro (Coulthard, McGarrity, Redfern; unpublished data), thus reducing the incorporation of dTG into DNA despite higher TGN levels. However, any interpretation of how dTG DNA relates to RBC TGNs has to be speculative given our limited understanding of thiopurine metabolism in white blood cells and transfer of thiopurine metabolites between cell and tissue compartments. Indeed, the lack of a relationship between dTG DNA and RBC-TGN has echoes of studies showing a lack of correlation between RBC-TGN levels and biological or clinical responses. [35][36][37] RBCs have low IMPDH activity 38,39 and RBC-TGNs must, therefore, be derived from liver thioguanine metabolized from azathioprine, or remain from immature erythrocytes. Further studies are needed to test the relationships between dTG DNA and biological and clinical responses in IBD, and how dTG DNA levels are affected by the dynamics of thiopurine and purine metabolism in different cell populations.
Although thiopurines are used extensively as immunosuppressants there is uncertainty of the mechanisms involved, a question that can be addressed by transcriptome analysis of peripheral blood with respect to immune system markers. Inhibition of GTP signaling, specifically the RAC1 GTPase which is elevated in expression in activated T-cells, 40 may be one mechanism, at least in vitro. 23,41,42 However, evidence that this underlies clinical responses is lacking 43 and interference in RAC1 signalling was not apparent from our RNA-seq data. We found a clear response with respect to the downregulation of genes involved in host-defense mechanisms against microorganisms in LDAA patients. This may result from gene downregulation in relevant immune cells, or from a thiopurine-induced reduction of specific immune-cell populations. Whatever the mechanism, this is consistent with a role of microorganisms and host-defense interactions in the etiology of IBD. The wide range of dTG DNA values despite a fairly homogenous LDAA dataset with respect to GO categories affected implies that dTG DNA may not be a good biological marker of response.
Despite the similar range of dTG DNA levels to the LDAA samples, the AZA-treated patients did not show a consistent transcriptome-response profile. While more extensive studies are needed to test for potentially confounding effects of bias from the treatment site, sample batch with respect to the LDAA/AZA comparisons or the effects of additional drugs on some patients in the AZA group, the Gene Ontology analysis, both for the shared and the full lists of differentially expressed genes, and the small number of shared genes up-or downregulated by either treatment, support the view that there are differences between LDAA and AZA in their underlying therapeutic mechanisms. An important caveat for the interpretation of the LDAA data is that allopurinol may have anti-inflammatory properties, independently of its effects on thiopurine metabolism, by scavenging free radicals, or reducing the T-cell activation 44 involved in the pathogenesis of IBD. 45 Thus, by providing complementary anti-inflammatory effects, the therapeutic value of allopurinol for IBD may extend beyond merely improving thiopurine pharmacokinetics.
Two differentially expressed genes may warrant further investigation. Firstly, the downregulation of GUCY2C was common to both treatments. This gene is well known for its role in endotoxin-mediated diarrhoea, 24 is a potential colorectal cancer marker, 46 and may be a link between obesity and colorectal cancer. 47 Although GUCY2C expression has been detected in CD34 + progenitor cells in peripheral blood, 48 a possible reflection of immune system status, the presence of GUCY2C transcripts in peripheral blood may also be a consequence of IBD disease processes. The detection and significant downregulation of GUCY2C transcripts in the LDAA and AZA patient groups raises questions about the potential of GUCY2C as a marker of IBD disease status and response to therapy with thiopurines or other drugs. Secondly, MENT is associated with lymphomagenesis in mice, and decreased methylation activity upregulates MENT expression. 49 Noting that thiopurines decrease DNMT activity, 50 and aside from the possibility that lower dTG DNA with LDAA may reduce secondary cancer risk, the evidence for downregulation of MENT by LDAA suggests that the elevated lymphoma risk associated with thiopurine therapy in elderly patients 51 may be avoided using LDAA rather than standard-dose azathioprine.
Whatever the mechanism of LDAA versus AZA, recent data from a small randomized clinical trial suggest that patients treated with an allopurinol and azathioprine combination have better response rates than those on AZA. 52 Lower dTG DNA and MMP levels in LDAA-treated patients suggest that LDAA as a first-line therapy may be the most appropriate for IBD patients, obviating the need for preselection on the basis of RBC metabolites. 12 Further studies of the potential beneficial role of LDAA in different clinical contexts, such as a reduction in cancer risk, and to identify the precise mechanisms of action, are badly needed, but there is increasing evidence that this could be the initial treatment of choice for patients with IBD for whom immunosuppression therapy is appropriate.